Next Article in Journal
Reproductive Outcomes in Adults with 22q11.2 Deletion Syndrome
Next Article in Special Issue
Interplay of Ecological Opportunities and Functional Traits Drives the Evolution and Diversification of Millettiod Legumes (Fabaceae)
Previous Article in Journal
Genomic Consequences of Fragmentation in the Endangered Fennoscandian Arctic Fox (Vulpes lagopus)
Previous Article in Special Issue
Potential of In Vitro Culture of Scutellaria baicalensis in the Formation of Genetic Variation Confirmed by ScoT Markers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 11
10.3390/genes13112125
genes-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 of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence

by
Dangwei Zhou
1,2,*,†,
Furrukh Mehmood
3,4,†,
Pengcheng Lin
1,
Tingfeng Cheng
2,
Huan Wang
2,
Shenbo Shi
2,
Jinkui Zhang
1,
Jing Meng
1,
Kun Zheng
1 and
Péter Poczai
5,*
1
The College of Pharmacy, Qinghai Nationalities University, Xining 810007, China
2
Key Laboratory of Adaptation and Evolution of Plateau Biota (AEPB), Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China
3
Department of Biochemistry, Faculty of Biological Sciences, Quaid-i-Azam University, Islamabad 45320, Pakistan
4
Department of Biochemistry, Faculty of Sciences, University of Sialkot, Daska Road, Punjab 51040, Pakistan
5
Faculty of Biological and Environmental Sciences, University of Helsinki, FI-00014 Helsinki, Finland
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2022, 13(11), 2125; https://doi.org/10.3390/genes13112125
Submission received: 3 October 2022 / Revised: 6 November 2022 / Accepted: 7 November 2022 / Published: 15 November 2022
(This article belongs to the Special Issue Phylogenetics, Genetics, and Breeding of Medicinal Plants)
Browse Figures
Review Reports Versions Notes

Abstract

:
Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai–Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai–Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

1. Introduction

Chloroplasts (cp) in plant cells are photosynthetic organelles that play a critical role in the production of essential energy and metabolites [1]. Cps contain a 120–160 kb genome with typical conserved quadripartite regions, including an LSC region, an SSC region, two identical copies of IR regions, and approximately 110–130 encoded genes distributed in a circular genome [2,3]. As most plant species exhibit maternal inheritance in their cps, the cp genome has undergone less recombination and evolution [3,4,5]. Therefore, cp sequences are considered a vital resource of molecular genetic markers to analyze the profiling of gene distribution and their molecular phylogenetic relationships [6]. Furthermore, because of the high level of expression, the absence of post-transcriptional gene silencing and site-specific transformation [7], cps are useful in the agriculture and pharmaceutical industries. In recent decades, rapid advances in next-generation sequencing technologies along with reduced costs have enabled a better understanding of the phylogenetic status and the cp transformation technique in traditional herbs using cp genome information. A. tanguticus is a rare and endangered Solanaceae species found on the Qinghai–Tibetan Plateau in grasslands or meadows (2200–4200 m) [8]. According to traditional Tibetan medicine, this species includes an abundance of hyoscyamine, atropine, and total alkaloids [9,10,11] and is widely used for the treatment of asthma, epilepsy, and inflammation [12,13]. Despite the development of a domesticated plant [11], the main medicinal resource of A. tanguticus still relies on wild plant exploration and collection, which has led to the species becoming endangered [14]. Furthermore, as populations are self-incompatible and interspecific hybridization occurs [15], the population structure and genetic diversity of A. tanguticus in the Qinghai–Tibetan region should be evaluated. DNA barcoding is thought to be an effective approach for medicinal component identification and quality control of produced goods in addition to protecting consumers [16,17]. Previous phylogenetic studies of Solanaceae plants have been based mostly on a few cp genes or ITS sequences [18,19]. There is currently limited cp genome information of A. tanguticus. Here, we assembled and de novo sequenced the whole cp genome of A. tanguticus with Illumina sequencing platforms. Additionally, cp genome annotation and structure analysis were performed and several oligonucleotide repeats, SSR, and substitutions were identified as mutation hotspots for future DNA barcoding. Further insights on A. tanguticus plastome evolution and its phylogenetic relationships have been revealed through comparative studies with reported cp genomes of Solanaceae species.

2. Materials and Methods

2.1. DNA Extraction

Fresh young leaves of A. tanguticus were collected from Haibei station (37.48° N, 101.2° E; Alt. 3200 m), in Qinghai province. According to the manufacturer’s instructions, a DNeasy Plant MiniKit (QIAGEN, Düsseldorf, Germany) was utilized to isolate DNA from the fresh young leaves. Spectrophotometry and electrophoresis on a 1% (w/v) agarose gel were used to evaluate the purity of the DNA. DNA of sufficient integrity and purity DNA was used for library construction and sequencing with an IlluminaHiseq2500 (San Diego, CA, USA).

2.2. Plastome Assembly and Annotation

