A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants

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.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

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.

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.

Autoactivation of mycorrhizal symbiosis signaling through gibberellin deactivation in orchid seed germination

Orchids parasitically depend on external nutrients from mycorrhizal fungi for seed germination. Previous findings suggest that orchids utilize a genetic system of mutualistic arbuscular mycorrhizal (AM) symbiosis, in which the plant hormone gibberellin (GA) negatively affects fungal colonization and development, to establish parasitic symbiosis. Although GA generally promotes seed germination in photosynthetic plants, previous studies have reported low sensitivity of GA in seed germination of mycoheterotrophic orchids where mycorrhizal symbiosis occurs concurrently. To elucidate the connecting mechanisms of orchid seed germination and mycorrhizal symbiosis at the molecular level, we investigated the effect of GA on a hyacinth orchid (Bletilla striata) seed germination and mycorrhizal symbiosis using asymbiotic and symbiotic germination methods. Additionally, we compared the transcriptome profiles between asymbiotically and symbiotically germinated seeds. Exogenous GA negatively affected seed germination and fungal colonization, and endogenous bioactive GA was actively converted to the inactive form during seed germination. Transcriptome analysis showed that B. striata shared many of the induced genes between asymbiotically and symbiotically germinated seeds, including GA metabolism- and signaling-related genes and AM-specific marker homologs. Our study suggests that orchids have evolved in a manner that they do not use bioactive GA as a positive regulator of seed germination and instead autoactivate the mycorrhizal symbiosis pathway through GA inactivation to accept the fungal partner immediately during seed germination.

Selection on Visual Opsin Genes in Diurnal Neotropical Frogs and Loss of the SWS2 Opsin in Poison Frogs

Amphibians are ideal for studying visual system evolution because their biphasic (aquatic and terrestrial) life history and ecological diversity expose them to a broad range of visual conditions. Here, we evaluate signatures of selection on visual opsin genes across Neotropical anurans and focus on three diurnal clades that are well-known for the concurrence of conspicuous colors and chemical defense (i.e., aposematism): poison frogs (Dendrobatidae), Harlequin toads (Bufonidae: Atelopus), and pumpkin toadlets (Brachycephalidae: Brachycephalus). We found evidence of positive selection on 44 amino acid sites in LWS, SWS1, SWS2, and RH1 opsin genes, of which one in LWS and two in RH1 have been previously identified as spectral tuning sites in other vertebrates. Given that anurans have mostly nocturnal habits, the patterns of selection revealed new sites that might be important in spectral tuning for frogs, potentially for adaptation to diurnal habits and for color-based intraspecific communication. Furthermore, we provide evidence that SWS2, normally expressed in rod cells in frogs and some salamanders, has likely been lost in the ancestor of Dendrobatidae, suggesting that under low-light levels, dendrobatids have inferior wavelength discrimination compared to other frogs. This loss might follow the origin of diurnal activity in dendrobatids and could have implications for their behavior. Our analyses show that assessments of opsin diversification in across taxa could expand our understanding of the role of sensory system evolution in ecological adaptation.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

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.

Impact of climate change on the geographical distribution and niche dynamics of Gastrodia elata

Background

Gastrodia elata is widely used in China as a valuable herbal medicine. Owing to its high medicinal and nutrient value, wild resources of G. elata have been overexploited and its native areas have been severely damaged. Understanding the impacts of climate change on the distribution of this endangered species is important for the conservation and sustainable use of G. elata.

Methods

We used the optimized maximum entropy model to simulate the potential distribution of G. elata under contemporary and future time periods (1970-2000, 2050s, 2070s, and 2090s) and different climate change scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5). Under these conditions, we investigated the key environmental factors influencing the distribution of G. elata as well as the spatial and temporal characteristics of its niche dynamics.

Results

