Transcriptome analyses of the Chinese white pine beetle-fungal symbiont Leptographium qinlingensis under terpene stress or growth on host pine sawdust

Comparative metabarcoding and biodiversity of gut-associated fungal assemblages of Dendroctonus species (Curculionidae: Scolytinae)

The genus Dendroctonus is a Holarctic taxon composed of 21 nominal species; some of these species are well known in the world as disturbance agents of forest ecosystems. Under the bark of the host tree, these insects are involved in complex and dynamic associations with phoretic ectosymbiotic and endosymbiotic communities. Unlike filamentous fungi and bacteria, the ecological role of yeasts in the bark beetle holobiont is poorly understood, though yeasts were the first group to be recorded as microbial symbionts of these beetles. Our aim was characterize and compare the gut fungal assemblages associated to 14 species of Dendroctonus using the internal transcribed spacer 2 (ITS2) region. A total of 615,542 sequences were recovered yielding 248 fungal amplicon sequence variants (ASVs). The fungal diversity was represented by 4 phyla, 16 classes, 34 orders, 54 families, and 71 genera with different relative abundances among Dendroctonus species. The α-diversity consisted of 32 genera of yeasts and 39 genera of filamentous fungi. An analysis of β-diversity indicated differences in the composition of the gut fungal assemblages among bark beetle species, with differences in species and phylogenetic diversity. A common core mycobiome was recognized at the genus level, integrated mainly by Candida present in all bark beetles, Nakazawaea, Cladosporium, Ogataea, and Yamadazyma. The bipartite networks confirmed that these fungal genera showed a strong association between beetle species and dominant fungi, which are key to maintaining the structure and stability of the fungal community. The functional variation in the trophic structure was identified among libraries and species, with pathotroph-saprotroph-symbiotroph represented at the highest frequency, followed by saprotroph-symbiotroph, and saprotroph only. The overall network suggested that yeast and fungal ASVs in the gut of these beetles showed positive and negative associations among them. This study outlines a mycobiome associated with Dendroctonus nutrition and provides a starting point for future in vitro and omics approaches addressing potential ecological functions and interactions among fungal assemblages and beetle hosts.

Establishment of RNA Interference Genetic Transformation System and Functional Analysis of FlbA Gene in Leptographium qinlingensis

Leptographium qinlingensis is a pathogenic fungus of Pinus armandii that is epidemic in the Qinling Mountains. However, an effective gene interference strategy is needed to characterize the pathogenic genes in this fungus on a functional level. Using the RNA silencing vector pSilent-1 as a template, we established an RNA interference genetic transformation system mediated by Agrobacterium tumefaciens GV3101, which is suitable for the gene study for Leptographium qinlingensis by homologous recombination and strain interference system screening. The LqFlbA gene was silenced using the RNA interference approach described above, and the resulting transformants displayed various levels of silencing with a gene silencing effectiveness ranging from 41.8% to 91.4%. The LqFlbA-RNAi mutant displayed altered colony morphology, sluggish mycelium growth, and diminished pathogenicity toward the host P. armandii in comparison to the wild type. The results indicate that this method provides a useful reverse genetic system for studying the gene function of L. qinlingensis, and that LqFlbA plays a crucial role in the growth, development, and pathogenicity of L. qinlingensis.

The effect of inactivation of aldehyde dehydrogenase on pheromone production by a gut bacterium of an invasive bark beetle, Dendroctonus valens

