Search

Abstract Background Castanopsis carlesii is a dominant tree species in subtropical evergreen broad-leaved forests and holds significant ecological value. It serves as an excellent timber tree species and raw material for cultivating edible fungi. Henry Chinquapin (Castanea henryi) wood is known for its hardness and resistance to water and moisture, making it an exceptional timber species. Additionally, its fruit has a sweet and fruity taste, making it a valuable food source. However, the mitogenomes of these species have not been previously reported. To gain a better understanding of them, this study successfully assembled high-quality mitogenomes of C. carlesii and Ca. henryi for the first time. Results Our research reveals that the mitochondrial DNA (mtDNA) of C. carlesii exhibits a unique multi-branched conformation, while Ca. henryi primarily exists in the form of two independent molecules that can be further divided into three independent molecules through one pair of long repetitive sequences. The size of the mitogenomes of C. carlesii and Ca. henryi are 592,702 bp and 379,929 bp respectively, which are currently the largest and smallest Fagaceae mitogenomes recorded thus far. The primary factor influencing mitogenome size is dispersed repeats. Comparison with published mitogenomes from closely related species highlights differences in size, gene loss patterns, codon usage preferences, repetitive sequences, as well as mitochondrial plastid DNA segments (MTPTs). Conclusions Our study enhances the understanding of mitogenome structure and evolution in Fagaceae, laying a crucial foundation for future research on cell respiration, disease resistance, and other traits in this family.

Abstract Simple sequence repeats (SSRs) are found in nonrandom distributions in genomes and are thought to impact gene expression. The distribution patterns of 48 295 SSRs of Paphiopedilum malipoense are mined and characterized based on the first full‐length transcriptome and comprehensive transcriptome dataset from 12 organs. Statistical genomics analyses are used to investigate how SSRs in transcripts affect gene expression. The results demonstrate the correlations between SSR distributions, characteristics, and expression level. Nine expression‐modulating motifs (expMotifs) are identified and a model is proposed to explain the effect of their key features, potency, and gene function on an intra‐transcribed region scale. The expMotif‐transcribed region combination is the most predominant contributor to the expression‐modulating effect of SSRs, and some intra‐transcribed regions are critical for this effect. Genes containing the same type of expMotif‐SSR elements in the same transcribed region are likely linked in function, regulation, or evolution aspects. This study offers novel evidence to understand how SSRs regulate gene expression and provides potential regulatory elements for plant genetic engineering.

Abstract Background Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation. Results We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnN GUU –rpl32 > rpoB–trnC GCA > psbK–psbI > psaC–rps15 > trnE UUC –trnT GGU > accD–psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47–56 SSRs and 11–14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support. Conclusion This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

The Homeodomain-Leucine Zipper (HD-ZIP) transcription factors play a pivotal role in governing various aspects of plant growth, development, and responses to abiotic stress. Despite the well-established importance of HD-ZIPs in many plants, their functions in Acoraceae, the basal lineage of monocots, remain largely unexplored. Using recently published whole-genome data, we identified 137 putative HD-ZIPs in two Acoraceae species, Acorus gramineus and Acorus calamus. These HD-ZIP genes were further classified into four subfamilies (I, II, III, IV) based on phylogenetic and conserved motif analyses, showcasing notable variations in exon-intron patterns among different subfamilies. Two microRNAs, miR165/166, were found to specifically target HD-ZIP III genes with highly conserved binding sites. Most cis-acting elements identified in the promoter regions of Acoraceae HD-ZIPs are involved in modulating light and phytohormone responsiveness. Furthermore, our study revealed an independent duplication event in Ac. calamus and a one-to-multiple correspondence between HD-ZIP genes of Ac. calamus and Ac. gramineus. Expression profiles obtained from qRT-PCR demonstrated that HD-ZIP I genes are strongly induced by salinity stress, while HD-ZIP II members have contrasting stress responses in two species. HD-ZIP III and IV genes show greater sensitivity in stress-bearing roots. Taken together, these findings contribute valuable insights into the roles of HD-ZIP genes in stress adaptation and plant resilience in basal monocots, illuminating their multifaceted roles in plant growth, development, and response to abiotic stress.

The TIFY family is a plant-specific gene family that is involved in regulating a variety of plant processes, including developmental and defense responses. The Cymbidium species have certain ornamental and ecological value. However, the characteristics and functions of TIFY genes in Cymbidium remain poorly understood. This study conducted a genome analysis of the TIFY gene family in Cymbidium goeringii, C. ensifolium, and C. sinense and investigated their physicochemical properties, phylogenetic relationships, gene structures, and expression patterns under heat stress in C. goeringii. C. goeringii (26), C. ensifolium (19), and C. sinense (21). A total of 66 TIFY genes were identified, and they were classified into four subfamilies (JAZ, ZML, PPD, and TIFY) based on their systematic evolutionary relationships. Sequence analysis showed that TIFYs contained a conserved TIFY domain and that genes within the same subfamily had structural similarity. Analysis of cis-regulatory elements revealed that these genes contain numerous light-responsive elements and stress-responsive elements. We subjected C. goeringii (16 h light/8 h dark) to 24 h of 38 °C high-temperature stress in a climate chamber. Additionally, results from RT-qPCR experiments showed that under heat stress, the expression levels of eight genes in C. goeringii show significant differences. Among them, the JAZ subfamily exhibited the strongest response to heat stress, initially showing upregulation followed by a downregulation trend. In conclusion, this study investigated the role of TIFY genes in three Cymbidium species, providing insights into C. goeringii under heat stress.

Apostasia fujianica belongs to the genus Apostasia and is part of the basal lineage in the phylogenetic tree of the Orchidaceae. Currently, there are only ten reported complete mitochondrial genomes in orchids, which greatly hinders the understanding of mitochondrial evolution in Orchidaceae. Therefore, we assembled and annotated the mitochondrial genome of A. fujianica, which has a length of 573,612 bp and a GC content of 44.5%. We annotated a total of 44 genes, including 30 protein-coding genes, 12 tRNA genes, and two rRNA genes. We also performed relative synonymous codon usage (RSCU) analysis, repeat sequence analysis, intergenomic transfer (IGT) analysis, and Ka/Ks analysis for A. fujianica and conducted RNA editing site analysis on the mitochondrial genomes of eight orchid species. We found that most protein-coding genes are under purifying selection, but nad6 is under positive selection, with a Ka/Ks value of 1.35. During the IGT event in A. fujianica’s mitogenome, the trnN-GUU, trnD-GUC, trnW-CCA, trnP-UGG, and psaJ genes were identified as having transferred from the plastid to the mitochondrion. Compared to other monocots, the family Orchidaceae appears to have lost the rpl10, rpl14, sdh3, and sdh4 genes. Additionally, to further elucidate the evolutionary relationships among monocots, we constructed a phylogenetic tree based on the complete mitogenomes of monocots. Our study results provide valuable data on the mitogenome of A. fujianica and lay the groundwork for future research on genetic variation, evolutionary relationships, and breeding of Orchidaceae.

