The Gastrodia elata genome provides insights into plant adaptation to heterotrophy

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera

PREMISE

Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants.

METHODS

Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities.

RESULTS

The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae.

CONCLUSIONS

Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Orchids acquire fungal carbon for seed germination: pathways and players

To germinate in nature, orchid seeds strictly rely on seed germination-promoting orchid mycorrhizal fungi (sgOMFs) for provision of carbon nutrients. The underlying delivery pathway, however, remains elusive. We develop here a plausible model for sugar transport from sgOMFs to orchid embryonic cells to fuel germination. Orchids exploit sgOMFs to induce the formation of pelotons, elaborate intracellular hyphal coils in orchid embryos. The colonized orchid cells then obtain carbon nutrients by uptake from living hyphae and peloton lysis, primarily as glucose derived from fungal trehalose hydrolyzed by orchid-specific trehalases. The uptake of massive fungally derived glucose is likely to be mediated by two classes of membrane proteins, namely, sugars will eventually be exported transporters (SWEETs) and H+-hexose symporters. The proposed model serves as a launch pad for further research to better understand and improve orchid seed germination and conservation.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Composition and diversity of soil microbial communities change by introducing Phallus impudicus into a Gastrodia elata Bl.-based soil

BACKGROUND

The Gastrodia elata Bl. is an orchid, and its growth demands the presence of Armillaria species. The strong competitiveness of Armillaria species has always been a concern of major threat to other soil organisms, thus disrupting the equilibrium of soil biodiversity. Introducing other species to where G. elata was cultivated, could possibly alleviate the problems associated with the disequilibrium of soil microenvironment; however, their impacts on the soil microbial communities and the underlying mechanisms remain unclear. To reveal the changes of microbial groups associated with soil chemical properties responding to different cultivation species, the chemical property measurements coupled with the next-generation pyrosequencing analyses were applied with soil samples collected from fallow land, cultivation of G. elata and Phallus impudicus, respectively.

RESULTS

The cultivation of G. elata induced significant increases (p < 0.05) in soil pH and NO3-N content compared with fallow land, whereas subsequent cultivation of P. impudicus reversed these G. elata-induced increases and was also found to significantly increase (p < 0.05) the content of soil NH4+-N and AP. The alpha diversities of soil microbial communities were significantly increased (p < 0.01) by cultivation of G. elata and P. impudicus as indicated with Chao1 estimator and Shannon index. The structure and composition of soil microbial communities differed responding to different cultivation species. In particular, the relative abundances of Bacillus, norank_o_Gaiellales, Mortierella and unclassified_k_Fungi were significantly increased (p < 0.05), while the abundances of potentially beneficial genera such as Acidibacter, Acidothermus, Cryptococcus, and Penicillium etc., were significantly decreased (p < 0.05) by cultivation of G. elata. It's interesting to find that cultivation of P. impudicus increased the abundances of these genera that G. elata decreased before, which contributed to the difference of composition and structure. The results of CCA and heatmap indicated that the changes of soil microbial communities had strong correlations with soil nutrients. Specifically, among 28 genera presented, 50% and 42.9% demonstrated significant correlations with soil pH and NO3-N in response to cultivation of G. elata and P. impudicus.

CONCLUSIONS

Our findings suggested that the cultivation of P. impudicus might have potential benefits as result of affecting soil microorganisms coupled with changes in soil nutrient profile.

Discovery of seed germinating fungi (Mycetinis scorodonius) from Gastrodia elata Bl. f. glauca S. chow in Changbai Mountain and examination of their germination ability

Multi-generational asexual reproduction of Gastrodia elata Bl. will cause seedling species degeneration. Sexual reproduction of Gastrodia elata Bl. seed is an effective method to solve the problem of degeneration. The development of Gastrodia elata Bl. seeds cannot be separated from the germination fungus. However, there are few strains of germination fungus in production, and there is also the problem of species degradation in application for many years. It is very important for the sexual reproduction of Gastrodia elata Bl. to isolate more new strains of excellent germination fungus from the origin. This study used the Gastrodia elata Bl. f. glauca S. chow seeds germination vegetative propagation corms capture method to isolate its symbiotic germination fungus, and comprehensively identified the species of germination fungus by colony morphology, ITS, sporocarps regeneration and germination function, and compared the growth characteristics and germination ability with other germination fungus (Mycena purpureofusca, Mycena dendrobii and Mycena osmundicola). The germination fungus was isolated from the vegetative propagation corms of Gastrodia elata Bl. f. glauca S. chow seeds and named GYGL-1. After comprehensive identification, GYGL-1 was Mycetinis scorodonius. Compared with other germination fungus, GYGL-1 has fast germination speed, vigorous growth, and high germination ability for Gastrodia elata Bl. f. glauca S. chow seeds. Innovated the isolation method of Gastrodia elata Bl. seeds germination fungus, obtained the regenerated sporocarps of the germination fungus, and discovered that Mycetinis scorodonius has a new function of germinating Gastrodia elata Bl. f. glauca S. chow seeds, enriching the resource library of Gastrodia elata Bl. germination fungus.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BACKGROUND

Orchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.

RESULTS

Here, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.

CONCLUSION

In this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

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.

Identification and Characterization of the AREB/ABF Gene Family in Three Orchid Species and Functional Analysis of DcaABI5 in Arabidopsis

AREB/ABF (ABA response element binding) proteins in plants are essential for stress responses, while our understanding of AREB/ABFs from orchid species, important traditional medicinal and ornamental plants, is limited. Here, twelve AREB/ABF genes were identified within three orchids' complete genomes and classified into three groups through phylogenetic analysis, which was further supported with a combined analysis of their conserved motifs and gene structures. The cis-element analysis revealed that hormone response elements as well as light and stress response elements were widely rich in the AREB/ABFs. A prediction analysis of the orchid ABRE/ABF-mediated regulatory network was further constructed through cis-regulatory element (CRE) analysis of their promoter regions. And it revealed that several dominant transcriptional factor (TF) gene families were abundant as potential regulators of these orchid AREB/ABFs. Expression profile analysis using public transcriptomic data suggested that most AREB/ABF genes have distinct tissue-specific expression patterns in orchid plants. Additionally, DcaABI5 as a homolog of ABA INSENSITIVE 5 (ABI5) from Arabidopsis was selected for further analysis. The results showed that transgenic Arabidopsis overexpressing DcaABI5 could rescue the ABA-insensitive phenotype in the mutant abi5. Collectively, these findings will provide valuable information on AREB/ABF genes in orchids.

