Next Article in Journal
Effects of Age and Lifelong Moderate-Intensity Exercise Training on Rats’ Testicular Function
Previous Article in Journal
TRIM26 Maintains Cell Survival in Response to Oxidative Stress through Regulating DNA Glycosylase Stability
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 19
10.3390/ijms231911617
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization and Comparative Analysis of Chloroplast Genomes in Five Uncaria Species Endemic to China

by
Min-Min Chen
1,
Miao Zhang
1,
Zong-Suo Liang
2,* and
Qiu-Ling He
1,*
1
Key Laboratory of Plant Secondary Metabolism and Regulation of Zhejiang Province, College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
2
Shaoxing Academy of Biomedicine, Zhejiang Sci-Tech University, Shaoxing 312366, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(19), 11617; https://doi.org/10.3390/ijms231911617
Submission received: 15 August 2022 / Revised: 22 September 2022 / Accepted: 27 September 2022 / Published: 1 October 2022
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Uncaria, a perennial vine from the Rubiaceae family, is a typical Chinese traditional medicine. Currently, uncertainty exists over the Uncaria genus’ evolutionary relationships and germplasm identification. The complete chloroplast genomes of four Uncaria species mentioned in the Chinese Pharmacopoeia and Uncaria scandens (an easily confused counterfeit) were sequenced and annotated. The findings demonstrated that the whole chloroplast genome of Uncaria genus is 153,780–155,138 bp in full length, encoding a total of 128–131 genes, containing 83–86 protein-coding genes, eight rRNAs and 37 tRNAs. These regions, which include eleven highly variable loci and 31–49 SSRs, can be used to create significant molecular markers for the Uncaria genus. The phylogenetic tree was constructed according to protein-coding genes and the whole chloroplast genome sequences of five Uncaria species using four methods. The topology of the two phylogenetic trees showed no difference. The sequences of U. rhynchophylla and U. scandens are clustered in one group, while the U. hirsuta and U. macrophylla are clustered in another group. U. sessilifructus is clustered together with the above two small clades. New insights on the relationship were revealed via phylogenetic research in five Uncaria species. This study will provide a theoretical basis for identifying U. rhynchophylla and its counterfeits, as well as the species of the Uncaria genus. This research provides the initial chloroplast genome report of Uncaria, contributes to elucidating the chloroplast genome evolution of Uncaria in China.

1. Introduction

The Rubiaceae is the fifth largest flowering plant family spread worldwide. Uncaria, a perennial vine from the Rubiaceae family, often climbs on top of other plants with its hook-like inflorescence stalk. It is indigenous to tropical areas of Asia, South America, and Africa and is mainly used for medicinal purposes and as a natural dye. The Pharmacopoeia of the People’s Republic of China contains five species, including Uncaria sessilifructus, Uncaria sinensis, Uncaria hirsuta, Uncaria macrophylla, and Uncaria rhynchophylla. In China, they are mainly distributed in Yunnan, Guangdong, Guangxi, and southeast Guizhou Provinces [1]. Wild resources often grow in mountain forests or hilly areas. There are differences in agronomic characters, morphological characteristics, effective components, and molecular characteristics among different germplasms.
The Uncaria plant contains various of chemical constituents, including alkaloids, flavonoids, triterpenes, aliphatic compounds, phenolic acids, coumarins, lignans, steroids, and anthraquinones [2,3,4,5,6,7]. Alkaloids are the main medicinal components in Uncaria spp. Their main types are monoterpene indole alkaloids and β-carboline. Previous studies have shown that rhynchophylline can regulate the expression of phosphoinositide-3-kinase (PI3Ks), protein kinase B and glycogen synthase kinase 3β, which are involved in regulating neuronal signaling pathways [8]. The racemic analogues of rhynchophylline (G2) and their stereoisomers can effectively treat microvascular dysfunction caused by diabetes. In addition, it is reported that isorhynchophylline can cause an effect by regulating multiple signal cascades in hepatocellular carcinoma cells [9].
Chloroplasts are semiautonomous organelles commonly found in higher plants, algae cells, and some protozoa. For photosynthesis and energy transformation in plants, chloroplasts are crucial structures. Plant chloroplast genes are usually composed of photosynthesis, transcriptional translation and expression, other protein-coding genes, and open reading frame (ORF) [10]. The chloroplast genome usually presents a four-segment structure, which is composed of two reverse repeat regions (IR), a large single copy region (LSC), and a small single copy region (SSC) [11]. Two IR regions separate the LSC region from the SSC. Compared to mitochondrial and nuclear genomes, chloroplast genomes have uniparental inheritance, a low rate of nucleotide substitution, and a simple structural design [12,13]. The matrilineal inheritance of the chloroplast genome is responsible for maintaining a relatively stable genetic structure with a rate of evolution intermediate to that of the nuclear and mitochondrial genomes. Therefore, the chloroplast genomes become an ideal resource and a versatile tool for phylogenetic studies at different levels for species identification [14]. However, the chloroplast genomes of Uncaria spp. had not been reported, leading to limitations in the mining of its genetic information and the development of phylogenetic studies.
Therefore, in this study, morphological characteristics and the chloroplast genome of four species of Uncaria spp. were first reported. A comparative genomic analysis among five species of Uncaria was performed to reconstruct the phylogeny and explore the high variability sites to provide pertinent details for potential molecular markers and lay the theoretical foundation for the identification of genuine and counterfeit herbs of Uncaria.

2. Results

2.1. Morphological Comparison

The five species of Uncaria used in this study are as follows: Uncaria rhynchophylla, Uncaria macrophylla, Uncaria hirsuta, Uncaria scandens, and Uncaria sessilifructus. The morphological study showed that the stipules of U. rhynchophylla were significantly different from the other four species of Uncaria plant, which were bifid linear stipules. The stipules of U. macrophylla, U. hirsuta, and U. sessilifructus are all deeply bifid triangular lobes, while the U. sessilifructus’ are slenderer. Furthermore, the size of stipules in U. scandens are larger than other species, with which broad strips or triangular lobes. The leaves of the five species of Uncaria are all simple and opposite leaves with cylindrical to square stems. The stem epidermal color of U. macrophylla and U. scandens showed yellowish green while U. rhynchophylla, U. sessilifructus, and U. hirsuta showed lavender to brown, yellow, and light gray to yellowish green colors, respectively. In addition, the leaf characteristics of these five species of Uncaria also show marked differences. The leaves of U. sessilifructus are leathery, smooth, and brittle, ovate to long elliptic, with gray abaxially. The leaves of U. scandens are subpapery and pilose on both surfaces. The leaves of U. hirsuta are ovate-lanceolate, almost glabrous on the leaf surface, with a sparse hirsute abaxial side. The leaves of U. macrophylla are subleathery, brittle, with hairs only on the veins, and the leaves color is yellowish brown with hirsute abaxially. The leaves of U. rhynchophylla are smooth and subleathery, and the surface of young leaves is reddish brown (Figure 1 and Figure S1, Supplementary Table S1).

2.2. Chloroplast Genome Sequencing and Assembly

Illumina paired-end (150 bp) sequencing was used to obtain the raw reads of U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens. All raw reads were processed to eliminate adaptor sequences, short reads, and low-quality bases. Then, 13,659,533–18,093,602 bp clean reads were yielded. Contigs were produced by the de novo assembly, which were subsequently joined to generate the final cp genome for four species. The lengths of chloroplast genome sequences were 154,605, 155,177, 155,138, 155,138, and 153,780 bp for U. rhynchophylla (MN723865), U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens, respectively (Table 1, Figure 2). Four Uncaria species of chloroplast genome sequences were uploaded to NCBI GenBank.

2.3. Chloroplast Genome Structure and Characteristics Analyses

The entire cp genomes showed the standard quadripartite structure found in most angiosperms, which includes the large single copy (LSC), the small single copy (SSC), and two inverted repeats (Ira and Irb) [15]. The length varied from 85,146 bp to 85,749 bp in the LSC region, from 16,989 bp to 18,145 bp in the SSC region, and from 25,657 bp to 25,690 bp in IR region. The percentages of the cp genome encoding protein, rRNAs, and tRNAs are, respectively, 49.88–52.08, 5.83–5.88, and 1.77–1.80, and the remaining 40.27–42.47% is composed of non-coding regions. In all five species, the overall GC content was 37.5–37.7%. The GC contents of the LSC, SSC, and IR regions were 35.4–35.6%, 31.5–31.7%, and 43.2%. The GC contents of IR region were significantly higher than LSC and SSC. It is highest in the IR regions, lowest in the SSC regions and moderate in the LSC regions. The GC contents of the protein-coding genes, rRNA and tRNA were 37.9–38.1%, 55.3–55.4%, and 53.3–53.5%. Meanwhile, the genomes of five Uncaria species contained 128–131 genes in total, including 83-86 protein-coding genes, 37 rRNA genes, and eight tRNA genes (Table 2).
In each Uncaria species, there were sixteen genes with two copies, which were comprised of five protein-coding genes (ndhB, rps7, ycf2, rpl2, rpl23), seven tRNA genes (trnI-GAU, trnA-UGC, trnL-CAA, trnI-CAU, trnR-ACG, trnV-GAC, trnN-GUU), and four rRNA genes (rrn16, rrn23, rrn4.5, rrn5). The three genes (clpP, rps12 and ycf3) had two introns, while nine protein-coding genes (atpF, rpl2, ndhB, ndhA, rps16, rpoC1, petD, petB, rpl16) and seven tRNA genes (trnI-GAU, trnG-UCC, trnV-UAC, trnK-UUU, trnL-UAA, trnL-UAG, trnA-UGC) only had one intron (Supplementary Table S2).

2.4. Codon Usage Analyses

We evaluated the codon usage of the protein-coding sequences in five species (Supplementary Table S3). “Count” represents the number of codons in one sample (used universal genetic code) [16]. There were 25,728 (U. rhynchophylla), 26,722 (U. sessilifructus), 26,578 (U. macrophylla), 26,562 (U. hirsuta), and 25,699 (U. scandens) counts, respectively. Except for the three stop codons, with an average of 1048.6 (varied from 1014 to 1069), the UUU codons for the amino acid phenylalanine (Phe) were the most numerous, while the CGC codon corresponding to arginine (Arg) was the smallest with an average of 116.6 (varied from 112 to 121). With the exception of the unique codons (Trp, Met), Leucine (Leu) was the amino acid with the most codons, with an average number of 2585, while cysteine (Cys) was the least (640.6) (Figure 3).
In addition, relative synonymous codon usage (RSCU) values are the relationship between the number of practical codon emergence and the number of anticipated codon emergence [17]. If the RSCU > 1, which means that this codon has the higher preference. The results showed that compared with A/T (U), the codons ending with C/G had fewer codons and RSCU < 1. In all protein-coding sequences of five species, 30 (46.875%) codons had RSCU > 1 (higher preference). Most of these codons had an ending of U or A. AGA codon corresponding to the amino acid Arginine (Arg) was the highest preferred across the studied five species with an average RSCU of 1.906. Methionine (Met) and threonine (Thr) had no bias (RSCU = 1). However, 32 (50.00%) codons had RSCU < 1 and low preference (Figure 4).

