Next Article in Journal
Deciphering the Molecular Recognition Mechanism of Multidrug Resistance Staphylococcus aureus NorA Efflux Pump Using a Supervised Molecular Dynamics Approach
Next Article in Special Issue
Plastid Genomes of Carnivorous Plants Drosera rotundifolia and Nepenthes × ventrata Reveal Evolutionary Patterns Resembling Those Observed in Parasitic Plants
Previous Article in Journal
miR526b and miR655 Induce Oxidative Stress in Breast Cancer
Previous Article in Special Issue
Development of Plastid Genomic Resources for Discrimination and Classification of Epimedium wushanense (Berberidaceae)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 16
10.3390/ijms20164040
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

Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species

by
Yingxian Cui
1,2,
Xinlian Chen
1,2,
Liping Nie
1,2,
Wei Sun
3,
Haoyu Hu
3,
Yulin Lin
1,2,
Haitao Li
4,
Xilong Zheng
5,
Jingyuan Song
1,2 and
Hui Yao
1,2,*
1
Key Lab of Chinese Medicine Resources Conservation, State Administration of Traditional Chinese Medicine of the People’s Republic of China, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100193, China
2
Engineering Research Center of Chinese Medicine Resources, Ministry of Education, Beijing 100193, China
3
Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing 100700, China
4
Yunnan Branch, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Jinghong 666100, China
5
Hainan Branch, Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences & Peking Union Medical College, Wanning 571533, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(16), 4040; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20164040
Submission received: 24 June 2019 / Revised: 15 August 2019 / Accepted: 17 August 2019 / Published: 19 August 2019
(This article belongs to the Special Issue Chloroplast 2.0)
Browse Figures
Versions Notes

Abstract

:
Amomum villosum is an important medicinal and edible plant with several pharmacologically active volatile oils. However, identifying A. villosum from A. villosum var. xanthioides and A. longiligulare which exhibit similar morphological characteristics to A. villosum, is difficult. The main goal of this study, therefore, is to mine genetic resources and improve molecular methods that could be used to distinguish these species. A total of eight complete chloroplasts (cp) genomes of these Amomum species which were collected from the main producing areas in China were determined to be 163,608–164,069 bp in size. All genomes displayed a typical quadripartite structure with a pair of inverted repeat (IR) regions (29,820–29,959 bp) that separated a large single copy (LSC) region (88,680–88,857 bp) from a small single copy (SSC) region (15,288–15,369 bp). Each genome encodes 113 different genes with 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. More than 150 SSRs were identified in the entire cp genomes of these three species. The Sanger sequencing results based on 32 Amomum samples indicated that five highly divergent regions screened from cp genomes could not be used to distinguish Amomum species. Phylogenetic analysis showed that the cp genomes could not only accurately identify Amomum species, but also provide a solid foundation for the establishment of phylogenetic relationships of Amomum species. The availability of cp genome resources and the comparative analysis is beneficial for species authentication and phylogenetic analysis in Amomum.

1. Introduction

Amomum villosum Lour., which belongs to the monophyletic Zingiberaceae family, is a valuable medicinal plant in China with a history of more than 1300 years [1]. As a plant with both medicinal and edible functions, its ripe fruits, called amomi fructus, are not only widely used clinically but also often used as a type of condiment. In China, amomi fructus is one of the “Four Major Southern Medicines” and plays important roles in clinical treatment, such as warming the spleen, eliminating dampness and diarrhea, promoting appetite, and preventing miscarriage [2]. Amomi fructus possesses potential therapeutic effect for inflammatory bowel disease [3], radical scavenging [4], and has analgesic and anti-inflammatory effects [5], which have recently gained increasing attention. Three authentic plant sources of amomi fructus have been recorded in the Chinese pharmacopoeia, namely, A. villosum Lour., A. villosum Lour. var. xanthioides T. L. Wu et Senjen and A. longiligulare T. L. Wu [2]. However, research has shown that the content of the active ingredient (borneol acetate) is significantly higher in A. villosum than in A. villosum var. xanthioides and A. longiligulare, which is consistent with the statement that the effects of A. villosum are better than those of the other two species [6,7], resulting in the high value of A. villosum in clinical use. Given the special flower structure of A. villosum, the plant exhibits difficulty in self-pollination and must rely on a few kinds of insects or artificial pollination. This pollination characteristic greatly reduces the yield of A. villosum [6]. Given their high economic and medicinal values, the ripe fruits of A. villosum var. xanthioides and A. longiligulare, which exhibit high similarity in morphological characters and confusing Chinese names with A. villosum, are often used as adulterants or contaminants of A. villosum in the market (Figure 1) [8,9]. In addition, the appearances of these three Amomum plants are extremely similar, and there are only a few differences on their leaves, which have been marked in Figure 1 with red boxes. Generally speaking, the ligules of A. villosum are purple and relatively short; the ligules sizes of A. villosum var. xanthioides are almost similar to those of A. villosum, but the color is green; and those of A. longiligulare are also purple in color, but their size is longer than the other two. The mixture of Amomum species severely impacts efficient market processing and drug safety. ITS [10], ITS2 [9], randomly amplified polymorphic DNA [11], and single nucleotide polymorphisms (SNPs) [12] have been used to identify these three species, and some progress has been made. In recent years, the chloroplast (cp) genomes have shown great potential for species authentication, especially between closely related species [13,14], breeding, and phylogenetic analysis. Thus, developing cp genomic resources of these three Amomum species is not only beneficial in accurately identifying the closely related species, but it can also greatly contribute to the improvement of cp genetic engineering. In addition, up to now, only two Amomum cp genome sequences—A. kravanh and A. compactum—have been reported. So, the availability of complete cp genome sequences of these three Amomum species will be helpful in revealing the evolutionary relationships and phylogenetic position of species in commelinids.
Cp is an organelle specialized for providing essential energy for growth and reproduction of plants, which converts solar energy into chemical energy and releases oxygen [15]. A typical angiosperm cp genome consists of two inverted repeats (IRs) which are separated by a large single copy (LSC) and a small single copy (SSC), and encodes key proteins for photosynthesis and other necessary metabolic processes for plants in response to environmental stresses, such as drought, salt, heat, and light [16]. Cp genomes generally encode 120–130 genes with lengths of 120–170 kb [17]. Compared with the nuclear genome, the cp genome is smaller and employs multiple copies that can improve the expression level of the target gene [18]. Cp genome sequences are known for their highly conserved gene order and content. However, large-scale genome rearrangements, gene transfers, and gene insertions and losses have been identified in previous studies [19,20,21]. Over the past years, cp genomes have been shown to be an efficient tool to reveal phylogenetic relationships [20], identify the related species as a super-barcode [13,14], and develop cp genetic engineering [22]. With the development of sequencing technologies, the number of cp genome sequences recorded in the National Center for Biotechnology Information (NCBI) has increased dramatically. However, for the whole plant community, the total number of plants that have determined cp genome sequences is still insufficient. As far as medicinal plants are concerned, according to the national survey on Chinese material medical resources, there are more than 11,146 known medicinal plants in China alone. However, up to now, approximately 3300 cp genome sequences of plants have been recorded in NCBI, and this proportion is still small. In order to further develop cp genetic resources, more plant cp genome data need to be developed.
In this study, eight complete cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligular were determined. Differences in their essential characteristics and repeat sequences were revealed. Intraspecific and interspecific comparative analyses among the Amomum genus were conducted to discover highly divergent regions for species authentication. Thirty-two samples of eight Amomum species were collected to assess the species discriminatory power for these highly divergent regions. Furthermore, a phylogenetic tree was constructed to identify Amomum species and reveal their phylogenetic positions.

