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

Chemistry, Biological Activities, and Pharmacological Properties of Gastrodin: Mechanism Insights

Gastrodin, a bioactive compound derived from the rhizome of the orchid Gastrodia elata, exhibits a diverse range of biological activities. With documented neuroprotective, anti-inflammatory, antioxidant, anti-apoptotic, and anti-tumor effects, gastrodin stands out as a multifaceted therapeutic agent. Notably, it has demonstrated efficacy in protecting against neuronal damage and enhancing cognitive function in animal models of Alzheimer's disease, Parkinson's disease, and cerebral ischemia. Additionally, gastrodin showcases immunomodulatory effects by mitigating inflammation and suppressing the expression of inflammatory cytokines. Its cytotoxic activity involves the inhibition of angiogenesis, suppression of tumor growth, and induction of apoptosis. This comprehensive review seeks to elucidate the myriad potential effects of Gastrodin, delving into the intricate molecular mechanisms underpinning its pharmacological properties. The findings underscore the therapeutic potential of gastrodin in addressing various conditions linked to neuroinflammation and cancer.

Perspective and challenges of mycorrhizal symbiosis in orchid medicinal plants

The family Orchidaceae is of the most diverse taxon in the plant kingdom, and most of its members are highly valuable herbal medicines. Orchids have a unique mycorrhizal symbiotic relationship with specific fungi for carbohydrate and nutrient supplies in their whole lifecycle. The large-scale cultivation of the medicinal plant Gastodia elata is a successful example of using mycorrhizal symbiotic technology. In this review, we adopted G. elata and Dendrobium officinale as examples to describe the characteristics of orchid mycorrhiza and mycorrhizal benefits for host plants' growth and health (e.g. biotic and abiotic stress and secondary metabolite accumulation). The challenges in applying mycorrhizal technology to the cultivation of orchid medicinal plants in the future were also discussed. This review aims to serve as a theoretical guide for the cultivation of mycorrhizal technology in medicinal orchid plants.

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.

Insights into the nutritional properties and molecular basis of biosynthesis of amino acids and vitamins of Gastrodia elata offered by metabolomic and transcriptomic analysis

Gastrodia elata Blume (GE), a traditional and precious Chinese medicinal material, has been approved as a functional food. However, understanding GE's nutritional properties and its molecular basis remains limited. Here, metabolomic and transcriptomic analyses were performed on young and mature tubers of G. elata.f.elata (GEEy and GEEm) and G. elata.f.glauca (GEGy and GEGm). A total of 345 metabolites were detected, including 76 different amino acids and their derivatives containing all human essential amino acids (e.g., l-(+)-lysine, l-leucine), 13 vitamins (e.g., nicotinamide, thiamine), and 34 alkaloids (e.g., spermine, choline). GEGm has higher amino acid accumulation than GEEy, GEEm and GEGy, and vitamin contents were also slightly different in all four samples. Implying that GE, especially GEGm, is a kind of excellent complementary food as amino acid nutrition provider. From assembled 21,513 transcripts (genes) based on the transcriptome, we identified many genes that encode enzymes (e.g., pfkA, bglX, tyrAa, lysA, his B, aroA), which are responsible for the biosynthesis of amino acids and enzymes (e.g., nadA, URH1, NAPRT1, punA, rsgA) that related to vitamins metabolism. A total of 16 pairs of the differentially expressed genes (DEG) and differentially accumulated metabolites (DAM) (e.g., gene-tia006709 coding GAPDH and l-(+)-arginine, and gene-tia010180 coding tyrA and l-(+)-arginine) and three DEG-DAM pairs (e.g., gene-tia015379 coding NadA and nicotinate d-ribonucleoside) show significant similar positive or negative correlation based on three, and two comparisons of GEEy vs. GEGy, GEGy vs. GEGm, GEEy vs. GEGy and GEEm vs. GEGm, which involved into amino acid biosynthesis, and nicotinate nicotinamide metabolism, respectively. These results prove that the enzyme coded by these DEG promotes (positive correlation) or inhibits (negative correlation) the biosynthesis of parallel DAM in GE. Overall, the data and corresponding analysis in this study provide new insights into the nutritional properties of GE and the related molecular basis.

One-pot biosynthesis of gastrodin using UDP-glycosyltransferase itUGT2 with an in situ UDP-glucose recycling system

Gastrodin, the major effective ingredient in Tianma (Gastrodia elata), is a p-hydroxybenzoic acid derivative with various activities. Gastrodin has been widely investigated for food and medical applications. The last biosynthetic step for gastrodin is UDP-glycosyltransferase (UGT)-mediated glycosylation with UDP-glucose (UDPG) as glycosyl donor. In this study, we performed a one-pot reaction both in vitro and in vivo to synthesize gastrodin from p-hydroxybenzyl alcohol (pHBA) by coupling UDP-glucosyltransferase from Indigofera tinctoria (itUGT2) to sucrose synthase from Glycine max (GmSuSy) for regeneration of UDPG. The in vitro results showed that itUGT2 transferred a glucosyl group to pHBA to generate gastrodin. After 37 UDPG regeneration cycles with 2.5% (molar ratio) UDP, the pHBA conversion reached 93% at 8 h. Furthermore, a recombinant strain with itUGT2 and GmSuSy genes was constructed. Through optimizing the incubation conditions, a 95% pHBA conversion rate (220 mg/L gastrodin titer) was achieved in vivo without addition of UDPG, which was 2.6-fold higher than that without GmSuSy. This in situ system for gastrodin biosynthesis provides a highly efficient strategy for both in vitro gastrodin synthesis and in vivo biosynthesis of gastrodin in E. coli with UDPG regeneration.

High-throughput sequencing-based analysis of the composition and diversity of endophytic bacteria community in tubers of Gastrodia elata f.glauca

Gastrodia elata f.glauca (G. elata) is a commonly used Chinese Medicinal Materials with great medicinal value. The medicinal plant and its endophytic bacteria are a symbiotic whole, and the endophytic bacteria are rich in species, and their metabolites are a treasure trove of natural compounds. However, there is a relative lack of analysis on the diversity, flora composition and network interactions of the endophytic bacteria of G. elata. In this study, high-throughput sequencing technology based on the Illumina Miseq platform was used to reveal the core microbiota by examining the diversity and community structures of tuber endophytic bacteria in G. elata grown under different regions and exploring the effect of region on its endophytic bacteria. Here, 1,265 endophytic ASVs were found to coexist with G. elata tuber in Guizhou and Hubei. At the phylum level, the dominant phyla were Proteobacteria, Actinobacteria and Acdobacteriota. At the family level, the dominant family were Comamonadaceae, Nocardicaece, Xanthobacteraceae, and Burkholderiaceae. At the genus level, Delftia and Rhodococcus were represented the core microbiota in G. elata tuber, which served as the dominant genera that coexisted in all samples tested. Moreover, we found that the beta diversity of endophytic bacteria in G. elata tuber was higher level in the Guizhou region than Hubei region. Overall, this study results to provide a reference for screening active strains and interaction between plants and endophytic bacteria.

Comparative transcriptome and metabolome analyses of four Panax species explore the dynamics of metabolite biosynthesis

Background

The genus Panax in the Araliaceae family has been used as traditional medicinal plants worldwide and is known to biosynthesize ginsenosides and phytosterols. However, genetic variation between Panax species has influenced their biosynthetic pathways is not fully understood.

Methods

Simultaneous analysis of transcriptomes and metabolomes obtained from adventitious roots of two tetraploid species (Panax ginseng and P. quinquefolius) and two diploid species (P. notoginseng and P. vietnamensis) revealed the diversity of their metabolites and related gene expression profiles.

Results

The transcriptome analysis showed that 2,3-OXIDOSQUALENE CYCLASEs (OSCs) involved in phytosterol biosynthesis are upregulated in the diploid species, while the expression of OSCs contributing to ginsenoside biosynthesis is higher in the tetraploid species. In agreement with these results, the contents of dammarenediol-type ginsenosides were higher in the tetraploid species relative to the diploid species.

Conclusion

These results suggest that a whole-genome duplication event has influenced the triterpene biosynthesis pathway in tetraploid Panax species during their evolution or ecological adaptation. This study provides a basis for further efforts to explore the genetic variation of the Panax genus.

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

BACKGROUND

Gastrodia elata Bl. f. glauca, a perennial herb of G.elata Bl. in Orchidaceae, is one of the most valuable traditional Chinese medicines. G. elata Bl. is a chlorophyll-free myco-heterotrophic plant, which must rely on the symbiotic growth of Armillaria, but not all Armillaria strains can play the symbiotic role. Additionally, Armillaria is easy to degenerate after multiple generations, and the compatibility between the strains from other areas and G. elata Bl. f. glauca in Changbai Mountain is unstable. Therefore, it is incredibly significant to isolate, identify and screen the symbiotic Armillaria suitable for the growth of G. elata Bl. f. glauca in Changbai Mountain, and to explore the mechanism by which Armillaria improves the production performance of G. elata Bl. f. glauca.

RESULTS

Firstly, G. elata Bl. f. glauca tubers, and rhizomorphs and fruiting bodies of Armillaria were used for the isolation and identification of Armillaria. Five Armillaria isolates were obtained in our laboratory and named: JMG, JMA, JMB, JMC and JMD. Secondly, Armillaria was selected based on the yield and the effective component content of G. elata Bl. f. glauca. It was concluded that the yield and quality of G. elata Bl. f. glauca co-planted with JMG is the highest. Finally, the mechanism of its high quality and yield was explored by investigating the effects of different Armillaria strains on the soil, its nutrition element contents and the soil microbial diversity around G. elata Bl. f. glauca in Changbai Mountain.

CONCLUSIONS

Compared with commercial strains, JMG significantly increased the content of Na, Al, Si, Mn, Fe, Zn, Rb and the absorption of C, Na, Mg, Ca, Cr, Cu, Zn and Rb in G. elata Bl. f. glauca; it improved the composition, diversity and metabolic functions of soil microbial communities around G. elata Bl. f. glauca at phylum, class and genus levels; it markedly increased the relative abundance of bacteria such as Chthoniobacter and Armillaria in the dominant populations, and enhanced such functions as Cell motility, amino acid metabolism and Lipid metabolism; it dramatically decreased the relative abundance of Bryobacter and other fungi in the dominant populations, and reduced such functions as microbial energy metabolism, translation and carbohydrate metabolism. This is the main reason why excellent Armillaria strains promote the high quality and yield of G. elata Bl. f. glauca in Changbai Mountain.

