Analysis of codon usage bias and evolution in the chloroplast genome of Mesona chinensis Benth

Analysis of codon usage patterns in complete plastomes of four medicinal Polygonatum species (Asparagaceae)

Polygonati Rhizoma and Polygonati odorati Rhizoma, known as "Huangjing" and "Yuzhu" in China, are medicinal Polygonatum species resources with top-grade medical and edible properties. The chloroplast (cp) genome has been used to study species diversity, evolution, and breeding of species for applications in genetic engineering. Codon usage bias (CUB), a common and complex natural phenomenon, is essential for studies of codon optimization of exogenous genes, genetic engineering, and molecular evolution. However, the CUB of medicinal Polygonatum species chloroplast genomes has not been systematically studied. In our study, a detailed analysis of CUB was performed in the medicinal Polygonatum species chloroplast genomes. We investigated the codon bias of 204 plastid protein-coding genes (PCGs) in 4 medicinal Polygonatum species using CodonW and CUSP online software. Through the analysis of the codon bias index, we found that the medicinal Polygonatum species chloroplast genomes had weak codon usage bias. In addition, our results also showed a high preference for AT bases in medicinal Polygonatum species chloroplast genomes, and the preference to use AT-ending codons was observed in these species chloroplast genomes. The neutrality plot, ENC plot, PR2-Bias plot, and correspondence analysis showed that compared with mutation pressure, natural selection was the most important factor of CUB. Based on the comparative analysis of high-frequency codons and high expression codons, we also determined the 10-11 optimal codons of investigative medicinal Polygonatum species. Furthermore, the result of RSCU-based cluster analysis showed that the genetic relationship between different medicinal Polygonatum species could be well reflected. This study provided an essential understanding of CUB and evolution in the medicinal Polygonatum species chloroplast genomes.

Comparative analysis of codon usage bias in the chloroplast genomes of eighteen Ampelopsideae species (Vitaceae)

BACKGROUND

The tribe Ampelopsideae plants are important garden plants with both medicinal and ornamental values. The study of codon usage bias (CUB) facilitates a deeper comprehension of the molecular genetic evolution of species and their adaptive strategies. The joint analysis of CUB in chloroplast genomes (cpDNA) offers valuable insights for in-depth research on molecular genetic evolution, biological resource conservation, and elite breeding within this plant family.

RESULTS

The base composition and codon usage preferences of the eighteen chloroplast genomes were highly similar, with the GC content of bases at all positions of their codons being less than 50%. This indicates that they preferred A/T bases. Their effective codon numbers were all in the range of 35-61, which indicates that the codon preferences of the chloroplast genomes of the 18 Ampelopsideae plants were relatively weak. A series of analyses indicated that the codon preference of the chloroplast genomes of the 18 Ampelopsideae plants was influenced by a combination of multiple factors, with natural selection being the primary influence. The clustering tree generated based on the relative usage of synonymous codons is consistent with some of the results obtained from the phylogenetic tree of chloroplast genomes, which indicates that the clustering tree based on the relative usage of synonymous codons can be an important supplement to the results of the sequence-based phylogenetic analysis. Eventually, 10 shared best codons were screened on the basis of the chloroplast genomes of 18 species.

CONCLUSION

The codon preferences of the chloroplast genome in Ampelopsideae plants are relatively weak and are primarily influenced by natural selection. The codon composition of the chloroplast genomes of the eighteen Ampelopsideae plants and their usage preferences were sufficiently similar to demonstrate that the chloroplast genomes of Ampelopsideae plants are highly conserved. This study provides a scientific basis for the genetic evolution of chloroplast genes in Ampelopsideae species and their suitable strategies.

Repetitive Sequences, Codon Usage Bias and Phylogenetic Analysis of the Plastome of Miliusa glochidioides

Annonaceae is the largest family in Magnoliales, exhibiting the greatest diversity among and within genera. In this study, we conducted an analysis of repetitive sequences and codon usage bias in the previously acquired plastome of Miliusa glochidioides. Using a concatenated dataset of shared genes, we constructed the phylogenetic relationships among 27 Annonaceae species. The results showed that the size of the plastomes in the Annonaceae ranged from 159 to 202 kb, with the size of the inverted repeat region ranging from 40 to 65 kb. Within the plastome of M. glochidioides, we identified 42 SSRs, 36 tandem repeats, and 9 dispersed repeats. These SSRs consist of three nucleotide types and eight motif types, with a preference for A/T bases, primarily located in the large single-copy regions and intergenic spacers. Tandem and dispersed repeat sequences were predominantly detected in the IR region. Through codon usage bias analysis, we identified 30 high-frequency codons and 11 optimal codons. The plastome of M. glochidioides demonstrated relatively weak codon usage bias, favoring codons with A/T endings, primarily influenced by natural selection. Phylogenetic analysis revealed that all four subfamilies formed monophyletic groups, with Cananga odorata (Ambavioideae) and Anaxagorea javanica (Anaxagoreoideae) successively nested outside Annonoideae + Malmeoideae. These findings improve our understanding of the plastome of M. glochidioides and provide additional insights for studying plastome evolution in Annonaceae.