2. Results and Discussion

2.1. Genome Length and Features

A total of eight complete cp genomes of Amomum species were determined and submitted to GenBank (accession number: MH161416-18 and MN067431-35), which comprised four (Sample numbers: Y17085, Y17089, Y19017, and Y19021), two (Sample numbers: Y17088 and Y19018), and two (Sample numbers: Y19019 and Y19020) individual plants of A. villosum, A. villosum var. xanthioides and A. longiligulare, respectively. All three species studied here included individuals from the main producing areas in China [23]. In total, 34,775,994–49,499,932 raw reads were generated, and 257,464–780,870 reads were finally assembled to generate complete cp genomes with 232–714× sequencing depth. The complete cp genomes of the eight species from A. villosum, A. villosum var. xanthioides and A. longiligulare were determined to be 163,608-164,069 bp in size (Table 1). All displayed a typical quadripartite structure with a pair of IR regions (29,820–29,959 bp) that separated an LSC region (88,680–88,857 bp) from an SSC region (15,288–15,369 bp) (Table 1 and Figure 2). The overall GC content of the three species was 36.0–36.1% (Table 1). The GC content was unevenly distributed in the cp genome of these species. The GC content was the highest (41.0%) in the IR regions, the lowest (29.9–30.1%) in the SSC region, and approximately 33.7% in the LSC region (Table 1). These values were similar to most other reported cp genomes [19,24]. The high GC content of the IR region could be due to the four ribosomal RNA (rRNA) genes with a reduced number of duplicated AT nucleotides [25] and may be one of the important factors that cause the IR region to be more conservative than the LSC and SSC regions [26]. This phenomenon is also evident in many other angiosperms, such as A. kravanh [27], A. compactum [28], Scutellaria baicalensis [25] and Schisandra chinensis [29]. The AT representation at the third codon position (71.2~71.3%) was higher than that at the first (55.4%) and second (62.6%) positions in the protein-coding regions (CDS) of these species (Table 1). This bias was used as one principle to discriminate cp DNA from nuclear and mitochondrial DNA.
The complete cp genomes of these three species encode 113 different genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes (Table 2). A total of 20 genes were duplicated in the IR regions, resulting in a total of 133 genes located in the complete cp genomes (Table 2). The gene distribution in these three cp genomes was exactly the same: the LSC regions encoded 58 protein-coding genes and 23 tRNA genes, and the SSC regions contained 11 protein-coding genes and one tRNA gene. Moreover, eight protein-coding genes, eight tRNA genes and all four rRNA genes were duplicated in the IR regions.
In eukaryotic and semi-prokaryotic systems, gene expression occurs in nuclear-cytosolic and organelle compartments, respectively. Introns contribute greatly to the regulation of gene expression, and previous studies proven that intron can improve expression levels of exogenous genes in eukaryotic genomes [30]. Introns can accumulate more mutations than exons, and they maybe contain “old code”, which is the part of a gene that has lost its function during evolution [31,32,33]. A total of 18 intron-containing genes were found in the cp genomes of these three species, 15 of which contained one intron, and two (ycf3 and clpP) contained two introns (Figure 1, Table 2 and Table S1). The rps12 gene was a special trans-spliced gene with the 5′ end located in the LSC region and the duplicated 3′ ends in the IR regions. Among the 18 intron-containing genes, 12 were present in LSC regions, only one was present in the SSC region, and five were duplicated in the IR regions. The matK gene was contained in the intron of trnK-UUU gene, which showed the largest intron with more than 2500 bp.

2.2. Codon Usage

Codon usage patterns of coding sequences for A. villosum, A. villosum var. xanthioides and A. longiligulare were calculated based on the relative synonymous codon usage (RSCU) value [34]. We defined that the codon whose RSCU values were >1.00 was used more frequently, and vice versa. All protein-coding genes of A. villosum were encoded by 27,730 codons, A. longiligulare was encoded by 27,720 codons, and that of A. villosum var. xanthioides was 27,726–27,732 codons (Table S2). Like most angiosperms, leucine is the most prevalent amino acid in the cp genomes of A. villosum (2853, 10.3%), A. villosum var. xanthioides (2844–2851, 10.3%), and A. longiligulare (2855, 10.3%). Conversely, cysteine, with 316–317 (1.1%) codons, was the least frequent amino acid in the cp genomes of these three species. Except for trnL-CAA encoded by UUG, amino acid codon (RSCU > 1) in the cp genomes of three species preferentially showed A- or U-endings, which corresponded to the mentioned results that were calculated based on the cp genome sequences.

2.3. SSRs Analyses and Repeat Structures