As an important genus in Orchidaceae, Cymbidium has rich ecological diversity and significant economic value. DNA binding with one zinc finger (Dof) proteins are pivotal plant-specific transcription factors that play crucial roles in the growth, development, and stress response of plants. Although the Dof genes have been identified and functionally analyzed in numerous plants, exploration in Orchidaceae remains limited. We conducted a thorough analysis of the Dof gene family in Cymbidium goeringii, C. ensifolium, and C. sinensis. In total, 91 Dof genes (27 CgDofs, 34 CeDofs, 30 CsDofs) were identified, and Dof genes were divided into five groups (I–V) based on phylogenetic analysis. All Dof proteins have motif 1 and motif 2 conserved domains and over half of the genes contained introns. Chromosomal localization and collinearity analysis of Dof genes revealed their evolutionary relationships and potential gene duplication events. Analysis of cis-elements in CgDofs, CeDofs, and CsDofs promoters showed that light-responsive cis-elements were the most common, followed by hormone-responsive elements, plant growth-related elements, and abiotic stress response elements. Dof proteins in three Cymbidium species primarily exhibit a random coil structure, while homology modeling exhibited significant similarity. In addition, RT-qPCR analysis showed that the expression levels of nine CgDofs changed greatly under heat stress. CgDof03, CgDof22, CgDof27, CgDof08, and CgDof23 showed varying degrees of upregulation. Most upregulated genes under heat stress belong to group I, indicating that the Dof genes in group I have great potential for high-temperature resistance. In conclusion, our study systematically demonstrated the molecular characteristics of Dof genes in different Cymbidium species, preliminarily revealed the patterns of heat stress, and provided a reference for further exploration of stress breeding in orchids.

ZF-HD transcription factors, which are unique to land plants, are involved in the regulation of abiotic stress response and related signaling pathways, and play a crucial role in plant growth and development. Dendrobium is one of the largest genera of orchids, with a high ornamental and ecological value. However, the specific functions of the ZF-HDs in Dendrobium remain unknown. In this study, we identified a total of 53 ZF-HDs from D. chrysotoxum (17), D. catenatum (23), and D. huoshanense (13), and analyzed their physicochemical properties, phylogenetic relationships, chromosomal locations, protein structures, conserved motifs, and expression patterns. The phylogenetic relationships revealed that 53 ZF-HDs were classified into six subfamilies (ZHDI–V and MIF), and all ZF-HD proteins contained motif 1 and motif 4 conserved domains, while a minority of these proteins had exons. The analysis of cis-elements in the promoters of ZF-HDs from three Dendrobium species showed that growth- and development-related elements were the most prevalent, followed by hormone response and abiotic stress response elements. Through collinearity analysis, 14 DchZF-HDs were found to be collinear with DhuZF-HDs, and 12 DchZF-HDs were found to be collinear with DcaZF-HDs. Furthermore, RT-qPCR analysis revealed that DchZF-HDs play a regulatory role in the development of lateral organs during the flowering process. The results indicated that DchZHD2 plays a role in the unpigmented bud stage, while DchMIF8 and DchZHD16 play significant roles during the pigmented bud and initial bloom stages. Hence, this study provides a crucial basis for further exploring ZF-HDs functions in regulating the floral organs of orchids.

The GRAS gene family, responsible for encoding transcription factors, serves pivotal functions in plant development, growth, and responses to stress. The exploration of the GRAS gene family within the Orchidaceae has been comparatively limited, despite its identification and functional description in various plant species. This study aimed to conduct a thorough examination of the GRAS gene family in Cymbidum goeringii, focusing on its physicochemical attributes, phylogenetic associations, gene structure, cis-acting elements, and expression profiles under heat stress. The results show that a total of 54 CgGRASs were pinpointed from the genome repository and categorized into ten subfamilies via phylogenetic associations. Assessment of gene sequence and structure disclosed the prevalent existence of the VHIID domain in most CgGRASs, with around 57.41% (31/54) CgGRASs lacking introns. The Ka/Ks ratios of all CgGRASs were below one, indicating purifying selection across all CgGRASs. Examination of cis-acting elements unveiled the presence of numerous elements linked to light response, plant hormone signaling, and stress responsiveness. Furthermore, CgGRAS5 contained the highest quantity of cis-acting elements linked to stress response. Experimental results from RT-qPCR demonstrated notable variations in the expression levels of eight CgGRASs after heat stress conditions, particularly within the LAS, HAM, and SCL4/7 subfamilies. In conclusion, this study revealed the expression pattern of CgGRASs under heat stress, providing reference for further exploration into the roles of CgGRAS transcription factors in stress adaptation.

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

Wolfberry, also known as goji berry or Lycium barbarum, is a highly valued fruit with significant health benefits and nutritional value. For more efficient and comprehensive usage of published L. barbarum genomic data, we established the Wolfberry database. The utility of the Wolfberry Genome Database (WGDB) is highlighted through the Genome browser, which enables the user to explore the L. barbarum genome, browse specific chromosomes, and access gene sequences. Gene annotation features provide comprehensive information about gene functions, locations, expression profiles, pathway involvement, protein domains, and regulatory transcription factors. The transcriptome feature allows the user to explore gene expression patterns using transcripts per kilobase million (TPM) and fragments per kilobase per million mapped reads (FPKM) metrics. The Metabolism pathway page provides insights into metabolic pathways and the involvement of the selected genes. In addition to the database content, we also introduce six analysis tools developed for the WGDB. These tools offer functionalities for gene function prediction, nucleotide and amino acid BLAST analysis, protein domain analysis, GO annotation, and gene expression pattern analysis. The WGDB is freely accessible at https://cosbi7.ee.ncku.edu.tw/Wolfberry/. Overall, WGDB serves as a valuable resource for researchers interested in the genomics and transcriptomics of L. barbarum. Its user-friendly web interface and comprehensive data facilitate the exploration of gene functions, regulatory mechanisms, and metabolic pathways, ultimately contributing to a deeper understanding of wolfberry and its potential applications in agronomy and nutrition.

