Ya Liu
Qinzheng Zhou
Zheng Wang
Huiming Wang
Guiheng Zheng1
Quan Lu- 1Key Laboratory of Forest Protection of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, China
- 2State Key Laboratory of Tree Genetics and Breeding, Institute of Ecological Conservation and Restoration, Chinese Academy of Forestry, Beijing, China
Ophiostomatoid fungi exhibit a complex relationship with bark beetles; exhausting of host tree defenses is traditionally regarded as one of the key benefits provided to beetle vectors. Ophiostoma bicolor is one of the dominant species of the mycobiota associated with Ips genus bark beetles which infect the spruce trees across the Eurasian continent. Host spruce trees resist fungal invasion through structural and inducible defenses, but the underlying mechanisms at the molecular level, particularly with respect to the interaction between bark beetle-associated fungi and host trees, remain unclear. The aim of this study was to observe the pathological physiology and molecular changes in Picea koraiensis seedlings after artificial inoculation with O. bicolor strains (TS, BH, QH, MX, and LWQ). This study showed that O. bicolor was a weakly virulent pathogen of spruce, and that the virulent of the five O. bicolor strains showed differentiation. All O. bicolor strains could induce monoterpenoid release. A positive correlation between fungal virulence and release of monoterpenoids was observed. Furthermore, the release rate of monoterpenoids peaked at 4 days post-inoculation (dpi) and then decreased from 4 to 90 dpi. Transcriptomic analysis at 4 dpi showed that many plant-pathogen interaction processes and mitogen-activated protein kinase (MAPK) metabolic processes were activated. The expression of monoterpenoid precursor synthesis genes and diterpenoid synthesis genes was upregulated, indicating that gene expression regulated the release rate of monoterpenoids at 4 dpi. The enriched pathways may reveal the immune response mechanism of spruce to ophiostomatoid fungi. The dominant O. bicolor possibly induces the host defense rather than defense depletion, which is likely the pattern conducted by the pioneers of beetle-associated mycobiota, such as Endoconidiophora spp.. Overall, these results facilitate a better understanding of the interaction mechanism between the dominant association of beetles and the host at the molecular level.
Introduction
Spruce is an evergreen woody plant that is distributed in cold temperate and subalpine regions. The genus consists of approximately 50 species, including 20 species and five varieties found in China (Chen, 2013). In recent years, pest outbreaks have been accelerated by climate change and increasing international trade. Spruce forests in the Northern Hemisphere, especially in Europe, have been severely affected by bark beetles of genus Ips (Coleoptera: Curculionidae: Scolytinae) (Hlásny et al., 2021). Bark beetles and fungi have formed complex associations during long-term co-evolution, resulting in the ability for the beetles to carry ophiostomatoid fungi (Six, 2012). The combined impact of the association between beetle pests and ophiostomatoid fungi has been studied for a long time (Kirisits and Konrad, 2004; Wingfield et al., 2016) and has received attention in China since the turn of the century, in view of the successful invasion of Dendroctonus valens from North America and the epidemic outbreak of indigenous Ips spp. and Dendroctonus armandii (Lu et al., 2008; Sun et al., 2013; Wang et al., 2022).
Both beetle and fungal partners have developed from occasional and sporadic to relatively stable and intimate relationship during their historical co-occurrence in the same habits (Bleiker and Six, 2009; Biedermann and Vega, 2020; Chakraborty and Roy, 2021; Singh et al., 2021). The benefits of beetle vectors to fungi are evident; for example, beetles can bring the fungi to new hosts and create inoculation holes. However, the benefits of fungi to beetles are unclear and sometimes disputed (Six and Wingfield, 2011). Nevertheless, substantial evidences from studies in nutrition, chemical ecology, and molecular biology have supported that fungi (particularly ophiostomatoid fungi) can enhance the colonization success of beetles in the following ways: by producing ergosterol, which is indispensable for beetle development (Bentz and Six, 2006), improving food quality and availability (Davis et al., 2019; Dzurenko and Hulcr, 2022), detoxifying tree defense compounds (DiGuistini et al., 2011; Wadke et al., 2016; Itoh et al., 2018; Davis et al., 2019; Zhao et al., 2019a), attracting or repelling other individuals through semiochemicals (Zhao et al., 2015, 2019b; Kandasamy et al., 2019), or competitive exclusion of beetle pathogens (Davis et al., 2019). In addition, many members of the ophiostomatoid fungi are pathogens of forests and crops worldwide. Notorious diseases caused by ophiostomatoid fungi include Dutch elm disease, laurel wilt, and oak wilt on broadleaf trees, and stem canker stain and black root disease on coniferous trees, leading to severe ecological and economic losses (Henry, 1944; Brasier, 1986; Harrington et al., 2008; Santini and Faccoli, 2015). Artificial inoculation of certain ophiostomatoid fungi can generate inner bark necrosis and lead to sapwood blue staining and drying within a few months, severely hindering water transportation, affecting tree growth, and even causing the death of host plants (Yamaoka et al., 1998; Davydenko et al., 2017; Devkota and Eckhardt, 2018).
Conifers have evolved a combination of structural and inducible defenses against bark beetles and their associated fungi (Franceschi et al., 2005; Jones and Dangl, 2006; Celedon and Bohlmann, 2019). Studies on the chemical defense of conifers against bark beetle-ophiostomatoid fungal symbionts have shown that terpenes are toxic to the symbionts (Phillips and Croteau, 1999; Raffa et al., 2005; Keeling and Bohlmann, 2006; Pan et al., 2018; Celedon and Bohlmann, 2019; Ullah et al., 2021). The phloem tissue of Pinus spp. contains constitutive monoterpenes that provide immediate resistance to bark beetle attacks (Boone et al., 2011; Raffa et al., 2013); local concentrations of terpenes increase rapidly after beetle attack. When the induced terpene concentration levels exceed the beetle’s physiological tolerance threshold, they inhibit or repel beetles that arrive subsequently and alter the growth of fungi associated with beetles that are already present (Raffa et al., 2005; Erbilgin et al., 2017). In vitro experiments also confirmed that monoterpenes, such as pinene and limonene, can inhibit the growth of fungi (Zhao et al., 2019b; Fang et al., 2020; Wang et al., 2020a). Monoterpenes, such as α-pinene, with stronger growth inhibitory ability against fungi, have a stronger repelling effect on bark beetles (Fang et al., 2020). Long-term production of monoterpenoids consumes a large amount of carbon, which is essential for tree growth (Trapp and Croteau, 2001; Dai et al., 2022), impeding the growth of pest-infested trees. In addition to terpenes, phenolics seem less important in the induced defense against bark beetle-associated fungal infestation (Erbilgin et al., 2017).
RNA-seq facilitates genome-wide assessment of gene expression profiles in plants under various conditions, such as biotic and abiotic stresses (Conesa et al., 2016). At present, more and more studies have used this method to analyze host plant resistance (Nibedita and Jolly, 2017). In a transcriptomic comparison between healthy Gastrodia elata and those infected with Penicillium oxalicum, 10 potential resistance genes were identified, involving plant hormone signal transduction, jasmonic acid signaling, and plant-pathogen interaction pathways, revealing that the immune response mechanism of G. elata to fungal disease is a complex biological process (Wang et al., 2020b). In poplars infected with canker disease in the early stage, most of the carbon metabolism and transportation genes, aquaporin genes, and genes related to secondary metabolites and phenylpropane biosynthesis pathway were inhibited, while the expression of resistance genes was promoted (Li et al., 2019). The transcript abundance of Norway spruce leucoanthocyanidin reductase (LAR) genes increased significantly during Endoconidiophora polonica infection (Hammerbacher et al., 2014). Anthocyanin monomers can mitigate the effects of various abiotic stresses such as ultraviolet radiation and ozone by reducing oxidative stress and inhibiting bacterial growth and fungal spore germination (Jaakola et al., 2004; Karonen et al., 2006).
It is interesting to understand whether the beetle-vectored fungi cause plant disease and whether the plants develop resistance against them. Ophiostoma bicolor is one of the dominant associate fungi with Ips spp. and Dendroctonus micans, which attacked various spruce species in Europe and north China (Giordano et al., 2013; Nève Repe et al., 2018; Chang et al., 2019, 2020; Chakraborty et al., 2020; Wang et al., 2020c,2021). Previous studies have confirmed that the pioneer species in the associated fungal community of Ips, like Endoconidiophora spp., showed strong virulence, which could cause severe host necrosis and induce a drastic resistance response of host conifers. However, O. bicolor could cause a small amount of spruce sapwood blue-staining and drying, and form small necrotic lesions on the inner bark near the inoculation site, accompanied by bits of resin outflow, showing moderate to weak levels of virulence from a series of pathogenicity tests in Europe and Japan (Christiansen and Solheim, 1994; Yamaoka et al., 2000; Sallé et al., 2005; Repe et al., 2015). The pathogenicity of O. bicolor isolated from China is unknown, and whether O. bicolor can induce the host to produce a resistance response like the pioneer species need to be further investigated. In addition, the molecular mechanisms underlying host resistance to ophiostomatoid fungi remain unclear.
In this study, the representative strains of O. bicolor were inoculated on 4-year-old Picea koraiensis to observe the differences in virulence between various strains, given that different strains of the same pathogen species may have different levels of virulence (Lieutier et al., 2003; Sallé et al., 2005). Changes in host defense-related metabolites after inoculation were also analyzed and the defense mechanisms of spruce against O. bicolor were investigated at the molecular level.
Materials and methods
Plant materials and fungal strains
To compare the virulence of O. bicolor strains on spruces, 210 healthy 4-year-old P. koraiensis Nakai clonal seedlings were used in this study. These seedlings, which were 31–50 cm in height and 5.9–8.4 mm in diameter at the base of the plants, were grown outdoors in 20 cm pots containing a mixture of turfy soil and perlite (v:v = 4:1) in campus of Chinese Academy of Forestry in Beijing. The seedlings were well watered throughout the experiments.
Five O. bicolor strains, TS, BH, LWQ, QH, and MX were inoculated into the spruces. These strains were originally isolated from three Ips bark beetles (Ips typographus Linnaeus from northeast and northwest China, Ips nitidus Eggers from Qinghai-Tibet plateau, and Ips hauseri Reitter from northwest China) (Table 1) and their breeding galleries in infested spruces during the period from 2016 to 2018. Before inoculation, fungi were incubated at 25°C in darkness for 7 days.
Table 1. The information of O. bicolor isolate.
Treatments
On July 15, 2020, the seedlings were randomly assigned to seven different treatment groups with each treatment group including 30 replicates. Among them, five treatment groups were inoculated with TS, BH, LWQ, QH, and MX each, one was inoculated with 2% MEA without fungus, and another was a healthy control. Stems were drilled with a sterile 5 mm cork borer into the surface of xylem to create a small hole, and a 5 mm mycelium plug was placed in the hole using a sterile toothpick. The wound was covered with bark and wrapped with parafilm and tape to protect the wound against contamination and drying.
Analysis of host monoterpenoids
At 4, 30, and 90 days post-inoculation (dpi), 10 seedlings per treatment were used to evaluate host response. The crown and stem of the seedlings were enclosed using a polyethylene (PE) film (48 cm × 60 cm; Reynolds Consumer Products, Lake Forest, IL, United States). Collection of volatile monoterpenoids from P. koraiensis was performed using the closed circulation dynamic headspace sampling system of Porapak-Q absorbent (i.e., 200 mg, 50–80 mesh, 0.6 mm diameter, 160 mm long glass tube, Merck KGaA, Darmstadt, Germany) through mini vacuum pumps (Atmospheric Sampling Instrument, QC-1S, Beijing Labor Protection Institute, China) at an airflow rate of 500 mL/min for 1 h. Sampling tubes containing volatiles trapped on the Porapak-Q absorbent were then sealed at both ends with aluminum foil and placed inside a warm box filled with dry ice. The host tree headspace volatiles were extracted from the sampling tubes with 2 ml of HPLC-grade n-hexane. Prior to gas chromatography (GC), each sample was concentrated to 50 μl with a mild nitrogen stream and stored at −20°C.
Sample analyses were performed using a GC equipped with a flame ionization detector (GC-FID) (Agilent Technologies, Palo Alto, CA, United States) and an automatic sampler for liquid sample injections to identify and quantify host monoterpenes. For each sample run, 1 μl of extract was injected into an HP-5 column (Agilent Technologies, 30 m × 0.25 mm i.d. × 0.25 μm film thickness). Analysis conditions were as follows: H2 as carrier gas at 15 psi column head pressure; flame ionization detector temperature 270°C, and injector temperature 250°C; the oven temperature program started at 45°C (kept isothermal, 1 min), and increased linearly to 105°C at 2°C min–1 (kept isothermal, 1 min), and then to 250°C at 15°C min–1 (kept isothermal, 30 min). Data acquisition and subsequent processing were performed using the Agilent ChemStation GC Systems software. Absolute amounts of monoterpenes were quantified according to the retention time using the following standards: (-)-α-pinene, camphene, β-pinene, myrcene, 3-carene, α-phellandrene, (-)-limonene, and γ-terpinene, and the monoterpene concentrations in the samples were calculated according to external standards.
Analysis of fungal pathogenicity
At 30 and 90 dpi, 10 seedlings per treatment were used to evaluate fungal infection. All barks were removed from the lesion area around the inoculation point with a knife, and the length and width of each necrotic lesion were measured. Three random replicates of every treatment were selected, to take approximately 1 cm2 samples from the edges of each necrotic lesion for re-isolation of the fungi.
Analysis of host transcriptome
The monoterpenoid release peak reached at 4 dpi and then distinctly decreased at later two sampling time points (30 and 90 dpi). Thus, the early molecular changes of P. koraiensis after infection by O. bicolor were analyzed. At 4 dpi, seedlings inoculated with the most and least virulent strains were used for transcriptome sequencing by Illumina. Phloem tissue was collected from a 1 × 4 cm square at the disease/health junction, wrapped in aluminum foil, and placed in dry ice. These samples were stored at −80°C in the laboratory prior to transcriptomic analysis.
Ribonucleic acid from each treatment group with three biological replicates was extracted using the RNAprep Pure Plant Kit (Polysaccharides and Polyphenolics-rich, TIANGEN, Beijing, China) by following the manufacturer’s protocol. RNA concentration and purity were measured using NanoDrop 2000 (Thermo Fisher Scientific, Wilmington, DE, United States). RNA integrity was assessed using the RNA Nano 6000 Assay Kit on the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, United States). The obtained RNA was stored at −80°C.
The sequencing library was constructed using 1 μg RNA from each sample, which was generated using the NEBNext Ultra RN Library Prep Kit from Illumina (New England Biolabs, Ipswich, MA, United States) by following the manufacturer’s protocol. Briefly, mRNA was enriched using magnetic oligo (dT) beads. Fragmentation was carried out using divalent cations under elevated temperatures in NEBNext First Strand Synthesis Reaction Buffer (5X). The short mRNA fragments were used as templates, and the first-strand cDNA was synthesized using a random hexamer primer and M-MuLV reverse transcriptase. Buffer solution, dNTPs, RNase H, and DNA polymerase I were added to synthesize the second-strand cDNA. The remaining overhangs were converted into blunt ends via exonuclease/polymerase activity. After adenylation of the 3′ ends of the DNA fragments, NEBNext Adaptor with a hairpin loop structure was ligated for hybridization. The AMPure XP system (Beckman Coulter, Beverly, MA, United States) was then used for fragment size selection. Suitable fragments were used as templates for PCR amplification to generate a final cDNA library. Finally, the cDNA library was sequenced using the Illumina HiSeq™ 2000 system (Biomarker Technologies Co., Ltd., Beijing, China). The entire set of raw reads was submitted to the Gene Expression Omnibus (GEO) at NCBI under the accession number PRJNA835255.
To obtain high-quality reads to ensure the accuracy of the subsequent analyses, reads containing adapters and poly-N, and low-quality reads (reads containing more than 50% of bases with a Q-value of ≤ 10%) were removed from the sequencing results. Cleaned reads were then mapped to the reference genome of the Norway spruce using HISAT2 software (Kim et al., 2015). Mapped reads were assembled and quantified using the StringTie software (Pertea et al., 2015). The fragments per kilobase of transcript per million fragments mapped (FPKM) measure was used to estimate gene expression levels to determine significant changes in gene expression under different treatments. Hierarchical cluster analysis was used to evaluate the consistency of the sequencing data.
Gene functions were annotated based on the KEGG Ortholog database (KO) and Gene Ontology (GO) databases. Differential expression analysis was performed in the DESeq2 software (Love et al., 2014) for pairwise comparisons using a model based on a negative binomial distribution. The Benjamini–Hochberg method was used to control for the false discovery rate (FDR), with a p-value < 0.01. KOBAS software (Mao et al., 2005) was used to test the statistical enrichment of differentially expressed genes (DEGs) in KEGG pathways.
Quantitative real-time PCR
The CFX96TM Real-Time system (Bio-Rad, Hercules, CA, United States) and SYBR Green FP205 Kit (TIANGEN, China) were used to quantify the expression levels of eight genes in the most and least virulent strain treatment and control. Primer 5.0 (Premier Biosoft) was used to design primer sequences (Supplementary Table 1), and the primers were synthesized by Sangon Biotech (Shanghai) Co., Ltd. Each sample contained three technical replicates. The relative expression levels were calculated using the 2–ΔΔ Ct method and normalized using the actin gene (NCBI accession number: AAF03692) as the internal reference (OuYang et al., 2015).
Data analysis
The area of the lesion was calculated by multiplying the length by the width (mm2) (Rajtar et al., 2021). The monoterpenoid release rate was analyzed based on the difference (Δ) between the treated and healthy groups. Transcriptomic analysis was performed using the BMK cloud platform.1 The lesion area generated by different strains in the same time period and monoterpenoid concentrations from different treatments at the same time or from same treatments at different times were analyzed by one-way analysis of variance (ANOVA). The relationship between the pathogenicity and monoterpenoid changes and the relationship between RNA-seq and quantitative reverse transcription PCR (RT-qPCR) were analyzed by Pearson correlation. IBM SPSS Statistics 19 was used to perform ANOVA (LSD test, α = 0.05) and Pearson correlation analysis.
Results
Lesion area and fungal re-isolation
During the experiment, the inoculated P. koraiensis seedlings had no visible disease symptoms or wilt on the whole seedlings. After removal of the outer bark, all treatment seedlings inoculated with O. bicolor strains showed evident lesions on the phloem tissues of their stems (Figure 1).
Figure 1. Disease development of P. koraiensis stems infected by O. bicolor at 30 and 90 dpi.
As shown in Table 2 and Figure 1, at 30 dpi, the lesion area caused by TS was the largest among the five strains treatments (371.00 ± 37.45 mm2), followed by BH, QH, and MX. The lesion area caused by LWQ was the smallest (143.90 ± 30.94 mm2). ANOVA results showed that the lesion area of O. bicolor inoculation was significantly larger than that of the control (85.63 ± 10.00 mm2), and the lesion area caused by TS was significantly larger than that of other strains. At 90 dpi, the lesion area caused by TS was consistently the largest (432.72 ± 37.38 mm2) among the five strains treatments (Table 2; Figure 1), followed by BH, QH, and MX. The lesion area caused by LWQ was the smallest (160.53 ± 21.20 mm2). ANOVA results showed that the lesion area of O. bicolor inoculation was significantly larger than that of the control (93.06 ± 13.26 mm2), and the lesion area caused by LWQ was significantly smaller than that of other strains. The above results indicated that the pathogenicity of the five strains was differentiated, among which TS was the strongest, BH was the second, QH and MX were the third, and LWQ was the least. The mean re-isolation rate of the inoculated fungi from each of the three random lesion samples was 88.7%, confirming that necrotic lesions were caused by O. bicolor inoculation.
Table 2. The area of necrotic lesions caused by different O. bicolor isolate.
Analysis of host monoterpenoids
Due to the unavailability of α-phellandrene and γ-terpinene for detection, we analyzed and quantified the release rates of six different monoterpenes in whole seedlings, including (-)-α-pinene, camphene, β-pinene, myrcene, 3-carene, and (-)-limonene. From 4 to 90 dpi, the monoterpenoid release rate of spruce was the highest at 4 dpi and decreased gradually thereafter as shown in Figure 2. At 4 dpi, the release rates of monoterpenoids were significantly different between various strain treatments, particularly between the strains that were highly pathogenic strains and those that were the lowly pathogenic. For example, the six monoterpenes release rate in TS treatment were significantly higher than those in LWQ and control treatment (Figure 2; Supplementary Table 2). After 30 dpi, the release rates of the six monoterpenes among the inoculation treatments varied marginally. For example, at 90 dpi, the release rates of (-)-α-pinene, β-pinene, and myrcene were almost the same under the treatment of five strains, and same as the control as well (Figure 2; Supplementary Table 2). In addition, the release rates of different monoterpenes during the three time periods of the control treatments showed no significant differences (Figure 2; Supplementary Table 2). These results indicated that the changes in monoterpenoid release rates were caused by fungal infection, and not by wounding.
Figure 2. Release rates of six monoterpenes in different treatments (μg/h, mean ± SE, N = 3).
The release rate of monoterpenoids among five strain treatments showed a similar trend with the virulence gradients of five strains. In order to further analyze the relationship between monoterpenoids release rate and fungal infection ability, we conducted correlation analysis for the pathogenicity and monoterpenoid release rate at 30 dpi. As shown in Table 3, (-)-α-pinene, β-pinene, 3-carene, and (-)-limonene release rate were highly correlated with fungi pathogenicity (Pearson correlation analysis, R-value ≥ 0.919, p ≤ 0.01). Therefore, the amount of monoterpenoids released can be used to indicate the virulence of fungal strains.
Table 3. Correlation analysis for the pathogenicity and monoterpenoid changes at 30 dpi.
Global review of transcriptome sequencing data
Based on the results of lesion area and monoterpenoid release analyses, TS, LWQ, and control samples were selected for transcriptomic analysis, of which TS was a highly virulent strain and LWQ was a weakly virulent one. A total of 56.97 Gb of clean reads with Q30 > 93.57% were obtained from nine RNA-seq samples. The clean reads of each sample were aligned with the Picea abies reference genome (GCA_900067695.1), with efficiencies ranging from 79.23 to 80.55%, indicating high-quality raw data (Table 4). In total, 37,965 functionally annotated genes were identified (Supplementary Figure 1). A total of 24,433 new genes were discovered in all samples, of which 14,227 were sequenced and annotated to supplement and improve the original genome annotation (Supplementary Table 3). Pearson clustering analysis of the gene expression results showed that all the biological repeats in TS, LWQ, and control were clustered together (Figure 3A). The Pearson coefficient of all three biological replications in all treatments were ≥ 0.987, indicating that our sequencing data had high reliability (Figure 3A).
Table 4. Overview of the transcriptome sequencing dataset and quality check.
Figure 3. Global evaluation of transcriptome sequencing data of spruce. (A) Hierarchical clustering analysis of gene expression shown the correlation among samples. (B) Venn diagram of DEGs in compared between TS, LWQ, and control treatments. DEGs were selected using fold change ≥ 4, FDR correction < 0.01. (C) Volcano plot of all detected genes. Red represents upregulation; green represents downregulation; black represents non-differentially expression. (D) TFs differentially expressed under O. bicolor. X-axis represents the number of DEGs, and Y-axis represents the number names of transcription factor family.
Differentially expressed genes in spruce infected by Ophiostoma bicolor
To explore the difference in the gene expression of P. koraiensis in response to various O. bicolor strains with different virulence, the transcriptomes of the seedlings in the TS, LWQ, and control treatment groups were compared. DEG parameters were set as follows: fold change ≥ 4 and FDR correction < 0.01. A total of 8,496 DEGs were detected in the two comparisons (TS vs. control, and LWQ vs. control), of which 5,300 DEGs were detected in the comparison of TS vs. control and LWQ vs. control, simultaneously (Figure 3B). In addition, 2,388 and 808 DEGs were detected in the comparison of TS vs. control and LWQ vs. control, respectively (Figure 3B). The number of upregulated genes accounted for more than 81.63% of the total DEGs, indicating that O. bicolor inoculation increased the overall gene expression of spruce at the genomic level (Figure 3C; Table 5, chi-square test, p-value < 0.0001). Annotation analysis suggested that the DEGs comprised many transcription factors (TFs), including AP2, RLK, NAC, C2H2, SNF2, PHD, and WRKY (Figure 3D; Supplementary Table 4). Compared to control, in TS- and LWQ-inoculated spruces, 281 and 353 TF-encoding genes were differentially expressed, among which 251 and 271 genes were upregulated (Table 5).
Table 5. Statistical of differentially expressed gene (DEG).
Gene Ontology enrichment analysis was used to assign the DEGs to 46 GO terms. In biological processes, DEGs involved in cellular, metabolite, and signal-organism processes were enriched (Supplementary Figure 2). In terms of molecular function, the DEGs involved in regulated binding, catalytic activity, and transporter activity were enriched (Supplementary Figure 2). In the cellular component category, DEGs involved in cells, organelles, and membranes were enriched (Supplementary Figure 2). These results suggested that binding activity and high enzymatic activity were involved in the defense response of P. koraiensis.
Kyoto encyclopedia of genes and genomes (KEGG) analysis enriched these DEGs in 13 metabolic pathways in all treatments (q < 0.05) (Figure 4; Supplementary Table 5; Kanehisa et al., 2021). Some of the P. koraiensis DEGs were mutually enriched in the TS and LWQ treatments, including RNA transport, starch and sucrose metabolism, aminoacyl-tRNA biosynthesis, mRNA surveillance pathway, ABC transporters, nucleotide excision repair, and base excision repair. Nevertheless, some DEGs were enriched only in one of the treatments. For example, DEGs involved in protein processing in the endoplasmic reticulum and other glycan degradation were enriched in the TS treatment, whereas DEGs involved in homologous recombination, ribosome biogenesis in eukaryotes, mismatch repair, and RNA degradation were enriched in the LWQ treatment (Figure 4 and Supplementary Table 5). DEGs involved in the mRNA surveillance and base excision repair pathways were upregulated, and DEGs in the remaining pathways were both upregulated and downregulated.
Figure 4. Top 20 KEGG pathway analysis of DEGs.
Differentially expressed genes involved in terpenoid biosynthesis
A total of 43 DEGs (24 upregulated and 19 downregulated) involved in terpenoid biosynthesis and metabolism were detected in TS vs. control and LWQ vs. control together (Supplementary Table 6). In both comparisons, from the viewpoint of terpenoid backbone biosynthesis pathway, key enzymes including acetyl-CoA C-acetyltransferase (ACAT) and hydroxymethylglutaryl-CoA reductase (HMGCR) involved in the isopentenyl pyrophosphate (IPP) synthesis, 4-hydroxy-3-methylbut-2-en-1-yl diphosphate reductase involved in the synthesis of dimethylallyl pyrophosphate (DMAPP), and farnesyl diphosphate synthase (FDPS) involved in the synthesis of monoterpenoid precursors geranyl pyrophosphate (GPP) from DMAPP, were activated, whereas isoprene synthase (ispS) involved in the synthesis of isoprene from DMAPP, were downregulated (Figure 5A; Supplementary Table 6). Enzyme genes in the monoterpene synthesis pathway were nearly unchanged, while ent-copalyl diphosphate synthase (ent-CPS) and miltiradiene synthase/copalyl diphosphate synthase (MDS) genes involved in the diterpenoid synthesis pathway were upregulated (Supplementary Table 6).
Figure 5. Transcriptional profiling of differentially expressed genes (DEGs) associated with Terpenoid backbone biosynthesis, flavonoid biosynthesis and phenylpropanoid biosynthesis pathway. The log2FC (fold change) values for the DEGs were used for each treatment (TS vs. Control and LWQ vs. Control). The progression of the color scale from blue to red represents an increase in the log2FC values. Sequences of all new genes were in Supplementary Table 7. (A) The DEGs involved in terpenoid backbone biosynthesis pathway. (B) The DEGs involved in flavonoid biosynthesis pathway. (C) The DEGs involved in phenylpropanoid biosynthesis pathway.
The DEGs related to terpenoid backbone biosynthesis were not identical in these two comparisons (Supplementary Table 6). For example, in LWQ treatment, the MA_1558499g0010 gene related to HMGCR synthesis was downregulated and MA_70145g0010 gene related to ispS synthesis was upregulated, while these two genes were not differentially expressed in TS treatment. In TS treatment the ispS synthesis related genes, Picea_abies_newGene_3028, Picea_abies_newGene_32817, Picea_abies_newGene_41386, MA_10435955g0010, and MA_17608g0020 genes were downregulated, but they were not differentially expressed in LWQ treatment.
Differentially expressed genes involved in plant-pathogen interaction and mitogen-activated protein kinase signaling pathway
Of the 241 DEGs in the plant-pathogen interaction pathway, 212 DEGs were upregulated and involved in all changed nodes except pathogenesis-related protein 1 (PR1) node (Supplementary Figure 3; Supplementary Table 6). In TS vs. control, TS treatment activated the nodes of enhanced disease susceptibility 1 protein (EDS1) and elongation factor Tu (elf18). In LWQ vs. control, LWQ treatment activated the nodes of the molecular chaperone HtpG (HSP90A) and basic helix-loop-helix TF Upa20 (UPA20). Notably, the expression of calcium-binding protein (CML) nodes was activated by TS treatment and inhibited by LWQ treatment (Supplementary Table 6).
In mitogen-activated protein kinase (MAPK) signaling pathway, most of the 109 DEGs at 26 nodes were upregulated (Supplementary Figure 4; Supplementary Table 6). In addition to PR1, the catalase 1 (CAT1) gene also showed downregulated expression. Compared with TS vs. control, LWQ vs. control showed more downregulated DEGs in the MAPK signaling pathway, including basic endochitinase B (ChiB), ethylene-responsive TF 1 (ERF1), serine/threonine-protein kinase (OXI1), protein phosphatase 2C (PP2C), and vegetative storage protein 2 (VSP2).
Differentially expressed genes involved in flavonoid and phenylpropanoid biosynthesis
In samples treated with the two O. bicolor strains, enzymes in the flavonoid biosynthesis pathway displayed various expression patterns. Chalcone isomerase (CHI) showed upregulated expression, trans-cinnamate 4-monooxygenase (CYP73A) showed downregulated expression, and flavonoid 3′-monooxygenase (CYP75B1) and bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR) were both upregulated and downregulated (Figure 5B; Supplementary Table 6). In addition, in the samples inoculated with the low-virulence strain LWQ, the DEGs that encoded chalcone synthase (CHS), flavonoid 3′,5′-hydroxylase (CYP75A), naringenin 3-dioxygenase (F3H), flavonol synthase (FLS), anthocyanidin synthase (ANS), leucoanthocyanidin reductase (LAR), and anthocyanidin reductase (ANR) were downregulated, whereas these DEGs were not detected in the samples that were inoculated with high-virulence strain TS (Supplementary Table 6).
The results also showed that O. bicolor inoculation altered the expression of most genes encoding enzymes related to phenylpropanoid biosynthesis in P. koraiensis (Figure 5C). For example, the expression of genes encoding 4-coumarate CoA ligase (4CL), caffeic acid 3-O-methyltransferase/acetylserotonin O-methyltransferase (COMT), and cinnamoyl-CoA reductase (CCR) was upregulated. Genes encoding cinnamyl-CoA dehydrogenase (CAD), peroxidase (POD), and CYP73A were also downregulated. However, the genes encoding scopoletin glucosyltransferase (TOGT1) and beta-glucosidase (bglx) were both upregulated and downregulated. While some DEGs of P. koraiensis were enriched in the TS and LWQ treatment groups, some were enriched in either TS or LWQ treatment groups. For example, genes encoding coniferyl-aldehyde dehydrogenase (REF1) were solely enriched and downregulated with TS treatment, whereas genes encoding caffeoylshikimate esterase (CSE) were only enriched and upregulated with LWQ treatment (Supplementary Table 6).
RT-qPCR
This study verified the expression of eight genes including calcium-dependent protein kinase (CDPK), ent-CPS, ispS, FLS2, LAR1, (E)-8-carboxylinalool synthase (CYP76F14), (3S,6E)-nerolidol synthase (NES1), and LAR2. Although the expression level of each gene in the transcriptome was higher than that of each gene in the RT-qPCR, both expression levels possessed consistent upregulated or downregulation as shown in Figure 6.
Figure 6. Validation of RNA-seq data by RT-qPCR. R was the correlation coefficient between the FPKM value of RNA-seq and the relative expression level of RT-qPCR, and p was the significance between the FPKM value of RNA-seq and the relative expression level of RT-qPCR. p < 0.01 means a significant correlation at the 0.01 level.
Discussion
With global climate change, forest disasters caused by parasitic and semi-parasitic fungi have become more severe, which was also the case with the beetle-vectored fungi (Kim et al., 2021; Klesse et al., 2021; Morrison et al., 2021; Sitz et al., 2021; Li et al., 2022). However, compared with the former, the diseases caused by beetle-vectored fungi, particularly in gymnosperms, are not well studied. This study combined pathophysiological and transcriptomic analyses to investigate the defense response of P. koraiensis to O. bicolor infection. The results showed that: (1) different O. bicolor strains showed pathogenicity differentiation in P. koraiensis; (2) O. bicolor can induce the release of monoterpenoids in spruce in a short period of time; the release rates are correlated with the fungal virulence. Subsequently, the release rates gradually decrease to a normal level comparable to the control; (3) O. bicolor can significantly change the overall gene expression of P. koraiensis in early infection, and induce the upregulated expression of terpenoid backbone biosynthesis, plant-pathogen interaction, and MAPK signaling pathway.
Some ophiostomatoid fungi associated with I. typographus were pathogenic to spruce, and could induce the host to produce terpenoids, phenols, and other defense substances that were related to disease resistance (Solheim and Safranyik, 1997; Yamaoka et al., 2000; Sallé et al., 2005; Repe et al., 2015). In inoculation experiments by ophiostomatoid fungi, fungal virulence to hosts was not necessarily fatal (Kolb et al., 2019; Stewart et al., 2020). In this study, all five strains of O. bicolor could cause lesions around the inoculation site in P. koraiensis but could not cause wilting or death of the host tree (Figure 1). The elicited lesion areas around the inoculation sites decreased in the order of TS, BH, QH, MX, and LWQ (Table 2). Devkota and Eckhardt (2018) showed that there were differences in host defense responses to different fungal strains inoculation. These results indicated that O. bicolor was a fungus that was weakly virulent to P. koraiensis, and there was pathogenicity differentiation among different strains.
Trees can generate and transfer terpenes, such as monoterpenes, diterpenes, and sesquiterpenes, which are toxic to beetles and fungi; these terpenes could prevent further attacks from beetles and fungi. Raffa and Berryman (1983) reported that lodgepole pines could resist mountain pine beetles by producing high concentrations of two monoterpenes, α-pinene, and limonene, which was further verified in the field by Erbilgin et al. (2017). Meanwhile, the resistance induced by the bark beetle in the host also inhibits the infection of the fungus (Lombardero et al., 2019). The results of this study showed that in addition to the pathogenicity differentiation of different fungi strains, the spruce monoterpenoids synthesis that was induced by O. bicolor, also showed differentiation (Supplementary Table 2). Furthermore, the release rate of monoterpenoids from the host increased at first and peaked at 4 dpi. From 4 to 90 dpi, the release rate of monoterpenoids showed a gradual downward trend (Figure 2). The release rates of monoterpenoids were correlated with the fungal virulence at 30 dpi and almost did not differ with each other under all treatments at 90 dpi. This phenomenon is consistent with the findings of a previous study, which showed that after MeJA was applied to P. abies, the host displayed an induction of terpenes, with concentrations peaking around 16 days after treatment and returning to near-normal levels within 32 days after treatment (Erbilgin et al., 2006). Mageroy et al. (2020a) further verified this phenomenon in Norway spruce, which might be the result of defense priming. Previous studies have shown that host trees inoculated with ophiostomatoid fungi, i.e., E. polonica and Endoconidiophora fujiensis, could generate a rapid and long-term release of monoterpenoids, which possibly deplete the host defense ultimately (Arango-Velez et al., 2018; Fang et al., 2020; Ott et al., 2021). In this study, the release of monoterpenoids did not drastically increase, which was likely due to the different fungal species and inoculation scenarios. The dominant species of ophiostomatoid fungi, which is less virulent, was selected to evaluate host response, whereas in the previous study, the pioneer species, which is more virulent, was selected. This result is in accordance with the conclusion that more virulent fungi induce large changes in host defense terpenoids (Cale et al., 2019; Fang et al., 2020). Therefore, we speculate that O. bicolor stimulates the host to produce defense priming; however, subsequent experiments using MeJA treatment as a positive control of priming defense and challenging stimulus, e.g., E. polonica inoculation or beetle infestation are needed for further verification.
To reveal the regulation of gene expression related to disease resistance and immunity in P. koraiensis, transcriptomic analysis was performed. The results showed that in the early stage of O. bicolor infection of spruce stems (4 dpi), most of the DEGs were upregulated (Figure 3), and these DEGs were mainly related to repair, binding, and signal transduction processes. In the transcriptomic analysis of Grosmannia clavigera hyphae treated with terpenoids for 12 h, it was observed that genes encoding DNA repair, recombination, stabilization, and replication proteins were induced (DiGuistini et al., 2011). This study expands our knowledge of the effects of ophiostomatoid fungi on hosts. Early defense signals in plants often lead to the activation of downstream TF genes, thereby enhancing the expression of defense-related genes (Eulgem and Somssich, 2007). WRKYs play important roles in plant immunity through the ET, JA, and SA pathways in response to various biotic stressors. Several WRKY TFs (mostly upregulated) were involved in the spruce response to O. bicolor. We focused on previously reported pathways involved in plant immune responses, such as plant-pathogen interactions, MAPK signaling, flavonoid biosynthesis, phenylpropane biosynthesis, and terpene biosynthesis and metabolism, to elucidate the molecular mechanisms underlying host defense responses (Gao et al., 2018; Su et al., 2018; Song et al., 2019).
In plants, the biosynthetic precursor of terpenoids can be synthesized via two pathways: the MVA and MEP pathways (Roberts, 2007). In this study, numerous DEGs involved in terpenoid biosynthesis by the MVA and MEP pathways were upregulated (Figure 5A; Supplementary Table 6). On the one hand, the genes of FDPS, which are involved in the synthesis of various monoterpene precursors, and GPP, which comes from DMAPP, were upregulated. On the other hand, the genes of ispS involved in the synthesis of isoprene from DMAPP, were downregulated and thus enhanced the monoterpenoids precursor biosynthesis (Figure 5A). These results suggest that O. bicolor inoculation in P. koraiensis promoted the biosynthesis of host monoterpene precursor, thereby increasing the release of six monoterpenoids. The expression of enzyme genes in the monoterpene synthesis pathway remained nearly unchanged, which is expected and consistent with the findings of previous research (Mageroy et al., 2020b). Moreover, the differences of DEGs related to ispS and HMGCR in TS and LWQ treatments may be the reason for the differences in the release rates of monoterpenoids in the two strain treatments.
Ophiostoma bicolor significantly altered the expression of genes related to plant-pathogen interactions pathways, involving 241 DEGs at 25 nodes. Upregulated genes accounted for 87.97% and were mainly involved in biological processes such as hypersensitive response (HR), cell wall reinforcement, stomatal closure, defense-related gene induction, phytoalexin accumulation, miRNA production, programmed cell death, suppression of plant HR, and defense response (Supplementary Figure 3). These processes are closely related to plant defense responses (Yang et al., 1997). Binding of plant PR proteins to pathogen effector proteins activates ion channels, oxidative bursts, and other signal transductions (Dixon et al., 1994). Ca2+ activates the plant early surveillance system to prevent microbial infection in plant defense signaling (Nürnberger and Scheel, 2001). In the present study, the host treated with the two strains (TS and LWQ treatments) showed differential expression pattern of CML, and LWQ treatment showed less downstream calcium-dependent process in defense response.
The upregulated genes in the MAPK signaling pathway accounted for 68.81% (75 of 109) of all upregulated genes and were mainly involved in biological processes such as cell death, accumulation of reactive oxygen species in plant defense, defense responses, stress adaptation, and wounding. The downregulated genes were mainly related to defense responses against pathogens and wounding (Supplementary Figure 4). MAPK cascades play an important role in plant defense against pathogens (Pitzschke et al., 2009). Its activation is one of the earliest responses for sensing pathogen-associated molecular patterns (PAMPs) in plants (Meng and Zhang, 2013). Some DEGs of LWQ treatment were downregulated and enriched in defense response against pathogens, cell death process and active oxygen production process, while these genes were not differentially expressed in TS treatment (Supplementary Table 6), which might indicate weaker induction of host defense response by the weakly virulent strain. These results suggest that P. koraiensis may use different MAPKs to deliver attack signals, and O. bicolor activated the earliest defense via the MAPK cascade, thus inhibiting further infection.
Through the phenylpropanoid pathway, many kinds of lignin (p-hydroxyphenyl lignin, guaiacyl lignin, 5-hydroxy-guaiacly lignin, and syringyl lignin) can be produced by plants. Lignin can help plants resist biotic and abiotic stresses by regulating secondary cell wall development and stomata (Chen et al., 2019). The synthesis process can be roughly divided into two steps: firstly, the lignin monomer is synthesized under catalysis by a series of enzymes including the 4CL and CAD, and secondly the lignin monomer is polymerized into bioactive lignin by a series of chemical reactions catalyzed by POD and other enzymes (Boerjan et al., 2003). Inhibition of CAD activity in vascular plants changed lignin content. For example, a decrease in CAD activity in sorghum results in a decrease in the total amount of lignin (Palmer et al., 2008). In the present study, the key enzymes involved in lignin monomer synthesis, including 4CL, COMT, and CCR, were activated (Supplementary Table 6). However, CAD required for the final step of lignin monomer biosynthesis was inhibited (Supplementary Table 6). In addition, the gene expression of POD was downregulated, which is consistent with the findings of previous studies on poplar responses to canker pathogens (Li et al., 2019). Although our experiment did not measure the change in lignin content in spruce, it was known from the literature that the content of lignin was reduced in hosts infected with ophiostomatoid fungi (Villari et al., 2012). As a result, O. bicolor may inhibit the synthesis and metabolism of lignin in the spruce stem.
The flavonoid biosynthesis pathway produces a variety of phenolic compounds such as procyanidin, catechin, gallocatechin, and epicatechin. When pathogens invade spruce, fluorescent inclusion bodies containing phenolic compounds appear in the phloem parenchyma cells of the bark, suggesting that phenolic compounds may play a key role in defense against herbivores and pathogens (Franceschi et al., 2005). After an E. polonica infection in Norway spruce, the content of flavane-3-ols, catechin, and gallocatechin in the bark increased (Hammerbacher et al., 2014). The results of this study showed that most DEGs (22 of 28) associated with the flavonoid biosynthesis pathway were downregulated at 4 dpi (Supplementary Table 6), and O. bicolor infection inhibited phenolics-related gene expression in spruce, which differed from the expression pattern of the virulent pioneer E. polonica. The downregulated expression pattern in this study was consistent with the transcriptomic results of Norway spruce infected with Chrysomyxa rhododendri at 4 dpi (Trujillo-Moya et al., 2020). The flavonoid contents of spruces inoculated by O. bicolor were not investigated in this study. Thus, further research is needed to clarify the correlation between the contents and gene expressions of flavonoids.
Conclusion
The ophiostomatoid fungal flora of bark beetles is rich and diverse. This study showed that a dominant fungus of the flora, O. bicolor, is a weakly virulent pathogen of spruce, and various strains showed differentiation in pathogenicity. The release of monoterpenoids from the host is positively correlated with the virulence of the inoculated fungus. In terms of amounts of monoterpenoids released by hosts after inoculation, the pioneer species, e.g., Endoconidiophora spp., can induce rapid and long-term defense responses in the host, while O. bicolor is likely to induce the host defense priming phenomenon. The response of spruce to O. bicolor is a complex process involving multiple biological processes, including plant-pathogen interaction process and MAPK metabolic process. The expression patterns of spruce defense-related genes during infection with different O. bicolor strains were different. However, it was still unclear how the host coordinates these different defense responses against different fungal strains. Further studies are necessary to investigate the underlying mechanism of the host defense response against the pathogen fungi. Bark beetles and their associated fungi usually attack adult trees, thus our inoculations conducted on seedlings have certain limitations. The experiments should attempt to conduct under natural conditions in the future.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
QL and YL conceptualized and designed the research. YL, QZ, and GZ performed the experiments. YL and QZ analyzed the data and wrote the manuscript. YL, QZ, ZW, HW, JZ, and QL participated in the discussion for experimental details. QL and JZ revised the manuscript. QL directed the project. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (Project nos. 32071769 and 31770682).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Ruhua Li of Chinese Academy of Forestry for support in GC-FID metabolite analysis, Fude Wang of Heilongjiang Academy of forestry, and Dongsheng Chen of Chinese Academy of Forestry for their support in logistical assistance. We also thank editage (https://www.editage.cn) for its linguistic assistance during the preparation of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.944336/full#supplementary-material
Supplementary Figure 1 | Functional annotation of the all genes detected in P. koraiensis from various databases. (A) NR database homologous species distribution analysis; (B) GO; (C) eggNOG; (D) KOG; and (E) COG.
Supplementary Figure 2 | Gene Ontology (GO) classification of DEGs.
Supplementary Figure 3 | Plant-pathogen interaction pathway map. Upregulation is highlighted in red; downregulation is highlighted in green; mixed regulation is highlighted in blue.
Supplementary Figure 4 | Mitogen-activated protein kinase (MAPK) signaling pathway map. Upregulation is highlighted in red; downregulation is highlighted in green; mixed regulation is highlighted in blue.
Supplementary Table 1 | Primer sequences used in RT-qPCR.
Supplementary Table 2 | Release rates of major monoterpenes in different treatments (μg/h, mean ± SE, N = 3).
Supplementary Table 3 | Statistics of annotation results for new genes in databases.
Supplementary Table 4 | All the different transcription factors.
Supplementary Table 5 | Statistics of KEGG pathway with significant enrichment.
Supplementary Table 6 | Differential gene expression levels in five KEGG pathways.
Supplementary Table 7 | Sequences of all differentially expressed new genes involved in five KEGG pathways.
Footnotes
References
Arango-Velez, A., Chakraborty, S., Blascyk, K., Phan, M. T., Barsky, J., and El Kayal, W. (2018). Anatomical and chemical responses of eastern white pine (Pinus strobus L.) to blue-stain (Ophiostoma minus) inoculation. Forests 9:690. doi: 10.3390/f9110690
Bentz, B. J., and Six, D. L. (2006). Ergosterol content of fungi associated with Dendroctonus ponderosae and Dendroctonus rufipennis (Coleoptera: Curculionidae, Scolytinae). Annu. Entomol. Soc. Am. 99, 189–194.
Biedermann, P. H., and Vega, F. E. (2020). Ecology and evolution of insect–fungus mutualisms. Annu. Rev. Entomol. 65, 431–455. doi: 10.1146/annurev-ento-011019-024910
Bleiker, K., and Six, D. L. (2009). Competition and coexistence in a multi-partner mutualism: interactions between two fungal symbionts of the mountain pine beetle in beetle-attacked trees. Microb. Ecol. 57, 191–202. doi: 10.1007/s00248-008-9395-6
Boerjan, W., Ralph, J., and Baucher, M. (2003). Lignin biosynthesis. Annu. Rev. Plant Biol. 54, 519–546. doi: 10.1146/annurev.arplant.54.031902.134938
Boone, C. K., Aukema, B. H., Bohlmann, J., Carroll, A. L., and Raffa, K. F. (2011). Efficacy of tree defense physiology varies with bark beetle population density: a basis for positive feedback in eruptive species. Can. J. Forest Res. 41, 1174–1188. doi: 10.1139/x11-041
Brasier, C. (1986). Comparison of pathogencity and cultural characteristics in the EAN and NAN aggressive subgroups of Ophiostoma ulmi. Trans. Br. Mycol. Soc. 87, 1–13. doi: 10.1016/S0007-1536(86)80001-8
Cale, J. A., Klutsch, J. G., Dykstra, C. B., Peters, B., and Erbilgin, N. (2019). Pathophysiological responses of pine defensive metabolites largely lack differences between pine species but vary with eliciting ophiostomatoid fungal species. Tree Physiol. 39, 1121–1135. doi: 10.1093/treephys/tpz012
Celedon, J. M., and Bohlmann, J. (2019). Oleoresin defenses in conifers: chemical diversity, terpene synthases and limitations of oleoresin defense under climate change. New Phytol. 224, 1444–1463. doi: 10.1111/nph.15984
Chakraborty, A., Modlinger, R., Ashraf, M. Z., Synek, J., Schlyter, F., and Roy, A. (2020). Core mycobiome and their ecological relevance in the gut of five Ips bark beetles (Coleoptera: Curculionidae: Scolytinae). Front. Microbiol. 11:568853. doi: 10.3389/fmicb.2020.568853
Chakraborty, A., and Roy, A. (2021). “Microbial Influence on Plant–Insect Interaction,” in Plant-Pest Interactions: From Molecular Mechanisms to Chemical Ecology, eds I. K. Singh and A. Singh (Singapore: Springer), doi: 10.1007/978-981-15-2467-7_14
Chang, R., Duong, T. A., Taerum, S. J., Wingfield, M. J., Zhou, X., and De Beer, Z. W. (2020). Ophiostomatoid fungi associated with mites phoretic on bark beetles in Qinghai, China. IMA Fungus 11, 1–18. doi: 10.1186/s43008-020-00037-9
Chang, R., Duong, T. A., Taerum, S. J., Wingfield, M. J., Zhou, X., Yin, M., et al. (2019). Ophiostomatoid fungi associated with the spruce bark beetle Ips typographus, including 11 new species from China. Persoonia 42, 50–74. doi: 10.3767/persoonia.2019.42.03
Chen, X., Wang, H., Li, X., Ma, K., Zhan, Y., and Zeng, F. (2019). Molecular cloning and functional analysis of 4-Coumarate: coA ligase 4 (4CL-like 1) from Fraxinus mandshurica and its role in abiotic stress tolerance and cell wall synthesis. BMC Plant Biol. 19:231. doi: 10.1186/s12870-019-1812-0
Christiansen, E., and Solheim, H. (1994). Pathogenicity of five species of Ophiostoma fungi to Douglas-fir. Medd. Nor. Inst. Skogforsk 47, 1–12.
Conesa, A., Madrigal, P., Tarazona, S., Gomez-Cabrero, D., Cervera, A., Mcpherson, A., et al. (2016). A survey of best practices for RNA-seq data analysis. Genome Biol. 17, 1–19. doi: 10.1186/s13059-016-0881-8
Dai, L., Li, H., Zheng, J., and Chen, H. (2022). Transcriptome analyses of the Chinese white pine beetle-fungal symbiont Leptographium qinlingensis under terpene stress or growth on host pine sawdust. Symbiosis 86, 17–31. doi: 10.1007/s13199-021-00822-z
Davis, T. S., Stewart, J. E., Mann, A., Bradley, C., and Hofstetter, R. W. (2019). Evidence for multiple ecological roles of Leptographium abietinum, a symbiotic fungus associated with the North American spruce beetle. Fungal Ecol. 38, 62–70. doi: 10.1016/j.funeco.2018.04.008
Davydenko, K., Vasaitis, R., and Menkis, A. (2017). Fungi associated with Ips acuminatus (Coleoptera: Curculionidae) in Ukraine with a special emphasis on pathogenicity of ophiostomatoid species. Eur. J. Entomol. 114, 77–85. doi: 10.14411/eje.2017.011
Devkota, P., and Eckhardt, L. G. (2018). Variation in pathogenicity of different Leptographium terebrantis isolates to Pinus taeda L. Forest Pathol. 48:e12469. doi: 10.1111/efp.12469
DiGuistini, S., Wang, Y., Liao, N. Y., Taylor, G., Tanguay, P., Feau, N., et al. (2011). Genome and transcriptome analyses of the mountain pine beetle-fungal symbiont Grosmannia clavigera, a lodgepole pine pathogen. Proc. Natl. Acad. Sci. U. S. A. 108, 2504–2509. doi: 10.1073/pnas.1011289108
Dixon, R., Harrison, M., and Lamb, C. (1994). Early events in the activation of plant defense responses. Annu. Rev. phytopathol. 32, 479–501. doi: 10.1146/annurev.py.32.090194.002403
Dzurenko, M., and Hulcr, J. (2022). Ambrosia beetles. Curr. Biol. 32, R61–R62. doi: 10.1016/j.cub.2021.11.043
Erbilgin, N., Cale, J. A., Hussain, A., Ishangulyyeva, G., Klutsch, J. G., Najar, A., et al. (2017). Weathering the storm: how lodgepole pine trees survive mountain pine beetle outbreaks. Oecologia 184, 469–478.
Erbilgin, N., Krokene, P., Christiansen, E., Zeneli, G., and Gershenzon, J. (2006). Exogenous application of methyl jasmonate elicits defenses in Norway spruce (Picea abies) and reduces host colonization by the bark beetle Ips typographus. Oecologia 148, 426–436. doi: 10.1007/s00442-006-0394-3
Eulgem, T., and Somssich, I. E. (2007). Networks of WRKY transcription factors in defense signaling. Curr. Opin. Plant Biol. 10, 366–371. doi: 10.1016/j.pbi.2007.04.020
Fang, J., Liu, M., Zhang, S., Liu, F., Zhang, Z., Zhang, Q., et al. (2020). Chemical signal interactions of the bark beetle with fungal symbionts, and host/non-host trees. J. Exp. Bot. 71, 6084–6091. doi: 10.1093/jxb/eraa296
Franceschi, V. R., Krokene, P., Christiansen, E., and Krekling, T. (2005). Anatomical and chemical defenses of conifer bark against bark beetles and other pests. New Phytol. 167, 353–376. doi: 10.1111/j.1469-8137.2005.01436.x
Gao, Y., Wang, W., Zhang, T., Gong, Z., Zhao, H., and Han, G.-Z. (2018). Out of water: the origin and early diversification of plant R-genes. Plant Physiol. 177, 82–89. doi: 10.1104/pp.18.00185
Giordano, L., Garbelotto, M., Nicolotti, G., and Gonthier, P. (2013). Characterization of fungal communities associated with the bark beetle Ips typographus varies depending on detection method, location, and beetle population levels. Mycol. Prog. 12, 127–140. doi: 10.1007/s11557-012-0822-1
Hammerbacher, A., Paetz, C., Wright, L. P., Fischer, T. C., Bohlmann, J., Davis, A. J., et al. (2014). Flavan-3-ols in Norway spruce: biosynthesis, accumulation, and function in response to attack by the bark beetle-associated fungus Ceratocystis polonica. Plant Physiol. 164, 2107–2122. doi: 10.1104/pp.113.232389
Harrington, T., Fraedrich, S., and Aghayeva, D. (2008). Raffaelea lauricola, a new ambrosia beetle symbiont and pathogen on the Lauracea. Mycotaxon 104, 399–404.
Hlásny, T., König, L., Krokene, P., Lindner, M., Montagné-Huck, C., Müller, J., et al. (2021). Bark beetle outbreaks in Europe: state of knowledge and ways forward for management. Curr. For. Rep. 7, 138–165. doi: 10.1007/s40725-021-00142-x
Itoh, H., Tago, K., Hayatsu, M., and Kikuchi, Y. (2018). Detoxifying symbiosis: microbe-mediated detoxification of phytotoxins and pesticides in insects. Nat. Prod. Rep. 35, 434–454. doi: 10.1039/C7NP00051K
Jaakola, L., Määttä-Riihinen, K., Kärenlampi, S., and Hohtola, A. (2004). Activation of flavonoid biosynthesis by solar radiation in bilberry (Vaccinium myrtillus L.) leaves. Planta 218, 721–728. doi: 10.1007/s00425-003-1161-x
Jones, J. D., and Dangl, J. L. (2006). The plant immune system. Nature 444, 323–329. doi: 10.1038/nature05286
Kandasamy, D., Gershenzon, J., Andersson, M. N., and Hammerbacher, A. (2019). Volatile organic compounds influence the interaction of the Eurasian spruce bark beetle (Ips typographus) with its fungal symbionts. ISME J. 13, 1788–1800. doi: 10.1038/s41396-019-0390-3
Kanehisa, M., Furumichi, M., Sato, Y., Ishiguro-Watanabe, M., and Tanabe, M. J. N. A. R. (2021). KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 49, D545–D551. doi: 10.1093/nar/gkaa970
Karonen, M., Ossipov, V., Ossipova, S., Kapari, L., Loponen, J., Matsumura, H., et al. (2006). Effects of elevated carbon dioxide and ozone on foliar proanthocyanidins in Betula platyphylla, Betula ermanii, and Fagus crenata seedlings. J. Chem. Ecol. 32, 1445–1458. doi: 10.1007/s10886-006-9061-8
Keeling, C. I., and Bohlmann, J. (2006). Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens. New Phytol. 170, 657–675. doi: 10.1111/j.1469-8137.2006.01716.x
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Kim, M.-S., Hanna, J. W., Stewart, J. E., Warwell, M. V., Mcdonald, G. I., and Klopfenstein, N. B. (2021). Predicting present and future suitable climate spaces (potential distributions) for an Armillaria root disease pathogen (Armillaria solidipes) and its host, Douglas-fir (Pseudotsuga menziesii), under changing climates. Front. For. Glob. Change 4:740994. doi: 10.3389/ffgc.2021.740994
Kirisits, T., and Konrad, H. (2004). Dutch elm disease in Austria. For. Syst. 13, 81–92. doi: 10.5424/816
Klesse, S., Abegg, M., Hopf, S. E., Gossner, M. M., Rigling, A., and Queloz, V. (2021). Spread and severity of ash dieback in Switzerland–tree characteristics and landscape features explain varying mortality probability. Front. For. Glob. Change 4:645920. doi: 10.3389/ffgc.2021.645920
Kolb, T., Keefover-Ring, K., Burr, S. J., Hofstetter, R., Gaylord, M., and Raffa, K. F. (2019). Drought-mediated changes in tree physiological processes weaken tree defenses to bark beetle attack. J. Chem. Ecol. 45, 888–900. doi: 10.1007/s10886-019-01105-0
Li, P., Liu, W., Zhang, Y., Xing, J., Li, J., Feng, J., et al. (2019). Fungal canker pathogens trigger carbon starvation by inhibiting carbon metabolism in poplar stems. Sci. Rep. 9, 1–14. doi: 10.1038/s41598-019-46635-5
Li, Y., Bateman, C., Skelton, J., Wang, B., Black, A., Huang, Y.-T., et al. (2022). Preinvasion Assessment of Exotic Bark Beetle-Vectored Fungi to Detect Tree-Killing Pathogens. Phytopathology 112, 261–270. doi: 10.1094/PHYTO-01-21-0041-R
Lieutier, F., Brignolas, F., Sauvard, D., Yart, A., Galet, C., Brunet, M., et al. (2003). Intra-and inter-provenance variability in phloem phenols of Picea abies and relationship to a bark beetle-associated fungus. Tree Physiol. 23, 247–256. doi: 10.1093/treephys/23.4.247
Lombardero, M. J., Solla, A., and Ayres, M. P. (2019). Pine defenses against the pitch canker disease are modulated by a native insect newly associated with the invasive fungus. For. Ecol. Manag. 437, 253–262. doi: 10.1016/j.foreco.2019.01.041
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21. doi: 10.1186/s13059-014-0550-8
Lu, Q., Zhang, X., Yang, Z., Henri, M., Yin, D., and Ren, X. (2008). Research progress on the fungi associated with dominus rubiformis. Sci. Silvae Sin. 44, 134–142.
Mageroy, M. H., Christiansen, E., Långström, B., Borg-Karlson, A. K., Solheim, H., Björklund, N., et al. (2020a). Priming of inducible defenses protects Norway spruce against tree-killing bark beetles. Plant Cell Environ. 43, 420–430. doi: 10.1111/pce.13661
Mageroy, M. H., Wilkinson, S. W., Tengs, T., Cross, H., Almvik, M., Pétriacq, P., et al. (2020b). Molecular underpinnings of methyl jasmonate-induced resistance in Norway spruce. Plant Cell Environ. 43, 1827–1843. doi: 10.1111/pce.13774
Mao, X., Cai, T., Olyarchuk, J. G., and Wei, L. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430
Meng, X., and Zhang, S. (2013). MAPK cascades in plant disease resistance signaling. Annu. Rev. Phytopathol. 51, 245–266. doi: 10.1146/annurev-phyto-082712-102314
Morrison, E. W., Kasson, M. T., Heath, J. J., and Garnas, J. R. (2021). Pathogen and endophyte assemblages co-vary with beech bark disease progression, tree decline, and regional climate. bioRxiv [Preprint]. doi: 10.1101/2021.02.24.432696
Nève Repe, A., Groot, M. D., and Jurc, M. (2018). Assemblages of ophiostomatoid fungi vectored by Ips amitinus (Coleoptera: Scolytinae) on norway spruce depend on colonization time, position on the host tree and development stage. Šumarski List 142, 171–178. doi: 10.31298/sl.142.3-4.5
Nibedita, C., and Jolly, B. (2017). Transcriptomics: a successful approach to unravel the molecular mechanism of plant-pathogen interaction in post-genomic era. Res. J. Biotechnol. 8, 90–100.
Nürnberger, T., and Scheel, D. (2001). Signal transmission in the plant immune response. Trends Plant Sci. 6, 372–379. doi: 10.1016/S1360-1385(01)02019-2
Ott, D. S., Davis, T. S., and Mercado, J. E. (2021). Interspecific variation in spruce constitutive and induced defenses in response to a bark beetle–fungal symbiont provides insight into traits associated with resistance. Tree Physiol. 41, 1109–1121. doi: 10.1093/treephys/tpaa170
OuYang, F., Mao, J.-F., Wang, J., Zhang, S., and Li, Y. (2015). Transcriptome analysis reveals that red and blue light regulate growth and phytohormone metabolism in Norway spruce [Picea abies (L.) Karst.]. PLoS One 10:e0127896. doi: 10.1371/journal.pone.0127896
Palmer, N. A., Sattler, S. E., Saathoff, A. J., Funnell, D., Pedersen, J. F., and Sarath, G. (2008). Genetic background impacts soluble and cell wall-bound aromatics in brown midrib mutants of sorghum. Planta 229, 115–127. doi: 10.1007/s00425-008-0814-1
Pan, Y., Zhao, T., Krokene, P., Yu, Z.-F., Qiao, M., Lu, J., et al. (2018). Bark beetle-associated blue-stain fungi increase antioxidant enzyme activities and monoterpene concentrations in Pinus yunnanensis. Front. Plant Sci. 9:1731. doi: 10.3389/fpls.2018.01731
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T.-C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Phillips, M. A., and Croteau, R. B. (1999). Resin-based defenses in conifers. Trends Plant Sci. 4, 184–190. doi: 10.1016/S1360-1385(99)01401-6
Pitzschke, A., Schikora, A., and Hirt, H. (2009). MAPK cascade signalling networks in plant defence. Curr. Opin. Plant Biol. 12, 421–426. doi: 10.1016/j.pbi.2009.06.008
Raffa, K., and Berryman, A. (1983). The role of host plant resistance in the colonization behavior and ecology of bark beetles (Coleoptera: Scolytidae). Ecol. Monogr. 53, 27–49. doi: 10.2307/1942586
Raffa, K. F., Aukema, B. H., Erbilgin, N., Klepzig, K. D., and Wallin, K. F. (2005). Interactions among conifer terpenoids and bark beetles across multiple levels of scale: an attempt to understand links between population patterns and physiological processes. Ecol. Monogr. 39, 79–118.
Raffa, K. F., Powell, E. N., and Townsend, P. A. (2013). Temperature-driven range expansion of an irruptive insect heightened by weakly coevolved plant defenses. Proc. Natl. Acad. Sci. U. S. A. 110, 2193–2198. doi: 10.1073/pnas.1216666110
Rajtar, N. N., Held, B. W., and Blanchette, R. A. (2021). Fungi from Galleries of the Emerald Ash Borer Produce Cankers in Ash Trees. Forests 12:1509. doi: 10.3390/f12111509
Repe, A., Bojović, S., and Jurc, M. (2015). Pathogenicity of ophiostomatoid fungi on Picea abies in Slovenia. For. Pathol. 45, 290–297. doi: 10.1111/efp.12170
Roberts, S. C. (2007). Production and engineering of terpenoids in plant cell culture. Nat. Chem. Biol. 3, 387–395. doi: 10.1038/nchembio.2007.8
Sallé, A., Monclus, R., Yart, A., Garcia, J., Romary, P., and Lieutier, F. (2005). Fungal flora associated with Ips typographus: frequency, virulence, and ability to stimulate the host defence reaction in relation to insect population levels. Can. J. For. Res. 35, 365–373. doi: 10.1139/x04-186
Santini, A., and Faccoli, M. (2015). Dutch elm disease and elm bark beetles: a century of association. iForest 8, 126–134. doi: 10.3832/ifor1231-008
Singh, S., Singh, A., Baweja, V., Roy, A., Chakraborty, A., and Singh, I. K. J. M. (2021). Molecular Rationale of Insect-Microbes Symbiosis—From Insect Behaviour to Mechanism. Microorganisms 9:2422. doi: 10.3390/microorganisms9122422
Sitz, R. A., Luna, E. K., Ibarra Caballero, J., Tisserat, N. A., Cranshaw, W. S., Mckenna, J. R., et al. (2021). Eastern Black Walnut (Juglans nigra L.) Originating From Native Range Varies in Their Response to Inoculation With Geosmithia morbida. Front. For. Glob. Change 4:627911. doi: 10.3389/ffgc.2021.627911
Six, D. L. (2012). Ecological and evolutionary determinants of bark beetle—fungus symbioses. Insects 3, 339–366. doi: 10.3390/insects3010339
Six, D. L., and Wingfield, M. J. (2011). The role of phytopathogenicity in bark beetle–fungus symbioses: a challenge to the classic paradigm. Annu. Rev. Entomol. 56, 255–272. doi: 10.1146/annurev-ento-120709-144839
Solheim, H., and Safranyik, L. (1997). Pathogenicity to Sitka spruce of Ceratocystis rufipenni and Leptographium abietinum blue-stain fungi associated with the spruce beetle. Can. J. For. Res. 27, 1336–1341. doi: 10.1139/x97-096
Song, N., Ma, L., Wang, W., Sun, H., Wang, L., Baldwin, I. T., et al. (2019). An ERF2-like transcription factor regulates production of the defense sesquiterpene capsidiol upon Alternaria alternata infection. J. Exp. Bot. 70, 5895–5908. doi: 10.1093/jxb/erz327
Stewart, J. E., Harris, F. L., Otto, K., and Davis, T. S. (2020). Responses of Engelmann spruce to inoculation with Leptographium abietinum, a symbiotic fungus of the North American spruce beetle. Can. J. For. Res. 50, 465–472. doi: 10.1139/cjfr-2019-0362
Su, X., Lu, G., Guo, H., Zhang, K., Li, X., and Cheng, H. (2018). The dynamic transcriptome and metabolomics profiling in Verticillium dahliae inoculated Arabidopsis thaliana. Sci. Rep. 8, 1–11. doi: 10.1038/s41598-018-33743-x
Sun, J., Lu, M., Gillette, N. E., and Wingfield, M. J. (2013). Red turpentine beetle: innocuous native becomes invasive tree killer in China. Annu. Rev. Entomol. 58, 293–311. doi: 10.1146/annurev-ento-120811-153624
Trapp, S., and Croteau, R. (2001). Defensive resin biosynthesis in conifers. Annu. Rev. Plant Biol. 52, 689–724. doi: 10.1146/annurev.arplant.52.1.689
Trujillo-Moya, C., Ganthaler, A., Stöggl, W., Kranner, I., Schüler, S., Ertl, R., et al. (2020). RNA-Seq and secondary metabolite analyses reveal a putative defence-transcriptome in Norway spruce (Picea abies) against needle bladder rust (Chrysomyxa rhododendri) infection. BMC Genomics 21:336. doi: 10.1186/s12864-020-6587-z
Ullah, A., Klutsch, J. G., and Erbilgin, N. (2021). Production of complementary defense metabolites reflects a co-evolutionary arms race between a host plant and a mutualistic bark beetle-fungal complex. Plant Cell Environ. 44, 3064–3077. doi: 10.1111/pce.14100
Villari, C., Battisti, A., Chakraborty, S., Michelozzi, M., Bonello, P., and Faccoli, M. (2012). Nutritional and pathogenic fungi associated with the pine engraver beetle trigger comparable defenses in Scots pine. Tree Physiol. 32, 867–879. doi: 10.1093/treephys/tps056
Wadke, N., Kandasamy, D., Vogel, H., Lah, L., Wingfield, B. D., Paetz, C., et al. (2016). The bark-beetle-associated fungus, Endoconidiophora polonica, utilizes the phenolic defense compounds of its host as a carbon source. Plant Physiol. 171, 914–931. doi: 10.1104/pp.15.01916
Wang, F., Cale, J. A., and Erbilgin, N. (2020a). Induced defenses of a novel host tree affect the growth and interactions of bark beetle-vectored fungi. Microb. Ecol. 80, 181–190. doi: 10.1007/s00248-020-01490-0
Wang, Y., Gao, Y., Zang, P., and Xu, Y. (2020b). Transcriptome analysis reveals underlying immune response mechanism of fungal (Penicillium oxalicum) disease in Gastrodia elata Bl. f. glauca S. chow (Orchidaceae). BMC Plant Biol. 20:445. doi: 10.1186/s12870-020-02653-4
Wang, Z., Liu, Y., Wang, H., Meng, X., Liu, X., Decock, C., et al. (2020c). Ophiostomatoid fungi associated with Ips subelongatus, including eight new species from northeastern China. IMA Fungus 11, 1–29. doi: 10.1186/s43008-019-0025-3
Wang, H., Wang, T., Liu, Y., Zeng, F., Zhang, H., Decock, C., et al. (2022). Diversity of Ophiostomatoid Fungi Associated with Dendroctonus armandi Infesting Pinus armandii in Western China. J. Fungi 8:214. doi: 10.3390/jof8030214
Wang, Z., Zhou, Q., Zheng, G., Fang, J., Han, F., Zhang, X., et al. (2021). Abundance and diversity of ophiostomatoid fungi associated with the Great Spruce Bark Beetle (Dendroctonus micans) in the Northeastern Qinghai-Tibet Plateau. Front. Microbiol. 12:721395. doi: 10.3389/fmicb.2021.721395
Wingfield, M. J., Garnas, J. R., Hajek, A., Hurley, B. P., De Beer, Z. W., and Taerum, S. J. (2016). Novel and co-evolved associations between insects and microorganisms as drivers of forest pestilence. Biol. Invasions 18, 1045–1056. doi: 10.1007/s10530-016-1084-7
Yamaoka, Y., Takahashi, I., and Iguchi, K. (2000). Virulence of ophiostomatoid fungi associated with the spruce bark beetle Ips typographus f. japonicus in Yezo spruce. J. For. Res. 5, 87–94. doi: 10.1007/BF02762525
Yamaoka, Y., Wingfield, M. J., Ohsawa, M., and Kuroda, Y. (1998). Ophiostomatoid fungi associated with Ips cembrae in Japan and their pathogenicity of Japanese larch. Mycoscience 39, 367–378. doi: 10.1007/BF02460897
Yang, Y., Shah, J., and Klessig, D. F. (1997). Signal perception and transduction in plant defense responses. Genes Dev. 11, 1621–1639.
Zhao, T., Axelsson, K., Krokene, P., and Borg-Karlson, A.-K. (2015). Fungal symbionts of the spruce bark beetle synthesize the beetle aggregation pheromone 2-methyl-3-buten-2-ol. J. Chem. Ecol. 41, 848–852. doi: 10.1007/s10886-015-0617-3
Zhao, T., Ganji, S., Schiebe, C., Bohman, B., Weinstein, P., Krokene, P., et al. (2019a). Convergent evolution of semiochemicals across Kingdoms: bark beetles and their fungal symbionts. ISME J. 13, 1535–1545. doi: 10.1038/s41396-019-0370-7
Zhao, T., Kandasamy, D., Krokene, P., Chen, J., Gershenzon, J., and Hammerbacher, A. (2019b). Fungal associates of the tree-killing bark beetle, Ips typographus, vary in virulence, ability to degrade conifer phenolics and influence bark beetle tunneling behavior. Fungal Ecol. 38, 71–79. doi: 10.1016/j.funeco.2018.06.003
Keywords: transcriptome, ophiostomatoid fungi, monoterpenoids, spruce, Ips
Citation: Liu Y, Zhou Q, Wang Z, Wang H, Zheng G, Zhao J and Lu Q (2022) Pathophysiology and transcriptomic analysis of Picea koraiensis inoculated by bark beetle-vectored fungus Ophiostoma bicolor. Front. Plant Sci. 13:944336. doi: 10.3389/fpls.2022.944336
Received: 15 May 2022; Accepted: 27 June 2022;
Published: 19 July 2022.
Edited by:
Lilin Zhao, Institute of Zoology (CAS), ChinaReviewed by:
Kateryna Davydenko, Swedish University of Agricultural Sciences, SwedenAmit Roy, Czech University of Life Sciences Prague, Czechia
Yonglin Wang, Beijing Forestry University, China
Copyright © 2022 Liu, Zhou, Wang, Wang, Zheng, Zhao and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiaping Zhao, zhaojiaping@caf.ac.cn; Quan Lu, luquan@caf.ac.cn