Simple sequence repeats (SSRs), also called microsatellites, are short and tandem repeat DNA sequences with sizes of 1–6 bp [35] that are widely distributed throughout the cp genome and are usually used as important molecular markers for species authentication [36,37]. Here, we analyzed the distribution and the type of SSRs contained in the cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare. A total of 157 SSRs were identified in the whole cp genome of A. villosum, including 91 mono-, 33 di-, 7 tri-, 18 tetra-, three penta-, and five hexa-nucleotide SSRs (Table S3 and Figure 3). Among these SSRs, 29 (18.47%) were contained in the CDS (Table S4). A total of 154 (Y17088) or 158 (Y19018) SSRs were present in the whole cp genomes of A. villosum var. xanthioides and included 89 mono-, 31 (Y17088) or 33 (Y19018) di-, eight tri-, 18 tetra-, three (Y17088) or five (Y19018) penta-, and five hexa-nucleotide SSRs (Table S3 and Figure 3). The total number of identified SSRs in the CDS was 29 (Y17088) or 31 (Y19018) (18.83%) (Table S4). For A. longiligulare, there were 152 SSRs that comprised 90 mono-, 33 di-, six tri-, 16 tetra-, five penta-, and two hexa-nucleotide SSRs (Table S3 and Figure 3). Our results, show in Table S3, were consistent with the finding [19,25] that SSRs appeared more frequently in the LSC regions than in the SSC and IR regions. The overwhelming majority of mononucleotide repeats were composed of A or T (Table S3). All SSR types are listed in Table S3. Among all SSR types, A and T were always the most frequently used bases, which was in line with previous findings [24,27] that cp SSRs are generally composed of short polyadenine or polythymine repeats, resulting from the bias toward A and T of cp genomes. Cp SSRs have been widely used as molecular markers for the population genetic structure and phylogeographic study of some species [38,39,40] due to high substitution rates. And cp SSRs markers can be applied as complementary tools of nuclear markers to determine species identification and genetic relationships among closely related species [41,42]. Therefore, the availability of cp SSRs is of great benefit to develop useful molecular markers for the studies of genetic diversity, population structure, evolutionary studies, molecular identification and other further investigations of Amomum species.
Except for SSRs, some repeat structures with a length ≥30 bp are known as long repeat sequences, including forward repeats (F), palindromic repeats (P), reverse repeats (R), and complement repeats (C). These repeat structures promote the rearrangement of the cp genome and increase the population’s genetic diversity [29]. We analyzed the repeat structures of these three Amomum species. A total of 74 long repeats were present in the cp genome of A. villosum, as follows: 29 forward repeats, 39 palindromic repeats, five reverse repeats, and 1 complement repeat (Figure 4). A total of 63 (Y17088) or 68 (Y19018) long repeats were identified in the cp genome of A. villosum var. xanthioides including 25 (Y17088) or 24 (Y19018) forward repeats, 36 (Y17088) or 32 (Y19018) palindromic repeats, two (Y17088) or 11 (Y19018) reverse repeats and zero (Y17088) or one (Y19018) compliment repeats. As for the cp genomes of A. longiligulare, there were 53 long repeats, including 16 forward repeats, 29 palindromic repeats, five reverse repeats, and three complement repeats. In all three species, the majority of these repeats showed lengths between 30 and 39 bp. Only A. villosum var. xanthioides (Y17088) exhibited eight long repeats more than 100 bp in size (Figure 4).

2.4. Interspecific Comparison

From the above analysis, we can see that intraspecific differences between the cp genomes of four individuals of A. villosum, two individuals of A. villosum var. xanthioides and two individuals of A. longiligulare were very small. Here, in each of the three species, one of the cp genomes was used for interspecies comparative analyses within Amomum species A. villosum (Y17085), A. villosum var. xanthioides (Y17088), and A. longiligulare (Y19019).
Interspecific comparisons between five Amomum species (A. villosum, A. villosum var. xanthioides, A. longiligulare, A. krervanh and A. compactum) were conducted using mVISTA software with the annotated cp genome of A. villosum as a reference (Figure 5). Among the five Amomum species, A. villosum exhibited the longest cp genome size and was closely followed by A. villosum var. xanthioides and A. longiligulare, and A. krervanh was the smallest one with the size of 162,766 bp [27]. The common features of most angiosperms in the IR regions were more conserved than in the LSC and SSC regions, and the four rRNA (rrn4.5, rrn5, rrn16, and rrn23) genes were the most conserved in the two IR regions. The non-coding regions were more divergent than the coding regions. The cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare showed little difference from one another. Comparing A. krervanh and A. compactum, the highly divergent regions among these five cp genomes mainly occurred in the intergenic regions, including atpH-atpI, trnD-trnY, accD-psaI, ycf4-cemA and trnI-ycf2, which may be potential molecular markers for species authentication.
DNA polymorphism analyses were executed to detect highly variable regions and to show divergence at the sequence level in the cp genomes of the five Amomum species. The average value of nucleotide variability (PI) of all the five Amomum species was 0.00306. Results shown in Figure 6 also suggest that IR regions were less divergent than the LSC and SSC regions. A total of four mutational hotspots that showed remarkably high values of PI (≥0.014) were observed. As shown in Figure 6, mutational hotspots within these Amomum species were commonly located in the LSC and the SSC regions, which was in line with the result from mVISTA.

2.5. Species Authentication Analyses Based on cp Highly Divergent Regions

In order to assess the species discriminatory power of the highly divergent regions (atpH-atpI, trnD-trnY, accD-psaI, ycf4-cemA, and trnI-ycf2), a total of 32 samples from eight Amomum species (A. villosum, A. villosum var. xanthioides, A. longiligulare, A. chinense, A. compactum, A. tsaoko, A. koenigii, and A. maximum) were collected. Information about sample numbers and locations are listed in Table S5. To develop identification markers, specific primers were designed against conserved regions of atpH-atpI, trnD-trnY, accD-psaI, ycf4-cemA, and trnI-ycf2 (Table S6). All highly divergent regions except atpH-atpI were amplified successfully (Figure S1). For atpH-atpI, we designed two pairs of specific primers, but amplification was still unsuccessful. Amplified fragments from all tested Amomum samples were sequenced, and high-quality sequences were used to construct an NJ phylogenetic tree (Figure S2) to exhibit the species discriminatory power of the four highly divergent regions. The results showed that trnD-trnY and trnI-ycf2 could not distinguish these eight Amomum species. accD-psaI could differentiate A. villosum, A. villosum var. xanthioides and A. longiligulare from their closely related species- but could not distinguish these three species from each other. The ycf4-cemA could not only distinguish these three species from the other related species, but it could also distinguish A. longiligulare from A. villosum and A. villosum var. xanthioides. The NJ trees showed that accD-psaI and ycf4-cemA were potential markers that could be used to distinguish the three species from their closely related species, but they all had relatively low support rates. All four highly divergent regions could not provide valuable discriminatory power to distinguish these species from each other. Although there are many studies [19,25,43,44,45] using cp genomes to screen regions with high variations which may be used as molecular markers for species authentication, the efficiency and availability of identification remains to be determined.