Identification and Analysis of PEPC Gene Family Reveals Functional Diversification in Orchidaceae and the Regulation of Bacterial-Type PEPC

Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating the PEPC gene family in CAM plants. In this study, a total of 33 PEPC genes were identified across 15 orchids. Specifically, one PEPC gene was found in Cymbidium goeringii and Platanthera guangdongensis; two in Apostasia shenzhenica, Dendrobium chrysotoxum, D. huoshanense, Gastrodia elata, G. menghaiensis, Phalaenopsis aphrodite, Ph. equestris, and Pl. zijinensis; three in C. ensifolium, C. sinense, D. catenatum, D. nobile, and Vanilla planifolia. These PEPC genes were categorized into four subgroups, namely PEPC-i, PEPC-ii, and PEPC-iii (PTPC), and PEPC-iv (BTPC), supported by the comprehensive analyses of their physicochemical properties, motif, and gene structures. Remarkably, PEPC-iv contained a heretofore unreported orchid PEPC gene, identified as VpPEPC4. Differences in the number of PEPC homolog genes among these species were attributed to segmental duplication, whole-genome duplication (WGD), or gene loss events. Cis-elements identified in promoter regions were predominantly associated with light responsiveness, and circadian-related elements were observed in each PEPC-i and PEPC-ii gene. The expression levels of recruited BTPC, VpPEPC4, exhibited a lower expression level than other VpPEPCs in the tested tissues. The expression analyses and RT-qPCR results revealed diverse expression patterns in orchid PEPC genes. Duplicated genes exhibited distinct expression patterns, suggesting functional divergence. This study offered a comprehensive analysis to unveil the evolution and function of PEPC genes in Orchidaceae.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Organelle Genomes of Epipogium roseum Provide Insight into the Evolution of Mycoheterotrophic Orchids

Epipogium roseum, commonly known as one of the ghost orchids due to its rarity and almost transparent color, is a non-photosynthetic and fully mycoheterotrophic plant. Given its special nutritional strategies and evolutionary significance, the mitogenome was first characterized, and three plastomes sampled from Asia were assembled. The plastomes were found to be the smallest among Orchidaceae, with lengths ranging from 18,339 to 19,047 bp, and exhibited high sequence variety. For the mitogenome, a total of 414,552 bp in length, comprising 26 circular chromosomes, were identified. A total of 54 genes, including 38 protein-coding genes, 13 tRNA genes, and 3 rRNA genes, were annotated. Multiple repeat sequences spanning a length of 203,423 bp (45.47%) were discovered. Intriguingly, six plastid regions via intracellular gene transfer and four plastid regions via horizontal gene transfer to the mitogenome were observed. The phylogenomics, incorporating 90 plastomes and 56 mitogenomes, consistently revealed the sister relationship of Epipogium and Gastrodia, with a bootstrap percentage of 100%. These findings shed light on the organelle evolution of Orchidaceae and non-photosynthetic plants.

The role of symbiotic fungi in the life cycle of Gastrodia elata Blume (Orchidaceae): a comprehensive review

Gastrodia elata Blume, a fully mycoheterotrophic perennial plant of the family Orchidaceae, is a traditional Chinese herb with medicinal and edible value. Interestingly, G. elata requires symbiotic relationships with Mycena and Armillaria strains for seed germination and plant growth, respectively. However, there is no comprehensive summary of the symbiotic mechanism between fungi and G. elata. Here, the colonization and digestion of hyphae, the bidirectional exchange of nutrients, the adaptation of fungi and G. elata to symbiosis, and the role of microorganisms and secondary metabolites in the symbiotic relationship between fungi and G. elata are summarized. We comprehensively and deeply analyzed the mechanism of symbiosis between G. elata and fungi from three perspectives: morphology, nutrition, and molecules. The aim of this review was to enrich the understanding of the mutualistic symbiosis mechanisms between plants and fungi and lay a theoretical foundation for the ecological cultivation of G. elata.

Effect of light conditions on trophic level and gene expression of partially mycoheterotrophic orchid, Cymbidium goeringii

Partial mycoheterotrophic i.e., mixotrophic, plants are the species which partially depend on mycorrhizal fungi for its nutrients. Although some of these plants are known to show plasticity in the degree of fungal dependence induced by the changes in light condition, the genetic background of this plasticity is largely unsolved. Here, we investigated the relationships between environmental conditions and nutrient sources based on ¹³C and ¹⁵N enrichment in mixotrophic orchid Cymbidium goeringii. We also shaded them for 2 months and evaluated the effect of light condition on the nutrient sources based on the abundance of ¹³C and ¹⁵N and the gene expressions by RNA-seq based de novo assembly. The shading had no effect on isotope enrichment, possibly because of the translocation of carbon and nitrogen from the storage organs. Gene expression analysis showed the upregulation of genes involved in jasmonic acid response in leaves of the shaded plants, which suggests that the jasmonic acid played an important role in regulation of degree of dependence against the mycorrhizal fungi. Our results suggest that mixotrophic plants might be controlling their dependency against the mycorrhizal fungi by a common mechanism with the autotrophic plants.

Mitochondrial genome comparison and phylogenetic analysis of Dendrobium (Orchidaceae) based on whole mitogenomes

BACKGROUND

Mitochondrial genomes are essential for deciphering the unique evolutionary history of seed plants. However, the rules of their extreme variation in genomic size, multi-chromosomal structure, and foreign sequences remain unresolved in most plant lineages, which further hindered the application of mitogenomes in phylogenetic analyses.

RESULTS

Here, we took Dendrobium (Orchidaceae) which shows the great divergence of morphology and difficulty in species taxonomy as the study focus. We first de novo assembled two complete mitogenomes of Dendrobium wilsonii and Dendrobium henanense that were 763,005 bp and 807,551 bp long with multichromosomal structures. To understand the evolution of Dendrobium mitogenomes, we compared them with those of four other orchid species. The results showed great variations of repetitive and chloroplast-derived sequences in Dendrobium mitogenomes. Moreover, the intergenic content of Dendrobium mitogenomes has undergone expansion during evolution. We also newly sequenced mitogenomes of 26 Dendrobium species and reconstructed phylogenetic relationships of Dendrobium based on genomic mitochondrial and plastid data. The results indicated that the existence of chloroplast-derived sequences made the mitochondrial phylogeny display partial characteristics of the plastid phylogeny. Additionally, the mitochondrial phylogeny provided new insights into the phylogenetic relationships of Dendrobium species.

CONCLUSIONS

Our study revealed the evolution of Dendrobium mitogenomes and the potential of mitogenomes in deciphering phylogenetic relationships at low taxonomic levels.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Genome-based identification of the CYP75 gene family in Orchidaceae and its expression patterns in Cymbidium goeringii

With a great diversity of species, Orchidaceae stands out as an essential component of plant biodiversity, making it a primary resource for studying angiosperms evolution and genomics. This study focuses on 13 published orchid genomes to identify and analyze the CYP75 gene family belonging to the cytochrome P450 superfamily, which is closely related to flavonoid biosynthetic enzymes and pigment regulation. We found 72 CYP75s in the 13 orchid genomes and further classified them into two classes: CYP75A and CYP75B subfamily, the former synthesizes blue anthocyanins, while the latter is involved in the production of red anthocyanins. Furthermore, the amount of CYP75Bs (53/72) greatly exceeds the amount of CYP75As (19/72) in orchids. Our findings suggest that CYP75B genes have a more important evolutionary role, as red plants are more common in nature than blue plants. We also discovered unique conserved motifs in each subfamily that serve as specific recognition features (motif 19 belong to CYP75A; motif 17 belong to CYP75B). Two diverse-colored varieties of C. goeringii were selected for qRT-PCR experiments. The expression of CgCYP75B1 was significantly higher in the purple-red variant compared to the yellow-green variant, while CgCYP75A1 showed no significant difference. Based on transcriptomic expression analysis, CYP75Bs are more highly expressed than CYP75As in floral organs, especially in colorful petals and lips. These results provide valuable information for future studies on CYP75s in orchids and other angiosperms.

Effect of symbiotic fungi-Armillaria gallica on the yield of Gastrodia elata Bl. and insight into the response of soil microbial community

Armillaria members play important roles in the nutrient supply and growth modulation of Gastrodia elata Bl., and they will undergo severe competition with native soil organisms before colonization and become symbiotic with G. elata. Unraveling the response of soil microbial organisms to symbiotic fungi will open up new avenues to illustrate the biological mechanisms driving G. elata's benefit from Armillaria. For this purpose, Armillaria strains from four main G. elata production areas in China were collected, identified, and co-planted with G. elata in Guizhou Province. The result of the phylogenetic tree indicated that the four Armillaria strains shared the shortest clade with Armillaria gallica. The yields of G. elata were compared to uncover the potential role of these A. gallica strains. Soil microbial DNA was extracted and sequenced using Illumina sequencing of 16S and ITS rRNA gene amplicons to decipher the changes of soil bacterial and fungal communities arising from A. gallica strains. The yield of G. elata symbiosis with the YN strain (A. gallica collected from Yunnan) was four times higher than that of the GZ strain (A. gallica collected from Guizhou) and nearly two times higher than that of the AH and SX strains (A. gallica collected from Shanxi and Anhui). We found that the GZ strain induced changes in the bacterial community, while the YN strain mainly caused changes in the fungal community. Similar patterns were identified in non-metric multidimensional scaling analysis, in which the GZ strain greatly separated from others in bacterial structure, while the YN strain caused significant separation from other strains in fungal structure. This current study revealed the assembly and response of the soil microbial community to A. gallica strains and suggested that exotic strains of A. gallica might be helpful in improving the yield of G. elata by inducing changes in the soil fungal community.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green 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.

Metabolomics analysis of mycelial exudates provides insights into fungal antagonists of Armillaria

The genus Armillaria has high edible and medical values, with zones of antagonism often occurring when different species are paired in culture on agar media, while the antagonism-induced metabolic alteration remains unclear. Here, the metabolome of mycelial exudates of two Chinese Armillaria biological species, C and G, co-cultured or cultured separately was analysed to discover the candidate biomarkers and the key metabolic pathways involved in Armillaria antagonists. A total of 2,377 metabolites were identified, mainly organic acids and derivatives, lipids and lipid-like molecules, and organoheterocyclic compounds. There were 248 and 142 differentially expressed metabolites between group C-G and C, C-G, and G, respectively, and fourteen common differentially expressed metabolites including malate, uracil, Leu-Gln-Arg, etc. Metabolic pathways like TCA cycle and pyrimidine metabolism were significantly affected by C-G co-culture. Additionally, 156 new metabolites (largely organic acids and derivatives) including 32 potential antifungal compounds, primarily enriched into biosynthesis of secondary metabolites pathways were identified in C-G co-culture mode. We concluded that malate and uracil could be used as the candidate biomarkers, and TCA cycle and pyrimidine metabolism were the key metabolic pathways involved in Armillaria antagonists. The metabolic changes revealed in this study provide insights into the mechanisms underlying fungal antagonists.

Hallmarks of Comparative Transcriptome between Rhizomorphs and Hyphae of Armillaria sp. 541 Participating in Fungal Symbiosis with Emphasis on LysM Domains

Armillaria sp. 541, a genus of root-infecting fungi, forms a symbiosis with traditional Chinese medicine Gastrodia elata (Orchid) and Polyporus umbellatus via extensive networks of durable rhizomorphs. It is not clear the hallmarks of comparative transcriptome between the rhizomorphs and hyphae of Armillaria sp. 541. In the present study, transcriptomic analysis of Armillaria sp. 541 identified 475 differentially expressed genes (DEGs) between Armillaria rhizomorphs (AR) and hyphae (AH). Of them, 285 genes were upregulated and 190 were downregulated. Bioinformatics analyses and tests demonstrated DEGs involved in oxidoreductase activity and peptidoglycan binding were significantly enriched in this process when rhizomorph formed from hyphae. We accordingly obtained 14 gene-encoding proteins containing the LysM domain, and further consensus pattern and phylogenetic analysis indicated that their amino acid sequences were conserved and their biological functions may be peptidoglycan binding for recognition between the fungus and host. Among these genes, one, named Armillaria LysM domain recognition gene (aLDRG), was expressed significantly when rhizomorphs were differentiated from hyphae. It was located in the cortical cells of the rhizomorph by in situ hybridization. Furthermore, biolayer interferometry (BLI) assay demonstrated that aLDRG can bind specifically to chitin oligosaccharide of the fungal cell wall, including N,N',N″-Triacetylchitotriose (CO3) and N,N',N″,N'″,N″″-Pentaacetylchitopentaose (CO5). Therefore, we deduced that Armillaria sp. 541 expressed higher levels of LysM protein aLDRG for better binding of oligosaccharide after rhizomorphs were generated. This study provides functional genes for further studies on the interaction between Armillaria sp. 541 and its host.

Ancient Horizontal Gene Transfers from Plastome to Mitogenome of a Nonphotosynthetic Orchid, Gastrodia pubilabiata (Epidendroideae, Orchidaceae)

Gastrodia pubilabiata is a nonphotosynthetic and mycoheterotrophic orchid belonging to subfamily Epidendroideae. Compared to other typical angiosperm species, the plastome of G. pubilabiata is dramatically reduced in size to only 30,698 base pairs (bp). This reduction has led to the loss of most photosynthesis-related genes and some housekeeping genes in the plastome, which now only contains 19 protein coding genes, three tRNAs, and three rRNAs. In contrast, the typical orchid species contains 79 protein coding genes, 30 tRNAs, and four rRNAs. This study decoded the entire mitogenome of G. pubilabiata, which consisted of 44 contigs with a total length of 867,349 bp. Its mitogenome contained 38 protein coding genes, nine tRNAs, and three rRNAs. The gene content of G. pubilabiata mitogenome is similar to the typical plant mitogenomes even though the mitogenome size is twice as large as the typical ones. To determine possible gene transfer events between the plastome and the mitogenome individual BLASTN searches were conducted, using all available orchid plastome sequences and flowering plant mitogenome sequences. Plastid rRNA fragments were found at a high frequency in the mitogenome. Seven plastid protein coding gene fragments (ndhC, ndhJ, ndhK, psaA, psbF, rpoB, and rps4) were also identified in the mitogenome of G. pubilabiata. Phylogenetic trees using these seven plastid protein coding gene fragments suggested that horizontal gene transfer (HGT) from plastome to mitogenome occurred before losses of photosynthesis related genes, leading to the lineage of G. pubilabiata. Compared to species phylogeny of the lineage of orchid, it was estimated that HGT might have occurred approximately 30 million years ago.

Genome-Wide Identification and Expression Analysis of the SPL Gene Family in Three Orchids

SPL transcription factors regulate important processes such as plant growth and development, metabolic regulation, and abiotic stress. They play crucial roles in the development of flower organs. However, little is known about the characteristics and functions of the SPLs in the Orchidaceae. In this study, Cymbidium goeringii Rchb. f., Dendrobium chrysotoxum Lindl., and Gastrodia elata BI. were used as research objects. The SPL gene family of these orchids was analyzed on a genome-wide scale, and their physicochemical properties, phylogenetic relationships, gene structures, and expression patterns were studied. Transcriptome and qRT-PCR methods were combined to investigate the regulatory effect of SPLs on the development of flower organs during the flowering process (bud, initial bloom, and full bloom). This study identifies a total of 43 SPLs from C. goeringii (16), D. chrysotoxum (17), and G. elata (10) and divides them into eight subfamilies according to the phylogenetic tree. Most SPL proteins contained conserved SBP domains and complex gene structures; half of the genes had introns longer than 10 kb. The largest number and variety of cis-acting elements associated with light reactions were enriched, accounting for about 45% of the total (444/985); 13/43 SPLs contain response elements of miRNA156. GO enrichment analysis showed that the functions of most SPLs were mainly enriched in the development of plant flower organs and stems. In addition, expression patterns and qRT-PCR analysis suggested the involvement of SPL genes in the regulation of flower organ development in orchids. There was little change in the expression of the CgoSPL in C. goeringii, but DchSPL9 and GelSPL2 showed significant expression during the flowering process of D. chrysotoxum and G. elata, respectively. In summary, this paper provides a reference for exploring the regulation of the SPL gene family in orchids.

Whole genome sequencing and analysis of Armillaria gallica Jzi34 symbiotic with Gastrodia elata

BACKGROUND

Armillaria species are plant pathogens, but a few Armillaria species can establish a symbiotic relationship with Gastrodia elata, a rootless and leafless orchid, that is used as a Chinese herbal medicine. Armillaria is a nutrient source for the growth of G. elata. However, there are few reports on the molecular mechanism of symbiosis between Armillaria species and G. elata. The genome sequencing and analysis of Armillaria symbiotic with G. elata would provide genomic information for further studying the molecular mechanism of symbiosis.

RESULTS

The de novo genome assembly was performed with the PacBio Sequel platform and Illumina NovaSeq PE150 for the A. gallica Jzi34 strain, which was symbiotic with G. elata. Its genome assembly contained ~ 79.9 Mbp and consisted of 60 contigs with an N50 of 2,535,910 bp. There were only 4.1% repetitive sequences in the genome assembly. Functional annotation analysis revealed a total of 16,280 protein coding genes. Compared with the other five genomes of Armillaria, the carbohydrate enzyme gene family of the genome was significantly contracted, while it had the largest set of glycosyl transferase (GT) genes. It also had an expansion of auxiliary activity enzymes AA3-2 gene subfamily and cytochrome P450 genes. The synteny analysis result of P450 genes reveals that the evolutionary relationship of P450 proteins between A. gallica Jzi34 and other four Armillaria was complex.

CONCLUSIONS

These characteristics may be beneficial for establishing a symbiotic relationship with G. elata. These results explore the characteristics of A. gallica Jzi34 from a genomic perspective and provide an important genomic resource for further detailed study of Armillaria. This will help to further study the symbiotic mechanism between A. gallica and G. elata.

Rethinking convergence in plant parasitism through the lens of molecular and population genetic processes

The autotrophic lifestyle of photosynthetic plants has profoundly shaped their body plan, physiology, and gene repertoire. Shifts to parasitism and heterotrophy have evolved at least 12 times in more than 4000 species, and this transition has consequently left major evolutionary footprints among these parasitic lineages. Features that are otherwise rare at the molecular level and beyond have evolved repetitively, including reduced vegetative bodies, carrion-mimicking during reproduction, and the incorporation of alien genetic material. Here, I propose an integrated conceptual model, referred to as the funnel model, to define the general evolutionary trajectory of parasitic plants and provide a mechanistic explanation for their convergent evolution. This model connects our empirical understanding of gene regulatory networks in flowering plants with classical theories of molecular and population genetics. It emphasizes that the cascading effects brought about by the loss of photosynthesis may be a major force constraining the physiological capacity of parasitic plants and shaping their genomic landscapes. Here I review recent studies on the anatomy, physiology, and genetics of parasitic plants that lend support to this photosynthesis-centered funnel model. Focusing on nonphotosynthetic holoparasites, I elucidate how they may inevitably reach an evolutionary terminal status (i.e., extinction) and highlight the utility of a general, explicitly described and falsifiable model for future studies of parasitic plants.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