2.5. Repeat Sequence Analysis

Repeat sequences occupy an important position in phylogenetic studies and genome rearrangements [18]. Repeat sequences are indispensable to make indels and substitutions [19]. Dispersed repeats were identified for the five cp genomes using REPuter, including forward repeats (F), reverse repeats (R), palindromic repeats (P), and complement repeats (C). The five Uncaria cp genomes contained 49–51 dispersed repeats, comprising 17–23 forward repeats, 14–19 palindromic repeats, 11–15 reverse repeats, and one-four complement repeat. Five cp genomes possessed four types of dispersed repeated sequences. Forward repeat was the most universal type (Figure 5A). Meanwhile, it largely consisted of dispersed repeats with a unit length of <20 bp and 20–30 bp in all five cp genomes (Figure 5B).
Tandem repeats finder was used to identify the tandem repeated sequences for the five cp genomes. The number of the tandem repeated sequences were different in five cp genomes, as follows: 28, 32, 30, 29, 20. It ranged from 20 to 32. U. sessilifructus had the most tandem repeats and U. scandens had the least. The period size (repeat unit size) for the tandem repeats ranged from 3 to 49 bp. Tandem repeated sequences with periods ranging from 11 to 30 bp predominated in the cp genomes (Figure 6A). The majority of the tandem repeats in the five cp genomes had a copy number of no more than 4 (Figure 6B). Five species’ sequence lengths of tandem ranged from 25 to 98 bp. Tandem repeats with sequence lengths ranging from 30 to 50 bp predominated in all cp genomes (Figure 6C).

2.6. SSR Analysis

The five cp genomes’ SSRs were discovered using MISA. There were 49, 40, 31, 31, and 37 SSRs in U. rhynchophylla, U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens, respectively. Most SSRs possess a mononucleotide repeat motif in five cp genomes, ranging from 26 to 36. Dinucleotides SSRs are the second most common. The five cp genomes did not have hexanucleotide SSRs. Only U. scandens did not find pentanucleotides SSRs. Meanwhile, no trinucleotides SSRs and tetranucleotides SSRs were discovered in U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens (Figure 7A).
Five cp genomes yielded a total of 188 SSRs, of which 121 (64.36%) were in IGS, 27 (14.36%) were in CDS, and 40 (21.28%) were in intron. Most loci were located in intergenic spacer (IGS) regions, and the fewest loci were located in CDS, probably due to the IGS regions having higher mutation rates than coding regions. In LSC, IR, and SSC regions, 163 (86.70%), 9 (4.79%), and 16 (8.51%) SSRs were found, respectively. Most SSRs existed in the LSC region (Figure 7BC). In five cp genomes, the mononucleotide A/T repeat units had the highest percentage occupancy (73.47–83.87%). Furthermore, the largest number of mononucleotides SSRs (A/T) was found in U. rhynchophylla (13A + 23T), while U. macrophylla and U. hirsuta had the least mononucleotide SSRs (8A + 18T, 7A + 19T) (Figure 7D).

2.7. Contraction and Expansion of Inverted Repeats

The IR region is well-conserved, but the IR/SC boundary regions differ among plant species. The IR/SC border positions and their adjacent genes were compared using IRscope (Figure 8). In five cp genomes, the gene rps19 crossed the LSC/IRb and was primarily found in the LSC region with a length of 247 bp and the remaining 32 bp in the IRb region. These overlapping sequences resulted in pseudogene sequences of rps19 gene at the IRa/SSC boundary. The ycf1 gene traversed the SSC and IRa regions, which was located at the SSC region with 4489–4504 bp (U. sessilifructus: 4504 bp, other: 4489 bp), while 1136 bp sequences located at IRa. These overlapping sequences resulted in the pseudogene sequences of ycf1 gene at the IRb/SSC boundary. The ndhF genes did not cross the IRb/SSC boundary. JLA was located between trnH and rpl2, the distances between trnH and JLA were 9–12 bp. The results indicated that there were differences at the boundaries in the five cp genomes. Overall, the length variance of the complete genome sequences among the five cp genomes was caused by the contraction and extension of IR/SC boundary regions. In addition, the results of the synteny analysis showed that the chloroplast genomes did not undergo genomic rearrangements (Figure 9).

2.8. RNA Editing Site Analysis and Selective Pressure Analyses

The five species were predicted to have identical RNA editing sites. A total of 123 RNA editing sites were predicted in 28 genes (Supplementary Tables S4 and S5). All predicted editing sites were C to U transitions. The majority of predicted editing sites were found in the ndhB gene (14), which is consistent with the results for other plants, such as Rhynchanthus beesianus, Pommereschea lackneri, and Hedychium coronarium [20]. In addition, psbB contained ten predicted editing sites; atpB and ndhD contained nine predicted editing sites; ndhF contained eight predicted editing sites; and atpA, psaB, and rpoB contained seven predicted editing sites. From one to five predicted editing sites were present in the other gene (rpoC2, 5; ndhA, atpI, ccsA, matK, ndhG, rpoA, 4; petB, petG, 3; accD, atpF, petD, rpoC1, rps14, rpl2, 2; psbE, psbF, rpl23, rps2, ycf3, 1). Of the total number of RNA editing sites, there were 49 (39.84%) RNA editing sites at the first base of the codon and 74 (60.16%) RNA editing sites at the second base of the codon. The transformation of amino acids included 53 hydrophilic amino acids to hydrophobic amino acids (S→L, S→F, H→Y, T→M, R→W, T→I), 15 hydrophobic amino acids to hydrophilic amino acids (P→S), and 55 hydrophobic amino acids to different hydrophobic amino acids (P→L, P→F, L→F, A→V).
We calculated the ratio of dN/dS to discover proteins undergoing selection in the protein-coding genes. The ratio of ω = dN/dS has become a measure of selection pressure, with ω = 1, >1, and <1 indicating neutral evolution, positive selection, and negative or purifying selection, respectively [21]. The dN/dS were between 0.0001–0.81066 in five Uncaria species chloroplast genomes, the dN/dS values are all less than 1, suggesting that almost all of the genes were undergoing strong purifying selection. The dN/dS values of ten genes, including photosynthetic device genes and ribosomal protein genes, such as ndhF, ndhJ, rbcL, rpl2, rpoC2, rps11, rps15, rps8, ycf1, and ycf2, were >0.4, demonstrating that they are evolving particularly fast at the protein level. It can be inferred that there may be the potential for positive selection loci (Figure 10) [22]. Furthermore, most of the dN/dS ratio values in the protein-coding genes of U. hirsuta vs. U. rhynchophylla and U. macrophylla vs. U. rhynchophylla were less than 1, except petA and petB, whose values were 1.207 and 1.206, respectively, indicating that both genes were undergoing positive selection. In the protein-coding genes of U. macrophylla vs. U. rhynchophylla and U. hirsuta vs. U. rhynchophylla, most genes had the lowest dN/dS values. Meanwhile, we found seven genes (ndhD, rbcL, petD, rpl20, rps3, ycf1, ycf2) (dN/dS > 0.4) had the higher ratios of dN and dS in the U. macrophylla vs. U. rhynchophylla set, while seven genes (petD, rbcL, rpl2, rpl20, rps3, ycf1, ycf2) were also identified with higher dN/dS values in the U. macrophylla vs. U. rhynchophylla set (Figure 11).

2.9. Sequence Divergence Analysis

Using Uncaria rhynchophylla (MN723865) as a reference, we used the mVISTA tool to compare the variations in five cp genomes and used DnaSP to identify the divergence hotspot regions (Figure 12 and Figure 13). The results revealed that the five species were proven to have high sequence identity by mVISTA. The pi value of nucleotide diversity ranged from 0 to 0.028. We found that nine intergenic spacers (IGSs) had relatively higher divergence values (Pi > 0.015) [23], including trnH-GUG-psbA, rps16-trnQ-UUG, trnQ-UUG-psbK, atpH-atpL, ndhC-trnV-UAC, petA-psbJ, rpl32-trnL-UAG, ndhF-rpl32, and rps15-ycf1. The partial sequences of the ycf1 and ccsA genes showed a comparatively high nucleotide diversity. Most of the higher divergence regions were found in LSC and SSC regions, while IR regions showed a lower sequence divergence in the cp genomes. The results from mVISTA were consistent with those from DnaSP. The higher divergence regions can be used to identify the divergence between closely related species.

2.10. Phylogenetic Analysis

The cp genomes and protein-coding genes have been successfully used to reveal phylogenetic relationships [24,25]. The available chloroplast genomes sequences of 13 species from Rubiaceae and Asclepias nivea were selected to accomplish the phylogenetic analysis, including five Uncaria species, Neolamarckia cadamba, Mitragyna speciosa, Cinchona officinalis, Galium aparine, Gardenia jasminoides, Paederia scandens, Coffea canephora, and Coffea arabica. Asclepias nivea was set as outgroups. The bayesian inference (BI), maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ) methods yielded nearly identical the tree topologies. Consequently, phylogenetic trees using four methods were integrated. From left to right, numbers above the branches are, respectively, represented Bayesian inference (BI) posterior probabilities (PPs) values, bootstrap support values of maximum likelihood (ML), neighbor-joining (NJ), and maximum parsimony (MP) from the complete cp genomes (Figure 14). The cp genome and protein-coding genes both produced phylogenetic trees with identical topological patterns. The phylogenetic trees all achieved relatively higher bootstrap support values and high BI (PPBI = 1.00), nearly all bootstrap values reached 100%. A strong monophyletic branch of five Uncaria species was formed, and they are clustered into two clades with 100% bootstrap values and one PPs values. One monophyletic clade was included two sub-clades with 97–100% bootstrap values and one PPs values: U. scandens, U. macrophylla, and U. rhynchophylla, as well as U. hirsuta. U. scandens shares a clade with U. rhynchophylla and U. macrophylla shares a clade with U. hirsuta, which indicated that U. scandens was closely related to U. rhynchophylla. It is suggested that U. hirsuta and U. macrophylla were closely related. The other monophyletic clade included U. sessilifructus and four varieties of Uncaria, and the phylogenetic tree revealed that U. sessilifructus was closely related to four other varieties of Uncaria in this study.

