Chromosome-scale assembly of the Dendrobium chrysotoxum genome enhances the understanding of orchid evolution

Pangeneric genome analyses reveal the evolution and diversity of the orchid genus Dendrobium

Orchids constitute one of the most diverse families of angiosperms, yet their genome evolution and diversity remain unclear. Here we construct and analyse chromosome-scale de novo assembled genomes of 17 representative accessions spanning 12 sections in Dendrobium, one of the largest orchid genera. These accessions represent a broad spectrum of phenotypes, lineages and geographical distributions. We first construct haplotype-resolved genomes for a Dendrobium hybrid and uncover haplotypic variations and allelic imbalance in the heterozygous genome, demonstrating the significance of diverse ancestry. At Dendrobium genus-wide scale, we further elucidate phylogenetic relationships, evolutionary dynamics, entire gene repertoire, and the mechanisms of preserving ancient genetic variants and rapid recent genome evolution for habitat adaption. We also showcase distinctive evolutionary trajectories in MADS-box and PEBP families over 28 Ma. These results considerably contribute to unearthing the mystery of orchid origin, evolution and diversification, laying the foundation for efficient use of genetic diversity in breeding.

Expression Analysis and Functional Validation of DcTPSb1 in Terpene Synthesis of Dendrobium chrysotoxum

Terpenes are critical components of the floral fragrance component in Dendrobium chrysotoxum, synthesized by terpene synthase (TPS). Analysis of the D. chrysotoxum genome and transcriptional data revealed that the gene DcTPSb1 was significantly up-regulated during flowering periods, showing a strong correlation with the accumulation of aromatic monoterpenes in the floral components of Dendrobium chrysotoxum. Consequently, the DcTPSb1 gene was selected for further analysis. DcTPSb1 exhibited elevated expression levels in flowers among four organs (roots, stems, leaves, flowers) of D. chrysotoxum, with the highest expression observed during the blooming phase, which aligned with the accumulation of volatile terpenes during flowering. DcTPSb1, located in the chloroplasts, was identified as a member of the TPS-b subfamily associated with monoterpenes synthesis, showing close phylogenetic relationships with homologous proteins in related plant species. An analysis of the promoter region of DcTPSb1 indicated that it may be regulated by methyl jasmonate (MeJA) responsiveness. Functionally, DcTPSb1 was shown to catalyze the conversion of geranyl diphosphate (GPP) to linalool, ocimene, and (-)-α-pinitol in vitro. Overexpression of DcTPSb1 in tobacco resulted in a significant increase in terpenoid release during the blooming stage; however, the up-regulated substances did not include their catalytic products. The classification of DcTPSb1 as a terpene synthase capable of producing multiple products provides valuable insights into the complex biosynthesis of terpenes in orchids. These findings enhance our understanding of the functional diversity of DcTPSb1 and the processes involved in terpene biosynthesis in orchids.

Plant AT-rich protein and zinc-binding protein (PLATZ) family in Dendrobium huoshanense: identification, evolution and expression analysis

PLATZ (plant A/T-rich protein and zinc-binding protein) transcription factors are essential for plant growth, development, and responses to abiotic stress. The regulatory role of PLATZ genes in the environmental adaptation of D. huoshanense is inadequately comprehended. The genome-wide identification of D. huoshanense elucidates the functions and regulatory processes of the gene family. Our investigation encompassed the examination of PLATZ gene structures and chromosome distribution, the construction of the phylogenetic tree with its relatives, and the analysis of the cis-acting elements and expression profiles potentially implicated in growth and stress responses. Eleven DhPLATZs were classified into three clades (I, II, and III) according to their evolutionary homology. The distribution of these genes over six chromosomes indicated that both whole genome duplication (WGD) and segmental duplication events have contributed to the expansion of this gene family. The Ka/Ks analysis revealed a pattern of purifying selection after duplication occurrences, suggesting little alterations in functional divergence. The collinearity and microsynteny results revealed that the three DhPLATZ genes shared the same conserved domains as the paralogs from D. huoshanense and D. chrysotoxum. Expression profiling and quantitative analysis demonstrated that DhPLATZ genes had unique expression patterns in response to phytohormones and cold stress. Subcellular localization indicated that three DhPLATZ genes were expressed in the nucleus, suggesting their role as transcription factors. These findings enhance our understanding of PLATZ genes' involvement in D. huoshanense species and underscore their significance as important areas for further research.

The Gongora gibba genome assembly provides new insights into the evolution of floral scent in male euglossine bee-pollinated orchids

Orchidaceae is one of the most prominent flowering plant families, with many species exhibiting highly specialized reproductive and ecological adaptations. An estimated 10% of orchid species in the American tropics are pollinated by scent-collecting male euglossine bees; however, to date, there are no published genomes of species within this pollination syndrome. In this study, we present the first draft genome of an epiphytic orchid from the genus Gongora, a representative of the male euglossine bee-pollinated subtribe Stanhopeinae. The 1.83-Gb de novo genome with a scaffold N50 of 1.7 Mb was assembled using short- and long-read sequencing and chromosome capture (Hi-C) information. Over 17,000 genes were annotated, and 82.95% of the genome was identified as repetitive content. Furthermore, we identified and manually annotated 26 terpene synthase genes linked to floral scent biosynthesis and performed a phylogenetic analysis with other published orchid terpene synthase genes. The Gongora gibba genome assembly will serve as the foundation for future research to understand the genetic basis of floral scent biosynthesis and diversification in orchids.

Chromosome-scale genome assembly of Docynia delavayi provides new insights into the α-farnesene biosynthesis

Docynia delavayi is an economically significant fruit species with a high market potential due to the special aroma of its fruit. Here, a 653.34 Mb high-quality genome of D. delavayi was first reported, of which 93.8 % of the sequences (612.98 Mb) could be anchored to 17 chromosomes, containing 48,325 protein-coding genes. Ks analysis proved that two whole genome duplication (WGD) events occurred in D. delavayi, resulting in the expansion of genes associated with terpene biosynthesis, which promoted its fruit-specific aroma production. Combined multi-omics analysis, α-farnesene was detected as the most abundant aroma substance emitted by D. delavayi fruit during storage, meanwhile one α-farnesene synthase gene (AFS) and 15 transcription factors (TFs) were identified as the candidate genes potentially involved in α-farnesene biosynthesis. Further studies for the regulation network of α-farnesene biosynthesis revealed that DdebHLH, DdeERF1 and DdeMYB could activate the transcription of DdeAFS. To our knowledge, it is the first report that MYB TF plays a regulatory role in α-farnesene biosynthesis, which will greatly facilitate future breeding programs for D. delavayi.