Apostasia Mitochondrial Genome Analysis and Monocot Mitochondria Phylogenomics

Apostasia shenzhenica belongs to the subfamily Apostasioideae and is a primitive group located at the base of the Orchidaceae phylogenetic tree. However, the A. shenzhenica mitochondrial genome (mitogenome) is still unexplored, and the phylogenetic relationships between monocots mitogenomes remain unexplored. In this study, we discussed the genetic diversity of A. shenzhenica and the phylogenetic relationships within its monocotyledon mitogenome. We sequenced and assembled the complete mitogenome of A. shenzhenica, resulting in a circular mitochondrial draft of 672,872 bp, with an average read coverage of 122× and a GC content of 44.4%. A. shenzhenica mitogenome contained 36 protein-coding genes, 16 tRNAs, two rRNAs, and two copies of nad4L. Repeat sequence analysis revealed a large number of medium and small repeats, accounting for 1.28% of the mitogenome sequence. Selection pressure analysis indicated high mitogenome conservation in related species. RNA editing identified 416 sites in the protein-coding region. Furthermore, we found 44 chloroplast genomic DNA fragments that were transferred from the chloroplast to the mitogenome of A. shenzhenica, with five plastid-derived genes remaining intact in the mitogenome. Finally, the phylogenetic analysis of the mitogenomes from A. shenzhenica and 28 other monocots showed that the evolution and classification of most monocots were well determined. These findings enrich the genetic resources of orchids and provide valuable information on the taxonomic classification and molecular evolution of monocots.

Fungal gut microbiota dysbiosis in systemic lupus erythematosus

Introduction

Despite recent developments in our comprehension of how the gut microbiota and systemic lupus erythematosus (SLE) are related. The mycobiome: which is a small but crucial part of the gut microbiota and is involved in hosts’ homeostasis and physiological processes, remained unexplored in SLE.

Methods

We profiled the gut fungal mycobiota based on internal transcribed spacer region 1 (ITS1) sequencing for the gut microbial DNA from the SLE individuals with lupus nephritis (LN) (n = 23), SLE without LN (n = 26) and healthy controls (n = 14) enrolled in Nanjing Drum Tower Hospital, The Affiliated Hospital of Nanjing University Medical School.

Results

The ITS sequencing generated a total of 4.63 million valid tags which were stratified into 4,488 operational taxonomic units (OTUs) and identified about 13 phyla and 262 genera. Patients with SLE were characterized with unique fungal flora feature. The fungal microbiomes of the three groups displayed distinct beta diversity from each other. Compared with HC group, the abundance of fungal dysbiosis was reflected in a higher ratio of opportunistic fungi in SLE or LN group, as well as the loss of Rhizopus and Malassezia. The main principal components of the flora between the SLE and LN group were generally consistent. The relative abundance of Vanrija in the fecal fungal community was higher in LN group, while the relative abundance of Fusarium was higher in SLE group. Moreover, our data revealed superior diagnostic accuracy for SLE with the fungal species (e.g. Candida, Meyerozyma). Correlations between gut fungi and clinical parameters were identified by Spearman’s correlation analysis. Interestingly, Aspergillus in SLE patients was positively correlated with ACR, 24 h proteinuria, proteinuria, anti-dsDNA, ANA, and SLEDAI, while Rhizopus was negatively correlated with lymphocytes and Hb. Finally, we successfully cultured the fungi and identified it as Candida glabrata by microscopic observation and mass spectrometry.

Discussion

We first explored the highly significant gut fungal dysbiosis and ecology in patients with SLE, and demonstrated the applicability of fungal species as SLE diagnostic tools, signifying that the gut fungal mycobiome-host interplay can potentially contribute in disease pathogenesis.

GelFAP v2.0: an improved platform for Gene functional analysis in Gastrodia elata

BACKGROUND

Gastrodia elata (tianma), a well-known medicinal orchid, is widely used to treat various kinds of diseases with its dried tuber. In recent years, new chromosome-level genomes of G.elata have been released in succession, which offer an enormous resource pool for understanding gene function. Previously we have constructed GelFAP for gene functional analysis of G.elata. As genomes are updated and transcriptome data is accumulated, collection data in GelFAP cannot meet the need of researchers.

RESULTS

Based on new chromosome-level genome and transcriptome data, we constructed co-expression network of G. elata, and then we annotated genes by aligning with sequences from NR, TAIR, Uniprot and Swissprot database. GO (Gene Ontology) and KEGG (Kyoto Encylopaedia of Genes and Genomes) annotations were predicted by InterProScan and GhostKOALA software. Gene families were further predicted by iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid. Finally, we developed an improved platform for gene functional analysis in G. elata (GelFAP v2.0) by integrating new genome, transcriptome data and processed functional annotation. Several tools were also introduced to platform including BLAST (Basic Local Alignment Search Tool), GSEA (Gene Set Enrichment Analysis), Heatmap, JBrowse, Motif analysis and Sequence extraction. Based on this platform, we found that the flavonoid biosynthesis might be regulated by transcription factors (TFs) such as MYB, HB and NAC. We also took C4H and GAFP4 as examples to show the usage of our platform.

CONCLUSION

An improved platform for gene functional analysis in G. elata (GelFAP v2.0, www.gzybioinformatics.cn/Gelv2 ) was constructed, which provides better genome data, more transcriptome resources and more analysis tools. The updated platform might be preferably benefit researchers to carry out gene functional research for their project.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

Evolutionary differences in gene loss and pseudogenization among mycoheterotrophic orchids in the tribe Vanilleae (subfamily Vanilloideae)

INTRODUCTION

Galeola lindleyana is a mycoheterotrophic orchid belonging to the tribe Vanilleae within the subfamily Vanilloideae.

METHODS

In this study, the G. lindleyana plastome was assembled and annotated, and compared with other Vanilleae orchids, revealing the evolutionary variations between the photoautotrophic and mycoheterotrophic plastomes.

RESULTS