The ARF gene family is a representative transcription factor that plays a crucial role in the regulation of various growth and development processes in plants. Although the ARF gene family has been identified in five Orchidaceae species, limited research has been conducted on the ARF gene within Dendrobium. To explore ARF family genes in different Dendrobium species, we selected chromosome-level genomic data from D. nobile, D. chrysotoxum, and D. huoshanense for genome-wide identification, and to analyze expression patterns in the D. nobile flower. In this study, 13, 18, and 23 ARF genes were identified in the genomes of D. chrysotoxum, D. huoshanense, and D. nobile, respectively. These genes were then subsequently classified into four classes (Classes I, II, III, and IV) based on our phylogenetic analysis. Additional protein sequence analysis found that 30 ARF proteins with three classically conserved structural domains (BDB, MR, and RD) were present in the three Dendrobium species. Our gene structure comparative analysis also found the same evolutionary branch with similar intron-exon structural features. Specifically, Class I and Class III display longer introns that potentially constitute a distinctive characteristic of Dendrobium. Gene expression patterns analysis showed the potential involvement of DnoARF 5, 7, 10, 11, 12, 13, 16, 17, and 23 in initial differentiation and flower bud formation. Likewise, transcriptomic analysis and RT-qPCR expression profiles indicated flower-specific expression patterns for four ARF genes (DnoARF5, DnoARF6, DnoARF16, and DnoARF22), which suggest an important impact on flower development and regulation. Ultimately, this study provides comprehensive data to elucidate the potential functions of ARF genes in three Dendrobium species and suggests new insights for further exploration concerning the function and regulatory mechanisms in Dendrobium during flower development.

The TIFY gene family (formerly known as the zinc finger proteins expressed in inflorescence meristem (ZIM) family) not only functions in plant defense responses but also are widely involved in regulating plant growth and development. However, the identification and functional analysis of TIFY proteins remain unexplored in Orchidaceae. Here, we identified 19 putative TIFY genes in the Phalaenopsis aphrodite genome. The phylogenetic tree classified them into four subfamilies: 14 members from JAZ, 3 members from ZML, and 1 each from PPD and TIFY. Sequence analysis revealed that all Phalaenopsis TIFY proteins contained a TIFY domain. Exon–intron analysis showed that the intron number and length of Phalaenopsis TIFY genes varied, whereas the same subfamily and subgroup genes had similar exon or intron numbers and distributions. The most abundant cis-elements in the promoter regions of the 19 TIFY genes were associated with light responsiveness, followed by MeJA and ABA, indicating their potential regulation by light and phytohormones. The 13 candidate TIFY genes screened from the transcriptome data exhibited two types of expression trends, suggesting their different roles in cell proliferation and cell expansion of floral organ growth during Phalaenopsis flower opening. Overall, this study serves as a background for investigating the underlying roles of TIFY genes in floral organ growth in Phalaenopsis.

Apostasioideae, the early divergent subfamily of Orchidaceae, comprises Apostasia and Neuwiedia genera with approximately 20 species. Despite extensive research on Apostasioideae, previous studies have struggled to resolve taxonomic issues, particularly concerning the position of species within this subfamily. Here, we sequenced and annotated plastomes of Apostasia fujianica and Neuwiedia malipoensis, unveiling their phylogenetic relationships and shared plastome features with the other five published plastomes. We identified and analyzed the length, GC content, repeat sequences, and RSCU values of the chloroplast genomes. It is noteworthy that the chloroplast genome of N. malipoensis stands out as the largest among all known chloroplast genomes within the Apostasioideae subfamily, primarily due to contributions from both the LSC and SSC regions. Furthermore, our analysis revealed three unique structural rearrangements located approximately 10k–47k bp (ycf3–trnS-GCU) and 58k–59k bp(accD) in the LSC region and 118k–119k (ndhI) bp in the SSC region of the chloroplast genomes across all five species within the Apostasia genus, which presents a potential avenue for identifying distinctive chloroplast genetic markers, setting them apart from other orchid plants. And a total of four mutational hotspots (rpoC2, atpH, rps4, ndhK, and clpP) were identified. Moreover, our study suggested that Apostasia and Neuwiedia formed a monophyletic group, with Apostasia being sister to Neuwiedia. Within the Apostasia genus, five species were classified into two major clades, represented as follows: (A. odorata (A. shenzhenica and A. fujianica) (A. ramifera and A. wallichii)). These findings hold significance in developing DNA barcoding of Apostasioideae and contribute to the further phylogenetic understanding of Apostasioideae species.

Cymbidium sinense, a type of orchid plant, is more drought-resistant and ornamental than other terrestrial orchids. Research has shown that many members of the NUCLEAR FACTOR Y (NF-Y) transcription factor family are responsive to plant growth, development, and abiotic stress. However, the mechanism of the NF-Y gene family’s response to abiotic stress in orchids has not yet been reported. In this study, phylogenetic analysis allowed for 27 CsNF-Y genes to be identified (5 CsNF-YAs, 9 CsNF-YBs, and 13 CsNF-YC subunits), and the CsNF-Ys were homologous to those in Arabidopsis and Oryza. Protein structure analysis revealed that different subfamilies contained different motifs, but all of them contained Motif 2. Secondary and tertiary protein structure analysis indicated that the CsNF-YB and CsNF-YC subfamilies had a high content of alpha helix structures. Cis-element analysis showed that elements related to drought stress were mainly concentrated in the CsNF-YB and CsNF-YC subfamilies, with CsNF-YB3 and CsNF-YC12 having the highest content. The results of a transcriptome analysis showed that there was a trend of downregulation of almost all CsNF-Ys in leaves under drought stress, while in roots, most members of the CsNF-YB subfamily showed a trend of upregulation. Additionally, seven genes were selected for real-time reverse transcription quantitative PCR (qRT-PCR) experiments. The results were generally consistent with those of the transcriptome analysis. The regulatory roles of CsNF-YB 1, 2, and 4 were particularly evident in the roots. The findings of our study may make a great contribution to the understanding of the role of CsNF-Ys in stress-related metabolic processes.

