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

Comparative transcriptome analysis reveals the role of sugar signaling in response to high temperature stress in Armillaria gallica

BACKGROUND

Armillaria gallica establish a unique symbiotic relationship with Gastrodia elata, which is a valuable Chinese herbal plant that relies on this symbiosis. High temperature (HT) seriously affected the growth of both A. gallica and G. elata. However, the molecular mechanisms and gene networks involved in the response of A. gallica to HT are not well understood.

RESULTS

In this study, transcriptomic analyses at multiple time point were conducted using rhizomorph from two A. gallica strains, GZ1 (tolerance to HT) and SX8 (sensitive to HT) under HT and normal temperature (NT) conditions. We identified 2,056 differentially expressed genes, primarily associated with carbohydrate metabolism and protein processing in endoplasmic reticulum (ER). Heat shock proteins (HSP20 and HSP90) were up-regulated in both GZ1 and SX8 under HT, indicating a conserved HT induced response in A. gallica. Additionally, genes related to glycolysis were up-regulated in GZ1 following HT treatment, while genes involved in the conversion of sugar to amino acids were down-regulated in GZ1 under HT treatment. GZ1 also showed an increased accumulation of soluble sugar and polysaccharide under HT, which were significantly correlated with pyruvate kinase and aldo/keto reductase genes. Furthermore, a sucrose concentration of 30 g/L enhanced HT resistance in A. gallica.

CONCLUSIONS

These results reveal that sugar signaling, particularly in carbohydrate metabolism, plays a critical role in the rhizomorph response to HT stress. Both conserved and species-specific transcriptome changes may contribute to the dynamic regulation of rhizomorph development during adaptation to HT stress in A. gallica.

CLINICAL TRIAL NUMBER

Not applicable.

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Centenary Progress on Orchidaceae Research: A Bibliometric Analysis

BACKGROUND

Research on orchids has experienced substantial growth since the early 20th century, reflecting their ecological and evolutionary significance.

METHODS

This paper provides a comprehensive bibliometric analysis of orchid-related literature published between 1902 and 2024, based on data retrieved from the Web of Science Core Collection™ (WoS).

RESULTS

The primary goal is to assess the global research landscape of orchids by identifying key authors, institutions, and journals, as well as major research themes in the field. A thorough analysis of publication trends, citation frequencies, and keyword co-occurrence networks was conducted to uncover significant research hotspots. The findings indicate that orchid research has evolved from foundational topics such as taxonomy and classification to more intricate subjects, including conservation strategies, orchid-pollinator dynamics, and the role of orchids in ecosystem functions. Additionally, biotechnology-related research is emerging as a dominant trend. This study also highlights that China has the highest publication output, while collaboration between the United States and Europe continues to grow. The co-word analysis of keywords suggests that future research is likely to continue to focus on orchid conservation, the impacts of climate change, pollination biology, and symbiotic relationships with mycorrhizal fungi.

CONCLUSIONS

This review offers valuable insights for researchers and conservationists, helping to identify future research priorities and strategies for the preservation and sustainable use of orchids.

Convergent reduction of immune receptor repertoires during plant adaptation to diverse special lifestyles and habitats

Plants deploy cell-surface pattern recognition receptors (PRRs) and intracellular nucleotide-binding site-leucine-rich repeat receptors (NLRs) to recognize pathogens. However, how plant immune receptor repertoires evolve in responding to changed pathogen burdens remains elusive. Here we reveal the convergent reduction of NLR repertoires in plants with diverse special lifestyles/habitats (SLHs) encountering low pathogen burdens. Furthermore, a parallel but milder reduction of PRR genes in SLH species was observed. The reduction of PRR and NLR genes was attributed to both increased gene loss and decreased gene duplication. Notably, pronounced loss of immune receptors was associated with the complete absence of signalling components from the enhanced disease susceptibility 1 (EDS1) and the resistance to powdery mildew 8 (RPW8)-NLR (RNL) families. In addition, evolutionary pattern analysis suggested that the conserved toll/interleukin-1 receptor (TIR)-only proteins might function tightly with EDS1/RNL. Taken together, these results reveal the hierarchically adaptive evolution of the two-tiered immune receptor repertoires during plant adaptation to diverse SLHs.

