The genome sequence of the orchid Phalaenopsis equestris

Diverse geotropic responses in the orchid family

In epiphytes, aerial roots are important to combat water-deficient, nutrient-poor, and high-irradiance microhabitats. However, whether aerial roots can respond to gravity and whether auxin plays a role in regulating aerial root development remain open-ended questions. Here, we investigated the gravitropic response of the epiphytic orchid Phalaenopsis aphrodite. Our data showed that aerial roots of P. aphrodite failed to respond to gravity, and this was correlated with a lack of starch granules/statolith sedimentation in the roots and the absence of the auxin efflux carrier PIN2 gene. Using an established auxin reporter, we discovered that auxin maximum was absent in the quiescent center of aerial roots of P. aphrodite. Also, gravity failed to trigger auxin redistribution in the root caps. Hence, loss of gravity sensing and gravity-dependent auxin redistribution may be the genetic factors contributing to aerial root development. Moreover, the architectural and functional innovations that achieve fast gravitropism in the flowering plants appear to be lost in both terrestrial and epiphytic orchids, but are present in the early diverged orchid subfamilies. Taken together, our findings provide physiological and molecular evidence to support the notion that epiphytic orchids lack gravitropism and suggest diverse geotropic responses in the orchid family.

The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application

Mycorrhizal symbioses between plants and fungi are vital for the soil structure, nutrient cycling, plant diversity, and ecosystem sustainability. More than 250 000 plant species are associated with mycorrhizal fungi. Recent advances in genomics and related approaches have revolutionized our understanding of the biology and ecology of mycorrhizal associations. The genomes of 250+ mycorrhizal fungi have been released and hundreds of genes that play pivotal roles in regulating symbiosis development and metabolism have been characterized. rDNA metabarcoding and metatranscriptomics provide novel insights into the ecological cues driving mycorrhizal communities and functions expressed by these associations, linking genes to ecological traits such as nutrient acquisition and soil organic matter decomposition. Here, we review genomic studies that have revealed genes involved in nutrient uptake and symbiosis development, and discuss adaptations that are fundamental to the evolution of mycorrhizal lifestyles. We also evaluated the ecosystem services provided by mycorrhizal networks and discuss how mycorrhizal symbioses hold promise for sustainable agriculture and forestry by enhancing nutrient acquisition and stress tolerance. Overall, unraveling the intricate dynamics of mycorrhizal symbioses is paramount for promoting ecological sustainability and addressing current pressing environmental concerns. This review ends with major frontiers for further research.

The road less taken: Dihydroflavonol 4-reductase inactivation and delphinidin anthocyanin loss underpins a natural intraspecific flower colour variation

Visual cues are of critical importance for the attraction of animal pollinators, however, little is known about the molecular mechanisms underpinning intraspecific floral colour variation. Here, we combined comparative spectral analysis, targeted metabolite profiling, multi-tissue transcriptomics, differential gene expression, sequence analysis and functional analysis to investigate a bee-pollinated orchid species, Glossodia major with common purple- and infrequent white-flowered morphs. We found uncommon and previously unreported delphinidin-based anthocyanins responsible for the conspicuous and pollinator-perceivable colour of the purple morph and three genetic changes underpinning the loss of colour in the white morph - (1) a loss-of-function (LOF; frameshift) mutation affecting dihydroflavonol 4-reductase (DFR1) coding sequence due to a unique 4-bp insertion, (2) specific downregulation of functional DFR1 expression and (3) the unexpected discovery of chimeric Gypsy transposable element (TE)-gene (DFR) transcripts with potential consequences to the genomic stability and post-transcriptional or epigenetic regulation of DFR. This is one of few known cases where regulatory changes and LOF mutation in an anthocyanin structural gene, rather than transcription factors, are important. Furthermore, if TEs prove to be a frequent source of mutation, the interplay between environmental stress-induced TE evolution and pollinator-mediated selection for adaptive colour variation may be an overlooked mechanism maintaining floral colour polymorphism in nature.

Crocus genome reveals the evolutionary origin of crocin biosynthesis

Crocus sativus (saffron) is a globally autumn-flowering plant, and its stigmas are the most expensive spice and valuable herb medicine. Crocus specialized metabolites, crocins, are biosynthesized in distant species, Gardenia (eudicot) and Crocus (monocot), and the evolution of crocin biosynthesis remains poorly understood. With the chromosome-level Crocus genome assembly, we revealed that two rounds of lineage-specific whole genome triplication occurred, contributing important roles in the production of carotenoids and apocarotenoids. According to the kingdom-wide identification, phylogenetic analysis, and functional assays of carotenoid cleavage dioxygenases (CCDs), we deduced that the duplication, site positive selection, and neofunctionalization of Crocus-specific CCD2 from CCD1 members are responsible for the crocin biosynthesis. In addition, site mutation of CsCCD2 revealed the key amino acids, including I143, L146, R161, E181, T259, and S292 related to the catalytic activity of zeaxanthin cleavage. Our study provides important insights into the origin and evolution of plant specialized metabolites, which are derived by duplication events of biosynthetic genes.

Breeding of ornamental orchids with focus on Phalaenopsis: current approaches, tools, and challenges for this century