The purpose of this study was to examine the antecedents of academic emotions among university students through real-time and multilevel analyses. We explored the link between state self-control and academic emotions, the influence of self-determination on the relationship between them, and the moderating role of trait self-control. The data was collected over seven consecutive days during which university students (N = 155) completed smartphone questionnaires. Data was organized into hierarchical two-level structures, where situations (Level 1) nested within individuals (Level 2). The results showed that positive emotions were negatively predicted by state self-control while negative emotions were positively predicted by state self-control. Moreover, state self-control under low self-determination was a significant negative predictor of positive emotions, whereas high self-determination had no predictive effect on positive emotions. The relationship between self-determination and negative emotions was further moderated by trait self-control. The limitations of this study and future research directions are also discussed.

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.

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

Abstract Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome‐level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex‐determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

The Tripterospermum, comprising 34 species, is a genus of Gentianaceae. Members of Tripterospermum are mostly perennial, entwined herbs with high medicinal value and rich in iridoids, xanthones, flavonoids, and triterpenes. However, our inadequate understanding of the differences in the plastid genome sequences of Tripterospermum species has severely hindered the study of their evolution and phylogeny. Therefore, we first analyzed the 86 Gentianae plastid genomes to explore the phylogenetic relationships within the Gentianae subfamily where Tripterospermum is located. Then, we analyzed six plastid genomes of Tripterospermum, including two newly sequenced plastid genomes and four previously published plastid genomes, to explore the plastid genomes’ evolution and phylogenetic relationships in the genus Tripterospermum. The Tripterospermum plastomes have a quadripartite structure and are between 150,929 and 151,350 bp in size. The plastomes of Tripterospermum encoding 134 genes were detected, including 86 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and three pseudogenes (infA, rps19, and ycf1). The result of the comparison shows that the Tripterospermum plastomes are very conserved, with the total plastome GC content ranging from 37.70% to 37.79%. In repeat sequence analysis, the number of single nucleotide repeats (A/T) varies among the six Tripterospermum species, and the identified main long repeat types are forward and palindromic repeats. The degree of conservation is higher at the SC/IR boundary. The regions with the highest divergence in the CDS and the intergenic region (IGS) are psaI and rrn4.5-rrn5, respectively. The average pi of the CDS and the IGS are only 0.071% and 0.232%, respectively, indicating that the Tripterospermum plastomes are highly conserved. Phylogenetic analysis indicated that Gentianinae is divided into two clades, with Tripterospermum as a sister to Sinogeniana. Phylogenetic trees based on CDS and CDS + IGS combined matrices have strong support in Tripterospermum. These findings contribute to the elucidation of the plastid genome evolution of Tripterospermum and provide a foundation for further exploration and resource utilization within this genus.

Bulbophyllum is one of the largest genera and presents some of the most intricate taxonomic problems in the family Orchidaceae, including species of ornamental and medical importance. The lack of knowledge regarding the characterization of Bulbophyllum chloroplast (cp) genomes has imposed current limitations on our study. Here, we report the complete cp genomes of seven Bulbophyllum species, including B. ambrosia, B. crassipes, B. farreri, B. hamatum, B. shanicum, B. triste, and B. violaceolabellum, and compared with related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny. A total of 28 Bulbophyllum cp genomes exhibit typical quadripartite structures with lengths ranging from 145,092 bp to 165,812 bp and a GC content of 36.60% to 38.04%. Each genome contained 125–132 genes, encompassing 74–86 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The genome arrangements, gene contents, and length were similar, with differences observed in ndh gene composition. It is worth noting that there were exogenous fragment insertions in the IR regions of B. crassipes. A total of 18–49 long repeats and 38–80 simple sequence repeats (SSRs) were detected and the single nucleotide (A/T) was dominant in Bulbophyllum cp genomes, with an obvious A/T preference. An analysis of relative synonymous codon usage (RSCU) revealed that leucine (Leu) was the most frequently used codon, while cysteine (Cys) was the least used. Six highly variable regions (rpl32-trnLUAG > trnTUGU-trnLUAA > trnFGAA-ndhJ > rps15-ycf1 > rbcL-accD > psbI-trnSGCU) and five coding sequences (ycf1 > rps12 > matK > psbK > rps15) were identified as potential DNA markers based on nucleotide diversity. Additionally, 31,641 molecular diagnostic characters (MDCs) were identified in complete cp genomes. A phylogenetic analysis based on the complete cp genome sequences and 68 protein-coding genes strongly supported that 28 Bulbophyllum species can be divided into four branches, sects. Brachyantha, Cirrhopetalum, and Leopardinae, defined by morphology, were non-monophyly. Our results enriched the genetic resources of Bulbophyllum, providing valuable information to illustrate the complicated taxonomy, phylogeny, and evolution process of the genus.

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.

Epipogium roseum, commonly known as one of the ghost orchids due to its rarity and almost transparent color, is a non-photosynthetic and fully mycoheterotrophic plant. Given its special nutritional strategies and evolutionary significance, the mitogenome was first characterized, and three plastomes sampled from Asia were assembled. The plastomes were found to be the smallest among Orchidaceae, with lengths ranging from 18,339 to 19,047 bp, and exhibited high sequence variety. For the mitogenome, a total of 414,552 bp in length, comprising 26 circular chromosomes, were identified. A total of 54 genes, including 38 protein-coding genes, 13 tRNA genes, and 3 rRNA genes, were annotated. Multiple repeat sequences spanning a length of 203,423 bp (45.47%) were discovered. Intriguingly, six plastid regions via intracellular gene transfer and four plastid regions via horizontal gene transfer to the mitogenome were observed. The phylogenomics, incorporating 90 plastomes and 56 mitogenomes, consistently revealed the sister relationship of Epipogium and Gastrodia, with a bootstrap percentage of 100%. These findings shed light on the organelle evolution of Orchidaceae and non-photosynthetic plants.