Establishment of an efficient Agrobacterium tumefaciens-mediated transformation system for an Armillaria species, a host of the fully mycoheterotrophic plant Gastrodia elata

The genus Armillaria (Basidiomycota, Agaricales, Physalacriaceae) comprises pathogenic fungi that cause root-rot disease in plants, as well as species with low pathogenicity, some of which are hosts of the fully mycoheterotrophic orchid plant Gastrodia elata (Orchidaceae). To investigate the mechanisms underlying such special interactions between Armillaria fungi and G. elata, it is crucial to establish genetic transformation platforms for the Armillaria fungi and G. elata. In this study, an Armillaria strain Arm37 was isolated from G. elata, which can form symbiosis with G. elata in axenic culture under laboratory conditions. A vector pYT-EV containing a cassette for hygromycin-resistance selection and a cassette for expressing or silencing target genes was constructed. An Agrobacterium tumefaciens-mediated transformation (ATMT) system for Arm37 was successfully developed and optimized to achieve a transformation efficiency of 32%. The ATMT system was successfully used to express the reporter genes eGFP encoding enhanced green fluorescent protein and GUS encoding β-glucuronidase and to effectively silence the endogenous gene URA3 encoding orotidine-5'-phosphate decarboxylase in Arm37. This ATMT system established for Arm37 provides an efficient genetic tool for exploring the Arm37 genes that are involved in the unique interaction between the Armillaria fungi and fully mycoheterotrophic plant G. elata.

Different roles of the phytohormone gibberellin in the wide-spread arbuscular mycorrhiza and in orchid mycorrhiza

Gibberellin (GA) is a classical plant hormone that regulates many physiological processes, such as plant growth, development, and environmental responses. GA inhibits arbuscular mycorrhizal (AM) symbiosis, the most ancient and widespread type of mycorrhizal symbiosis. Knowledge about mycorrhizal symbioses at the molecular level has been obtained mainly in model plants such as legumes and rice. In contrast, molecular mechanisms in non-model plants are still unclear. Recent studies have revealed the novel roles of GA in mycorrhizal symbioses: its positive effect in Paris-type AM symbiosis in Eustoma grandiflorum and its negative effect on both seed germination and mycorrhizal symbiosis in orchids. This review focuses on the recent data on GA function in AM and orchid mycorrhizal symbioses.

The fall armyworm converts maize endophytes into its own probiotics to detoxify benzoxazinoids and promote caterpillar growth

BACKGROUND

The fall armyworm (FAW, Spodoptera frugiperda) threatens maize production worldwide, and benzoxazinoids (Bxs) are known as the main secondary metabolites produced by maize to defend against FAW. However, we do not yet know whether and in what ways certain endophytes in the digestive system of FAW can metabolize Bxs, thus enhancing the fitness of FAW when feeding on maize.

RESULTS

Using Bxs as the sole carbon and nitrogen source, we isolated Pantoea dispersa from the guts of FAW. P. dispersa can colonize maize roots and leaves as indicated by GFP-labeling and further successfully established itself as an endophyte in the Malpighian tubules and the gut of FAW after FAW feeding activities. Once established, it can be vertically transmitted through FAW eggs, suggesting the potential that FAW can convert maize-derived endophytes into symbiotic bacteria for intergenerational transmission. The prevalence of P. dispersa in FAW guts and maize leaves was also confirmed over large geographic regions, indicating its evolutionary adaptation in fields. Bxs determination in the gut and frass of FAW combined with bioassays performance on maize bx2 mutants revealed that the colonization of P. dispersa can promote FAW growth by metabolizing Bxs rather than other metabolites. Additionally, genome and transcriptome analyses identified plasmid-borne genes, rather than chromosomes of this species, were crucial for Bxs metabolism. This was further validated through in vitro prokaryotic expression assays by expressing two candidate genes form the plasmid.

CONCLUSIONS

FAW can convert maize endophytes into its own probiotics to detoxify Bxs and thus enhance caterpillar growth. This represents a novel strategy for lepidopteran pests-transforming allies of the host into its own-thereby shedding light on the rapid spread of FAW and enhancing our understanding of ecological and evolutionary mechanisms underlying the pest-microbe-plant interactions. Video Abstract.