Pangenome Identification and Analysis of Terpene Synthase Gene Family Members in Gossypium

Terpene synthases (TPSs), key gatekeepers in the biosynthesis of herbivore-induced terpenes, are pivotal in the diversity of terpene chemotypes across and within plant species. Here, we constructed a gene-based pangenome of the Gossypium genus by integrating the genomes of 17 diploid and 10 tetraploid species. Within this pangenome, 208 TPS syntelog groups (SGs) were identified, comprising 2 core SGs (TPS5 and TPS42) present in all 27 analyzed genomes, 6 softcore SGs (TPS11, TPS12, TPS13, TPS35, TPS37, and TPS47) found in 25 to 26 genomes, 131 dispensable SGs identified in 2 to 24 genomes, and 69 private SGs exclusive to a single genome. The mutational load analysis of these identified TPS genes across 216 cotton accessions revealed a great number of splicing variants and complex splicing patterns. The nonsynonymous/synonymous Ka/Ks value for all 52 analyzed TPS SGs was less than one, indicating that these genes were subject to purifying selection. Of 208 TPS SGs encompassing 1795 genes, 362 genes derived from 102 SGs were identified as atypical and truncated. The structural analysis of TPS genes revealed that gene truncation is a major mechanism contributing to the formation of atypical genes. An integrated analysis of three RNA-seq datasets from cotton plants subjected to herbivore infestation highlighted nine upregulated TPSs, which included six previously characterized TPSs in G. hirsutum (AD1_TPS10, AD1_TPS12, AD1_TPS40, AD1_TPS42, AD1_TPS89, and AD1_TPS104), two private TPSs (AD1_TPS100 and AD2_TPS125), and one atypical TPS (AD2_TPS41). Also, a TPS-associated coexpression module of eight genes involved in the terpenoid biosynthesis pathway was identified in the transcriptomic data of herbivore-infested G. hirsutum. These findings will help us understand the contributions of TPS family members to interspecific terpene chemotypes within Gossypium and offer valuable resources for breeding insect-resistant cotton cultivars.

Identification and Analysis of Aluminum-Activated Malate Transporter Gene Family Reveals Functional Diversification in Orchidaceae and the Expression Patterns of Dendrobium catenatum Aluminum-Activated Malate Transporters

Aluminum-activated malate transporter (ALMT) genes play an important role in aluminum ion (Al3+) tolerance, fruit acidity, and stomatal movement. Although decades of research have been carried out in many plants, there is little knowledge about the roles of ALMT in Orchidaceae. In this study, 34 ALMT genes were identified in the genomes of four orchid species. Specifically, ten ALMT genes were found in Dendrobium chrysotoxum and D. catenatum, and seven were found in Apostasia shenzhenica and Phalaenopsis equestris. These ALMT genes were further categorized into four clades (clades 1-4) based on phylogenetic relationships. Sequence alignment and conserved motif analysis revealed that most orchid ALMT proteins contain conserved regions (TM1, GABA binding motif, and WEP motif). We also discovered a unique motif (19) belonging to clade 1, which can serve as a specifically identified characteristic. Comparison with the gene structure of AtALMT genes (Arabidopsis thaliana) showed that the gene structure of ALMT was conserved across species, but the introns were longer in orchids. The promoters of orchid ALMT genes contain many light-responsive and hormone-responsive elements, suggesting that their expression may be regulated by light and phytohormones. Chromosomal localization and collinear analysis of D. chrysotoxum indicated that tandem duplication (TD) is the main reason for the difference in the number of ALMT genes in these orchids. D. catenatum was chosen for the RT-qPCR experiment, and the results showed that the DcaALMT gene expression pattern varied in different tissues. The expression of DcaALMT1-9 was significantly changed after ABA treatment. Combining the circadian CO2 uptake rate, titratable total acid, and RT-qPCR data analysis, most DcaALMT genes were highly expressed at night and around dawn. The result revealed that DcaALMT genes might be involved in photosynthate accumulation. The above study provides more comprehensive information for the ALMT gene family in Orchidaceae and a basis for subsequent functional analysis.

Multi-omics analyses and botanical perfumer hypothesis provide insights into the formation and maintenance of aromatic characteristics of Dendrobium loddigesii flowers

Dendrobium loddigesii, a member of the Orchidaceae family, is a valuable horticultural crop known for its aromatic qualities. However, the mechanisms responsible for the development of its aromatic characteristics remain poorly understood. To elucidate these underlying mechanisms, we assembled the first chromosome-level reference genome of D. loddigesii using PacBio HiFi-reads, Illumina short-reads, and Hi-C data. The assembly comprises 19 pseudochromosomes with N50 contig and N50 scaffold sizes of 55.15 and 89.94 Mb, respectively, estimating the genome size to be 1.68 Gb, larger than that of other sequenced Dendrobium species. During the flowering stages, we conducted a comprehensive analysis combining volatilomics and transcriptomics to understand the characteristics and biosynthetic mechanisms pathways of the floral scent. Our findings emphasize the significant contribution of aromatic terpenoids, especially monoterpenoids, in defining the floral aroma. Furthermore, we identified two crucial terpene synthase (TPS) genes that play a key role in maintaining the aroma during flowering. Through the integration volatilomics data with catalytic assays of DlTPSbs proteins, we identified specific compounds responsible for the aromatic characteristics of D. loddigesii. This integrated analysis of the genome, transcriptome, and volatilome, offers valuable insights into the development and preservation of D. loddigesii's aromatic characteristics, setting the stage for further exploration of the botanical perfumer hypothesis.

