Abstract
Pine wood nematode (PWN) causes serious diseases in conifers, especially pine species. To investigate the transcriptomic profiles of genes involved in pine-PWN interactions, two different pine species, namely, Pinus thunbergii and P. massoniana, were selected for this study. Weighted gene coexpression network analysis (WGCNA) was used to determine the relationship between changes in gene expression and the PWN population after PWN infection. PWN infection negatively affects the expression of most genes in pine trees, including plant defense-related genes such as genes related to plant hormone signal transduction, plant-pathogen interactions, and the MAPK signaling pathway in plants. However, the expression of chalcone synthase genes and their related genes were proportional to the changes in nematode populations, and chalcone synthase genes were dominant within the coexpression module enriched by genes highly correlated with the nematode population. Many genes that were closely related to chalcone synthase genes in the module were related to flavonoid biosynthesis, flavone and flavonol biosynthesis, and phenylpropanoid biosynthesis. Pine trees could actively adjust their defense strategies in response to changes in the number of invasive PWNs, but the sustained expression of chalcone synthase genes should play an important role in the inhibition of PWN infection.