Pan-Chloroplast Genomes Reveal the Accession-Specific Marker for Gastrodia elata f. glauca

Gastrodia elata rhizomes have been applied as traditional medicinal materials for thousands of years. In China, G. elata f. elata (red flower and stem, Ge), G. elata f. viridis (green, Gv), and G. elata f. glauca (black, Gg) represent the primary cultivars in artificial cultivation. Although the annual output of G. elata amounts to 150,000 tons, only 20% is Gg. The long production period, low yield, and high quality of Gg led to its extremely high market prices. However, an effective method to identify this crude drug based solely on its morphological or chemical characteristics is lacking. In this study, the complete chloroplast genomes of three G. elata variants were sequenced using the Illumina HiSeq 2500 platform. Another 21 chloroplast genomes from Gastrodia species, which have published in previous reports, were combined and analyzed together. Our results showed that larger genomic sizes, fewer long tandem repeats, and more simple sequence repeats were the major features of the Gg chloroplast genomes. Phylogenetic analysis showed that the Gg samples were separately clustered in a subclade. Moreover, an accession-specific marker was successfully developed and validated for distinguishing additional rhizome samples. Our study provides investigations of the taxonomic relationships of Gastrodia species. The molecular marker will be useful for differentiating Gastrodia products on the market.

Unraveling the key mechanisms of Gastrodia elata continuous cropping obstacles: soil bacteria Massilia, Burkholderia-Caballeronia-Paraburkholderia, and Dyella along with soil metabolites 4-hydroxy-benzenemethanol and N-(2-butyl)-N-octadecyl-, ethyl ester as crucial indicators

Background

Tian-ma (Gastrodia elata) is a traditional medicinal herb found in China. It is used in healthy food and to treat various diseases, therefore cultivated extensively in southwest China. However, continuous cropping of this species has led to various obstacles, such as microbial disease and pest infestation, significantly affecting the production and development of valuable medicinal and food resources. As per the growth habit, soil is presumed to be the primary factor contributing to these obstacles, despite the known issues of continuous cropping obstacles in Gastrodia elata, such as microbial disease, there is a lack of comprehensive understanding of the specific soil bacterial communities and metabolites involved in these processes.

Methods

We analyzed soil samples collected during the year of Tian-ma cultivation (0 Year), after the Tian-ma harvest (1 Year), after two years (2 Year), and three years (3 Year) of fallowing post-cultivation using soil 16S rRNA metabarcoding sequencing by illumina platform and metabolomics (GC-MS/MS). Soil sample collected from the uncultivated field was used as the control (CK).

Results

Metabarcoding sequencing showed high bacterial alpha diversity during the cultivation of Tian-ma (0 Year) and the period of deterioration of soil bacterial community. (1 Year), with decreased anaerobic bacterial abundance and increased copiotrophic bacterial abundance. Bacteria associated with sulfur metabolism also showed increased abundance during the year of cropping obstacles. Further metabolomics approach identified 4-hydroxy-benzenemethanol as an indicator of Tian-ma continuous cropping obstacles. Besides, metabolites of the carbohydrate class were found to be the most abundant during the occurrence of continuous cropping obstacles of Gastrodia elata, suggesting that regulation of soil microbial diversity may be a critical factor in addressing these obstacles. Finally, the correlation analysis indicated a positive association between the abundance of some metabolite, e.g., carbamic acid, N-(2-butyl)-N-octadecyl-, ethyl ester detected after Tian-ma cultivation and the abundance of bacteria capable of degrading toxic metabolites, such as Massilia, Burkholderia-Caballeronia-Paraburkholderia, and Dyella.

Conclusion

This study has revealed the specific soil bacteria and metabolic factors related to the continuous cropping obstacles of Gastrodia elata. These findings not only deepen our understanding of the continuous cropping issues but also pave the way for developing effective strategies to overcome them.

Mitochondrial genome insights into the spatio-temporal distribution and genetic diversity of Dendrobium hancockii Rolfe (Orchidaceae)

Introduction

With its distinctive evolutionary rate and inheritance patterns separate from the nuclear genome, mitochondrial genome analysis has become a prominent focus of current research. Dendrobium hancockii Rolfe, a species of orchid with both medicinal and horticultural value, will benefit from the application of the fully assembled and annotated mitochondrial genome. This will aid in elucidating its phylogenetic relationships, comparative genomics, and population genetic diversity.

