The Complete Mitogenome of Apostasia fujianica Y.Li & S.Lan and Comparative Analysis of Mitogenomes across Orchidaceae.

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.

Genome-Based Identification of the Dof Gene Family in Three Cymbidium Species and Their Responses to Heat Stress in Cymbidium goeringii.

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.

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

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

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree).

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Genome-Wide Identification and Expression Analysis of the GRAS Gene Family and Their Responses to Heat Stress in Cymbidium goeringii.

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.

On the mechanism of wogonin against acute monocytic leukemia using network pharmacology and experimental validation.

Wogonin is a natural flavone compound from the plant Scutellaria baicalensis, which has a variety of pharmacological activities such as anti-cancer, anti-virus, anti-inflammatory, and immune regulation. However, the potential mechanism of wogonin remains unknown. This study was to confirm the molecular mechanism of wogonin for acute monocytic leukemia treatment, known as AML-M5. The potential action targets between wogonin and acute monocytic leukemia were predicted from databases. The compound-target-pathway network and protein-protein interaction network (PPI) were constructed. The enrichment analysis of related targets and molecular docking were performed. The network pharmacological results of wogonin for AML-M5 treatment were verified using the THP-1 cell line. 71 target genes of wogonin associated with AML-M5 were found. The key genes TP53, SRC, AKT1, RELA, HSP90AA1, JUN, PIK3R1, and CCND1 were preliminarily found to be the potential central targets of wogonin for AML-M5 treatment. The PPI network analysis, GO analysis and KEGG pathway enrichment analysis demonstrated that the PI3K/AKT signaling pathway was the significant pathway in the wogonin for AML-M5 treatment. The antiproliferative effects of wogonin on THP-1 cells of AML-M5 presented a dose-dependent and time-dependent manner, inducing apoptosis, blocking the cell cycle at the G2/M phase, decreasing the expressions of CCND1, CDK2, and CyclinA2 mRNA, as well as AKT and p-AKT proteins. The mechanisms of wogonin on AML-M5 treatment may be associated with inhibiting cell proliferation and regulating the cell cycle via the PI3K/AKT signaling pathway.

Genome-Wide Identification and Expression Pattern Analysis of TIFY Family Genes Reveal Their Potential Roles in Phalaenopsis aphrodite Flower Opening.

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.

Genome-wide characterization and expression profiling of the HD-ZIP gene family in Acoraceae under salinity and cold stress.

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 Evolution of the WUSCHEL-Related Homeobox Gene Family in Dendrobium Species and Its Role in Sex Organ Development in D. chrysotoxum.

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

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae).

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Statistical Genomics Analysis of Simple Sequence Repeats from the Paphiopedilum Malipoense Transcriptome Reveals Control Knob Motifs Modulating Gene Expression.

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.

Comparative Analysis of Six Complete Plastomes of Tripterospermum spp.

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.

The Complete Chloroplast Genomes of Bulbophyllum (Orchidaceae) Species: Insight into Genome Structure Divergence and Phylogenetic Analysis.

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.

Selection and Validation of qRT-PCR Internal Reference Genes to Study Flower Color Formation in Camellia impressinervis.

Real-time quantitative PCR (qRT-PCR) is a pivotal technique for gene expression analysis. To ensure reliable and accurate results, the internal reference genes must exhibit stable expression across varied experimental conditions. Currently, no internal reference genes for Camellia impressinervis have been established. This study aimed to identify stable internal reference genes from eight candidates derived from different developmental stages of C. impressinervis flowers. We employed geNorm, NormFinder, and BestKeeper to evaluate the expression stability of these candidates, which was followed by a comprehensive stability analysis. The results indicated that CiTUB, a tubulin gene, exhibited the most stable expression among the eight reference gene candidates in the petals. Subsequently, CiTUB was utilized as an internal reference for the qRT-PCR analysis of six genes implicated in the petal pigment synthesis pathway of C. impressinervis. The qRT-PCR results were corroborated by transcriptome sequencing data, affirming the stability and suitability of CiTUB as a reference gene. This study marks the first identification of stable internal reference genes within the entire genome of C. impressinervis, establishing a foundation for future gene expression and functional studies. Identifying such stable reference genes is crucial for advancing molecular research on C. impressinervis.

Genome-Wide Identification and Drought Stress Response Pattern of the NF-Y Gene Family in Cymbidium sinense.

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.

Wolfberry genome database: integrated genomic datasets for studying molecular biology.

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.

Organelle Genomes of Epipogium roseum Provide Insight into the Evolution of Mycoheterotrophic Orchids.

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.

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

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