Genome-wide identification of Ankyrin (ANK) repeat gene families in three Dendrobium species and the expression of ANK genes in D. officinale under gibberellin and abscisic acid treatments

BACKGROUND

Dendrobium Sw. represents one of the most expansive genera within the Orchidaceae family, renowned for its species' high medicinal and ornamental value. In higher plants, the ankyrin (ANK) repeat protein family is characterized by a unique ANK repeat domain, integral to a plethora of biological functions and biochemical activities. The ANK gene family plays a pivotal role in various plant physiological processes, including stress responses, hormone signaling, and growth. Hence, investigating the ANK gene family and identifying disease-resistance genes in Dendrobium is of paramount importance.

RESULTS

This research identified 78 ANK genes in Dendrobium officinale Kimura et Migo, 77 in Dendrobium nobile Lindl., and 58 in Dendrobium chrysotoxum Lindl. Subsequently, we conducted comprehensive bioinformatics analyses on these ANK gene families, encompassing gene classification, chromosomal localization, phylogenetic relationships, gene structure and motif characterization, cis-acting regulatory element identification, collinearity assessment, protein-protein interaction network construction, and gene expression profiling. Concurrently, three DoANK genes (DoANK14, DoANK19, and DoANK47) in D. officinale were discerned to indirectly activate the NPR1 transcription factor in the ETI system via SA, thereby modulating the expression of the antibacterial PR gene. Hormonal treatments with GA3 and ABA revealed that 17 and 8 genes were significantly up-regulated, while 4 and 8 genes were significantly down-regulated, respectively. DoANK32 was found to localize to the ArfGAP gene in the endocytosis pathway, impacting vesicle transport and the polar movement of auxin.

CONCLUSION

Our findings provide a robust framework for the taxonomic classification, evolutionary analysis, and functional prediction of Dendrobium ANK genes. The three highlighted ANK genes (DoANK14, DoANK19, and DoANK47) from D. officinale may prove valuable in disease resistance and stress response research. DoANK32 is implicated in the morphogenesis and development of D. officinale through its role in vesicular transport and auxin polarity, with subcellular localization studies confirming its presence in the nucleus and cell membrane. ANK genes displaying significant expression changes in response to hormonal treatments could play a crucial role in the hormonal response of D. officinale, potentially inhibiting its growth and development through the modulation of plant hormones such as GA3 and ABA.

Recent Advances and New Insights in Genome Analysis and Transcriptomic Approaches to Reveal Enzymes Associated with the Biosynthesis of Dendrobine-Type Sesquiterpenoid Alkaloids (DTSAs) from the Last Decade

Dendrobium species, which are perennial herbs widely distributed in tropical and subtropical regions, are notable for their therapeutic properties attributed to various bioactive compounds, including dendrobine-type sesquiterpenoid alkaloids (DTSAs). The objective of this review article is to provide a comprehensive overview of recent advances in the biosynthesis of DTSAs, including their extraction from Dendrobium species and endophytes, elucidation of associated genes through genomic and transcriptomic sequencing in both Dendrobium spp. and endophytes, exploration of the biosynthetic pathways of DTSAs, and drawing conclusions and outlining future perspectives in this field. Alkaloids, predominantly nitrogen-containing compounds found in medicinal orchids, include over 140 types discovered across more than 50 species. DTSAs, identified in 37 picrotoxane alkaloids, have a distinctive five-membered nitrogen heterocyclic ring. This review highlights endophytic fungi as alternative sources of DTSAs, emphasizing their potential in pharmaceutical applications when plant-derived compounds are scarce or complex. Genomic and transcriptomic sequencing of Dendrobium spp. and their endophytes has identified key genes involved in DTSAs biosynthesis, elucidating pathways such as the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways. Genes encoding enzymes, such as acetyl-CoA C-acetyltransferase and diphosphomevalonate decarboxylase, are positively associated with dendrobine production. Despite significant advancements, the complexity of terpenoid biosynthesis in different subcellular compartments remains a challenge. Future research should focus on leveraging high-quality genomic data and omics technologies to further understand and manipulate the biosynthetic pathways of DTSAs and enhance their medicinal use.

Genome-Wide Identification and Expression Pattern Analysis of GATA Gene Family in Orchidaceae

The GATA transcription factors play crucial roles in plant growth, development, and responses to environmental stress. Despite extensive studies of GATA genes in many plants, their specific functions and mechanisms in orchids remain unexplored. In our study, a total of 149 GATA genes were identified in the genomes of seven sequenced orchid species (20 PeqGATAs, 23 CgGATAs, 24 CeGATAs, 23 DcaGATAs, 20 DchGATAs, 27 DnoGATAs, and 12 GelGATAs), classified into four subfamilies. Subfamily I typically contains genes with two exons, while subfamily II contains genes with two or three exons. Most members of subfamilies III and IV have seven or eight exons, with longer introns compared to subfamilies I and II. In total, 24 pairs (CgGATAs-DchGATAs), 27 pairs (DchGATAs-DnoGATAs), and 14 pairs (DnoGATAs-GelGATAs) of collinear relationships were identified. Cis-acting elements in GATA promoters were mainly enriched in abscisic acid (ABA) response elements and methyl jasmonate (MeJA) elements. Expression patterns and RT-qPCR analysis revealed that GATAs are involved in the regulation of floral development in orchids. Furthermore, under high-temperature treatment, GL17420 showed an initial increase followed by a decrease, GL18180 and GL17341 exhibited a downregulation followed by upregulation and then a decrease, while GL30286 and GL20810 displayed an initial increase followed by slight inhibition and then another increase, indicating diverse regulatory mechanisms of different GATA genes under heat stress. This study explores the function of GATA genes in orchids, providing a theoretical basis and potential genetic resources for orchid breeding and stress resistance improvement.

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.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

Metabolic and Transcriptomic Profile Revealing the Differential Accumulating Mechanism in Different Parts of Dendrobium nobile