3. Discussion

Currently, only a few numbers of the whole chloroplast genome sequences for the Uncaria genus were released into GenBank. Thus, the complete cp genomes of five Uncaria species were reported in this study. The cp genomes of five Uncaria species comprised 128–131 genes, including 83-86 protein-coding genes, 37 rRNA genes, and eight tRNA genes, other Rubiaceae plants have comparable gene patterns as well [26]. Chloroplast genomes are largely preserved in terms of gene composition and structure, compared to the nuclear and mitochondrial genomes. Additionally, recent research has demonstrated that gene or intron loss occurs in chloroplast genomes [27,28,29]. In addition, the deletion of the petD intron was observed in U. macrophylla and U. hirsuta. petD intron loss had been reported in numerous angiosperms, such as a few species of the families Quercus and Brassicaceae [30,31]. Introns play a key role in gene expression regulation. They can enhance gene expression in a given location and at a specific time [32]. Some introns have been shown to improve or necessitate proper amounts of mRNA transcription, processing, and transport. Intron loss appears to be under selection pressure in several unicellular eukaryotes [33]. The loss of the ndhD gene was only found in U. hirsuta, consistent with the cp genome of Vigna angularis, Lathyrus sativus, Trifolium subterraneum, and Robinia pseudoacacia [34]. Natural selection or chromosomal rearrangements during evolution may account for their absence in the species [35]. Pseudogenization and gene deletions were common and lineage-specific, with the ndh genes involved in electron recycling being particularly often and independently lost [29]. It is suggested that the ndh pathway is redundant in some species of Uncaria spp., and it may even be selected against [36]. Further investigation into the nature of cp genome modification should yield a new understanding of which genetic (and thus physiological) functions are lost. Thus, we infer that the deletion of some genes encoded in chloroplasts might not have an impact on the plants’ life cycle and homologous genes from other sources may be capable of performing the functions of deleted genes, although experiments are still needed to support this hypothesis.
The herbal medicine industry is rapidly expanding in China and the authenticity of herbs are becoming a major concern. U. scandens is a confusing counterfeit of U. rhynchophylla, so we sought to find the differences between these two species at the molecular level. Compared to the other four Uncaria species, two genes (matk and ccsA) are missing in U. scandens. The ccsA gene produces a protein called cytochrome C biogenesis, which mediates the binding of heme to c-type cytochromes [37,38]. Pedicularis cyathophylloides and Pedicularis superba also had ccsA gene deletion [39]. It is suggested that changes in the selection pressure on CCSA genes located at the junction of the IR and SSC regions may cause more frequent parallel loss or reversal of these genes relative to other genes in the LSC and IR regions. According to previous reports, the matK gene is a highly variable sequence frequently employed in many plant taxa phylogenetic studies [40,41]. Additionally, the loss of matK genes also occurs in the Cuscuta genus, Arachnitis genus, and Thismia genus [10,42]. The loss of ccsA and matK happened in other U. scandens individuals, which we had also uploaded on NCBI (GenBank ON086916). The current study will contribute to the identifying of U. rhynchophylla and its counterfeits, providing a specific experimental foundation for future studies.
The chloroplast of five Uncaria species was 153,780–155,177 bp with a quadripartite structure (LSC, SSC, and two IR regions), a characteristic of higher plants [43]. Whether in LSC regions, SSC regions, IR regions, or the complete cp genome, the five cp genomes of the Uncaria displayed a high degree of conservation in terms of GC content. The GC content was identical to the other Rubiaceae family members [26,44,45]. The IR areas contained four rRNA genes with high GC content, which supported the hypothesis of previous research [46,47,48,49]. A previous study found that the GC content was substantially higher, most likely due to the abundance of rRNA in the IR regions [50]. However, the exact specifics are still unknown. GC content is an important property of plastid genomes, which are likely produced following endosymbiosis by DNA replication and repair [51]. In viruses, GC content is not determined by the structure of the gene, but it is connected to its location [52]. More research is needed to determine if or not this exists in the chloroplast genome. It is supposed that one of the factors contributing to the low diversity in tRNA sequences and IR regions may be the high GC content, the maintenance of sequence stability may be greatly influenced by GC content. More evidence is needed from subsequent experiments to provide further.
During the evolution of terrestrial plants, the chloroplast genome’s IR region frequently experiences length variations, which results in the formation of a variety of boundary features [32,49,53]. Furthermore, the primary causes of the size change in cp genomes are the contraction and expansion of IR regions and the single-copy (SC) border regions, which influences the cp genome’s rate of evolution [54,55]. This research proved that boundary genes in the Uncaria species were mainly rpl22, rps19, rpl2, ndhF, trnH, and ycf1. Contraction and expansion as well as variation at the SC-IR boundary were identified, suggesting that gene location information in the SC-IR region can reveal the distance between species to some extent. During the evolution of terrestrial plants, length variation in the IR region of the chloroplast genome was a frequent occurrence, which resulted in the formation of a variety of boundary features [32,49,53]. Furthermore, the contraction and expansion of IR regions and the single-copy (SC) border regions are also the main contributors to the size variation in cp genomes, which influences the evolutionary rate of the cp genome [54,55]. This study demonstrated that boundary genes in the species of the Uncaria were mainly rpl22, rps19, rpl2, ndhF, trnH, and ycf1. Variation at the SC-IR boundary and contraction and expansion were observed in five Uncaria species, which suggested that the position information of the genes in the SC-IR region can, to some extent, reveal the distance between species.
It is a well-known fact that bias in codon usage occurs in wide different varieties of organisms. Codon usage preferences are closely linked to gene expression and have an impact on the level of mRNA and proteins in the genome [56,57,58,59]. The most abundant amino acid in the Uncaria was Leucine (Leu) with an average number of 2585, which had also been regularly observed in the other angiosperms. Meanwhile, in accordance with results for the chloroplast genomes of other angiosperms, our research showed that the majority of codons ending in A/U have RSCU values higher than 1, which may have been caused by a bias in composition toward a high A/T ratio [60,61,62]. Further studies will provide a deeper understanding of Uncaria’s gene expression and molecular evolutionary mechanisms. In the five cp genomes, the largest similarities in codon usage were discovered. These findings may indicate that these five Uncaria species underwent comparable environmental pressure during their evolutionary process. Identifying the RNA editing sites is critical for understanding the proper translation process and gene mutations in the cp genome [63]. In our study, potential RNA-editing sites in five Uncaria species were explored using the PREP tool. Twenty-eight genes were found to include a total of 123 predicted RNA editing sites, with the highest number found in the ndhB genes, which is consistent with the general features of RNA editing of chloroplast genes in higher plants [64].
The plastid genome contains many oligonucleotide repeats, which are assumed to be biomarkers for mutational hotspots [65,66,67,68]. In the present study, we detected forward repeats (F), reverse repeats (R), palindromic repeats (P), and complement repeats (C). Most of the dispersed repeated sequences were less than or equal to 30 bp. Tandem repeated sequences with periods ranging from 11 to 30 bp predominated in the cp genomes. The sequence length of tandem repeats of five Uncaria species ranged from 25 to 98 bp, which were similar to those found in Salvia [16]. In higher plants, SSRs are recognized as significant molecular markers for population variation investigations and are widely used to assess genetic diversity, population structure and evolutionary studies [69,70,71,72]. In addition, due to their distinct uniparental inheritance, SSRs have been widely used as molecular markers [73]. In total, five Uncaria species cp genomes with 31–49 perfect SSRs were found, with strong A/T bias. The majority of SSR types were discovered to be mono-nucleotide repeats in this study, which also discovered that the non-coding region included the largest number of SSRs. These repeats provide a crucial starting point for creating genetic markers in Uncaria species, which may be used in phylogenetic and ecological research. We will be able to assess the variation at the intraspecific level after the SSRs has been found in the cp genome of five Uncaria species.
In the current study, we identified that nine intergenic spacers (IGSs) had relatively high divergence values (Pi > 0.015), including trnH-GUG-psbA, rps16-trnQ-UUG, trnQ-UUG-psbK, atpH-atpL, ndhC-trnV-UAC, petA-psbJ, rpl32-trnL-UAG, ndhF-rpl32, and rps15-ycf1. Moreover, the partial sequences of the ccsA and ycf1 genes showed a relatively high nucleotide diversity. These variable regions could be exploited to provide molecular markers for DNA barcoding and genealogical research in Uncaria. Over the last decade, nearly all taxonomic levels of phylogenetic connections have been successfully resolved using whole chloroplast genomes and protein-coding genes [24,25,74]. Our research provides the first phylogeny of the Uncaria genus based on a significant sampling of both chloroplast genomes and protein-coding genes. The topology structures of phylogenetic trees constructed with the whole chloroplast genome and protein-coding regions were almost identical. Phylogenetic analysis strongly demonstrated that U. scandens had been resolved as a sister relationship with U. rhynchophylla, and the same relationship was determined between U. hirsuta, and U. macrophylla. Furthermore, U. sessilifructus formed a clade with other four species. It is believed that U. scandens and U. rhynchophylla are more closely related because U. scandens is a commonly confused counterfeit of U. rhynchophylla. However, the results of the dN/dS analysis indicated that several different genes had the highest ratios of dN and dS, indicating that they were selected for sequence diversity and affect the sequence evolution rate of the chloroplast genomes of U. hirsuta and U. machophylla. This research will not only help shed lighter on the evolutionary position of Uncaria but also offer useful chloroplast genomic data for future investigations into the genesis and diversification of Rubiaceae. Overall, our phylogenomic investigations using chloroplast genomes have made the first successful effort to clarify intrageneric connections within the Uncaria genus [49].

4. Materials and Methods

4.1. Phenotype Measurement

The morphological characteristics of five Uncaria species were observed according to the taxonomic characteristics of botany and the characteristics of Uncaria spp. in Chinese flora. All the morphologic records based on our field observations.

4.2. Chloroplast DNA Extraction, Sequencing, Assembly and Annotation