Trimmomatic V0.36 was used to remove adapters and poor-quality reads from the raw reads [20]. Furthermore, clean reads were then mapped to the database, which was constructed from all cp genome sequences reported in the NCBI based on coverage depth and similarity. The mapped reads were then assembled into contigs [21] using SOAPdenovo2. The reads were filtered with careful error correction and different k-mers (55, 87, and 121) [22]. The scaffold of the cp genome was constructed with SSPACE [23] and GapCloser was used to fill the gaps [21]. PCR amplification and Sanger sequencing were used to validate the assembly by verifying the four boundaries between single copy (SC) and inverted repeat (IR) regions of the assembled sequences. Dual Organellar GenoMe Annotator (DOGMA, http://dogma.ccbb.utexas.edu/, accessed on 15 March 2021), CPGAVAS, and manual corrections were used to annotate the full cp genomes [24,25]. The software tRNAscan-SE was used for the prediction of tRNA genes. Circular cp genome map was drawn by Organellar Genome DRAW (OGDRAW) V1.2 [26]. The correct and complete cp genome was deposited in the National Center for Biotechnology Information (NCBI) under accession number MW246825 (A. tanguticus). Web-based REPuter (https://bibiserv.cebitec.uni-bielefeld.de/reputer/ accessed on 15 March 2021) was used for the analysis of repeat sequences, including forward, palindrome, reverse, and complement repeats. Sequence identity was set to >90% with a minimum repeat size of 30 base pairs. Micro Satellite Identification Tool (MISA) [27] was used to identify simple sequence repeats (SSRs), with minimum repeats of mono-, di-, tri-, tetra-, penta-, and hexa-nucleotides set to 10, 5, 4, 3, and 3, respectively. Using the tool CodonW1.4.2 (http://downloads.fyxm.net/CodonW-76666.html accessed on 15 March 2021), only 53 protein-coding genes with lengths more than 300 bp were chosen for analysis of synonymous codon usage. The relative synonymous codon usage (RSCU) and overall codon usage were both analyzed.

2.3. Comparative Analysis of cp Genomes

The plastome of A. tanguticus was compared with the cp genomes of Prewalskia tangutica, Scopolia parviflora, H. niger, Datura stramonium, and Nicotiana sylvestris in Solanaceae using Mauve software to identify evolutionary events such as gene loss, duplication, and rearrangement [28]. As a reference, the annotation of N. sylvestris was selected.

2.4. Molecular Evolution Analysis

The protein-coding genes of A. tanguticus were evaluated for synonymous (Ks) and non-synonymous (Ka) substitution rates in comparison with the six closely related Solanaceae species (P. tangutica, S. parviflora, H. niger, D. stramonium, and N. tabacum). Between the compared species, the corresponding functional protein-coding gene was separately aligned using MAFFT [29]. The Ka/Ks value was then determined using the NG method of computation, created by Nei and Gojobori and implemented in KaKs calculator 2.0 [30], with the settings genetic code Table S1 (bacterial and plant plastid code). In the results, some genes had Ka/Ks values of “NA”, which meant they were not relevant. When Ks = 0, this indicates no substitutions in the alignment, or 100% match. Concatenated sequences of 67 common protein-coding genes among the studied species were used in a phylogenetic analysis based on the complete cp genomes. The two sets of sequences were aligned by MAFFT [29], and the alignments were then adjusted by the Gblocks program [31]. Phylogenetic trees were created by the RAxMLversion 8.0 program using the GTRGAMMA model by the maximum likelihood (ML) approach [32]. Each branch’s bootstrap analysis was calculated using 1000 replications.

3. Results and Discussion

3.1. Chloroplast Genome Features of A. tanguticus

The cp genome of A. tanguticus is 155,765 nucleotides and has a quadripartite architecture. The genome consists of a pair of IRs (25,881 bp), a large single copy (LSC) region (86,516 bp), and a small single copy (SSC) region (17,487 bp) (Figure 1). These results are similar to other species of Solanaceae (Table 1). The cp genome’s total GC content was 37.63% and the GC content of the IR regions was 42.87%, which is higher than those of the LSC region (35.65%) and SSC region (31.93%). Similar results have also been observed in other species [32,33], which may be due to the four rRNAs and tRNAs [34,35]. The A. tanguticus plastome has 132 unique genes, including 87 protein-coding genes and 37 tRNA and eight rRNA genes. The cp genome contained 131 genes altogether, including 21 functional gene duplications in IR areas (Table 1). The plastome of A. tanguticus (87) has more protein-coding genes than P. tangutica (85) and H. niger (80) (Table 1). Among these, eight protein-coding genes (atpF, ndhA, ndhB, petD, rpoC1, rpl16, rpl2, and rps16) and six tRNAs had one intron, while three genes (ycf3, rps12, and clp2) contained two introns (Table 2). Nineteen genes, including four ribosomal (rrn4.5, rrn5, rrn16, and rrn23), seven trn (trn-UGC, trnI-CAU, trnI-GAU, trnL-CAA, trnN-GUU, trnR-ACG, and trnV-GAC), ndhB, rps7, rps12, rpl2, rpl23, ycf1, ycf2, and ycf15 had two repeats in the A. tanguticus plastome (Table 2). The largest intron (2502 bp) was found in the trnK-UUU gene, and it contained matK. These results are consistent with the mustard family (Brassicaceae) cp genome [35].

3.2. Codon Usage Analysis

RSCU frequency plays an important role in reflecting mutation bias during evolution [36]. RSCU values > 1 show the use of a codon more frequently than expected and vice versa [37]. The A. tanguticus plastome has a total of 80,910 bp of protein-coding genes, which encodes 26,970 codons. Among these, 2867 (6%) encoded leucine, which is the most abundant amino acid, whereas only 483 (1.0%) encoded tyrosine, which is the least prevalent amino acid in the A. tanguticus cp genome. According to the codon-anticodon recognition patterns, 29 tRNAs had codons that matched all 20 amino acids (Figure 2, Supplementary Materials, Table S1). AT content at the first, second, and third codon positions covered 54.3%, 61.9%, and 67.8%, respectively. This phenomenon was also higher in other plants [38,39]. There were 31 codons that exhibited clearly biased usage (RSCU > 1). Except for two codons (UUG for Leu and AUG for Met), all biased codon usage bases were A/U at the third position. Similarly, the codon for Trp did not exhibit biased usage (RSCU = 1) (Figure 2, Supplementary Materials, Table S1). This result is consistent with the N. officinale cp genome [35]. However, unlike these species, ATG showed the highest RSCU (1.9936) value of all codons in the A. tanguticus cp genome. The significant evolutionary character of the genome is evident in the codon usage bias [37,38]. Our results suggest that the A. tanguticus cp genome has experienced a different evolution process.

3.3. Repeat Sequence and SSR Analysis

Repeat sequences are a constant feature of illegitimate recombination and slipped-strand mispairing, and they can be utilized in plant phylogenetic studies [39,40]. The A. tanguticus cp genome contains 39 repeats, of which 18 are direct repeat sequences, 18 are palindromic sequences, and three are reverse repeat sequences. There were no complementary repeat sequences (Figure 3A). Analyses of oligonucleotide repeats revealed that all sequences were ≥30 bp, and direct and palindromic repeats covered 46% of the total dispersed repeat sequences in the cp genome. Interestingly, there was one palindromic sequence of 25,881 bp (Figure 3A). The number of forward and palindromic sequences was close among the different types of repeats, although they still differed among the six species (Figure 3B). Further analysis of repeat length among the four species of Hyoscyameae tribe revealed that the length of 40–49 bp in forward and palindromic sequences was greater than the other two species (Figure 3C,D), and the length of 50–59 bp may have played an important role in H. niger cp genome rearrangement (Figure 3C,D). However, the length of 30–40 bp for forward and palindromic sequences did not differ significantly among six species. Surprisingly, the number of reverse sequences with lengths of 30–39 bp was substantially greater in the cp genomes of A. tanguticus and P. tangutica, and there were no reverse sequences in the cp genomes of H. niger and D. stramonium (Figure 3E). Previous research has discovered long repetitions of >30 bp, which may contribute to cp genome rearrangement and increase population genetic diversity [41,42]. Our findings revealed that forward and palindromic sequences with lengths of 40–49 bp may play an essential role in the Hyoscyameae tribe, while reverse sequences with lengths of 30–39 bp are more abundant in A. tanguticus and P. tangutica than in the Hyoscyameae tribe. This may suggest that the reverse sequences of A. tanguticus and P. tangutica may play a more important role in cp genome rearrangement to adapt to the alpine environment on the Qinghai–Tibet Plateau. SSRs were 1–6 bp in the cp genome and were widely dispersed in intergenic, intron, or coding regions. SSRs had a high mutation rate and diverse copy number, which confers their important value in phylogenetic and population analyses [43,44]. Using MISA, 199 SSRs were identified in the A. tanguticus plastome (Supplementary Materials, Table S2). Among these, there were 120 mononucleotide, eight dinucleotide, 72 trinucleotide, seven tetranucleotide, and one pentanucleotide repeats. In addition, 60.3% mononucleotide repeats were A/T, which is consistent with results from Physalis species [45]. There were 61 trinucleotide SSRs in the A. tangutica cp genome, which is similar to four species of Hyoscyameae (Figure 3F) and more than that observed in four Physalis species [45]. Interestingly, dinucleotide SSRs were considerably lower than that of P. tangutica, S. parviflora, and H. niger, which contained 43, 43, and 42 dinucleotide SSRs, respectively. Although there was one pentanucleotide repeat in A. tanguticus, S. parviflora, and H. niger, there were three such repeats in the P. tangutica cp genome (Figure 3F). Therefore, these SSR differences may be helpful for phylogenetic and evolutionary studies. Furthermore, the LSC regions had 63.3% SSR markers, which was considerably greater than that of the SSC and IR regions. Similarly, nearly all SSRs in the A. tanguticus cp genome were comprised of A/T, which added 62.3% to the genome’s base composition.

3.4. IR Contraction and Expansion

The expansion and contraction of the IR boundaries were supposed to be an evolutionary trait of plastomes, and contraction and expansion at the IR region borders were believed to explain plastome differences [46,47]. To characterize the variation at the IR boundary, we selected the nine phylogenetically related Solanaceae species (P. peruviana, Iochroma ellipticum, Capsicum frutescens, Solanum lycopersicum, D. stramonium, N. sylvestris, H. niger, S. parviflora, and P. tangutica). Genomic structure and gene order are conserved in these species, especially in the Hyoscyamus group (Figure 4). The IRa/SSC border for the A. tanguticus cp genome was found in the 3’ region of the whole ycf1 gene. This is similar to H. niger, P. tangutica, and the other species except for S. parviflora. Similarly, except for S. parviflora and N. sylvestris, A. tanguticus and the other seven species contained the ycf1 pseudogene in the IRb region. This suggests that the ycf1 pseudogene in the IRb was lost in the evolution of the two species. In the cp genomes of A. tanguticus and the Solanaceae species mentioned above, the rps19 gene is located within the LSC/IRb border, except P. tangutica, where rps19 exists only in the LSC region. In some species, the 3′ region truncated rps19 pseudogene is in the IRa/LSC border [48]. Similarly, rps19 only appeared in A. tanguticus and S. parviflora in the IRa region. Our results indicate that the rps19 gene duplication is lost in the IRa/LSC border in the cp genomes of eight Solanaceae species (Figure 4). In the P. tangutica cp genome, there was also no rpl22 gene in the LSC region and rpl2 expanded to the IRb/LSC border (Figure 4). Similarly, the rpl2 gene is duplicated in A. tanguticus and five other species, while one duplication is lost in P. tangutica, H. niger, I. ellipticum, and P. peruviana. In general, our results showed that the IR boundary expansion and contraction differed from the Solanaceae species.

3.5. Comparative Chloroplast Genome Analysis

The whole cp genome of A. tanguticus was compared with five other species (P. tangutica, H. niger, S. parviflora, N. tabacum, and D. stramonium) using the Mauve program [28]. All cp genomes appeared as locally collinear blocks (Figure 5). The coding proteins and the tRNA and rRNA collinearly blocks are relatively conserved in structures and regions. However, gene loss occurred in the six species, such as 150,000–155,000, in which the P. tangutica and H. niger cp genomes lost two genes compared with the other four species. In addition, there was no gene rearrangement in Solanaceae genomes (Figure 5); this result is consistent with the other species [48,49]. The CD regions in the cp genomes of these six species were analyzed with MAFFT software and VCFtools was used to calculate the nucleotide variability (π) of each gene (Figure 6, Supplementary Materials, Table S3). The mean Pi value was 0.006936 among the six species. Similarly, IR regions showed less variability than that of the LSC and SSC regions. Among 113 genes, 26 genes showed no change in the six species, whereas 30 mutation hotspot regions were discovered (Pi > 0.01). These highly divergent proteins included matK, clpP, rps16, rps12, and ycf1, which had the highest Pi value (0.04246) (Figure 6, Supplementary Materials, Table S3). These regions might be going through rapid nucleotide substitution at the species level, suggesting the possible use of molecular markers for phylogenetic studies and plant identification.

3.6. Synonymous (Ks) and Non-Synonymous (Ka) Substitution Rate Analysis

Synonymous (Ks) and non-synonymous (Ka) were parameters for assessing the rate of gene divergence in a molecular evolution study [30]. To characterize the molecular evolution pattern, we calculated the Ka/Ks ratios of the A. tanguticus cp genome and the results were compared with those of five closely related species (P. tangutica, H. niger, S. parviflora, D. stramonium, and N. tabacum) (Figure 7; Supplementary Materials, Table S4). Of the 86 protein-coding genes, most had Ka/Ks values < 1. The average ratio was 0.24–0.32 for the five comparison groups. This shows that these species were mostly the outcome of vigorous purifying selection. In the Hyoscyameae tribe, cemA, rpoC2, and clpP exhibited positive selection in S. parviflora and P. tangutica. In a previous study, the clpP gene showed positive selection in various plant lineages, such as Oenothera [50], and cemA and rpoC2 also appeared to have undergone adaptation evolution in Gossypium and Populus, respectively [51,52]. Our results are consistent with these observations. In a recent study, clpP was shown to degrade misfolded proteins and protected cps from abiotic and biotic stress [53]; loss of this protein may lead to inappropriate structural changes in cps [54]. cemA proteins play an important role in protein sorting signals [55] and rpoC2 is an important subunit of PEP, which is required for transcription of photosynthesis genes and is involved in cp biogenesis [56]. In Nicotiana, positive selection acting on the ATP synthase and NAD(P)H dehydrogenase encoding genes was believed to be an adaptation to novel ecological conditions [57]. Here, our results indicate that protein metabolism may be an important adaptation strategy for A. tanguticus and other species in Hyoscyameae. As an endemic species on the Qinghai–Tibetan Plateau, A. tanguticus is often exposed to low temperatures, high UV-B radiation, and other abiotic stresses. Degrading misfolded proteins and a sorting signal system may help maintain plasmid function. Therefore, our results suggest that these genes are also important for Hyoscyameae tribe species evolution and adaptation to adverse environments.

3.7. Phylogenetic Analysis of the A. tanguticus cp Genome

The tribe Hyoscyameae is distributed throughout Eurasia and is considered as a monophyly in Solanaceae [18,19]. However, its genera relationship within this tribe remains unclear. Cp genomes provide whole phylogenetic information, which is suitable for establishing phylogenetic relationships in plants. Whole cp genomes and protein-coding genes of 21 species were utilized to create an ML phylogenetic tree to understand the phylogenetic position of A. tanguticus and to further elucidate the evolutionary relationship in the Solanaceae (Supplementary Materials, Table S5). Most nodes were strongly supported by 100% bootstrap value (Figure 8A). Moreover, A. tanguticus was clustered with three species (P. tangutica, H. niger, and S. parviflora) with 100% bootstrap support. The three Hyoscyameae tribe species were clustered together more closely. D. stramonium, N. sylvestris, and other species were clustered outside this branch with 100% bootstrap support. These findings are consistent with previous studies [18,19,57,58]. In a Hyoscyameae tribe study, A. tanguticus grouped with H. niger when several chloroplast markers were combined [18,19]. Our results are consistent with this observation. However, when whole genome data are used, P. tangutica and S. parviflora exhibited a closer relationship in the Hyoscyameae tribe (Figure 8B). This may be due to a data information shortage, such as in Physochlaina. Since A. tanguticus is an endemic species that grows on the Qinghai–Tibetan Plateau, sequencing, and characterization of additional cp genomes of Hyoscyameae tribe species may provide further insight on the phylogenetic relationships and evolution among these species.

4. Conclusions

High throughput sequencing was used to obtain the whole cp genome of A. tanguticus. The cp genome is 152,714 bp and encodes 126 functional genes. Analysis of repeat sequences revealed that forward and palindromic repeats of 40–49 bp were more common in the Hyoscyameae tribe, and 30–39 bp reverse repeats may play a more significant role in A. tanguticus and P. tangutica cp genome rearrangement and in adapting to extreme conditions on the Qinghai–Tibetan Plateau. The A. tanguticus cp genome contained a total of 199 SSRs and can be utilized for the development of molecular markers for phylogenetic and population studies. The genome structure and gene organization were similar to three previously reported species in the Hyoscyameae tribe, while gene contraction and expansion were clearly different in the A. tanguticus cp genome. Evolutionary analysis indicated that clpP, cemA, and ropC2 showed positive selection, which suggests that protein metabolism may be an important adaptation strategy for A. tanguticus and other species in Hyoscyameae on the Qinghai–Tibetan Plateau. Phylogenetic analysis indicated that A. tanguticus has a closer relationship with P. tangutica. These results provide significant information for forthcoming studies on phylogenetic and evolutionary adaptation in A. tanguticus and its related species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13112125/s1, Table S1: The codon usage and codon-anticodon recognition pattern for the A. tanguticus cp genome; Table S2a: Comparison of repeat sequences on A. tanguticus and five related species; Table S2b: Distribution of SSRs in the A. tanguticus cp genome; Table S3: Comparison of nucleotide variability (Pi) among A. tanguticus and related species; Table S4: The Ka/Ks ratios of protein-coding genes of A. tanguticus and related species; Table S5 The chloroplast genome of plant species for construction phylogenetic tree.

Author Contributions

D.Z. designed and managed the project; P.L., T.C. and H.W. performed the molecular experiment; F.M., S.S., J.Z., performed the bioinformation analysis; J.M. and K.Z. provided the sequencing data; D.Z., F.M. and P.P. wrote and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Training Qualified People Plan “Hope of West China” (Y229151211), the Foundation of Science in Qinghai (2021-ZJ-734), Construction Project for Innovation Plat form of Qinghai province, and the Innovation Group on Medicinal Resource Utilization (2021XJP102) and Institute fund of medicinal technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the corresponding author.

Acknowledgments

We acknowledge Nanjing Genepioneer Biotechnologies Inc. for sequencing of A. tanguticus and appreciate Chuanbei Jiang, Peipei Cao, Xiaolong Luo and Mingxing Yang for assistance in data analysis.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

chloroplast genome: cP; inverted repeated: IR; large single copy: LSC; small single copy: SSC; simple sequence repeats: SSR; intergenic spacer regions: IGS; relative synonymous codon usage: RSCU; maximum likelihood: ML; locally collinear blocks: LCBs; synonymous Ks; non-synonymous Ka.

References

  1. Rodríguez-Ezpeleta, N.; Brinkmann, H.; Burey, S.C.; Roure, B.; Burger, G.; Löffelhardt, W.; Bohnert, H.J.; Philippe, H.; Lang, B.F. Monophyly of primary photosynthetic eukaryotes: Green plants, red algae, and glaucophytes. Curr. Biol. 2005, 15, 1325–1330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Palmer, J.D. Comparative organization of chloroplast genomes. Annu. Rev. Genet. 1985, 19, 325–354. [Google Scholar] [CrossRef] [PubMed]
  3. Smith, D.R. Mutation rates in plastid genomes: They are lower than you might think. Genome Biol. Evol. 2015, 7, 1227–1234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wolfe, K.H.; Li, W.H.; Sharp, P.M. Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs. Proc. Natl. Acad. Sci. USA 1987, 84, 9054–9058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Duchene, D.; Bromham, L. Rates of molecular evolution and diversification in plants: Chloroplast substitution rates correlate with species-richness in the Proteaceae. BMC Evol. Biol. 2013, 13, 65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Shaw, J.; Lickey, E.B.; Schilling, E.E.; Small, R.L. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. Am. J. Bot. 2007, 94, 275–288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Maliga, P. Plastid transformation in higher plants. Annu. Rev. Plant Biol. 2004, 55, 289–313. [Google Scholar] [CrossRef]
  8. Liu, S.W. Flora of Qinghai; Qinghai People’s Publishing House: Xining, China, 1999. [Google Scholar]
  9. Yang, Y.C. Tibetan Medcine; Qinghai People’s Publishing House: Xining, China, 1990. [Google Scholar]
  10. Xiao, P.G.; Xia, G.C.; He, L.Y. The occurrence of some important tropane alkaloids in Chinese solanaceous plants. J. Integr. Plant Biol. 1973, 15, 187–194.4. [Google Scholar]
  11. Wang, H.; Zhang, X.F.; Chen, G.C.; Li, T.C.; Zhou, G.Y.; Shen, J.W. Comparative study of contents of four tropane alkaloids in cultural and wild Anisodus tanguticus. Acta Bot. Boreal.-Occid. Sin. 2005, 25, 575–577. [Google Scholar]
  12. Spinks, A.; Wasiak, J.; Villanueva, E.; Bernath, V. Scopolamine (hyoscine) for preventing and treating motion sickness. Cochrane Database Syst. Rev. 2007, 18, CD002851. [Google Scholar]
  13. Lei, T.X.; Cai, X.J.; Wang, H.; Li, S.L.; Shen, J.W.; Zhou, D.W. Progress on molecular mechanism of tropane alkaloids synthesis and plant bioengineering research. Acta Bot. Boreal. Occid. Sin. 2016, 36, 0204–0214. [Google Scholar]
  14. Zheng, W.; Wang, L.Y.; Meng, L.H.; Liu, J.Q. Genetic variation in the endangered Anisodus tanguticus (Solanaceae), an alpine perennial endemic to the Qinghai-Tibetan Plateau. Genetica 2008, 132, 123–129. [Google Scholar] [CrossRef] [PubMed]
  15. Duan, Y.W.; Zhang, T.F.; Liu, J.Q. Pollination biology of Anisodus tanguticus (Solanaceae). Biodivers. Sci. 2007, 15, 584–591. [Google Scholar]
  16. Li, M.; Cao, H.; But, P.P.H.; Shaw, P.C. Identification of herbal medicinal materials using DNA barcodes. J. Syst. Evol. 2011, 49, 271–283. [Google Scholar] [CrossRef]
  17. Raclariu, A.C.; Heinrich, M.; Ichim, M.C.; de Boer, H. Benefits and limitations of DNA barcoding and metabarcoding in herbal product authentication. Phytochem. Anal. 2018, 29, 123–128. [Google Scholar] [CrossRef] [PubMed]
  18. Tu, T.Y. Phylogeny and Biogeography of the tribes Nolaneae, Hyoscyameae and Mandragoreae of Solanaceae. Ph.D. Thesis, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China, 2008. [Google Scholar]
  19. Sanchez-Puerta, M.V.; Abbona, C.C. The Chloroplast genome of Hyoscyamus niger and a phylogenetic study of the tribe Hyoscyameae (Solanaceae). PLoS ONE 2014, 9, e98353. [Google Scholar] [CrossRef] [Green Version]
  20. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  21. Luo, R.; Liu, B.; Xie, Y.; Li, Z.; Huang, W.; Yuan, J.; He, G.; Chen, Y.; Qi, P.; Liu, Y. SOAPdenovo2: An empirically improved memory-efficient short-read de novo assembler. Gigascience 2012, 1, 18. [Google Scholar] [CrossRef]
  22. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  23. Boetzer, M.; Henkel, C.V.; Jansen, H.J.; Butler, D.; Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics 2011, 27, 578–579. [Google Scholar] [CrossRef] [Green Version]
  24. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Liu, C.; Shi, L.; Zhu, Y.; Chen, H.; Zhang, J.; Lin, X.; Guan, X. CPGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences. BMC Genom. 2012, 13, 715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lohse, M.; Drechsel, O.; Bock, R. OrganellarGenomeDRAW (OGDRAW): A tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr. Genet. 2007, 52, 7–274. [Google Scholar] [CrossRef] [PubMed]
  27. Thiel, T.; Michalek, W.; Varshney, R.; Graner, A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor. Appl. Genet. 2003, 106, 411–422. [Google Scholar] [CrossRef]
  28. 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] [PubMed]
  29. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, D.; Zhang, Y.; Zhang, Z.; Zhu, J.; Yu, J. KaKs_calculator 2.0: A toolkit incorporating gamma-series methods and sliding window strategies. Genom. Proteom. Bioinform. 2010, 8, 77–80. [Google Scholar] [CrossRef] [Green Version]
  31. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [Green Version]
  32. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, X.; Rong, C.X.; Qin, L.; Mo, C.Y.; Fan, L.; Yan, J.; Zhang, M.R. Complete chloroplast genome sequence of Malus hupehensis: Genome Structure, comparative analysis, and phylogenetic relationships. Molecules 2018, 23, 2917. [Google Scholar] [CrossRef] [Green Version]
  34. Shen, X.; Wu, M.; Liao, B.; Liu, Z.; Bai, R.; Xiao, S.; Li, X.; Zhang, B.; Xu, J.; Chen, S. Complete chloroplast genome sequence and phylogenetic analysis of the medicinal plant Artemisia annua. Molecules 2017, 22, 1330. [Google Scholar] [CrossRef] [PubMed]
  35. Yan, C.; Du, J.C.; Gao, L.; Li, Y.; Hou, X.L. The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae. Gene 2019, 699, 24–36. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, Q.; Xue, Q. Comparative studies on codon usage pattern of chloroplasts and their host nuclear genes in four plant species. J. Genet. 2005, 84, 55–62. [Google Scholar] [PubMed]
  37. Sharp, P.M.; Li, W.H. The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15, 1281–1295. [Google Scholar] [CrossRef] [Green Version]
  38. Redwan, R.M.; Saidin, A.; Kumar, S.V. Complete chloroplast genome sequence of MD-2 pineapple and its comparative analysis among nine other plants from the subclass Commelinidae. BMC Plant Biol. 2015, 15, 196. [Google Scholar] [CrossRef]
  39. Cavaliersmith, T. Chloroplast evolution: Secondary symbiogenesis and multiple losses. Curr. Biol. 2002, 12, R62–R64. [Google Scholar] [CrossRef] [Green Version]
  40. 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] [Green Version]
  41. Nie, X.; Lv, S.; Zhang, Y.; Du, X.; Wang, L.; Biradar, S.S.; Tan, X.; Wan, F.; Weining, S. Complete chloroplast genome sequence of a major invasive species, Croftonweed (Ageratina adenophora). PLoS ONE 2012, 7, e36869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Timme, R.E.; Kuehl, J.V.; Boore, J.L.; Jansen, R.K. A comparative analysis of the Lactuca and Helianthus (Asteraceae) plastid genomes: Identification of divergent regions and categorization of shared repeats. Am. J. Bot. 2007, 94, 302–312. [Google Scholar] [CrossRef]
  43. Huang, H.; Shi, C.; Liu, Y.; Mao, S.Y.; Gao, L.Z. Thirteen Camellia chloroplast genome sequences determined by high-throughput sequencing: Genome structure and phylogenetic relationships. BMC Evol. Biol. 2014, 14, 151. [Google Scholar] [CrossRef] [Green Version]
  44. Zhao, Y.B.; Yin, J.L.; Guo, H.Y.; Zhang, Y.Y.; Xiao, W.; Sun, C.; Wu, J.Y.; Qu, X.B.; Yu, J.; Wang, X.M.; et al. The complete chloroplast genome provides insight into the evolution and polymorphism of Panax ginseng. Front. Plant Sci. 2015, 5, 696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Feng, S.G.; Zheng, K.X.; Jiao, K.L.; Cai, Y.C.; Chen, C.L.; Mao, Y.Y.; Wang, L.Y.; Zhan, X.R.; Ying, Q.C.; Wang, H.Z. Complete chloroplast genomes of four Physalis species (Solanaceae): Lights into genome structure, comparative analysis, and phylogenetic relationships. BMC Plant Biol. 2020, 20, 242. [Google Scholar] [CrossRef] [PubMed]
  46. Raubeson, L.A.; Peery, R.; Chumley, T.W.; Dziubek, C.; Fourcade, H.M.; Boore, J.L.; Jansen, R.K. Comparative chloroplast genomics: Analyses including new sequences from the angiosperms Nuphar advena and Ranunculus macranthus. BMC Genom. 2007, 8, 174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Wang, R.J.; Cheng, C.L.; Chang, C.C.; Wu, C.L.; Su, T.M.; Chaw, S.M. Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots. BMC Evol. Biol. 2008, 8, 36. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, J.; Yue, M.; Niu, C.; Ma, X.F.; Li, Z.H. Comparative analysis of the complete chloroplast genome of four endangered herbals of Notopterygium. Genes 2017, 8, 124. [Google Scholar] [CrossRef]
  49. Zhao, K.H.; Li, L.Q.; Lu, Y.Z.; Yang, J.B.; Zhang, Z.R.; Zhao, F.Y.; Quan, H.; Ma, X.J.; Liao, Z.H.; Lan, X.Z. Characterization and comparative analysis of two Rheum complete chloroplast genomes. BioMed. Res. Int. 2020, 2020, 6490164. [Google Scholar] [CrossRef] [PubMed]
  50. Oxelman, B.; Erixon, P. Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene. PLoS ONE 2008, 3, e1386. [Google Scholar]
  51. Zong, D.; Zhou, A.P.; Zhang, Y.; Zou, X.L.; Li, D.; Duan, A.N.; He, C.Z. Characterization of the complete chloroplast genomes of five Populus species from the western Sichuan plateau, southwest China: Comparative and phylogenetic analyses. Peer J. 2019, 7, e6386. [Google Scholar]
  52. Wu, Y.; Liu, F.; Yang, D.G.; Li, W.; Zhou, X.J.; Pei, X.Y.; Liu, Y.G.; He, K.L.; Zhang, W.S.; Ren, Z.Y.; et al. Comparative chloroplast genomics of Gossypium species: Insights into repeat sequence variations and phylogeny. Front. Plant Sci. 2018, 9, 376. [Google Scholar] [CrossRef] [Green Version]
  53. Ali, M.D.; Baek, K.H. Protective roles of cytosolic and plastidal proteasomes on abiotic stress and pathogene invasion. Plants 2020, 9, 832. [Google Scholar] [CrossRef]
  54. Ali, M.D.; Baek, K.H. Co-Suppression of NbClpC1 and NbClpC2, encoding Clp protease chaperons, elicits significant changes in the metabolic profile of Nicotiana benthamiana. Plants 2020, 9, 259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zoschke, R.; Barkan, A. Genome-wide analysis of thylakoid-bound ribosomes in maize reveals principles of cotranslational targeting to the thylakoid membrane. Proc. Natl. Acad. Sci. USA 2015, 112, E1678–E1687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Börner, T.; Yu, A.; Yan, O.; Zubo, A.; Victor, V.; Kusnetsov, V.V. Chloroplast RNA polymerases: Role in chloroplast biogenesis. Biochim. Biophys. Acta 2015, 1874, 761–769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Mehmood, F.; Abdullah.; Ubaid, Z.; Shahzadi, I.; Ahmed, I.; Waheed, M.T.; Poczai, P.; Mirza, B. Plastid genomics of Nicotiana (Solanaceae): Insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica). Peer J. 2020, 8, e9552. [Google Scholar] [CrossRef]
  58. Zhang, D.; Fan, L.Q.; Zhang, L.; Wang, T.J.; Guo, X.Y.; Hu, Q.J. Characterization of the complete chloroplast genome of the endangered Przewalskia tangutica Maxim. Conserv. Genet. Resour. 2017, 9, 409–413. [Google Scholar] [CrossRef]