Comparative genomic and phylogenetic analyses of mitochondrial genomes of hawthorn (Crataegus spp.) in Northeast China

Hawthorn (Crataegus spp.) plants are major sources of health food and medicines. Twenty species and seven variations of Crataegus are present in China. A variety of unique Crataegus species was found in their natural distribution in northeast China. In the present study, we assembled and annotated the mitochondrial genomes of five Crataegus species from northeastern China. The sizes of the newly sequenced mitochondrial genomes ranged from 245,907 bp to 410,837 bp. A total of 45-55 genes, including 12-19 transfer RNA genes, three ribosomal RNA genes, and 29-33 protein-coding genes (PCGs) were encoded by these mitochondrial genomes. Seven divergent hotspot regions were identified by comparative analyses: atp6, nad3, ccmFN, matR, nad1, nad5, and rps1. The most conserved genes among the Crataegus species, according to the whole-genome correlation analysis, were nad1, matR, nad5, ccmFN, cox1, nad4, trnQ-TTG, trnK-TTT, trnE-TTC, and trnM-CAT. Horizontal gene transfer between organellar genomes was common in Crataegus plants. Based on the phylogenetic trees of mitochondrial PCGs, C. maximowiczii, C. maximowiczii var. ninganensis, and C. bretschneideri shared similar maternal relationships. This study improves Crataegus mitochondrial genome resources and offers important insights into the taxonomy and species identification of this genus.

Preparation technologies, structural features, and biological activities of polysaccharides from Mesona chinensis Benth.: A review

ETHNOPHARMACOLOGICAL RELEVANCE

Mesona chinensis Benth. (or Platostoma palustre (Blume) A. J. Paton) is an important medicinal and edible plant also known as the Hsian-tsao in China and Southeast Asian countries. It is cold in nature and sweet in taste, with the effects of clearing heat, relieving heatstroke and diuretic, and traditionally used to treat heatstroke, erysipelas, hypertension, joint pain and other diseases in folk medicine. It is also a popular supplement with the function of detoxifying and heat-clearing use in Asia. It is used to be processed into the popular tea, Bean jelly, and so on. Published studies have demonstrated that polysaccharides from M. chinensis (MCPs) are one of the principal bioactive ingredients with a variety of health-promoting effects in the prevention and treatment of diseases, including antioxidant, immunomodulation, anti-inflammatory, hepatoprotective, anti-tumor, hypoglycemic, regulation of gut microbiota, and other pharmacological properties.

AIM OF THE REVIEW

This review aims to compile the extraction and purification methods, structural characteristics, pharmacological activities including the mechanism of action of MCPs, and to further understand the applications of M. chinensis in order to lay the foundation for the development of MCPs.

MATERIALS AND METHODS

By inputting the search term "Mesona chinensis polysaccharides", relevant research information was obtained from databases such as PubMed, Google Scholar, Web of Science, and China National Knowledge Infrastructure (CNKI).

RESULTS

More than 40 polysaccharides have been extracted from M. chinensis, different extraction and purification methods have been described, as well as the structural features and pharmacological activities of MCPs have been systematically reviewed. Polysaccharides, as important components of M. chinensis, were mainly extracted by methods such as hot water dipping method, hot alkali extraction method, enzyme-assisted extraction method and ultrasonic-assisted extraction method, subsequently obtained by decolorization, deproteinization, removal of other small molecules and separation on various chromatographic columns. The chemical composition and structure of MCPs show diversity and have a variety of pharmacological activities, including antioxidant, immunomodulation, anti-inflammatory, hepatoprotective, anti-tumor, hypoglycemic, regulation of gut microbiota, and so on.

CONCLUSIONS

This article systematically reviews the research progress of MCPs in terms of extraction and purification, structural characteristics, rheological gel properties, pharmacological properties, and safety assessment. The potentials and roles of M. chinensis in the field of medicine, functional food, and materials are further highlighted to provide references and bases for the high-value processing and utilization of MCPs.

Volatilome and flavor analyses based on e-nose combined with HS-GC-MS provide new insights into ploidy germplasm diversity in Platostoma palustre

Platostoma palustre (Mesona chinensis Benth or Hsian-tsao, also known as "Xiancao" in China), is an edible and medicinal plant native to India, Myanmar, and Indo-China. It is the main ingredient in the popular desserts Hsian-tsao tea, herbal jelly, and sweet herbal jelly soup. P. palustre is found abundantly in nutrient-rich substances and possesses unique aroma compounds. Variations in the contents of volatile compounds among different germplasms significantly affect the quality and flavor of P. palustre, causing contradiction in demand. This study investigates the variation in the volatile compound profiles of distinct ploidy germplasms of P. palustre by utilising headspace gas chromatography-mass spectrometry (HS-GC-MS) and an electronic nose (e-nose). The results showed significant differences in the aroma characteristics of stem and leaf samples in diverse P. palustre germplasms. A total of sixty-seven volatile compounds have been identified and divided into ten classes. Six volatile compounds (caryophyllene, α-bisabolol, benzaldehyde, β-selinene, β-elemene and acetic acid) were screened as potential marker volatile compounds to discriminate stems and leaves of P. palustre. In this study, leaves of P. palustre showed one odor pattern and stems showed two odor patterns under the influence of α-bisabolol, acetic acid, and butyrolactone. In addition, a correlation analysis was conducted on the main volatile compounds identified by HS-GC-MS and e-nose. This analysis provided additional insight into the variations among samples resulting from diverse germplasms. The present study provides a valuable volatilome, and flavor, and quality evaluation for P. palustre, as well as new insights and scientific basis for the development and use of P. palustre germplasm resources.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

