- Research article
- Open Access
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 Biology volume 15, Article number: 40 (2015)
Plastids originated from cyanobacteria and the majority of the ancestral genes were lost or functionally transferred to the nucleus after endosymbiosis. Comparative genomic investigations have shown that gene transfer from plastids to the nucleus is an ongoing evolutionary process but molecular evidence for recent functional gene transfers among seed plants have only been documented for the four genes accD, infA, rpl22, and rpl32.
The complete plastid genome of Thalictrum coreanum, the first from the subfamily Thalictroideae (Ranunculaceae), was sequenced and revealed the losses of two genes, infA and rpl32. The functional transfer of these two genes to the nucleus in Thalictrum was verified by examination of nuclear transcriptomes. A survey of the phylogenetic distribution of the rpl32 loss was performed using 17 species of Thalictrum and representatives of related genera in the subfamily Thalictroideae. The plastid-encoded rpl32 gene is likely nonfunctional in members of the subfamily Thalictroideae (Aquilegia, Enemion, Isopyrum, Leptopyrum, Paraquilegia, and Semiaquilegia) including 17 Thalictrum species due to the presence of indels that disrupt the reading frame. A nuclear-encoded rpl32 with high sequence identity was identified in both Thalictrum and Aquilegia. The phylogenetic distribution of this gene loss/transfer and the high level of sequence similarity in transit peptides suggest a single transfer of the plastid-encoded rpl32 to the nucleus in the ancestor of the subfamily Thalictroideae approximately 20–32 Mya.
The genome sequence of Thalictrum coreanum provides valuable information for improving the understanding of the evolution of plastid genomes within Ranunculaceae and across angiosperms. Thalictrum is unusual among the three sequenced Ranunculaceae plastid genomes in the loss of two genes infA and rpl32, which have been functionally transferred to the nucleus. In the case of rpl32 this represents the third documented independent transfer from the plastid to the nucleus with the other two transfers occurring in the unrelated angiosperm families Rhizophoraceae and Salicaceae. Furthermore, the transfer of rpl32 provides additional molecular evidence for the monophyly of the subfamily Thalictroideae.
Massive transfer of genes from the plastid to the nucleus occurred following the endosymbiotic origin of the plastid from cyanobacteria . Photosynthetic land plant plastid genomes (plastomes) only encode 101–118 genes, most of which represent genetic system and photosynthetic genes [2,3]. A considerable number of organelle-targeted genes in the nucleus are translated in the cytosol and imported into the plastids and mitochondria where they perform essential functions. Many studies have revealed that gene transfer from organelles to the nucleus is an ongoing process [1,4], however subsequent molecular characterization of these events has been limited. Transferred plastid genes must obtain nuclear expression elements as well as transit peptides for import of gene products into the plastids [5,6]. Successful functional gene transfers from the plastid to the nucleus in seed plants have been documented for only four genes: infA in multiple lineages , rpl22 in Fabaceae and Fagaceae [8,9], rpl32 in Rhizophoraceae and Salicaceae [10,11] and accD in Trifolium [12,13]. Transferred plastid genes have either adopted a transit peptide from an existing nuclear gene or acquired a novel transit peptide [9,10,13]. In addition to functional gene transfers, movement of DNA fragments from the plastid to the nucleus is common among flowering plants (referred to as NUPTs; nuclear plastid DNA) [1,14], and the proportion of NUPTs differs considerably among species [15,16].
The angiosperm family Ranunculaceae (buttercups) exhibits enormous ecological, anatomical, biochemical, and morphological diversity and comprises approximately 2,500 species in 59 genera and five subfamilies distributed throughout the world . Ranunculaceae have two chromosome types: R (Ranunculus)-type with large chromosomes, and T (Thalictrum)-type with small chromosomes [17,18]. Although there are several different classification systems for Ranunculaceae [17,19-23], multiple lines of evidence suggest that genera with the T-type chromosome (excluding Hydrastis) form a monophyletic group [22-24]. Thalictrum, a member of the subfamily Thalictroideae, is one of the most diverse genera of Ranunculaceae in terms of number of species and morphological variation . Recent studies have estimated phylogenetic relationships of Thalictrum using molecular data to understand the evolution of sexual systems and polyploidy [25,26]. This genus has great medicinal value because it contains high levels of Thaliblastin (Thalicarpine), which has anticancer properties [27,28]. Thalictrum coreanum is a popular, economically important endemic plant native to Korea and it is used widely in horticulture and medicine. Its natural habitat is restricted to small areas in Korea and it is often confused with a species of Berberidaceae, Epimedium koreanum, which is used in traditional Chinese and Korean herbal medicine as a potent enhancer of erectile function.
Previous studies performed restriction site mapping of the plastid genome of Ranunculaceae and identified several phylogenetically informative rearrangements, including inversions, the loss of the rps16 gene and loss of the rps12 cis-spliced intron [29,30]. The complete plastid genome sequences of only two species of Ranunculaceae have been reported [31,32] and neither of these are members of the subfamily Thalictroideae.
In this study the plastome sequence of T. coreanum is presented, which represents the first sequenced member of the subfamily Thalictroideae. Genome organization is examined, including identification of transfers of two genes, infA and rpl32, from the plastid to the nucleus. In addition, the phylogenetic distribution of the rpl32 gene loss in the Ranunculaceae is examined. The plastome sequence of T. coreanum provides valuable additional information about variation within the Ranunculaceae.
Plastome of Thalictrum coreanum
The Thalictrum coreanum plastome is 155,088 bp with a pair of inverted repeats (IRs) of 26,403 bp separated by a small single copy (SSC) region of 17,549 bp and a large single copy (LSC) region of 84,733 bp (Figure 1A and Table 1). The genome encodes 112 different genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes and consists of 58.23% genes (i.e. protein-coding, tRNA, and rRNA genes) (Table 1). The translation initiation factor A (infA) is a pseudogene due to the presence of frameshift mutations. The ribosomal protein L32 (rpl32), which is usually located between ndhF and trnL-UAG (Figure 1A), is a pseudogene because deletions near the 5’ end generate two internal stop codons.
General features of the plastomes of three Ranunculaceae are summarized in Table 1. Compared with two other sequenced Ranunculaceae plastomes [31,32], Megaleranthis saniculifolia and Ranunculus macranthus, changes in genome organization reflect shifts of the IRs at the LSC/IR boundary relative to Nicotiana tabacum (Figure 1B). For example, IRb of T. coreanum and M. saniculifolia extends into the LSC to include the N-terminal portion of rps19, generating a truncated rps19 fragment in IRa. However, in R. macranthus, IRa extends into the LSC to include the C-terminal portion of trnH-GUG, generating a trnH-GUG fragment in IRb. In terms of gene losses, the infA loss is shared by T. coreanum and R. macranthus, whereas M. saniculifolia contains an intact infA gene in its plastome. The presence of rpl32 as a pseudogene is unique to T. coreanum among all three Ranuculaceae analyzed.
Identification of functional gene transfers to the nucleus
To determine if the plastid-encoded rpl32 gene in Thalictrum has been transferred to the nucleus, the transcriptome database (1KP project) for T. thalictroides was queried with the rpl32 coding sequence of M. saniculifolia and R. macranthus. A transcript with high sequence identity to rpl32 is present and has an extended sequence of 417 bp upstream from the conserved ribosomal protein L32 domain (CHL00152). The first 66 amino acids of the open reading frame (ORF) is predicted by both TargetP and Predotar to be a transit peptide that is targeted to the plastid (Table 2). The extended region including the transit peptide had no significant hits with BlastN to any sequences in the NCBI databases and Phytozome genomics portal. Extensive searching of the Phytozome genomics portal revealed the presence of a nuclear-encoded rpl32 ORF in Aquilegia coerulea, which is also a member of the subfamily Thalictroideae. The sequence upstream from the conserved domain also has a transit peptide (66 amino acids; Table 2). However, an rpl32-like gene sequence was not detected in the Hydrastis canadensis transcriptome. Alignment of the nuclear-encoded rpl32 from Thalictrum and Aquilegia revealed a pairwise nucleotide sequence identity of 94.2% and 93.2% for the extended region and the conserved domain, respectively (Figure 2A). Amino acid alignment of four nuclear-encoded rpl32 copies (Aquilegia, Thalictrum, Populus [AB302219], and Bruguiera [AM711843]) shows that the extended region of Thalictrum is highly similar to Aquilegia with 89.9% identity, whereas Populus and Bruguiera are highly divergent with very low identities (19.3% and 16.5%) to Thalictrum (Figure 2B). The conserved ribosomal protein L32 domain of nuclear and plastid copies has pairwise identities ranging from 61.4% to 100% (Figure 2B).
Phylogenetic analyses of the nuclear-encoded rpl32 copies (Aquilegia, Bruguiera, Thalictrum, and Populus) and the plastid-encoded copies from 48 other angiosperms show that the Thalictrum and Aquilegia nuclear copies are nested within a clade with the plastid copies of the two Ranunculaceae Ranunculus and Megaleranthis, and the Populus and Bruguiera nuclear-encoded copies are grouped with the rosid Cucumis (Additional file 1: Figure S1). The nuclear copies of Thalictrum and Aquilegia group together with high bootstrap support (100%). The branch lengths on the tree indicate that the four nuclear-encoded copies have much higher substitution rates compared to plastid-encoded copies of closely related species. However, bootstrap support across the angiosperms is weak because the tree is based on only a single, short gene sequence.
To examine rate variation further, pairwise analysis of nonsynonymous (d N ) and synonymous (d S ) substitutions for plastid and nuclear rpl32 homologs was performed (Figure 3). The analysis shows higher divergence in both Aquilegia and Thalictrum nuclear-encoded genes compared to other species of Ranunculaceae. Higher sequence divergence in the Populus nuclear-encoded copy is also evident. The synonymous substitution rate of Thalictrum and Aquilegia clade is 2.5 and 8.8 times higher than their closest relatives Megaleranthis and Ranunculus, respectively. The branch lengths on the tree indicate that the Thalictrum copy has experienced much higher synonymous substitution rates than Aquilegia (Figure 3A). The correlation of d N and d S was moderate (P < 1 x 10−15, r = 0.7547). The d N /d S ratio among plastid copies shows similar patterns with d S larger than d N , which is also the case for the three nuclear copies (Figure 3B).
In addition, a Blast search of the T. thalictroides transcriptome from the 1KP database identified one or more transcripts of the translation initiation factor IF1 (cd04451) domain that has a transit peptide for targeting back to the plastid (Table 2). The Aquilegia transcriptome databases from Phytozome v.10 were queried with the infA domain sequence from the Thalictrum nuclear copy, confirming an infA-like ORF acquired a transit peptide (Table 2). Examination of the Aquilegia coerulea v1.1 nuclear genome (Phytozome; scaffold_1) showed the presence of the nuclear-encoded infA gene containing two exons totaling 1,171 bp separated by a 105 bp intron (Additional file 1: Figure S2). Nuclear-like infA sequences were not detected in the Hydrastis transcriptome.
Characterization of rpl32 gene in the subfamily Thalictroideae
The plastid-encoded rpl32 is a pseudogene in T. coreanum (Figure 1A). Seventeen additional species of Thalictrum representing two subgenera were surveyed for the presence of a pseudogene using PCR and Sanger sequencing (Figure 4). In T. thalictroides, PCR failed to amplify a product, which may be due to variation in primer binding sites. The product sequence sizes for the other 16 species of Thalictrum range from 745 bp in T. alpinum to 1,198 bp in T. rochebrunianum (median size of 16 Thalictrum species was 1,104 bp). Blast searches using intact rpl32 from M. saniculifolia (174 bp) and R. macranthus (162 bp) revealed that 15 examined Thalictrum species have remnant sequences of rpl32, ranging from 164 to 210 bp (Figure 5). However, one species, T. alpinum, lacks any detectable rpl32-like sequences, suggesting a loss of the entire gene. Nucleotide alignment of rpl32 revealed a consistent pattern, the majority of indel events are shared by members of the T. coreaum clade (Figures 4 and 5).
To further investigate the rpl32 gene loss, six other genera (Aquilegia, Enemion, Isopyrum, Leptopyrum, Paraquilegia, and Semiaquilegia) were examined in the subfamily Thalictroideae. The results show frameshift mutations due to insertions and deletions (indels) in five of the genera (Figure 5), and the sixth genus Leptopyrum has entirely lost rpl32. Maximum likelihood (ML) analysis of a concatenated data set resolves phylogenetic relationships among members of the subfamily Thalictroideae with bootstrap values of 98% for the monophyly of Thalictrum and 100% for the monophyly of subfamily Thalictroideae (Figure 4). Overall the rpl32 gene in the plastid genome of subfamily Thalictroideae is likely nonfunctional due to indels that disrupt the reading frame.
Correlation between reduction of ndhF-trnL intergenic spacer and rpl32 gene loss
The ndhF-trnL intergenic spacer (IGS) including rpl32 gene, which is either a pseudogene or absent within the subfamily Thalictroideae, shows considerable length variation (1.6-5.5 fold reduction compared to a full length IGS with rpl32, Figure 6A). This IGS region in the subfamily Thalictroideae is nearly two times shorter than in other angiosperms (Figure 6B). Both t-test and Wilcoxon signed rank test estimates indicated that the mean size of IGS between the two groups is significantly different (t-test; P < 1 × 10−13 and Wilcoxon signed-rank test; P < 1 × 10−8).
Functional gene transfer to the nucleus
Two protein coding genes, translation initiation factor A (infA) and ribosomal protein L32 (rpl32), are pseudogenes in the T. coreanum plastome. In case of infA, multiple independent losses have been reported across angiosperms including Caltha from the Ranunculaceae [7,33]. This previous report, combined with the phylogenetic distribution of infA loss from the sequenced Ranunculaceae genomes, indicates that this gene has been lost multiple times in the family. In order for a gene transfer event to be successful, transferred genes must acquire a transit peptide to shuttle the product back into plastids. Nuclear-encoded infA copies from Thalictrum and Aquilegia were identified in the transcriptome and they have high levels of sequence identity. In view of the high nucleotide sequence identity of both infA (94.1%) and the transit peptide (85.1%), it is likely that there has been a single transfer of this gene to the nucleus within the subfamily Thalictroideae, although expanded sampling is needed to confirm this hypothesis.
Most ribosomal protein subunits have been transferred to the nuclear genome since the endosymbiotic origin of plastids; however, land plant plastid genomes still retain a set of 12 small ribosomal protein subunits (rps) and 9 large ribosomal protein subunits (rpl) . Among the remaining plastid-encoded rps and rpl genes, several examples of gene losses across seed plants have been demonstrated [3,33]. Comparative analysis of the three sequenced Ranunculaceae plastomes (Megaleranthis, Ranunculus, and Thalictrum) indicates that the loss of the plastid-encoded rpl32 gene is unique to the Thalictrum plastome. However, comparisons of 17 additional Thalictrum species suggest that pseudogenization of the plastid encoded rpl32 gene occurred within the entire genus (Figure 4). Alignment of rpl32 pseudogenes from the sequenced Thalictrum species with intact rpl32 genes from M. saniculifolia and R. macranthus reveals that the majority of indel events are shared by members of the T. coreaum clade (Figure 5), indicating that the deletions occurred in the ancestor of this clade. Examination of the transcriptome sequences of Thalictrum and Aquilegia reveals that rpl32 has been transferred to the nucleus and acquired a target peptide for transport back to the plastid (Figure 2). The nuclear copies from Aquilegia and Thalictrum have high sequence identity at both nucleotide and amino acid levels (93.9% and 92.8%). The transferred genes have significantly elevated synonymous substitution rates and have experienced purifying selection (Figure 3). Phylogenetic analysis provided strong support for monophyly of the nuclear-encoded rpl32 copies (Figure 4), suggesting a single transfer of rpl32 to the nucleus. Plastid-encoded rpl32 gene losses have also been reported from Bruguiera, Populus, Yucca, and some parasitic plants [33,34]. There is evidence in only two of these cases, Bruguiera and Populus, that rpl32 has been functionally transferred to the nucleus [10,11]. In the case of Bruguiera and Populus rpl32 fused to an existing nuclear gene (Cu-Zn superoxide dismutase) to acquire a transit peptide, whereas Thalictrum and Aquilegia have acquired a novel transit peptide.
Loss of plastid-encoded rpl32 gene in the subfamily Thalictroideae
The high level of conservation of genome organization among the three sequenced Ranunculaceae plastomes enabled a PCR and sequencing survey of the ndhF and trnL-UAG region, which contains the rpl32 gene. The absence of intact rpl32 gene was identified for seven genera of the subfamily Thalictroideae (Aquilegia, Enemion, Isopyrum, Leptopyrum, Paraquilegia, Semiaquilegia, and Thalictrum) and the evolutionary fate of the plastid-encoded rpl32 differed among the genera or species examined; the gene is completely absent in Leptopyrum and T. alpinum and pseudogenes of varying length are present in the remaining species (Figure 6A). This suggests that rpl32 was transferred to the nucleus in the ancestor of subfamily Thalictroideae. Previous studies have shown that reductions of IGS regions are caused by gene loss, which has led to a more compact genome [35,36]. Although most examined Thalictroideae have a portion of rpl32 remaining, the ndhF-trnL intergenic spacer is significantly shorter in the subfamily Thalictroideae than in other angiosperms (Figure 6B) due to extreme degradation of the IGS. This finding indicates that the reduction of the ndhF and trnL-UAG IGS region is associated with the loss or pseudogenization of rpl32.
Two different types of chromosomes based on size have been characterized in Ranunculaceae, R-type and T-type [17,18]. The subfamilies Thalictroideae and Hydrastidoideae belong to T-type chromosome group, however, phylogenetic analyses have shown that these two subfamilies are polyphyletic [23,24]. The distribution of the transfer of rpl32 to the nucleus in Thalictrum and Aquilegia but not in Hydrastis indicates that this transfer does not represent a synapomorphy for the lineages with the T-type chromosomes.
Fior et al.  used the rbcL, matK and 26S nuclear ribosomal DNA (nrDNA) sequences generated by Wang et al.  to infer divergence times for the main clades of the Ranunculaceae. The divergence time of the subfamily Thalictroideae was estimated at 26.2 Mya (95% highest posterior density, HPD = 20.3-32.3 Mya). Another estimate indicated slightly later divergence times with the shorter interval for the subfamily Thalictroideae at 27.61 Mya (95% HPD = 26.6-28.6 Mya) . Thus, the transfer of rpl32 to the nucleus at the base of the subfamily Thalictroideae occurred approximately 20–32 Mya.
The monophyly of subfamily Thalictroideae has been confirmed based on phylogenetic analyses of multiple DNA markers: rbcL, matK, trnL-F spacer, and 26S nrDNA , 26S nrDNA , and atpB, rbcL, and 18S nrDNA . The rpl32 gene transfer event, combined with divergence time estimates, provides valuable phylogenetic data in support of the monophyly of subfamily Thalictroideae. Although there are multiple examples of plastid gene losses that exhibit homoplasy [e.g., 7, 9], the loss of rpl32 by all sampled members of subfamily Thalictroideae provides an excellent example of a genomic change that supports the monophyly of this subfamily.
The plastome sequence of Thalictrum coreanum, the first genome completed from the subfamily Thalictroideae, provides new insights into the evolution of plastomes within Ranunculaceae. The T. coreanum plastome is highly conserved with gene order identical to the ancestral organization of angiosperms  and at 155 kb it has the median genome size for photosynthetic land plants . The only unusual feature of the plastome is the loss of two genes, infA and rpl32. Examination of nuclear transcriptomes indicates that both of these genes have been transferred to the nucleus. Comparing the plastome sequence of Thalictrum with the two other Ranunculaceae and the survey of the rpl32 gene loss resolve the phylogenetic distribution and timing of this gene loss/transfer event in Ranunculaceae.
Plant material, plastid isolation, and RCA
Fresh leaf tissue of Thalictrum coreanum was sampled from a single individual from a natural population in Gangwon-do, Korea. Intact plastids were isolated from 1.45 g of tissue using the sucrose step gradient method of Jansen et al. . Isolated plastids were used to amplify the plastid genome by rolling circle amplification (RCA) using REPLI-g midi Kit (cat. No. 150043, Qiagen, Valencia, CA, USA) following the protocol described in Jansen et al. . RCA products were digested with EcoRI and the resulting fragments were separated by gel electrophoresis in a 1% agarose gel to verify the purity and quantity of plastid DNA.
Genome sequencing, assembly, annotation, and analyses
Plastid DNA (538.9 ng/ul) was sheared by nebulization, subjected to library preparation and sequencing on a Roche 454 Genome Sequencer (GS) FLX Titanium platform at Solgent Co. (Deajeon, Korea). The Roche 454 sequencing produced approximately 80 Mb of sequence with an average read length of 357 bp.
The quality filtered sequence reads were assembled using the GS de novo sequence assembler v.2.5.3 (Roche 454 Life Sciences, Branford, CT, USA) and multiple assemblies were performed with modified parameters (i.e. adjusting minimum overlap length). Three long contigs representing a nearly complete plastid genome sequence were generated and the contigs were mapped against two complete plastid genomes of Ranunculaceae, Megaleranthis saniculifolia (NC_012615) and Ranunculus macranthus (NC_008796), in Geneious R6 v.6.1.6 . The presence of gaps between the junctions of LSC, SSC, and IR regions were filled by polymerase chain reaction (PCR) and Sanger sequencing. The Roche 454 pyrosequencing platform is known to have a high error rate in long homopolymer regions [44,45]. There were 36 homopolymers > 7 bp in protein-coding genes, five of which were nonsense mutations, and these regions were corrected by PCR and Sanger sequencing. All primers for PCR were designed by Primer3  in Geneious R6 (Additional file 2: Table S1).
Annotation of plastid genome was done in DOGMA  and all tRNA genes were verified by their predicted secondary structures using tRNAscan-SE 1.3.1 . A genome map was drawn with OGDRAW . The plastome sequence of T. coreanum was deposited in GenBank (accession number KM206568).
Two published plastomes of Ranunculaceae [31,32], M. saniculifolia and R. macranthus, were used for genomic comparisons with T. coreanum. Whole genome alignment was performed under ‘progressiveMauve algorithm’  in Geneious R6. Repetitive sequences were identified by performing BLASTN v.2.2.28+ (word size = 11) searches of each plastome against itself with an e-value cutoff of 1e−10 and at least 90% sequence identity. The analysis was performed on Lonestar Dell Linux Cluster of the Texas Advanced Computing Center (TACC).
Identification of gene transfers to the nucleus
Three genera of Ranunculaceae (Aquilegia, Hydrastis, and Thalictrum) with T-type chromosomes were surveyed for gene transfers to the nucleus. Thalictrum thalictroides and H. canadensis transcriptomes from the 1KP project database  and A. coerulea transcriptome from the genomics portal Phytozome v.10  were searched. Transferred genes were identified using BlastN of the infA and rpl32 sequences from the M. saniculifolia and R. macranthus plastomes against the transcriptomes. The NCBI Conserved Domain Database (CDD) was used for functional domain annotation . TargetP v.1.1  and Predotar v.1.03  were used to predict transit peptides. Putative ORFs were searched for using Phytozome with BLASTX and ‘angiosperms’ as a reference sequence source to identify plant gene families. Nucleotide and amino acid sequences of nuclear and plastid genes were aligned with MUSCLE  in Geneious R6.
Survey for loss of rpl32 gene in the subfamily Thalictroideae
Seventeen Thalictrum species from all major clades of the phylogenetic tree of the genus (S. Park, unpublished) and six other genera of the subfamily Thalictroideae were sampled (Additional file 2: Table S2). Total genomic DNA was isolated from either fresh leaves or herbarium specimens using the methods of Allen et al.  with the following modifications to the extraction buffer: Cetyl trimethylammonium bromide (CTAB) was increased to 3%; and 1% polyvinylpyrrolidone (PVP, w/v, MW 4,000) and 2% beta-mercaptoethanol (Sigma, St. Louis, MO) were added. To detect the rpl32 gene, the intergenic spacer (IGS) region between ndhF and trnL-UAG genes was amplified by PCR using the Shaw et al.  primers (ndhF: GAAAGGTATKATCCAYGMATATT and trnL-UAG: CTGCTTCCTAAGAGCAGCGT). PCR products were purified by using Solg™ Gel & PCR Purification System Kit (Solgent Co., Daejeon, Korea) following the manufacturer’s protocol. All sequencing of PCR products was performed using an ABI 3730XL DNA Analyzer (Applied Biosystems, California, USA) at Solgent Co., and nucleotide sequences were aligned with MUSCLE in Geneious R6. Statistical analysis was conducted by using R v.2.1.5  to test whether gene loss/transfer was associated with the size of intergenic spacer.
Phylogenetic analyses were performed on two data sets. The first included 39 species with nucleotide sequence of five plastid genes/regions (rbcL, ndhF, ndhA intron, trnL intron, and trnL-F intergenic spacer), including 31 Thalictrum species and a single species from six other genera of the subfamily Thalictroideae (Additional file 2: Table S2). Megaleranthis saniculifolia and R. macranthus were used as outgroups by extracting nucleotide sequences of the five genes/regions from the published plastomes. The second data set included sequences of the plastid-encoded rpl32 gene for 48 taxa and four nuclear-encoded copies (Additional file 2: Table S3). The data sets were aligned with MUSCLE in Geneious R6. Maximum likelihood (ML) analyses were performed with RAxML v.7.2.8  using the ‘GTRGAMMA’ model under the rapid bootstrap algorithm with 1000 replicates at TACC.
Estimating nucleotide substitution rates
To analyze rates of nucleotide substitution, photosystem I (psaA, B, and C) and II (psbB, C, D, E, F, H, J, L, M, N, T, and Z) genes and rbcL were sampled from selected angiosperms (Additional file 2: Table S3). The data were concatenated into a single data set and a phylogenetic tree was generated using the ML method (see phylogenetic analyses section) and used as a constraint tree (Additional file 1: Figure S3) for all rate comparisons. Nonsynonymous (d N ) and synonymous (d S ) substitution rates for 48 plastid-encoded rpl32 sequences and three nuclear-encoded sequences (rpl32 from Bruguiera was not used in the rate variation estimation because there were insufficient plastid data to generate a constraint tree) were calculated in PAML v.4.8  using codeml option with codon frequencies estimated with an F3 × 4 model.
Availability of supporting data
The data sets supporting the results of this article are included within additional files. The phylogenetic data sets (including amino acid sequence) supporting the results of this article are available in Dryad Digital Repository (http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.g84g5).
Nuclear DNAs of plastid origin
Open reading frame
- d N :
Number of substitutions per nonsynonymous site
- d S :
Number of substitutions per synonymous site
Rolling circle amplification
Million years ago
Texas advanced computing center
Timmis JN, Ayliffe MA, Huang CY, Martin W. Endosymbiotic gene transfer: organelle genomes forge eukaryotic chromosomes. Nature Rev Genet. 2004;5:123–35.
Bock R. Structure, function, and inheritance of plastid genomes. In: Bock R, editor. Cell and Molecular Biology of Plastids, vol. 19. Berlin Heidelberg: Springer; 2007. p. 29–63.
Jansen RK, Ruhlman TA. Plastid genomes of seed plants. In: Bock R, Knoop V, editors. Genomics of Chloroplasts and Mitochondria. New York: Springer; 2012. p. 103–26.
Stegemann S, Hartmann S, Ruf S, Bock R. High-frequency gene transfer from the chloroplast genome to the nucleus. Proc Natl Acad Sci USA. 2003;100:8828–33.
Glaser E, Soll J. Targeting signals and import machinery of plastids and plant mitochondria. In: Daniell H, Chase C, editors. Molecular Biology and Biotechnology of Plant Organelles: Chloroplasts and Mitochondria. New York: Springer; 2004. p. 385–418.
Bock R, Timmis JN. Reconstructing evolution: gene transfer from plastids to the nucleus. BioEssays. 2008;30:556–66.
Millen RS, Olmstead RG, Adams KL, Palmer JD, Lao NT, Heggie L, et al. Many parallel losses of infA from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus. Plant Cell. 2001;13:645–58.
Gantt JS, Baldauf SL, Calie PJ, Weeden NF, Palmer JD. Transfer of rpl22 to the nucleus greatly preceded its loss from the chloroplast and involved the gain of an intron. EMBO J. 1991;10:3073–8.
Jansen RK, Saski C, Lee SB, Hansen AK, Daniell H. Complete plastid genome sequences of three rosids (Castanea, Prunus, Theobroma): evidence for at least two independent transfers of rpl22 to the nucleus. Mol Biol Evol. 2011;28:835–47.
Cusack BP, Wolfe KH. When gene marriages don't work out: divorce by subfunctionalization. Trends Genet. 2007;23:270–2.
Ueda M, Fujimoto M, Arimura S, Murata J, Tsutsumi N, Kadowaki K. 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:51–6.
Magee AM, Aspinall S, Rice DW, Cusack BP, Sémon M, Perry AS, et al. Localized hypermutation and associated gene losses in legume chloroplast genomes. Genome Res. 2010;20:1700–10.
Sabir J, Schwarz E, Ellison N, Zhang J, Baeshen NA, Mutwakil M, Jansen RK, Ruhlman TA. Evolutionary and biotechnology implications of plastid genome variation in the inverted-repeat-lacking clade of legumes. Plant Biotechnol J 2014, doi:10.1111/pbi.12179.
Noutsos C, Richly E, Leister D. Generation and evolutionary fate of insertions of organelle DNA in the nuclear genomes of flowering plants. Genome Res. 2005;15:616–28.
Richly E, Leister D. NUPTs in sequenced eukaryotes and their genomic organization in relation to NUMTs. Mol Biol Evol. 2004;21:1972–80.
Michalovova M, Vyskot B, Kejnovsky E. Analysis of plastid and mitochondrial DNA insertions in the nucleus (NUPTs and NUMTs) of six plant species: size, relative age and chromosomal localization. Heredity. 2013;111:314–20.
Tamura M. Angiospermae. Ordnung Ranunculales. Fam. Ranunculaceae. II. Systematic Part. In: Hiepko P, editor. Die natürliche Pflanzenfamilien, 17aIV. 2nd ed. Berlin: Duncker and Humblot; 1995. p. 223–519.
Langlet OFI. Beitrage zur zytologie der Ranunculazeen. Sven Bot Tidskr. 1927;21:1–17.
Hsiao PK. Ranunculaceae. In: Institutum Botanicum Academiae Sinicae, editor. Flora Reipublicae Popularis Sinicae, vol. 27. Beijing: Science Press; 1979. p. 465–601.
Fu DZ. Phylogenetic considerations on the subfamily Thalictroideae (Ranunculaceae). Cathaya. 1990;2:181–90.
Jensen U, Hoot SB, Johansson JT, Kosuge K. Systematics and phylogeny of the Ranunculaceae - a revised family concept on the basis of molecular data. Plant Syst Evol. 1995;9 Suppl 9:273–80.
Wang W, Chen ZD. Generic level phylogeny of Thalictroideae (Ranunculaceae) - Implications for the taxonomic status of Paropyrum and petal evolution. Taxon. 2007;56:811–21.
Wang W, Lu AM, Ren Y, Endress ME, Chen ZD. Phylogeny and classification of Ranunculales: evidence from four molecular loci and morphological data. Perspect Plant Ecol Evol Syst. 2009;11:81–110.
Ro KE, Keener CS, McPheron BA. Molecular phylogenetic study of the Ranunculaceae: utility of the nuclear 26S ribosomal DNA in inferring intrafamilial relationships. Mol Phylogenet Evol. 1997;8:117–27.
Soza VL, Brunet J, Liston A, Smith PS, Di Stilio VS. Phylogenetic insights into the correlates of dioecy in meadow-rues (Thalictrum, Ranunculaceae). Mol Phylogenet Evol. 2012;63:180–92.
Soza VL, Haworth KL, Di Stilio VS. Timing and consequences of recurrent polyploidy in meadow-rues (Thalictrum, Ranunculaceae). Mol Biol Evol. 2013;30:1940–54.
Chen G, Ramachandran C, Krishan A. Thaliblastine, a plant alkaloid, circumvents multidrug resistance by direct binding to P-glycoprotein. Cancer Res. 1993;53:2544–7.
Pajeva I, Todorov DK, Seydel J. Membrane effects of the antitumor drugs doxorubicin and thaliblastine: comparison to multidrug resistance modulators verapanil and trans-flupentixol. Eur J Pharm Sci. 2004;21:243–50.
Hoot SB, Palmer JD. Structural rearrangements, including parallel inversions, within the chloroplast genome of Anemone and related genera. J Mol Evol. 1994;38:274–81.
Johansson JT. There large inversions in the chloroplast genomes and one loss of the chloroplast gene rps16 suggest an early evolutionary split in the genus Adonis (Ranunculaceae). Plant Syst Evol. 1999;218:133–43.
Raubeson LA, Peery R, Chumley TW, Dziubek C, Fourcade HM, Boore JL, et al. Comparative chloroplast genomics: analyses including new sequences from the angiosperms Nuphar advena and Ranunculus macranthus. BMC Genomics. 2007;8:174.
Kim YK, Park CW, Kim KJ. Complete chloroplast DNA sequence from a Korean endemic genus, Megaleranthis saniculifolia, and its evolutionary implications. Mol Cells. 2009;27:365–81.
Jansen RK, Cai Z, Raubeson LA, Daniell H, de Pamphilis CW, Leebens-Mack J, et al. 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–74.
Wicke S, Schneeweiss GM, de Pamphilis 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:273–97.
Wu CS, Lai YT, Lin CP, Wang YN, Chaw SM. Evolution of reduced and compact chloroplast genomes (cpDNAs) in gnetophytes: selection toward a lower-cost strategy. Mol Phylogenet Evol. 2009;52:115–24.
Wu CS, Chaw SM. Highly rearranged and size-variable chloroplast genomes in conifers II clade (cupressophytes): evolution towards shorter intergenic spacers. Plant Biotechnol J. 2014;12:344–53.
Fior S, Li M, Oxelman B, Viola R, Hodges SA, Ometto L, et al. Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions. New Phytol. 2013;198:579–92.
Bastida JM, Alcantara JM, Rey PJ, Vargas P, Herrera CM. Extended phylogeny of Aquilegia: the biogeographical and ecological patterns of two simultaneous but contrasting radiations. Plant Syst Evol. 2010;284:171–85.
Hoot SB. Phylogeny of the Ranunculaceae based on preliminary atpB, rbcL and 18S nuclear ribosomal DNA sequence data. Plant Syst Evol. 1995;9 Suppl 9:241–51.
Ruhlman TA, Jansen RK. The plastid genomes of flowering plants. In: Maliga P, editor. Chloroplast Biotechnology: Methods and Protocols Vol. 1132. New York: Springer Science and Business Media, LLC; 2014. p. 3–38.
Weng ML, Blazier JC, Govindu M, Jansen RK. Reconstruction of the ancestral plastid genome in Geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Mol Biol Evol. 2014;31:645–59.
Jansen RK, Boore JL, de Pamphilis CW, Chumley TW, Haberle RC, Wyman SK, et al. Methods for obtaining and analyzing whole chloroplast genome sequences. Methods Enzymol. 2005;395:348–84.
Biomatters: Geneious R6 v.6.1.6 [http://www.geneious.com/]
Moore MJ, Dhingra A, Soltis PS, Shaw R, Farmerie WG, Folta KM, et al. Rapid and accurate pyrosequencing of angiosperm plastid genomes. BMC Plant Biol. 2006;6:17.
Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM. Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007;8:R143.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3 - new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.
Wyman SK, Jansen RK, Boore JL. Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004;20:3252–5.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.
Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW–a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41:W575–81.
Darling AC, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.
1KP project database [http://onekp.com/project.html]
Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, et al. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012;40:D1178–86.
Marchler-Bauer A, Zheng C, Chitsaz F, Derbyshire MK, Geer LY, Geer RC, et al. CDD: conserved domains and protein three-dimensional structure. Nucleic Acids Res. 2013;41:D348–52.
Emanuelsson O, Brunak S, von Heijne G, Nielsen H. Locating proteins in the cell using TargetP, SignalP, and related tools. Nature Protocols. 2007;2:953–71.
Small I, Peeters N, Legeai F, Lurin C. Predotar: a tool for rapidly screening proteomes for N-terminal targeting sequences. Proteomics. 2004;4:1581–90.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Allen GC, Flores-Vergara MA, Krasynanski S, Kumar S, Thompson WF. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nature Protocols. 2006;1:2320–5.
Shaw J, Lickey EB, Schilling EE, Small RL. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. Am J Bot. 2007;94:275–88.
R Development Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: The R Foundation for Statistical Computing; 2014.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.
This paper represents a portion of SP’s Ph.D. thesis in Department of Biology at Yeungnam University, Korea. The authors thank the Texas Advanced Computing Center (TACC) at the University of Texas at Austin for access to supercomputers, Boram An for assistance on colleting material and plastid isolation. Support was provided by the National Science Foundation (IOS-1027259) and the S. F. Blake Centennial Professorship to RJK. We also thank Chris Blazier, Mao-Lun Weng, Jin Zhang and three anonymous reviewers for valuable comments on an earlier version of the manuscript.
The authors declare that they have no competing interests.
SP performed all analyses, prepared the figures and tables, and drafted the manuscript. RKJ and SJP contributed to the design of the project and assisted with manuscript preparation. All authors read and approved the final draft of the manuscript.
Maximum likelihood phylogenetic tree of 52 taxa based on rpl32 gene sequence. Figure S2. Alignment of the infA gene. Figure S3. Maximum likelihood phylogenetic tree inferred from 16 genes from 51 taxa of angiosperms used as a constraint tree for rate comparisons.
Primers used for amplification. Table S2. NCBI accession numbers for species included in the phylogenetic analysis based on five plastid regions. Table S3. Taxon sampling and NCBI accession numbers for phylogenetic analyses, estimating nucleotide substitution rates, and examination of the correlation between reduction of ndhF-trnL intergenic spacer and rpl32 gene loss.
About this article
Cite this article
Park, S., Jansen, R.K. & 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 15, 40 (2015) doi:10.1186/s12870-015-0432-6
- Gene loss
- Intracellular gene transfer
- Plastid genome