Dendrobium nobile is an important orchid plant that has been used as a traditional herb for many years. For the further pharmaceutical development of this resource, a combined transcriptome and metabolome analysis was performed in different parts of D. nobile. First, saccharides, organic acids, amino acids and their derivatives, and alkaloids were the main substances identified in D. nobile. Amino acids and their derivatives and flavonoids accumulated strongly in flowers; saccharides and phenols accumulated strongly in flowers and fruits; alkaloids accumulated strongly in leaves and flowers; and a nucleotide and its derivatives and organic acids accumulated strongly in leaves, flowers, and fruits. Simultaneously, genes for lipid metabolism, terpenoid biosynthesis, and alkaloid biosynthesis were highly expressed in the flowers; genes for phenylpropanoids biosynthesis and flavonoid biosynthesis were highly expressed in the roots; and genes for other metabolisms were highly expressed in the leaves. Furthermore, different members of metabolic enzyme families like cytochrome P450 and 4-coumarate-coA ligase showed differential effects on tissue-specific metabolic accumulation. Members of transcription factor families like AP2-EREBP, bHLH, NAC, MADS, and MYB participated widely in differential accumulation. ATP-binding cassette transporters and some other transporters also showed positive effects on tissue-specific metabolic accumulation. These results systematically elucidated the molecular mechanism of differential accumulation in different parts of D. nobile and enriched the library of specialized metabolic products and promising candidate genes.

The Evolution of the WUSCHEL-Related Homeobox Gene Family in Dendrobium Species and Its Role in Sex Organ Development in D. chrysotoxum

The WUSCHEL-related homeobox (WOX) transcription factor plays a vital role in stem cell maintenance and organ morphogenesis, which are essential processes for plant growth and development. Dendrobium chrysotoxum, D. huoshanense, and D. nobile are valued for their ornamental and medicinal properties. However, the specific functions of the WOX gene family in Dendrobium species are not well understood. In our study, a total of 30 WOX genes were present in the genomes of the three Dendrobium species (nine DchWOXs, 11 DhuWOXs, and ten DnoWOXs). These 30 WOXs were clustered into ancient clades, intermediate clades, and WUS/modern clades. All 30 WOXs contained a conserved homeodomain, and the conserved motifs and gene structures were similar among WOXs belonging to the same branch. D. chrysotoxum and D. huoshanense had one pair of fragment duplication genes and one pair of tandem duplication genes, respectively; D. nobile had two pairs of fragment duplication genes. The cis-acting regulatory elements (CREs) in the WOX promoter region were mainly enriched in the light response, stress response, and plant growth and development regulation. The expression pattern and RT-qPCR analysis revealed that the WOXs were involved in regulating the floral organ development of D. chrysotoxum. Among them, the high expression of DchWOX3 suggests that it might be involved in controlling lip development, whereas DchWOX5 might be involved in controlling ovary development. In conclusion, this work lays the groundwork for an in-depth investigation into the functions of WOX genes and their regulatory role in Dendrobium species' floral organ development.

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

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

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

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

Genome-Wide Identification and Expression Profile Analysis of Sugars Will Eventually Be Exported Transporter (SWEET) Genes in Zantedeschia elliottiana and Their Responsiveness to Pectobacterium carotovora subspecies Carotovora (Pcc) Infection

SWEET, sugars will eventually be exported transporter, is a novel class of sugar transporter proteins that can transport sugars across membranes down a concentration gradient. It plays a key role in plant photosynthetic assimilates, phloem loading, nectar secretion from nectar glands, seed grouting, pollen development, pathogen interactions, and adversity regulation, and has received widespread attention in recent years. To date, systematic analysis of the SWEET family in Zantedeschia has not been documented, although the genome has been reported in Zantedeschia elliottiana. In this study, 19 ZeSWEET genes were genome-wide identified in Z. elliottiana, and unevenly located in 10 chromosomes. They were further clustered into four clades by a phylogenetic tree, and almost every clade has its own unique motifs. Synthetic analysis confirmed two pairs of segmental duplication events of ZeSWEET genes. Heatmaps of tissue-specific and Pectobacterium carotovora subsp. Carotovora (Pcc) infection showed that ZeSWEET genes had different expression patterns, so SWEETs may play widely varying roles in development and stress tolerance in Zantedeschia. Moreover, quantitative reverse transcription-PCR (qRT-PCR) analysis revealed that some of the ZeSWEETs responded to Pcc infection, among which eight genes were significantly upregulated and six genes were significantly downregulated, revealing their potential functions in response to Pcc infection. The promoter sequences of ZeSWEETs contained 51 different types of the 1380 cis-regulatory elements, and each ZeSWEET gene contained at least two phytohormone responsive elements and one stress response element. In addition, a subcellular localization study indicated that ZeSWEET07 and ZeSWEET18 were found to be localized to the plasma membrane. These findings provide insights into the characteristics of SWEET genes and contribute to future studies on the functional characteristics of ZeSWEET genes, and then improve Pcc infection tolerance in Zantedeschia through molecular breeding.

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.

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.

Approaches to increase the validity of gene family identification using manual homology search tools

Identifying homologs is an important process in the analysis of genetic patterns underlying traits and evolutionary relationships among species. Analysis of gene families is often used to form and support hypotheses on genetic patterns such as gene presence, absence, or functional divergence which underlie traits examined in functional studies. These analyses often require precise identification of all members in a targeted gene family. Manual pipelines where homology search and orthology assignment tools are used separately are the most common approach for identifying small gene families where accurate identification of all members is important. The ability to curate sequences between steps in manual pipelines allows for simple and precise identification of all possible gene family members. However, the validity of such manual pipeline analyses is often decreased by inappropriate approaches to homology searches including too relaxed or stringent statistical thresholds, inappropriate query sequences, homology classification based on sequence similarity alone, and low-quality proteome or genome sequences. In this article, we propose several approaches to mitigate these issues and allow for precise identification of gene family members and support for hypotheses linking genetic patterns to functional traits.

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.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

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.

Genome-Wide Identification of TCP Gene Family in Dendrobium and Their Expression Patterns in Dendrobium chrysotoxum

The TCP gene family are plant-specific transcription factors that play important roles in plant growth and development. Dendrobium chrysotoxum, D. nobile, and D. huoshanense are orchids with a high ornamental value, but few studies have investigated the specific functions of TCPs in Dendrobium flower development. In this study, we used these three Dendrobium species to analyze TCPs, examining their physicochemical properties, phylogenetic relationships, gene structures, and expression profiles. A total of 50 TCPs were identified across three Dendrobium species; they were divided into two clades-Class-I (PCF subfamily) and Class-II (CIN and CYC/TB1 subfamilies)-based on their phylogenetic relationships. Our sequence logo analysis showed that almost all Dendrobium TCPs contain a conserved TCP domain, as well as the existence of fewer exons, and the cis-regulatory elements of the TCPs were mostly related to light response. In addition, our transcriptomic data and qRT-PCR results showed that DchTCP2 and DchTCP13 had a significant impact on lateral organs. Moreover, changes in the expression level of DchTCP4 suggested its important role in the phenotypic variation of floral organs. Therefore, this study provides a significant reference for the further exploration of TCP gene functions in the regulation of different floral organs in Dendrobium orchids.

Genome-Wide Identification and Analysis of TCP Gene Family among Three Dendrobium Species

Dendrobium orchids, which are among the most well-known species of orchids, are appreciated for their aesthetic appeal across the globe. Furthermore, due to their strict living conditions, they have accumulated high levels of active ingredients, resulting not only in their medicinal value but also in their strong ability to respond to harsh environments. The TCP gene family plays an important role in plant growth and development, and signal transduction. However, these genes have not been systematically investigated in Dendrobium species. In this study, we detected a total of 24, 23, and 14 candidate TCP members in the genome sequences of D. officinale, D. nobile, and D. chrysotoxum, respectively. These genes were classified into three clades on the basis of a phylogenetic analysis. The TCP gene numbers among Dendrobium species were still highly variable due to the independent loss of genes in the CIN clade. However, only three gene duplication events were detected, with only one tandem duplication event (DcTCP9/DcTCP10) in D. chrysotoxum and two pairs of paralogous DoTCP gene duplication events (DoTCP1/DoTCP23 and DoTCP16/DoTCP24) in D. officinale. A total of 25 cis-acting elements of TCPs related to hormone/stress and light responses were detected. Among them, the proportions of hormone response, light response, and stress response elements in D. officinale (100/421, 127/421, and 171/421) were similar to those in D. nobile (83/352, 87/352, and 161/352). Using qRT-PCR to determine their expression patterns under MeJA treatment, four DoTCPs (DoTCP2, DoTCP4, DoTCP6, and DoTCP14) were significantly upregulated under MeJA treatment, which indicates that TCP genes may play important roles in responding to stress. Under ABA treatment, seven DoTCPs (DoTCP3, DoTCP7, DoTCP9, DoTCP11, DoTCP14, DoTCP15, and DoTCP21) were significantly upregulated, indicating that TCP genes may also play an important role in hormone response. Therefore, these results can provide useful information for studying the evolution and function of TCP genes in Dendrobium species.

Genome-wide identification, characterization and transcriptional profile of the SWEET gene family in Dendrobium officinale

BACKGROUND

Dendrobium officinale Kimura et Migo (D. officinale) is a well-known traditional Chinese medicine with high content polysaccharides in stems. The SWEET (Sugars Will Eventually be Exported Transporters) family is a novel class of sugar transporters mediating sugar translocation among adjacent cells of plants. The expression patterns of SWEETs and whether they are associated with stress response in D. officinale remains uncovered.

RESULTS

Here, 25 SWEET genes were screened out from D. officinale genome, most of which typically contained seven transmembrane domains (TMs) and harbored two conserved MtN3/saliva domains. Using multi-omics data and bioinformatic approaches, the evolutionary relationship, conserved motifs, chromosomal location, expression patterns, correlationship and interaction network were further analyzed. DoSWEETs were intensively located in nine chromosomes. Phylogenetic analysis revealed that DoSWEETs were divided into four clades, and conserved motif 3 specifically existed in DoSWEETs from clade II. Different tissue-specific expression patterns of DoSWEETs suggested the division of their roles in sugar transport. In particular, DoSWEET5b, 5c, and 7d displayed relatively high expression levels in stems. DoSWEET2b and 16 were significantly regulated under cold, drought, and MeJA treatment, which were further verified using RT-qPCR. Correlation analysis and interaction network prediction discovered the internal relationship of DoSWEET family.

CONCLUSIONS

Taken together, the identification and analysis of the 25 DoSWEETs in this study provide basic information for further functional verification in D. officinale.

Comprehensive chemical profiling of two Dendrobium species and identification of anti-hepatoma active constituents from Dendrobium chrysotoxum by network pharmacology

BACKGROUND

Dendrobium nobile and Dendrobium chrysotoxum are important species of the genus Dendrobium and have great economic and medicinal value. However, the medicinal properties of these two plants remain poorly understood. This study aimed to investigate the medical properties of D. nobile and D. chrysotoxum by conducting a comprehensive chemical profiling of the two plants. Additionally, active compounds and predictive targets for anti-hepatoma activity in D. chrysotoxum extracts were identified using Network Pharmacology.

RESULTS

Chemical profiling showed that altogether 65 phytochemicals were identified from D. nobile and D. chrysotoxum, with major classes as alkaloids, terpenoids, flavonoids, bibenzyls and phenanthrenes. About 18 compounds were identified as the important differential metabolites in D. nobile and D. chrysotoxum. Furtherly, CCK-8 results showed that the extracts of stems and leaves of D. nobile and D. chrysotoxum could inhibit the growth of Huh-7 cells, and the anti-hepatoma activity of extracts were dose-dependent. Among the extracts, the extract of D. chrysotoxum showed significant anti-hepatoma activity. In order to find the potential mechanism of anti-hepatoma activity of D. chrysotoxum, five key compounds and nine key targets were obtained through constructing and analyzing the compound-target-pathway network. The five key compounds were chrysotobibenzyl, chrysotoxin, moscatilin, gigantol and chrysotoxene. Nine key targets, including GAPDH, EGFR, ESR1, HRAS, SRC, CCND1, HIF1A, ERBB2 and MTOR, could be considered as the core targets of the anti-hepatoma activity of D. chrysotoxum.