Functionality and Health-Promoting Properties of Polysaccharide and Plant-Derived Substances from Mesona chinensis

Mesona chinensis, in Thai called Chao Kuay and in Chinese Hsian-tsao, belongs to the Lamiaceae family. This herbal plant grows widely in Southern China, Taiwan (China), Malaysia, the Philippines, Indonesia, Vietnam, and Thailand. The Mesona plant is used to make functional products such as drinks and soft textured sweet treats, and also traditional medicine, to treat heat stroke, high blood pressure, heart attack, high blood sugar, hepatic diseases, colon diseases, inflammatory conditions, and to alleviate myalgia. The proximate composition of M. chinensis is a mixture of protein, fat, fiber, ash, and minerals. The main biological compounds in M. chinensis extracts are polysaccharides, terpenoids, flavonoids, and polyphenols, with wide-ranging pharmacological properties including antioxidant, antidiabetic, antilipidemic, carcinoma-inhibitory, renal-protective, antihypertensive, DNA damage-protective, and anti-inflammatory effects. This review investigated the proximate composition, polysaccharide type, and pharmacological properties of M. chinensis extracts. Phytochemical properties enhance the actions of the gut microbiota and improve health benefits. This review assessed the functional and medicinal activities of M. chinensis extracts. Future studies should further elucidate the in vitro/in vivo mechanisms of this plant extract and its impact on gut health.

Codon usage characterization and phylogenetic analysis of the mitochondrial genome in Hemerocallis citrina

BACKGROUND

Hemerocallis citrina Baroni is a traditional vegetable crop widely cultivated in eastern Asia for its high edible, medicinal, and ornamental value. The phenomenon of codon usage bias (CUB) is prevalent in various genomes and provides excellent clues for gaining insight into organism evolution and phylogeny. Comprehensive analysis of the CUB of mitochondrial (mt) genes can provide rich genetic information for improving the expression efficiency of exogenous genes and optimizing molecular-assisted breeding programmes in H. citrina.

RESULTS

Here, the CUB patterns in the mt genome of H. citrina were systematically analyzed, and the possible factors shaping CUB were further evaluated. Composition analysis of codons revealed that the overall GC (GCall) and GC at the third codon position (GC3) contents of mt genes were lower than 50%, presenting a preference for A/T-rich nucleotides and A/T-ending codons in H. citrina. The high values of the effective number of codons (ENC) are indicative of fairly weak CUB. Significant correlations of ENC with the GC3 and codon counts were observed, suggesting that not only compositional constraints but also gene length contributed greatly to CUB. Combined ENC-plot, neutrality plot, and Parity rule 2 (PR2)-plot analyses augmented the inference that the CUB patterns of the H. citrina mitogenome can be attributed to multiple factors. Natural selection, mutation pressure, and other factors might play a major role in shaping the CUB of mt genes, although natural selection is the decisive factor. Moreover, we identified a total of 29 high-frequency codons and 22 optimal codons, which exhibited a consistent preference for ending in A/T. Subsequent relative synonymous codon usage (RSCU)-based cluster and mt protein coding gene (PCG)-based phylogenetic analyses suggested that H. citrina is close to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in evolutionary terms, reflecting a certain correlation between CUB and evolutionary relationships.

CONCLUSIONS

There is weak CUB in the H. citrina mitogenome that is subject to the combined effects of multiple factors, especially natural selection. H. citrina was found to be closely related to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in terms of their evolutionary relationships as well as the CUB patterns of their mitogenomes. Our findings provide a fundamental reference for further studies on genetic modification and phylogenetic evolution in H. citrina.

Analysis of codon usage patterns in 48 Aconitum species

BACKGROUND

The Aconitum genus is a crucial member of the Ranunculaceae family. There are 350 Aconitum species worldwide, with about 170 species found in China. These species are known for their various pharmacological effects and are commonly used to treat joint pain, cold abdominal pain, and other ailments. Codon usage bias (CUB) analysis contributes to evolutionary relationships and phylogeny. Based on protein-coding sequences (PCGs), we selected 48 species of Aconitum for CUB analysis.

RESULTS

The results revealed that Aconitum species had less than 50% GC content. Furthermore, the distribution of GC content was irregular and followed a trend of GC1 > GC2 > GC3, indicating a bias towards A/T bases. The relative synonymous codon usage (RSCU) heat map revealed the presence of conservative codons with slight variations within the genus. The effective number of codons (ENC)-Plot and the parity rule 2 (PR2)-bias plot analysis indicate that natural selection is the primary factor influencing the variation in codon usage. As a result, we screened various optimal codons and found that A/T bases were preferred as the last codon. Furthermore, our Maximum Likelihood (ML) analysis based on PCGs among 48 Aconitum species yielded results consistent with those obtained from complete chloroplast (cp.) genome data. This suggests that analyzing mutation in PCGs is an efficient method for demonstrating the phylogeny of species at the genus level.