Auxin Response Factors (ARFs) mediate auxin signaling and govern diverse biological processes. However, a comprehensive analysis of the ARF gene family and identification of their key regulatory functions have not been conducted in Melastoma dodecandrum, leading to a weak understanding of further use and development for this functional shrub. In this study, we successfully identified a total of 27 members of the ARF gene family in M. dodecandrum and classified them into Class I–III. Class II–III showed more significant gene duplication than Class I, especially for MedARF16s. According to the prediction of cis-regulatory elements, the AP2/ERF, BHLH, and bZIP transcription factor families may serve as regulatory factors controlling the transcriptional pre-initiation expression of MedARF. Analysis of miRNA editing sites reveals that miR160 may play a regulatory role in the post-transcriptional expression of MeARF. Expression profiles revealed that more than half of the MedARFs exhibited high expression levels in the stem compared to other organs. While there are some specific genes expressed only in flowers, it is noteworthy that MedARF16s, MedARF7A, and MedARF9B, which are highly expressed in stems, also demonstrate high expressions in other organs of M. dodecandrum. Further hormone treatment experiments revealed that these MedARFs were sensitive to auxin changes, with MedARF6C and MedARF7A showing significant and rapid changes in expression upon increasing exogenous auxin. In brief, our findings suggest a crucial role in regulating plant growth and development in M. dodecandrum by responding to changes in auxin. These results can provide a theoretical basis for future molecular breeding in Myrtaceae.

Heat shock factors (HSFs) are the key regulators of heat stress responses and play pivotal roles in tissue development and the temperature-induced regulation of secondary metabolites. In order to elucidate the roles of HSFs in Cymbidium ensifolium, we conducted a genome-wide identification of CeHSF genes and predicted their functions based on their structural features and splicing patterns. Our results revealed 22 HSF family members, with each gene containing more than one intron. According to phylogenetic analysis, 59.1% of HSFs were grouped into the A subfamily, while subfamily HSFC contained only two HSFs. And the HSF gene families were differentiated evolutionarily between plant species. Two tandem repeats were found on Chr02, and two segmental duplication pairs were observed on Chr12, Chr17, and Chr19; this provided evidence for whole-genome duplication (WGD) events in C. ensifolium. The core region of the promoter in most CeHSF genes contained cis-acting elements such as AP2/ERF and bHLH, which were associated with plant growth, development, and stress responses. Except for CeHSF11, 14, and 19, each of the remaining CeHSFs contained at least one miRNA binding site. This included binding sites for miR156, miR393, and miR319, which were responsive to temperature and other stresses. The HSF gene family exhibited significant tissue specificity in both vegetative and floral organs of C. ensifolium. CeHSF13 and CeHSF15 showed relatively significant expression in flowers compared to other genes. During flower development, CeHSF15 exhibited markedly elevated expression in the early stages of flower opening, implicating critical regulatory functions in organ development and floral scent-related regulations. During the poikilothermic treatment, CeHSF14 was upregulated over 200-fold after 6 h of heat treatment. CeHSF13 and CeHSF14 showed the highest expression at 6 h of low temperature, while the expression of CeHSF15 and CeHSF21 continuously decreased at a low temperature. The expression patterns of CeHSFs further confirmed their role in responding to temperature stress. Our study may help reveal the important roles of HSFs in plant development and metabolic regulation and show insight for the further molecular design breeding of C. ensifolium.

Though conserved in higher plants, the WOX transcription factors play crucial roles in plant growth and development of Melastoma dodecandrum Lour., which shows pioneer position in land ecosystem formation and produces nutritional fruits. Identifying the WOX family genes in M. dodecandrum is imperative for elucidating its growth and development mechanisms. However, the WOX genes in M. dodecandrum have not yet been characterized. In this study, by identification 22 WOX genes in M. dodecandrum based on current genome data, we classified family genes into three clades and nine types with homeodomains. We highlighted gene duplications of MedWOX4, which offered evidences of whole-genome duplication events. Promoter analysis illustrated that cis-regulatory elements related to light and stress responses and plant growth were enriched. Expression pattern and RT-qPCR results demonstrated that the majority of WOX genes exhibited expression in the stem. MedWOX13s displayed highest expression across various tissues. MedWOX4s displayed a specific expression in the stem. Collectively, our study provided foundations for elucidating WOX gene functions and further molecular design breeding in M. dodecandrum.

The PEBP gene family plays a significant role in regulating flower development and formation. To understand its function in Dendrobium chrysotoxum and D. nobile flowering, we identified 22 PEBP genes (11 DchPEBPs and 11 DnoPEBPs) from both species. We conducted analyses on their conserved domains and motifs, phylogenetic relationships, chromosome distribution, collinear correlation, and cis elements. The classification results showed that the 22 PEBPs were mainly divided into three clades, as follows: FT, MFT, and TFL1. A sequence analysis showed that most PEBP proteins contained five conserved domains, while a gene structure analysis revealed that 77% of the total PEBP genes contained four exons and three introns. The promoter regions of the 22 PEBPs contained several cis elements related to hormone induction and light response. This suggests these PEBPs could play a role in regulating flower development by controlling photoperiod and hormone levels. Additionally, a collinearity analysis revealed three pairs of duplicate genes in the genomes of both D. chrysotoxum and D. nobile. Furthermore, RT-qPCR has found to influence the regulatory effect of DchPEBPs on the development of flower organs (sepals, petals, lip, ovary, and gynostemium) during the flowering process (bud, transparent stage, and initial bloom). The results obtained imply that DchPEBP8 and DchPEBP9 play a role in the initial bloom and that DchPEBP7 may inhibit flowering processes. Moreover, DchPEBP9 may potentially be involved in the development of reproductive functionality. PEBPs have regulatory functions that modulate flowering. FT initiates plant flowering by mediating photoperiod and temperature signals, while TFL1 inhibits flowering processes. These findings provide clues for future studies on flower development in Dendrobium.