CONCLUSIONS

In this study, the chemical composition difference and anti-hepatoma activity of stems and leaves of D. nobile and D. chrysotoxum were compared, and the potential anti-hepatoma mechanism of D. chrysotoxum was revealed in a multi-target and multi-pathway manner.

Aromatic Terpenes and Their Biosynthesis in Dendrobium, and Conjecture on the Botanical Perfumer Mechanism

This review presents a systematic analysis of the studies on volatiles in Dendrobium. Among the various components, aromatic terpenes are a crucial component in the development of the aromatic characteristics of Dendrobium and other plants. Recent advancements in detection and sequencing technology have resulted in a considerable rise in research on the biosynthetic processes of aromatic terpenes in Dendrobium and other flowering plants. Nevertheless, the inquiry into the precise means by which plants regulate the proportion of diverse aromatic terpenes in their floral scent, thereby preserving their olfactory traits, requires further investigation. A conjecture on the botanical perfumer mechanism, which condensed the findings of earlier studies, was put forward to address this area of interest. Specific transcription factors likely govern the coordinated expression of multiple key terpene synthase (TPS) genes during the flowering stage of plants, thereby regulating the proportional biosynthesis of diverse aromatic terpenes and sustaining the distinctive aromatic properties of individual plants. This review serves as a significant theoretical reference for further investigations into aromatic volatile compounds in Dendrobium.

Genome-Wide Identification of the YABBY Gene Family in Dendrobium Orchids and Its Expression Patterns in Dendrobium chrysotoxum

The small plant-specific YABBY gene family plays key roles in diverse developmental processes in plants. Dendrobium chrysotoxum, D. huoshanense, and D. nobile are perennial herbaceous plants belonging to Orchidaceae with a high ornamental value. However, the relationships and specific functions of the YABBY genes in the Dendrobium species remain unknown. In this study, six DchYABBYs, nine DhuYABBYs, and nine DnoYABBYs were identified from the genome databases of the three Dendrobium species, which were unevenly distributed on five, eight, and nine chromosomes, respectively. The 24 YABBY genes were classified into four subfamilies (CRC/DL, INO, YAB2, and FIL/YAB3) based on their phylogenetic analysis. A sequence analysis showed that most of the YABBY proteins contained conserved C2C2 zinc-finger and YABBY domains, while a gene structure analysis revealed that 46% of the total YABBY genes contained seven exons and six introns. All the YABBY genes harbored a large number of Methyl Jasmonate responsive elements, as well as anaerobic induction cis-acting elements in the promoter regions. Through a collinearity analysis, one, two, and two segmental duplicated gene pairs were identified in the D. chrysotoxum, D. huoshanense, and D. nobile genomes, respectively. The Ka/Ks values of these five gene pairs were lower than 0.5, indicating that the Dendrobium YABBY genes underwent negative selection. In addition, an expression analysis revealed that DchYABBY2 plays a role in ovary and early-stage petal development, while DchYABBY5 is essential for lip development and DchYABBY6 is crucial for early sepal formation. DchYABBY1 primarily regulates sepals during blooming. Furthermore, there is the potential involvement of DchYABBY2 and DchYABBY5 in gynostemium development. The results of a comprehensive genome-wide study would provide significant clues for future functional investigations and pattern analyses of YABBY genes in different flower parts during flower development in the Dendrobium species.

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.

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.

Molecular and metabolic insights into floral scent biosynthesis during flowering in Dendrobium chrysotoxum

Dendrobium chrysotoxum is considered as an important ornamental dendrobium because of its strong and long-lasting floral scent. Nevertheless, few information is known about the dynamic changes and related formation mechanism of dendrobium floral scent at different flowering stages. In this study, the characteristics and biosynthetic mechanism of floral scent in D. chrysotoxum during flowering was revealed by using widely-targeted volatilomics (WTV) combined with transcriptome analysis. Over 500 kinds of volatile organic compounds (VOCs) were detected in the floral scents of D. chrysotoxum, which improved the knowledge about floral scent components of dendrobium. A total of 153 differential VOCs and 4,487 differentially expressed genes (DEGs) were identified between flowers of different flowering stages, respectively. The results for both volatilomics and transcriptomics data indicated that terpenes and related genes played an important role in the formation of floral characteristics of D. chrysotoxum. But in general, the expression of genes showed an opposite trend to the accumulation of metabolites during flowering, suggesting that the regulation of floral scent biosynthesis might have started at the budding stage in D. chrysotoxum. Additionally, a transcriptional metabolic regulatory network consisting of terpenes, terpene synthases and candidate transcription factors was established. This research is the first systematic and comprehensive exploration of floral characteristics and related mechanisms during flowering in D. chrysotoxum. It provides basis for exploration of mechanisms on the floral scents and the breeding of aromatic dendrobium.

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.

Evolution patterns of NBS genes in the genus Dendrobium and NBS-LRR gene expression in D. officinale by salicylic acid treatment

BACKGROUND

Dendrobium officinale Kimura et Migo, which contains rich polysaccharides, flavonoids and alkaloids, is a Traditional Chinese Medicine (TCM) with important economic benefits, while various pathogens have brought huge losses to its industrialization. NBS gene family is the largest class of plant disease resistance (R) genes, proteins of which are widely distributed in the upstream and downstream of the plant immune systems and are responsible for receiving infection signals and regulating gene expression respectively. It is of great significance for the subsequent disease resistance breeding of D. officinale to identify NBS genes by using the newly published high-quality chromosome-level D. officinale genome.

RESULTS