CONCLUSIONS

The CUB analysis of 48 species of Aconitum was mainly influenced by natural selection. This study reveals the CUB pattern of Aconitum and lays the foundation for future genetic modification and phylogenetic analyses.

Comparative analysis of the medicinal plant Polygonatum kingianum (Asparagaceae) with related verticillate leaf types of the Polygonatum species based on chloroplast genomes

Background

Polygonatum kingianum has been widely used as a traditional Chinese medicine as well as a healthy food. Because of its highly variable morphology, this medicinal plant is often difficult to distinguish from other related verticillate leaf types of the Polygonatum species. The contaminants in P. kingianum products not only decrease the products' quality but also threaten consumer safety, seriously inhibiting the industrial application of P. kingianum.

Methods

Nine complete chloroplast (cp) genomes of six verticillate leaf types of the Polygonatum species were de novo assembled and systematically analyzed.

Results

The total lengths of newly sequenced cp genomes ranged from 155,437 to 155,977 bp, including 86/87 protein-coding, 38 tRNA, and 8 rRNA genes, which all exhibited well-conserved genomic structures and gene orders. The differences in the IR/SC (inverted repeats/single-copy) boundary regions and simple sequence repeats were detected among the verticillate leaf types of the Polygonatum cp genomes. Comparative cp genomes analyses revealed that a higher similarity was conserved in the IR regions than in the SC regions. In addition, 11 divergent hotspot regions were selected, providing potential molecular markers for the identification of the Polygonatum species with verticillate leaf types. Phylogenetic analysis indicated that, as a super barcode, plastids realized a fast and efficient identification that clearly characterized the relationships within the verticillate leaf types of the Polygonatum species. In brief, our results not only enrich the data on the cp genomes of the genus Polygonatum but also provide references for the P. kingianum germplasm resource protection, herbal cultivation, and drug production.

Conclusion

This study not only accurately identifies P. kingianum species, but also provides valuable information for the development of molecular markers and phylogenetic analyses of the Polygonatum species with verticillate leaf types.

Comparative Analysis of Codon Usage Patterns in Nuclear and Chloroplast Genome of Dalbergia (Fabaceae)

The Dalbergia plants are widely distributed across more than 130 tropical and subtropical countries and have significant economic and medicinal value. Codon usage bias (CUB) is a critical feature for studying gene function and evolution, which can provide a better understanding of biological gene regulation. In this study, we comprehensively analyzed the CUB patterns of the nuclear genome, chloroplast genome, and gene expression, as well as systematic evolution of Dalbergia species. Our results showed that the synonymous and optimal codons in the coding regions of both nuclear and chloroplast genome of Dalbergia preferred ending with A/U at the third codon base. Natural selection was the primary factor affecting the CUB features. Furthermore, in highly expressed genes of Dalbergia odorifera, we found that genes with stronger CUB exhibited higher expression levels, and these highly expressed genes tended to favor the use of G/C-ending codons. In addition, the branching patterns of the protein-coding sequences and the chloroplast genome sequences were very similar in the systematic tree, and different with the cluster from the CUB of the chloroplast genome. This study highlights the CUB patterns and features of Dalbergia species in different genomes, explores the correlation between CUB preferences and gene expression, and further investigates the systematic evolution of Dalbergia, providing new insights into codon biology and the evolution of Dalbergia plants.

Codon Usage Analyses Reveal the Evolutionary Patterns among Plastid Genes of Saxifragales at a Larger-Sampling Scale

Saxifragales is a 15-family order of early-divergent Eudicots with a rich morphological diversity and an ancient rapid radiation. Codon usage bias (CUB) analyses have emerged as an essential tool for understanding the evolutionary dynamics in genes. Thus far, the codon utilization patterns had only been reported in four separate genera within Saxifragales. This study provides a comprehensive assessment of the codon manipulation based on 50 plastid genes, covering 11 constituent families at a larger sampling scale. Our results first showed a high preference for AT bases and AT-ending codons. We then used effective number of codons (ENC) to assess a range of codon bias levels in the plastid genes. We also detected high-informative intrafamilial differences of ENC in three families. Subsequently, parity rule 2 (PR2) plot analyses revealed both family-unique and order-shared bias patterns. Most importantly, the ENC plots and neutrality analyses collectively supported the dominant roles of selection in the CUB of Saxifragales plastid genes. Notably, the phylogenetic affinities inferred by both ML and BI methods were consistent with each other, and they all comprised two primary clades and four subclades. These findings significantly enhance our understanding of the evolutionary processes of the Saxifrage order, and could potentially inspire more CUB analyses at higher taxonomic levels.

Codon usage bias and genetic diversity in chloroplast genomes of Elaeagnus species (Myrtiflorae: Elaeagnaceae)

Codon usage bias (CUB) reveals the characteristics of species and can be utilized to understand their evolutionary relationship, increase the target genes' expression in the heterologous receptor plants, and further provide theoretic assistance for correlative study on molecular biology and genetic breeding. The chief aim of this work was to analyze the CUB in chloroplast (cp.) genes in nine Elaeagnus species to provide references for subsequent studies. The codons of Elaeagnus cp. genes preferred to end with A/T bases rather than with G/C bases. Most of the cp. genes were prone to mutation, while the rps7 genes were identical in sequences. Natural selection was inferred to have a powerful impact on the CUB in Elaeagnus cp. genomes, and their CUB was extremely strong. In addition, the optimal codons were identified in the nine cp. genomes based on the relative synonymous codon usage (RSCU) values, and the optimal codon numbers were between 15 and 19. The clustering analyses based on RSCU were contrasted with the maximum likelihood (ML)-based phylogenetic tree derived from coding sequences, suggesting that the t-distributed Stochastic Neighbor Embedding clustering method was more appropriate for evolutionary relationship analysis than the complete linkage method. Moreover, the ML-based phylogenetic tree based on the conservative matK genes and the whole cp. genomes had visible differences, indicating that the sequences of specific cp. genes were profoundly affected by their surroundings. Following the clustering analysis, Arabidopsis thaliana was considered the optimal heterologous expression receptor plant for the Elaeagnus cp. genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01289-6.

Complete chloroplast genome molecular structure, comparative and phylogenetic analyses of Sphaeropteris lepifera of Cyatheaceae family: a tree fern from China

Sphaeropteris lepifera is a tree fern in the Cyatheaceae, a family that has played an important role in the evolution of plant systems. This study aimed to analyze the complete chloroplast genome of S. lepifera and compared it with previously published chloroplast genomes Cyatheaceae family. The chloroplast genome of S. lepifera comprised 162,114 bp, consisting of a large single copy (LSC) region of 86,327 bp, a small single copy (SSC) region of 27,731 bp and a pair of inverted repeats (IRa and IRb) of 24,028 bp each. The chloroplast genome encoded 129 genes, comprising 32 transfer RNAs, 8 ribosomal RNAs, and 89 protein-coding genes. Comparison of the genomes of 7 Cyatheaceae plants showed that the chloroplast genome of S. lepifera was missing the gene trnV-UAC. Expansion of the SSC region led to the difference in the chloroplast genome size of S. lepifera. Eight genes, atpI, ccsA, petA, psaB, rpl16, rpoA, rpoC1, and ycf2 have high nucleic acid diversity and can be regarded as potential molecular markers. The genes trnG-trnR and atpB were suitable for DNA barcodes between different communities of S. lepifera. The S. lepifera groups in Zhejiang Province probably diffused from Pingtan and Ningde, Fujian. The results will provide a basis for species identification, biological studies, and endangerment mechanism of S. lepifera.

Mitochondrial genome characteristics and phylogenetic analysis of the medicinal and edible plant Mesona chinensis Benth

Mesona chinensis Benth (MCB) (or Platostoma palustre or Platostoma chinense) is an important edible and medicinal plant in China. However, the mitochondrial genome (mitogenome, or mtDNA) of MCB has not been characterized or reported yet. In this study, we first sequenced and characterized the complete mitogenome of MCB. The MCB mitogenome was 494,599 bp in length and encoded 59 genes containing 37 protein-coding genes (PCGs), 19 tRNAs, and 3 rRNAs. Gene transfer analysis revealed that a total of 12 transfer segments with more than 93% identity (total length of 25,427 bp) were detected in the MCB mitogenome. Simple sequence repeats (SSR) analysis showed that 212 simple sequence repeats (SSR) were identified. Repeat sequence analysis revealed 305 repeat sequences (158 forward and 147 palindromic repeats) ranging from 30 bp to 48,383 bp and the 30-39 bp repeats were the majority type. Relative synonymous codon usage (RSCU) analysis uncovered that in total, 9,947 codons were encoding the protein-coding genes (PCGs). Serine (909, 9.1%) and leucine (879, 8.8%) were the two most abundant amino acids, while terminator (32, .3%) was the least abundant amino acid. Ka/Ks analysis indicated that almost all genes were subject to purification selection, except ccmB. Analysis of Lamiaceae mitogenomes constitution revealed that atpB and atpE were unique to the Rotheca serrata and Salvia miltiorrhiza mitogenomes. mttB gene loss was unique to the Boea hygrometrica mitogenome. The core fragments of the Lamiaceae mitogenomes harbored a higher GC content than the specific and variable fragments. In addition, phylogenetic analysis revealed that MCB was closely related to Salvia miltiorrhiza based on the mitogenomes. The current study provided valuable genomic resources for understanding and utilizing this important medicinal plant in the future.

Codon Usage Provide Insights into the Adaptation of Rice Genes under Stress Condition

Plants experience different stresses, i.e., abiotic, or biotic, and to combat them, plants re-program the expression of growth-, metabolism-, and resistance-related genes. These genes differ in their synonymous codon usage frequency and show codon usage bias. Here, we investigated the correlation among codon usage bias, gene expression, and underlying mechanisms in rice under abiotic and biotic stress conditions. The results indicated that genes with higher expression (up- or downregulated) levels had high GC content (≥60%), a low effective number of codon usage (≤40), and exhibited strong biases towards the codons with C/G at the third nucleotide position, irrespective of stress received. TTC, ATC, and CTC were the most preferred codons, while TAC, CAC, AAC, GAC, and TGC were moderately preferred under any stress (abiotic or biotic) condition. Additionally, downregulated genes are under mutational pressure (R² ≥ 0.5) while upregulated genes are under natural selection pressure (R² ≤ 0.5). Based on these results, we also identified the possible target codons that can be used to design an optimized set of genes with specific codons to develop climate-resilient varieties. Conclusively, under stress, rice has a bias towards codon usage which is correlated with GC content, gene expression level, and gene length.

The complete chloroplast genome of Elsholtzia fruticosa (D. Don) Rehd. (Labiatae), an ornamental plant with high medicinal value

Elsholtzia fruticosa is an ornamental plant with high medicinal value. In this study, we sequenced and analyzed the complete chloroplast (cp) genome of the species. The complete cp sequence is 151,550 bp, including the large single-copy (LSC) region of 82,778 bp, the small single-copy (SSC) region of 17,492 bp, and a pair of invert repeats (IRs) regions of 25,640 bp. It encodes 132 unique genes in total, including 87 protein-coding genes, 37 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). The comparative analysis of complete cp genomes showed that the genomic structure and gene order of E. fruticosa cps were conserved. The sequences of rps15, rps19, ycf1, ycf3, ycf15, psbL, psaI, trnG-UCC, trnS-GCU, trnR-UCU, trnL-UAG, trnP-UG, and trnL-UAA serve as hotspots for developing the DNA barcoding of Elsholtzia species. There are 49 SSR loci in the cp genome of E. fruticosa, among which the repeat numbers of mononucleotide, dinucleotide, trinucleotide, tetranucleotide, and pentanucleotides SSR are 37, 9, 3, 0, and 0, respectively. A total of 50 repeats were detected, including 15 forward repeats, seven reverse repeats, 26 palindromic repeats, and two complementary repeats. Phylogenetic analysis based on the complete cp genome and protein-coding DNA sequences of 26 plants indicates that E. fruticosa has a dose relationship with E. splendens and E. byeonsanensis.

A comparative study of the chloroplast genomes of five Lepidium species with high medicinal value

Plantgenomics is a rapidly developing field in medicinal plant research. This study analysed the relevant information of chloroplasts genome sequences of five medicinal plants from the genus Lepidium . We sequenced the complete chloroplast (cp) genomes of Lepidium apetalum Willd. and Lepidium perfoliatum Linnaeus., and assessed their genetic profiles against the reported profiles of Lepidium sativum Linnaeus., Lepidium meyenii Walp., and Lepidium virginicum Linn. We found that L. apetalum and L. perfoliatum possessed 130 distinct genes that included 85 protein-coding, 37 transfer RNA (tRNA), and eight ribosomal RNA (rRNA) genes. Our repeat analyses revealed that L. apetalum harboured 20 direct repeats, 16 palindrome repeats, 30 tandem repeats, and 87 simple sequence repeats, whereas, L. perfoliatum had 15 direct repeats, 20 palindrome repeats, four reverse repeats, 21 tandem repeats, and 98 simple sequence repeats. Using syntenic analysis, we also revealed a high degree of sequence similarity within the coding regions of Lepidium medicinal plant cp genomes, and a high degree of divergence among the intergenic spacers. Pairwise alignment and single-nucleotide polymorphism (SNP) examinations further revealed certain Lepidium -specific gene fragments. Codon usage analysis showed that codon 14 was the most frequently used codon in the Lepidium coding sequences. Further, correlation investigations suggest that L. apetalum and L. perfoliatum originate from similar genetic backgrounds. Analysis of codon usage bias of Lepidium cp genome was strongly influenced by mutation and natural selection. We showed that L. apetalum and L. perfoliatum will likely enhance breeding, species recognition, phylogenetic evolution, and cp genetic engineering of the Lepidium medicinal plants.

Complete Chloroplast Genome Sequences of Four Species in the Caladium Genus: Comparative and Phylogenetic Analyses

Caladiums are promising colorful foliage plants due to their dazzling colors of the leaves, veins, stripes, and patches, which are often cultivated in pots or gardens as decorations. Four wild species, including C. bicolor, C. humboldtii, C. praetermissum, and C. lindenii, were employed in this study, where their chloroplast (cp) genomes were sequenced, assembled, and annotated via high-throughput sequencing. The whole cp genome size ranged from 162,776 bp to 168,888 bp, and the GC contents ranged from 35.09% to 35.91%. Compared with the single large copy (LSC) and single small copy (SSC) regions, more conserved sequences were identified in the inverted repeat regions (IR). We further analyzed the different region borders of nine species of Araceae and found the expansion or contraction of IR/SSC regions might account for the cp genome size variation. Totally, 131 genes were annotated in the cp genomes, including 86 protein-coding genes (PCGs), 37 tRNAs, and eight rRNAs. The effective number of codons (ENC) values and neutrality plot analyses provided the foundation that the natural selection pressure could greatly affect the codon preference. The GC₃ content was significantly lower than that of GC₁ and GC₂, and codons ending with A/U had higher usage preferences. Finally, we conducted phylogenetic relationship analysis based on the chloroplast genomes of twelve species of Araceae, in which C. bicolor and C. humboldtii were grouped together, and C. lindenii was furthest from the other three Caladium species occupying a separate branch. These results will provide a basis for the identification, development, and utilization of Caladium germplasm.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Plastome evolution of Aeonium and Monanthes (Crassulaceae): insights into the variation of plastomic tRNAs, and the patterns of codon usage and aversion