The gastrodin biosynthetic pathway in Pholidota chinensis Lindl. revealed by transcriptome and metabolome profiling

Pholidota chinensis Lindl. is an epiphytic or lithophytic perennial herb of Orchidaceae family used as a garden flower or medicinal plant to treat high blood pressure, dizziness and headache in traditional Chinese medicine. Gastrodin (GAS) is considered as a main bioactive ingredient of this herb but the biosynthetic pathway remains unclear in P. chinensis. To elucidate the GAS biosynthesis and identify the related genes in P. chinensis, a comprehensive analysis of transcriptome and metabolome of roots, rhizomes, pseudobulbs and leaves were performed by using PacBio SMART, Illumina Hiseq and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS/MS). A total of 1,156 metabolites were identified by UPLC-MS/MS, of which 345 differential metabolites were mainly enriched in phenylpropanoid/phenylalanine, flavone and flavonol biosynthesis. The pseudobulbs make up nearly half of the fresh weight of the whole plant, and the GAS content in the pseudobulbs was also the highest in four tissues. Up to 23,105 Unigenes were obtained and 22,029 transcripts were annotated in the transcriptome analysis. Compared to roots, 7,787, 8,376 and 9,146 differentially expressed genes (DEGs) were identified in rhizomes, pseudobulbs and leaves, respectively. And in total, 80 Unigenes encoding eight key enzymes for GAS biosynthesis, were identified. Particularly, glycosyltransferase, the key enzyme of the last step in the GAS biosynthetic pathway had 39 Unigenes candidates, of which, transcript28360/f2p0/1592, was putatively identified as the most likely candidate based on analysis of co-expression, phylogenetic analysis, and homologous searching. The metabolomics and transcriptomics of pseudobulbs versus roots showed that 8,376 DEGs and 345 DEMs had a substantial association based on the Pearson's correlation. This study notably enriched the metabolomic and transcriptomic data of P. chinensis, and it provides valuable information for GAS biosynthesis in the plant.

Chromosome-level genome assembly of the fully mycoheterotrophic orchid Gastrodia elata

Gastrodia elata, an obligate mycoheterotrophic orchid, requires complete carbon and mineral nutrient supplementation from mycorrhizal fungi during its entire life cycle. Although full mycoheterotrophy occurs most often in the Orchidaceae family, no chromosome-level reference genome from this group has been assembled to date. Here, we report a high-quality chromosome-level genome assembly of G. elata, using Illumina and PacBio sequencing methods with Hi-C technique. The assembled genome size was found to be 1045 Mb, with an N50 of 50.6 Mb and 488 scaffolds. A total of 935 complete (64.9%) matches to the 1440 embryophyte Benchmarking Universal Single-Copy Orthologs were identified in this genome assembly. Hi-C scaffolding of the assembled genome resulted in 18 pseudochromosomes, 1008 Mb in size and containing 96.5% of the scaffolds. A total of 18,844 protein-coding sequences (CDSs) were predicted in the G. elata genome, of which 15,619 CDSs (82.89%) were functionally annotated. In addition, 74.92% of the assembled genome was found to be composed of transposable elements. Phylogenetic analysis indicated a significant contraction of genes involved in various biosynthetic processes and cellular components and an expansion of genes for novel metabolic processes and mycorrhizal association. This result suggests an evolutionary adaptation of G. elata to a mycoheterotrophic lifestyle. In summary, the genomic resources generated in this study will provide a valuable reference genome for investigating the molecular mechanisms of G. elata biological functions. Furthermore, the complete G. elata genome will greatly improve our understanding of the genetics of Orchidaceae and its mycoheterotrophic evolution.

Armillaria mellea Symbiosis Drives Metabolomic and Transcriptomic Changes in Polyporus umbellatus Sclerotia

Sclerotia, the medicinal part of Polyporus umbellatus, play important roles in diuresis and renal protection, with steroids and polysaccharides as the main active ingredients. The sclerotia grow and develop only after symbiosis with Armillaria sp. In this study, a systematic metabolomics based on non-targeted UPLC-MS method was carried out between the infected part of the separated cavity wall of the sclerotia (QR) and the uninfected part (the control group, CK) to find and identify differential metabolites. The biosynthetic pathway of characteristic steroids in sclerotia of P. umbellatus was deduced and the content of ergosterol, polyporusterone A and B in the QR and CK groups were detected with the High Performance Liquid Chromatography (HPLC). Furthermore, the expression patterns of putative genes associated with steroid biosynthesis pathway were also performed with quantitative real-time PCR. The results showed that a total of 258 metabolites originated from fungi with the fragmentation score more than 45 and high resolution mass were identified, based on UPLC-MS metabolomic analysis, and there were 118 differentially expressed metabolites (DEMs) between both groups. The metabolic pathways indicated that steroids, fatty acid and carbohydrate were active and enriched during P. umbellatus sclerotia infected by A. mellea. The content of ergosterol, polyporusterone A and B in the QR group increased by 32.2, 75.0, and 20.0%, in comparison to that of the control group. The qRT-PCR analysis showed that series of enzymes including C-8 sterol isomerase (ERG2), sterol C-24 methyltransferase (ERG6) and sterol 22-desaturase (ERG5), which played important roles in the final steps of ergosterol biosynthesis, all presented up-regulated patterns in the QR group in P. umbellatus. The comprehensive metabolomic and transcriptomic information will contribute to further study concerning the mechanisms of P. umbellatus sclerotial formation infected by A. mellea in the future.