To characterize the feature of among endemic species of Uncaria in China, five members of the Uncaria family were collected from Guangxi Province in China in July 2021 and subsequently deposited in the College of Life Sciences and Medicine, Zhejiang Sci-Tech University (Hangzhou, China), with herbarium numbers ZSTU70128 (Uncaria rhynchophylla), ZSTU09554 (Uncaria macrophylla), ZSTU01022 (Uncaria hirsuta), ZSTU09751 (Uncaria sessilifructus), and ZSTU70322 (Uncaria scandens). Following the manufacturer’s directions, DNA Plantzol Reagent (Invitrogen, Carlsbad, CA, USA) was used to extract the total DNA from the plant’s leaves [75,76,77]. Utilizing Illumina’s TruSeq Nano DNA Library Preparation kit (350 bp median insert), sequencing libraries were created in accordance with the manufacturer’s instructions [78,79,80]. The Illumina HiSeq 2500 platform was used to produce the chloroplast sequences (Illumina Inc., San Diego, CA, USA). Trimmomatic v0.39 software (Jülich, Germany) was used to treat all raw readings in order to eliminate adaptor sequences, short reads (length 75 bp), and low-quality bases (Q-value 20). Following that, adaptors were cut to produce high-quality clean reads (150 bp PE read length). To complete the chloroplast genome, these clear data were assembled utilizing GetOrganelle (-R 15-k 105,121) (Kunming, China) [81]. The sequence of Uncaria rhynchophylla (MN723865) was used as seeds for reference. Plastid-homologous sequences were selected from the total cellular DNA reads, and then assembled with a de novo assembly pipeline to obtain many longer fragment sequences (Contigs). The Contigs are then positioned on the chloroplast reference genomes of closely related species to obtain longer chloroplast genome fragment sequences. Finally, the full-length framework of the chloroplast genome sequence was constructed. The chloroplast genomes annotation was assembled based on the comparison by CPGAVAS2 (Beijing, China) [82] and downloaded the publicly available chloroplast genome of Uncaria rhynchophylla (GenBank accession number MN723865) from NCBI (Bethesda, MD, USA) as a reference for annotation. The annotation results were checked and corrected manually by Geneious v11.0.5 software (Biomatters, Auckland, New Zealand) [83]. The structural characteristics of the cp genomes were demonstrated using the CPGview-RSG (Beijing, China).

4.3. Codon Usage Analyses

The genetic code is degenerate, that means one amino acid of multiple codons is called the synonymous codons of this amino acid [84]. In different conditions, the synonymous codons selected by the amino acids on the protein during translation are not the same, and they prefer to select some specific codons for different species [85]. The selection phenomenon is called codon preference [86]. CodonW (Nottingham, UK) software was used to calculate the relative synonymous codon usage (RSCU).

4.4. IR/SC Boundary Region