This study reported 13 new plastomes from Aeonium and Monanthes, and observed new markers for phylogeny and DNA barcoding, such as novel tRNA structures and codon usage bias and aversion. The Macaronesian clade of Crassulaceae consists of three genera: Aichryson, with about 15 species; Monanthes, with about 10 species; Aeonium, with about 40 species. Within this clade, Aeonium, known as "the botanical equivalent of Darwin's finches", is regarded as an excellent model plant for researching adaptive evolution. Differing from the well-resolved relationships among three genera of the Macaronesian clade, the internal branching patterns within the genus Aeonium are largely unclear. In this study, we first reported 13 new plastomes from genus Aeonium and the closely related genus Monanthes. We further performed comprehensive analyses of the plastomes, with focuses on the secondary structures of pttRNAs and the patterns of codon usage and aversion. With a typical circular and quadripartite structure, the 13 plastomes ranged from 149,900 to 151,030 bp in size, and the unique pattern in IR junctions might become a family-specific marker for Crassulaceae species. Surprisingly, the π values of plastomes from Monanthes were almost twice those from Aeonium. Most importantly, we strongly recommend that highly polymorphic regions, novel putative pttRNA structures, patterns of codon usage bias and aversion derived from plastomes might have phylogenetic implications, and could act as new markers for DNA barcoding of plants. The results of phylogenetic analyses strongly supported a clear internal branching pattern in Macaronesian clade (represented by Aeonium and Monanthes), with higher nodal support values. The findings reported here will provide new insights into the variation of pttRNAs, and the patterns of codon usage and aversion of the family Crassulaceae.

Comparative Analysis of the Chloroplast Genome for Aconitum Species: Genome Structure and Phylogenetic Relationships

Aconitum is an important medicinal group of the Ranunculaceae family and has been used as conventional medicine in Bai, Yi, and other ethnic groups of China. There are about 350 Aconitum species globally and about 170 species in China. It is challenging to identify the species in morphology, and the lack of molecular biology information hinders the identification and rational utilization of the germplasm of this genus. Therefore, it is necessary to increase the molecular data of Aconitum species. This paper acquired the complete chloroplast (CP) genome sequence of ten medicinal plants of Aconitum species from Yunnan by Illumina paired-end (PE) sequencing technology and compared it with other species in the same family and genus. These CP genomes exhibited typical circular quadripartite structure, and their sizes ranged from 155,475 (A. stylosum) to 155,921 bp (A. vilmoinianum), including a large single-copy region (LSC), a small single-copy region (SSC), and two inverted repeat regions (IRs). Their gene content, order, and GC content (38.1%) were similar. Moreover, their number of genes ranged from 129 (A. vilmoinianum) to 132 (A. ramulosum), including 83 to 85 protein-coding genes (PCGs), 37 tRNA genes (tRNAs), eight rRNA genes (rRNAs), and two pseudogenes. In addition, we performed repeated sequence analysis, genomic structure, and comparative analysis using 42 Aconitum chloroplast genomes, including ten Aconitum chloroplast genomes and other sequenced Aconitum species. A total of 48-79 simple sequence repeats (SSRs) and 17 to 77 long repeat sequences were identified. IR regions showed higher variability than the SSC region and LSC region. Seven mutational hotspots were screened out, including trnK-UUU-trnQ-UGG, psbD, ndhJ-ndhK, clpP, psbH-petB, ycf1, and trnA-UGC-trnI-GAU, respectively. The phylogenetic trees of ten Aconitum species and other Aconitum species revealed that the complete CP genome was beneficial in determining the complex phylogenetic relationships among Aconitum species. This study provides a potential molecular marker and genomic resource for phylogeny and species identification of Aconitum species and an important reference and basis for Ranunculaceae species identification and phylogeny.

Physio-Morphological, Biochemical and Transcriptomic Analyses Provide Insights Into Drought Stress Responses in Mesona chinensis Benth

