Abstract
As one of the largest families of angiosperms, the Orchidaceae family is diverse. Dendrobium represents the second largest genus of the Orchidaceae. However, an assembled high-quality genome of species in this genus is lacking. Here, we report a chromosome-scale reference genome of Dendrobium chrysotoxum, an important ornamental and medicinal orchid species. The assembled genome size of D. chrysotoxum was 1.37 Gb, with a contig N50 value of 1.54 Mb. Of the sequences, 95.75% were anchored to 19 pseudochromosomes. There were 30,044 genes predicted in the D. chrysotoxum genome. Two whole-genome polyploidization events occurred in D. chrysotoxum. In terms of the second event, whole-genome duplication (WGD) was also found to have occurred in other Orchidaceae members, which diverged mainly via gene loss immediately after the WGD event occurred; the first duplication was found to have occurred in most monocots (tau event). We identified sugar transporter (SWEET) gene family expansion, which might be related to the abundant medicinal compounds and fleshy stems of D. chrysotoxum. MADS-box genes were identified in D. chrysotoxum, as well as members of TPS and Hsp90 gene families, which are associated with resistance, which may contribute to the adaptive evolution of orchids. We also investigated the interplay among carotenoid, ABA, and ethylene biosynthesis in D. chrysotoxum to elucidate the regulatory mechanisms of the short flowering period of orchids with yellow flowers. The reference D. chrysotoxum genome will provide important insights for further research on medicinal active ingredients and breeding and enhances the understanding of orchid evolution.
Similar content being viewed by others
Introduction
With more than 25,000 species, Orchidaceae is the largest angiosperm family1 and comprises 8–10% of flowering plants. Orchids are renowned for their specialized flowers, which have a very wide variety of growth forms, and have been successful colonizers of a wide variety of different habitats2. As one of the largest genera of Orchidaceae, Dendrobium encompasses ~1450 species with fleshy stems3. Many species of Dendrobium have high medicinal and commercial value, and the main medicinal active ingredients are in the stems4,5,6,7,8,9. Therefore, studying the molecular mechanism of these active ingredients and breeding cultivars with increased contents of natural products are the main objectives in Dendrobium scientific research and industrialization10.
Guchui Shihu (鼓槌石斛) Dendrobium chrysotoxum, a medicinal species, is listed in the Chinese Pharmacopoeia (2020, 2015, and 2010 edition) and contains an abundance of erianin, gigantol, polysaccharides, and fluorenones, among other compounds11,12,13,14,15,16,17,18,19 (Fig. 1). These compounds show antipyretic, analgesic, antihyperglycemic, and antioxidant effects and enhance immune function11,12,13,14,15,16,17,18,19. Recently, preliminary clinical study results suggested that gigantol could delay lens turbidity through the inhibition of aldose reductase and aldose reductase mRNA expression, which have good effects on diabetic cataracts16,17. Erianin has been demonstrated to exhibit metabolic inhibition20 and antitumor21, antiproliferative22, and antiangiogenic activity23. Moreover, it also inhibits high glucose-induced retinal angiogenesis12. Polysaccharides isolated from D. chrysotoxum have potential utility in enhancing antioxidation, immune function, and/or hypoglycemic activity11. The stems of D. chrysotoxum are fusiform and rich in active medicinal substances, which makes it a suitable species for scientific research and industrial applications in Dendrobium.
D. chrysotoxum is lithophytic on rocks or epiphytic on tree trunks with bright yellow flowers and fusiform, fleshy stems
With the improvement of sequencing technology and cost reduction, genome sequencing has become a necessary method for obtaining comprehensive genetic information and an effective method for screening candidate genes for specific traits, especially for identifying candidate genes involved in the biosynthesis pathways of medicinal compounds24,25,26,27,28. To date, only two Dendrobium spp. genomes have been sequenced, and some candidate genes involved in polysaccharide metabolic pathways have been identified in those two species24,29,30. However, these studies were largely limited due to their low-quality genome assemblies. Therefore, high-quality reference genomes and additional Dendrobium species need to be sequenced to better understand the molecular mechanisms underlying the production of medicinal compounds and enable the breeding of new varieties.
In this study, we used PacBio sequencing and Hi-C technology to generate a chromosome-level genome assembly. The specific genes of D. chrysotoxum were identified, which lays a foundation for further research on the functions of medicinal active ingredients, provides a reference for breeding new varieties and enhances the understanding of orchid evolution.
Results and discussion
Genome sequencing and characteristics
D. chrysotoxum has a karyotype of 2N = 2X = 38, with uniform chromosomes31. To completely sequence the D. chrysotoxum genome, 138.15 Gb of clean reads were generated by BGISEQ sequencing system (Supplementary Table 1). The estimated genome size was 1.38 Gb with 1.84% heterozygosity, as determined by K-mer analysis (Supplementary Fig. 1). To obtain a better assembly, PacBio technology was employed, and 132.64 Gb of PacBio sequencing data were generated (Supplementary Table 1). The assembly size was 1.37 Gb with a corresponding contig N50 value of 1.54 Mb (Supplementary Table 2). The BUSCO32 assessment indicated that the completeness of the gene set of the assembled genome was 90.3% (Supplementary Table 3). This indicates that the D. chrysotoxum genome assembly was complete and could be used for subsequent analysis. We further used 125.96 Gb of reads from the Hi-C library. The assembled scaffolds were ultimately clustered into 19 pseudomolecules, which represented the 19 chromosomes in the haploid genome of D. chrysotoxum (Fig. 2a). The lengths of the 19 pseudochromosomes ranged from 38.28 to 100.49 Mb with a scaffold N50 value of 67.80 Mb (Supplementary Tables 4 and 5). In addition, contigs with a length of 1.31 Gb were mapped onto the 19 pseudochromosomes at a 95.75% anchor rate (Supplementary Tables 4 and 5). The chromatin interaction data suggest that our Hi-C assembly is of high quality (Fig. 2b). Compared with those of other orchid genome assemblies, the contig N50 and scaffold N50 values of the D. chrysotoxum genome were much higher (Table 1), and the assembly completeness was higher than 90% (Table 1), suggesting high genome quality and completeness.
a From inside outward: chromosome (purple), gene density (red), DNA type repeat sequence density (green), copy density (blue), and gypsy density (orange). All the data are shown with sliding windows of 500 kb, and the inner lines (green indicates the positive direction, and red indicates the opposite direction) represent syntenic blocks on homologous chromosomes. b Heat map of the intensity of the Hi-C chromosome. The heat map represents the contact matrices generated by aligning the Hi-C data to the chromosome-scale assembly of the D. chrysotoxum genome. A higher value on the scale bar indicates a higher contact frequency
Gene prediction and annotation
In D. chrysotoxum genome, 30,044 protein-coding genes were annotated (see Materials and methods; Supplementary Table 6). The completeness of the genome was 95.64%, indicating that the D. chrysotoxum genome annotation was relatively complete (Supplementary Table 7).
In addition to a high number of genes, the average length of genes and introns was also larger in D. chrysotoxum than in Phalaenopsis equestris, Gastrodia elata, and D. catenatum24,33,34 and much higher than that in most other angiosperms (Supplementary Table 8). The average length of the coding DNA sequences (CDSs) in D. chrysotoxum was longer than those in other angiosperms, and a greater average intron length was also previously observed for P. equestris, G. elata, and D. catenatum24,33,34; thus, a relatively long CDS might be a unique characteristic of Orchidaceae (Supplementary Fig. 2; Supplementary Table 8). Regulatory elements are frequently present in introns, and alternative splicing events often occur among different introns and exons, diversifying the protein-coding aspect of the genome. All these factors might contribute to genome structure evolution, genome size, gene function diversification, and gene expression patterns35,36,37,38. For example, intron transcriptional delay in Drosophila is particularly important for proper development of the embryo39,40. Thus, this characteristic of orchids needs to be further analyzed and researched. Moreover, 80 microRNAs, 1281 transfer RNAs, 2275 ribosomal RNAs, and 882 small nuclear RNAs were identified in the D. chrysotoxum genome (Supplementary Table 9).
We estimated that the D. chrysotoxum genome comprised 62.81% repetitive sequences (Supplementary Figs. 3 and 4; Supplementary Table 10), the percentage of which was higher than 62% in P. equestris but lower than 78.1% in D. catenatum24,33. Transposable elements (TEs) are important forms of repeats and constitute a substantial part of the D. chrysotoxum genome (61.22%); TEs are the most abundant repeat subtypes in this species. In addition, repeats predicted de novo were much larger than those obtained based on Repbase11 database, suggesting that, compared with other plants species whose genome has been sequenced, D. chrysotoxum has many specific repeats (Supplementary Table 10). Long terminal repeats (LTRs) represented the highest proportion among all subtypes of repeats, accounting for ~53.15% of the genome, which was higher than the 46% for D. catenatum24 (Supplementary Table 11).
In addition, 27,575 (91.78%) predicted genes were functionally annotated (Supplementary Table 12). Among them, 27,268 (90.76%) and 26,808 (89.23%) genes were annotated to the TrEMBL and Nr databases, respectively (Supplementary Table 12). The numbers of annotated genes were 22,735 (75.67%), 19,185 (63.86%), and 18,666 (62.13%) in the InterPro, SwissProt, and KEGG databases, respectively (Supplementary Table 12).
Evolution of gene families
A high-confidence phylogenetic tree was constructed, and the divergence times were estimated based on 274 single-copy genes from 17 different plant species (Supplementary Fig. 5 and Supplementary Table 6). As expected, D. chrysotoxum was sister to D. catenatum, forming an Epidendroideae clade together with P. equestris, G. elata, and A. shenzhenica located at the bases of Orchidaceae branches (Supplementary Fig. 6). The Orchidaceae divergence was estimated to have occurred 123 Mya; the divergence between subfamily Apostasioideae and subfamily Epidendroideae occurred 80 Mya; the divergence between D. chrysotoxum and D. catenatum occurred 11 Mya; and the divergence between Dendrobium and Phalaenopsis occurred 38 Mya (Fig. 3). Then, the expansion and contraction of orthologous gene families were analyzed. According to the results, 140 and 1112 gene families expanded and contracted, respectively, in the lineage leading to Orchidaceae. In D. chrysotoxum, 953 gene families were expanded, as opposed to 783 in D. catenatum, 853 in P. equestris, 358 in G. elata, and 562 in A. shenzhenica. At the same time, 1335 gene families were contracted in D. chrysotoxum, as opposed to 644 in D. catenatum, 1009 in P. equestris, 2748 in G. elata, and 1615 in A. shenzhenica. A greater number of expanded gene families in D. chrysotoxum may lead to a larger genome size than that in other sequenced orchid species2,24,33,34.
The green and red numbers represent the numbers of expanded gene families and contracted gene families, respectively. The blue portions of the pie charts show that the copy numbers of gene families are constant. Divergence times are represented by light blue bars at internodes; the divergence time is at a 95% confidence interval and shown by the range of these bars. The expansion and contraction of gene families are represented by members of branches (see Materials and methods and Supplementary Fig. 6). The orange part of the pie chart at the top left represents the ratio of 11,252 gene families found in the most recent common ancestor (MRCA) that expanded or contracted during recent differentiation events
The ancestral clade of Dendrobium had 464 expanded gene families and 216 contracted gene families. The D. chrysotoxum clade had 953 expanded gene families and 1335 contracted gene families. In the ancestral clade of Dendrobium, there were 19 significantly expanded gene families, including 236 genes from D. chrysotoxum. In the D. chrysotoxum clade, 107 gene families were significantly expanded, including 1048 genes, and 43 gene families were significantly contracted, including 59 genes. We also conducted Gene Ontology (GO) enrichment analysis for the expanded gene families, and the GO terms “cytoplasmic part” and “intracellular organelle” were found to be enriched (Supplementary Table 13). In addition, the bidirectional sugar transporter gene SWEET was identified (Supplementary Fig. 7), whose product plays important roles in sugar translocation between compartments41, phloem loading for long-distance translocation42, pollen nutrition43, and seed filling44. Further phylogenetic analysis showed that 17 genes were expanded in clade II (Supplementary Fig. 7), suggesting that these SWEET genes might be associated with a fleshy stem that is abundant in polysaccharides and other medicinal compounds.
Synteny analysis and whole-genome duplication (WGD)
Both the loss of a substantial fraction of genes and the increase in substitution rate complications were indicated by WGD in D. chrysotoxum, which is thought to have occurred among different orchid species2. WGD is evident in many lineages and is a practical method for genome expansion45. To determine the occurrence of WGDs in D. chrysotoxum, JCVI v0.9.1446 was used to analyze the protein sequences of D. chrysotoxum, P. equestris, P. aphrodite, and D. catenatum with the default parameters and obtain collinear gene pairs. There were 21,881 collinear gene pairs between D. chrysotoxum and P. equestris, 21,592 between D. chrysotoxum and P. aphrodite, 24,550 between D. chrysotoxum and D. catenatum, and 2800 between D. chrysotoxum and itself (Supplementary Table 14). Although D. chrysotoxum was assembled to the chromosome level, its self-collinearity was still very low compared to that of other sequenced orchid species. The collinearity between D. catenatum and D. chrysotoxum was fragmented, which may be the result of the quality of the D. catenatum genome, which was not at the chromosome level. The chromosomes of Dendrobium and Phalaenopsis showed a good corresponding relationship, indicating that after the divergence of Dendrobium and Phalaenopsis, the chromosomes were conserved, with few rearrangements. Syntenic figures show that the collinearity blocks were mainly in a 1:1 pattern, indicating that after the differentiation of D. chrysotoxum, no species-specific WGD events had occurred (Fig. 4; Supplementary Figs. 8–12).
The values on the X- and Y-axes are the numbers of cumulative genes on the 19 chromosomes
The distributions of synonymous substitutions per synonymous site (Ks) were estimated to infer polyploidization events that occurred in the D. chrysotoxum genome. There were two peaks in the distribution of Ks for paralogous D. chrysotoxum genes: Ks = 1.0 and 1.7–1.8 (Fig. 5). These results suggested that two polyploidization events occurred in D. chrysotoxum. To further verify the polyploidization events in D. chrysotoxum, its genome was compared with that of P. equestris, A. shenzhenica, and D. catenatum. The peaks in Ks values between both D. chrysotoxum/P. equestris and D. chrysotoxum/D. catenatum were less than 1.0, suggesting that the events occurred before the differentiation of these three species. There was a diverging peak in the Ks distribution of D. chrysotoxum and A. shenzhenica at Ks = 0.7–0.8, which was smaller than but close to the Ks peak of the Orchidaceae (Ks = 1), indicating that extant orchid species differentiated immediately after experiencing a shared WGD event. Based on the evolution of gene families, species differentiation mainly occurred through gene loss with little gene expansion, which confirmed that the WGD event occurred in the most recent common ancestor of extant orchid species. The second peak in the Ks distributions within D. chrysotoxum (1.7–1.8) indicated that the τ WGD had occurred in most monocot species45. Furthermore, the peak of the Ks distribution in D. chrysotoxum was smaller than 0.2, suggesting that it originated from background (tandem) duplications and likely did not signify additional recent WGDs2. Therefore, this study found that D. chrysotoxum experienced two polyploidization events: an early WGD event was shared among all extant orchid species, and a later event that was shared among most monocot species.
The Ks distributions of paralogs using the reciprocal best hit (RBH) and anchor are shown as light blue and red histograms, respectively. The Ks distribution for the one-to-one orthologs of D. chrysotoxum–G. elata, D. chrysotoxum–A. shenzhenica, D. chrysotoxum–D. catenatum, and D. chrysotoxum–P. equestris are shown as green, gray, blue, and cyan curves, respectively. RBH reciprocal best blast hit; “anchor” refers to colinear regions
MADS-box genes and the evolution of flowers
MADS-box genes are among the most important regulators of plant floral development and compose major class of regulators mediating floral transition. In total, the D. chrysotoxum genome encodes 58 putative functional MADS-box genes and 1 pseudogene (Table 2; Supplementary Table 15). Interestingly, the number of MADS-box genes was similar to that in other sequenced orchid species but smaller than that in most sequenced angiosperms2,24,33. D. chrysotoxum has 31 type II MADS-box genes, which is higher than that found in P. equestris (29) and A. shenzhenica (27), but smaller than that of D. catenatum (35)2,24,33. Phylogenetic analysis (Supplementary Fig. 13) showed that, except for those in the MIKC*, Bs, and OsMADS32 clades, most genes in the type II MADS-box clade were contracted. Bs genes are involved in the differentiation and development of ovules47. In D. chrysotoxum, there are four Bs members, more than the number found in other sequenced orchid species. The Bs genes had duplicated, as evidenced by higher seed production in D. chrysotoxum than in other sequenced orchid species. This must have been accompanied by duplication of the type I MADS-box gene Mα, as D. chrysotoxum has more Mα genes (19) than other sequenced orchid species (Table 2), ensuring seed development. In addition, there were no genes from the FLOWERING LOCUS C (FLC), AGL12, or AGL15 clades in the D. chrysotoxum genome or other sequenced orchid genomes. In comparison with genes in the AGL12 and AGL15 clades, which are present in both rice and Arabidopsis, orthologous genes of FLC, AGL12, and AGL15 might have been specifically lost in orchids. Although AGL12-like genes (XAL1 in A. thaliana) are necessary for root development and flowering48, D. chrysotoxum and P. equestris have varying mechanisms that perform the same function2, showing that D. chrysotoxum is not a terrestrial orchid but is an epiphytic orchid.
The D. chrysotoxum genome has 26 putative functional type I genes and 1 pseudogene (Table 2), which might have resulted in a lower expansion rate or a higher contraction rate compared with those of type II MADS-box genes in D. chrysotoxum (31 functional genes). Tandem gene duplication might play an important role in the increasing number of type I genes in the α group (type I Mα), suggesting that the type I genes have mainly been duplicated on a smaller scale from more-recent duplications49. Although members of the β group of type I MADS-box genes (type I Mβ) do exist in A. thaliana, poplar, and rice, they are absent in the D. chrysotoxum genome. Interactions among these type I MADS-box genes are essential for initiating endosperm development50; therefore, like in other sequenced orchids2,24,33, endosperm is also absent in D. chrysotoxum.
Floral color regulatory pathway in D. chrysotoxum
The flowering time of a single flower of D. chrysotoxum was ~10 days51,52, the limit of which might be associated with yellow flower color. All photosynthetic tissues in each of the biological kingdoms can produce carotenoids53. More than 1100 naturally occurring carotenoids (http://carotenoiddb.jp/) are involved in many of the red, orange, and yellow colors of flowers53. These compounds also play important roles in photosynthesis. Interestingly, carotenoids function as precursors for the biosynthesis of abscisic acid (ABA)53. Moreover, ethylene plays a role in senescing flowers54. Ethylene and ABA regulate plant growth and development55 synergistically or antagonistically. We therefore analyzed the network involving carotenoid, ABA, and ethylene biosynthesis and regulation (Fig. 6).
Increased gene expression during the four developmental stages of flowers is shown by thick green arrows, while decreased gene expression is shown by thick red arrows. The enzymes or genes are indicated alongside the arrows. The dashed lines suggest that there are multiple steps. Carotenes are shown in orange, and light gold indicates xanthophylls. ABA-related genes are shown in yellow, and ethylene-related genes are shown in light green. CCDs (orange arrows) regulate carotenoid accumulation. Plastoglobule-localized metallopeptidase 48 (PGM48) is hypothesized to regulate CCD4. The orange protein regulates PSY activity, acting as a chaperone. ClpB3 (red arrow) regulates the activity of DXS, acting as an enzyme in the MEP pathway
Eighteen genes or gene family members in the carotenoid biosynthesis pathway and related regulatory mechanisms were identified (Supplementary Table 16). These genes encode phytoene synthase (PSY), orange protein, casein lytic proteinase B3 (ClpB3), deoxy-D-xylulose 5-phosphate synthase (DXS), phytoene desaturase, ζ-carotene isomerase, ζ-carotene desaturase, carotenoid isomerase (CRTISO), β-lycopene cyclase, ε-lycopene cyclase (LYCE), β-carotene hydroxylase, carotene ε-hydroxylase (CYP97C), zeaxanthin epoxidase (ZEP), violaxanthin de-epoxidase, neoxanthin synthase (NXS), xanthophyll acyl-transferase (XAT), plastoglobule-localized metallopeptidase 48, and carotenoid cleavage dioxygenase (CCD). The expression of these genes increased with flower development, except for LYCE, which is targeted for downregulation during biofortification, ZEP, NXS, and XAT (Supplementary Table 16; Fig. 6), suggesting that more carotenoids and fewer xanthophylls were produced during flowering to senescence.
The substrates used to produce ABA were neoxanthin and violaxanthin, and the process was regulated by nine-cis-epoxy carotenoid dioxygenases (NCEDs). The biosynthesis of ABA is catalyzed by the short-chain dehydrogenase/reductase-like (SDR1) enzyme abscisic aldehyde oxidase (AAO) and molybdenum cofactor (MoCo). The expression of NCED and SDR genes increased gradually with the development of flowers, while the expression of AAO and MoCo genes decreased gradually (Supplementary Table 17; Fig. 6). Furthermore, there were four AAO gene members detected in Arabidopsis, while there was only one gene detected in D. chrysotoxum (Supplementary Fig. 14). Taken together, these findings might indicate that relatively low amounts of ABA (C15) were produced, which might improve ethylene biogenesis.
For ethylene biogenesis, genes encoding three kinds of enzymes were identified. The expression of Maker79017, encoding S-AdoMet synthetase (SAMS), Maker75695 and Maker66290, encoding ACC synthase (ACS), and Maker29641, encoding ACC oxidase (ACO), increased gradually, suggesting that increased amounts of ethylene were produced during the development of flowers (Supplementary Table 18; Fig. 6). CONSTITUTIVE TRIPLE RESPONSE 1 and ETHYLENE INSENSITIVE 2 regulate the interaction between ethylene and the ABA pathway and are partially dependent on the MHZ5/CRTISO-mediated ABA pathway in rice55. Therefore, we also analyzed the expression patterns of the two genes in D. chrysotoxum, but there were no obvious differences in any of the four stages of flower development (Supplementary Table 18; Fig. 6).
In conclusion, carotenoid production increased gradually, and the content of xanthophylls decreased gradually in yellow D. chrysotoxum flowers during flowering to senescence. Less xanthophyll was degraded into less ABA, and less ABA led to more ethylene being produced. As a result, yellow flowers of D. chrysotoxum generally have a relatively short flowering period.
Identification of the terpene synthases (TPS) and Hsp90 gene families and adaptive evolution
Dendrobium spp., with epiphytic or lithophytic lifestyles, frequently experience adverse environmental conditions, such as chilling and water deficit56. During plant responses to environmental stresses, volatile terpenes play critical roles56. Moreover, terpenes also play an important role in the formation of orchid floral scents56. TPSs are the key enzymes involved in terpene biosynthesis57. Different sizes of TPS families and subfamilies in plant species have evolved to synthesize a specific set of terpene compounds57. There are seven subfamilies in the TPS family: TPS-a, TPS-b, TPS-c, TPS-d, TPS-e/f, TPS-g, and TPS-h57. Among them, TPS-a encodes a sesquiterpene synthase that is found in both dicotyledonous and monocotyledonous plants. Angiosperm-specific TPS-b encodes a monoterpene synthase with an R(R)X8W motif that catalyzes the isomerization cyclization reaction. TPS-c belongs to the ancestral clade and catalyzes the activity of copalyl diphosphate synthase. Gymnosperm-specific TPS-d performs several functions, such as those of diterpene, monoterpene, and sesquiterpene synthases. TPS-e/f encodes copalyl diphosphate/kaurene synthases, which are critical enzymes for gibberellic acid production. Another angiosperm-specific TPS, TPS-g, encodes monoterpene synthase enzymes that lack the R(R)X8W motif. TPS-h has been observed only in Selaginella moellendorffii58,59,60,61. Phylogenetic analysis of the TPS gene family members and their expression in bud formation and initial flower opening, blooming, and withering are shown in Fig. 7. In this study, the TPS gene number in D. chrysotoxum was 48, which was greater than that in D. catenatum (42) (Fig. 7a). Moreover, there were 14 and 21 genes in A. shenzhenica and P. equestris, respectively. The TPS-b subfamily can be divided into monocot and eudicot clades. More D. chrysotoxum TPS genes than D. catenatum ones clustered in the monocot A clade—14 (red gene ID) and 4 (blue gene ID), respectively (Fig. 7a). Fewer TPS genes were found in D. chrysotoxum than in D. catenatum in the monocot B clade—7 (red gene ID) and 10 (blue gene ID), respectively (Fig. 7a). The different distribution patterns might contribute to the difference in terpenoid compositions between these two species, which needs further validation.
a Phylogenetic analysis of TPS genes in D. chrysotoxum, D. catenatum, A. shenzhenica, A. thaliana, and P. equestris. Ash A. shenzhenica, Maker D. chrysotoxum, Dca D. catenatum, Peq P. equestris, AT A. thaliana. b Expression patterns of TPS genes in buds and in the initial flower opening, blooming, and withering stages of D. chrysotoxum
To explore heat stress-related genes in D. chrysotoxum, we also analyzed heat stress-related gene families across orchid species. Only two Hsp90 genes (red gene ID) were identified (clustering in group III), with high expression during bud formation (Supplementary Fig. 15a, b). This number was lower than that for the other four species (six were identified in D. catenatum, seven in P. equestris, six in A. shenzhenica, and seven in A. thaliana). This large gene loss might be related to resistance to heat stress.
Conclusion
Although D. chrysotoxum has high ornamental and medicinal value, further molecular mechanism research and development of medicinal compounds have been limited by a lack of omics data. In this study, a chromosome-level reference genome of D. chrysotoxum with an assembled genome size of 1.37 Gb and 30,044 annotated protein-coding genes was obtained. Ks analysis suggested that two polyploidization events occurred in D. chrysotoxum: a recent WGD shared among other orchid species and an ancient polyploidization event shared among most monocots (τ event). Phylogenetic analysis of the SWEET gene family in D. chrysotoxum showed that gene expansion occurred in clade II of the SWEET gene family, which might be related to fleshy stems containing an abundance of polysaccharides. Floral color regulation analysis showed that fewer xanthophylls degraded into ABA, which led to more ethylene production, thus accelerating the senescence of D. chrysotoxum flowers. The analysis of D. chrysotoxum helped elucidate the mechanism through which fleshy stems produce an abundance of polysaccharides and other medicinal compounds, as well as flowering time regulation, which is critical for industrial development. Our results provide the first high-quality genome of Dendrobium and give important insights into the molecular mechanism underlying the production of medicinal active ingredients, breeding, and orchid evolution.
Materials and methods
DNA preparation and sequencing
Fresh leaves of wild D. chrysotoxum were collected for genome sequencing. A modified cetyltrimethylammonium bromide protocol was used to extract the genomic DNA. To estimate genome size and heterozygosity, 143.78 Gb of raw data from paired-end libraries (PE150) constructed from a MGISEQ-2000 sequencer were generated. After data filtering was carried out by SOAPnuke v1.6.5 software with the parameters -n 0.02 -l 20 -q 0.4 -Q 2 -i -G --seqType 0 –rmdup, clean data (138.15 Gb) were obtained (Supplementary Table 1). Then, a SMRTbell Template Prep Kit 1.0 (PacBio, Menlo Park, CA, USA) and a PacBio Sequel system were used to construct and sequence the DNA libraries, respectively, for PacBio long-read sequencing. A total of 132.64 Gb of sequencing data (coverage of 96.12%) were generated, with an N50 read length of 19.5 kb (Supplementary Table 1). Furthermore, all libraries with a 500 bp insert size were sequenced on a NovaSeq platform (2 × 150 bp). We ultimately produced 169.25 Gb of data and 125.96 Gb of clean data for Hi-C analysis. The transcriptomes of flowers of D. chrysotoxum were obtained from Huang’s doctoral thesis62 to assist gene annotation.
Genome assembly
Genome size and heterozygosity were measured using Jellyfish v.2.2.6 and GenomeScope (http://qb.cshl.edu/genomescope)63 based on a 17-K-mer distribution. Canu64 was used to assemble the PacBio sequencing reads, with the following parameters: minOverlapLength = 700; minReadLength = 1000; and corOutCoverage = 50. Then, Arrow software was used to polish the assembly, and Pilon v1.2365 was further used for correction of the assembly based on short reads, with the following parameters: fixed bases; mindepth 10; minqual 20; and diploid. Finally, the completeness and quality of the final assembled genome were evaluated with BUSCO v332.
Hi-C library construction and chromosome assembly
The raw reads produced by the NovaSeq sequencing platform were filtered by SOAPnuke66 (v1.6.5, https://github.com/BGI-flexlab/SOAPnuke) software with the following parameters: -n 0.02 -l 20 -q 0.4 -i –rmdup. Then, the obtained clean reads were compared with the preassembled contigs using Juicer67 software. After filtering the results and removing the misaligned reads, 3D-DNA68 software was used to preliminarily cluster, sequence, and direct the pseudochromosomes. Juicer-box was used to adjust, reset, and cluster the pseudochromosomes to improve the chromosome assembly quality. For the evaluation of the Hi-C assembly results, the final pseudochromosome assemblies were divided into 100 kb bins of equal lengths, and a heat map was used to visualize the interaction signals generated by the valid mapped read pairs between each bin.
Genome annotation
Repetitive sequences are an important part of a genome and are divided into two types, namely, tandem repeats and interspersed repeats. Two methods, de novo prediction and homology-based searches, were used to annotate repeat sequences in the genome. RepeatMasker v4.0.7 and RepeatProteinMask v4.0.7 software69 (http://www.repeatmasker.org) were used to identify repetitive sequences based on the Repbase v21.12 database69 (http://www.girinst.org/repbase). For de novo prediction, a repetitive sequence database was constructed using RepeatModeler v1.0.870 and LTR_FINDER v1.0671. RepeatMasker software and Tandem Repeats Finder v4.0972 were subsequently used to predict repeat sequences and identify tandem repeats in the genome, respectively. The annotation of high-quality protein-coding genes was carried out by integrating homology-based, de novo and transcriptome-based predictions. For homology-based prediction, protein sequences from six species (Arabidopsis thaliana, Oryza sativa, Sorghum bicolor, Zea mays, G. elata, and P. equestris) were used to align D. chrysotoxum genome sequences via Exonerate v2.2.073. Then, the complete sequences of 3000 genes from the homology-based prediction method were used to produce a training model through Augustus v3.2.374 and SNAP v2006-07-2875 software. The RNA-seq data of D. chrysotoxum were mapped to genome sequences through HISAT2 and StringTie software76,77. Finally, Maker v2.31.878 was used to annotate and integrate the results generated by the above methods. BUSCO v332 was then used to evaluate the completeness and quality of the gene models.
Functional annotation of the predicted gene models was carried out by BLAST v2.2.3179 software and aligned against the contents of the SwissProt80, TrEMBL (http://www.uniprot.org/), KEGG (http://www.genome.jp/kegg/), InterPro81, Nr, and GO (The Gene Ontology Consortium) databases82. For noncoding RNA annotations, tRNAscan-SE 1.3.1 (http://lowelab.ucsc.edu/tRNAscan-SE/)83 was used to annotate tRNA sequences. BLASTN79 was used to search for rRNA, and miRNA and snRNA sequences were predicted by Infernal 1.1 (http://infernal.janelia.org/) software84.
WGD analysis
Ks distribution analysis was used to infer the occurrence of WGD events in D. chrysotoxum and those between D. chrysotoxum and A. shenzhenica, D. catenatum, and P. equestris. BLASTP79 was used to search for putative paralogous and orthologous genes within and between genomes by alignment of each genome pair. MCScanX v1.5.285 was used to identify colinear regions, and then CodeML in the PAML package86 was used to calculate the Ks value of each salicoid duplicated gene pair. We used CAFÉ87 to evaluate the significance of each expanded and contracted gene family (P < 0.01).
SWEET gene family analyses
To identify SWEET proteins, proteomic datasets of four orchid species (D. chrysotoxum, A. shenzhenica, D. catenatum, and P. equestris) and A. thaliana were constructed. The MtN3_slv domain PF03083 model profile from the Pfam database88 was used for performing local searches of proteome datasets containing five species via the HMMER program89. The SWEET protein sequences were aligned with MAFFT90. The alignment was then used for phylogenetic tree reconstruction by PhyML 3.091,92,93 with the default parameters.
MADS-box gene family analysis
The sequences of the MADS-box proteins of A. thaliana and the HMM profile (PF00319) were obtained from the Arabidopsis information resource (TAIR) (https://www.arabidopsis.org/) and the Pfam database88, respectively. Then, the sequences of the MADS-box gene family members in the D. chrysotoxum genome were obtained using HMMER 3.2.1 software89 and BLASTP83 methods. The obtained amino acid sequences were used for TBLASTN79 analysis of the D. chrysotoxum transcriptomic assemblies. SMART94 was subsequently used to confirm the obtained sequences by domain analysis. MEGA X95 was then used for the alignment of the candidate genes, and the CIPRES website (https://www.phylo.org/portal2/) was used for phylogenetic tree construction. iTOL (https://itol.embl.de) was subsequently used to visualize the phylogenetic trees.
Identification of genes involved in the carotenoid, ABA, and ethylene biosynthesis pathways and regulatory mechanisms in D. chrysotoxum
The sequences of all 17 genes or gene family members involved in the carotenoid biosynthesis pathway and regulatory mechanisms in A. thaliana, Triticum aestivum, and Pantoea ananatis53 were used as queries to search against the protein database of D. chrysotoxum. The obtained amino acid sequences were aligned using MAFFT90. We then manually inspected the aligned sequences and removed any obviously inconsistent sequences.
Four genes or gene family members involved in the ABA biosynthesis pathway or regulatory mechanisms in Arabidopsis were obtained96. BLASTP79 was used to search for homologous genes by querying the protein database of D. chrysotoxum. After aligning the amino acid sequences with MAFFT90 software, we removed any obviously inconsistent sequences.
The sequences of genes encoding SAMS, ACS, and ACO, all of which are involved in the ethylene biosynthesis pathway, in Arabidopsis97 were used as queries for searching proteins by BLASTP79 software.
For gene families, a phylogenetic tree was constructed with PhyML98 based on the alignment of sequences from D. chrysotoxum, D. catenatum, A. shenzhenica, P. equestris, and A. thaliana. The tree was generated by the maximum likelihood method based on the Jones–Taylor–Thornton (JTT) matrix-based model99, and the fast likelihood-based method was used for phylogenetic tests with SH-like branch supports.
Gene expression analysis
Transcriptome data from flowers at four developmental stages (flower buds, initial flowering stage, blooming period, and withering flowers), stems, and leaves were obtained (BioProject PRJNA691441), and Salmon v1.3.0100 was used to quantify gene expression, with the default settings.
TPS and Hsp90 gene family identification
The HMM profiles for PF01397 (Terpene_synth) and PF03936 (Terpene_synth_C) were downloaded from the Pfam database (pfam.xfam.org/), and both profiles were used to carry out HMM searches against the information of the protein databases for five species (D. chrysotoxum, Dendrobium catenatum, P. equestris, Apostasia shenzhenica, and A. thaliana). The sequences aligned with MAFFT90 were used for phylogenetic tree construction through PhyML79. The tree was generated by the maximum likelihood method based on the JTT matrix-based model99 and the bootstrap method for phylogenetic tests with 1000 replications. Similarly, the HMM profile for PF00183 (Hsp90) was downloaded from the Pfam database (pfam.xfam.org/), and the subsequent steps were the same as those for TPS gene family identification.
Data availability
All the data from this study have been deposited in the NCBI database under BioProject ID PRJNA664445.
References
Christenhusz, M. J. M. & Byng, J. W. The number of known plants species in the world and its annual increase. Phytotaxa 261, 201–217 (2016).
Zhang, G. et al. The Apostasia genome and the evolution of orchid. Nature 549, 379–383 (2017).
Pridgeon, A. M. et al. Genera Orchidacearum: Epidendroideae (Part three) (Oxford University, 2014).
Luo, J.-P., Deng, Y.-Y. & Zha, X.-Q. Mechanism of polysaccharides from Dendrobium huoshanense. on streptozotocin-induced diabetic cataract. Pharm. Biol. 46, 243–249 (2008).
Leitch, I. J. et al. Genome size diversity in orchids: consequences and evolution. Ann. Bot. 104, 469–481 (2009).
Ng, T. B. et al. Review of research on Dendrobium, a prized folk medicine. Appl. Microbiol. Biotechnol. 93, 1795–1803 (2012).
Pan, L. H. et al. Comparison of hypoglycemic and antioxidative effects of polysaccharides from four different Dendrobium species. Int. J. Biol. Macromol. 64, 420–427 (2014).
Tian, C.-C., Zha, X.-Q. & Luo, J.-P. A polysaccharide from Dendrobium huoshanense prevents hepatic inflammatory response caused by carbon tetrachloride. Biotechnol. Biotechnol. Equip. 29, 132–138 (2014).
Huang, K. et al. Purification, characterization and biological activity of polysaccharides from Dendrobium officinale. Molecules 21, 701 (2016).
Lu, J. et al. High-density genetic map construction and stem total polysaccharide content-related QTL exploration for Chinese endemic Dendrobium (Orchidaceae). Front. Plant Sci. 9, 398 (2018).
Zhao, Y. et al. Antioxidant and anti-hyperglycemic activity of polysaccharide isolated from Dendrobium chrysotoxum Lindl. J. Biochem. Mol. Biol. 40, 670–677 (2007).
Li, S. et al. Elution-extrusion counter-current chromatography separation of five bioactive compounds from Dendrobium chrysototxum Lindl. J. Chromatogr. A 1218, 3124–3128 (2011).
Hu, J., Fan, W., Dong, F., Miao, Z. & Zhou, J. Chemical components of Dendrobium chrysotoxum. Chin. J. Chem. 30, 1327–1330 (2012).
Yu, Z. et al. Dendrobium chrysotoxum Lindl. alleviates diabetic retinopathy by preventing retinal inflammation and tight junction protein decrease. J. Diabetes Res. 2015, 518317 (2015).
Yu, Z. et al. Erianin inhibits high glucose-induced retinal angiogenesis via blocking ERK1/2-regulated HIF-1alpha-VEGF/VEGFR2 signaling pathway. Sci. Rep. 6, 34306 (2016).
Wu, J. et al. Gigantol from Dendrobium chrysotoxum Lindl. binds and inhibits aldose reductase gene to exert its anti-cataract activity: an in vitro mechanistic study. J. Ethnopharmacol. 198, 255–261 (2017).
Lim, V., Schneider, E., Wu, H. & Pang, I. H. Cataract preventive role of isolated phytoconstituents: findings from adecade of research. Nutrients 10, 1580 (2018).
Ren, Z. Y. et al. Functional analysis of a novel C-glycosyltransferase in the orchid Dendrobium catenatum. Hortic. Res. 7, 111 (2020).
Robles-Rivera, R. R. et al. Adjuvant therapies in diabetic retinopathy as an early approach to delay its progression: the importance of oxidative stress and inflammation. Oxid. Med. Cell Longev. 2020, 1–23 (2020).
Gong, Y. et al. Erianin induces a jnk/sapk-dependent metabolic inhibition in human umbilical vein endothelial cells. In Vivo 18, 223–228 (2004).
Cushman, M. et al. Synthesis and evaluation of stilbene and dihydrostilbene derivatives as potential anticancer agents that inhibit tubulin polymerization. J. Med. Chem. 34, 2579–2588 (1991).
Ma, G. et al. Iuhibitory effects of Dendrobium chrysotoxutm and its constiuents on the mouse HePA and ESC. J. China Pharm. Univ. 25, 188–189 (1994).
Gong, Y. Q. et al. In vivo and in vitro evaluation of erianin, a novel anti-angiogenic agent. Eur. J. Cancer 40, 1554–1565 (2004).
Zhang, G. et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution. Sci. Rep. 6, 19029 (2016).
Hu, M. J. et al. Chromosome-scale assembly of the Kandelia obovata genome. Hortic. Res. 7, 75 (2020).
Kang, M. et al. A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine: an Isatis genome. Hortic. Res. 7, 18 (2020).
Chen, F. et al. The sequenced angiosperm genomes and genome databases. Front. Plant Sci. 9, 418 (2018).
Chen, F. et al. Genome sequences of horticultural plants: past, present, and future. Hortic. Res. 6, 112 (2019).
Yan, L. et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese orchid herb. Mol. Plant 8, 922–934 (2015).
Si, C. et al. DoRWA3 from Dendrobium officinale plays an essential role in acetylation of polysaccharides. Int. J. Mol. Sci. 21, 6250 (2020).
Zheng, S. G. et al. Genome-wide researches and applications on Dendrobium. Planta 248, 769–784 (2018).
Simão, F. A. et al. BUSCO online supplementary information: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212 (2015).
Cai, J. et al. The genome sequence of the orchid Phalaenopsis equestris. Nat. Genet. 47, 65–72 (2015).
Yuan, Y. et al. The Gastrodia elata genome provides insights into plant adaptation to heterotrophy. Nat. Commun. 9, 1615 (2018).
Sena, J. S. et al. Evolution of gene structure in the conifer Picea glauca: a comparative analysis of the impact of intron size. BMC Plant Biol. 14, 95 (2014).
De La Torre, A. R. et al. Insights into conifer giga-genomes. Plant Physiol. 166, 1724–1732 (2014).
Castillo-Davis, C. I. et al. Selection for short introns in highly expressed genes. Nat. Genet. 31, 415 (2002).
Keane, P. A. & Seoighe, C. Intron length coevolution across mammalian genomes. Mol. Biol. Evol. 33, 2682–2691 (2016).
Swinburne, I. A. & Silver, P. A. Intron delays and transcriptional timing during development. Dev. Cell 14, 324–330 (2008).
Artieri, C. G. & Fraser, H. B. Transcript length mediates developmental timing of gene expression across Drosophila. Mol. Biol. Evol. 31, 2879–2889 (2014).
Lin, I. W. et al. Nectar secretion requires sucrose phosphate synthases and the sugar transporter SWEET9. Nature 508, 546–549 (2014).
Chen, L. Q. et al. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science 335, 207–211 (2012).
Sun, M. X. et al. Arabidopsis RPG1 is important for primexine deposition and functions redundantly with RPG2 for plant fertility at the late reproductive stage. Plant Reprod. 26, 83–91 (2013).
Chen, L. Q. et al. A cascade of sequentially expressed sucrose transporters in the seed coat and endosperm provides nutrition for the Arabidopsis embryo. Plant Cell 27, 607–619 (2015).
Jiao, Y., Li, J., Tang, H. & Paterson, A. H. Integrated syntenic and phylogenomic analyses reveal an ancient genome duplication in monocots. Plant Cell 26, 2792–2802 (2014).
Tang, H. et al. jcvi: JCVI Utility Libraries (Zenodo, 2015). https://doi.org/10.5281/zenodo.31631.
Sang, X. et al. CHIMERIC FLORAL ORGANS1, encoding a monocot-specific MADS box protein, regulates floral organ identity in rice. Plant Physiol. 160, 788–807 (2012).
Tapia-López, R. et al. An AGAMOUS-related MADS-box gene, XAL1 (AGL12), regulates root meristem cell proliferation and flowering transition in Arabidopsis. Plant Physiol. 146, 1182–1192 (2008).
Par̆enicová, L. et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: new openings to the MADS world. Plant Cell 15, 1538–1551 (2003).
Masiero, S. et al. The emerging importance of type I MADS box transcription factors for plant reproduction. Plant Cell 23, 865–872 (2011).
Niu, S. C. et al. Morphological type identification of self-incompatibility in Dendrobium and its phylogenetic evolution pattern. Int. J. Mol. Sci. 19, 2595 (2018).
Wu, R. H. Studies on the Germplasm Resources of Dendrobium in China and their Genetic Relationships. Doctoral dissertation, Chin. Acad. For. (2007).
Watkins, J. L. & Pogson B. J. Prospects for carotenoid biofortification targeting retention and catabolism. Trends Plant Sci. 25, 501–512 (2020).
Johnson, P. R. & Ecker, J. R. The ethylene gas signal transduction pathway: a molecular perspective. Annu. Rev. Genet. 32, 227 (1998).
Yin, C. C. et al. Ethylene responses in rice roots and coleoptiles are differentially regulated by a carotenoid isomerase-mediated abscisic acid pathway. Plant Cell 27, 1061–1081 (2015).
Yu, Z. et al. Genome-wide identification and expression profile of TPS gene family in Dendrobium officinale and the role of DoTPS10 in linalool biosynthesis. Int. J. Mol. Sci. 21, 5419 (2020).
Jiang, S. Y., Jin, J. J., Sarojam, R. & Ramachandran, S. A comprehensive survey on the terpene synthase gene family provides new insight into its evolutionary patterns. Genome Biol. Evol. 11, 2078–2098 (2019).
Chen, F. The family of terpene synthases in plants: A mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 66, 212–229 (2011).
Tholl, D. Biosynthesis and biological functions of terpenoids in plants. Adv. Biochem. Eng. Biotechnol. 148, 63–106 (2015).
Bohlmann, J., Meyer-Gauen, G. & Croteau, R. Plant terpenoid synthases: molecular biology and phylogenetic analysis. Proc. Natl Acad. Sci. USA 95, 4126–4133 (1998).
Jiang, S. Y. et al. Comprehensive survey on the terpene synthase gene family provides new insight into its evolutionary patterns. Genome Biol. Evol. 11, 2078–2098 (2019).
Huang, X. L. Research on molecular regulation mechanism of formation of floral color and floral fragrance of Dendrobium chrysotoxum based on transcriptome sequencing. Doctoral dissertation, Chin. Acad. For. (2019).
Ranallo-Benavidez, T. R., Jaron, K. S. & Schatz, M. C. GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat. Commun. 11, 1–10 (2020).
Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res .27, 722–736 (2017).
Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).
Chen, Y. et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. GigaScience 7, 120 (2017).
Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 3, 95–98 (2016).
Dudchenko, O. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92–95 (2017).
Jurka, J. et al. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 110, 462–467 (2005).
Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, i351–i358 (2005).
Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265–W268 (2007).
Benson, G. et al. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573–580 (1999).
Slater, G. S. & Birney, E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 6, 31 (2005).
Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, W435–W439 (2006).
Johnson, A. D. et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 24, 2938–2939 (2008).
Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295 (2015).
Kim, D. et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915 (2019).
Holt, C. & Yandell, M. MAKER2: an annotation pipeline and genome database management tool for second generation genome projects. BMC Bioinform. 12, 491 (2011).
Kent, W. J. BLAT—the BLAST-like alignment tool. Genome Res. 12, 656–664 (2002).
Boeckmann, B. et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 31, 365–370 (2003).
Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240 (2014).
Ashburner, M. et al. Gene Ontology: tool for the unification of biology. Nat. Genet. 25, 25–29 (2000).
Lowe, T. M. & Eddy, S. R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964 (1997).
Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933–2935 (2013).
Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49 (2012).
Yang, Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Bioinformatics 13, 555–556 (1997).
De Bie, T. et al. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269–1271 (2006).
El-Gebali, S. et al. The Pfam protein families database in 2019. Nucleic Acids Res. 47, D427–D432 (2019).
Eddy, S. R. Accelerated profile HMM searches. PLoS Comput. Biol. 7, e1002195 (2011).
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).
Guindon, S., Lethiec, F., Duroux, P. & Gascuel, O. PHYML Online—a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 33, W557–W559 (2005).
Wu, S. S. et al. The genome sequence of star fruit (Averrhoa carambola). Hortic. Res. 7, 95 (2020).
Chen, S. P. et al. The Phoebe genome sheds light on the evolution of magnoliids. Hortic. Res. 7, 146 (2020).
Letunic, I., Doerks, T. & Bork, P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, D257–D260 (2015).
Kumar, S. et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).
Finkelstein, R. Abscisic Acid synthesis and response. Arabidopsis Book 11, e0166–e0166 (2013).
Lin, Z., Zhong, S. & Grierson, D. Recent advances in ethylene research. J. Exp. Bot. 60, 3311–3336 (2009).
Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).
Jones, D. T., Taylor, W. R. & Thornton, J. M. The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282 (1992).
Patro, R. et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419 (2017).
Nisar, N. et al. Carotenoid metabolism in plants. Mol. Plant 8, 68–82 (2015).
Sun, T. et al. Carotenoid metabolism in plants: the role of plastids. Mol. Plant 11, 58–74 (2017).
Acknowledgements
This project was supported by the Guangdong Innovation Research Team Fund (2014ZT05S078); National Natural Science Foundation of China (grants 31571252 and 31772322); Guangdong Special Support Program for Young Talents in Innovation Research of Science and Technology (2019TQ05N940); Shenzhen Peacock Grant (827/000189); Science and Technology Program of Guangdong Province, China (2019B121202006); Program of Forestry Administration of Guangdong Province (E036011002); Department for Wildlife and Forest Plant Protection of the National Forest and Grassland Administration (2019073010); National Key Research and Development Program of China (No. 2018YFD1000400); Special Research Foundation of Hebei Agricultural University (YJ201848); and Natural Science Foundation of Hebei Province (C2019204295).
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, Y., Zhang, GQ., Zhang, D. et al. Chromosome-scale assembly of the Dendrobium chrysotoxum genome enhances the understanding of orchid evolution. Hortic Res 8, 183 (2021). https://doi.org/10.1038/s41438-021-00621-z
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41438-021-00621-z
This article is cited by
-
The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii
BMC Plant Biology (2024)
-
Integrated metabolome and transcriptome analysis provides insight into the leaf color change of Cymbidium ensifolium
Acta Physiologiae Plantarum (2024)
-
Genome-wide identification, characterization and transcriptional profile of the SWEET gene family in Dendrobium officinale
BMC Genomics (2023)
-
Comprehensive chemical profiling of two Dendrobium species and identification of anti-hepatoma active constituents from Dendrobium chrysotoxum by network pharmacology
BMC Complementary Medicine and Therapies (2023)
-
Approaches to increase the validity of gene family identification using manual homology search tools
Genetica (2023)