With high Maxent model accuracy (AUCmean = 0.947 ± 0.012, and the Kappa value is 0.817), our analysis revealed that annual precipitation, altitude, and mean temperature of driest quarter are the most important environmental factors influencing the distribution of G. elata. Under current bioclimatic conditions, the potentially suitable area for G. elata in China is 71.98 × 10⁴ km², while the highly suitable region for G. elata growth is 7.28 × 10⁴ km². Our models for three future periods under four climate change scenarios indicate that G. elata can maintain stable distributions in southern Shaanxi, southwestern Hubei, and around the Sichuan basin, as these areas are highly suitable for its growth. However, the center of the highly suitable areas of G. elata shift depending on different climatic scenarios. The values of niche overlap for G. elata show a decreasing trend over the forecasted periods, of which the niche overlap under the SSP3-7.0 scenario shows the greatest decrease.

Discussions

Under the condition of global climate change in the future, our study provides basic reference data for the conservation and sustainable utilization of the valuable and endangered medicinal plant G. elata. It is important to carefully choose the protection area of G. elata wild resources according the suitable area conditions modeled. Moreover, these findings will be valuable for providing insights into the breeding and artificial cultivation of this plant, including the selection of suitable areas for planting.

The Sapria himalayana genome provides new insights into the lifestyle of endoparasitic plants

BACKGROUND

Sapria himalayana (Rafflesiaceae) is an endoparasitic plant characterized by a greatly reduced vegetative body and giant flowers; however, the mechanisms underlying its special lifestyle and greatly altered plant form remain unknown. To illustrate the evolution and adaptation of S. himalayasna, we report its de novo assembled genome and key insights into the molecular basis of its floral development, flowering time, fatty acid biosynthesis, and defense responses.

RESULTS

The genome of S. himalayana is ~ 1.92 Gb with 13,670 protein-coding genes, indicating remarkable gene loss (~ 54%), especially genes involved in photosynthesis, plant body, nutrients, and defense response. Genes specifying floral organ identity and controlling organ size were identified in S. himalayana and Rafflesia cantleyi, and showed analogous spatiotemporal expression patterns in both plant species. Although the plastid genome had been lost, plastids likely biosynthesize essential fatty acids and amino acids (aromatic amino acids and lysine). A set of credible and functional horizontal gene transfer (HGT) events (involving genes and mRNAs) were identified in the nuclear and mitochondrial genomes of S. himalayana, most of which were under purifying selection. Convergent HGTs in Cuscuta, Orobanchaceae, and S. himalayana were mainly expressed at the parasite-host interface. Together, these results suggest that HGTs act as a bridge between the parasite and host, assisting the parasite in acquiring nutrients from the host.

CONCLUSIONS

Our results provide new insights into the flower development process and endoparasitic lifestyle of Rafflesiaceae plants. The amount of gene loss in S. himalayana is consistent with the degree of reduction in its body plan. HGT events are common among endoparasites and play an important role in their lifestyle adaptation.

Habitat-related plastome evolution in the mycoheterotrophic Neottia listeroides complex (Orchidaceae, Neottieae)

BACKGROUND

Mycoheterotrophs, acquiring organic carbon and other nutrients from mycorrhizal fungi, have evolved repeatedly with substantial plastid genome (plastome) variations. To date, the fine-scale evolution of mycoheterotrophic plastomes at the intraspecific level is not well-characterized. A few studies have revealed unexpected plastome divergence among species complex members, possibly driven by various biotic/abiotic factors. To illustrate evolutionary mechanisms underlying such divergence, we analyzed plastome features and molecular evolution of 15 plastomes of Neottia listeroides complex from different forest habitats.

RESULTS

These 15 samples of Neottia listeroides complex split into three clades according to their habitats approximately 6 million years ago: Pine Clade, including ten samples from pine-broadleaf mixed forests, Fir Clade, including four samples from alpine fir forests and Fir-willow Clade with one sample. Compared with those of Pine Clade members, plastomes of Fir Clade members show smaller size and higher substitution rates. Plastome size, substitution rates, loss and retention of plastid-encoded genes are clade-specific. We propose to recognized six species in N. listeroides complex and slightly modify the path of plastome degradation.

CONCLUSIONS

Our results provide insight into the evolutionary dynamics and discrepancy of closely related mycoheterotrophic orchid lineages at a high phylogenetic resolution.

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.

ZmMYC2s play important roles in maize responses to simulated herbivory and jasmonate

Both herbivory and jasmonic acid (JA) activate the biosynthesis of defensive metabolites in maize, but the mechanism underlying this remains unclear. We generated maize mutants in which ZmMYC2a and ZmMYC2b, two transcription factor genes important in JA signaling, were individually or both knocked out. Genetic and biochemical analyses were used to elucidate the functions of ZmMYC2 proteins in the maize response to simulated herbivory and JA. Compared with the wild-type (WT) maize, the double mutant myc2ab was highly susceptible to insects, and the levels of benzoxazinoids and volatile terpenes, and the levels of their biosynthesis gene transcripts, were much lower in the mutants than in the WT maize after simulated insect feeding or JA treatment. Moreover, ZmMYC2a and ZmMYC2b played a redundant role in maize resistance to insects and JA signaling. Transcriptome and Cleavage Under Targets and Tagmentation-Sequencing (CUT&Tag-Seq) analysis indicated that ZmMYC2s physically targeted 60% of the JA-responsive genes, even though only 33% of these genes were transcriptionally ZmMYC2-dependent. Importantly, CUT&Tag-Seq and dual luciferase assays revealed that ZmMYC2s transactivate the benzoxazinoid and volatile terpene biosynthesis genes IGPS1/3, BX10/11/12/14, and TPS10/2/3/4/5/8 by directly binding to their promoters. Furthermore, several transcription factors physically targeted by ZmMYC2s were identified, and these are likely to function in the regulation of benzoxazinoid biosynthesis. This work reveals the transcriptional regulatory landscapes of both JA signaling and ZmMYC2s in maize and provides comprehensive mechanistic insight into how JA signaling modulates defenses in maize responses to herbivory through ZmMYC2s.

A proposed metric set for evaluation of genome assembly quality

Quality control is essential for genome assemblies; however, a consensus has yet to be reached on what metrics should be adopted for the evaluation of assembly quality. N50 is widely used for contiguity measurement, but its effectiveness is constantly in question. Prevailing metrics for the completeness evaluation focus on gene space, yet challenging areas such as tandem repeats are commonly overlooked. Achieving correctness has become an indispensable dimension for quality control, while prevailing assembly releases lack scores reflecting this aspect. We propose a metric set with a set of statistic indexes for effective, comprehensive evaluation of assemblies and provide a score of a finished assembly for each metric, which can be utilized as a benchmark for achieving high-quality genome assemblies.

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.

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.

Genome-wide identification and adaptive evolution of CesA/Csl superfamily among species with different life forms in Orchidaceae

Orchidaceae, with more than 25,000 species, is one of the largest flowering plant families that can successfully colonize wide ecological niches, such as land, trees, or rocks, and its members are divided into epiphytic, terrestrial, and saprophytic types according to their life forms. Cellulose synthase (CesA) and cellulose synthase-like (Csl) genes are key regulators in the synthesis of plant cell wall polysaccharides, which play an important role in the adaptation of orchids to resist abiotic stresses, such as drought and cold. In this study, nine whole-genome sequenced orchid species with three types of life forms were selected; the CesA/Csl gene family was identified; the evolutionary roles and expression patterns of CesA/Csl genes adapted to different life forms and abiotic stresses were investigated. The CesA/Csl genes of nine orchid species were divided into eight subfamilies: CesA and CslA/B/C/D/E/G/H, among which the CslD subfamily had the highest number of genes, followed by CesA, whereas CslB subfamily had the least number of genes. Expansion of the CesA/Csl gene family in orchids mainly occurred in the CslD and CslF subfamilies. Conserved domain analysis revealed that eight subfamilies were conserved with variations in orchids. In total, 17 pairs of CesA/Csl homologous genes underwent positive selection, of which 86%, 14%, and none belonged to the epiphytic, terrestrial, and saprophytic orchids, respectively. The inter-species collinearity analysis showed that the CslD genes expanded in epiphytic orchids. Compared with terrestrial and saprophytic orchids, epiphytic orchids experienced greater strength of positive selection, with expansion events mostly related to the CslD subfamily, which might have resulted in strong adaptability to stress in epiphytes. Experiments on stem expression changes under abiotic stress showed that the CslA might be a key subfamily in response to drought stress for orchids with different life forms, whereas the CslD might be a key subfamily in epiphytic and saprophytic orchids to adapt to freezing stress. This study provides the basic knowledge for the further systematic study of the adaptive evolution of the CesA/Csl superfamily in angiosperms with different life forms, and research on orchid-specific functional genes related to life-history trait evolution.

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.