2.6. Phylogenetic Analysis

To discuss the discriminatory power of cp genomes and determine the phylogenetic positions and evolutionary relationships of Amomum species among angiosperm, the maximum likelihood (ML) tree (Figure 7) of 17 representatives from order Zingiberales based on 64 common protein-coding genes, was constructed with Nicotiana tabacum and Salvia miltiorrhiza as outgroups. On the basis of the topologic structure, each of the four selected families (Zingiberales, Commelinales, Arecales and Poales) was clustered into a monophyletic branch respectively. All five Amomum species formed a strongly supported monophyletic clade sister to Alpinia oxyphylla, which is sometimes used as an adulterant of amomi fructus. This topological structure thereby hinted at a close relationship between these taxa. In addition, among the three Amomum species studied here, branches of four A. villosum and two A. villosum var. xanthioides clustered into a clade, thereby reflecting a closer relationship of these species. Two A. longiligulare clustered into a monophyletic branch separately from A. villosum and A. villosum var. xanthioides with high support rates. This showed that common protein-coding genes based on cp genomes could be considered an efficient resource for species authentication in Amomum. Results were consistent with traditional phylogenetic theory and provide a valuable basis for the establishment of phylogenetic relationships of species in Zingiberales.

3. Materials and Methods

3.1. Plant Material, DNA Extraction, and Sequencing

Fresh leaves of 8 individuals of A. villosum, A. villosum var. xanthioides and A. longiligulare were collected from their main producing areas. A total of four samples of A. villosum were collected respectively, from Yangchun City, Guangdong Province; Nanning City, Guangxi Province; Haikou City, Hainan Province; and Xishuangbanna City, Yunnan Province. Two samples of A. villosum var. xanthioides were collected from Xishuangbanna City, Yunnan Province; and two samples of A. longiligulare were collected from Haikou City, Hainan Province and Guangzhou City, Guangdong Province. A total of 32 dry leaf samples of eight Amomum species (A. villosum, A. villosum var. xanthioides, A. longiligulare, A. chinense, A. compactum, A. tsaoko, A. koenigii, and A. maximum) were collected for verification. The origin information of all samples is listed in Table S4. All samples and vouchers were stored in the herbarium of the Institute of Medicinal Plant Development (IMPLAD), Chinese Academy of Medical Sciences, and Peking Union Medical College (PUMC), and they were identified by Professor Yu-lin Lin. Leaf samples were cleansed by 70% ethanol, and total DNA was extracted using a DNeasy Plant Mini Kit (Qiagen Co., Hilden, Germany) using the standard protocol. DNA concentration and quality were respectively assessed through Nanodrop 2000C spectrophotometry and electrophoresis in 1% (w/v) agarose gel. The pure DNA from dry leaves was used to perform polymerase chain reaction (PCR) and Sanger sequencing. The pure DNA from fresh leaves was used to construct a shotgun library with an average insert size of 500 bp and carry out next generation sequencing (NGS) by an Illumina HiSeq X 10 platform in accordance with standard protocol.

3.2. Cp Genome Assembly

Raw reads obtained from NGS were first checked by FastQC and trimmed using Trimmomatic software [27,46]. The trimmed reads were composed of the data from nuclear and organelle genomes [47]. In order to extract cp-like reads, we mapped the trimmed reads to the reference database, which was constructed by all cp genomes recorded in NCBI, on the basis of their coverage and similarity. The extracted reads were then assembled into contigs by SOAPdenovo [48]. Scaffolds of the cp genome were constructed using SSPACE [49] and gaps were filled using GapCloser [50]. The assembly strategy of raw reads sequenced by Illumina HiSeq X was based on the method of Zhou et al. [19]. The accuracy of assembly of four junctions between the IRs and SSC/LSC regions was verified by PCR amplification and Sanger sequencing using newly designed primers listed in Table S7 and Figure S3.

3.3. Cp genome Annotation and Structure Analysis

The complete cp genome was annotated using CPGAVAS [51] and DOGMA [52] with default settings and was checked manually. The tRNA genes were identified by tRNAscan-SE [53]. A BLAST search was used to annotate boundaries of genes, introns/exons and coding regions versus reference sequences. A circular cp genome map was drawn using OGDRAW v1.2 (Organellar Genome DRAW) [54].
The GC content of the cp genome was analyzed using MEGA 6.0 [55]. The distribution of codon usage was investigated using the software CodonW with RSCU ratios. REPuter [56] was used to identify the size and location of repeat sequences, including forward (F), palindromic (P), reverse (R), and complement (C) repeats in the cp genomes. The minimal size for all repeat types was 30 bp, and the two repeat copies had at least 90% similarity. MISA software [57] was used to detect SSRs.

3.4. Cp genome Comparisons and Species Authentication

To detect variations within the Amomum cp genomes, we compared the cp genomes of A. villosum, A. villosum var. xanthioides, A. longiligulare and two other Amomum species, namely, A. kravanh and A. compactum, which were downloaded from GenBank by mVISTA [58]. All cp genome sequences were aligned by MAFFT software [59], and nucleotide diversity of the cp genomes was analyzed through the sliding window by DNA Sequence Polymorphism (DnaSP) [60] software. The step size was set to 200 bp with an 800-bp window length.
In order to evaluate the species discriminatory power for the highly divergent regions, we designed the specific primers for each region (Table S5) and carried out PCR amplification based on some Amomum species. PCR amplification was performed in 25-μL reaction mixtures containing 20 ng to ~100 ng of genomic DNA template, 12.5 μL of 2 × EasyTaq PCR SuperMix, and 1 μL of each primer (2.5μM). Purified PCR products were sequenced using the specific primers in the ABI3730XL sequencer. Proofreading and contig assembly of sequencing peak diagrams were performed using CodonCode Aligner 3.7.1 and the sequences were aligned using ClustalW. To evaluate the capability of these sequences to authenticate the studied Amomum species, a phylogenetic tree was constructed using the neighbor-joining algorithm (NJ tree), and bootstrap tests were conducted using 1000 resamples to assess the confidence of the phylogenetic relationships using MEGA 6.0.

3.5. Phylogenetic Analyses