Ornamental orchid breeding programs have been conducted to develop commercially valuable cultivars with improved characteristics of commercial interest, such as size, flower color, pattern, shape, and resistance to pathogens. Conventional breeding, including sexual hybridization followed by selection of desirable characteristics in plants, has so far been the main method for ornamental breeding, but other techniques, including mutation induction by polyploidization and gamma irradiation, and biotechnological techniques, such as genetic transformation, have also been studied and used in ornamental breeding programs. Orchids are one of the most commercially important families in floriculture industry, having very particular reproductive biology characteristics and being a well-studied group of ornamentals in terms of genetic improvement. The present review focuses on the conventional and biotechnological techniques and approaches specially employed in breeding Phalaenopsis orchids, the genus with highest worldwide importance as an ornamental orchid, highlighting the main limitations and strengths of the approaches. Furthermore, new opportunities and future prospects for ornamental breeding in the CRISPR/Cas9 genome editing era are also discussed. We conclude that conventional hybridization remains the most used method to obtain new cultivars in orchids. However, the emergence of the first biotechnology-derived cultivars, as well as the new biotechnological tools available, such as CRISPR-Cas9, rekindled the full potential of biotechnology approaches and their importance for improve ornamental orchid breeding programs.

Bringing CAM photosynthesis to the table: Paving the way for resilient and productive agricultural systems in a changing climate

Modern agricultural systems are directly threatened by global climate change and the resulting freshwater crisis. A considerable challenge in the coming years will be to develop crops that can cope with the consequences of declining freshwater resources and changing temperatures. One approach to meeting this challenge may lie in our understanding of plant photosynthetic adaptations and water use efficiency. Plants from various taxa have evolved crassulacean acid metabolism (CAM), a water-conserving adaptation of photosynthetic carbon dioxide fixation that enables plants to thrive under semi-arid or seasonally drought-prone conditions. Although past research on CAM has led to a better understanding of the inner workings of plant resilience and adaptation to stress, successful introduction of this pathway into C3 or C4 plants has not been reported. The recent revolution in molecular, systems, and synthetic biology, as well as innovations in high-throughput data generation and mining, creates new opportunities to uncover the minimum genetic tool kit required to introduce CAM traits into drought-sensitive crops. Here, we propose four complementary research avenues to uncover this tool kit. First, genomes and computational methods should be used to improve understanding of the nature of variations that drive CAM evolution. Second, single-cell 'omics technologies offer the possibility for in-depth characterization of the mechanisms that trigger environmentally controlled CAM induction. Third, the rapid increase in new 'omics data enables a comprehensive, multimodal exploration of CAM. Finally, the expansion of functional genomics methods is paving the way for integration of CAM into farming systems.

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

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

The Calvin-Benson-Bassham cycle in C4 and Crassulacean acid metabolism species

The Calvin-Benson-Bassham (CBB) cycle is the ancestral CO2 assimilation pathway and is found in all photosynthetic organisms. Biochemical extensions to the CBB cycle have evolved that allow the resulting pathways to act as CO2 concentrating mechanisms, either spatially in the case of C4 photosynthesis or temporally in the case of Crassulacean acid metabolism (CAM). While the biochemical steps in the C4 and CAM pathways are known, questions remain on their integration and regulation with CBB cycle activity. The application of omic and transgenic technologies is providing a more complete understanding of the biochemistry of C4 and CAM species and will also provide insight into the CBB cycle in these plants. As the global population increases, new solutions are required to increase crop yields and meet demands for food and other bioproducts. Previous work in C3 species has shown that increasing carbon assimilation through genetic manipulation of the CBB cycle can increase biomass and yield. There may also be options to improve photosynthesis in species using C4 photosynthesis and CAM through manipulation of the CBB cycle in these plants. This is an underexplored strategy and requires more basic knowledge of CBB cycle operation in these species to enable approaches for increased productivity.

Drought Stress Induced Different Response Mechanisms in Three Dendrobium Species under Different Photosynthetic Pathways

Many Dendrobium species, which hold a high status and value in traditional Chinese medicine, grow on barks and rocks in the wild, often encountering harsh environments and facing droughts. However, the molecular mechanisms underlying the shift in the photosynthetic pathway induced by drought remain unclear. To address this issue, three Dendrobium species with different photosynthetic pathways were selected for sequencing and transcriptome data analysis after drought treatment. The findings included 134.43 GB of sequencing data, with numerous Differentially Expressed Genes (DEGs) exhibiting different response mechanisms under drought stress. Gene Ontology (GO)-KEGG-based enrichment analysis of DEGs revealed that metabolic pathways contributed to drought tolerance and alterations in photosynthetic pathways. Phosphoenolpyruvate Carboxylase (PEPC) was subjected to phylogenetic tree construction, sequence alignment, and domain analysis. Under drought stress, variations were observed in the PEPC gene structure and expression among different Dendrobium species; the upregulation of Dc_gene2609 expression may be caused by dof-miR-384, which resulted in the shift from C3 photosynthesis to CAM, thereby improving drought tolerance in Dendrobium. This study revealed the expression patterns and roles of PEPC genes in enhancing plant drought tolerance and will provide an important basis for in-depth research on Dendrobium's adaptation mechanisms in arid environments.

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

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

Genome-wide identification and characterization of TCP gene family in Dendrobium nobile and their role in perianth development