WNK (With No Lysine) kinases are members of serine/threonine protein kinase family, which lack conserved a catalytic lysine (K) residue in protein kinase subdomain II and this residue is replaced by either asparagine, serine, or glycine residues. They are involved in various physiological regulations of flowering time, circadian rhythms, and abiotic stresses in plants. In this study, we identified the WNK gene family in two species of Acorus, and analyzed their phylogenetic relationship, physiochemical properties, subcellular localization, collinearity, and cis-elements. The results showed twenty-two WNKs in two Acorus (seven in Ac. gramineus and fifteen in Ac. calamus) have been identified and clustered into five main clades phylogenetically. Gene structure analysis showed all WNKs possessed essential STKc_WNK or PKc_like superfamily domains, and the gene structures and conserved motifs of the same clade were similar. All the WNKs harbored a large number of light response elements, plant hormone signaling elements, and stress resistance elements. Through a collinearity analysis, two and fourteen segmental duplicated gene pairs were identified in the Ac. gramineus and Ac. calamus, respectively. Moreover, we observed tissue-specificity of WNKs in Acorus using transcriptomic data, and their expressions in response to salt stress and cold stress were analyzed by qRT-PCR. The results showed WNKs are involved in the regulation of abiotic stresses. There were significant differences in the expression levels of most of the WNKs in the leaves and roots of Acorus under salt stress and cold stress, among which two members in Ac. gramineus (AgWNK3 and AgWNK4) and two members in Ac. calamus (AcWNK8 and AcWNK12) were most sensitive to stress. In summary, this paper will significantly contribute to the understanding of WNKs in monocots and thus provide a set up for functional genomics studies of WNK protein kinases.

Angraecum, commonly known as Darwin’s orchid, is the largest genus of Angraecinae (Orchidaceae). This genus exhibits a high morphological diversity, making it as a good candidate for macroevolutionary studies. In this study, four complete plastomes of Angraecum were firstly reported and the potential variability hotspots were explored. The plastomes possessed the typical quadripartite structure and ranged from 150,743 to 151,818 base pair (bp), with a guanine–cytosine (GC) content of 36.6–36.9%. The plastomes all contained 120 genes, consisting of 74 protein-coding genes (CDS), 38 transfer RNA (tRNA) genes and 8 ribosomal RNA (rRNA) genes; all ndh genes were pseudogenized or lost. A total of 30 to 46 long repeats and 55 to 63 SSRs were identified. Relative synonymous codon usage (RSCU) analysis indicated a high degree of conservation in codon usage bias. The Ka/Ks ratios of most genes were lower than 1, indicating that they have undergone purifying selection. Based on the ranking of Pi (nucleotide diversity) values, five regions (trnSGCU-trnGGCC, ycf1-trnNGGU, trnNGUU-rpl32, psaC-ndhE and trnSGCU-trnGGCC) and five protein-coding genes (rpl32, rps16, psbK, rps8, and ycf1) were identified. The consistent and robust phylogenetic relationships of Angraecum were established based on a total of 40 plastomes from the Epidendroideae subfamily. The genus Angraecum was strongly supported as a monophyletic group and sister to Aeridinae. Our study provides an ideal system for investigating molecular identification, plastome evolution and DNA barcoding for Angraecum.

Jasmine, a recently domesticated shrub, is renowned for its use as a key ingredient in floral tea and its captivating fragrance, showcasing significant ornamental and economic value. When cultivated to subtropical zone, a significant abiotic stress adaptability occurs among different jasmine varieties, leading to huge flower production changes and plantlet survival. The bZIP transcription factors (TFs) are reported to play indispensable roles in abiotic stress tolerance. Here, we performed a genome-level comparison of bZIPs using three-type jasmine genomes. Based on their physicochemical properties, conserved motif analysis and phylogenetic analysis, about 63 bZIP genes were identified and clustered in jasmine genomes, noting a difference of one member compared to the other two types of jasmines. The HTbZIP genes were categorized into 12 subfamilies compared with A. thaliana. In cis-acting element analysis, all genes contained light-responsive elements. The abscisic acid response element (ABRE) was the most abundant in HTbZIP62 promoter, followed by HTbZIP33. Tissue-specific genes of the bZIPs may play a crucial role in regulating the development of jasmine organs and tissues, with HTbZIP36 showing the most significant expressions in roots. Combined with complicated protein interactions, HTbZIP62 and HTbZIP33 might play a crucial role in the ABA signaling pathway and stress tolerance. Combined with RT-qPCR analysis, SJbZIP37/57/62 were more sensitive to ABA response genes compared with other bZIPs in DJ amd HT genomes. Our findings provide a useful resource for further research on the regulation of key genes to improve abiotic stress tolerance in jasmine.

Unlike Panax species, which can produce dammarane-type saponins, Aralia elata can only synthesize oleananetype saponins. Here, the authors reveal that the loss of the dammaranediol synthase-encoding gene and tandem duplication of triterpene biosynthetic genes drive structural divergences of saponins between the two genera.

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

AP2/ERF transcription factors play crucial roles in various biological activities, including plant growth, development, and responses to biotic and abiotic stressors. However, limited research has been conducted on the AP2/ERF genes of Melastoma dodecandrum for breeding of this potential fruit crop. Leveraging the recently published whole genome sequence, we conducted a comprehensive assessment of this superfamily and explored the expression patterns of AP2/ERF genes at a genome-wide level. A significant number of genes, totaling 218, were discovered to possess the AP2 domain sequence and displayed notable structural variations among five subfamilies. An uneven distribution of these genes was observed on 12 pseudochromosomes as the result of gene expansion facilitated by segmental duplications. Analysis of cis-acting elements within promoter sites and 87.6% miRNA splicing genes predicted their involvement in multiple hormone responses and abiotic stresses through transcriptional and post-transcriptional regulations. Transcriptome analysis combined with qRT-PCR results indicated that certain candidate genes are involved in tissue formation and the response to developmental changes induced by IAA hormones. Overall, our study provides valuable insights into the evolution of ERF genes in angiosperms and lays a solid foundation for future breeding investigations aimed at improving fruit quality and enhancing adaptation to barren land environments.

The Elsholtzieae, comprising ca. 7 genera and 70 species, is a small tribe of Lamiaceae (mint family). Members of Elsholtzieae are of high medicinal, aromatic, culinary, and ornamentals value. Despite the rich diversity and value of Elsholtzieae, few molecular markers or plastomes are available for phylogenetics. In the present study, we employed high-throughput sequencing to assemble two Mosla plastomes, M. dianthera and M. scabra, for the first time, and compared with other plastomes of Elsholtzieae. The plastomes of Elsholtzieae exhibited a quadripartite structure, ranging in size from 148,288 bp to 152,602 bp. Excepting the absence of the pseudogene rps19 in Elsholtzia densa, the exhaustive tally revealed the presence of 132 genes (113 unique genes). Among these, 85 protein-coding genes (CDS), 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1) were annotated. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. Notably, the E. eriostchya plastid genome exhibited increased GC content regions in the LSC and SSC, resulting in an increased overall GC content of the entire plastid genome. The E. densa plastid genome displayed modified boundaries due to inverted repeat (IR) contraction. The sequences of CDS and intergenic regions (IGS) with elevated variability were identified as potential molecular markers for taxonomic inquiries within Elsholtzieae. Phylogenetic analysis indicated that four genera formed monophyletic entities, with Mosla and Perilla forming a sister clade. This clade was, in turn, sister to Collinsonia, collectively forming a sister group to Elsholtzia. Both CDS, and CDS + IGS could construct a phylogenetic tree with stronger support. These findings facilitate species identification and DNA barcoding investigations in Elsholtzieae and provide a foundation for further exploration and resource utilization within this tribe.