To discuss the discriminatory power of the cp genomes and reveal the phylogenetic position of Amomum species, we downloaded 35 complete cp genome sequences from NCBI’s Organelle Genome and Nucleotide Resources database (Table S8). The sequences of 64 common protein-coding genes (atpA, atpB, atpE, atpF, atpH, atpI, ccsA, cemA, clpP, matK, ndhA, ndhB, ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK, petA, petG, petL, petN, psaA, psaB, psaC, psaI, psaJ, psbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, rbcL, rpl14, rpl2, rpl22, rpl32, rpl33, rpl36, rpoA, rpoB, rpoC1, rps11, rps16, rps19, rps2, rps3, rps4, rps7, rps8, ycf3, and ycf4) shared in these species were extracted and aligned separately using MAFFT v7 [61], and the alignment was manually verified. Phylogenetic trees were reconstructed based on the 64 concatenated protein-coding gene sequences by ML methods. ML analysis was conducted with a bootstrap of 1000 repetitions based on the GTR+F+R3 nucleotide substitution model using IQ-TREE software [62]. This adopted best-fit model was also determined by IQ-TREE. N. tabacum and S. miltiorrhiza were set as outgroups.

4. Conclusions

This is the first study to sequence and determine the complete cp genomes of eight individual species of A. villosum, A. villosum var. xanthioides and A. longiligulare via high-throughput sequencing, which were collected from the different main producing areas in China. First of all, we obtained and compared eight complete cp genomes, and we found that the cp genomes of the same species from different producing areas showed little difference. Then, five divergent regions (atpH-atpI, trnD-trnY, accD-psaI, ycf4-cemA and trnI-ycf2), which were detected by comparative analyses using Sanger sequencing based on 32 samples, proved to not be useful as effective molecular markers to identify species in Amomum. Further, phylogenetic analyses revealed that cp genomes could be used to identify Amomum species. The availability of cp genome sequences is expected to improve species identification and phylogenetic analyses in Amomum.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/16/4040/s1.

Author Contributions

Conceptualization, H.Y.; Data curation, Y.C., X.C. and L.N.; Formal analysis, Y.C. and W.S.; Funding acquisition, H.L., J.S. and H.Y.; Investigation, Y.L.; Methodology, Y.C., X.C. and L.N.; Project administration, H.Y.; Resources, Y.L., H.L. and X.Z.; Software, Y.C. and H.H.; Supervision, H.Y.; Validation, X.C.; Writing—original draft, Y.C.; Writing—review & editing, J.S. and H.Y.

Funding

This research was funded by National Key R&D Program of China (2017YF C1701102), Chinese Academy of Medical Sciences (CAMS) Innovation Fund for Medical Sciences (CIFMS) (NO. 2016-I2M-3-016), and Yunnan Provincal Science and Technology Major Projects (2018ZF011).

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

LSCLarge single copy
SSCSmall single copy
IRInverted repeat
MLMaximum likelihood
SSRSimple sequence repeats
ATPAdenosine triphosphate
NADHNicotinamide adenine dinucleotide