TCP is a widely distributed, essential plant transcription factor that regulates plant growth and development. An in-depth study of TCP genes in Dendrobium nobile, a crucial parent in genetic breeding and an excellent model material to explore perianth development in Dendrobium, has not been conducted. We identified 23 DnTCP genes unevenly distributed across 19 chromosomes and classified them as Class I PCF (12 members), Class II: CIN (10 members), and CYC/TB1 (1 member) based on the conserved domain and phylogenetic analysis. Most DnTCPs in the same subclade had similar gene and motif structures. Segmental duplication was the predominant duplication event for TCP genes, and no tandem duplication was observed. Seven genes in the CIN subclade had potential miR319 and -159 target sites. Cis-acting element analysis showed that most DnTCP genes contained many developmental stress-, light-, and phytohormone-responsive elements in their promoter regions. Distinct expression patterns were observed among the 23 DnTCP genes, suggesting that these genes have diverse regulatory roles at different stages of perianth development or in different organs. For instance, DnTCP4 and DnTCP18 play a role in early perianth development, and DnTCP5 and DnTCP10 are significantly expressed during late perianth development. DnTCP17, 20, 21, and 22 are the most likely to be involved in perianth and leaf development. DnTCP11 was significantly expressed in the gynandrium. Specially, MADS-specific binding sites were present in most DnTCP genes putative promoters, and two Class I DnTCPs were in the nucleus and interacted with each other or with the MADS-box. The interactions between TCP and the MADS-box have been described for the first time in orchids, which broadens our understanding of the regulatory network of TCP involved in perianth development in orchids.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Orchid fruit and root movement analyzed using 2D photographs and a bioinformatics pipeline for processing sequential 3D scans

Premise

Most studies of the movement of orchid fruits and roots during plant development have focused on morphological observations; however, further genetic analysis is required to understand the molecular mechanisms underlying this phenomenon. A precise tool is required to observe these movements and harvest tissue at the correct position and time for transcriptomics research.

Methods

We utilized three-dimensional (3D) micro-computed tomography (CT) scans to capture the movement of fast-growing Erycina pusilla roots, and built an integrated bioinformatics pipeline to process 3D images into 3D time-lapse videos. To record the movement of slowly developing E. pusilla and Phalaenopsis equestris fruits, two-dimensional (2D) photographs were used.

Results

The E. pusilla roots twisted and resupinated multiple times from early development. The first period occurred in the early developmental stage (77-84 days after germination [DAG]) and the subsequent period occurred later in development (140-154 DAG). While E. pusilla fruits twisted 45° from 56-63 days after pollination (DAP), the fruits of P. equestris only began to resupinate a week before dehiscence (133 DAP) and ended a week after dehiscence (161 DAP).

Discussion

Our methods revealed that each orchid root and fruit had an independent direction and degree of torsion from the initial to the final position. Our innovative approaches produced detailed spatial and temporal information on the resupination of roots and fruits during orchid development.

Research advances on the gene regulation of floral development and color in orchids

Orchidaceae is one of the largest monocotyledon families and contributes significantly to worldwide biodiversity, with value in the fields of landscaping, medicine, and ecology. The diverse phenotypes and vibrant colors of orchid floral organs make them excellent research objects for investigating flower development and pigmentation. In recent years, a number of orchid genomes have been published, laying the molecular foundation for revealing flower development and color presentation. In this article, we review transcription factors, the structural genes responsible for the floral pigment synthesis pathways, the molecular mechanisms of flower morphogenesis, and the potential relationship between flower type and flower color. This study provides a theoretical reference for the research on molecular mechanisms related to flower morphogenesis and color presentation, genetic improvement, and new variety creation in orchids.

Genome-wide identification and expression analysis of growth-regulating factors in Dendrobium officinale and Dendrobium chrysotoxum

Background

Dendrobium, one of the largest genera in Orchidaceae, is popular not only for its aesthetic appeal but for its significant medicinal value. Growth-regulating factors (GRFs) play an essential role in plant growth and development. However, there is still a lack of information about the evolution and biological function analysis of the GRF gene family among Dendrobiumspecies.

Methods

Growth-regulating factors from Dendrobium officinale Kimura et Migo and Dendrobium chrysotoxum Lindl. were identified by HMMER and BLAST. Detailed bioinformatics analysis was conducted to explore the evolution and function of GRF gene family in D. officinale and D. chrysotoxum using genomic data, transcriptome data and qRT-PCR technology.

Results

Here, we evaluated the evolution of the GRF gene family based on the genome sequences of D. officinale and D. chrysotoxum. Inferred from phylogenetic trees, the GRF genes were classified into two clades, and each clade contains three subclades. Sequence comparison analysis revealed relatively conserved gene structures and motifs among members of the same subfamily, indicating a conserved evolution of GRF genes within Dendrobiumspecies. However, considering the distribution of orthologous DoGRFs and DcGRFs, and the differences in the number of GRFs among species, we suggest that the GRF gene family has undergone different evolutionary processes. A total of 361 cis-elements were detected, with 33, 141, and 187 related to plant growth and development, stress, and hormones, respectively. The tissue-specific expression of GRFs showed that DoGRF8 may have a significant function in the stem elongation of D. officinale. Moreover, four genes were up-regulated under Methyl-jasmonic acid/methyl jasmonate (MeJA) treatment, showing that DoGRFs and DcGRFs play a crucial role in stress response. These findings provide valuable information for further investigations into the evolution and function of GRF genes in D. officinale and D. chrysotoxum.

Transcriptome analysis of flower colour reveals the correlation between SNP and differential expression genes in Phalaenopsis

BACKGROUND

Phalaenopsis is an important ornamental plant that has great economic value in the world flower market as one of the most popular flower resources.

OBJECTIVE

To investigate the flower colour formation of Phalaenopsis at the transcription level, the genes involved in flower color formation were identified from RNA-seq in this study.

METHODS

In this study, white and purple petals of Phalaenopsis were collected and analyzed to obtained (1) differential expression genes (DEGs) between white and purple flower color and (2) the association between single nucleotide polymorphisms (SNP) mutations and DEGs at the transcriptome level.

RESULTS

