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.

Melampsora larici-populina causes serious poplar foliar diseases called rust worldwide. Salicylic acid (SA) and jasmonic acid (JA) are important phytohormones that are related to plant defence responses. To investigate the transcriptome profiles of SA- and JA-related genes involved in poplar rust interaction, two tolerant poplars and one intolerant poplar were selected for this study. Weighted gene coexpression network analysis (WGCNA) was applied to characterize the changes in the transcriptome profiles and contents of SA and JA after infection with the virulent E4 race of M. larici-populina. In response to infection with the E4 race of M. larici-populina, tolerant symptoms were correlated with the expression of genes related to SA and JA biosynthesis, the levels of SA and JA, and the expression of defence-related genes downstream of SA and JA. Tolerant poplars could promptly regulate the occurrence of defence responses by activating or inhibiting SA or JA pathways in a timely manner, including regulating the expression of genes related to programmed cell death, such as Kunitz-type trypsin inhibitor (KTI), to limit the growth of E4 and protect themselves. WGCNA suggested that KTI might be regulated by a Cytochrome P450 family (CYP) gene. Some CYPs should play an important role in both JA- and SA-related pathways. In contrast, in intolerant poplar, the inhibition of SA-related defence signalling through increasing JA levels in the early stage led to continued inhibition of a large number of plant–pathogen interaction-related and signalling-related genes, including NBS-LRRs, EDS1, NDR1, WRKYs, and PRs. Therefore, timely activation or inhibition of the SA or JA pathways is the key difference between tolerant and intolerant poplars.

Background: Recently, pine wood nematode (PWN, Bursaphelenchus xylophilus) has been found in the extreme cold area of northeast China. The third-stage dispersal juvenile (DJ3) of PWN, which is a long-lived stress-resistant stage, plays an important role in the process of PWN spreading to low-temperature areas, as this stage can survive under unfavorable conditions. Results: Weighted correlation network analysis (WGCNA) was used to analyze the expression patterns of 15,889 genes included in 21 RNA-Seq results of PWN at DJ3 and the other 6 different stages, and a total of 12 coexpression modules were obtained. Among them, the magenta module has the highest correlation with DJ3, which included a total of 652 genes. KEGG enrichment analysis showed that most of the genes in the magenta module were involved in metabolic processes, which were related to autophagy and longevity regulation. These pathways included starch and sucrose metabolism, which contains trehalose metabolism. To explore the function of trehalose in DJ3 formation and survival under − 20 °C, a trehalose-6-phosphate synthase encoding gene (Bx-tps), a trehalose-6-phosphate phosphatase encoding gene (Bx-tpp) and 7 trehalase encoding genes (Bx-tres) were identified and investigated. The expression of these 9 genes was related to the formation of DJ3. A treatment under − 20 °C induced the accumulation of trehalose. The survival rate of DJ3 at -20 °C reduced after silencing of any of these trehalose metabolism genes. Further analysis suggested that two trehalose synthesis genes were highly correlated with DJ3 and might be involved in autophagy by regulating with energy conversion related genes. Conclusions: The above results indicated that trehalose metabolism promotes DJ3 formation and helps DJ3 survive at -20 °C. Although trehalose accumulation is favorable for DJ3 to cope with low-temperature stress, multiple trehalose metabolism genes need to work together. There may be a multi-path regulated physiological process involving trehalose synthesis genes under low-temperature stress resistance. This physiological process may regulate the formation and maintenance of DJ3 through autophagy and energy conversion.

Neofusicoccum laricinum is the causal agent of larch shoot blight, a fungal disease affecting several species of larch. It causes severe damage, including stunting and mortality. This study aims to address the severe impact of larch shoot blight by evaluating the effect of farrerol on the inhibition of Neofusicoccum laricinum in Larix olgensis. We used LC-MS/MS and weighted gene co-expression network analysis to investigate farrerol’s effects on Neofusicoccum laricinum and identify associated genes in resistant and susceptible larch. Our study identified significant differences in metabolite profiles between resistant and susceptible cultivars, with higher concentrations of farrerol showing complete inhibition of N. laricinum. Additionally, specific genes associated with farrerol content were up-regulated in resistant larch. Farrerol at higher concentrations completely inhibited N. laricinum, showing a strong correlation with increased disease resistance. This research suggests that farrerol enhances disease resistance in larch and provides a foundation for developing disease-resistant larch varieties based on antifungal metabolite traits.

Larch shoot blight, caused by Neofusicoccum laricinum, threatens global larch resources, while conventional chemical control is constrained by pollution and resistance. To address this gap, we integrated metabolomics, transcriptomics, and antifungal efficacy assays to identify Fraxetin, a disease-induced phytoalexin, and to elucidate its antifungal activity and mechanism. Metabolomics showed infection-triggered accumulation of Fraxetin in resistant Larix olgensis shoots. Antifungal experiments showed that within the range of 68–1088 μg/mL, the optimal antifungal concentration was 1088 μg/mL. When inoculated larches were treated with 1088 μg/mL Fraxetin, the maximum inhibition rate of pathogen growth reached 66.67% within 12 days, and the symptoms of the treated plants were alleviated. Transcriptomics revealed activation of damage responses, disruption of oxidative homeostasis, and compromised membrane integrity in the pathogen under Fraxetin treatment. Physiological measurements confirmed increased lipid peroxidation, redox collapse, membrane leakage, and reduced fungal viability. These findings indicate a lipid peroxidation–mediated oxidative–membrane mode of action and support the potential of plant-derived Fraxetin for more sustainable management of larch shoot blight.

Neofusicoccum laricinum, the pathogen responsible for larch shoot blight, is a hemibiotrophic pathogen. A hypersensitive reaction of plants does not inhibit the growth of the pathogen, while phytoalexin is an antifungal compound secreted by plants, which can directly destroy the cells of pathogens and help plants achieve resistance. This study aimed to investigate the chemical defense mechanisms in resistant larch, identify key biocontrol agents, and assess their potential for field application. By integrating multi-omics technologies with time-dependent models and dose–response curve analysis, the accumulation and antimicrobial properties of rutaevin were examined. Its application potential was verified through cross-species resource screening and field trials. The results revealed that rutaevin acts as a phytoalexin against larch shoot blight. It began to accumulate linearly 0.62 days after pathogen exposure, with its antifungal activity demonstrating a dose-dependent response, achieving 100% inhibition at 0.5 mg/mL. The activation of terpene metabolic pathways in disease-resistant plants resulted in a significant increase in rutaevin content compared to susceptible plants. Cross-species screening showed that the highest concentration of this compound is found in the fruit of Evodia rutaecarpa var. rutaecarpa, with its crude extract exhibiting strong field efficacy. The findings provide theoretical and technical support for disease-resistant breeding and the development of plant-derived fungicides.