References

  1. Duan, L.S.; Zhang, L.X.; Peng, J.M.; Jie, M.A. Original Report on Investigation of Ammomum Villosum Germplasm Resources in Xishuangbanna. Lishizhen Med. Mater. Med. Res. 2009, 20, 627–628. [Google Scholar]
  2. Commission, C.P. Pharmacopoeia of the Peoples Republic of China; China Medical Science Press, 2015. Available online: https://www.tsoshop.co.uk/Medical/Pharmacopoeia/Chinese-Pharmacopoeia/ (accessed on 18 August 2019).
  3. Chen, Z.; Ni, W.; Yang, C.; Zhang, T.; Lu, S.; Zhao, R.; Mao, X.; Yu, J. Therapeutic Effect of Amomum villosum on Inflammatory Bowel Disease in Rats. Front. Pharmacol. 2018, 9, 639. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lu, C.L.; Zhao, H.Y.; Jiang, J.G. Evaluation of multi-activities of 14 edible species from Zingiberaceae. Int. J. Food Sci. Nutr. 2013, 64, 28–35. [Google Scholar] [CrossRef] [PubMed]
  5. Wu, X.; Li, X.; Xiao, F.; Zhang, Z.; Xu, Z.; Wang, H. Studies on the analgesic and anti-inflammatory effect of bornyl acetate in volatile oil from Amomum villosum. J. Chin. Med. Mater. 2004, 27, 438. [Google Scholar]
  6. Tang, L.; He, G.; Su, J.; Xu, H. The Strategy to Promote the Development of Industry of Genuine Medicinal Material of Amomum villosum. Chin. Agric. Sci. Bull. 2012, 28, 94–99. [Google Scholar]
  7. Quan, Q.; Lin, J.C.; Liang, J.M.; Hu, Z.L.; Huang, J.X. Compare on the Content of Bronyl Acetate and Total Volatile Oil in Fructus Amomi from Different Producing Area. Guid. J. Tradit. Chin. Med. Pharm. 2017, 23, 70–73. [Google Scholar]
  8. Chen, J.; Ding, P.; Xu, X.; Xu, H. A resource investigation and commodity identification of Fructus Amomi. J. Chin. Med. Mater. 2001, 24, 18. [Google Scholar]
  9. Han, J.P.; Li, M.N.; Shi, L.C. Identification of Amomi fructus and its adulterants based on ITS2 sequences. Glob. Tradit. Chin. Med. 2011, 4, 99–102. [Google Scholar]
  10. Zhou, L.; Wang, P.; Huang, F.; Cao, L.; Liang, R. ITS sequence analysis of Amomum villosum. Chin. Tradit. Herb. Drugs 2002, 33, 72–75. [Google Scholar]
  11. Wang, P.; Huang, F.; Zhou, L.; Cao, L.; Liang, S.; Xu, H.; Liu, J. Analysis of Amomun villosum species and some adulterants of zingiberaceae by RAPD. J. Chin. Med. Mater. 2000, 23, 71–74. [Google Scholar]
  12. Huang, Q.; Duan, Z.; Yang, J.; Ma, X.; Zhan, R.; Xu, H.; Chen, W. SNP typing for germplasm identification of Amomum villosum Lour. Based on DNA barcoding markers. PLoS ONE 2014, 9, e114940. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, X.; Zhou, J.; Cui, Y.; Wang, Y.; Duan, B.; Yao, H. Identification of Ligularia Herbs Using the Complete Chloroplast Genome as a Super-Barcode. Front. Pharmacol. 2018, 9, 695. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Park, I.; Yang, S.; Kim, W.; Noh, P.; Lee, H.; Moon, B. Authentication of Herbal Medicines Dipsacus asper and Phlomoides umbrosa Using DNA Barcodes, Chloroplast Genome, and Sequence Characterized Amplified Region (SCAR) Marker. Molecules 2018, 23, 1748. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, J.; Qi, Z.C.; Zhao, Y.P.; Fu, C.X.; Jenny Xiang, Q.Y. Complete cpDNA genome sequence of Smilax china and phylogenetic placement of Liliales--influences of gene partitions and taxon sampling. Mol. Phylogenetics Evol. 2012, 64, 545–562. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, X.; Li, Y.; Yang, H.; Zhou, B. Chloroplast Genome of the Folk Medicine and Vegetable Plant Talinum paniculatum (Jacq.) Gaertn.: Gene Organization, Comparative and Phylogenetic Analysis. Molecules 2018, 23, 857. [Google Scholar]
  17. Green, B.R. Chloroplast genomes of photosynthetic eukaryotes. Plant J. 2011, 66, 34–44. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, H.-H.; Yin, W.-B.; Hu, Z.-M. Advances in chloroplast engineering. J. Genet. Genom. 2009, 36, 387–398. [Google Scholar] [CrossRef]
  19. Zhou, J.; Chen, X.; Cui, Y.; Sun, W.; Li, Y.; Wang, Y.; Song, J.; Yao, H. Molecular Structure and Phylogenetic Analyses of Complete Chloroplast Genomes of Two Aristolochia Medicinal Species. Int. J. Mol. Sci. 2017, 18, 1839. [Google Scholar] [CrossRef] [PubMed]
  20. Jansen, R.K.; Cai, Z.; Raubeson, L.A.; Daniell, H.; Depamphilis, C.W.; Leebensmack, J.; Müller, K.F.; Guisingerbellian, M.; Haberle, R.C.; Hansen, A.K. Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns. Proc. Natl. Acad. Sci. USA 2007, 104, 19369–19374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. He, L.; Qian, J.; Li, X.; Sun, Z.; Xu, X.; Chen, S. Complete Chloroplast Genome of Medicinal Plant Lonicera japonica: Genome Rearrangement, Intron Gain and Loss, and Implications for Phylogenetic Studies. Molecules 2017, 22, 249. [Google Scholar] [CrossRef]
  22. Daniell, H.; Lin, C.S.; Ming, Y.; Chang, W.J. Chloroplast genomes: Diversity, evolution, and applications in genetic engineering. Genome Biol. 2016, 17, 134. [Google Scholar] [CrossRef] [PubMed]
  23. Peter, H.R.; Zhang, L.B.; Ihsan, A.; Nicholas, J.; Zhu, G.H. Flora of China; Science Press: Beijing, China; Missouri Botanical Garden Press: St. Louis, MO, USA, 2000; Volume 24, pp. 347–356. [Google Scholar]
  24. Xiang, B.; Li, X.; Qian, J.; Wang, L.; Ma, L.; Tian, X.; Wang, Y. The Complete Chloroplast Genome Sequence of the Medicinal Plant Swertia mussotii Using the PacBio RS II Platform. Molecules 2016, 21, 1029. [Google Scholar] [CrossRef] [PubMed]
  25. Jiang, D.; Zhao, Z.; Zhang, T.; Zhong, W.; Liu, C.; Yuan, Q.; Huang, L. The Chloroplast Genome Sequence of Scutellaria baicalensis Provides Insight into Intraspecific and Interspecific Chloroplast Genome Diversity in Scutellaria. Genes 2017, 8, 227. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, Y.; Dang, Y.; Li, Q.; Lu, J.; Li, X.; Wang, Y. Complete Chloroplast Genome Sequence of Poisonous and Medicinal Plant Datura stramonium: Organizations and Implications for Genetic Engineering. PLoS ONE 2014, 9, e110656. [Google Scholar] [CrossRef] [PubMed]
  27. Wu, M.; Li, Q.; Hu, Z.; Li, X.; Chen, S. The Complete Amomum kravanh Chloroplast Genome Sequence and Phylogenetic Analysis of the Commelinids. Molecules 2017, 22, 1875. [Google Scholar] [CrossRef] [PubMed]
  28. Wu, M.L.; Li, Q.; Xu, J.; Li, X.W. Complete chloroplast genome of the medicinal plant Amomum compactum: Gene organization, comparative analysis, and phylogenetic relationships within Zingiberales. Chin. Med. 2018, 13, 10. [Google Scholar] [CrossRef]
  29. Guo, H.; Liu, J.; Luo, L.; Wei, X.; Zhang, J.; Qi, Y.; Zhang, B.; Liu, H.; Xiao, P. Complete chloroplast genome sequences of Schisandra chinensis: Genome structure, comparative analysis, and phylogenetic relationship of basal angiosperms. Sci. China Life Sci. 2017, 60, 1286–1290. [Google Scholar] [CrossRef]
  30. Xu, J.; Feng, D.; Song, G.; Wei, X.; Chen, L.; Wu, X.; Li, X.; Zhu, Z. The first intron of rice EPSP synthase enhances expression of foreign gene. Sci. China 2003, 46, 561–569. [Google Scholar] [CrossRef]
  31. Kelchner, S.A. Group II introns as phylogenetic tools: Structure, function, and evolutionary constraints. Am. J. Bot. 2002, 89, 1651–1669. [Google Scholar] [CrossRef]
  32. Mattick, J.S. Introns: Evolution and function. Curr. Opin. Genet. Dev. 1994, 4, 823–831. [Google Scholar] [CrossRef]
  33. Chorev, M.; Carmel, L. The Function of Introns. Front. Genet. 2012, 3, 55. [Google Scholar] [CrossRef] [Green Version]
  34. Butt, A.M.; Nasrullah, I.; Qamar, R.; Tong, Y. Evolution of codon usage in Zika virus genomes is host and vector specific. Emerg. Microbes Infect. 2016, 5, e107. [Google Scholar] [CrossRef]
  35. Powell, W.; Morgante, M.; Mcdevitt, R.; Vendramin, G.G.; Rafalski, J.A. Polymorphic simple sequence repeat regions in chloroplast genomes: Applications to the population genetics of pines. Proc. Natl. Acad. Sci. USA 1995, 92, 7759–7763. [Google Scholar] [CrossRef]
  36. Ai-Hong, Y.; Jin-Ju, Z.; Xiao-Hong, Y.; Hong-Wen, H. Chloroplast microsatellite markers in Liriodendron tulipifera (Magnoliaceae) and cross-species amplification in L. chinense. Am. J. Bot. 2011, 98, 123–126. [Google Scholar]
  37. Jiao, Y. Development of simple sequence repeat (SSR) markers from a genome survey of Chinese bayberry (Myrica rubra). Bmc Genom. 2012, 13, 201. [Google Scholar] [CrossRef]
  38. Echt, C.S.; Deverno, L.L.; Anzidei, M.; Vendramin, G.G. Chloroplast microsatellites reveal population genetic diversity in red pine, Pinus resinosa Ait. Mol. Ecol. 2010, 7, 307–316. [Google Scholar] [CrossRef]
  39. Navascués, M.; Vaxevanidou, Z.; Climent, J.; Gil, L.; Emerson, B.C. Chloroplast microsatellites reveal colonization and metapopulation dynamics in the Canary Island pine. Mol. Ecol. 2010, 15, 2691–2698. [Google Scholar] [CrossRef]
  40. Baneh, H.D.; Mohammadi, S.A.; Labra, M.; Nazemieh, A.; Mardi, M. Chloroplast microsatellites markers to assess genetic diversity in wild and cultivated grapevines of Iran. Pak. J. Biol. Sci. 2007, 10, 1855–1859. [Google Scholar]
  41. Park, H.; Kim, C.; Lee, Y.M.; Kim, J.H.J.A.i.P.S. Development of Chloroplast Microsatellite Markers for the Endangered Maianthemum bicolor (Asparagaceae s.l.). Appl. Plant. Sci. 2016, 4, apps.1600032. [Google Scholar] [CrossRef]
  42. Yan, Y.D.; Li, X.Y.; Worth, J.R.P.; Lin, X.Y.; Ruhsam, M.; Chen, L.; Wu, X.T.; Wang, M.Q.; Thomas, P.I.; Wen, Y.F. Development of chloroplast microsatellite markers for Glyptostrobus pensilis (Cupressaceae). Appl. Plant. Sci. 2019, 7, e11277. [Google Scholar] [CrossRef] [Green Version]
  43. Zhou, Y.; Nie, J.; Xiao, L.; Hu, Z.; Wang, B. Comparative Chloroplast Genome Analysis of Rhubarb Botanical Origins and the Development of Specific Identification Markers. Molecules 2018, 23, 2811. [Google Scholar] [CrossRef]
  44. Niu, Z.; Xue, Q.; Zhu, S.; Sun, J.; Liu, W.; Ding, X. The Complete Plastome Sequences of Four Orchid Species: Insights into the Evolution of the Orchidaceae and the Utility of Plastomic Mutational Hotspots. Front. Plant Sci. 2017, 8, 715. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, J.B.; Tang, M.; Li, H.T.; Zhang, Z.R.; Li, D.Z. Complete chloroplast genome of the genus Cymbidium: Lights into the species identification, phylogenetic implications and population genetic analyses. Bmc Evol. Biol. 2013, 13, 84. [Google Scholar] [CrossRef]
  46. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  47. Li, Y.; Zhou, J.-g.; Chen, X.-l.; Cui, Y.-x.; Xu, Z.-c.; Li, Y.-h.; Song, J.-y.; Duan, B.-z.; Yao, H. Gene losses and partial deletion of small single-copy regions of the chloroplast genomes of two hemiparasitic Taxillus species. Sci. Rep. 2017, 7, 12834. [Google Scholar] [CrossRef]
  48. Gogniashvili, M.; Naskidashvili, P.; Bedoshvili, D.; Kotorashvili, A.; Kotaria, N.; Beridze, T. Complete chloroplast DNA sequences of Zanduri wheat (Triticum spp.). Genet. Resour. Crop. Evol. 2015, 62, 1269–1277. [Google Scholar] [CrossRef]
  49. 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]
  50. Acemel, R.D.; Tena, J.J.; Irastorza-Azcarate, I.; Marletaz, F.; Gomez-Marin, C.; De, l.C.-M.E.; Bertrand, S.; Diaz, S.G.; Aldea, D.; Aury, J.M. A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation. Nat. Genet. 2016, 48, 336–341. [Google Scholar] [CrossRef] [Green Version]
  51. Liu, C.; Shi, L.; Zhu, Y.; Chen, H.; Zhang, J.; Lin, X.; Guan, X.J.B.G. 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]
  52. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [Green Version]
  53. Schattner, P.; Brooks, A.N.; Lowe, T.M. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005, 33, 686–689. [Google Scholar] [CrossRef]
  54. 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, 267–274. [Google Scholar] [CrossRef]
  55. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  56. 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]
  57. Yang, X.M.; Sun, J.T.; Xue, X.F.; Zhu, W.C.; Hong, X.Y. Development and characterization of 18 novel EST-SSRs from the western flower thrips, Frankliniella occidentalis (Pergande). Int. J. Mol. Sci. 2012, 13, 2863–2876. [Google Scholar] [CrossRef]
  58. Frazer, K.A.; Pachter, L.; Poliakov, A.; Rubin, E.M.; Dubchak, I. VISTA: Computational tools for comparative genomics. Nucleic Acids Res. 2004, 32, W273. [Google Scholar] [CrossRef]
  59. Nakamura, T.; Yamada, K.D.; Tomii, K.; Katoh, K. Parallelization of MAFFT for large-scale multiple sequence alignments. Bioinformatics 2018, 34, 2490–2492. [Google Scholar] [CrossRef] [Green Version]
  60. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef]
  61. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2017. [Google Scholar] [CrossRef]
  62. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; 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]