In this study, a total of 655 NBS genes were uncovered from the genomes of D. officinale, D. nobile, D. chrysotoxum, V. planifolia, A. shenzhenica, P. equestris and A. thaliana. The phylogenetic results of CNL-type protein sequences showed that orchid NBS-LRR genes have significantly degenerated on branches a and b. The Dendrobium NBS gene homology analysis showed that the Dendrobium NBS genes have two obvious characteristics: type changing and NB-ARC domain degeneration. Because the NBS-LRR genes have both NB-ARC and LRR domains, 22 D. officinale NBS-LRR genes were used for subsequent analyses, such as gene structures, conserved motifs, cis-elements and functional annotation analyses. All these results suggested that D. officinale NBS-LRR genes take part in the ETI system, plant hormone signal transduction pathway and Ras signaling pathway. Finally, there were 1,677 DEGs identified from the salicylic acid (SA) treatment transcriptome data of D. officinale. Among them, six NBS-LRR genes (Dof013264, Dof020566, Dof019188, Dof019191, Dof020138 and Dof020707) were significantly up-regulated. However, only Dof020138 was closely related to other pathways from the results of WGCNA, such as pathogen identification pathways, MAPK signaling pathways, plant hormone signal transduction pathways, biosynthetic pathways and energy metabolism pathways.

CONCLUSION

Our results revealed that the NBS gene degenerations are common in the genus Dendrobium, which is the main reason for the diversity of NBS genes, and the NBS-LRR genes generally take part in D. officinale ETI system and signal transduction pathways. In addition, the D. officinale NBS-LRR gene Dof020138, which may have an important breeding value, is indirectly activated by SA in the ETI system.

Triploid cultivars of Cymbidium act as a bridge in the formation of polyploid plants

Triploid is considered a reproductive barrier and also a bridge in the formation of polyploids. However, few reports are available in Cymbidium. In this study, diploid 'Xiaofeng', sexual triploid 'Yuchan' and 'Huanghe' of Cymbidium were used to evaluate hybridization compatibility of the triploids. Results showed that the sexual triploids were fertile whether they were used as male or female parents. 'Yuchan' produced male gametes of 1x, 1x~2x, 2x, 2x~3x, and 3x at frequencies of 8.89%, 77.78%, 6.67%, 3.33%, and 3.33%, respectively; while 'Huanghe' produced 3.33% 1x, 80.00% 1x~2x, 8.89% 2x, 5.56% 2x~3x, and 2.22% 3x male gametes. The cross of 'Xiaofeng' with 'Yuchan' produced progenies with a wide range of ploidy levels, including one diploid, 34 2×~3× aneuploids, 12 triploids, and one tetraploid, indicating that male gametes produced by sexual triploid were fertile and could be transmitted and fused with egg cells. On the other hand, 10 progenies obtained from the cross of 'Yuchan' × 'Xiaofeng' were all aneuploids. The cross of 'Yuchan' with 'Huanghe' produced 40 progenies including three 2×~3× aneuploids, nine 3×~4× aneuploids, 21 tetraploids, six 4×~5× aneuploids, and one pentaploid, suggesting that 2x gametes, instead of the unreduced ones played a more important role in the formation of tetraploids. The survival rates of the hybrids were all above 80.00%, with the tetraploids at 96.67%. Cytological analysis revealed that during meiosis of sexual polyploids, two chromosome sets of the 2n gamete were inclined to enter into the same daughter cell, resulting in the production of 2x gametes. Our results indicate that the triploid cymbidiums are not reproductive barrier but serve as a bridge in the formation of polyploid plants.

Identification of Two GDSL-Type Esterase/Lipase Genes Related to Tissue-Specific Lipolysis in Dendrobium catenatum by Multi-Omics Analysis

Dendrobium catenatum is an important herb and widely cultivated in China. GDSL-Type Esterase/Lipase proteins (GELPs) are widely distributed in plants and play crucial roles in stress responses, plant growth, and development. However, no identification or functional analysis of GELPs was reported in D. catenatum. This study identifies 52 GELPs in D. catenatum genome, which is classified into four groups by phylogenetic analysis. Four conservative blocks (Ser-Gly-Asn-His) are found in most GELP domains. Transcriptome analysis reveals the expression profiles of GELPs in different organs and flowering phases. Co-expression analysis of the transcriptome and lipidome identifies a GELP gene, Dca016600, that positively correlates with 23 lipids. The purified Dca016600 protein shows the optimum pH is active from 8.0 to 8.5, and the optimum temperature is active from 30 °C to 40 °C. The kinetic study provides V_max (233.43 μmol·min^-1·mg^-1) and K_m (1.49 mM) for substrate p-nitrophenyl palmitate (p-NPP). Integrated analysis of the transcriptome and proteome identifies a GELP gene, Dca005399, which is specially induced by freezing. Interestingly, Dca005399 shows high expression in symbiotic germination seeds and sepals. This study provides new insights into the function of D. catenatum GELPs in plant development and stress tolerance.

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.

A chromosome-level, haplotype-phased Vanilla planifolia genome highlights the challenge of partial endoreplication for accurate whole-genome assembly

Vanilla planifolia, the species cultivated to produce one of the world's most popular flavors, is highly prone to partial genome endoreplication, which leads to highly unbalanced DNA content in cells. We report here the first molecular evidence of partial endoreplication at the chromosome scale by the assembly and annotation of an accurate haplotype-phased genome of V. planifolia. Cytogenetic data demonstrated that the diploid genome size is 4.09 Gb, with 16 chromosome pairs, although aneuploid cells are frequently observed. Using PacBio HiFi and optical mapping, we assembled and phased a diploid genome of 3.4 Gb with a scaffold N50 of 1.2 Mb and 59 128 predicted protein-coding genes. The atypical k-mer frequencies and the uneven sequencing depth observed agreed with our expectation of unbalanced genome representation. Sixty-seven percent of the genes were scattered over only 30% of the genome, putatively linking gene-rich regions and the endoreplication phenomenon. By contrast, low-coverage regions (non-endoreplicated) were rich in repeated elements but also contained 33% of the annotated genes. Furthermore, this assembly showed distinct haplotype-specific sequencing depth variation patterns, suggesting complex molecular regulation of endoreplication along the chromosomes. This high-quality, anchored assembly represents 83% of the estimated V. planifolia genome. It provides a significant step toward the elucidation of this complex genome. To support post-genomics efforts, we developed the Vanilla Genome Hub, a user-friendly integrated web portal that enables centralized access to high-throughput genomic and other omics data and interoperable use of bioinformatics tools.

Genome-Wide Identification and Expression Patterns of the SWEET Gene Family in Bletilla striata and its Responses to Low Temperature and Oxidative Stress