Epidendrum, one of the three largest genera of Orchidaceae, exhibits significant horticultural and ornamental value and serves as an important research model in conservation, ecology, and evolutionary biology. Given the ambiguous identification of germplasm and complex evolutionary relationships within the genus, the complete plastome of this genus (including five species) were firstly sequenced and assembled to explore their characterizations. The plastomes exhibited a typical quadripartite structure. The lengths of the plastomes ranged from 147,902 bp to 150,986 bp, with a GC content of 37.16% to 37.33%. Gene annotation revealed the presence of 78–82 protein-coding genes, 38 tRNAs, and 8 rRNAs. A total of 25–38 long repeats and 130–149 SSRs were detected. Analysis of relative synonymous codon usage (RSCU) indicated that leucine (Leu) was the most and cysteine (Cys) was the least. The consistent and robust phylogenetic relationships of Epidendrum and its closely related taxa were established using a total of 43 plastid genomes from the tribe Epidendreae. The genus Epidendrum was supported as a monophyletic group and as a sister to Cattleya. Meanwhile, four mutational hotspots (trnCGCA–petN, trnDGUC–trnYGUA, trnSGCU–trnGUCC, and rpl32–trnLUAG) were identified for further phylogenetic studies. Our analysis demonstrates the promising utility of plastomes in inferring the phylogenetic relationships of Epidendrum.

Trichoglottis exhibits a range of rich variations in colors and shapes of flower and is a valuable ornamental orchid genus. The genus Trichoglottis has been expanded by the inclusion of Staurochilus, but this Trichoglottis sensu lato (s.l.) was recovered as a non-monophyletic genus based on molecular sequences from one or a few DNA regions. Here, we present phylogenomic data sets, incorporating complete plastome sequences from seven species (including five species sequenced in this study) of Trichoglottis s.l. (including two species formerly treated as Staurochilus), to compare plastome structure and to reconstruct the phylogenetic relationships of this genus. The seven plastomes possessed the typical quadripartite structure of angiosperms and ranged from 149,402 bp to 149,841 bp with a GC content of 36.6–36.7%. These plastomes contain 120 genes, which comprise 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes, all ndh genes were pseudogenized or lost. A total of 98 (T. philippinensis) to 134 (T. ionosma) SSRs and 33 (T. subviolacea) to 46 (T. ionosma) long repeats were detected. The consistent and robust phylogenetic relationships of Trichoglottis were established using a total of 25 plastid genomes from the Aeridinae subtribe. The genus Trichoglottis s.l. was strongly supported as a monophyletic group, and two species formerly treated as Staurochilus were revealed as successively basal lineages. In addition, five mutational hotspots (trnNGUU-rpl32, trnLUAA, trnSGCU-trnGUCC, rbcL-accD, and trnTGGU-psbD) were identified based on the ranking of PI values. Our research indicates that plastome data is a valuable source for molecular identification and evolutionary studies of Trichoglottis and its related genera.

WUSCHEL-related homeobox (WOX) is a plant-specific transcription factor (TF), which plays an essential role in the regulation of plant growth, development, and abiotic stress responses. However, little information is available on the specific roles of WOX TFs in sacred lotus (Nelumbo nucifera), which is a perennial aquatic plant with important edible, ornamental, and medicinal values. We identified 15 WOX TFs distributing on six chromosomes in the genome of N. nucifera. A total of 72 WOX genes from five species were divided into three clades and nine subclades based on the phylogenetic tree. NnWOXs in the same subclades had similar gene structures and conserved motifs. Cis-acting element analysis of the promoter regions of NnWOXs found many elements enriched in hormone induction, stress responses, and light responses, indicating their roles in growth and development. The Ka/Ks analysis showed that the WOX gene family had been intensely purified and selected in N. nucifera. The expression pattern analysis suggested that NnWOXs were involved in organ development and differentiation of N. nucifera. Furthermore, the protein–protein interaction analysis showed that NnWOXs might participate in the growth, development, and metabolic regulation of N. nucifera. Taken together, these findings laid a foundation for further analysis of NnWOX functions.

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.

Cymbidium, which includes approximately 80 species, is one of the most ornamental and cultivated orchid genera. However, a lack of markers and sparse sampling have posed great challenges to resolving the phylogenetic relationships within the genus. In the present study, we reconstructed the phylogenetic relationships by utilizing one nuclear DNA (nrITS) and seven plastid genes (rbcL, trnS, trnG, matK, trnL, psbA, and atpI) from 70 species (varieties) in Cymbidium. We also examined the occurrence of phylogenetic conflict between nuclear (nrITS) and plastid loci and investigated how phylogenetic conflict bears on taxonomic classification within the genus. We found that phylogenetic conflict and low support values may be explained by hybridization and a lack of informative characteristics. Our results do not support previous classification of the subgenera and sections within Cymbidium. Discordance between gene trees and network analysis indicate that reticulate evolution occurred in the genus Cymbidium. Overall, our study indicates that Cymbidium has undergone a complex evolution.

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

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.

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.