Comparative transcriptome analysis of tubers, stems, and flowers of Gastrodia elata Blume reveals potential genes involved in the biosynthesis of phenolics

Orchidaceae, well known for its fascinating flowers, is one of the largest and most diverse families of flowering plants. There are many kinds of plants in this family; these are distributed practically globally and have high ornamental and medicinal values. Gastrodia elata Blume, a traditional Chinese medicinal herb, is a rootless and leafless achlorophyllous orchid. Phenolic compounds are considered to be the major bioactive constituents in G. elata, with antioxidant, antiangiogenic, neuroprotective, antidepressant, anxiolytic, and sedative activities. In this study, we determined the contents of six main phenolic components in tubers, stems and flowers from G. elata. Meanwhile, the transcriptomes of the tuber, stem and flower tissues of G. elata were obtained using the BGISEQ-500 platform. A total of 58.29 Gb of data and 113,067 unigenes were obtained, of which 74,820 unigenes were functionally annotated against seven public databases. Differentially expressed genes between tuber, stem and flower tissues were identified. A total of 76 DEGs encoding eight key enzymes were identified as candidate genes involved in the biosynthesis of phenolics in G. elata. For further validation, the expression levels of unigenes were measured using quantitative real-time PCR. Our results greatly enrich the transcriptomic data of G. elata and provide valuable information for the identification of candidate genes involved in the biosynthesis of secondary metabolites.

Transcriptome Reveals Roles of Lignin-Modifying Enzymes and Abscisic Acid in the Symbiosis of Mycena and Gastrodia elata

Gastrodia elata is a well-known medicinal and heterotrophic orchid. Its germination, limited by the impermeability of seed coat lignin and inhibition by abscisic acid (ABA), is triggered by symbiosis with fungi such as Mycena spp. However, the molecular mechanisms of lignin degradation by Mycena and ABA biosynthesis and signaling in G. elata remain unclear. In order to gain insights into these two processes, this study analyzed the transcriptomes of these organisms during their dynamic symbiosis. Among the 25 lignin-modifying enzyme genes in Mycena, two ligninolytic class II peroxidases and two laccases were significantly upregulated, most likely enabling Mycena hyphae to break through the lignin seed coats of G. elata. Genes related to reduced virulence and loss of pathogenicity in Mycena accounted for more than half of annotated genes, presumably contributing to symbiosis. After coculture, upregulated genes outnumbered downregulated genes in G. elata seeds, suggesting slightly increased biological activity, while Mycena hyphae had fewer upregulated than downregulated genes, indicating decreased biological activity. ABA biosynthesis in G. elata was reduced by the downregulated expression of 9-cis-epoxycarotenoid dioxygenase (NCED-2), and ABA signaling was blocked by the downregulated expression of a receptor protein (PYL12-like). This is the first report to describe the role of NCED-2 and PYL12-like in breaking G. elata seed dormancy by reducing the synthesis and blocking the signaling of the germination inhibitor ABA. This study provides a theoretical basis for screening germination fungi to identify effective symbionts and for reducing ABA inhibition of G. elata seed germination.

Armillaria root rot fungi host single-stranded RNA viruses

Species of Armillaria are distributed globally and include some of the most important pathogens of forest and ornamental trees. Some of them form large long-living clones that are considered as one of the largest organisms on earth and are capable of long-range spore-mediated transfer as well as vegetative spread by drought-resistant hyphal cords called rhizomorphs. However, the virus community infecting these species has remained unknown. In this study we used dsRNA screening and high-throughput sequencing to search for possible virus infections in a collection of Armillaria isolates representing three different species: Armillaria mellea from South Africa, A. borealis from Finland and Russia (Siberia) and A. cepistipes from Finland. Our analysis revealed the presence of both negative-sense RNA viruses and positive-sense RNA viruses, while no dsRNA viruses were detected. The viruses included putative new members of virus families Mymonaviridae, Botourmiaviridae and Virgaviridae and members of a recently discovered virus group tentatively named "ambiviruses" with ambisense bicistronic genomic organization. We demonstrated that Armillaria isolates can be cured of viruses by thermal treatment, which enables the examination of virus effects on host growth and phenotype using isogenic virus-infected and virus-free strains.

GeSUT4 mediates sucrose import at the symbiotic interface for carbon allocation of heterotrophic Gastrodia elata (Orchidaceae)