Genes 13 02125 g001 550
Figure 1. Gene map of the A. tanguticus chloroplast genome. (Genes shown outside the outer circle are transcribed clockwise and those insides are transcribed counterclockwise. The colored bars indicated different functional groups).
Figure 1. Gene map of the A. tanguticus chloroplast genome. (Genes shown outside the outer circle are transcribed clockwise and those insides are transcribed counterclockwise. The colored bars indicated different functional groups).
Genes 13 02125 g001
Genes 13 02125 g002 550
Figure 2. Codon content of 20 amino acids and the stop codon of 81 coding genes of the A. tanguticus.
Figure 2. Codon content of 20 amino acids and the stop codon of 81 coding genes of the A. tanguticus.
Genes 13 02125 g002
Genes 13 02125 g003 550
Figure 3. Repeat sequences and SSR analysis of A. tanguticus and related species of Solanaceae. (A) Profiling of repeat sequences of A. tanguticus; (B) number of four types of repeat of A. tanguticus and five related species; (C) comparison of a forward repeat of A. tanguticus and five related species; (D) comparison of a palindromic repeat of A. tanguticus and five related species; (E) comparison of a reverse repeat of A. tanguticus and five related species; (F) SSR number of A. tanguticus and three related species.
Figure 3. Repeat sequences and SSR analysis of A. tanguticus and related species of Solanaceae. (A) Profiling of repeat sequences of A. tanguticus; (B) number of four types of repeat of A. tanguticus and five related species; (C) comparison of a forward repeat of A. tanguticus and five related species; (D) comparison of a palindromic repeat of A. tanguticus and five related species; (E) comparison of a reverse repeat of A. tanguticus and five related species; (F) SSR number of A. tanguticus and three related species.
Genes 13 02125 g003
Genes 13 02125 g004 550
Figure 4. Comparison of the border position of LSC, SSC, and LR regions among A. tanguticus and nine species’ chloroplast genomes in Solanaceae. Genes are represented by boxes above (negative strand) and below (positive strand) the proportionate line. The size of each gene and their relative position at the junctions are shown in base pairs (bp).
Figure 4. Comparison of the border position of LSC, SSC, and LR regions among A. tanguticus and nine species’ chloroplast genomes in Solanaceae. Genes are represented by boxes above (negative strand) and below (positive strand) the proportionate line. The size of each gene and their relative position at the junctions are shown in base pairs (bp).
Genes 13 02125 g004
Genes 13 02125 g005 550
Figure 5. Sequence alignment of four chloroplast genomes in Solanaceae with A. tanguticus by Mauve software. Within each of the alignments, local collinear blocks are presented by blocks of the same color connected by lines. Annotations of rRNA, protein-coding, and tRNA genes are shown in red, white, and green boxes, respectively.
Figure 5. Sequence alignment of four chloroplast genomes in Solanaceae with A. tanguticus by Mauve software. Within each of the alignments, local collinear blocks are presented by blocks of the same color connected by lines. Annotations of rRNA, protein-coding, and tRNA genes are shown in red, white, and green boxes, respectively.
Genes 13 02125 g005
Genes 13 02125 g006 550
Figure 6. Comparison of nucleotide variability (Pi) among A. tanguticus and related species.
Figure 6. Comparison of nucleotide variability (Pi) among A. tanguticus and related species.
Genes 13 02125 g006
Genes 13 02125 g007 550
Figure 7. The Ka/Ks ratios of protein-coding genes of A. tanguticus and related species.
Figure 7. The Ka/Ks ratios of protein-coding genes of A. tanguticus and related species.
Genes 13 02125 g007
Genes 13 02125 g008a 550Genes 13 02125 g008b 550
Figure 8. The phylogenetic tree base on 20 cp chloroplast species was constructed with ML methods. Bootstrap values were shown at the nodes. ((A). ML tree of coding genes; (B). ML tree of whole cp genomes).
Figure 8. The phylogenetic tree base on 20 cp chloroplast species was constructed with ML methods. Bootstrap values were shown at the nodes. ((A). ML tree of coding genes; (B). ML tree of whole cp genomes).
Genes 13 02125 g008aGenes 13 02125 g008b
Table
Table 1. Comparison analyses of A. tanguticus with six close species.
Table 1. Comparison analyses of A. tanguticus with six close species.
Genome FeatureA. tanguticusP. tanguticaS. parvifloraH. nigerD. stramoniumN. tabacumS. tuberosum
Genome Size (bp)155,765155,569156,193155,720155,871155, 943155,298
LSC (bp)86,51686,70786,36486,10586,30186, 68685,749
SSC (bp)17,48718,28825,90517,86318,36618, 57318,373
IR (bp)25,88125,28725,87625,87625,60225,34225,595
GC content (%)37.6337.637.637.637.937.937.9
Total number
of genes		132138131118134156130
Protein-coding gene878586808811181
tRNA37443730383737
rRNA8884888
Table
Table 2. List of genes in the cp genome of A. tanguticus.
Table 2. List of genes in the cp genome of A. tanguticus.
Category for GenesGroups of GenesName of Genes
PhotosynthesisSubunits of photosystem ⅠpsaA, psaB, psaC, psaI, psaJ
Subunits of photosystemⅡpsbA, psbB, psbC, psbD, psbE, psbF, psbH
psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ
Subunits of cytochrome b/f complexpetA, petB, petD *, petG, petL, petN
Large subunit of RubiscorbcL
Subunits of ATP synthase atpA, atpB, atpE, atpF *, atpH, atpI
Subunits of NADH-dehydrogenasendhA *, ndhB a*, ndhC, ndhD, ndhE, ndhF, ndhG, ndhHndhI, ndhJ, ndhK
Self-replicationRibosomal RNA genes rrn16 a, rrn23 a, rrn5 a, rrn4.5 a
Transfer RNA genestrnA-UGC *a, trnC-GCA, trnD-GUC, trnE-UUC *a, trnF-GAA, trnfM-CAU, trnG-GCC, trnG-UCC *, trnH-GUG, trnI-CAU a, trnI-GAU a*, trnK-UUU *, trnL-CAA a, trnL-UAG, trnL-UAA *, trnM-CAU, trnN-GUU a, trnP-UGG, trnQ-UUG, trnR-ACG a,
trnR-UCU, trnS-GCU, trnS-GGA, trnS-UGA, trnY-GUA, trnT-GGU, trnT-UGU, trnV-GAC a, trnV-UAC *, trnW-CCA
Small subunit of ribosomerps2, rps3, rps4, rps7 a, rps8, rps11, rps12 a**, rps14, rps15, rps16, rps18, rps19
Large subunit of ribosomerpl2 a*, rpl14, rpl16 *, rpl20, rpl22, rpl23 a, rpl32, rpl33, rpl36
DNA-dependent RNA polymeraserpoA, rpoB, rpoC1 *, rpoC2
Other genesMaturasematK
Envelope membrane proteincemA
Subunit of acetyl-CoAaccD
C-type cytochrome synthesis geneccsA
ProteaseclpP **
Genes of unknown functionConserved Open reading framesycf1 a, ycf2 a, ycf3 **, ycf4, ycf15 a
Note: * means the gene contained one intron; ** means the gene contained two introns; a indicates the number of the repeat unit is 2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, D.; Mehmood, F.; Lin, P.; Cheng, T.; Wang, H.; Shi, S.; Zhang, J.; Meng, J.; Zheng, K.; Poczai, P. Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence. Genes 2022, 13, 2125. https://doi.org/10.3390/genes13112125

AMA Style

Zhou D, Mehmood F, Lin P, Cheng T, Wang H, Shi S, Zhang J, Meng J, Zheng K, Poczai P. Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence. Genes. 2022; 13(11):2125. https://doi.org/10.3390/genes13112125

Chicago/Turabian Style

Zhou, Dangwei, Furrukh Mehmood, Pengcheng Lin, Tingfeng Cheng, Huan Wang, Shenbo Shi, Jinkui Zhang, Jing Meng, Kun Zheng, and Péter Poczai. 2022. "Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence" Genes 13, no. 11: 2125. https://doi.org/10.3390/genes13112125

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