Methods

Based on sequencing results from Illumina combined with PacBio and Nanopore, the mitochondrial genome map of D. hancockii was constructed. Comparative analysis was conducted from the perspectives of phylogeny across multiple species, selection pressure on protein-coding genes, and homologous segments. The population diversity of D. hancockii was analyzed using single nucleotide polymorphism (SNP) data from the mitochondrial genome and single-copy nuclear genes.

Results and discussion

This research constructed a circular mitochondrial map for D. hancockii, spanning 523,952 bp, containing 40 unique protein-coding genes, 37 transfer RNA genes, and 4 ribosomal RNA genes. Comparative analysis of mitochondrial genes from 26 land plants revealed a conserved gene cluster, "rpl16-ccmFn-rps3-rps19," particularly within the Dendrobium genus. The mitochondrial genome of D. hancockii exhibits a lower point mutation rate but significant structural variation. Analysis of 103 resequencing samples identified 19,101 SNP sites, dividing D. hancockii into two major groups with limited gene flow between them, as supported by population diversity, genetic structure analysis, principal component analysis, and phylogenetic trees. The geographical distribution and genetic differentiation of D. hancockii into two major groups suggest a clear phytogeographical division, likely driven by ancient geological or climatic events. The close alignment of mitochondrial data with nuclear gene data highlights the potential of the mitochondrial genome for future studies on genetic evolution in this species.

Cytonuclear evolution in fully heterotrophic plants: lifestyles and gene function determine scenarios

BACKGROUND

Evidence shows that full mycoheterotrophs and holoparasites often have reduced plastid genomes with rampant gene loss, elevated substitution rates, and deeply altered to conventional evolution in mitochondrial genomes, but mechanisms of cytonuclear evolution is unknown. Endoparasitic Sapria himalayana and mycoheterotrophic Gastrodia and Platanthera guangdongensis represent different heterotrophic types, providing a basis to illustrate cytonuclear evolution. Here, we focused on nuclear-encoded plastid / mitochondrial (N-pt / mt) -targeting protein complexes, including caseinolytic protease (ClpP), ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), oxidative phosphorylation system (OXPHOS), DNA recombination, replication, and repair (DNA-RRR) system, and pentatricopeptide repeat (PPR) proteins, to identify evolutionary drivers for cytonuclear interaction.

RESULTS

The severity of gene loss of N-pt PPR and pt-RRR genes was positively associated with increased degree of heterotrophy in full mycoheterotrophs and S. himalayana, while N-mt PPR and mt-RRR genes were retained. Substitution rates of organellar and nuclear genes encoding N-pt/mt subunits in protein complexes were evaluated, cytonuclear coevolution was identified in S. himalayana, whereas disproportionate rates of evolution were observed in the OXPHOS complex in full mycoheterotrophs, only slight accelerations in substitution rates were identified in N-mt genes of full mycoheterotrophs.

CONCLUSIONS

Nuclear compensatory evolution was identified in protein complexes encoded by plastid and N-pt genes. Selection shaping codon preferences, functional constraint, mt-RRR gene regulation, and post-transcriptional regulation of PPR genes all facilitate mito-nuclear evolution. Our study enriches our understanding of genomic coevolution scenarios in fully heterotrophic plants.

What we know so far and what we can expect next: A molecular investigation of plant parasitism

The review explores parasitic plants' evolutionary success and adaptability, highlighting their widespread occurrence and emphasizing the role of an invasive organ called haustorium in nutrient acquisition from hosts. It discusses the genetic and physiological adaptations that facilitate parasitism, including horizontal gene transfer, and the impact of environmental factors like climate change on these relationships. It addresses the need for further research into parasitic plants' genomes and interactions with their hosts to better predict environmental changes' impacts.

Progress in systematics and biogeography of Orchidaceae