The G. lindleyana plastome was found to include 32 protein-coding genes, 16 tRNA genes and four ribosomal RNA genes, including 11 pseudogenes. Almost all of the genes encoding photosynthesis have been lost physically or functionally, with the exception of six genes encoding ATP synthase and psaJ in photosystem I. The length of the G. lindleyana plastome has decreased to 100,749 bp, while still retaining its typical quadripartite structure. Compared with the photoautotrophic Vanilloideae plastomes, the inverted repeat (IR) regions and the large single copy (LSC) region of the mycoheterotrophic orchid's plastome have contracted, while the small single copy (SSC) region has expanded significantly. Moreover, the difference in length between the two ndhB genes was found to be 682 bp, with one of them spanning the IRb/SSC boundary. The Vanilloideae plastomes were varied in their structural organization, gene arrangement, and gene content. Even the Cyrtosia septentrionalis plastome which was found to be closest in length to the G. lindleyana plastome, differed in terms of its gene arrangement and gene content. In the LSC region, the psbA, psbK, atpA and psaB retained in the G. lindleyana plastome were missing in the C. septentrionalis plastome, while, the matK, rps16, and atpF were incomplete in the C. septentrionalis plastome, yet still complete in that of the G. lindleyana. Lastly, compared with the G. lindleyana plastome, a 15 kb region located in the SSC area between ndhB-rrn16S was found to be inverted in the C. septentrionalis plastome. These changes in gene content, gene arrangment and gene structure shed light on the polyphyletic evolution of photoautotrophic orchid plastomes to mycoheterotrophic orchid plastomes.

DISCUSSION

Thus, this study's decoding of the mycoheterotrophic G. lindleyana plastome provides valuable resource data for future research and conservation of endangered orchids.

Effects of Steam and Water Blanching on Drying Characteristics, Water Distribution, Microstructure, and Bioactive Components of Gastrodia Elata

In the current work, the effects of steam and boiling water blanching on the drying characteristics, water distribution, microstructure, and contents of bioactive substances of Gastrodia elata (G. elata) were explored. Results showed that the degree of steaming and blanching was related to the core temperature of G. elata. The steaming and blanching pretreatment increased the drying time of the samples by more than 50%. The low-field nuclear magnetic resonance (LF-NMR) of treated samples showed that the relaxation time corresponded to water molecule states (bound, immobilized, and free) and G. elata became shorter, which indicated a reduction in free moisture and increased resistance of water diffusion in the solid structure during drying. Hydrolysis of polysaccharides and gelatinization of starch granules was observed in the microstructure of treated samples, which was consistent with changes in water status and drying rates. Steaming and blanching increased gastrodin and crude polysaccharide contents and decreased p-hydroxybenzyl alcohol content. These findings will contribute to a better understanding of the effect of steaming and blanching on the drying behavior and quality attributes of G. elata.

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.

Gnawing pressure led to the expansion of JAZ genes in angiosperms

A long-standing problem in evolutionary biology is why some populations differentiate into many species while the majority do not. Angiosperms is an excellent group for investigating this problem because their diversity is unevenly distributed in space and phylogeny. Plant hormone participates in growth, development and defense. However, jasmonic acid (JA) was the only hormone response to bites. We first searched jasmonate ZIM-domain (JAZ), AUXIN/INDOLE ACETIC ACID (IAA / aux), PYR/PYL/RCAR (PYL), DELLA, and SUPPRESSOR OF MAX2 1-like (SMAX) in 272 plant species. We found the gene number change trends were consistent with origination rates and species numbers of angiosperms. So, 26 representative species were selected as an example for further analysis. The results showed JAZ had experienced two lineage-specific gene expansion events in angiosperms, which coincided with increases in mammalian body size and dental diversity. The proliferation of large herbivores as a results of mammalian prosperity after dinosaur extinction may be related to angiosperm evolution and bursting. The proliferation of large herbivores as the result of mammalian prosperity after the extinction of the dinosaurs was related to angiosperm evolution and bursting. Overall, our study uncovered a previously unknown co-evolution mechanism in terrestrial plants exposed to extreme environmental conditions.

Multichromosomal Mitochondrial Genome of Paphiopedilum micranthum: Compact and Fragmented Genome, and Rampant Intracellular Gene Transfer

Orchidaceae is one of the largest families of angiosperms. Considering the large number of species in this family and its symbiotic relationship with fungi, Orchidaceae provide an ideal model to study the evolution of plant mitogenomes. However, to date, there is only one draft mitochondrial genome of this family available. Here, we present a fully assembled and annotated sequence of the mitochondrial genome (mitogenome) of Paphiopedilum micranthum, a species with high economic and ornamental value. The mitogenome of P. micranthum was 447,368 bp in length and comprised 26 circular subgenomes ranging in size from 5973 bp to 32,281 bp. The genome encoded for 39 mitochondrial-origin, protein-coding genes; 16 tRNAs (three of plastome origin); three rRNAs; and 16 ORFs, while rpl10 and sdh3 were lost from the mitogenome. Moreover, interorganellar DNA transfer was identified in 14 of the 26 chromosomes. These plastid-derived DNA fragments represented 28.32% (46,273 bp) of the P. micranthum plastome, including 12 intact plastome origin genes. Remarkably, the mitogenome of P. micranthum and Gastrodia elata shared 18% (about 81 kb) of their mitochondrial DNA sequences. Additionally, we found a positive correlation between repeat length and recombination frequency. The mitogenome of P. micranthum had more compact and fragmented chromosomes compared to other species with multichromosomal structures. We suggest that repeat-mediated homologous recombination enables the dynamic structure of mitochondrial genomes in Orchidaceae.

Evolutionary Patterns of the Chloroplast Genome in Vanilloid Orchids (Vanilloideae, Orchidaceae)

The Vanilloideae (vanilloids) is one of five subfamilies of Orchidaceae and is composed of fourteen genera and approximately 245 species. In this study, the six new chloroplast genomes (plastomes) of vanilloids (two Lecanorchis, two Pogonia, and two Vanilla species) were decoded, and then the evolutionary patterns of plastomes were compared to all available vanilloid plastomes. Pogonia japonica has the longest plastome, with 158,200 bp in genome size. In contrast, Lecanorchis japonica has the shortest plastome with 70,498 bp in genome size. The vanilloid plastomes have regular quadripartite structures, but the small single copy (SSC) region was drastically reduced. Two different tribes of Vanilloideae (Pogonieae and Vanilleae) showed different levels of SSC reductions. In addition, various gene losses were observed among the vanilloid plastomes. The photosynthetic vanilloids (Pogonia and Vanilla) showed signs of stage 1 degradation and had lost most of their ndh genes. The other three species (one Cyrotsia and two Lecanorchis), however, had stage 3 or stage 4 degradation and had lost almost all the genes in their plastomes, except for some housekeeping genes. The Vanilloideae were located between the Apostasioideae and Cypripedioideae in the maximum likelihood tree. A total of ten rearrangements were found among ten Vanilloideae plastomes when compared to the basal Apostasioideae plastomes. The four sub-regions of the single copy (SC) region shifted into an inverted repeat (IR) region, and the other four sub-regions of the IR region shifted into the SC regions. Both the synonymous (dS) and nonsynonymous (dN) substitution rates of IR in-cooperated SC sub-regions were decelerated, while the substitution rates of SC in-cooperated IR sub-regions were accelerated. A total of 20 protein-coding genes remained in mycoheterotrophic vanilloids. Almost all these protein genes show accelerated base substitution rates compared to the photosynthetic vanilloids. Two of the twenty genes in the mycoheterotrophic species faced strong "relaxed selection" pressure (p-value < 0.05).