The extremely reduced, diverged and reconfigured plastomes of the largest mycoheterotrophic orchid lineage

BACKGROUND

Plastomes of heterotrophic plants have been greatly altered in structure and gene content, owing to the relaxation of selection on photosynthesis-related genes. The orchid tribe Gastrodieae is the largest and probably the oldest mycoheterotrophic clade of the extant family Orchidaceae. To characterize plastome evolution across members of this key important mycoheterotrophic lineage, we sequenced and analyzed the plastomes of eleven Gastrodieae members, including representative species of two genera, as well as members of the sister group Nervilieae.

RESULTS

The plastomes of Gastrodieae members contain 20 protein-coding, four rRNA and five tRNA genes. Evolutionary analysis indicated that all rrn genes were transferred laterally and together, forming an rrn block in the plastomes of Gastrodieae. The plastome GC content of Gastrodia species ranged from 23.10% (G. flexistyla) to 25.79% (G. javanica). The plastome of Didymoplexis pallens contains two copies each of ycf1 and ycf2. The synonymous and nonsynonymous substitution rates were very high in the plastomes of Gastrodieae among mycoheterotrophic species in Orchidaceae and varied between genes.

CONCLUSIONS

The plastomes of Gastrodieae are greatly reduced and characterized by low GC content, rrn block formation, lineage-specific reconfiguration and gene content, which might be positively selected. Overall, the plastomes of Gastrodieae not only serve as an excellent model for illustrating the evolution of plastomes but also provide new insights into plastome evolution in parasitic plants.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

Comparative genomics of orobanchaceous species with different parasitic lifestyles reveals the origin and stepwise evolution of plant parasitism

Orobanchaceae is the largest family of parasitic plants, containing autotrophic and parasitic plants with all degrees of parasitism. This makes it by far the best family for studying the origin and evolution of plant parasitism. Here we provide three high-quality genomes of orobanchaceous plants, the autotrophic Lindenbergia luchunensis and the holoparasitic plants Phelipanche aegyptiaca and Orobanche cumana. Phylogenomic analysis of these three genomes together with those previously published and the transcriptomes of other orobanchaceous species created a robust phylogenetic framework for Orobanchaceae. We found that an ancient whole-genome duplication (WGD; about 73.48 million years ago), which occurred earlier than the origin of Orobanchaceae, might have contributed to the emergence of parasitism. However, no WGD events occurred in any lineage of orobanchaceous parasites except for Striga after divergence from their autotrophic common ancestor, suggesting that, in contrast with previous speculations, WGD is not associated with the emergence of holoparasitism. We detected evident convergent gene loss in all parasites within Orobanchaceae and between Orobanchaceae and dodder Cuscuta australis. The gene families in the orobanchaceous parasites showed a clear pattern of recent gains and expansions. The expanded gene families are enriched in functions related to the development of the haustorium, suggesting that recent gene family expansions may have facilitated the adaptation of orobanchaceous parasites to different hosts. This study illustrates a stepwise pattern in the evolution of parasitism in the orobanchaceous parasites and will facilitate future studies on parasitism and the control of parasitic plants in agriculture.

Systematics and Plastome Evolution in Schizaeaceae

