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.

Molecular insights into self-incompatibility systems: From evolution to breeding.

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes.

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 (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > 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.

Diploid and tetraploid genomes of Acorus and the evolution of monocots.

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

Apostasia Mitochondrial Genome Analysis and Monocot Mitochondria Phylogenomics.

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.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism.

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

Comparative Analysis of Luisia (Aeridinae, Orchidaceae) Plastomes Shed Light on Plastomes Evolution and Barcodes Investigation.

Luisia, a genus of the subtribe Aeridinae of Orchidaceae, comprises ca. 40 species. Members of Luisia exhibit unique morphological characteristics and represent a valuable ornamental orchid genus. However, due to the scarcity of distinct morphological characters, species identification within this genus is ambiguous and controversial. In the present study, next-generation sequencing (NGS) methods were used to assemble the plastomes of five Luisia species and compare them with one publicly available Luisia plastid genome data. The plastomes of Luisia possessed a quadripartite structure, with sizes ranging from 146,243 bp to 147,430 bp. The plastomes of six Luisia species contained a total of 120 genes, comprising 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. Notably, all ndh genes were pseudogenized or lost. An analysis of codon usage bias showed that leucine (Leu) exhibited the highest frequency, while cysteine (Cys) exhibited the lowest frequency. A total of 57 to 64 SSRs and 42 to 49 long repeats were identified. Five regions and five coding sequences were identified for DNA barcodes, based on the nucleotide diversity (Pi) analysis. The species of Luisia constituted a monophyletic group and were sister to Paraphalaenopsis with strong support. Our study deepens the understanding of species identification, plastome evolution and the phylogenetic positions of Luisia.

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

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

OrchidBase 5.0: updates of the orchid genome knowledgebase.

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

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

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

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Cymbidium goeringii.

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.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy.

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata.

Araliaceae species produce various classes of triterpene and triterpenoid saponins, such as the oleanane-type triterpenoids in Aralia species and dammarane-type saponins in Panax, valued for their medicinal properties. The lack of genome sequences of Panax relatives has hindered mechanistic insight into the divergence of triterpene saponins in Araliaceae. Here, we report a chromosome-level genome of Aralia elata with a total length of 1.05 Gb. The loss of 12 exons in the dammarenediol synthase (DDS)-encoding gene in A. elata after divergence from Panax might have caused the lack of dammarane-type saponin production, and a complementation assay shows that overexpression of the PgDDS gene from Panax ginseng in callus of A. elata recovers the accumulation of dammarane-type saponins. Tandem duplication events of triterpene biosynthetic genes are common in the A. elata genome, especially for AeCYP72As, AeCSLMs, and AeUGT73s, which function as tailoring enzymes of oleanane-type saponins and aralosides. More than 13 aralosides are de novo synthesized in Saccharomyces cerevisiae by overexpression of these genes in combination. This study sheds light on the diversity of saponins biosynthetic pathway in Araliaceae and will facilitate heterologous bioproduction of aralosides.

The Melastoma dodecandrum genome and the evolution of Myrtales.

Melastomataceae has abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size is 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis, and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. The SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

[History and application of foreign medicinal resources serving healthcare industry].

Foreign medicinal resources have always been an important part of Chinese medicine and have made great contributions to the development of the Chinese medical and health industry. Since the opening of the Silk Road in the Han Dynasty, foreign medicinal resources have been introduced for different purposes, some of which have become Chinese medicine in clinical practice and are still in use today.Today, foreign medicinal resources also serve the Big Health industry in China. They are introduced and applied to the fields in the Big Health industry, such as food, cosmetics, health products, decoction pieces, and daily chemical products. With the integration and development of the "Healthy China" initiative, more foreign resources will enter the Big Health industry. This paper retrospectedthe history of foreign medicinal resources serving the ancient medical and health industry, reviewedits current development under the Big Health industry, summarizedthe experience of foreign medicinal resources serving the ancient medical and health industry, as well as the development and problems of new foreign medicinal resources, and put forward some suggestions to provide ideas for the development and application of foreign medicinal resources under the Big Health industry.

The complete chloroplast genome sequence of Phalaenopsis wilsoniii (Orchidaceae).

In the present study, we reported and characterized the complete chloroplast genome of a moth orchid, Phalaenopsis wilsonii, which is endemic to South China. Its plastid genome size is 145,373 bp, consisting of a large single copy (LSC) region (84,996 bp), a small single-copy region (10,668 bp), and two inverted repeats (IRs) regions (24,855 bp). A total of 122 plastid genes were annotated, comprising 76 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. The phylogenetic tree further revealed that P. wilsonii showed a sister relationship with P. lowii within subgenus Parishianae.

Plastid phylogenomics improves resolution of phylogenetic relationship in the Cheirostylis and Goodyera clades of Goodyerinae (Orchidoideae, Orchidaceae).

Goodyerinae are one of phylogenetically unresolved groups of Orchidaceae. The lack of resolution achieved through the analyses of previous molecular sequences from one or a few markers has long confounded phylogenetic estimation and generic delimitation. Here, we present large-scale phylogenomic data to compare the plastome structure of the two main clades (Goodyera and Cheirostylis) in this subtribe and further adopt two strategies, combining plastid coding sequences and the whole plastome, to investigate phylogenetic relationships. A total of 46 species in 16 genera were sampled, including 39 species in 15 genera sequenced in this study. The plastomes of heterotrophic species are not drastically reduced in overall size, but display a pattern congruent with a loss of photosynthetic function. The plastomes of autotrophic species ranged from 147 to 165 kb and encoded from 132 to 137 genes. Three unusual structural features were detected: a 1.0-kb inversion in the large single-copy region of Goodyera schlechtendaliana; the loss and/or pseudogenization of ndh genes only in two species, Cheirostylis chinensis and C. montana; and the expansion of inverted repeat regions and contraction of small single-copy region in Hetaeria oblongifolia. Phylogenomic analyses provided improved resolution for phylogenetic relationships. All genera were recovered as monophyletic, except for Goodyera and Hetaeria, which were each recovered as non-monophyletic. Nomenclatural changes are needed until the broader sampling and biparental inherited markers. This study provides a phylogenetic framework of Goodyerinae and insight into plastome evolution of Orchidaceae.

The complete plastid genome of Thrixspermum centipeda (Orchidaceae, Aeridinae).

The complete plastid genome of the type species of Thrixspermum, Th. centipeda, was determined and analyzed in this work. The plastome was 147,888 bp in length with 85,899 bp of the large single-copy (LSC) region, 11,055 bp of the small single-copy (SSC) region and 25,467 bp of the invert repeats (IR) regions. The genome contained 120 genes, including 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis divided 18 Aeridinae plastomes into four groups, and Th. centipeda was sister to Th. tsii.