SWEETs (sugars will eventually be exported transporters), a well-known class of sugar transporters, are involved in plant growth and development, sugar transport, biotic and abiotic stresses, etc. However, to date, there have been few investigations of SWEETs in Orchidaceae. In this study, 23 SWEET genes were identified in Bletilla striata for the first time, with an MtN3/saliva conserved domain, and were divided into four subgroups by phylogenetic tree. The same subfamily members had similar gene structures and motifs. Multiple cis-elements related to sugar and environmental stresses were found in the promoter region. Further, 21 genes were localized on 11 chromosomes and 2 paralogous pairs were found via intraspecific collinearity analysis. Expression profiling results showed that BsSWEETs were tissue-specific. It also revealed that BsSWEET10 and BsSWEET18 were responsive to low temperature and oxidative stresses. In addition, subcellular localization study indicated that BsSWEET15 and BsSWEET16 were localized in the cell membrane. This study provided important clues for the in-depth elucidation of the sugar transport mechanism of BsSWEET genes and their functional roles in response to abiotic stresses.

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.

First De novo whole genome sequencing and assembly of mutant Dendrobium hybrid cultivar 'Emma White'

The Dendrobium hybrid cultivar 'Emma White' is an ornamental, successfully commercialised orchid. We used a gamma ray-induced early flowering mutant and the Illumina HiSeqX10 sequencing platform to generate the first draft de novo whole genome sequence and assembly. The draft sequence was 678,650,699 bp in length, comprising 447,500 contigs with an N50 of 1423 and 33.48% GC content. Comparing 95,529 predicted genes against the Uniprot database revealed 60,741 potential genes governing molecular functions, biological processes and cellular components. We identified 216,232 simple sequence repeats and 138,856 microsatellite markers. Chromosome-level genome assembly of Dendrobium huoshanense was used to RagTag-scaffold available contigs of the mutant, revealing a total length of 687,254,899 bp with an N50 of 2096. The longest final contiguous length was 18,000,059 bp from 30,571 bp. BUSCO genome completeness was 93.6%. This study is valuable for investigating the mechanisms of mutation, and developing Dendrobium hybrid cultivars using mutation breeding.

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.

Orchid Phylotranscriptomics: The Prospects of Repurposing Multi-Tissue Transcriptomes for Phylogenetic Analysis and Beyond

The Orchidaceae is rivaled only by the Asteraceae as the largest plant family, with the estimated number of species exceeding 25,000 and encompassing more than 700 genera. To gain insights into the mechanisms driving species diversity across both global and local scales, well-supported phylogenies targeting different taxonomic groups and/or geographical regions will be crucial. High-throughput sequencing technologies have revolutionized the field of molecular phylogenetics by simplifying the process of obtaining genome-scale sequence data. Consequently, there has been an explosive growth of such data in public repositories. Here we took advantage of this unprecedented access to transcriptome data from predominantly non-phylogenetic studies to assess if it can be repurposed to gain rapid and accurate phylogenetic insights across the orchids. Exhaustive searches revealed transcriptomic data for more than 100 orchid species spanning 5 subfamilies, 13 tribes, 21 subtribes, and 50 genera that were amendable for exploratory phylotranscriptomic analysis. Next, we performed re-assembly of the transcriptomes before strategic selection of the final samples based on a gene completeness evaluation. Drawing on these data, we report phylogenetic analyses at both deep and shallow evolutionary scales via maximum likelihood and shortcut coalescent species tree methods. In this perspective, we discuss some key outcomes of this study and conclude by highlighting other complementary, albeit rarely explored, insights beyond phylogenetic analysis that repurposed multi-tissue transcriptome can offer.

Characterization of Three SEPALLATA-Like MADS-Box Genes Associated With Floral Development in Paphiopedilum henryanum (Orchidaceae)

Paphiopedilum (Orchidaceae) is one of the world's most popular orchids that is found in tropical and subtropical forests and has an enormous ornamental value. SEPALLATA-like (SEP-like) MADS-box genes are responsible for floral organ specification. In this study, three SEP-like MADS-box genes, PhSEP1, PhSEP2, and PhSEP3, were identified in Paphiopedilum henryanum. These genes were 732-916 bp, with conserved SEPI and SEPII motifs. Phylogenetic analysis revealed that PhSEP genes were evolutionarily closer to the core eudicot SEP3 lineage, whereas none of them belonged to core eudicot SEP1/2/4 clades. PhSEP genes displayed non-ubiquitous expression, which was detectable across all floral organs at all developmental stages of the flower buds. Furthermore, subcellular localization experiments revealed the localization of PhSEP proteins in the nucleus. Yeast two-hybrid assays revealed no self-activation of PhSEPs. The protein-protein interactions revealed that PhSEPs possibly interact with B-class DEFICIENS-like and E-class MADS-box proteins. Our study suggests that the three SEP-like genes may play key roles in flower development in P. henryanum, which will improve our understanding of the roles of the SEP-like MADS-box gene family and provide crucial insights into the mechanisms underlying floral development in orchids.

Natural Composition and Biosynthetic Pathways of Alkaloids in Medicinal Dendrobium Species

Dendrobium is the second biggest genus in the Orchidaceae family, some of which have both ornamental and therapeutic values. Alkaloids are a group of active chemicals found in Dendrobium plants. Dendrobine has emerged specific pharmacological and therapeutic properties. Although Dendrobium alkaloids have been isolated and identified since the 1930s, the composition of alkaloids and their biosynthesis pathways, including metabolic intermediates, alkaloid transporters, concrete genes involved in downstream pathways, and associated gene clusters, have remained unresolved scientific issues. This paper comprehensively reviews currently identified and tentative alkaloids from the aspect of biogenic pathways or metabolic genes uncovered based on the genome annotations. The biosynthesis pathways of each class of alkaloids are highlighted. Moreover, advances of the high-throughput sequencing technologies in the discovery of Dendrobium alkaloid pathways have been addressed. Applications of synthetic biology in large-scale production of alkaloids are also described. This would serve as the basis for further investigation into Dendrobium alkaloids.

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.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.