Gastrodia elata, a fully mycoheterotrophic orchid without photosynthetic ability, only grows symbiotically with the fungus Armillaria. The mechanism of carbon distribution in this mycoheterotrophy is unknown. We detected high sucrose concentrations in all stages of Gastrodia tubers, suggesting sucrose may be the major sugar transported between fungus and orchid. Thick symplasm-isolated wall interfaces in colonized and adjacent large cells implied involvement of sucrose importers. Two sucrose transporter (SUT)-like genes, GeSUT4 and GeSUT3, were identified that were highly expressed in young Armillaria-colonized tubers. Yeast complementation and isotope tracer experiments confirmed that GeSUT4 functioned as a high-affinity sucrose-specific proton-dependent importer. Plasma-membrane/tonoplast localization of GeSUT4-GFP fusions and high RNA expression of GeSUT4 in symbiotic and large cells indicated that GeSUT4 likely functions in active sucrose transport for intercellular allocation and intracellular homeostasis. Transgenic Arabidopsis overexpressing GeSUT4 had larger leaves but were sensitive to excess sucrose and roots were colonized with fewer mutualistic Bacillus, supporting the role of GeSUT4 in regulating sugar allocation. This is not only the first documented carbon import system in a mycoheterotrophic interaction but also highlights the evolutionary importance of sucrose transporters for regulation of carbon flow in all types of plant-microbe interactions.

Metabolic engineering of Saccharomyces cerevisiae for high-level production of gastrodin from glucose

BACKGROUND

The natural phenolic glycoside gastrodin is the major bioactive ingredient in the well-known Chinese herb Tianma and is widely used as a neuroprotective medicine in the clinic. Microbial production from sustainable resources is a promising method to replace plant extraction and chemical synthesis which were currently used in industrial gastrodin production. Saccharomyces cerevisiae is considered as an attractive host to produce natural plant products used in the food and pharmaceutical fields. In this work, we intended to explore the potential of S. cerevisiae as the host for high-level production of gastrodin from glucose.

RESULTS

Here, we first identified the plant-derived glucosyltransferase AsUGT to convert 4-hydroxybenzyl alcohol to gastrodin with high catalytic efficiency in yeast. Then, we engineered de novo production of gastrodin by overexpressing codon-optimized AsUGT^syn, the carboxylic acid reductase gene CAR^syn from Nocardia species, the phosphopantetheinyl transferase gene PPTcg-1^syn from Corynebacterium glutamicum, the chorismate pyruvate-lyase gene UbiC^syn from Escherichia coli, and the mutant ARO4^K229L. Finally, we achieved an improved product titer by a chromosomal multiple-copy integration strategy and enhancement of metabolic flux toward the aglycon 4-hydroxybenzyl alcohol. The best optimized strain produced 2.1 g/L gastrodin in mineral medium with glucose as the sole carbon source by flask fermentation, which was 175 times higher than that of the original gastrodin-producing strain.

CONCLUSIONS

The de novo high-level production of gastrodin was first achieved. Instead of chemical synthesis or plants extraction, our work provides an alternative strategy for the industrial production of gastrodin by microbial fermentation from a sustainable resource.

GelFAP: Gene Functional Analysis Platform for Gastrodia elata

Gastrodia elata, also named Tianma, is a valuable traditional Chinese herbal medicine. It has numerous important pharmacological roles such as in sedation and lowering blood pressure and as anticonvulsant and anti-aging, and it also has effects on the immune and cardiovascular systems. The whole genome sequencing of G. elata has been completed in recent years, which provides a strong support for the construction of the G. elata gene functional analysis platform. Therefore, in our research, we collected and processed 39 transcriptome data of G. elata and constructed the G. elata gene co-expression networks, then we identified functional modules by the weighted correlation network analysis (WGCNA) package. Furthermore, gene families of G. elata were identified by tools including HMMER, iTAK, PfamScan, and InParanoid. Finally, we constructed a gene functional analysis platform for G. elata . In our platform, we introduced functional analysis tools such as BLAST, gene set enrichment analysis (GSEA), and cis-elements (motif) enrichment analysis tool. In addition, we analyzed the co-expression relationship of genes which might participate in the biosynthesis of gastrodin and predicted 19 mannose-binding lectin antifungal proteins of G. elata. We also introduced the usage of the G. elata gene function analysis platform (GelFAP) by analyzing CYP51G1 and GFAP4 genes. Our platform GelFAP may help researchers to explore the gene function of G. elata and make novel discoveries about key genes involved in the biological processes of gastrodin.

Draft genomic sequence of Armillaria gallica 012m: insights into its symbiotic relationship with Gastrodia elata

Armillaria species (Basidiomycota, Physalacriaceae) are well known as plant pathogens related to serious root rot disease on various trees in forests and plantations. Interestingly, some Armillaria species are essential symbionts of the rare Chinese medicinal herb Gastrodia elata, a rootless and leafless orchid used for over 2000 years. In this work, an 87.3-M draft genome of Armillaria gallica 012m strain, which was symbiotic with G. elata, was assembled. The genome includes approximately 23.6% repetitive sequences and encodes 26,261 predicted genes. In comparison with other four genomes of Armillaria, the following gene families related to pathogenicity/saprophytic phase, including cytochrome P450 monooxygenases, carbohydrate-active enzyme AA3, and hydrophobins, were significantly contracted in A. gallica 012m. These characteristics may be beneficial for G. elata to get less injuries. The genome-guided analysis of differential expression between rhizomorph (RH) and vegetative mycelium (VM) showed that a total of 2549 genes were differentially expressed, including 632 downregulated genes and 1917 upregulated genes. In the RH, most differentially expressed genes (DEGs) related to pathogenicity were significantly upregulated. To further elucidate gene function, Gene Ontology enrichment analysis showed that the upregulated DEGs significantly grouped into monooxygenase activity, hydrolase activity, glucosidase activity, extracellular region, fungal cell wall, response to xenobiotic stimulus, response to toxic substance, etc. These phenomena indicate that RH had better infection ability than VM. The infection ability of RH may be beneficial for G. elata to obtain nutrition, because the rhizomorph constantly infected the nutritional stems of G. elata and formed the hyphae that can be digested by G. elata. These results clarified the characteristics of A. gallica 012m and the reason why the strain 012m can establish a symbiotic relationship with G. elata in some extent from the perspective of genomics.