Comparative transcriptome analysis of Armillaria gallica 012m in response to ethephon treatment

Background

Gastrodia elata, known as a rootless, leafless, achlorophyllous and fully mycoheterotrophic orchid, needs to establish symbionts with particular Armillaria species to acquire nutrition and energy. Previous research findings had approved that ethylene (ET) played an important role in plant-fungi interaction and some receptors of ET had been discovered in microorganisms. However, the molecular mechanisms underlying the role of ET in the interaction between G. elata and Armillaria species remain unknown.

Methods

Exiguous ethephon (ETH) was added to agar and liquid media to observe the morphological features of mycelium and count the biomass respectively. Mycelium cultured in liquid media with exiguous ETH (0.1 ppm, 2.0 ppm, 5.0 ppm) were chosen to perform whole-transcriptome profiling through the RNA-seq technology (Illumina NGS sequencing). The DEGs of growth-related genes and candidate ET receptor domains were predicted on SMART.

Results

ETH-0.1 ppm and ETH-2 ppm could significantly improve the mycelium growth of A. gallica 012m, while ETH-5 ppm inhibited the mycelium growth in both solid and liquid media. The number of up-regulated or down-regulated genes increased along with the concentrations of ETH. The growth of mycelia might benefit from the up-regulated expression of Pyr_redox (Pyridine nucleotide-disulphide oxidoreductase), GAL4 (C6 zinc finger) and HMG (High Mobility Group) genes in the ETH-0.1 ppm and ETH-2 ppm. Therefore, the growth of mycelia might be impaired by the down-regulated expression of ZnF_C2H2 and ribosomal protein S4 proteins in the ETH-5 ppm. Seven ET receptor domains were predicted in A. gallica 012m. Based on cluster analysis and comparative studies of proteins, the putative ETH receptor domains of A. gallica 012m have a higher homologous correlation with fungi.

Conclusions

The responses of A. gallica 012m to ETH had a concentration effect similar to the plants’ responses to ET. Therefore, the number of up-regulated or down-regulated genes are increased along with the concentrations of ETH. Seven ET receptor protein domains were predicted in the genome and transcriptome of A. gallica 012m. We speculate that ETH receptors exist in A. gallica 012m and ethylene might play an important role in the plant-fungi interaction.

Genome sequencing and resequencing identified three horizontal gene transfers and uncovered the genetic mechanism on the intraspecies adaptive evolution of Gastrodia elata Blume

Horizontal gene transfer is a rare and useful genetic mechanism in higher plants. Gastrodia elata Blume (GE) (Orchidaceae), well known as traditional medicinal material in East Asia, adopts a heterotrophic lifestyle, thus being considered to be more prone to horizontal gene transfer (HGT). GE is a "polytypic species" that currently comprised of five recognized forms according to the plant morphology. G. elata Blume forma elata (GEE) and G. elata Bl.f.glauca (GEG) are two common forms that naturally grow in different habitats with difference in altitude and latitude. G. elata Bl.f.viridis (GEV) often occurs sporadically in cultivated populations of GEE and GEG. However, the genetic relationships and genetic mechanism underpinned the divergent ecological adaptations of GEE and GEG have not been revealed. Here, we assembled a chromosome-level draft genome of GEE with 1.04 Gb. Among predicted 17,895 protein coding genes, we identified three HGTs. Meanwhile, we resequenced 10 GEE accessions, nine GEG accessions, and 10 GEV accessions, and identified two independent genetic lineages: GEG_pedigree (GEG individuals and GEV individuals collected from GEG populations) and GEE_pedigree (GEE individuals and GEV individuals collected from GEE populations), which strongly support the taxonomic status of GEE and GEG as subspecies, not as different forms. In highly differentiated genomic regions of GEE_pedigree and GEG_pedigree, three chalcone synthase-encoding genes and one Phox/Bem1p (PB1) domain of encoding Auxin (AUX)/Indoleacetic acid (IAA) were identified in selection sweeping genome regions, which suggested that differentiation between GEE_pedigree and GEG_pedigree was promoted by the selection of genes related to photoresponse and growth and development. Overall, this new genome would be helpful for breeding and utilization of GE and the new findings would deepen the understanding about ecological adaptation and evolution of GE.

Transcriptome analysis of Armillaria gallica 012 m in response to auxin

Gastrodia elata is an achlorophyllous and fully mycoheterotrophic orchid which obtains carbon and other nutrients from Armillaria species in its life cycle. Many researchers suggested that plant hormones, as signing molecules, play a central role in the plant-fungi interaction. In the process of Armillaria gallica 012 m cultivation, both exogenous indole-3-acetic acid (IAA) and indole-3-butyric acid (IBA) distinctly stimulated the growth of mycelia in solid media. The differential expression genes (DEGs) of A. gallica 012 m with IAA versus blank control (BK) and IBA versus BK were investigated. The results showed that more than 80% of DEGs of the IAA group were coincident with the DEGs of the IBA group, and more than half of upregulated DEGs and most of the downregulated DEGs of the IAA group coincided with those DEGs of the IBA group. Above research implied that A. gallica 012 m could perceive IAA and IBA, and possess similar responses and signaling pathways to IAA and IBA. The overlapping differential genes of the IAA group and IBA group were analyzed by GO term, and the results showed that several DEGs identified were related to biological processes including positive regulation of the biological process and biological process. The downregulated NmrA-like and FKBP_C genes might be benefit to the growth of mycelia. Those results can explain that exiguous IAA and IBA improved the growth of A. gallica to some extent. We speculate that IAA and IBA are signaling molecules, and regulate the expression of growth-related genes of A. gallica 012 m by the same signaling pathway.

