Highly diversified fungi are associated with the achlorophyllous orchid Gastrodia flavilabella

Metagenomics and metatranscriptomics as potential driving forces for the exploration of diversity and functions of micro-eukaryotes in soil

Micro-eukaryotes are ubiquitous and play vital roles in diverse ecological systems, yet their diversity and functions are scarcely known. This may be due to the limitations of formerly used conventional culture-based methods. Metagenomics and metatranscriptomics are enabling to unravel the genomic, metabolic, and phylogenetic diversity of micro-eukaryotes inhabiting in different ecosystems in a more comprehensive manner. The in-depth study of structural and functional characteristics of micro-eukaryote community residing in soil is crucial for the complete understanding of this major ecosystem. This review provides a deep insight into the methodologies employed under these approaches to study soil micro-eukaryotic organisms. Furthermore, the review describes available computational tools, pipelines, and database sources and their manipulation for the analysis of sequence data of micro-eukaryotic origin. The challenges and limitations of these approaches are also discussed in detail. In addition, this review summarizes the key findings of metagenomic and metatranscriptomic studies on soil micro-eukaryotes. It also highlights the exploitation of these methods to study the structural as well as functional profiles of soil micro-eukaryotic community and to screen functional eukaryotic protein coding genes for biotechnological applications along with the future perspectives in the field.

An exploration of mechanism of high quality and yield of Gastrodia elata Bl. f. glauca by the isolation, identification, and evaluation of Mycena

Gastrodia elata Bl. f. glauca is an important traditional Chinese medicinal plant. The yield and quality of Gastrodia elata Bl. have significantly decreased due to multigenerational asexual reproduction. Therefore, it is necessary to have sexual reproduction of Gastrodia elata Bl. to supplement the market supply. Seeds of G. elata Bl. have no endosperm, and their sexual reproduction depends on the nutrients provided by the embryo cells infected by Mycena fungi to complete seed germination. However, Mycena fungi are small and have many species, and not all Mycena fungi can promote the germination of G. elata Bl. seeds. Therefore, it is of great significance to isolate and identify suitable germination fungi and explore the mechanism for improving the production performance and yield, and quality of G. elata Bl. Six closely related Mycena isolates, JFGL-01, JFGL-02, JFGL-03, JFGL-04, JFGL-05, and JFGL-06, were isolated from the leaves and protocorms of G. elata Bl. f. glauca and were identified as Mycena purpureofusca. The mycelial state and number of germinating protocorms were used as indicators to preferentially select Mycena fungi, and it was concluded that JFGL-06 had the best mycelial state and ability to germinate G. elata Bl. seeds. Finally, a mechanism to increase the yield of G. elata Bl. was explored by comparing the changes in nutrient elements and microbial diversity in the soil around G. elata Bl. with different strains. JFGL-06 proved to be an excellent Mycena fungal strain suitable for G. elata Bl. f. glauca. Compared with the commercial strain, JFGL-06 significantly increased the C, N, Na, Mg, S, Cl, K, Ca, and Fe contents of the soil surrounding the protocorms of G. elata Bl. f. glauca. JFGL-06 improved the composition, diversity, and metabolic function of the surrounding soil microbial community of G. elata Bl. f. glauca protocorms at the phylum, class, and genus levels, significantly increased the relative abundance of bacteria such as Acidobacteria and fungi such as Trichoderma among the dominant groups, and increased the abundance of functional genes in metabolic pathways such as nucleotide metabolism and energy metabolism. There was a significant reduction in the relative abundance of bacteria, such as Actinomycetes, and fungi, such as Fusarium, in the dominant flora, and a reduced abundance of functional genes, such as amino acid metabolism and xenobiotic biodegradation and metabolism. This is the main reason why the JFGL-06 strain promoted high-quality and high-yield G. elata Bl. f. glauca in Changbai Mountain.

Characterization of the transcriptional responses of Armillaria gallica 012m to GA3

Gastrodia elata needs to establish a symbiotic relationship with Armillaria strains to obtain nutrients and energy. However, the signaling cross talk between G. elata and Armillaria strains is still unclear. During our experiment, we found that the vegetative mycelium of Armillaria gallica 012m grew significantly better in the media containing gibberellic acid (GA3) than the blank control group (BK). To explore the response mechanism, we performed an RNA-sequencing experiment to profile the transcriptome changes of A. gallica 012m cultured in the medium with exogenous GA3. The transcriptome-guided differential expression genes (DEGs) analysis of GA3 and BK showed that a total of 1309 genes were differentially expressed, including 361 upregulated genes and 948 downregulated genes. Some of those DEGs correlated with the biological process, including positive regulation of chromosome segregation, mitotic metaphase/anaphase transition, attachment of mitotic spindle microtubules to kinetochore, mitotic cytokinesis, and nuclear division. These analyses explained that GA3 actively promoted the growth of A. gallica to some extent. Further analysis of protein domain features showed that the deduced polypeptide contained 41 candidate genes of GA receptor, and 27 of them were expressed in our samples. We speculate that GA receptors exist in A. gallica 012m. Comparative studies of proteins showed that the postulated GA receptor domains of A. gallica 012m have a higher homologous correlation with fungi than others based on cluster analysis.

Structure and diversity of mycorrhizal fungi communities of different part of Bulbophyllum tianguii in three terrestrial environments

Mycorrhizal fungi plays important roles in the seed germination and subsequent growth of orchids. The research of fungi in orchid roots, especially dominant mycorrhizal fungi is critical for orchids protection. In this study, the fungal community and composition of mycorrhizal fungi in roots, rhizomes and rhizosphere soil of Bulbophyllum tianguii grown in three terrestrial environments were analyzed by the second generation sequencing technology. The results of OTU clustering and α and β diversity analysis showed that there were significant differences in fungal communities in roots, rhizomes and rhizosphere soil of B. tianguii. The total number of OTUs in rhizomes was much less than that in roots and rhizosphere soil. The number of OTUs in rhizosphere soil and the diversity of mycorrhizal fungi were the highest. Meanwhile, the species and abundance of mycorrhizal fungi in roots and rhizomes of B. tianguii were different from those in rhizosphere soil. For different elevations, compared with B. tianguii that grow in middle of Tiankeng and top of Tiankeng, the OTUs number of B. tianguii in orchid garden is richest, and the diversity of mycorrhizal fungi in orchid garden was significantly higher than other locations. Among the three different habitats of B. tianguii, the number of OTUs in humus soil and stone habitats was notably higher than tree habitats, and the diversity of mycorrhizal fungi in humus soil was the highest. The analysis of mycorrhizal fungi in different habitats and altitudes of B. tianguii showed that Sebacina and Exophiala were the dominant mycorrhizal fungi in B. tianguii. The results of species annotation, phylogenetic tree and co-occurrence network analysis showed the dominant mycorrhizal fungi of B. tianguii mainly included Sebacina, Cladosporium, Exophiala, Fusarium. This study reveals the symbiotic relationship between Sebacina, Exophiala, Cladosporium and the B. Tianguii. It will provide a theoretical basis for the protection and biological function study of B. Tianguii.

The selective anti-fungal metabolites from Irpex lacteus and applications in the chemical interaction of Gastrodia elata, Armillaria sp., and endophytes

The investigation of the metabolites from endophyte Irpex lacteus cultured in host "tian ma" (Gastrodia elata) revealed five new tremulane sesquiterpenes (1-5), and a new tetrahydrofuran derivative (6). Compound 1 was the first tremulane glucoside, and 6 possessed a rare tetrahydropyran-tetrahydrofuran scaffold. Main metabolite (2,3-dihydroxydodacane-4,7-dione, 14) from I. lacteus showed significant selectivity for antifungal activity against phytopathogen and endophytes associated with G. elata rather than against Armillaria sp. providing nutrition for the host G. elata. 14 accounted for 27.4% of isolated compounds from G. elata medium, and 69.3% by co-culturing with Armillaria sp. So the I. lacteus tended to promote the growth of Armillaria sp. in co-culture by producing 2,3-dihydroxydodacane-4,7-dione (14) to selective inhibit the phytopathogen and endophyte existed in host G. elata for the benefit of G. elata-Armillaria symbiosis. And the results were in accord with the real environment of G. elata depending on the nutrition of Armillaria. Some metabolites had anti-fungal activities against phytopathogens of G. elata with MICs ≤8 μg/mL.

Progress and Prospects of Mycorrhizal Fungal Diversity in Orchids