Angiosperm chloroplast genomes’ length variation is commonly due to the expansion and contraction of the IR regions and the single-copy boundary regions [87,88]. We used the IRscope (Helsinki, Finland) (https://irscope.shinyapps.io/irapp/, accessed on 20 July 2022) to compare IR/SC boundary region and Mauve (Madison, WI, USA) [89] to compare the chloroplast genomes for synteny analyses.

4.5. Characterization of Repeat Sequences and SSRs

We used the REPuter (Bielefeld, Germany) (https://bibiserv.cebitec.uni-bielefeld.de/reputer, accessed on 20 July 2022) to identify dispersed repeats with the cp genome, including forward repeats, reverse repeats, palindromic repeats, and complement repeats [90]. The following settings were used: (1) Minimal Repeat Size of 8 bp; (2) Hamming Distance of 3; (3) Maximum Computed Repeats of 50; and (4) 90% or greater sequence identity. To explore tandem repeats with the cp genome using Tandem Repeats Finder (New York, NY, USA) (https://tandem.bu.edu/trf/trf.html, accessed on 20 July 2022) [91]. The parameter settings were used: (1) Match of 2; (2) Mismatch of 7; (3) Indel of 7; (4) Alignment score of 50; and (5) Maximum period size of 500. Simple sequence repeats (SSRs) in the cp genome were found using MISA (Seeland, Germany) (http://pgrc.ipk-gatersleben.de/misa/, accessed on 20 July 2022) [92], with the parameters set at ten repeat units ≥10 for mononucleotides, eight repeat units ≥8 for dinucleotides, four repeat units ≥4 for trinucleotides, and three repeat units ≥3 for tetranucleotides, pentanucleotide, and hexanucleotide.

4.6. Predict RNA Editing Sites and Selective Pressure Analyses

The alteration of nucleotides at the RNA level by RNA editing is frequent in higher organisms [93]. The nucleotide substitutions, deletions, or insertions occur in mRNA after gene transcription, resulting in a change in the original genetic information [94,95]. According to previous research, chloroplast RNA editing controls the genes’ expression, allowing the genes to produce a variety of protein products and increase genetic information [96]. In order to ensure the accuracy of the predictions, PREP-cp (Lincoln, NE, USA) (http://prep.unl.edu/, accessed on 20 July 2022) [97] was used to predict RNA editing sites, using a parameter threshold (cut-off value) of 0.8. To calculate the rates of synonymous ratio (dS) and nonsynonymous ratio (dN) substitution and their ratio (ω = dN/dS), the protein-coding sequences of five Uncaria species were aligned using the MAFFT v7.3 [98]. Then, stop codons and gaps between the comparative sequences were deleted. dS/dN ratios were finally calculated using the CODEML program in PAML (London, UK) [99].

4.7. Sequence Divergence Analysis

The genome comparison among the five Uncaria species in the chloroplast genomes was analyzed using the mVISTA program (Berkeley, CA, USA) (https://genome.lbl.gov/vista/mvista/submit.shtml, accessed on 20 July 2022) as a point of comparison. The Shuffle-LAGAN mode was used with default parameters to align the five chloroplast genomes [100]. The nucleotide diversity of coding sections, introns, and intergenic spacers was evaluated using DnaSP v6 software (Barcelona, Spain) to determine the divergence hotspot locations in the five cp genomes [101]. The regions mentioned above were selected based on two factors: the total number of mutations (Eta) was greater than zero, and the aligned length was more than 200 bp. The following settings for sliding window analysis were used: (1) windows length of 600 bp and (2) step size of 200 bp.

4.8. Phylogenetic Analysis

The available chloroplast genomes sequences and protein-coding genes of 13 species from Rubiaceae were selected to study for the phylogenetic analysis, including Uncaria rhynchophylla, Uncaria sessilifructus, Uncaria macrophylla, Uncaria hirsuta, Uncaria scandens, Neolamarckia cadamba, Mitragyna speciosa, Cinchona officinalis, Coffea arabica, Coffea canephora, Gardenia jasminoides, Galium aparine, and Paederia scandens. Using Asclepias nivea of Asclepias as an outgroup, the phylogenetic tree was built. Geneious v11.0.5 (Biomatters, Auckland, New Zealand) [83] software was used to extract protein-coding genes from the chloroplast genome. The alignment of chloroplast genomes sequences and protein-coding genes was generated by MAFFT v7.3 (Suita, Osaka, Japan) with the default parameters [98]. We employed four methods to structure phylogenetic trees. Bayesian inference (BI) analyses were conducted in PhyloSuite v1.2.2 software (Wuhan, China) [102] with the best-fit GTR+F+I+G4 model. The maximum likelihood (ML) phylogenetic analyses were performed using IQTREE v1.6.7 (Vienna, Austria) [103]. Utilizing the best TVM+F+R2 model with 5000 bootstrap replicates, the maximum likelihood (ML) phylogenetic analyses were executed using IQTREE v1.6.7 [103]. Maximum parsimony (MP) and Neighbor-joining (NJ) analyses used MEGA-X (State College, PA, USA) [104] with the 1000 bootstrap replicates.

5. Conclusions

The whole cp genomes of four Uncaria species were firstly revealed in present study. They are conserved by comparing the complete chloroplast genome of five Uncaria, but some highly variable loci in coding regions and intergenic regions occurred. Nine highly variable intergenic regions (trnH-GUG-psbA, rps16-trnQ-UUG, trnQ-UUG-psbK, atpH-atpL, ndhC-trnV-UAC, petA-psbJ, rpl32-trnL-UAG, ndhF-rpl32, and rps15-ycf1) and two highly variable genic regions (ccsA and ycf1) were detected, these regions could serve as significant molecular markers for the Uncaria genus. We also found that the ndhB, psbB, ndhD, and atpB genes are characterized by a large number of potential RNA editing sites and 31-49 SSRs in five Uncaria species cp genomes. Compared to U. rhynchophylla, the other four Uncaria species showed gene or intron loss. These results will provide a theoretical basis for identifying U. rhynchophylla and its counterfeits, as well as the species of the Uncaria genus. Phylogenetic analysis brought unprecedented perspectives on the relationship among five Uncaria species. The availability of these findings will supplement important genetic knowledge and serve as a helpful guide for future studies on plastid genetic engineering research in Uncaria.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms231911617/s1.

Author Contributions

Methodology, software, validation, formal analysis, and writing—original draft preparation, M.-M.C.; data curation, investigation, conceptualization, and resources, M.Z.; funding acquisition and project administration, Z.-S.L.; funding acquisition, writing—review and editing, visualization, and supervision, Q.-L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Zhejiang Province, grant number LGN21H280004, and was funded by Zhejiang Sci-Tech University and Tasly Group Modern Chinese Medicine Resources Development Key Technology Research Center, grant number TZYHT201903250246.

Institutional Review Board Statement

We comply with the Convention on Biological Diversity and the Convention on the Trade in Endangered Species of Wild Fauna and Flora. This article does not involve any endangered or protected species. The research and collection of plant material were conducted according to the guidelines provided by Zhejiang Sci-Tech University. Uncaria plants were used in this study. Uncaria hirsuta, Uncaria macrophylla, and Uncaria scandens cultivars selected from Guangxi University of Chinese Medicine (Nanning, China) were kindly provided by Mei Ru (Guangxi University of Chinese Medicine, Nanning, China). Uncaria rhynchophylla and Uncaria sessilifructus cultivars selected from Guangxi Botanical Garden of Medicinal Plants (Nanning, China) were kindly provided by Yu-De Peng (Guangxi Botanical Garden of Medicinal Plants, Nanning, China).

Informed Consent Statement

The studies not involving humans.

Data Availability Statement

The chloroplast genome sequences of this study are openly available at the NCBI database (https://www.ncbi.nlm.nih.gov, accessed on 10 September 2022).

Acknowledgments

We thank the Guangxi University of Chinese Medicine (Mei Ru) and Guangxi Botanical Garden of Medicinal Plants (Yu-Dei Peng) for providing plant materials. At the same time, we thank the Tasly Group (Xue-Min Zhang) for coming up with this research’s initial idea and for this research being completed with the financial support of the Tasly Group.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, M.; Chen, M.M.; Zhang, X.M.; Chen, S.N.; Liang, Z.S. The complete chloroplast genome sequence of traditional Chinese medicine Uncaria macrophylla (Rubiaceae). Mitochondrial DNA Part B 2022, 7, 694–695. [Google Scholar] [CrossRef] [PubMed]
  2. Laus, G. Advances in Chemistry and Bioactivity of the Genus Uncaria. Phytother. Res. 2004, 18, 259–274. [Google Scholar] [CrossRef] [PubMed]
  3. Hao, H.F.; Liu, L.M.; Liu, Y.Y.; Liu, J.; Yan, L.; Pan, C.S.; Wang, M.X.; Wang, C.S.; Fan, J.Y.; Gao, Y.S.; et al. Inhibitory effect of rhynchophylline on contraction of cerebral arterioles to endothelin 1: Role of rho kinase—ScienceDirect. J. Ethnopharmacol. 2014, 155, 147–153. [Google Scholar] [CrossRef]
  4. Zhang, Q.; Chen, L.; Hu, L.J.; Liu, W.Y.; Feng, F.; Qu, W. Two new ortho benzoquinones from Uncaria rhynchophylla. Chin. J. Nat. Med. 2016, 14, 235. [Google Scholar] [CrossRef]
  5. Guo, Q.; Si, X.L.; Shi, Y.T.; Yang, H.S.; Liu, X.Y.; Liang, H.; Tu, P.F.; Zhang, Q.Y. Glucoconjugated Monoterpene Indole Alkaloids from Uncaria rhynchophylla. J. Nat. Prod. 2019, 82, 3288–3301. [Google Scholar] [CrossRef] [PubMed]
  6. Li, C.; Jiang, F.; Li, Y.L.; Jiang, Y.H.; Yang, W.Q.; Sheng, J.; Xu, W.J.; Zhu, Q.J. Rhynchophylla total alkaloid rescues autophagy, decreases oxidative stress and improves endothelial vasodilation in spontaneous hypertensive rats. Acta Pharmacol. Sin. 2018, 39, 345–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Sakakibara, I.; Terabayashi, S.; Kubo, M.; Higuchi, M.; Komatsu, Y.; Okada, M.; Taki, K.; Kamei, J. Effect on locomotion of indole alkaloids from the hooks of Uncaria plants. Phytomedicine 1999, 6, 163–168. [Google Scholar] [CrossRef]
  8. Hu, S.Q.; Mak, S.H.; Zuo, X.L.; Li, H.T.; Wang, Y.Q.; Han, Y.F. Neuroprotection Against MPP+-Induced Cytotoxicity Through the Activation of PI3-K/Akt/GSK3β/MEF2D Signaling Pathway by Rhynchophylline, the Major Tetracyclic Oxindole Alkaloid Isolated From Uncaria rhynchophylla. Front. Pharmacol. 2018, 9, 768. [Google Scholar] [CrossRef] [Green Version]
  9. Hanwool, L.; Seung, B.; Jong, L.; Chulwon, K.; Jeong Hyeon, K.; Seok Geun, L.; Arunachalam, C.; Sulaiman, A.; Yang, W.; Jae Young, U. Isorhynchophylline, a Potent Plant Alkaloid, Induces Apoptotic and Anti-Metastatic Effects in Human Hepatocellular Carcinoma Cells through the Modulation of Diverse Cell Signaling Cascades. Int. J. Mol. Sci. 2017, 18, 1095. [Google Scholar]
  10. Wicke, S.; Schneeweiss, G.M.; De Pamphilis, C.W.; Kai, F.M.; Quandt, D. The evolution of the plastid chromosome in land plants: Gene content, gene order, gene function. Plant Mol. Biol. 2011, 76, 273–297. [Google Scholar] [CrossRef] [Green Version]
  11. Ruhlman, T.A.; Jansen, R.K. Plastid Genomes of Flowering Plants: Essential Principles. Methods Mol. Biol. 2021, 2317, 3–47. [Google Scholar] [PubMed]
  12. Yang, Z.; Wang, G.; Ma, Q.; Ma, W.; Liang, L.; Zhao, T. The complete chloroplast genomes of three Betulaceae species: Implications for molecular phylogeny and historical biogeography. PeerJ 2019, 7, e6320. [Google Scholar] [CrossRef] [PubMed]
  13. Ruhfel, B.; Gitzendanner, M.A.; Soltis, P.S.; Soltis, D.E.; Burleigh, J.G. From algae to angiosperms–inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes. BMC Evol. Biol. 2014, 14, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chen, X.L.; Zhou, J.G.; Cui, Y.X.; Wang, Y.; Duan, B.Z.; Yao, H. Identification of Ligularia Herbs Using the Complete Chloroplast Genome as a Super-Barcode. Front. Pharmacol. 2018, 9, 695. [Google Scholar] [CrossRef] [Green Version]
  15. Li, L.; Hu, Y.; He, M.; Zhang, B.; Hong, Y. Comparative chloroplast genomes: Insights into the evolution of the chloroplast genome of Camellia sinensis and the phylogeny of Camellia. BMC Genom. 2020, 22, 138. [Google Scholar] [CrossRef]
  16. Ning, C.U.I.; Bao-Sheng, L.I.A.O.; Liang, C.L.; Shi-Feng, L.I.; Zhang, H.; Jiang, X.U.; Xi-Wen, L.I.; Chen, S.L. Complete chloroplast genome of Salvia plebeia:organization, specific barcode and phylogenetic analysis. Chin. J. Nat. Med. 2020, 18, 563–572. [Google Scholar]
  17. Zhou, J.H.; Zhang, J.; Chen, H.T.; Ma, L.N.; Liu, Y.S. Analysis of synonymous codon usage in foot-and-mouth disease virus. Vet. Res. Commun. 2010, 34, 393–404. [Google Scholar] [CrossRef]
  18. Nie, X.; Lv, S.; Zhang, Y.; Du, X.; Wang, L.; Biradar, S.S.; Tan, X.; Wan, F.; Song, W.; Sergios, O.K. Complete Chloroplast Genome Sequence of a Major Invasive Species, Crofton Weed (Ageratina adenophora). PLoS ONE 2012, 7, e36869. [Google Scholar] [CrossRef] [Green Version]
  19. Yi, X.; Lei, G.; Wang, B.; Su, Y.J.; Wang, T. The Complete Chloroplast Genome Sequence of Cephalotaxus oliveri (Cephalotaxaceae): Evolutionary Comparison of Cephalotaxus Chloroplast DNAs and Insights into the Loss of Inverted Repeat Copies in Gymnosperms. Genome Biol. Evol. Dev. 2013, 5, 688–698. [Google Scholar] [CrossRef] [Green Version]
  20. Yang, Q.; Fu, G.F.; Wu, Z.Q.; Li, L.; Zhao, J.L.; Li, Q.J. Chloroplast Genome Evolution in Four Montane Zingiberaceae Taxa in China. Front. Plant Sci. 2021, 12, 774482. [Google Scholar] [CrossRef]
  21. Hu, Y.; Xing, W.; Song, H.; Zhu, H.; Liu, G.; Hu, Z. Evolutionary Analysis of Unicellular Species in Chlamydomonadales Through Chloroplast Genome Comparison With the Colonial Volvocine Algae. Front. Microbiol. 2019, 10, 1351. [Google Scholar] [CrossRef] [PubMed]
  22. Ollivier, M.; Rispe, C. Evolutionary Dynamics in the Aphid Genome: Search for Genes Under Positive Selection and Detection of Gene Family Expansions. In Evolutionary Biology–Concepts, Molecular and Morphological Evolution; Springer: Berlin/Heidelberg, Germany, 2010; pp. 133–142. [Google Scholar]
  23. Wei, X.P.; Li, H.J.; Che, P.; Guo, H.J.; Zhang, B.G.; Liu, H.T.; Qi, Y.D. Comparing chloroplast genomes of traditional Chinese herbs Schisandra sphenanthera and S. chinensis. Chin. Herb. Med. 2020, 12, 48–57. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, Y.; Du, L.; Ao, L.; Chen, J.; Li, W.; Hu, W.; Zhang, W.; Kyunghee, K.; Sang-Choon, L.; Yang, T.J. The Complete Chloroplast Genome Sequences of Five Epimedium Species: Lights into Phylogenetic and Taxonomic Analyses. Front. Plant Sci. 2016, 7, 306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Moore, M.J.; Bell, C.D.; Soltis, P.S.; Soltis, D.E. Using plastid genome-scale data to resolve enigmatic relationships among basal angiosperms. Proc. Natl. Acad. Sci. USA 2007, 104, 19363–19368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zhang, Y.; Zhang, J.W.; Yang, Y.; Li, X.N. Structural and Comparative Analysis of the Complete Chloroplast Genome of a Mangrove Plant: Scyphiphora hydrophyllacea Gaertn. f. and Related Rubiaceae Species. Forests 2019, 10, 1000. [Google Scholar] [CrossRef] [Green Version]
  27. Gao, L.; Su, Y.J.; Wang, T. Plastid genome sequencing, comparative genomics, and phylogenomics:Current status and prospects. J. Syst. Evol. 2010, 48, 77–93. [Google Scholar] [CrossRef]
  28. Oyebanji, O.; Zhang, R.; Chen, S.Y.; Yi, T.S. New Insights Into the Plastome Evolution of the Millettioid Phaseoloid Clade (Papilionoideae, Leguminosae). Front. Plant Sci. 2021, 12, 652483. [Google Scholar] [CrossRef]
  29. Frailey, D.C.; Chaluvadi, S.R.; Vaughn, J.N.; Coatney, C.G.; Bennetzen, J.L. Gene loss and genome rearrangement in the plastids of five Hemiparasites in the family Orobanchaceae. BMC Plant Biol. 2018, 18, 30. [Google Scholar] [CrossRef] [Green Version]
  30. Li, X.; Li, Y.F.; Zang, S.P.S.M.; El-Kassaby, Y.A.; Fang, Y.M. Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus. Ecol. Evol. 2021, 11, 13401–13414. [Google Scholar] [CrossRef]
  31. Li, Y.; LÜ, G.H.; Zhang, X.N.; He, X.M. Chloroplast Genome Structure and Variation Analysis of Brassicaceae Species. Acta Bot. Boreali-Occident. Sin. 2017, 37, 1090–1101. [Google Scholar]
  32. Wang, W.; Yu, H.; Wang, J.; Lei, W.; Gao, J.; Qiu, X.; Wang, J. The Complete Chloroplast Genome Sequences of the Medicinal Plant Forsythia suspensa (Oleaceae). Int. J. Mol. Sci. 2017, 18, 2288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yi, D.K.; Lee, H.L.; Sun, B.Y.; Chung, M.Y.; Kim, K.J. The complete chloroplast DNA sequence of Eleutherococcus senticosus (Araliaceae); Comparative evolutionary analyses with other three asterids. Mol. Cells 2012, 33, 497–508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lei, W.J.; Ni, D.P.; Wang, Y.J.; Shao, J.J.; Liu, C. Intraspecific and heteroplasmic variations, gene losses and inversions in the chloroplast genome of Astragalus membranaceus. Sci. Rep. 2016, 6, 21669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Daniell, H.; Lin, C.S.; Yu, M.; Chang, W.J. Chloroplast genomes: Diversity, evolution, and applications in genetic engineering. Genome Biol. 2016, 17, 134. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, R.; Jiao, W.; Kai, H.; Ting, R.; Zeng, S.; Edward, B.; Liu, Z. Complete chloroplast genome sequence of Pedicularis cheilanthifolia, an alpine plant in China. Conserv. Genet. Resour. 2017, 9, 619–621. [Google Scholar] [CrossRef]
  37. Bock, R. Structure, function, and inheritance of plastid genomes. In Cell and Molecular Biology of Plastids; Springer: Berlin/Heidelberg, Germany, 2007; pp. 29–63. [Google Scholar]
  38. Ruhlman, T.A.; Jansen, R. The Plastid Genomes of Flowering Plants. Chloroplast Biotechnol. 2014, 1132, 3–38. [Google Scholar]
  39. Li, X.; Yang, J.B.; Wang, H.; Song, Y.; Yu, W.B. Plastid NDH pseudogenization and gene loss in a recently derived lineage from the largest hemiparasitic plant genus Pedicularis (Orobanchaceae). Plant Cell Physiol. 2021, 62, 971–984. [Google Scholar] [CrossRef]
  40. Hilu, K.W.; Borsch, T.; Müller, K.; Soltis, D.E.; Soltis, P.S.; Savolainen, V.; Chase, M.W.; Powell, M.P.; Alice, L.A.; Evans, R.; et al. Angiosperm phylogeny based on matK sequence information. Am. J. Bot. 2003, 90, 1758–1776. [Google Scholar] [CrossRef] [Green Version]
  41. Tamura, M.N.; Yamashita, J.; Fuse, S.; Haraguchi, M. Molecular phylogeny of monocotyledons inferred from combined analysis of plastid matK and rbcL gene sequences. J. Plant Res. 2004, 117, 109–120. [Google Scholar] [CrossRef]
  42. Graham, S.W.; Lam, V.K.Y.; Merckx, V.S.F.T. Plastomes on the edge: The evolutionary breakdown of mycoheterotroph plastid genomes. New Phytol. 2017, 214, 48–55. [Google Scholar] [CrossRef] [Green Version]
  43. Saldaña, C.L.; Grados, P.R.; Galarza, J.C.C.; Feijoo, S.; Abad, J.C.G.; Vásquez, H.V.; Maicelo, J.L.; Jhoncon, J.H.; Arbizu, C.I. Unlocking the Complete Chloroplast Genome of a Native Tree Species from the Amazon Basin, Capirona (Calycophyllum spruceanum, Rubiaceae), and Its Comparative Analysis with Other Ixoroideae Species. Genes 2022, 13, 113. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, R.; Li, Q.; Gao, J.; Qu, M.; Ding, P. The complete chloroplast genome sequence of the medicinal plant Morinda officinalis (Rubiaceae), an endemic to China. Mitochondrial DNA Part B 2016, 27, 4324–4325. [Google Scholar] [CrossRef] [PubMed]
  45. Duan, R.Y.; Huang, M.Y.; Yang, L.M.; Liu, Z.W. Characterization of the complete chloroplast genome of Emmenopterys henryi (Gentianales: Rubiaceae), an endangered relict tree species endemic to China. Conserv. Genet. Resour. 2017, 9, 459–461. [Google Scholar] [CrossRef]
  46. Mardanov, A.V.; Ravin, N.V.; Kuznetsov, B.B.; Samigullin, T.H.; Antonov, A.S.; Kolganova, T.V.; Skyabin, K.G. Complete Sequence of the Duckweed (Lemna minor) Chloroplast Genome: Structural Organization and Phylogenetic Relationships to Other Angiosperms. J. Mol. Evol. 2008, 66, 555–564. [Google Scholar] [CrossRef]
  47. Gao, L.; Yi, X.; Yang, Y.X.; Su, Y.J.; Wang, T. Complete chloroplast genome sequence of a tree fern Alsophila spinulosa: Insights into evolutionary changes in fern chloroplast genomes. BMC Evol. Biol. 2009, 9, 130. [Google Scholar] [CrossRef] [Green Version]
  48. Wanga, V.O.; Dong, X.; Oulo, M.A.; Mkala, E.M.; Wang, Q.F. Complete Chloroplast Genomes of Acanthochlamys bracteata (China) and Xerophyta (Africa) (Velloziaceae): Comparative Genomics and Phylogenomic Placement. Front. Plant Sci. 2021, 12, 1135. [Google Scholar] [CrossRef]
  49. Ren, J.; Tian, J.; Jiang, H.; Zhu, X.; Mutie, F.M.; Wanga, V.O.; Ding, S.; Yang, J.; Dong, X.; Chen, L.; et al. Comparative and Phylogenetic Analysis Based on the Chloroplast Genome of Coleanthus subtilis (Tratt.) Seidel, a Protected Rare Species of Monotypic Genus. Front. Plant Sci. 2022, 13, 828467. [Google Scholar] [CrossRef]
  50. Guo, S.; Liao, X.J.; Chen, S.Y.; Liao, B.S.; Guo, Y.M.; Cheng, R.Y.; Xiao, S.M.; Hu, H.Y.; Chen, J.; Pei, J.; et al. A Comparative Analysis of the Chloroplast Genomes of Four Polygonum Medicinal Plants. Front. Genet. 2022, 13, 764534. [Google Scholar] [CrossRef]
  51. Howe, C.J.; Barbrook, A.C.; Koumandou, V.L.; Nisbet, R.; Symington, H.A.; Wightman, T.F. Evolution of the chloroplast genome. Philos. Trans. R. Soc. B Biol. Sci. 2003, 358, 99–106. [Google Scholar] [CrossRef] [Green Version]
  52. Khrustalev, V.V.; Barkovsky, E.V. Mutational pressure is a cause of inter- and intragenomic differences in GC-content of simplex and varicello viruses. Comput. Biol. Chem. 2009, 33, 295–302. [Google Scholar] [CrossRef]
  53. Ding, S.X.; Dong, X.; Yang, J.X.; Guo, C.; Cao, B.B.; Guo, Y.; Hu, G.W. Complete Chloroplast Genome of Clethra fargesii Franch. an Original Sympetalous Plant from Central China: Comparative Analysis, Adaptive Evolution, and Phylogenetic Relationships. Forests 2021, 12, 441. [Google Scholar] [CrossRef]
  54. Kim, K.J.; Lee, H.L. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants. DNA Res. Int. J. Rapid Publ. Rep. Genes Genomes 2004, 11, 247–261. [Google Scholar] [CrossRef] [PubMed]
  55. Zhang, H.; Li, C.; Miao, H.; Xiong, S. Insights from the Complete Chloroplast Genome into the Evolution of Sesamum indicum L. PLoS ONE 2013, 8, e80508. [Google Scholar] [CrossRef] [PubMed]
  56. Zhou, M.; Guo, J.; Cha, J.; Chae, M.; She, C.; Barral, J.M.; Sachs, M.S.; Liu, Y. Non-optimal codon usage affects expression, structure and function of clock protein FRQ. Nature 2013, 495, 111–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Lyu, X.L.; Liu, Y. Non-optimal codon usage is critical for protein structure and function of the master general amino acid control regulator CPC-1. mBio 2020, 11, e2605–e2620. [Google Scholar] [CrossRef] [PubMed]
  58. Jian, H.Y.; Zhang, Y.H.; Yan, H.J.; Qiu, X.Q.; Wang, Q.G.; Li, S.B.; Zhang, S.D. The Complete Chloroplast Genome of a Key Ancestor of Modern Roses, Rosa chinensis var. spontanea, and a Comparison with Congeneric Species. Molecules 2018, 23, 389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Somaratne, Y.; Guan, D.L.; Wang, W.Q.; Zhao, L.; Xu, S.Q. The Complete Chloroplast Genomes of Two Lespedeza Species: Insights into Codon Usage Bias, RNA Editing Sites, and Phylogenetic Relationships in Desmodieae (Fabaceae: Papilionoideae). Plants 2020, 9, 51. [Google Scholar] [CrossRef] [Green Version]
  60. Kuang, D.Y.; Wu, H.; Wang, Y.L.; Gao, L.M.; Lu, L. Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): Implication for DNA barcoding and population. Genome 2011, 54, 663–673. [Google Scholar] [CrossRef] [Green Version]
  61. Dong, F.; Lin, Z.; Lin, J.; Ming, R.; Zhang, W. Chloroplast Genome of Rambutan and Comparative Analyses in Sapindaceae. Plants 2021, 10, 283. [Google Scholar] [CrossRef]
  62. Yu, X.; Zuo, L.; Lu, D.; Lu, B.; Yang, M.S.; Wang, J. Comparative analysis of chloroplast genomes of five Robinia species: Genome comparative and evolution analysis. Gene 2018, 689, 141–151. [Google Scholar] [CrossRef]
  63. Wang, M.; Liu, H.; Ge, L.; Xing, G.; Wang, M.; Weining, S.; Nie, X. Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L. Genes 2017, 8, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Josphat, S.; Li, Z.Z.; Andrew, G.; Liao, Y.Y. The Complete Chloroplast Genome Sequence of Tree of Heaven (Ailanthus altissima (Mill.) (Sapindales: Simaroubaceae), an Important Pantropical Tree. Int. J. Mol. Sci. 2018, 19, 929. [Google Scholar]
  65. Ibrar, A.; Biggs, P.J.; Matthews, P.J.; Collins, L.J.; Hendy, M.D.; Lockhart, P. Mutational Dynamics of Aroid Chloroplast Genomes. Genome Biol. Evol. Dev. 2012, 4, 1316–1323. [Google Scholar]
  66. Jungeun, L.; Kang, Y.; Chul, S.S.; Hyun, P.; Hyoungseok, L.; Szabolcs, S. Combined Analysis of the Chloroplast Genome and Transcriptome of the Antarctic Vascular Plant Deschampsia antarctica Desv. PLoS ONE 2014, 9, e92501. [Google Scholar]
  67. Liu, Q.; Li, X.Y.; Li, M.Z.; Xu, W.K.; Heslop, H.; John, S. Comparative chloroplast genome analyses of Avena: Insights into evolutionary dynamics and phylogeny. BMC Plant Biol. 2020, 20, 406. [Google Scholar] [CrossRef]
  68. Abdullah; Mehmood, F.; Shahzadi, I.; Ali, Z.; Islam, M.; Naeem, M.; Mirza, B.; Lockhart, P.J.; Ahmed, I.; Waheed, M.T. Correlations among oligonucleotide repeats, nucleotide substitutions, and insertion–deletion mutations in chloroplast genomes of plant family Malvaceae. J. Syst. Evol. Dev. 2020, 59, 388–402. [Google Scholar] [CrossRef]
  69. Liu, H.Y.; Yu, Y.; Deng, Y.Q.; Li, J.; Huang, Z.X.; Zhou, S.D. The Chloroplast Genome of Lilium henrici: Genome Structure and Comparative Analysis. Molecules 2018, 23, 1276. [Google Scholar] [CrossRef] [Green Version]
  70. Bejaoui, F.; Salas, J.J.; Nouairi, I.; Smaoui, A.; Abdelly, C.; Force, E.M.; Youssef, N.B. Changes in chloroplast lipid contents and chloroplast ultrastructure in Sulla carnosa and Sulla coronaria leaves under salt stress. J. Plant Physiol. 2016, 198, 32–38. [Google Scholar] [CrossRef]
  71. Dong, W.; Liu, H.; Xu, C.; Zuo, Y.; Chen, Z.; Zhou, S. A chloroplast genomic strategy for designing taxon specific DNA mini-barcodes: A case study on ginsengs. BMC Genet. 2014, 15, 138. [Google Scholar] [CrossRef] [Green Version]
  72. Jim, P.; Powell, W.; Hollingsworth, P.M. Chloroplast microsatellites: New tools for studies in plant ecology and evolution. Trends Ecol. Evol. 2001, 16, 142–147. [Google Scholar]
  73. Varshney, R.K.; Sigmund, R.; Börner, A.; Korzun, V.; Stein, N.; Sorrells, M.E.; Langridge, P.; Graner, A. Interspecific transferability and comparative mapping of barley EST-SSR markers in wheat, rye and rice. Plant Sci. 2005, 168, 195–202. [Google Scholar] [CrossRef]
  74. Rivas, J.D.L.; Lozano, J.; Ortiz, A.R. Comparative Analysis of Chloroplast Genomes: Functional Annotation, Genome-Based Phylogeny, and Deduced Evolutionary Patterns. Genome Res. 2002, 12, 567–583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Chen, M.M.; Wang, R.H.; Sha, H.K.; Liu, M.Z.; Tong, J.Q.; He, Q.L. The complete chloroplast genome sequence of Spiraea × vanhouttei (Briot) Zabel (Rosaceae). Mitochondrial DNA Part B 2022, 7, 505–506. [Google Scholar] [CrossRef] [PubMed]
  76. Wang, R.H.; Chen, M.M.; Wu, Y.Z.; Wang, R.R.; Xiao, X.F.; Qi, Z.C.; Yan, X.L. The complete chloroplast genome sequence of wild Japanese pepper Tubocapsicum anomalum Makino (Solanaceae). Mitochondrial DNA Part B 2021, 6, 2322–2323. [Google Scholar] [CrossRef] [PubMed]
  77. Zhou, X.; Huang, Y.B.; Zhang, Z.C.; Xu, X.Y.; Wang, R.H.; Xua, L.; Qi, Z.C.; Wei, Y.K. The complete chloroplast genome of endangered Zhangjiajie sage Salvia daiguii Y. K. Wei & Y. B. Huang (Lamiaceae). Mitochondrial DNA Part B 2020, 5, 3833–3834. [Google Scholar] [PubMed]
  78. Qi, Z.C.; Shen, C.; Han, Y.W.; Shen, W.; Yang, M.; Liu, J.L.; Liang, Z.S.; Li, P.; Fu, C.X. Development of microsatellite loci in Mediterranean sarsaparilla (Smilax aspera; Smilacaceae) using transcriptome data. Appl. Plant Sci. 2017, 5, 1700005. [Google Scholar] [CrossRef]
  79. Lu, Q.X.; Ren, L.J.; Wang, Y.N.; Xiao, Y.; Li, Y.D.; Ying, J.N.; Zhang, Y.; Wang, R.H.; Qi, Z.C.; You, Z.Y. Characterization of the complete chloroplast genome of Sparganium glomeratum (Typhaceae) from Jilin Province, China and phylogenetic analysis. Mitochondrial DNA Part B 2021, 6, 3253–3254. [Google Scholar] [CrossRef]
  80. Feng, J.Y.; Jin, X.J.; Zhang, S.L.; Yang, J.W.; Fei, S.P.; Huang, Y.S.; Liu, Y.; Qi, Z.C.; Li, P. Smilax weniae, a New Species of Smilacaceae from Limestone Areas Bordering Guizhou and Guangxi, China. Plants 2022, 11, 1032. [Google Scholar] [CrossRef]
  81. Jin, J.J.; Yu, W.B.; Yang, J.B.; Song, Y.; Li, D.Z. GetOrganelle: A fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020, 21, 241. [Google Scholar] [CrossRef]
  82. Shi, L.; Chen, H.; Jiang, M.; Wang, L.; Wu, X.; Huang, L.; Liu, C. CPGAVAS2, an integrated plastome sequence annotator and analyzer. Nuclc Acids Res. 2019, 47, W65–W73. [Google Scholar] [CrossRef]
  83. Drummond, A.J.B. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar]
  84. Uddin, A.; Choudhury, M.N.; Chakraborty, S. Factors influencing codon usage of mitochondrial ND1 gene in pisces, aves and mammals. Mitochondrion 2017, 37, 17–26. [Google Scholar] [CrossRef]
  85. Yu, T.H.; Li, J.S.; Yang, Y.; Qi, L.; Chen, B.B.; Zhao, F.Q.; Bao, Q.Y.; Wu, J.Y. Codon usage patterns and adaptive evolution of marine unicellular cyanobacteria Synechococcus and Prochlorococcus. Mol. Phylogenetics Evol. 2012, 62, 206–213. [Google Scholar] [CrossRef] [PubMed]
  86. Fei, L.; Gao, L.; Cheng, A.; Wang, M.; Han, Q. Analysis of Codon Usage Patterns of Tianfu Goose Interferon Alpha. Procedia Environ. Sci. 2011, 8, 730–736. [Google Scholar]
  87. Downie, S.R.; Jansen, R.K. A Comparative Analysis of Whole Plastid Genomes from the Apiales: Expansion and Contraction of the Inverted Repeat, Mitochondrial to Plastid Transfer of DNA, and Identification of Highly Divergent Noncoding Regions. Syst. Bot. 2015, 40, 336–351. [Google Scholar] [CrossRef]
  88. Lu, Q.; Ye, W.; Lu, R.; Xu, W.; Qiu, Y. Phylogenomic and Comparative Analyses of Complete Plastomes of Croomia and Stemona (Stemonaceae). Int. J. Mol. Sci. 2018, 19, 2383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Darling, A.C.E.; Mau, B.; Blattner, F.R.; Perna, N.T. Mauve: Multiple Alignment of Conserved Genomic Sequence with Rearrangements. Genome Res. 2004, 14, 1394–1403. [Google Scholar] [CrossRef] [Green Version]
  90. Kurtz, S.; Choudhuri, J.V.; Ohlebusch, E.; Schleiermacher, C.; Stoye, J.; Giegerich, R. REPuter: The manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001, 29, 4633–4642. [Google Scholar] [CrossRef] [Green Version]
  91. Gary, B. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar]
  92. Beier, S.; Thiel, T.; Münch, T.; Scholz, U.; Mascher, M. MISA-web: A web server for microsatellite prediction. Bioinformatics 2017, 33, 2583–2585. [Google Scholar] [CrossRef] [Green Version]
  93. Peng, X.X.; Xu, X.Y.; Wang, Y.M.; Hawke, D.H.; Yu, S.X.; Han, L.; Zhou, Z.C.; Mojumdar, K.; Jeong, K.J.; Labrie, M. A-to-I RNA Editing Contributes to Proteomic Diversity in Cancer. Cancer Cell 2018, 33, 817–828. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Fang, J.; Jiang, X.; Wang, T.; Deng, Z.; Zhang, A.; Zhang, X. Dynamic landscape of mitochondrial Cytidine-to-Uridine RNA editing in tobacco (Nicotiana tabacum) shows its tissue specificity. Plant Cell Tissue Organ Cult. 2022, 148, 363–376. [Google Scholar] [CrossRef]
  95. Brauer, N.L.; Alon, S.; Porath, H.T.; Elstein, B.; Unger, R.; Ziv, T.; Admon, A.; Levanon, E.Y.; Rosenthal, J.J.C.; Eisenberg, E. Trade-off between Transcriptome Plasticity and Genome Evolution in Cephalopods. Cell 2017, 169, 191–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Hanson, M.R.; Sutton, C.; Luis, B. Plant organelle gene expression: Altered by RNA editing. Trends Plant Sci. 1996, 1, 57–64. [Google Scholar] [CrossRef]
  97. Mower, J.P. The PREP suite: Predictive RNA editors for plant mitochondrial genes, chloroplast genes and user-defined alignments. Nucleic Acids Res. 2009, 37, W253–W259. [Google Scholar] [CrossRef]
  98. Katoh, K.; Standley, D. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 737–751. [Google Scholar] [CrossRef] [Green Version]
  99. Yang, Z.H. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef] [Green Version]
  100. Frazer, K.A.; Lior, P.; Alexander, P.; Rubin, E.M.; Inna, D. VISTA: Computational tools for comparative genomics. Nucleic Acids Res. 2004, 32, W273–W279. [Google Scholar] [CrossRef] [Green Version]
  101. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  102. Zhang, D.H.; Gao, F.L.; Jakovlić, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef]
  103. Nguyen, L.T.; Schmidt, H.A.; Haeseler, A.V.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  104. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 11617 g001 550
Figure 1. The leaf characteristics of five species in Uncaria. (A): U. rhynchophylla; (B): U. sessilifructus; (C): U. macrophylla; (D): U. hirsuta; (E): U. scandens.
Figure 1. The leaf characteristics of five species in Uncaria. (A): U. rhynchophylla; (B): U. sessilifructus; (C): U. macrophylla; (D): U. hirsuta; (E): U. scandens.
Ijms 23 11617 g001
Ijms 23 11617 g002 550
Figure 2. The complete chloroplast maps of five species in Uncaria. The gene in the circle is transcribed clockwise while the gene outside is transcribed counterclockwise. Genes are color-coded according to their roles. The AT content is shown in the lighter grey, while the GC content is shown in the deeper grey in the center circle.
Figure 2. The complete chloroplast maps of five species in Uncaria. The gene in the circle is transcribed clockwise while the gene outside is transcribed counterclockwise. Genes are color-coded according to their roles. The AT content is shown in the lighter grey, while the GC content is shown in the deeper grey in the center circle.
Ijms 23 11617 g002
Ijms 23 11617 g003 550
Figure 3. Amino acids and stop codons of codon numbers for the protein-coding regions in five Uncaria species. The histogram from left to right is U. rhynchophylla, U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens.
Figure 3. Amino acids and stop codons of codon numbers for the protein-coding regions in five Uncaria species. The histogram from left to right is U. rhynchophylla, U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens.
Ijms 23 11617 g003
Ijms 23 11617 g004 550
Figure 4. RSCU value comparison plots for amino acids in the protein-coding regions of five Uncaria species. The codon colors and the histogram colors match each other, their heights represent the RSCU value. The histogram from left to right is U. rhynchophylla, U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens.
Figure 4. RSCU value comparison plots for amino acids in the protein-coding regions of five Uncaria species. The codon colors and the histogram colors match each other, their heights represent the RSCU value. The histogram from left to right is U. rhynchophylla, U. sessilifructus, U. macrophylla, U. hirsuta, and U. scandens.
Ijms 23 11617 g004
Ijms 23 11617 g005 550
Figure 5. Dispersed repeated sequences analyses in cp genomes for five Uncaria species. (A) Statistics of four types of dispersed repeated sequences in five cp genomes; (B) Number of the dispersed repeated sequences with different period sizes.
Figure 5. Dispersed repeated sequences analyses in cp genomes for five Uncaria species. (A) Statistics of four types of dispersed repeated sequences in five cp genomes; (B) Number of the dispersed repeated sequences with different period sizes.
Ijms 23 11617 g005
Ijms 23 11617 g006 550
Figure 6. Tandem repeated sequence analyses for five Uncaria species in cp genomes. The number of tandem repeats with various period sizes, copy counts, and sequence lengths are shown in (A,B,C), respectively.
Figure 6. Tandem repeated sequence analyses for five Uncaria species in cp genomes. The number of tandem repeats with various period sizes, copy counts, and sequence lengths are shown in (A,B,C), respectively.
Ijms 23 11617 g006
Ijms 23 11617 g007 550
Figure 7. SSR analyses for five Uncaria species in cp genomes. (A) Number of different SSR types; (B) Frequency of identified SSRs in LSC, SSC, and IR regions; (C) Frequency of identified SSRs in CDS, IGS (intergenic spacer region), and intron; (D) Frequency of identified SSR in different repeat class types.
Figure 7. SSR analyses for five Uncaria species in cp genomes. (A) Number of different SSR types; (B) Frequency of identified SSRs in LSC, SSC, and IR regions; (C) Frequency of identified SSRs in CDS, IGS (intergenic spacer region), and intron; (D) Frequency of identified SSR in different repeat class types.
Ijms 23 11617 g007
Ijms 23 11617 g008 550
Figure 8. Comparison of LSC, IR, and SSC junction positions among five Uncaria species in cp genomes. JLB denotes the LSC/IRb junction, JSB denotes the SSC/IRb junction, JSA denotes the SSC/IRa junction, and JLA denotes the LSC/IRa junction.
Figure 8. Comparison of LSC, IR, and SSC junction positions among five Uncaria species in cp genomes. JLB denotes the LSC/IRb junction, JSB denotes the SSC/IRb junction, JSA denotes the SSC/IRa junction, and JLA denotes the LSC/IRa junction.
Ijms 23 11617 g008
Ijms 23 11617 g009 550
Figure 9. Synteny analyses for five Uncaria species in cp genomes.
Figure 9. Synteny analyses for five Uncaria species in cp genomes.
Ijms 23 11617 g009
Ijms 23 11617 g010 550
Figure 10. The dN/dS ratio values of common protein-coding genes in five Uncaria species chloroplast genomes.
Figure 10. The dN/dS ratio values of common protein-coding genes in five Uncaria species chloroplast genomes.
Ijms 23 11617 g010
Ijms 23 11617 g011 550
Figure 11. The dN/dS ratio values of common protein-coding genes in U. macrophylla, U. hirsuta, and U. rhynchophylla.
Figure 11. The dN/dS ratio values of common protein-coding genes in U. macrophylla, U. hirsuta, and U. rhynchophylla.
Ijms 23 11617 g011
Ijms 23 11617 g012 550
Figure 12. Sequence identity plots among five Uncaria species. The x-axis represents aligned base sequences, and the y-axis represents percent pairwise identity within 50–100%. (Coding regions were shown in blue, and non-coding regions were shown in pink. The locations and directions of each gene were indicated by gray arrows above the alignment. CNS: conserved noncoding sequences; UTR: untranslated region.).
Figure 12. Sequence identity plots among five Uncaria species. The x-axis represents aligned base sequences, and the y-axis represents percent pairwise identity within 50–100%. (Coding regions were shown in blue, and non-coding regions were shown in pink. The locations and directions of each gene were indicated by gray arrows above the alignment. CNS: conserved noncoding sequences; UTR: untranslated region.).
Ijms 23 11617 g012
Ijms 23 11617 g013 550
Figure 13. The nucleotide variability (Pi) values were compared among five Uncaria species. Position of the window midpoint on the X-axis and the nucleotide diversity within each window on the Y-axis (window length: 600 bp, step size: 200 bp).
Figure 13. The nucleotide variability (Pi) values were compared among five Uncaria species. Position of the window midpoint on the X-axis and the nucleotide diversity within each window on the Y-axis (window length: 600 bp, step size: 200 bp).
Ijms 23 11617 g013
Ijms 23 11617 g014 550
Figure 14. Phylogenetic trees based on 14 complete chloroplast genomes (A) and protein-coding genes (B) using different methods. From left to right, numbers above the branches are, respectively, represented Bayesian inference (BI) posterior probabilities (PPs) values, bootstrap support values of maximum likelihood (ML), neighbor-joining (NJ), and maximum parsimony (MP).
Figure 14. Phylogenetic trees based on 14 complete chloroplast genomes (A) and protein-coding genes (B) using different methods. From left to right, numbers above the branches are, respectively, represented Bayesian inference (BI) posterior probabilities (PPs) values, bootstrap support values of maximum likelihood (ML), neighbor-joining (NJ), and maximum parsimony (MP).
Ijms 23 11617 g014
Table
Table 1. Statistics of chloroplast genome sequencing of four species in Uncaria.
Table 1. Statistics of chloroplast genome sequencing of four species in Uncaria.
SampleU. sessilifructusU. macrophyllaU. hirsutaU. scandens
Clean reads18,093,60218,325,40114,436,56113,659,533
Clean bases2,300,585,4782,759,589,8672,174,504,8182,031,770,495
Clean data1.27 G2.06 G1.69 G1.68 G
Read length/bp150150150150
Q30(%)92.88%95.21%93.50%92.58%
Length/bp155,177155,138155,138153,780
GenBankON243635ON243636ON243634ON243637
Table
Table 2. Comparison of chloroplast genomic characteristics among five Uncaria species.
Table 2. Comparison of chloroplast genomic characteristics among five Uncaria species.
U. rhynchophyllaU. sessilifructusU. macrophyllaU. hirsutaU. scandens
Length (bp)
Total154,605155,177155,138155,138153,780
LSC85,146(55.07%)85,669(55.21%)85,749(55.27%)85,749(55.27%)85,465(55.58%)
SSC18,145(11.74%)18,128(11.68%)18,057(11.64%)18,057(11.64%)16,989(11.05%)
IR25,657(16.60%)25,690(16.56%)25,666(16.54%)25,666(16.54%)25,663(16.69%)
CDS80,517(52.08%)80,027(51.57%)80,069(51.61%)77,390(49.88%)77,526(50.41%)
tRNA2789(1.80%)2789(1.80%)2789(1.80%)2793(1.80%)2717(1.77%)
rRNA9048(5.85%)9048(5.83%)9048(5.83%)9070(5.85%)9048(5.88%)
GC content (%)
Overall37.737.637.537.537.7
LSC35.635.535.435.535.5
SSC31.731.531.531.731.7
IR43.243.243.243.243.2
CDS37.937.937.937.938.1
tRNA53.353.553.553.553.4
rRNA55.455.455.455.355.4
Number of genes
Overall131130130130128
CDS8685858583
tRNA3737373737
rRNA88888
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, M.-M.; Zhang, M.; Liang, Z.-S.; He, Q.-L. Characterization and Comparative Analysis of Chloroplast Genomes in Five Uncaria Species Endemic to China. Int. J. Mol. Sci. 2022, 23, 11617. https://doi.org/10.3390/ijms231911617

AMA Style

Chen M-M, Zhang M, Liang Z-S, He Q-L. Characterization and Comparative Analysis of Chloroplast Genomes in Five Uncaria Species Endemic to China. International Journal of Molecular Sciences. 2022; 23(19):11617. https://doi.org/10.3390/ijms231911617

Chicago/Turabian Style

Chen, Min-Min, Miao Zhang, Zong-Suo Liang, and Qiu-Ling He. 2022. "Characterization and Comparative Analysis of Chloroplast Genomes in Five Uncaria Species Endemic to China" International Journal of Molecular Sciences 23, no. 19: 11617. https://doi.org/10.3390/ijms231911617

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop