Transcriptomic response of Pinus massoniana to infection stress from the pine wood nematode Bursaphelenchus xylophilus

Abstract The pinewood nematode (PWN) Bursaphelenchus xylophilus is a forestry quarantine pest and causes an extremely dangerous forest disease that is spreading worldwide. Due to the complex pathogenic factors of pine wood nematode disease, the pathogenesis is still unknown. B. xylophilus ultimately invades a host and causes death. However, little is known about the defence-regulating process of host pine after infection by B. xylophilus at the molecular level. Therefore, we wanted to understand how Pinus massoniana regulates its response to invasion by B. xylophilus. P. massoniana were artificially inoculated with B. xylophilus solution, while those without B. xylophilus solution were used as controls. P. massoniana inoculated with B. xylophilus solution for 0 h, 6 h, 24 h, and 120 h was subjected to high-throughput sequencing to obtain transcriptome data. At various time points (0 h, 6 h, 24 h, 120 h), gene transcription was measured in P. massoniana inoculated with PWN. At different time points, P. massoniana gene transcription differed significantly, with a response to early invasion by PWN. According to Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, P. massoniana response to PWN invasion involves a wide range of genes, including plant hormone signal transformation, flavonoid biosynthesis, amino sugar and nucleoside sugar metabolism, and MAPK signalling pathways. Among them, inoculation for 120 hours had the greatest impact on differential genes. Subsequently, weighted gene coexpression network analysis (WGCNA) was used to analyse transcriptional regulation of P. massoniana after PWN infection. The results showed that the core gene module of P. massoniana responding to PWN was “MEmagenta”, enriched in oxidative phosphorylation, amino sugar and nucleotide sugar metabolism, and the MAPK signalling pathway. MYB family transcription factors with the highest number of changes between infected and healthy pine trees accounted for 20.4% of the total differentially expressed transcription factors. To conclude, this study contributes to our understanding of the molecular mechanism of initial PWN infection of P. massoniana. Moreover, it provides some important background information on PWN pathogenic mechanisms