Transcriptomic analysis reveals differentially expressed genes associated with pine wood nematode resistance in resistant Pinus thunbergii

Pine wilt disease (caused by the nematode Bursaphelenchus xylophilus) is extremely harmful to pine forests in East Asia. As a low-resistance pine species, Pinus thunbergii is more vulnerable to pine wood nematode (PWN) than Pinus densiflora and Pinus massoniana. Field inoculation experiments were conducted on PWN-resistant and -susceptible P. thunbergii, and the difference in transcription profiles 24 h after inoculation was analyzed. We identified 2603 differentially expressed genes (DEGs) in PWN-susceptible P. thunbergii, while 2559 DEGs were identified in PWN-resistant P. thunbergii. Before inoculation, DEGs between PWN-resistant and PWN-susceptible P. thunbergii were enriched in the REDOX (Oxidation-Reduction) activity pathway (152 DEGs), followed by the oxidoreductase activity pathway (106 DEGs). After inoculation with PWN, however, the opposite was observed; DEGs were enriched in the oxidoreductase activity pathway (119 DEGs), followed by the REDOX activity pathway (84 DEGs). Before inoculation, according to the metabolic pathway analysis results, we found more genes upregulated in phenylpropanoid metabolic pathways and enriched in lignin synthesis pathways; cinnamoyl-CoA reductase-coding genes related to lignin synthesis were upregulated in PWN-resistant P. thunbergii and downregulated in PWN-susceptible P. thunbergii, and the lignin content was always higher in resistant than in susceptible P. thunbergii. These results reveal distinctive strategies of resistant and susceptible P. thunbergii in dealing with PWN infections.