Ijms 20 04040 g001 550
Figure 1. Commercial dried samples and plant materials of A. villosum (A), A. villosum var. xanthioides (B) and A. longiligulare (C). In picture 1, differences in the size and color of ligules are highlighted in red boxes. In picture 2, fresh fruits were removed from soil and photographed. In picture 3, amomi fructus (dried fruits) were collected from herb markets and photographed.
Figure 1. Commercial dried samples and plant materials of A. villosum (A), A. villosum var. xanthioides (B) and A. longiligulare (C). In picture 1, differences in the size and color of ligules are highlighted in red boxes. In picture 2, fresh fruits were removed from soil and photographed. In picture 3, amomi fructus (dried fruits) were collected from herb markets and photographed.
Ijms 20 04040 g001
Ijms 20 04040 g002 550
Figure 2. Gene map of the complete cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare. Only one map is shown here because the differences among the three species are negligible, and the gene maps of the three cp genomes are almost identical. Genes of different functional groups are separated by color. Genes inside the circle are transcribed clockwise, whereas those on the outside are transcribed counter-clockwise. The dark grey area in the inner circle corresponds to GC content, whereas the light grey area corresponds to AT content.
Figure 2. Gene map of the complete cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare. Only one map is shown here because the differences among the three species are negligible, and the gene maps of the three cp genomes are almost identical. Genes of different functional groups are separated by color. Genes inside the circle are transcribed clockwise, whereas those on the outside are transcribed counter-clockwise. The dark grey area in the inner circle corresponds to GC content, whereas the light grey area corresponds to AT content.
Ijms 20 04040 g002
Ijms 20 04040 g003 550
Figure 3. Analysis of simple sequence repeats (SSRs) in the cp genomes of three Amomum species.
Figure 3. Analysis of simple sequence repeats (SSRs) in the cp genomes of three Amomum species.
Ijms 20 04040 g003
Ijms 20 04040 g004 550
Figure 4. Repeat sequences in the three cp genomes. F, P, R, and C indicate the repeat types F (forward), P (palindrome), R (reverse), and C (complement), respectively.
Figure 4. Repeat sequences in the three cp genomes. F, P, R, and C indicate the repeat types F (forward), P (palindrome), R (reverse), and C (complement), respectively.
Ijms 20 04040 g004
Ijms 20 04040 g005 550
Figure 5. Sequence identity plot comparing the five cp genomes in Amomum species with A. villosum as a reference by using mVISTA. Grey arrows and thick black lines above the alignment indicate genes with their orientation and the position of the inverted repeats (IRs), respectively. A cut-off of 70% identity was used for the plots, and the Y-scale represents the percent identity ranging from 50% to 100%.
Figure 5. Sequence identity plot comparing the five cp genomes in Amomum species with A. villosum as a reference by using mVISTA. Grey arrows and thick black lines above the alignment indicate genes with their orientation and the position of the inverted repeats (IRs), respectively. A cut-off of 70% identity was used for the plots, and the Y-scale represents the percent identity ranging from 50% to 100%.
Ijms 20 04040 g005
Ijms 20 04040 g006 550
Figure 6. Sliding window analysis based on the cp genomes of five Amomum species. Window length: 800 bp; step size: 200 bp. X-axis: position of the midpoint of a window. Y-axis: nucleotide diversity of each window.
Figure 6. Sliding window analysis based on the cp genomes of five Amomum species. Window length: 800 bp; step size: 200 bp. X-axis: position of the midpoint of a window. Y-axis: nucleotide diversity of each window.
Ijms 20 04040 g006
Ijms 20 04040 g007 550
Figure 7. Phylogenetic tree of A. villosum, A. villosum var. xanthioides and A. longiligulare inferred by maximum likelihood (ML) analyses based on 64 common protein-coding genes.
Figure 7. Phylogenetic tree of A. villosum, A. villosum var. xanthioides and A. longiligulare inferred by maximum likelihood (ML) analyses based on 64 common protein-coding genes.
Ijms 20 04040 g007
Table
Table 1. Comparisons among the cp genome characteristics of A. villosum, A. villosum var. xanthioides and A. longiligulare.
Table 1. Comparisons among the cp genome characteristics of A. villosum, A. villosum var. xanthioides and A. longiligulare.
TypeA. VillosumA. Villosum var. XanthioidesA. Longiligulare
Accession NumberMH161416
MH161418
MN067431
MN067432		MH161417
MN067433		MN067434
MN067435
Total Length (bp)164,068–164,069163,981–163,985163,608
LSC Length (bp)88,797–88,79888,720–88,85788,680
IR Length (bp)29,95929,886–29,94829,820
SSC Length (bp)135,352–135,35315,352–15,36915,288
CDS Length (bp)83,19083,178–83,19683,160
Total GC content36.0%36.0%36.1%
GC content of LSC33.7%33.7%33.7%
GC content of IRa41.1%41.1%41.1%
GC content of IRb41.1%41.1%41.1%
GC content of SSC30.0%30.0%30.1%
AT content at 1st position55.4%55.4%55.4%
AT content at 2nd position62.6%62.6%62.6%
AT content at 3rd position71.2%71.2%71.2%
Table
Table 2. Gene contents in the cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare.
Table 2. Gene contents in the cp genomes of A. villosum, A. villosum var. xanthioides and A. longiligulare.
Group of GenesGene NamesNumber of Genes
Photosystem IpsaA, psaB, psaC, psaI, psaJ5
Photosystem IIpsbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ15
Cytochrome b/f complexpetA, petB *, petD *, petG, petL, petN6
ATP synthaseatpA, atpB, atpE, atpF *, atpH, atpI6
NADH dehydrogenasendhA *, ndhB(×2) *, ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK12
RubisCO large subunitrbcL1
RNA polymeraserpoA, rpoB, rpoC1 *, rpoC24
Ribosomal proteins (SSU)rps2, rps3, rps4, rps7(×2), rps8, rps11, rps12(×2) **, rps14, rps15, rps16 *, rps18, rps19(×2)15
Ribosomal proteins (LSU)rpl2(×2) *, rpl14, rpl16 *, rpl20, rpl22, rpl23(×2), rpl32, rpl33, rpl3611
Proteins of unknown functionycf1(×2), ycf2(×2), ycf3 **, ycf46
Transfer RNAs38 tRNAs (8 in the IRs(×2), 6 contain one intron)38
Ribosomal RNAsrrn4.5(×2), rrn5(×2), rrn16(×2), rrn23(×2)8
Other genesaccD, clpP **, matK, ccsA, cemA, infA6
(×2) indicates the gene sequence is repeated twice. * indicates genes containing one intron; while ** indicates gene containing two introns. The rps12 gene is a trans-spliced gene.

Share and Cite

MDPI and ACS Style

Cui, Y.; Chen, X.; Nie, L.; Sun, W.; Hu, H.; Lin, Y.; Li, H.; Zheng, X.; Song, J.; Yao, H. Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species. Int. J. Mol. Sci. 2019, 20, 4040. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20164040

AMA Style

Cui Y, Chen X, Nie L, Sun W, Hu H, Lin Y, Li H, Zheng X, Song J, Yao H. Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species. International Journal of Molecular Sciences. 2019; 20(16):4040. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20164040

Chicago/Turabian Style

Cui, Yingxian, Xinlian Chen, Liping Nie, Wei Sun, Haoyu Hu, Yulin Lin, Haitao Li, Xilong Zheng, Jingyuan Song, and Hui Yao. 2019. "Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species" International Journal of Molecular Sciences 20, no. 16: 4040. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20164040

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