Drought stress affects the normal growth and development of Mesona chinensis Benth (MCB), which is an important medicinal and edible plant in China. To investigate the physiological and molecular mechanisms of drought resistance in MCB, different concentrations of polyethylene glycol 6000 (PEG6000) (0, 5, 10, and 15%) were used to simulate drought conditions in this study. Results showed that the growth of MCB was significantly limited under drought stress conditions. Drought stress induced the increases in the contents of Chla, Chlb, Chla + b, soluble protein, soluble sugar, and soluble pectin and the activities of superoxide dismutase (SOD), catalase (CAT), total antioxidant capacity (TAC), hydrogen peroxide (H2O2), and malondialdehyde (MDA). Transcriptome analysis revealed 3,494 differentially expressed genes (DEGs) (1,961 up-regulated and 1,533 down-regulated) between the control and 15% PEG6000 treatments. These DEGs were identified to be involved in the 10 metabolic pathways, including "plant hormone signal transduction," "brassinosteroid biosynthesis," "plant-pathogen interaction," "MAPK signaling pathway-plant," "starch and sucrose metabolism," "pentose and glucuronate interconversions," "phenylpropanoid biosynthesis," "galactose metabolism," "monoterpenoid biosynthesis," and "ribosome." In addition, transcription factors (TFs) analysis showed 8 out of 204 TFs, TRINITY_DN3232_c0_g1 [ABA-responsive element (ABRE)-binding transcription factor1, AREB1], TRINITY_DN4161_c0_g1 (auxin response factor, ARF), TRINITY_DN3183_c0_g2 (abscisic acid-insensitive 5-like protein, ABI5), TRINITY_DN28414_c0_g2 (ethylene-responsive transcription factor ERF1b, ERF1b), TRINITY_DN9557_c0_g1 (phytochrome-interacting factor, PIF3), TRINITY_DN11435_c1_g1, TRINITY_DN2608_c0_g1, and TRINITY_DN6742_c0_g1, were closely related to the "plant hormone signal transduction" pathway. Taken together, it was inferred that these pathways and TFs might play important roles in response to drought stress in MCB. The current study provided important information for MCB drought resistance breeding in the future.

Analysis of codon usage bias of chloroplast genomes in Gynostemma species

Gynostemma plants are important Chinese medicinal material and economic crops. Codon usage analysis is a good way to understand organism evolution and phylogeny. There is no report yet about analysis of codon usage bias of chloroplast genomes in Gynostemma species. In this study, the chloroplast genomes in nine Gynostemma species were analyzed systematically to explore the factors affecting the formation of codon usage bias. The codon usage indicators were analyzed. Multivariate statistical analysis including analysis of neutrality plot, effective number of codons plot, parity rule 2 plot and correspondence were performed. Composition analysis of codons showed that the frequency of GC in chloroplast genes of all nine Gynostemma species was less than 50%, and the protein-coding sequences of chloroplast genes preferred to end with A/T at the third codon position. The chloroplast genes had an overall weak codon usage bias. A total of 29 high frequency codons and 12 optimal codons were identified. These could provide useful information in optimizing and modifying codons thus improving the gene expression of Gynostemma species. The results of multivariate analysis showed that the codon usage patterns were not only affected by single one factor but multiple factors. Mutation pressure, natural selection and base composition might have an influence on the codon usage patterns while natural selection might be the main determinant. The study could provide a reference for organism evolution and phylogeny of Gynostemma species and help to understand the patterns of codons in chloroplast genomes in other plant species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01105-z.

Identification of Differentially Expressed Genes and Pathways Involved in Growth and Development of Mesona chinensis Benth Under Red- and Blue-Light Conditions

Mesona chinensis Benth (MCB) is an important Chinese herbal medicine. The plant factories might be one of the ways to solve the shortage of MCB supply. In this study, the MCB seedlings were treated under the red (R) and blue (B) lights in the plant factory. Results showed that the red light promoted the growth and development of MCB in comparison with the blue light. Under the red-light condition, the biomass, plant height, and root characteristics were significantly higher than those under blue-light condition, while the soil and plant analyzer development (SPAD) under the red-light treatment was significantly lower than that under the blue-light treatment. Red light also significantly promoted the content of soluble sugar and pectin of MCB compared with blue light. Transcriptome analysis showed that a total of 4,165 differentially expressed genes (DEGs) were detected including 2,034 upregulated and 2,131 downregulated. Of these, 1,112 DEGs including 410 upregulated and 702 downregulated genes were associated with 111 pathways. Moreover, a total of 8,723 differentially expressed transcription factors (TFs) were identified in R vs. B, and these TFs were distributed in 56 gene families. Metabonomic results revealed that a total of 184 metabolites and 99 differentially expressed metabolites (DEMs) (42 upregulated and 57 downregulated) were identified in the red- and blue-light treatments. Integrative analysis of transcriptome and metabolome unveiled that a total of 24 pathways included 70 compounds (metabolites) and were associated with 28 unigenes. In particular, these pathways included starch and sucrose metabolism, phenylpropanoid biosynthesis, cysteine and methionine metabolism, glycolysis/gluconeogenesis, and pentose and glucuronate interconversions. The unigenes included asparagine synthetase (AS), thymidine kinase (TK), alpha, alpha-trehalose-phosphate synthase (TPS), phosphatase IMPL1 (IMPL1), dihydroflavonol 4-reductase (D4R), and 4-coumarate-CoA ligase-like 6 (4CL6), bifunctional aspartokinase-homoserine dehydrogenase 1 (thrA), and abscisic acid 8'-hydroxylase 2 isoform X1 (ABA8). It was indicated that these pathways and genes might play important roles in the growth and development of MCB. This study laid a foundation for the future research of MCB.