Orchids form mycorrhizal symbioses with fungi in natural habitats that affect their seed germination, protocorm growth, and adult nutrition. An increasing number of studies indicates how orchids gain mineral nutrients and sometime even organic compounds from interactions with orchid mycorrhizal fungi (OMF). Thus, OMF exhibit a high diversity and play a key role in the life cycle of orchids. In recent years, the high-throughput molecular identification of fungi has broadly extended our understanding of OMF diversity, revealing it to be a dynamic outcome co-regulated by environmental filtering, dispersal restrictions, spatiotemporal scales, biogeographic history, as well as the distribution, selection, and phylogenetic spectrum width of host orchids. Most of the results show congruent emerging patterns. Although it is still difficult to extend them to all orchid species or geographical areas, to a certain extent they follow the "everything is everywhere, but the environment selects" rule. This review provides an extensive understanding of the diversity and ecological dynamics of orchid-fungal association. Moreover, it promotes the conservation of resources and the regeneration of rare or endangered orchids. We provide a comprehensive overview, systematically describing six fields of research on orchid-fungal diversity: the research methods of orchid-fungal interactions, the primer selection in high-throughput sequencing, the fungal diversity and specificity in orchids, the difference and adaptability of OMF in different habitats, the comparison of OMF in orchid roots and soil, and the spatiotemporal variation patterns of OMF. Further, we highlight certain shortcomings of current research methodologies and propose perspectives for future studies. This review emphasizes the need for more information on the four main ecological processes: dispersal, selection, ecological drift, and diversification, as well as their interactions, in the study of orchid-fungal interactions and OMF community structure.

Evolutionary histories and mycorrhizal associations of mycoheterotrophic plants dependent on saprotrophic fungi

Mycoheterotrophic plants (MHPs) are leafless, achlorophyllous, and completely dependent on mycorrhizal fungi for their carbon supply. Mycorrhizal symbiosis is a mutualistic association with fungi that is undertaken by the majority of land plants, but mycoheterotrophy represents a breakdown of this mutualism in that plants parasitize fungi. Most MHPs are associated with fungi that are mycorrhizal with autotrophic plants, such as arbuscular mycorrhizal (AM) or ectomycorrhizal (ECM) fungi. Although these MHPs gain carbon via the common mycorrhizal network that links the surrounding autotrophic plants, some mycoheterotrophic lineages are associated with saprotrophic (SAP) fungi, which are free-living and decompose leaf litter and wood materials. Such MHPs are dependent on the forest carbon cycle, which involves the decomposition of wood debris and leaf litter, and have a unique biology and evolutionary history. MHPs associated with SAP fungi (SAP-MHPs) have to date been found only in the Orchidaceae and likely evolved independently at least nine times within that family. Phylogenetically divergent SAP Basidiomycota, mostly Agaricales but also Hymenochaetales, Polyporales, and others, are involved in mycoheterotrophy. The fungal specificity of SAP-MHPs varies from a highly specific association with a single fungal species to a broad range of interactions with multiple fungal orders. Establishment of symbiotic culture systems is indispensable for understanding the mechanisms underlying plant-fungus interactions and the conservation of MHPs. Symbiotic culture systems have been established for many SAP-MHP species as a pure culture of free-living SAP fungi is easier than that of biotrophic AM or ECM fungi. Culturable SAP-MHPs are useful research materials and will contribute to the advancement of plant science.

Root-Associated Fungal Communities in Two Populations of the Fully Mycoheterotrophic Plant Arachnitis uniflora Phil. (Corsiaceae) in Southern Chile

The microbiological interactions of the roots of non-photosynthetic plants in South America have been scarcely explored. This study analyzes culturable fungal diversity associated with the mycoheterotrophic plant Arachnitis uniflora Phil. (Corsiaceae) in southern Chile, growing in two different understoreys of native (Nothofagus-dominated) and mixed forest (native, Cupressus sempervirens, and Pinus radiata). Rhizospheric and endophytic fungi were isolated, cultured, and purified to identify microorganisms associated with A. uniflora roots. We showed the different fungi associated with the plant, and that these distributions are influenced by the sampling site. We isolated 410 fungal strains (144 endophytic and 266 from the rhizosphere). We identified 13 operative taxonomical units from plants sampled in the mixed forest, while 15 were from the native forest. Rhizospheric microorganisms were mainly related to Penicillium spp., whereas some pathogenic and saprophytic strains were more frequent inside the roots. Our results have also shown that the fungal strains are weak for phosphate solubilization, but other pathways such as organic acid exudation and indole acetic acid production can be considered as major mechanisms to stimulate plant growth. Our results point to new fungal associates of A. uniflora plants reported in Andean ecosystems, identifying new beneficial endophytic fungi associated with roots of this fully mycoheterotrophic plant.

Dynamics of fungal communities during Gastrodia elata growth

BACKGROUND

Gastrodia elata is a widely distributed achlorophyllous orchid and is highly valued as both medicine and food. Gastrodia elata produces dust-like seeds and relies on mycorrhizal fungi for its germination and growth. In its life cycle, G. elata is considered to switch from a specific single-fungus relationship (Mycena) to another single-fungus relationship (Armillaria). However, no studies have investigated the changes in the plant-fungus relationship during the growth of G. elata in the wild. In this study, high-throughput sequencing was used to characterize the fungal community of tubers in different growth phases as well as the soils surrounding G. elata.

RESULTS

The predominant fungi were Basidiomycota (60.44%) and Ascomycota (26.40%), which exhibited changes in abundance and diversity with the growth phases of G. elata. Diverse basidiomycetes in protocorms (phase P) were Hyphodontia, Sistotrema, Tricholoma, Mingxiaea, Russula, and Mycena, but the community changed from a large proportion of Resinicium bicolor (40%) in rice-like tubers (phase M) to an unidentified Agaricales operational taxonomic unit 1(OTU1,98.45%) in propagation vegetation tubers (phase B). The soil fungi primarily included Simocybe, Psathyrella, Conocybe, and Subulicystidium. Three Mycena OTUs obtained in this study were differentially distributed among the growth phases of G. elata, accounting for less than 1.0% of the total reads, and were phylogenetically close to Mycena epipterygia and M. alexandri.

CONCLUSIONS

Our data indicated that G. elata interacts with a broad range of fungi beyond the Mycena genus. These fungi changed with the growth phases of G. elata. In addition, these data suggested that the development of the fungal community during the growth of G. elata was more complex than previously assumed and that at least two different fungi could be involved in development before the arrival of Armillaria.

Highly diverse endophytes in roots of Cycas bifida (Cycadaceae), an ancient but endangered gymnosperm

As an ancient seed plant, cycads are one of the few gymnosperms that develop a root symbiosis with cyanobacteria, which has allowed cycads to cope with harsh geologic and climatic conditions during the evolutionary process. However, the endophytic microbes in cycad roots remain poorly identified. In this study, using next-generation sequencing techniques, we investigated the microbial diversity and composition of both the coralloid and regular roots of Cycas bifida (Dyer) K.D. Hill. Highly diverse endophytic communities were observed in both the coralloid and regular roots. Of the associated bacteria, the top five families were the Nostocaceae, Sinobacteraceae, Bradyrhizobiaceae, Bacillaceae, and Hyphomicrobiaceae. The Nectriaceae, Trichocomaceae, and Incertae sedis were the predominant fungal families in all root samples. A significant difference in the endophytic bacterial community was detected between coralloid roots and regular roots, but no difference was observed between the fungal communities in the two root types. Cyanobacteria were more dominant in coralloid roots than in regular roots. The divergence of cycad root structures and the modified physiological processes may have contributed to the abundance of cyanobionts in coralloid roots. Consequently, the colonization of cyanobacteria inhibits the assemblage of other endophytes. Our results contribute to an understanding of the species diversity and composition of the cycad-endophyte microbiome and provide an abbreviated list of potential ecological roles of the core microbes present.

Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes

BACKGROUND

Gastrodia elata Blume (Orchidaceae) is an important Chinese medicine with several functional components. In the life cycle of G. elata, the orchid develops a symbiotic relationship with two compatible mycorrhizal fungi Mycena spp. and Armillaria mellea during seed germination to form vegetative propagation corm and vegetative growth to develop tubers, respectively. Gastrodin (p-hydroxymethylphenol-beta-D-glucoside) is the most important functional component in G. elata, and gastrodin significantly increases from vegetative propagation corms to tubers. To address the gene regulation mechanism in gastrodin biosynthesis in G. elata, a comparative analysis of de novo transcriptome sequencing among the vegetative propagation corms and tubers of G. elata and A. mellea was conducted using deep sequencing.

RESULTS

Transcriptome comparison between the vegetative propagation corms and juvenile tubers of G. elata revealed 703 differentially expressed unigenes, of which 298 and 405 unigenes were, respectively up-regulated (fold-change ≥ 2, q-value < 0.05, the trimmed mean of M-values (TMM)-normalized fragments per kilobase of transcript per Million mapped reads (FPKM) > 10) and down-regulated (fold-change ≤ 0.5, q-value <0.05, TMM-normalized FPKM > 10) in juvenile tubers. After Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 112 up-regulated unigenes with KEGG Ortholog identifiers (KOids) or enzyme commission (EC) numbers were assigned to 159 isogroups involved in seventy-eight different pathways, and 132 down-regulated unigenes with KOids or EC numbers were assigned to 168 isogroups, involved in eighty different pathways. The analysis of the isogroup genes from all pathways revealed that the two unigenes TRINITY_DN54282_c0_g1 (putative monooxygenases) and TRINITY_DN50323_c0_g1 (putative glycosyltransferases) might participate in hydroxylation and glucosylation in the gastrodin biosynthetic pathway.

CONCLUSIONS

The gene expression of the two unique unigenes encoding monooxygenase and glycosyltransferase significantly increases from vegetative propagation corms to tubers, and the molecular basis of gastrodin biosynthesis in the tubers of G. elata is proposed.

Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.