Volatile Organic Compounds from Orchids: From Synthesis and Function to Gene Regulation

Orchids are one of the most significant plants that have ecologically adapted to every habitat on earth. Orchids show a high level of variation in their floral morphologies, which makes them popular as ornamental plants in the global market. Floral scent and color are key traits for many floricultural crops. Volatile organic compounds (VOCs) play vital roles in pollinator attraction, defense, and interaction with the environment. Recent progress in omics technology has led to the isolation of genes encoding candidate enzymes responsible for the biosynthesis and regulatory circuits of plant VOCs. Uncovering the biosynthetic pathways and regulatory mechanisms underlying the production of floral scents is necessary not only for a better understanding of the function of relevant genes but also for the generation of new cultivars with desirable traits through molecular breeding approaches. However, little is known about the pathways responsible for floral scents in orchids because of their long life cycle as well as the complex and large genome; only partial terpenoid pathways have been reported in orchids. Here, we review the biosynthesis and regulation of floral volatile compounds in orchids. In particular, we focused on the genes responsible for volatile compounds in various tissues and developmental stages in Cymbidium orchids. We also described the emission of orchid floral volatiles and their function in pollination ecology. Taken together, this review will provide a broad scope for the study of orchid floral scents.

Transcriptome analysis of Clinopodium gracile (Benth.) Matsum and identification of genes related to Triterpenoid Saponin biosynthesis

BACKGROUND

Clinopodium gracile (Benth.) Matsum (C. gracile) is an annual herb with pharmacological properties effective in the treatment of various diseases, including hepatic carcinoma. Triterpenoid saponins are crucial bioactive compounds in C. gracile. However, the molecular understanding of the triterpenoid saponin biosynthesis pathway remains unclear.

RESULTS

In this study, we performed RNA sequencing (RNA-Seq) analysis of the flowers, leaves, roots, and stems of C. gracile plants using the BGISEQ-500 platform. The assembly of transcripts from all four types of tissues generated 128,856 unigenes, of which 99,020 were mapped to several public databases for functional annotation. Differentially expressed genes (DEGs) were identified via the comparison of gene expression levels between leaves and other tissues (flowers, roots, and stems). Multiple genes encoding pivotal enzymes, such as squalene synthase (SS), or transcription factors (TFs) related to triterpenoid saponin biosynthesis were identified and further analyzed. The expression levels of unigenes encoding important enzymes were verified by quantitative real-time PCR (qRT-PCR). Different chemical constituents of triterpenoid saponins were identified by Ultra-Performance Liquid Chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC/Q-TOF-MS).

CONCLUSIONS

Our results greatly extend the public transcriptome dataset of C. gracile and provide valuable information for the identification of candidate genes involved in the biosynthesis of triterpenoid saponins and other important secondary metabolites.

Orchids and their mycorrhizal fungi: an insufficiently explored relationship

Orchids are associated with diverse fungal taxa, including nonmycorrhizal endophytic fungi as well as mycorrhizal fungi. The orchid mycorrhizal (OM) symbiosis is an excellent model for investigating the biological interactions between plants and fungi due to their high dependency on these symbionts for growth and survival. To capture the complexity of OM interactions, significant genomic, numerous transcriptomic, and proteomic studies have been performed, unraveling partly the role of each partner. On the other hand, several papers studied the bioactive metabolites from each partner but rarely interpreted their significance in this symbiotic relationship. In this review, we focus from a biochemical viewpoint on the OM dynamics and its molecular interactions. The ecological functions of OM in plant development and stress resistance are described first, summarizing recent literature. Secondly, because only few studies have specifically looked on OM molecular interactions, the signaling pathways and compounds allowing the establishment/maintenance of mycorrhizal association involved in arbuscular mycorrhiza (AM) are discussed in parallel with OM. Based on mechanistic similarities between OM and AM, and recent findings on orchids' endophytes, a putative model representing the different molecular strategies that OM fungi might employ to establish this association is proposed. It is hypothesized here that (i) orchids would excrete plant molecule signals such as strigolactones and flavonoids but also other secondary metabolites; (ii) in response, OM fungi would secrete mycorrhizal factors (Myc factors) or similar compounds to activate the common symbiosis genes (CSGs); (iii) overcome the defense mechanism by evasion of the pathogen-associated molecular patterns (PAMPs)-triggered immunity and by secretion of effectors such as small inhibitor proteins; and (iv) finally, secrete phytohormones to help the colonization or disrupt the crosstalk of plant defense phytohormones. To challenge this putative model, targeted and untargeted metabolomics studies with special attention to each partner's contribution are finally encouraged and some technical approaches are proposed.

Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome

Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.