While the family Schizaeaceae (Schizaeales) represents only about 0.4% of the extant fern species diversity, it differs from other ferns greatly in gross morphologies, niche preferences, and life histories. One of the most notable features in this family is its mycoheterotrophic life style in the gametophytic stage, which appears to be associated with extensive losses of plastid genes. However, the limited number of sequenced plastomes, and the lack of a well-resolved phylogenetic framework of Schizaeaceae, makes it difficult to gain any further insight. Here, with a comprehensive sampling of ~77% of the species diversity of this family, we first inferred a plastid phylogeny of Schizaeaceae using three DNA regions. To resolve the deep relationships within this family, we then reconstructed a plastome-based phylogeny focusing on a selection of representatives that covered all the major clades. From this phylogenomic backbone, we traced the evolutionary histories of plastid genes and examined whether gene losses were associated with the evolution of gametophytic mycoheterotrophy. Our results reveal that extant Schizaeaceae is comprised of four major clades-Microschizaea, Actinostachys, Schizaea, and Schizaea pusilla. The loss of all plastid NADH-like dehydrogenase (ndh) genes was confirmed to have occurred in the ancestor of extant Schizaeaceae, which coincides with the evolution of mycoheterotrophy in this family. For chlorophyll biosynthesis genes (chl), the losses were interpreted as convergent in Schizaeaceae, and found not only in Actinostachys, a clade producing achlorophyllous gametophytes, but also in S. pusilla with chlorophyllous gametophytes. In addition, we discovered a previously undescribed but phylogenetically distinct species hidden in the Schizaea dichotoma complex and provided a taxonomic treatment and morphological diagnostics for this new species-Schizaea medusa. Finally, our phylogenetic results suggest that the current PPG I circumscription of Schizaea is non-monophyletic, and we therefore proposed a three-genus classification moving a subset of Schizaea species sensu PPG I to a third genus-Microschizaea.

The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions

BACKGROUND

To illustrate the molecular mechanism of mycoheterotrophic interactions between orchids and fungi, we assembled chromosome-level reference genome of Gastrodia menghaiensis (Orchidaceae) and analyzed the genomes of two species of Gastrodia.

RESULTS

Our analyses indicated that the genomes of Gastrodia are globally diminished in comparison to autotrophic orchids, even compared to Cuscuta (a plant parasite). Genes involved in arbuscular mycorrhizae colonization were found in genomes of Gastrodia, and many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids. The highly expressed genes for fatty acid and ammonium root transporters suggest that fungi receive material from orchids, although most raw materials flow from the fungi. Many nuclear genes (e.g. biosynthesis of aromatic amino acid L-tryptophan) supporting plastid functions are expanded compared to photosynthetic orchids, an indication of the importance of plastids even in totally mycoheterotrophic species.

CONCLUSION

Gastrodia menghaiensis has the smallest proteome thus far among angiosperms. Many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

Genome-wide identification of Cymbidium sinense WRKY gene family and the importance of its Group III members in response to abiotic stress

Transcription factors (TFs) of the WRKY family play pivotal roles in defense responses and secondary metabolism of plants. Although WRKY TFs are well documented in numerous plant species, no study has performed a genome-wide investigation of the WRKY gene family in Cymbidium sinense. In the present work, we found 64 C. sinense WRKY (CsWRKY) TFs, and they were further divided into eight subgroups. Chromosomal distribution of CsWRKYs revealed that the majority of these genes were localized on 16 chromosomes, especially on Chromosome 2. Syntenic analysis implied that 13 (20.31%) genes were derived from segmental duplication events, and 17 orthologous gene pairs were identified between Arabidopsis thaliana WRKY (AtWRKY) and CsWRKY genes. Moreover, 55 of the 64 CsWRKYs were detectable in different plant tissues in response to exposure to plant hormones. Among them, Group III members were strongly induced in response to various hormone treatments, indicating their potential essential roles in hormone signaling. We subsequently analyzed the function of CsWRKY18 in Group III. The CsWRKY18 was localized in the nucleus. The constitutive expression of CsWRKY18 in Arabidopsis led to enhanced sensitivity to ABA-mediated seed germination and root growth and elevated plant tolerance to abiotic stress within the ABA-dependent pathway. Overall, our study represented the first genome-wide characterization and functional analysis of WRKY TFs in C. sinense, which could provide useful clues about the evolution and functional description of CsWRKY genes.