Members of the YABBY gene family play significant roles in lamina development in cotyledons, floral organs, and other lateral organs. The Orchidaceae family is one of the largest angiosperm groups. Some YABBYs have been reported in Orchidaceae. However, the function of YABBY genes in Cymbidium is currently unknown. In this study, 24 YABBY genes were identified in Cymbidium ensifolium, C. goeringii, and C. sinense. We analyzed the conserved domains and motifs, the phylogenetic relationships, chromosome distribution, collinear correlation, and cis-elements of these three species. We also analyzed expression patterns of C. ensifolium and C. goeringii. Phylogenetic relationships analysis indicated that 24 YABBY genes were clustered in four groups, INO, CRC/DL, YAB2, and YAB3/FIL. For most YABBY genes, the zinc finger domain was located near the N-terminus and the helix-loop-helix domain (YABBY domain) near the C-terminus. Chromosomal location analysis results suggested that only C. goeringii YABBY has tandem repeat genes. Almost all the YABBY genes displayed corresponding one-to-one relationships in the syntenic relationships analysis. Cis-elements analysis indicated that most elements were clustered in light-responsive elements, followed by MeJA-responsive elements. Expression patterns showed that YAB2 genes have high expression in floral organs. RT-qPCR analysis showed high expression of CeYAB3 in lip, petal, and in the gynostemium. CeCRC and CeYAB2.2 were highly expressed in gynostemium. These findings provide valuable information of YABBY genes in Cymbidium species and the function in Orchidaceae.

The MYB gene family plays a vital regulatory role in plant metabolism, stress response, and floral color. The R2R3-MYB gene family of C. goeringii was identified, and its expression was analyzed using bioinformatics in this article. The R2R3-MYB genes of Arabidopsis thaliana were used as a reference to determine 104 CgMYB genes and categorize them into 22 subfamilies. Exon/intron organizations and conserved motif analysis revealed that the majority of CgMYB genes were highly conserved, and chromosome localization and collinearity analysis provided evidence of tandem duplication and segmental duplication events, indicating the phenomenon of gene family expansion and contraction. The function of CgMYB genes was analyzed by cis-acting element and gene ontology (GO) enrichment. In addition, we selected CgMYB91 and CgMYB32 for RT–qPCR, suggesting that CgMYB91 and CgMYB32 are associated with anthocyanin formation. In short, this study provides a comprehensive and specific function of the R2R3-MYB transcription factors (TFs) in orchids.

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.

Abstract Background The Orchid family is the largest families of the monocotyledons and an economically important ornamental plant worldwide. Given the pivotal role of this plant to humans, botanical researchers and breeding communities should have access to valuable genomic and transcriptomic information of this plant. Previously, we established OrchidBase, which contains expressed sequence tags (ESTs) from different tissues and developmental stages of Phalaenopsis as well as biotic and abiotic stress-treated Phalaenopsis. The database includes floral transcriptomic sequences from 10 orchid species across all the five subfamilies of Orchidaceae. Description Recently, the whole-genome sequences of Apostasia shenzhenica, Dendrobium catenatum, and Phalaenopsis equestris were de novo assembled and analyzed. These datasets were used to develop OrchidBase 4.0, including genomic and transcriptomic data for these three orchid species. OrchidBase 4.0 offers information for gene annotation, gene expression with fragments per kilobase of transcript per millions mapped reads (FPKM), KEGG pathways and BLAST search. In addition, assembled genome sequences and location of genes and miRNAs could be visualized by the genome browser. The online resources in OrchidBase 4.0 can be accessed by browsing or using BLAST. Users can also download the assembled scaffold sequences and the predicted gene and protein sequences of these three orchid species. Conclusions OrchidBase 4.0 is the first database that contain the whole-genome sequences and annotations of multiple orchid species. OrchidBase 4.0 is available at http://orchidbase.itps.ncku.edu.tw/

The wood sorrel family, Oxalidaceae, is mainly composed of annual or perennial herbs, a few shrubs, and trees distributed from temperate to tropical zones. Members of Oxalidaceae are of high medicinal, ornamental, and economic value. Despite the rich diversity and value of Oxalidaceae, few molecular markers or plastomes are available for phylogenetic analysis of the family. Here, we reported four new whole plastomes of Oxalidaceae and compared them with plastomes of three species in the family, as well as the plastome of Rourea microphylla in the closely related family Connaraceae. The eight plastomes ranged in length from 150,673 bp (Biophytum sensitivum) to 156,609 bp (R. microphylla). Genome annotations revealed a total of 129–131 genes, including 83–84 protein-coding genes, eight rRNA genes, 37 tRNA genes, and two to three pseudogenes. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. The smaller plastomes of herbs B. sensitivum and three Oxalis species are associated with variations in IR region sizes, intergenic region variation, and gene or intron loss. We identified sequences with high variation that may serve as molecular markers in taxonomic studies of Oxalidaceae. The phylogenetic trees of selected superrosid representatives based on 76 protein-coding genes corroborated the Oxalidaceae position in Oxalidales and supported it as a sister to Connaraceae. Our research also supported the monophyly of the COM (Celastrales, Oxalidales, and Malpighiales) clade.

Cao, Li, et al. sequence 13 perennial woody plant species of Lycium, and specifically provide a draft assembly of L. ruthenicum and a chromosome-level assembly of L. barbarum, the wolfberry or Goji berry. From a phylogenetic tree the authors identify an ancient hexaploidization event, and report the evolution of gene families including fruit ripening, fruit coloration, polysaccharide synthesis and self-incompatibility within Solanaceae and the general biogeography of L. barbarum.

Abstract Background Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum. Results Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 – 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates. Conclusions We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

Apostasia shenzhenica belongs to the subfamily Apostasioideae and is a primitive group located at the base of the Orchidaceae phylogenetic tree. However, the A. shenzhenica mitochondrial genome (mitogenome) is still unexplored, and the phylogenetic relationships between monocots mitogenomes remain unexplored. In this study, we discussed the genetic diversity of A. shenzhenica and the phylogenetic relationships within its monocotyledon mitogenome. We sequenced and assembled the complete mitogenome of A. shenzhenica, resulting in a circular mitochondrial draft of 672,872 bp, with an average read coverage of 122× and a GC content of 44.4%. A. shenzhenica mitogenome contained 36 protein-coding genes, 16 tRNAs, two rRNAs, and two copies of nad4L. Repeat sequence analysis revealed a large number of medium and small repeats, accounting for 1.28% of the mitogenome sequence. Selection pressure analysis indicated high mitogenome conservation in related species. RNA editing identified 416 sites in the protein-coding region. Furthermore, we found 44 chloroplast genomic DNA fragments that were transferred from the chloroplast to the mitogenome of A. shenzhenica, with five plastid-derived genes remaining intact in the mitogenome. Finally, the phylogenetic analysis of the mitogenomes from A. shenzhenica and 28 other monocots showed that the evolution and classification of most monocots were well determined. These findings enrich the genetic resources of orchids and provide valuable information on the taxonomic classification and molecular evolution of monocots.