Identification and Analysis of Potential Genes Regulated by an Alphasatellite (TYLCCNA) that Contribute to Host Resistance against Tomato Yellow Leaf Curl China Virus and Its Betasatellite (TYLCCNV/TYLCCNB) Infection in Nicotiana benthamiana

Recently, begomovirus/betasatellite disease complexes were found to be associated with alphasatellites, and their presence modulated disease symptoms and/or viral DNA accumulation in infected plants. However, the biological functions of alphasatellites during begomovirus/betasatellite infections remain unclear. Tomato yellow leaf curl China virus (TYLCCNV) associated with a betasatellite (TYLCCNB) is a widespread monopartite begomovirus in China. In the Yunnan province of China, the TYLCCNV/TYLCCNB disease complex is found in association with an alphasatellite (TYLCCNA). In this study, in order to explain the mechanisms underlying TYLCCNV/TYLCCNB infection and reductions in viral DNA accumulation caused by TYLCCNA, we analyzed the transcriptome profiles of Nicotiana benthamiana seedlings challenged by TYLCCNV/TYLCCNB or TYLCCNV/TYLCCNB/TYLCCNA using RNA sequencing. In total, 2272 and 1207 differentially expressed genes (DEGs) were identified to respond to TYLCCNV/TYLCCNB and TYLCCNV/TYLCCNB/TYLCCNA infections, respectively. Compared with the DEGs in the TYLCCNV/TYLCCNB-infected N. benthamiana seedlings, the number of DEGs in plants co-infected with TYLCCNV/TYLCCNB + TYLCCNA was significantly reduced. Additionally, 36 DEGs were identified to be regulated by TYLCCNA, six of which were further analyzed using the virus-induced gene silencing (VIGS) approach. Silencing of these six TYLCCNA responsive DEGs caused more severe disease symptoms and higher viral DNA accumulation levels, suggesting that TYLCCNA responsive DEGs may attenuate TYLCCNV/TYLCCNB infection.

Inducing Secondary Metabolite Production by Co-culture of the Endophytic Fungus Phoma sp. and the Symbiotic Fungus Armillaria sp

Co-culturing the endophytic fungus Phoma sp. YUD17001 from Gastrodia elata with Armillaria sp. in liquid nutrient medium resulted in the production of five new secondary metabolites, including two phenolic compounds, phexandiols A and B (1 and 2), three aliphatic ester derivatives, phomesters A-C (3-5), and a known fatty acid (6). The structures and absolute configurations of these compounds were elucidated by the interpretation of data from detailed spectroscopic analysis, Mosher's method, and electronic circular dichroism spectra, together with consideration of the biogenetic origins. None of the five new compounds were detected in single-strain cultures under identical fermentation conditions. The results of this work indicated that the production of 1-5 involved a complicated interaction process. None of the new compounds possessed significant cytotoxicity or antimicrobial activities.

New Insights into the Symbiotic Relationship between Orchids and Fungi

Mycorrhizas play an important role in plant growth and development. In mycorrhizal symbioses, fungi supply soil mineral nutrients, such as nitrogen and phosphorus, to their host plants in exchange for carbon resources. Plants gain as much as 80% of mineral nutrient requirements from mycorrhizal fungi, which form associations with the roots of over 90% of all plant species. Orchid seeds lack endosperms and contain very limited storage reserves. Therefore, the symbiosis with mycorrhizal fungi that form endomycorrhizas is essential for orchid seed germination and protocorm development under natural conditions. The rapid advancement of next-generation sequencing contributes to identifying the orchid and fungal genes involved in the orchid mycorrhizal symbiosis and unraveling the molecular mechanisms regulating the symbiosis. We aim to update and summarize the current understanding of the mechanisms on orchid-fungus symbiosis, and the main focus will be on the nutrient exchange between orchids and their fungal partners.

The Gastrodia elata genome provides insights into plant adaptation to heterotrophy

We present the 1.06 Gb sequenced genome of Gastrodia elata, an obligate mycoheterotrophic plant, which contains 18,969 protein-coding genes. Many genes conserved in other plant species have been deleted from the G. elata genome, including most of those for photosynthesis. Additional evidence of the influence of genome plasticity in the adaptation of this mycoheterotrophic lifestyle is evident in the large number of gene families that are expanded in G. elata, including glycoside hydrolases and urease that likely facilitate the digestion of hyphae are expanded, as are genes associated with strigolactone signaling, and ATPases that may contribute to the atypical energy metabolism. We also find that the plastid genome of G. elata is markedly smaller than that of green plant species while its mitochondrial genome is one of the largest observed to date. Our report establishes a foundation for studying adaptation to a mycoheterotrophic lifestyle.

The Gastrodia elata genome provides insights into plant adaptation to heterotrophy

We present the 1.06 Gb sequenced genome of Gastrodia elata, an obligate mycoheterotrophic plant, which contains 18,969 protein-coding genes. Many genes conserved in other plant species have been deleted from the G. elata genome, including most of those for photosynthesis. Additional evidence of the influence of genome plasticity in the adaptation of this mycoheterotrophic lifestyle is evident in the large number of gene families that are expanded in G. elata, including glycoside hydrolases and urease that likely facilitate the digestion of hyphae are expanded, as are genes associated with strigolactone signaling, and ATPases that may contribute to the atypical energy metabolism. We also find that the plastid genome of G. elata is markedly smaller than that of green plant species while its mitochondrial genome is one of the largest observed to date. Our report establishes a foundation for studying adaptation to a mycoheterotrophic lifestyle.

Revealing proteins associated with symbiotic germination of Gastrodia elata by proteomic analysis

BACKGROUND

Gastrodia elata, a mycoheterotrophic orchid, is a well-known medicinal herb. In nature, the seed germination of G. elata requires proper fungal association, because of the absence of endosperm. To germinate successfully, G. elata obtains nutrition from mycorrhizal fungi such as Mycena. However, Mycena is not able to supply nutrition for the further development and enlargement of protocorms into tubers, flowering and fruit setting of G. elata. To date, current genomic studies on this topic are limited. Here we used the proteomic approach to explore changes in G. elata at different stages of symbiotic germination.

RESULTS

Using mass spectrometry, 3787 unique proteins were identified, of which 599 were classified as differentially accumulated proteins. Most of these differentially accumulated proteins were putatively involved in energy metabolism, plant defense, molecular signaling, and secondary metabolism. Among them, the defense genes (e.g., pathogenesis-/wound-related proteins, peroxidases, and serine/threonine-protein kinase) were highly expressed in late-stage protocorms, suggesting that fungal colonization triggered the significant defense responses of G. elata.

CONCLUSIONS

The present study indicated the metabolic change and defensive reaction could disrupt the balance between Mycena and G. elata during mycorrhizal symbiotic germination.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Transcriptomic analyses reveal clathrin-mediated endocytosis involved in symbiotic seed germination of Gastrodia elata

BACKGROUND

Gastrodia elata is a well-known medicinal orchid. In nature, the germination rate of G. elata is extremely poor, because there is no endosperm within the mature seed. It is crucial for G. elata to obtain nutrition from mycorrhizal fungi (Mycena) at the early-stage of germination. After germination, the seed gives rise to a protocorm. However, there are no "omic" studies on understanding the interaction between Gastrodia and Mycena. Here, we used transcriptomic approaches to explore changes in seed germination of G. elata.

RESULTS

Based on RNA-Seq, a total of ~221 million clean reads were assembled denovo into 139,756 unigenes, including 42,140 unigenes that were annotated in public databases. Meanwhile, 1750 unigenes were identified as differentially expressed genes. Most of these differentially expressed genes were putatively involved in energy metabolism, plant defense, molecular signaling, and secondary metabolism. Additionally, numerous genes involved in clathrin-mediated endocytosis were identified from our data. Most of these genes (e.g., clathrin, adaptor protein, dynamin, HSC70) were basally expressed in seeds and highly expressed in protocorms.

CONCLUSIONS

Our data suggested that clathrin-mediated endocytosis could play important roles in symbiotic seed germination of G. elata with Mycena infections.

Transcriptome Analysis of Salicylic Acid Treatment in Rehmannia glutinosa Hairy Roots Using RNA-seq Technique for Identification of Genes Involved in Acteoside Biosynthesis

Rehmannia glutinosa is a common bulk medicinal material that has been widely used in China due to its active ingredients. Acteoside, one of the ingredients, has antioxidant, antinephritic, anti-inflammatory, hepatoprotective, immunomodulatory, and neuroprotective effects, is usually selected as a quality-control component for R. glutinosa herb in the Chinese Pharmacopeia. The acteoside biosynthesis pathway in R. glutinosa has not yet been clearly established. Herein, we describe the establishment of a genetic transformation system for R. glutinosa mediated by Agrobacterium rhizogenes. We screened the optimal elicitors that markedly increased acteoside accumulation in R. glutinosa hairy roots. We found that acteoside accumulation dramatically increased with the addition of salicylic acid (SA); the optimal SA dose was 25 μmol/L for hairy roots. RNA-seq was applied to analyze the transcriptomic changes in hairy roots treated with SA for 24 h in comparison with an untreated control. A total of 3,716, 4,018, and 2,715 differentially expressed transcripts (DETs) were identified in 0 h-vs.-12 h, 0 h-vs.-24 h, and 12 h-vs.-24 h libraries, respectively. KEGG pathway-based analysis revealed that 127 DETs were enriched in "phenylpropanoid biosynthesis." Of 219 putative unigenes involved in acteoside biosynthesis, 54 were found to be up-regulated at at least one of the time points after SA treatment. Selected candidate genes were analyzed by quantitative real-time PCR (qRT-PCR) in hairy roots with SA, methyl jasmonate (MeJA), AgNO₃ (Ag⁺), and putrescine (Put) treatment. All genes investigated were up-regulated by SA treatment, and most candidate genes were weakly increased by MeJA to some degree. Furthermore, transcription abundance of eight candidate genes in tuberous roots of the high-acteoside-content (HA) cultivar QH were higher than those of the low-acteoside-content (LA) cultivar Wen 85-5. These results will pave the way for understanding the molecular basis of acteoside biosynthesis in R. glutinosa, and can serve as a basis for future validation studies.