- Research article
- Open Access
Analysis of six chloroplast genomes provides insight into the evolution of Chrysosplenium (Saxifragaceae)
BMC Genomics volume 21, Article number: 621 (2020)
Chrysosplenium L. (Saxifragaceae) is a genus of plants widely distributed in Northern Hemisphere and usually found in moist, shaded valleys and mountain slopes. This genus is ideal for studying plant adaptation to low light conditions. Although some progress has been made in the systematics and biogeography of Chrysosplenium, its chloroplast genome evolution remains to be investigated.
To fill this gap, we sequenced the chloroplast genomes of six Chrysosplenium species and analyzed their genome structure, GC content, and nucleotide diversity. Moreover, we performed a phylogenetic analysis and calculated non-synonymous (Ka) /synonymous (Ks) substitution ratios using the combined protein-coding genes of 29 species within Saxifragales and two additional species as outgroups, as well as a pair-wise estimation for each gene within Chrysosplenium. Compared with the outgroups in Saxifragaceae, the six Chrysosplenium chloroplast genomes had lower GC contents; they also had conserved boundary regions and gene contents, as only the rpl32 gene was lost in four of the Chrysosplenium chloroplast genomes. Phylogenetic analyses suggested that the Chrysosplenium separated to two major clades (the opposite group and the alternate group). The selection pressure estimation (Ka/Ks ratios) of genes in the Chrysosplenium species showed that matK and ycf2 were subjected to positive selection.
This study provides genetic resources for exploring the phylogeny of Chrysosplenium and sheds light on plant adaptation to low light conditions. The lower average GC content and the lacking gene of rpl32 indicated selective pressure in their unique habitats. Different from results previously reported, our selective pressure estimation suggested that the genes related to photosynthesis (such as ycf2) were under positive selection at sites in the coding region.
Challenging environments may impose selective pressure on genes, which could leave a footprint of natural selection in genes involved in adaptation to the environment. The chloroplast genome is typically quadripartite in structure, containing a large single copy (LSC) and a small single copy (SSC) separated by a pair of inverted repeats (IR). It is widely used for chloroplast inheritance, domestication studies, phylogeny and adaptative evolution [1,2,3]. Adaptive evolution is considered as the improved adaptability of species for changing environmental conditions during their evolutionary processes. And it is driven by evolutionary processes such as natural selection, which act on genetic variations produced by genetic recombination, gene mutations, and gene flow . Strong purifying selection detected rather than expected positive selection in chloroplast genome of the green alga (Ostreobium quekettii) facilitated its extremely low light adaptation . This study contributed to our understanding of plant adaptive evolution. However, to our knowledge, a comparative analysis of chloroplast genomes of angiosperms with low light requirements has not been conducted.
Chrysosplenium L. is a genus of Saxifragaceae and belongs to the subfamily Saxifragoideae according to the APG IV . The genus plays an important role in the phylogeny of Saxifragaceae, comprises about 79 perennial herbs , and mainly occurs in the northern hemisphere, with the highest species diversity in East Asia; only two species, Chrysosplenium valdivicum Hook and Chrysosplenium macranthum Hook, are found in the Southern Hemisphere [8, 9]. Thirty-six species and fifteen variants have been recorded in China [10, 11]. Chrysosplenium is divided into two subgenera according to the leaf arrangement: Alternifolia Franchet with alternate leaves and Oppositifolia Franchet with opposite leaves [8, 10, 12]. The genus was regarded as a typical group with floristic disjunction and is important for studying speciation [13, 14]. Chrysosplenium (Saxifragaceae) is usually found in shaded valleys and mountain slopes and, compared with other genera within Saxifragaceae, this genus has the lowest light requirement [8, 15]. Therefore, examination of the chloroplast genomes of Chrysosplenium species may provide insight into the impacts of low light in angiosperms.
Nuclear markers (such as the internal transcribed spacer, ITS) of the ribosomal DNA, and chloroplast markers were employed to determine the molecular phylogeny of Chrysosplenium [16, 17]. Compared to nuclear markers, the chloroplast genome possesses highly conserved DNA sequences and a lower substitution level (especially in inverted repeat regions). Therefore, the chloroplast genome is ideal for phylogenetic inference at the species and higher levels [9, 17,18,19,20].
In this study, we aimed to provide a comprehensive insight into the evolution of the chloroplast genomes of several Chrysosplenium species. First, we sequenced the chloroplast genomes of six Chrysosplenium species in addition to a previous study . Next, we conducted comparative chloroplast genome analyses for these six genomes, plus four Saxifragaceae chloroplast genomes from GenBank. Then, we constructed a phylogeny of Chrysosplenium using chloroplast genomes of 29 species within Saxifragales and two outgroups. Finally, we estimated selective pressures to investigate whether the genes related to photosynthesis in Chrysosplenium are under purifying selection or positive selection.
Organization of the Chloroplast Genomes of Chrysosplenium species
The chloroplast genomes of the Chrysosplenium species contain the typical quadripartite structures (Fig. 1), which include a large single copy (LSC), two inverted repeats (IR) and a small single copy (SSC) region. The seven Chrysosplenium chloroplast genomes ranged from 151,679 bp to 153,460 bp in length (see Table 1 for details), with IRs 25,974–26,224 bps, LSCs 82,771–83,752 bps and SSCs 16,960–17,342 bps. Each Chrysosplenium chloroplast genome encoded 30 transfer RNAs (tRNAs) and 4 ribosomal RNAs (rRNAs). Each genome also includes 79 functional proteins (Table 2) encoding genes except for C. macrophyllum, C. flagelliferum, C. alternifolium, C. ramosum, which lacked rpl32. By homolog and expression analysis of chloroplast rpl32 in C. sinicum (Cp_rpl32), another homolog (Nu_rpl32) was also identified nuclear genome of C. sinicum, and the expression value (4.86 FPKMs) of Nu_rpl32 is much lower than that (20,835.9 FPKMs) of Cp_rpl32 (Additional file 1: Supplementary Figure S1). In total, each chloroplast genome includes 113 (rpl32 present) or 112 (rpl32 loss) genes. The rps12 gene in Chrysosplenium was recognized as a trans-spliced gene, with the first exon located in the LSC region and the other one or two exons distributed in the IR regions. In addition, 17 intron-containing genes were also detected (Additional file 2: Supplementary Table S4). The chloroplast genome size and gene content neither significantly diverged between Oppositifolia and Alternifolia subgenera (Table 1) nor significantly diverged between Chrysosplenium and other genera of Saxifragaceae.
GC content, nucleotide diversity, and repeat analysis
When we compared the total GC content of the chloroplast genomes of Chrysosplenium species with that of the chloroplast genomes of the three non-Chrysosplenium Saxifragaceae species (S. stolonifera, B. scopulosa, and O. rupifraga), we found the Chrysosplenium species have the lowest total GC contents (< 37.5%) (Fig. 2 and Additional file 2: Supplementary Table S5). In addition, Chrysosplenium has the lowest GC contents (< 29.7%) at the third codon position (GC3). Within the Chrysosplenium species, the GC contents in subgenus Oppositifolia were slightly lower than those in subgenus Alternifolia, regardless of the total GC contents or those in GC3.
The IR regions were more conserved than the LSC and SSC regions, with average Pi values of 0.00586 in IR regions, 0.01760 in the LSC region, and 0.01900 in the SSC region (Additional file 2: Supplementary Table S6 and Additional file 3: Supplementary Figure S2). In the LSC region, psbT has the highest Pi value of 0.22159, followed by trnG-GCC with Pi value of 0.10369.
Among the mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide categories of SSRs in the chloroplast genomes of the Chrysosplenium species, mono-nucleotide repeats were the most common (Additional File 4: Supplementary Table S7 and Additional File 5: Supplementary Figure S3A) ranging from 42.42% (C. sinicum) to 61.29% (C. flagelliferum). Hexa-nucleotide repeats account for the lowest proportion of SSRs in C. ramosum, C. sinicum, and C. flagelliferum. Chrysosplenium species contained fewer SSRs than B. scopulosa and O. rupifraga. Among the four repeat types, the most common repeat type was palindromic repeats, which ranged from 53.13% in C. aureobracteatum to 42.42% in C. macrophyllum (Additional File 4: Supplementary Table S8 and Additional File 5: Supplementary Figure S3B).
Boundary regions and comparative analysis
When comparing the chloroplast genomes of Chrysosplenium species, we found that IR/LSC junctions of IRb are largely located between rpl2 and rps19 (Fig. 3). Moreover, the overlap of ycf1 pseudogenes and ndhF appeared in different locations among the Chrysosplenium species: in the region of the SSC for C. ramosum, and at the IRb/SSC border for the other five species. C. alternifolium did not contain the ycf1 pseudogene. The ycf1 genes were sited at the SSC/IRa boundary and the length of ycf1 ranged from 5402 to 5546 bps. The trnH genes of the seven Chrysosplenium species were located in the LSC region, 2–19 bp away from the IRa–LSC border.
When comparing the genome boundaries of the Chrysosplenium species to the other three non-Chrysosplenium species of Saxifragaceae, ndhF was at the IRb/SSC boundary in most species of Chrysosplenium, except for C. ramosum, which showed contraction of the SSC and expansion of IRb. In addition, S. stolonifera was slightly different from the other two non-Chrysosplenium species in Saxifragaceae. In S. stolonifera, the contraction of the LSC region resulted in the rpl22 gene being at the IRb/LSC junction, which placed the whole rps19 gene in the IRb region. The rps19 pseudogenes were also found in the IRa region in S. stolonifera and B. scopulosa. When these data were combined with the phylogenetic tree of the three clades (S. stolonifera, B. scopulosa and O. rupifraga, and Chrysosplenium) inferred from whole-chloroplast protein-coding genes (Fig. 3), we found that the chloroplast genome structure within Chrysosplenium species is not strongly conserved, although the gene content is conserved.
LAGAN and Shuffle-LAGAN gave very similar results in the genetic divergence among the chloroplast genomes of Saxifragaceae species (Fig. 4 and Additional File 6: Supplementary Figure S4). The chloroplast genomes of the Chrysosplenium species were more conserved when compared with the three non-Chrysosplenium species of Saxifragaceae, and the intergenic spacer (IGS) regions had the highest levels of divergence: trnK–rps16, rps16–trnQ, rpoB–trnC, petN–psbM, trnT–psbD, psbZ–trnG, trnT–trnL, accD–psaI, ycf4–cemA, ndhF–rpl32, and rps15–ycf1. In addition, we found some highly variable coding sequences (ndhD, ycf2, ndhA, and ycf1), and the IR regions were more conserved than LSC and SSC regions in all the species tested. We also found slight difference for rpoC2, ycf2, and ycf1, which correspond to the difference between the Alternifolia and Oppositifolia subgenera.
Selective pressure analyses
We calculated the Ka/Ks ratios, the ratios of the rate of non-synonymous substitutions (Ka) to the rate of synonymous substitutions (Ks), at the species level by concatenating all of the 79 genes into a super-matrix. In Chrysosplenium species, the Ka/Ks ratios were around 0.2. This result suggested that at the whole-chloroplast protein level, the Chrysosplenium species have been subjected to a stronger purifying selection (Fig. 5, Additional File 4: Supplementary Table S9 and Additional File 4: Supplementary Table S10).
The Ka/Ks ratios were also calculated for all of the 79 protein-coding genes of the ten chloroplast genomes of Chrysosplenium separately (Fig. 6 and Additional File 4: Supplementary Table S7). Two genes (matK, ycf2) had Ka/Ks ratios around 1.0 in most species, implying possible positive selection. Specially, matK showed an average Ka/Ks ratio of 0.74 when compared with C. ramosum. Among the Chrysosplenium species, ycf2 often had a ratio higher than 0.8. Most of the other genes had a Ka/Ks ratio range from 0.1–0.3, implying strong purification (Additional File 4: Supplementary Table S10).
Sixty-six single-copy genes were used for selective pressure estimation with the branch-site model (Additional File 4: Supplementary Table S9 and Supplementary Table S11). We found that matK was positively selected in Chrysosplenium with the p-value = 0.022 and the Bayes Empirical Bayes (BEB) posterior probability for one amino acid site (117S, from polar Ser to non-polar Val) larger than 0.972. And the gene of ycf2 was also positively selected in Chrysosplenium with p-value = 0.00003 and the BEB posterior probability for 0.953 in 1028 K (from Lys to Leu). In addition, positively selected sites were detected for 18 genes (atpB, atpE, atpF, atpI, cemA, clpP, matK, ndhC, ndhE, ndhF, ndhH, ndhK, petA, psaB, psbH, psbJ, psbN, rps14, rps16) (Fig. 7 and Additional File 4: Supplementary Table S9).
Phylogenetic analyses yielded a well-supported phylogeny of Saxifrageles with most of the nodes having maximum likelihood (ML) bootstrap support values > 95 and bayesian inference (BI) posterior probabilities =1 (Fig. 8). The topologies yielded from ML analysis and BI analysis were completely identical. The topology of Saxifragales in our study was similar to the APG IV system  with Saxifragaceae closer to Iteaceae, phylogenetically. And, Chrysosplenium was divided into two clades corresponding to the two subgenera (Alternifolia and Oppositifolia) in our phylogenetic tree.
Gene numbers were slightly different due to the loss of rpl32 or transfer of the gene to the nucleus [22,23,24,25,26]. The gene transfer and its dormant expression in Alternifolia of Chrysosplenium could be explained by the decreased demand on photosynthesis and plastid translational capacity, which increased the success rate of gene transfer to the nucleus . Based on the previous study and our results, we speculate the rpl32 gene might transfer to the nucleus due to the adaptation for the specific habitat of Chrysosplenium, although this hypothesis remains to be verified by experiments.
Lower GC contents in Chrysosplenium chloroplast genomes compared to other Saxifragaceae members can be explained by the natural selection . DNA sequences of closely related species from different environments show marked differences in GC content, which has a direct impact on the amino acid sequences of the proteins in the respective environments . Genes with low GC contents are more prone to be transcribed than those with high GC content as GC pairs have three hydrogen bonds, making them more stable than AT pairs with two hydrogen bonds . Therefore, the selective pressure of the unique habitat of Chrysosplenium species (insufficient light energy) resulted in the lower overall GC contents and GC3 contents in their chloroplast genomes.
Compared with other genera within Saxifragaceae, the genus has the lowest light requirement. Among the 79 chloroplast genes in plants, 46 are related to photosynthesis pathway (Table 2) . Genes related to a specific environment are normally assumed to be under positive selection . This assumption was widely used to detect genes related to environmental adaptation . Our expectation was that the 46 genes were under positive selection. However, the lower Ka/Ks ratios at the chloroplast genome level within the Chrysosplenium species compared to non-Chrysosplenium species indicated that most genes were subjected to purifying selection to retain conserved functions in the Chrysosplenium. In the opposite environment, sunlight, including UV radiation, induces DNA damage, mutations and rearrangements [33, 34], which may contribute to an increase in mutation rates. Moreover, it has been proposed that more solar radiation and higher temperatures increase metabolism and growth rates, shortening generation times and increasing mutation rates . Meanwhile, low light is also a stress to the plants which can increase the mutation rates. As one of the most prevalent form in natural selection, purifying selection constantly sweeps away deleterious mutations in population. Therefore, the purifying selection on most chloroplast genes within Chrysosplenium would be evolutionary result of the preservation of the adaptive characteristics of Chrysosplenium species.
Plants have a variety of strategies to adapt to the environment; therefore, multiple genes may have been subjected to positive selection in Chrysosplenium during its adaptation. The gene matK is transcribed from the sole intact plastid group IIA intron ORF localized between the exons coding for the lysine-tRNA (trnK-UUU). In contrast to other group IIA ORFs, matK has lost domains assigned to a reverse transcriptase and endonuclease function [31, 36]. matK is usually used as a phylogenetic signal that can resolve evolutionary relationships because of its the high nucleotide and amino acid substitution rates [37, 38]. However, an extremely significant positive selection site at the 117S loci in matK of Chrysosplenium suggests this positive selection fixes beneficial variations within Chrysosplenium. The gene ycf2 is the largest chloroplast gene reported in angiosperms . ycf2 has become a useful gene for assessing sequence variation and evolution in plants . Positive selection of ycf2 was also found to be involved in adaptation in other species . However, due to the extremely high Ka/Ks value (ѡ = 231.51421) and unknown function, ycf2 is a valuable resource for future research of the adaptive evolution of Chrysosplenium. Further functional studies on the adaptive amino acid sites of chloroplast genes in Chrysosplenium are needed.
In this study, we sequenced the chloroplast genomes of six Chrysosplenium species and revealed the chloroplast genomic features between the Oppositifolia (C. macrophyllum, C. flagelliferum, C. alternifolium, and C. ramosum) and Alternifolia (C. ramosum, C. kamtschaticum, and C. sinicum) subgenera. In addition, we combined these six sequences with the previously reported chloroplast genomes of C. aureobracteatum (Oppositifolia), S. stolonifera (Saxifragaceae), B. scopulosa (Saxifragaceae), and O. rupifraga (Saxifragaceae). We discussed the comprehensive features of the chloroplast genomes, such as gene content and GC content, in these seven species of Chrysosplenium. All the species of Saxifragaceae shared similar genome structures, whereas the seven species of Chrysosplenium showed a lower average GC content, indicating selective pressure in their unique habitats. At the chloroplast genome level, the Ka/Ks ratios of the individual sequences showed that Chrysosplenium species were subjected to purifying selection compared to the non-Chrysosplenium species. At the level of amino acid sites, we found that matK and ycf2 were under positive selection with both high posterior probability and statistical significance. Other 15 genes with lower posterior probability involved in photosynthesis out of 19 were also possibly subjected to positive selection via a change of amino acid sites, which may be the adaptive response to its moist and shaded habitat. Using the protein-coding sequences from the whole chloroplast genome of 31 species, the robust consensus of phylogenetic trees reconstructed with both ML and BI algorithms suggested that Chrysosplenium species are sister to B. scopulosa and O. rupifraga within Saxifragaceae of Saxifragales. Also, our results supported the classification of the genus into two subgenera based on the morphology of opposite leaves or alternate leaves. These findings will be valuable for further study of the chloroplast genomes of Chrysosplenium species and provide valuable resources for studies of plant adaptation to low light conditions.
Sampling and sequencing
To represent the Chrysosplenium, six species were selected based on their morphological characteristics: C. macrophyllum, C. flagelliferum, and C. alternifolium belonging to the Alternifolia subgenus, and C. ramosum, C. kamtschaticum, and C. sinicum belonging to the Oppositifolia subgenus. Among the six species from wild, three were collected from China and three were from Japan (See details in Additional File 2: Supplementary Table S1). Due to a high content of secondary metabolites, the chloroplast DNA of C. macrophyllum was extracted using a high-salt method . To get a complete chloroplast genome, which can be used as a reference in assembling chloroplast genome of the other five species, C. macrophyllum was sequenced using the PacBio Sequel I platform at Frasergen (Wuhan, China) and the Illumina Hiseq 2500 at the Novogene Company (Beijing, China). The total genomic DNA for the other five Chrysosplenium species was extracted using a modified cetyltrimethylammonium bromide (CTAB) method  and sequenced using the Illumina Hiseq 2500 platform at the Novogene Company (Beijing, China).
Chloroplast genome assembly and annotation
The sequencing of the chloroplast DNA of C. macrophyllum with the PacBio Sequel I platform generated 218,330 reads with the N50 of 4452 bp. De novo genome assembly was conducted using Canu (v1.5) , which produced 4028 contigs with an N50 of 5011 bp. To discard nuclear DNA sequences, we aligned the contigs to a whole-chloroplast reference genome with the Burrow-Wheeler Aligner bwa . Then the contigs were polished with Arrow implemented in SMRT Link v6.0.0. Finally, the draft chloroplast genome was manually adjusted based on the two inverted repeats and scaffolds assembled from the Illumina Hiseq 2500 platform.
Sequence data for the other five species generated from the Illumina Hiseq 2500 platform were processed to remove the low-quality reads and adaptors. The clean reads were aligned to the complete chloroplast genome of C. macrophyllum with bwa-0.7.12 . The aligned reads were then selected for de novo assembly with ABYSS-2.0.2  after the optimal K-mer was chosen with the software kmergenie . Then, the contigs were connected with Sequencher 5.4.6 and scaffolded again with the original data by the software SSPACE_Standard_v3.0 . Last, the assembled scaffolds were manually adjusted based on the two inverted repeats and verified by Sanger sequencing (Additional File 2: Supplementary Table S2). We also assembled these scaffolds with GetOrganelle  to validate the ones assembled with ABYSS-2.0.2.
Gene annotation was performed using CPGAVAS2  and PGA . The different annotations of protein-coding sequences were confirmed using BLASTx. The tRNAs were checked with tRNAscan-SE v2.0.3 . Final chloroplast genome maps were drawn using OGDRAW .
To identify the nuclear homologs of chloroplast rpl32 (Cp_rpl32) in C. sinicum, the protein sequence of Cp_rpl32 was subjected to a BLASTp analysis against our the whole protein sequences from our nuclear genome of C. sinicum with the threshold of 30% identity, and the nuclear rpl32 genes were further annotated via the BLASTp of these potential nuclear homologs against NCBI non-redundant protein database with the threshold of 30% identity and an e-value <1E− 10.
Analysis of GC content, nucleotide diversity, and repeat content
We accessed chloroplast genome sequences of Chrysosplenium aureobracteatum (MG878089; Saxifragaceae; Chrysosplenium), Saxifraga stolonifera (MH191389; Saxifragaceae; Saxifraga), Bergenia scopulosa (KY412195; Saxifragaceae; Bergenia), and Oresitrophe rupifraga (MF774190; Saxifragaceae; Oresitrophe) from GenBank to compare the features among these chloroplast genomes in Saxifragaceae. The nucleotide diversity (Pi) among the seven species of Chrysosplenium was calculated using the software DnaSP v6.12.03 . The GC content of the whole chloroplast genome and the third position GC content of codons for all ten species were calculated using an in-house Python script. The simple sequence repeats (SSRs) were detected using the MIcroSAtellite (MISA) identification tool with the minimum repeat number set at 10, 5, 4, 3, 3, and 3 for mono-, di-, tri-, tetra-, penta-, and hexanucleotides, respectively. We also identified tandem repeat sequences using REPuter  with minimal repeats of more than 30 bp and hamming distances of less than 3 bp.
Boundary regions and comparative analysis
The contraction or expansion between boundary regions of the chloroplast genome in each species was drawn by IRscope . To compare the conservation of each gene, we visualized the results with mVISTA through two alignment programs: LAGAN, which produces true multiple alignments regardless of whether they contain inversions or not, and Shuffle-LAGAN, which can detect rearrangements and inversions in sequences .
Selective pressure estimation
We carried out selective pressure estimation for the 6 species of Chrysosplenium and 25 species of non-Chrysosplenium with two strategies: calculation based on pairwise comparison and calculation based on the branch-site model:
All protein-coding sequences (CDSs) from each of the 31 species were concatenated into a super matrix for inferring phylogenetic tree. Then, species vs. species Ka/Ks ratio was estimated. In addition, the Ka/Ks ratio was estimated for each of the 79 genes within Chrysosplenium separately. The CDS for each gene was translated to amino acid sequences, which were aligned with MEGA7 . Then, the corresponding CDS were aligned according to the amino acid sequences. Lastly, Ka/Ks ratios were calculated using the KaKs-calculator v 2.0 . Genes with 1 < Ka/Ks ratio < 45 were considered as under positive selection; genes with Ka/Ks ratio < 1 were considered as under purifying selection. The ratio > =45 or NA indicates that the gene has few nonsynonymous sites/substitutions, and was not considered in our analysis.
A total of 66 CDSs presented in all the analysed species, and were used for identification of positive selection using the branch-site model . CDSs of each gene were aligned according to their amino acid sequences with MEGA7 . The branch-site model in the program codeml of the PAML v4.9 package  was used to assess potential positive selection in Chrysosplenium that was set as the foreground branch. Selective pressure is measured by the ratio (ω) of the nonsynonymous substitution rate (dN) to the synonymous substitutions rate (dS). A neutral branch-site model (Model = 2, NSsites = 2, Fix = 1, and Fix ω = 1) and an alternative branch-site model (Model = 2, NSsites = 2, and Fix = 0) were applied separately. The right-tailed Chi-square test was used to compute p-values based on the difference of log-likelihood values between the two models with one degree of freedom. Moreover, BEB method  was implemented to calculate the posterior probabilities for amino acid sites that are potentially under positive selection. A gene with a p-value < 0.05 and ω > 1 was considered as a positively selected gene. An amino acid site with posterior probabilities > 0.95 was considered as positively selected.
To construct a phylogeny of Chrysosplenium, 29 species of Saxifrageles and two Buxaceae species (as outgroups) were selected (see Additional file 2: Supplementary Table S3 for details). The whole-chloroplast protein-coding genes of these 31 species were aligned with MUSCLE v3.8.31 . The best-fitting nucleotide substitution model was determined using the Akaike Information Criterion in the model-finder IQ-TREE . ML analysis was performed using IQ-TREE with the best model of GTR + F + R4 and 1000 bootstrap replicates, and BI analysis was performed in MrBayes 3.2.6  using the Markov Chain Monte Carlo method with 200,000 generations and sampling trees every 100 generations. The first 20% of trees were discarded as burn-in with the remaining trees being used for generating a consensus tree.
Availability of data and materials
All chloroplast genomes used in this study can be found in Genbank and their Genbank accessions can be found in Additional File 2: Supplementary Table S1 and Additional File 2: Supplementary Table S3. The whole protein sequences of the nuclear genome of C. sinicum in our study are archived at the Dryad Digital Repository (https://doi.org/10.5061/dryad.jdfn2z38m). The other data sets generated in this study are included within the article and additional files. All materials used or generated during the study are kept in our laboratory. Most figures were completed by the open source software, R package; some figures were done by the relevant software with the proper citations.
Large single copy region
Inverted repeat region
Small single copy region
Polymerase chain reaction
The rate of non-synonymous substitutions to the rate of synonymous substitutions
Simple sequence repeats
Bayes empirical bayes
Intergenic spacer region
Fragments per kilobase of exon per million reads
Wu Z, Gui S, Quan Z, Pan L, Wang S, Ke W, Liang D, Ding Y. A precise chloroplast genome of Nelumbo nucifera (Nelumbonaceae) evaluated with sanger, Illumina MiSeq, and PacBio RS II sequencing platforms: insight into the plastid evolution of basal eudicots. BMC Plant Biol. 2014;14:289.
Alzahrani DA, Yaradua SS, Albokhari EJ, Abba A. Complete chloroplast genome sequence of Barleria prionitis, comparative chloroplast genomics and phylogenetic relationships among Acanthoideae. BMC Genomics. 2020;21(1):393.
Cui L, Leebens-Mack J, Wang L, Tang J, Rymarquis L, Stern DB, Claude WD. Adaptive evolution of chloroplast genome structure inferred using a parametric bootstrap approach. BMC Evol Biol. 2006;6:13.
Scott-Phillips TC, Laland KN, Shuker DM, Dickins TE, West SA. The niche construction perspective: a critical appraisal. Evolution. 2014;68(5):1231–43.
R. Marcelino V, MCM C, Jackson CJ, AAW L, Verbruggen H. Evolutionary dynamics of chloroplast genomes in low light: a case study of the endolithic green alga Ostreobium quekettii. Genome Biol Evol. 2016;8(9):2939–51.
The Angiosperm Phylogeny G. An update of the angiosperm phylogeny group classification for the orders and families of flowering plants: APG IV. Bot J Linn Soc. 2016;181(1):1–20.
POWO. Plants of the world online. In: facilitated by the Royal Botanic Gardens. Kew: Published on the Internet; 2019. http://powo.science.kew.org/taxon/urn:lsid:ipni.org:names:30129301-30129302#children.
Hara H. Synopsis of the genus Chrysosplenium L. J Fac Sci Univ Tokyo. 1957;3(7):1–90.
Soltis DE, Tago-Nakazawa M, Xiang Q-Y, Kawano S, Murata J, Wakabayashi M, Hibsch-Jetter C. Phylogenetic relationships and evolution in Chrysosplenium (Saxifragaceae) based on matK sequence data. Am J Bot. 2001;88(5):883–93.
Pan J, Ohba H. Chrysosplenium. In: Wu ZY, Raven PH, editors. Flora of China, vol. 8. Beijing and St. Louis: Science press and Missouri botanical garden press; 2001. p. 346–58.
Liu H, Luo J, Liu Q, Lan D, Qin R, Yu X. A new species of Chrysosplenium (Saxifragaceae) from Zhangjiajie, Hunan, Central China. Phytotaxa. 2016;277:287–92.
Franchet A. Monographie du genere Chrysosplenium tourn. Nouv Arch duMuséum Hist Nat. 1891;3(2):1–32.
Guo Q, Ricklefs R, Cody M. Vascular plant diversity in eastern Asia and North America: historical and ecological explanations. Bot J Linn Soc. 1998;128(2):123–36.
Xiang QY, Zhang WH, Ricklefs RE, Qian H, Chen ZD, Wen J, Li JH. Regional differences in rates of plant speciation and molecular evolution: a comparison between eastern Asia and eastern North America. Evolution. 2004;58(10):2175–84.
Ibrahim K, De Luca V, He K, Latchinian L, Brisson L, Charest PM. Enzymology and compartmentation of polymethylated flavonol glucosides in Chrysosplenium americanum. Phytochemistry. 1987;26:1237–45.
Han J, Yang S, Kim H, Jang C, Park J, Kang S. Phylogenetic study of Korean Chrysosplenium based on nrDNA ITS sequences. Korean J Plant Resour. 2011;24(4):358–69.
Kim Y, Kim Y. Molecular systematic study of Chrysosplenium series pilosa (Saxifragaceae) in Korea. J Plant Biol. 2011;54(6):396.
Ravi V, Khurana JP, Tyagi AK, Khurana P. An update on chloroplast genomes. Plant Syst Evol. 2008;271(1):101–22.
Parks M, Cronn R, Liston A. Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes. BMC Biol. 2009;7:84.
Xie D, Yu H, Price M, Xie C, Deng Y, Chen J, Yu Y, Zhou S, He X. Phylogeny of chinese Allium species in section daghestanica and adaptive evolution of Allium (Amaryllidaceae, Allioideae) species revealed by the chloroplast complete genome. Front Plant Sci. 2019;10:460.
Kim Y, Lee J, Kim Y. The complete chloroplast genome of a Korean endemic plant Chrysosplenium aureobracteatum Y.I. Kim & Y.D. Kim (Saxifragaceae). Mitochondrial DNA Part B. 2018;3:380–1.
Huang CY, Ayliffe MA, Timmis JN. Direct measurement of the transfer rate of chloroplast DNA into the nucleus. Nature. 2003;422(6927):72–6.
Stegemann S, Hartmann S, Ruf S, Bock R. High-frequency gene transfer from the chloroplast genome to the nucleus. Proc Natl Acad Sci. 2003;100(15):8828.
Park S, Jansen RK, Park S. Complete plastome sequence of Thalictrum coreanum (Ranunculaceae) and transfer of the rpl32 gene to the nucleus in the ancestor of the subfamily Thalictroideae. BMC Plant Biol. 2015;15(1):40.
Ueda M, Fujimoto M. Arimura S-i, Murata J, Tsutsumi N, Kadowaki K-i. loss of the rpl32 gene from the chloroplast genome and subsequent acquisition of a preexisting transit peptide within the nuclear gene in Populus. Gene. 2007;402(1):51–6.
Cusack BP, Wolfe KH. When gene marriages don't work out: divorce by subfunctionalization. Trends Genet. 2007;23(6):1208–13.
Fleischmann TT, Scharff LB, Alkatib S, Hasdorf S, Schöttler MA, Bock R. Nonessential plastid-encoded ribosomal proteins in tobacco: a developmental role for plastid translation and implications for reductive genome evolution. Plant Cell. 2011;23(9):3137–55.
Charlesworth B. Genetic recombination: patterns in the genome. Curr Biol. 1994;4(2):182–4.
Foerstner KU, von Mering C, Hooper SD, Bork P. Environments shape the nucleotide composition of genomes. EMBO Rep. 2005;6(12):1208–13.
Jia Q, Wu H, Zhou X, Gao J, Zhao W, Aziz J, Wei J, Hou L, Wu S, Zhang Y, et al. A "GC-rich" method for mammalian gene expression: a dominant role of non-coding DNA GC content in regulation of mammalian gene expression. Sci China Life Sci. 2010;53(1):94–100.
Wicke S, Schneeweiss GM. dePamphilis CW, Müller KF, Quandt D. the evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011;76(3):273–97.
Yang Z, Wong WSW, Nielsen R. Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005;22(4):1107–18.
Raven JA, Beardall J, AWD L, Sánchez-Baracaldo P. Interactions of photosynthesis with genome size and function. Philos Trans R Soc Lond B Biol Sci. 2013;368(1622):20120264.
Kumar RA, Oldenburg DJ, Bendich AJ. Changes in DNA damage, molecular integrity, and copy number for plastid DNA and mitochondrial DNA during maize development. J Exp Bot. 2014;65(22):6425–39.
Klaus R. Latitudinal gradients in species diversity: the search for the primary cause. Oikos. 1992;65:514–27.
Liere K, Link G. RNA-binding activity of the matK protein encodecd by the chloroplast trnk intron from mustard ( Sinapis alba L.). Nucleic Acids Res. 1995;23(6):917–21.
Hilu KW, Borsch T, Müller K, Soltis DE, Soltis PS, Savolainen V, Chase MW, Powell MP, Alice LA, Evans R, et al. Angiosperm phylogeny based on matK sequence information. Am J Bot. 2003;90(12):1758–76.
Hilu KW, Liang G. The matK gene: sequence variation and application in plant systematics. Am J Bot. 1997;84(6):830–9.
Drescher A, Ruf S, Calsa T Jr, Carrer H, Bock R. The two largest chloroplast genome-encoded open reading frames of higher plants are essential genes. Plant J. 2000;22(2):97–104.
Huang J, Sun G, Zhang D. Molecular evolution and phylogeny of the angiosperm ycf2 gene. J Syst Evol. 2010;48(4):240–8.
Zhong Q, Yang S, Sun X, Wang L, Li Y. The complete chloroplast genome of the Jerusalem artichoke (Helianthus tuberosus L.) and an adaptive evolutionary analysis of the ycf2 gene. PeerJ. 2019;7:e7596.
Shi C, Hu N, Huang H, Gao J, Zhao Y, Gao L. An improved chloroplast DNA extraction procedure for whole plastid genome sequencing. PLoS One. 2012;7(2):e31468.
Doyle J. DNA protocols for plants-CTAB total DNA isolation. Mol Tech Taxonomy. 1991;57:283–93.
Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017;27(5):722–36.
Li H, Durbin R. Fast and accurate long-read alignment with burrows–wheeler transform. Bioinformatics. 2010;26(5):589–95.
Jackman SD, Vandervalk BP, Mohamadi H, Chu J, Yeo S, Hammond SA, Jahesh G, Khan H, Coombe L, Warren RL, et al. ABySS 2.0: resource-efficient assembly of large genomes using a bloom filter. Genome Res. 2017;27(5):768–77.
Chikhi R, Medvedev P. Informed and automated k-mer size selection for genome assembly. Bioinformatics. 2013;30(1):31–7.
Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2010;27(4):578–9.
Jin J, Yu W, Yang J, Song Y, Yi T, Li D. GetOrganelle: a simple and fast pipeline for de novo assembly of a complete circular chloroplast genome using genome skimming data. bioRxiv. 2018. https://doi.org/10.1101/256479.
Shi L, Chen H, Jiang M, Wang L, Wu X, Huang L, Liu C. CPGAVAS2, an integrated plastome sequence annotator and analyzer. Nucleic Acids Res. 2019;47(W1):W65–73.
Qu X, Moore MJ, Li D, Yi T. PGA: a software package for rapid, accurate, and flexible batch annotation of plastomes. Plant Methods. 2019;15(1):50.
Schattner P, Brooks AN, Lowe TM. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005;33(suppl_2):W686–9.
Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 2019;47(W1):W59–64.
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, Sánchez-Gracia A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34(12):3299–302.
Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R. REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001;29(22):4633–42.
Amiryousefi A, Hyvönen J, Poczai P. IRscope: an online program to visualize the junction sites of chloroplast genomes. Bioinformatics. 2018;34(17):3030–1.
Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Program NCS, Green ED, Sidow A, Batzoglou S. Lagan and multi-Lagan: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003;13(4):721–31.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Wang D, Zhang Y, Zhang Z, Zhu J, Yu J. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010;8(1):77–80.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
Yang Z, dos Reis M. Statistical properties of the branch-site test of positive selection. Mol Biol Evol. 2010;28(3):1217–28.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Nguyen L, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2014;32(1):268–74.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.
This work was supported by the Construction Plan of Hubei Province Science and Technology basic conditions platform (No. 2017BEC014), Fund for Key Laboratory Construction of Hubei Province (No. 2018BFC360) and Hubei Provincial Natural Science Foundation of China (No. 2019CFB214). All funders we mentioned provide financial support for our study.
Ethics approval and consent to participate
Not applicable. The plant was collected in non protected area; no any legal authorization/license is required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The protein sequences (A) and expression values in leaf (B) of both chloroplast rpl32 (Cp_rpl32) and its nuclear homolog (Nu_rpl32) in C. sinicum.
Supplementary Tables S1-S6.
Comparison of nucleotide diversity (Pi) between the chloroplast genomes of Chrysosplenium.
Supplementary Tables S7-S11.
Analyses of repeat sequences in the ten species of Saxifragaceae. (A) The analysis of simple sequence repeats (SSRs) in chloroplast genomes of Saxifragaceae. (B) The repeat types in Saxifragaceae.
The comparative analysis with Shuffle-LAGAN program of the whole chloroplast genome of seven different species from the family of Saxifragaceae.
About this article
Cite this article
Wu, Z., Liao, R., Yang, T. et al. Analysis of six chloroplast genomes provides insight into the evolution of Chrysosplenium (Saxifragaceae). BMC Genomics 21, 621 (2020). https://doi.org/10.1186/s12864-020-07045-4
- Chloroplast genome
- Opposite leaves
- Alternate leaves