Semiochemical-based management strategies are important for controlling bark beetles, such as invasive Red Turpentine Beetle (Denroctonus valens), the causal agent for mass mortality of pine trees (Pinus spp.) in China. It has been previously shown that the pheromone verbenone regulates the attack density of this beetle in a dose-dependent manner and that the gut bacteria of D. valens are involved in verbenone production. However, molecular functional verification of the role of gut bacteria in the pheromone production of D. valens is still lacking. To better understand the molecular function of gut bacterial verbenone production, we chose a facultative anaerobic gut bacterium (Enterobacter xiangfangensis) of D. valens based on its strong ability to convert cis-verbenol to verbenone, as shown in our previous study, and investigated its transcriptomics in the presence or absence of cis-verbenol under anaerobic conditions (simulating the anoxic environment in the beetle's gut). Based on this transcriptome analysis, aldehyde dehydrogenase (ALDH1) was identified as a putative key gene responsible for verbenone production and was knocked-down by homologous recombination to obtain a mutant E. xiangfangensis strain. Our results show that these mutants had significantly decreased the ability to convert the monoterpene precursor to verbenone compared with the wild-type bacteria, indicating that ALDH1 is primarily responsible for verbenone conversion for this bacterium species. These findings provide further mechanistic evidence of bacterially mediated pheromone production by D. valens, add new perspective for functional studies of gut bacteria in general, and may aid the development of new gene silencing-based pest management strategies.

Transcriptomic analysis of the fungus Graphilbum sp. in response to the pine wood nematode

Graphilbum species are important blue stain fungi associated with pine trees and are widely distributed throughout Asia, Australia, and North Africa. Pine wood nematode (PWN) primarily feed on ophiostomatoid fungi such as Graphilbum sp. in wood, the population of PWNs was increased, and incomplete organelle structures were observed in Graphilbum sp. hyphal cells following exposure to PWNs. In this study, we showed that Rho and Ras were involved in the MAPK pathway, SNARE binding and small GTPase-mediated signal transduction, and their expression was upregulated in the treatment group. However, the expression of the Rab7 involved in MAPK and small GTPase-mediated signal pathway was downregulated in the treatment group. Thus, further research is needed to study the MAPK pathway and related Ras and Rho genes in Graphilbum sp. associated with the PWN population. Overall, transcriptomic analysis clarified the basic mechanisms of mycelial growth in Graphilbum sp. fungus used as a food source by PWNs.

Metabarcoding of mycetangia from the Dendroctonus frontalis species complex (Curculionidae: Scolytinae) reveals diverse and functionally redundant fungal assemblages

Dendroctonus-bark beetles are associated with microbes that can detoxify terpenes, degrade complex molecules, supplement and recycle nutrients, fix nitrogen, produce semiochemicals, and regulate ecological interactions between microbes. Females of some Dendroctonus species harbor microbes in specialized organs called mycetangia; yet little is known about the microbial diversity contained in these structures. Here, we use metabarcoding to characterize mycetangial fungi from beetle species in the Dendroctonus frontalis complex, and analyze variation in biodiversity of microbial assemblages between beetle species. Overall fungal diversity was represented by 4 phyla, 13 classes, 25 orders, 39 families, and 48 genera, including 33 filamentous fungi, and 15 yeasts. The most abundant genera were Entomocorticium, Candida, Ophiostoma-Sporothrix, Ogataea, Nakazawaea, Yamadazyma, Ceratocystiopsis, Grosmannia-Leptographium, Absidia, and Cyberlindnera. Analysis of α-diversity indicated that fungal assemblages of D. vitei showed the highest richness and diversity, whereas those associated with D. brevicomis and D. barberi had the lowest richness and diversity, respectively. Analysis of β-diversity showed clear differentiation in the assemblages associated with D. adjunctus, D. barberi, and D. brevicomis, but not between closely related species, including D. frontalis and D. mesoamericanus and D. mexicanus and D. vitei. A core mycobiome was not statistically identified; however, the genus Ceratocystiopsis was shared among seven beetle species. Interpretation of a tanglegram suggests evolutionary congruence between fungal assemblages and species of the D. frontalis complex. The presence of different amplicon sequence variants (ASVs) of the same genus in assemblages from species of the D. frontalis complex outlines the complexity of molecular networks, with the most complex assemblages identified from D. vitei, D. mesoamericanus, D. adjunctus, and D. frontalis. Analysis of functional variation of fungal assemblages indicated multiple trophic groupings, symbiotroph/saprotroph guilds represented with the highest frequency (∼31% of identified genera). These findings improve our knowledge about the diversity of mycetangial communities in species of the D. frontalis complex and suggest that minimal apparently specific assemblages are maintained and regulated within mycetangia.

Phylogeny of Regulators of G-Protein Signaling Genes in Leptographium qinlingensis and Expression Levels of Three RGSs in Response to Different Terpenoids

Leptographium qinlingensis is a bark beetle-vectored pine pathogen in the Chinese white pine beetle (Dendroctonus armandi) epidemic in Northwest China. L. qinlingensis colonizes pines despite the trees’ massive oleoresin terpenoid defenses. Regulators of G-protein signaling (RGS) proteins modulate heterotrimeric G-protein signaling negatively and play multiple roles in the growth, asexual development, and pathogenicity of fungi. In this study, we have identified three L. qinlingensis RGS genes, and the phylogenetic analysis shows the highest homology with the regulators of G-protein signaling proteins sequence from Ophiostoma piceae and Grosmannia clavigera. The expression profiles of three RGSs in the mycelium of L. qinlingensis treated with six different terpenoids were detected, as well as their growth rates. Under six terpenoid treatments, the growth and reproduction in L. qinlingensis were significantly inhibited, and the growth inflection day was delayed from 8 days to 12–13 days. By analyzing the expression level of three RGS genes of L. qinlingensis with different treatments, results indicate that LqFlbA plays a crucial role in controlling fungal growth, and both LqRax1 and LqRgsA are involved in overcoming the host chemical resistances and successful colonization.

Pathophysiology and transcriptomic analysis of Picea koraiensis inoculated by bark beetle-vectored fungus Ophiostoma bicolor

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.