Orchidaceae are one of the largest families of angiosperms in terms of species richness. In the last decade, numerous studies have delved into reconstructing the phylogenetic framework of Orchidaceae, leveraging data from plastid, mitochondrial and nuclear sources. These studies have provided new insights into the systematics, diversification and biogeography of Orchidaceae, establishing a robust foundation for future research. Nevertheless, pronounced controversies persist regarding the precise placement of certain lineages within these phylogenetic frameworks. To address these discrepancies and deepen our understanding of the phylogenetic structure of Orchidaceae, we provide a comprehensive overview and analysis of phylogenetic studies focusing on contentious groups within Orchidaceae since 2015, delving into discussions on the underlying reasons for observed topological conflicts. We also provide a novel phylogenetic framework at the subtribal level. Furthermore, we examine the tempo and mode underlying orchid species diversity from the perspective of historical biogeography, highlighting factors contributing to extensive speciation. Ultimately, we delineate avenues for future research aimed at enhancing our understanding of Orchidaceae phylogeny and diversity.

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.

The Armillaria response to Gastrodia elata is partially mediated by strigolactone-induced changes in reactive oxygen species

Armillaria root diseases, caused by Armillaria spp., pose a significant threat to woody plants worldwide and result in substantial economic losses. However, certain species in the genus Armillaria can establish a unique symbiotic relationship with Gastrodia elata, which is the only known example of a plant benefiting from Armillaria. Although various plant signals that play a role in this interaction have been identified, the mechanism remains largely unknown from the Armillaria's perspective. In this study, we performed whole-genome sequencing of an Armillaria gallica strain named NRC001 isolated from G. elata. Comparative genomic analysis showed it is low-pathogenic Armillaria spp., which possesses 169 expanded gene families compared to high-pathogenic Armillaria spp. Among these expanded families, transcriptomic analysis revealed a significant increase in expression levels of four reactive oxygen species (ROS)-related gene families in A. gallica on G. elata compared to A. gallica on wood. Thus, a systematic survey of ROS-related gene families was carried out, and a total of 218 genes belonging to 44 ROS-related gene families in A. gallica were identified. Physiological experiments and transcriptome analysis showed that strigolactones (SLs) released by G. elata have a mediation impact on ROS, particularly enhancing the ROS scavenging activities by increasing the expression level and activity of several enzymes, such as catalase and glutathione reductase. Among the ROS-related genes, the aquaporin (AQP) is crucial as it is responsible for transporting hydrogen peroxide (H2O2) across the cell membrane. Five orthologs of AQP genes in A. gallica were identified and overexpressed in yeast. Only AgAQPA from the so-called 'other aquaglyceroporin' subfamily was demonstrated to be capable of mediating H2O2 transport in A. gallica. To our best knowledge, this is the first 'other aquaglyceroporins' gene in fungi to be identified as having transporter capacity. This study not only provides new insights into the mechanisms by which SL signaling regulates interactions between Armillaria and G. elata, but also sheds light on the function of fungal AQPs.

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.

LncRNA evolution and DNA methylation variation participate in photosynthesis pathways of distinct lineages of Populus

During the independent process of evolution in plants, photosynthesis appears to have been under convergent evolution to adapt to specific selection pressure in their geographical regions. However, it is unclear how lncRNA regulation and DNA methylation are involved in the phenotypic convergence in distinct lineages. Here, we present a large-scale comparative study of lncRNA transcription profile and whole-genome bisulfite sequencing (WGBS) data in two unrelated Populus species, selected from three relatively overlapping geographical regions. The results indicated that 39.75% lncRNAs of Populus tomentosa were shown to have homologous sequences in the 46.99% lncRNA of Populus simonii. Evolutionary analysis revealed that lncRNAs showed a rapid gain rate in the Populus lineage. Furthermore, co-expression networks in two Populus species identified eight lncRNAs that have the potential to simultaneously cis- or trans-regulate eight photosynthetic-related genes. These photosynthetic lncRNAs and genes were predominantly expressed in accessions from the southern region, indicating a conserved spatial expression in photosynthetic pathways in Populus. We also detected that most lncRNA targeted photosynthetic genes hypomethylated in promoter regions of Southern accessions compared with Northern accessions. Geographical DMRs correlated with genetic SNP variations in photosynthetic genes among Populus from the three geographic regions, indicating that DNA methylation coordinated with lncRNAs in convergent evolution of photosynthesis in Populus. Our results shed light on the evolutionary forces acting on patterns of lncRNA and DNA methylation, and provide a better understanding of the genetic and epigenetic mechanism in photosynthetic convergence evolution.

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.

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.

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.