Transcriptomic heterochrony and completely cleistogamous flower development in the mycoheterotrophic orchid Gastrodia

Cleistogamy, in which plants can reproduce via self-fertilization within permanently closed flowers, has evolved in > 30 angiosperm lineages; however, consistent with Darwin's doubts about its existence, complete cleistogamy - the production of only cleistogamous flowers - has rarely been recognized. Thus far, the achlorophyllous orchid genus, Gastrodia, is the only known genus with several plausible completely cleistogamous species. Here, we analyzed the floral developmental transcriptomes of two recently evolved, completely cleistogamous Gastrodia species and their chasmogamous sister species to elucidate the possible changes involved in producing common cleistogamous traits. The ABBA-BABA test did not support introgression and protein sequence convergence as evolutionary mechanisms leading to cleistogamy, leaving convergence in gene expression as a plausible mechanism. Regarding transcriptomic differentiation, the two cleistogamous species had common modifications in the expression of developmental regulators, exhibiting a gene family-wide signature of convergent expression changes in MADS-box genes. Our transcriptomic pseudotime analysis revealed a prolonged juvenile state and eventual maturation, a heterochronic pattern consistent with partial neoteny, in cleistogamous flower development. These findings indicate that transcriptomic partial neoteny, arising from changes in the expression of conserved developmental regulators, might have contributed to the rapid and repeated evolution of cleistogamous flowers in Gastrodia.

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.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Genome-wide analysis of the TCP gene family and their expression pattern in Cymbidium goeringii

TCP gene family are specific transcription factors for plant, and considered to play an important role in development and growth. However, few related studies investigated the TCP gene trait and how it plays a role in growth and development of Orchidaceae. In this study, we obtained 14 TCP genes (CgTCPs) from the Spring Orchid Cymbidium goeringii genome. The classification results showed that 14 CgTCPs were mainly divided into two clades as follows: four PCF genes (Class I), nine CIN genes and one CYC gene (Class II). The sequence analysis showed that the TCP proteins of C. goeringii contain four conserved regions (basic Helix-Loop-Helix) in the TCP domain. The exon-intron structure varied in the clade according to a comparative investigation of the gene structure, and some genes had no introns. There are fewer CgTCP homologous gene pairs compared with Dendrobium catenatum and Phalaenopsis equestris, suggesting that the TCP genes in C. goeringii suffered more loss events. The majority of the cis-elements revealed to be enriched in the function of light responsiveness, followed by MeJA and ABA responsiveness, demonstrating their functions in regulating by light and phytohormones. The collinearity study revealed that the TCPs in D. catenatum, P. equestris and C. goeringii almost 1:1. The transcriptomic data and real-time reverse transcription-quantitative PCR (RT-qPCR) expression profiles showed that the flower-specific expression of the TCP class II genes (CgCIN2, CgCIN5 and CgCIN6) may be related to the regulation of florescence. Altogether, this study provides a comprehensive analysis uncovering the underlying function of TCP genes in Orchidaceae.

OrchidBase 5.0: updates of the orchid genome knowledgebase

Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii

PREMISE

Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution.

METHODS

We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome.

RESULTS

All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis.

CONCLUSIONS

Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

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.

Advances and prospects of orchid research and industrialization

Orchidaceae is one of the largest, most diverse families in angiosperms with significant ecological and economical values. Orchids have long fascinated scientists by their complex life histories, exquisite floral morphology and pollination syndromes that exhibit exclusive specializations, more than any other plants on Earth. These intrinsic factors together with human influences also make it a keystone group in biodiversity conservation. The advent of sequencing technologies and transgenic techniques represents a quantum leap in orchid research, enabling molecular approaches to be employed to resolve the historically interesting puzzles in orchid basic and applied biology. To date, 16 different orchid genomes covering four subfamilies (Apostasioideae, Vanilloideae, Epidendroideae, and Orchidoideae) have been released. These genome projects have given rise to massive data that greatly empowers the studies pertaining to key innovations and evolutionary mechanisms for the breadth of orchid species. The extensive exploration of transcriptomics, comparative genomics, and recent advances in gene engineering have linked important traits of orchids with a multiplicity of gene families and their regulating networks, providing great potential for genetic enhancement and improvement. In this review, we summarize the progress and achievement in fundamental research and industrialized application of orchids with a particular focus on molecular tools, and make future prospects of orchid molecular breeding and post-genomic research, providing a comprehensive assemblage of state of the art knowledge in orchid research and industrialization.

In-depth analysis of genomes and functional genomics of orchid using cutting-edge high-throughput sequencing

High-throughput sequencing technology has been facilitated the development of new methodologies and approaches for studying the origin and evolution of plant genomes and subgenomes, population domestication, and functional genomics. Orchids have tens of thousands of members in nature. Many of them have promising application potential in the extension and conservation of the ecological chain, the horticultural use of ornamental blossoms, and the utilization of botanical medicines. However, a large-scale gene knockout mutant library and a sophisticated genetic transformation system are still lacking in the improvement of orchid germplasm resources. New gene editing tools, such as the favored CRISPR-Cas9 or some base editors, have not yet been widely applied in orchids. In addition to a large variety of orchid cultivars, the high-precision, high-throughput genome sequencing technology is also required for the mining of trait-related functional genes. Nowadays, the focus of orchid genomics research has been directed to the origin and classification of species, genome evolution and deletion, gene duplication and chromosomal polyploidy, and flower morphogenesis-related regulation. Here, the progressing achieved in orchid molecular biology and genomics over the past few decades have been discussed, including the evolution of genome size and polyploidization. The frequent incorporation of LTR retrotransposons play important role in the expansion and structural variation of the orchid genome. The large-scale gene duplication event of the nuclear genome generated plenty of recently tandem duplicated genes, which drove the evolution and functional divergency of new genes. The evolution and loss of the plastid genome, which mostly affected genes related to photosynthesis and autotrophy, demonstrated that orchids have experienced more separate transitions to heterotrophy than any other terrestrial plant. Moreover, large-scale resequencing provide useful SNP markers for constructing genetic maps, which will facilitate the breeding of novel orchid varieties. The significance of high-throughput sequencing and gene editing technologies in the identification and molecular breeding of the trait-related genes in orchids provides us with a representative trait-improving gene as well as some mechanisms worthy of further investigation. In addition, gene editing has promise for the improvement of orchid genetic transformation and the investigation of gene function. This knowledge may provide a scientific reference and theoretical basis for orchid genome studies.