The results indicated that a total of 1,175 DEGs were identified, and 718 and 457 of them were up- and down-regulated genes, respectively. Gene Ontology and pathway enrichment showed that the biosynthesis of the secondary metabolites pathway was key to color formation, and the expression of 12 crucial genes (C4H, CCoAOMT, F3'H, UA3'5'GT, PAL, 4CL, CCR, CAD, CALDH, bglx, SGTase, and E1.11.17) that are involved in the regulation of flower color in Phalaenopsis.

CONCLUSION

This study reported the association between the SNP mutations and DEGs for color formation at RNA level, and provides a new insight to further investigate the gene expression and its relationship with genetic variants from RNA-seq data in other species.

The CAM lineages of planet Earth

BACKGROUND AND SCOPE

The growth of experimental studies of crassulacean acid metabolism (CAM) in diverse plant clades, coupled with recent advances in molecular systematics, presents an opportunity to re-assess the phylogenetic distribution and diversity of species capable of CAM. It has been more than two decades since the last comprehensive lists of CAM taxa were published, and an updated survey of the occurrence and distribution of CAM taxa is needed to facilitate and guide future CAM research. We aimed to survey the phylogenetic distribution of these taxa, their diverse morphology, physiology and ecology, and the likely number of evolutionary origins of CAM based on currently known lineages.

RESULTS AND CONCLUSIONS

We found direct evidence (in the form of experimental or field observations of gas exchange, day-night fluctuations in organic acids, carbon isotope ratios and enzymatic activity) for CAM in 370 genera of vascular plants, representing 38 families. Further assumptions about the frequency of CAM species in CAM clades and the distribution of CAM in the Cactaceae and Crassulaceae bring the currently estimated number of CAM-capable species to nearly 7 % of all vascular plants. The phylogenetic distribution of these taxa suggests a minimum of 66 independent origins of CAM in vascular plants, possibly with dozens more. To achieve further insight into CAM origins, there is a need for more extensive and systematic surveys of previously unstudied lineages, particularly in living material to identify low-level CAM activity, and for denser sampling to increase phylogenetic resolution in CAM-evolving clades. This should allow further progress in understanding the functional significance of this pathway by integration with studies on the evolution and genomics of CAM in its many forms.

The Agavoideae: an emergent model clade for CAM evolutionary biology

Crassulacean acid metabolism - or CAM photosynthesis - was described in the early to mid-20th century, and our understanding of this metabolic pathway was later expanded upon through detailed biochemical analyses of carbon balance. Soon after, scientists began to study the ecophysiological implications of CAM, and a large part of this early work was conducted in the genus Agave, in the subfamily Agavoideae of the family Asparagaceae. Today, the Agavoideae continues to be important for the study of CAM photosynthesis, from the ecophysiology of CAM species, to the evolution of the CAM phenotype and to the genomics underlying CAM traits. Here we review past and current work on CAM in the Agavoideae, in particular highlighting the work of Park Nobel in Agave, and focusing on the powerful comparative system the Agavoideae has become for studying the origins of CAM. We also highlight new genomics research and the potential for studying intraspecific variation within species of the Agavoideae, particularly species in the genus Yucca. The Agavoideae has served as an important model clade for CAM research for decades, and undoubtedly will continue to help push our understanding of CAM biology and evolution in the future.

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.

Characterizing the chloroplast genome of Mammillaria elongata DC. 1828 in the Cactaceae family and unveiling its phylogenetic affinities within the genus Mammillaria

With its nearly 200 species, the Mammillaria genus is the most species-rich within the Cactaceae family, yet surprisingly, few of its chloroplast genomes have been studied. We focused on the species Mammillaria elongata DC. 1828, a petite cactus native to Mexico and favored by horticulturists, yet whose phylogenetic relationships remain uncertain due to a lack of genomic data. We extracted the DNA from a sample obtained in China, sequenced it using the NovaSeq 6000 platform, and assembled the chloroplast genome using GetOrganelle software. Our assembly resulted in a chloroplast genome of 110,981 base pairs with an overall GC content of 36.28%, which included 100 genes (95 unique). Notably, several protein-coding genes were absent. Phylogenetic analysis using 59 shared genes across nine Mammillaria species and one Obregonia species revealed that M. elongata and M. gracilis are closely related, suggesting a recent common ancestor and possible shared evolutionary pressures or ecological niches. This study provides crucial genomic data for M. elongata and hints at intriguing phylogenetic relationships within the Mammillaria genus.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

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.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

Recurrent allopolyploidizations diversify ecophysiological traits in marsh orchids (Dactylorhiza majalis s.l.)

Whole-genome duplication has shaped the evolution of angiosperms and other organisms, and is important for many crops. Structural reorganization of chromosomes and repatterning of gene expression are frequently observed in allopolyploids, with physiological and ecological consequences. Recurrent origins from different parental populations are widespread among polyploids, resulting in an array of lineages that provide excellent models to uncover mechanisms of adaptation to divergent environments in early phases of polyploid evolution. We integrate here transcriptomic and ecophysiological comparative studies to show that sibling allopolyploid marsh orchid species (Dactylorhiza, Orchidaceae) occur in different habitats (low nutrient fens vs. meadows with mesic soils) and are characterized by a complex suite of intertwined, pronounced ecophysiological differences between them. We uncover distinct features in leaf elemental chemistry, light-harvesting, photoprotection, nutrient transport and stomata activity of the two sibling allopolyploids, which appear to match their specific ecologies, in particular soil chemistry differences at their native sites. We argue that the phenotypic divergence between the sibling allopolyploids has a clear genetic basis, generating ecological barriers that maintain distinct, independent lineages, despite pervasive interspecific gene flow. This suggests that recurrent origins of polyploids bring about a long-term potential to trigger and maintain functional and ecological diversity in marsh orchids and other groups.

Defining Mechanisms of C3 to CAM Photosynthesis Transition toward Enhancing Crop Stress Resilience

Global climate change and population growth are persistently posing threats to natural resources (e.g., freshwater) and agricultural production. Crassulacean acid metabolism (CAM) evolved from C3 photosynthesis as an adaptive form of photosynthesis in hot and arid regions. It features the nocturnal opening of stomata for CO2 assimilation, diurnal closure of stomata for water conservation, and high water-use efficiency. To cope with global climate challenges, the CAM mechanism has attracted renewed attention. Facultative CAM is a specialized form of CAM that normally employs C3 or C4 photosynthesis but can shift to CAM under stress conditions. It not only serves as a model for studying the molecular mechanisms underlying the CAM evolution, but also provides a plausible solution for creating stress-resilient crops with facultative CAM traits. This review mainly discusses the recent research effort in defining the C3 to CAM transition of facultative CAM plants, and highlights challenges and future directions in this important research area with great application potential.

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

Acclimation to a combination of water deficit and nutrient deprivation through simultaneous increases in abscisic acid and bioactive jasmonates in the succulent plant Sempervivum tectorum L

Activation of hormonal responses defines the drought acclimation ability of plants and may condition their survival. However, aside ABA, little is known about the possible contribution of other phytohormones, such as jasmonates and salicylates, in the response of CAM plants to water deficit. Here, we aimed to study the physiological mechanisms underlying the stress tolerance of house leek (Sempervivum tectorum L.), a CAM plant adapted to survive harsh environments, to a combination of water deficit and nutrient deprivation. We exposed plants to the combination of these two abiotic stresses by withholding nutrient solution for 10 weeks and monitored their physiological response every two weeks by measuring various stress makers together with the accumulation of stress-related phytohormones and photoprotective molecules, such as tocopherols (vitamin E). Results showed that ABA content increased by 4.2-fold after four weeks of water deficit to keep later constant up to 10 weeks of stress, variations that occurred concomitantly with reductions in the relative leaf water content, which decreased by up to 20% only. The bioactive jasmonate, jasmonoyl-isoleucine was the other stress-related phytohormone that simultaneously increased under stress together with ABA. While contents of salicylic acid and the jasmonoyl-isoleucine precursors, 12-oxo-phytodienoic acid and jasmonic acid decreased with water deficit, those of jasmonoyl-isoleucine increased 3.6-fold at four weeks of stress. The contents of ABA and jasmonoyl-isoleucine correlated positively between them and with the content of α-tocopherol per unit of chlorophyll, thus suggesting a photoprotective activation role. It is concluded that S. tectorum not only withstands a combination of water deficit and nutrient deprivation for 10 weeks without any symptom of damage but also activates effective defense strategies through the simultaneous accumulation of ABA and the bioactive jasmonate form, jasmonoyl-isoleucine.

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

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

Complete chloroplast genome structural characterization of two Phalaenopsis (Orchidaceae) species and comparative analysis with their alliance

BACKGROUND

The taxonomy and infrageneric delimitation of Phalaenopsis Blume has been significantly disputed due to some overlapping morphological features between species related, which needed further evidence for clarification. The structural characterization of complete chloroplast genomes of P. storbatiana and P. wilsonii were analyzed and compared with those of related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny.

RESULTS

It was shown that chloroplast genomes of Phalaenopsis storbatiana and P. wilsonii had a typical quadripartite structure with conserved genome arrangements and moderate divergence. The chloroplast genomes of P. storbatiana and P. wilsonii were 145,885 bp and 145,445 bp in length, respectively, and shared a similar GC content of 36.8%. Gene annotations of two species revealed 109 single-copy genes consistently. In addition, 20 genes duplicated in the inverted regions, 16 genes each possessed one or more introns, and five ndh (NA (D)H dehydrogenase) genes were observed in both. Comparative analysis of the total cp genomes of P. storbatiana and P. wilsonii with those of other six related Phalaenopsis species confirmed the stable sequence identity for coding and non-coding regions and higher sequence variation in SC regions than IR regions. Most of their protein-coding genes had a high degree of codon preference. Moreover, 45 genes were discovered with significantly positive selection. However, different amplifications in IR regions were observed in these eight species. Phylogenetic analysis based on CDS from 60 species representing main clades in Orchidaceae indicated that Phalaenopsis species including P. stobartiana and P. wilsonii formed a monophyletic clade with high bootstrap nested in tribe Vandeae of Epidendroideae, which was consistent with those from previous studies.

CONCLUSIONS

The results could provide insight into understanding the plastome evolution and phylogenetic relationships of Phalaenopsis.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

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

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

Construction of a de novo assembly pipeline using multiple transcriptome data sets from Cypripedium macranthos (Orchidaceae)

The family Orchidaceae comprises the most species of any monocotyledonous family and has interesting characteristics such as seed germination induced by mycorrhizal fungi and flower morphology that co-adapted with pollinators. In orchid species, genomes have been decoded for only a few horticultural species, and there is little genetic information available. Generally, for species lacking sequenced genomes, gene sequences are predicted by de novo assembly of transcriptome data. Here, we devised a de novo assembly pipeline for transcriptome data from the wild orchid Cypripedium (lady slipper orchid) in Japan by mixing multiple data sets and integrating assemblies to create a more complete and less redundant contig set. Among the assemblies generated by combining various assemblers, Trinity and IDBA-Tran yielded good assembly with higher mapping rates and percentages of BLAST hit contigs and complete BUSCO. Using this contig set as a reference, we analyzed differential gene expression between protocorms grown aseptically or with mycorrhizal fungi to detect gene expressions required for mycorrhizal interaction. A pipeline proposed in this study can construct a highly reliable contig set with little redundancy even when multiple transcriptome data are mixed, and can provide a reference that is adaptable to DEG analysis and other downstream analysis in RNA-seq.

The spatiotemporal profile of Dendrobium huoshanense and functional identification of bHLH genes under exogenous MeJA using comparative transcriptomics and genomics

Introduction

Alkaloids are one of the main medicinal components of Dendrobium species. Dendrobium alkaloids are mainly composed of terpene alkaloids. Jasmonic acid (JA) induce the biosynthesis of such alkaloids, mainly by enhancing the expression of JA-responsive genes to increase plant resistance and increase the content of alkaloids. Many JA-responsive genes are the target genes of bHLH transcription factors (TFs), especially the MYC2 transcription factor.

Methods

In this study, the differentially expressed genes involved in the JA signaling pathway were screened out from Dendrobium huoshanense using comparative transcriptomics approaches, revealing the critical roles of basic helix-loop-helix (bHLH) family, particularly the MYC2 subfamily.

Results and discussion

Microsynteny-based comparative genomics demonstrated that whole genome duplication (WGD) and segmental duplication events drove bHLH genes expansion and functional divergence. Tandem duplication accelerated the generation of bHLH paralogs. Multiple sequence alignments showed that all bHLH proteins included bHLH-zip and ACT-like conserved domains. The MYC2 subfamily had a typical bHLH-MYC_N domain. The phylogenetic tree revealed the classification and putative roles of bHLHs. The analysis of cis-acting elements revealed that promoter of the majority of bHLH genes contain multiple regulatory elements relevant to light response, hormone responses, and abiotic stresses, and the bHLH genes could be activated by binding these elements. The expression profiling and qRT-PCR results indicated that bHLH subgroups IIIe and IIId may have an antagonistic role in JA-mediated expression of stress-related genes. DhbHLH20 and DhbHLH21 were considered to be the positive regulators in the early response of JA signaling, while DhbHLH24 and DhbHLH25 might be the negative regulators. Our findings may provide a practical reference for the functional study of DhbHLH genes and the regulation of secondary metabolites.

Identification and transcriptome data analysis of ARF family genes in five Orchidaceae species

The Orchidaceae is a large family of perennial herbs especially noted for the exceptional diversity of specialized flowers. Elucidating the genetic regulation of flowering and seed development of orchids is an important research goal with potential utility in orchid breeding programs. Auxin Response Factor (ARF) genes encode auxin-responsive transcription factors, which are involved in the regulation of diverse morphogenetic processes, including flowering and seed development. However, limited information on the ARF gene family in the Orchidaceae is available. In this study, 112 ARF genes were identified in the genomes of 5 orchid species (Apostasia shenzhenica, Dendrobium catenatum, Phalaenopsis aphrodite, Phalaenopsis equestris and Vanilla planifolia,). These genes were grouped into 7 subfamilies based on their phylogenetic relationships. Compared with the ARF family in model plants, such as Arabidopsis thaliana and Oryza sativa, one group of ARF genes involved in pollen wall synthesis has been lost during evolution of the Orchidaceae. This loss corresponds with absence of the exine in the pollinia. Through mining of the published genomic and transcriptomic data for the 5 orchid species: the ARF genes of subfamily 4 may play an important role in flower formation and plant growth, whereas those of subfamily 3 are potentially involved in pollen wall development. the study results provide novel insights into the genetic regulation of unique morphogenetic phenomena of orchids, which lay a foundation for further analysis of the regulatory mechanisms and functions of sexual reproduction-related genes in orchids.

A plastome phylogeny of Rumex (Polygonaceae) illuminates the divergent evolutionary histories of docks and sorrels

The genus Rumex L. (Polygonaceae) provides a unique system for investigating the evolutionary development of sex determination and molecular rate evolution. Historically, Rumex has been divided, both taxonomically and colloquially into two groups: 'docks' and 'sorrels'. A well-resolved phylogeny can help evaluate a genetic basis for this division. Here we present a plastome phylogeny for 34 species of Rumex, inferred using maximum likelihood criteria. The historical 'docks' (Rumex subgenus Rumex) were resolved as monophyletic. The historical 'sorrels' (Rumex subgenera Acetosa and Acetosella) were resolved together, though not monophyletic due to the inclusion of R. bucephalophorus (Rumex subgenus Platypodium). Emex is supported as its own subgenus within Rumex, instead of resolved as sister taxa. We found remarkably low nucleotide diversity among the docks, consistent with recent diversification in that group, especially as compared to the sorrels. Fossil calibration of the phylogeny suggested that the common ancestor for Rumex (including Emex) has origins in the lower Miocene (22.13 MYA). The sorrels appear to have subsequently diversified at a relatively constant rate. The origin of the docks, however, was placed in the upper Miocene, but with most speciation occurring in the Plio-Pleistocene.

Origin and diversification of a Himalayan orchid genus Pleione

Pleione is an orchid endemically distributed in high mountain areas across the Hengduan Mountains (HDM), Himalayas, Southeast Asia and South of China. The unique flower shapes, rich colors and immense medicinal importance of Pleione are valuable ornamental and economic resources. However, the phylogenetic relationships and evolutionary history of the genus have not yet been comprehensively resolved. Here, the evolutionary history of Pleione was investigated using single-copy gene single nucleotide polymorphisms and chloroplast genome datasets. The data revealed that Pleione could be divided into five clades. Discordance in topology between the two phylogenetic trees and network and D-statistic analyses indicated the occurrence of reticulate evolution in the genus. The evolution could be attributed to introgression and incomplete lineage sorting. Ancestral area reconstruction suggested that Pleione was originated from the HDM. Uplifting of the HDM drove rapid diversification by creating conditions favoring rapid speciation. This coincided with two periods of consolidation of the Asian monsoon climate, which caused the first rapid diversification of Pleione from 8.87 to 7.83 Mya, and a second rapid diversification started at around 4.05 Mya to Pleistocene. The interaction between Pleione and climate changes, especially the monsoons, led to the current distribution pattern and shaped the dormancy characteristic of the different clades. In addition to revealing the evolutionary relationship of Pleione with orogeny and climate changes, the findings of this study provide insights into the speciation and diversification mechanisms of plants in the East Asian flora.

The climate changes promoted the chloroplast genomic evolution of Dendrobium orchids among multiple photosynthetic pathways

Dendrobium orchids have multiple photosynthetic pathways, which can be used as a model system for studying the evolution of crassulacean acid metabolism (CAM). In this study, based on the results of the net photosynthetic rates (P_n), we classified Dendrobium species into three photosynthetic pathways, then employed and compared their chloroplast genomes. The Dendrobium chloroplast genomes have typical quartile structures, ranging from 150,841-153,038 bp. The apparent differences in GC content, sequence variability, and IR junctions of SSC/IR_B junctions (J_SBs) were measured within chloroplast genomes among different photosynthetic pathways. The phylogenetic analysis has revealed multiple independent CAM origins among the selected Dendrobium species. After counting insertions and deletions (InDels), we found that the occurrence rates and distribution densities among different photosynthetic pathways were inconsistent. Moreover, the evolution patterns of chloroplast genes in Dendrobium among three photosynthetic pathways were also diversified. Considering the diversified genome structure variations and the evolution patterns of protein-coding genes among Dendrobium species, we proposed that the evolution of the chloroplast genomes was disproportional among different photosynthetic pathways. Furthermore, climatic correlation revealed that temperature and precipitation have influenced the distribution among different photosynthetic pathways and promoted the foundation of CAM pathway in Dendrobium orchids. Based on our study, we provided not only new insights into the CAM evolution of Dendrobium but also provided beneficial genetic data resources for the further systematical study of Dendrobium.

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

PeGRF6-PeGIF1 complex regulates cell proliferation in the leaf of Phalaenopsis equestris

Phalaenopsis equestris is an ornamental plant with very large leaves. In this study, we identified genes related to the regulation of leaf development in Phalaenopsis and explored their mechanism of action. Sequence alignment and phylogenetic analyses revealed that PeGRF6 in the PeGRF family of P. equestris has similarities with the Arabidopsis genes AtGRF1 and AtGRF2, which are known to be involved in the regulation of leaf development. Among the PeGRFs, PeGRF6 was continuously and stably expressed at various stages of leaf development. The functions of PeGRF6 and of its complex formed with PeGIF1 in leaf development were verified by virus-induced gene silencing (VIGS) technology. The results show that the PeGRF6-PeGIF1 complex forms in the nucleus and positively regulates leaf cell proliferation via influencing cell size. Interestingly, VIGS suppression of PeGRF6 resulted in anthocyanin accumulation in Phalaenopsis leaves. Analyses of the regulatory mechanism of the miR396-PeGRF6 model based on the P. equestris small RNA library constructed here suggested that PeGRF6 transcripts are cleaved by Peq-miR396. These results show that, compared with PeGRF6 or PeGIF1 alone, the PeGRF6-PeGIF1 complex plays a more important role in the leaf development of Phalaenopsis, possibly by regulating the expression of cell cycle-related genes.

Genome-wide identification and analysis of the PAL genes from the orchids Apostasia shenzhenica, Dendrobium catenatum and Phalaenopsis equestris

Phenylalanine ammonia-lyase (PAL) is a key gateway enzyme that connects the phenylpropanoid pathway to primary metabolism. The phenylpropanoid pathway plays a vital role in the growth and environmental adaptation of many plants leading to the production of valuable bioactive compounds with industrial and medical applications. In the present study, nine putative PAL genes from three orchids were identified; five in Apostasia shenzhenica and two each in Dendrobium catenatum and Phalaenopsis equestris. Eighteen motifs and four major conserved functional domains were identified as reported in PAL proteins of other species. All the nine PALs were stable based on their computed physicochemical properties and localized in the cytoplasm. The three-dimensional structures of PALs revealed a homo-tetrameric structure consisting of four identical subunits. A total of 21 cis-regulatory elements with known functions were identified from the promoter regions of all PALs which are responsible for various plant responses to light, stress and growth regulators like auxins, gibberellins and abscisic acid. Phylogenetic analysis showed that the studied PAL proteins clustered in two major clades (clade I and II), placing dicot and monocot PALs in two separate monophyletic clades. In silico gene expression of the identified PALs in different vegetative and reproductive tissues revealed the differential expressions based on tissue type and disclosed that the expression of PAL genes was upregulated in all the tissues examined with an exception of PePAL leaf samples where no expression was detected, however, the same being highly expressed in reproductive tissues (PePAL1-labellum; PePAL2-sepal). In case of AsPALs, the expression was found to be highest in reproductive tissues (AsPAL4-maximum in inflorescence). On the other hand, the expression of DcPALs was found to be highest in vegetative tissues (DcPAL2-maximum in root). Based on the medicinal importance of orchids and the significant role of PAL genes in synthesis of bioactive compounds, the functional characterization of PAL genes can be further exploited in genetic improvement of medicinal orchids.Communicated by Ramaswamy H. Sarma.

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

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

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

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

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

Orchid NAC Transcription Factors: A Focused Analysis of CUPULIFORMIS Genes

Plant transcription factors are involved in different developmental pathways. NAC transcription factors (No Apical Meristem, Arabidopsis thaliana Activating Factor, Cup-shaped Cotyledon) act in various processes, e.g., plant organ formation, response to stress, and defense mechanisms. In Antirrhinum majus, the NAC transcription factor CUPULIFORMIS (CUP) plays a role in determining organ boundaries and lip formation, and the CUP homologs of Arabidopsis and Petunia are involved in flower organ formation. Orchidaceae is one of the most species-rich families of angiosperms, known for its extraordinary diversification of flower morphology. We conducted a transcriptome and genome-wide analysis of orchid NACs, focusing on the No Apical Meristem (NAM) subfamily and CUP genes. To check whether the CUP homologs could be involved in the perianth formation of orchids, we performed an expression analysis on the flower organs of the orchid Phalaenopsis aphrodite at different developmental stages. The expression patterns of the CUP genes of P. aphrodite suggest their possible role in flower development and symmetry establishment. In addition, as observed in other species, the orchid CUP1 and CUP2 genes seem to be regulated by the microRNA, miR164. Our results represent a preliminary study of NAC transcription factors in orchids to understand the role of these genes during orchid flower formation.

OrchidBase 5.0: updates of the orchid genome knowledgebase

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

Genome-wide identification and expression analysis of the GRAS gene family in Dendrobium chrysotoxum

The GRAS gene family encodes transcription factors that participate in plant growth and development phases. They are crucial in regulating light signal transduction, plant hormone (e.g. gibberellin) signaling, meristem growth, root radial development, response to abiotic stress, etc. However, little is known about the features and functions of GRAS genes in Orchidaceae, the largest and most diverse angiosperm lineage. In this study, genome-wide analysis of the GRAS gene family was conducted in Dendrobium chrysotoxum (Epidendroideae, Orchidaceae) to investigate its physicochemical properties, phylogenetic relationships, gene structure, and expression patterns under abiotic stress in orchids. Forty-six DchGRAS genes were identified from the D. chrysotoxum genome and divided into ten subfamilies according to their phylogenetic relationships. Sequence analysis showed that most DchGRAS proteins contained conserved VHIID and SAW domains. Gene structure analysis showed that intronless genes accounted for approximately 70% of the DchGRAS genes, the gene structures of the same subfamily were the same, and the conserved motifs were also similar. The Ka/Ks ratios of 12 pairs of DchGRAS genes were all less than 1, indicating that DchGRAS genes underwent negative selection. The results of cis-acting element analysis showed that the 46 DchGRAS genes contained a large number of hormone-regulated and light-responsive elements as well as environmental stress-related elements. In addition, the real-time reverse transcription quantitative PCR (RT-qPCR) experimental results showed significant differences in the expression levels of 12 genes under high temperature, drought and salt treatment, among which two members of the LISCL subfamily (DchGRAS13 and DchGRAS15) were most sensitive to stress. Taken together, this paper provides insights into the regulatory roles of the GRAS gene family in orchids.

A draft genome of the medicinal plant Cremastra appendiculata (D. Don) provides insights into the colchicine biosynthetic pathway

Cremastra appendiculata (D. Don) Makino is a rare terrestrial orchid with a high market value as an ornamental and Chinese traditional medicinal herb with a wide range of pharmacological properties. The pseudobulbs of C. appendiculata are one of the primary sources of the famous traditional Chinese medicine "Shancigu", which has been clinically used for treating many diseases, especially, as the main component to treat gout. The lack of genetic research and genome data restricts the modern development and clinical use of C. appendiculata. Here, we report a 2.3 Gb chromosome-level genome of C. appendiculata. We identify a series of candidates of 35 candidate genes responsible for colchicine biosynthesis, among which O-methyltransferase (OMT) gene exhibits an important role in colchicine biosynthesis. Co-expression analysis reveal purple and green-yellow module have close relationships with pseudobulb parts and comprise most of the colchicine pathway genes. Overall, our genome data and the candidate genes reported here set the foundation to decipher the colchicine biosynthesis pathways in medicinal plants.

WGDTree: a phylogenetic software tool to examine conditional probabilities of retention following whole genome duplication events

BACKGROUND

Multiple processes impact the probability of retention of individual genes following whole genome duplication (WGD) events. In analyzing two consecutive whole genome duplication events that occurred in the lineage leading to Atlantic salmon, a new phylogenetic statistical analysis was developed to examine the contingency of retention in one event based upon retention in a previous event. This analysis is intended to evaluate mechanisms of duplicate gene retention and to provide software to generate the test statistic for any genome with pairs of WGDs in its history.

RESULTS

Here a software package written in Python, 'WGDTree' for the analysis of duplicate gene retention following whole genome duplication events is presented. Using gene tree-species tree reconciliation to label gene duplicate nodes and differentiate between WGD and SSD duplicates, the tool calculates a statistic based upon the conditional probability of a gene duplicate being retained after a second whole genome duplication dependent upon the retention status after the first event. The package also contains methods for the simulation of gene trees with WGD events. After running simulations, the accuracy of the placement of events has been determined to be high. The conditional probability statistic has been calculated for Phalaenopsis equestris on a monocot species tree with a pair of consecutive WGD events on its lineage, showing the applicability of the method.

CONCLUSIONS

A new software tool has been created for the analysis of duplicate genes in examination of retention mechanisms. The software tool has been made available on the Python package index and the source code can be found on GitHub here: https://github.com/cnickh/wgdtree .