Plastid Genome Evolution Research Articles

Markov-Modulated Continuous-Time Markov Chains to Identify Site- and Branch-Specific Evolutionary Variation in BEAST

Markov models of character substitution on phylogenies form the foundation of phylogenetic inference frameworks. Early models made the simplifying assumption that the substitution process is homogeneous over time and across sites in the molecular sequence alignment. While standard practice adopts extensions that accommodate heterogeneity of substitution rates across sites, heterogeneity in the process over time in a site-specific manner remains frequently overlooked. This is problematic, as evolutionary processes that act at the molecular level are highly variable, subjecting different sites to different selective constraints over time, impacting their substitution behavior. We propose incorporating time variability through Markov-modulated models (MMMs), which extend covarion-like models and allow the substitution process (including relative character exchange rates as well as the overall substitution rate) at individual sites to vary across lineages. We implement a general MMM framework in BEAST, a popular Bayesian phylogenetic inference software package, allowing researchers to compose a wide range of MMMs through flexible XML specification. Using examples from bacterial, viral, and plastid genome evolution, we show that MMMs impact phylogenetic tree estimation and can substantially improve model fit compared to standard substitution models. Through simulations, we show that marginal likelihood estimation accurately identifies the generative model and does not systematically prefer the more parameter-rich MMMs. To mitigate the increased computational demands associated with MMMs, our implementation exploits recent developments in BEAGLE, a high-performance computational library for phylogenetic inference. [Bayesian inference; BEAGLE; BEAST; covarion, heterotachy; Markov-modulated models; phylogenetics.]

Plastid Genome Evolution in the Subtribe Calypsoinae (Epidendroideae, Orchidaceae)

Calypsoinae is a small subtribe in Orchidaceae (Epidendroideae) characterized by diverse trophic strategies and morphological characters. Calypsoinae includes 13 genera, four of which are leafless and mycoheterotrophic. Mycoheterotrophic species in the leafless genus Corallorhiza are well suited to studies of plastome evolution. However, the lack of plastome sequences for other genera in Calypsoinae limits the scope of comparative and phylogenetic analyses, in particular our understanding of plastome evolution. To understand plastid genome evolution in Calypsoinae, we newly sequenced the plastomes of 12 species in the subtribe, including representatives of three mycoheterotrophic genera as well as five autotrophic genera. We detected two parallel photosynthetic losses in Corallorhiza. Evolutionary analyses indicated that the transition to obligate mycoheterotrophy leads to the relaxation of selection in a highly gene-specific pattern.

Extensive plastome reduction and loss of photosynthesis genes in Diphelypaea coccinea, a holoparasitic plant of the family Orobanchaceae.

BackgroundParasitic plants have the ability to obtain nutrients from their hosts and are less dependent on their own photosynthesis or completely lose this capacity. The reduction in plastid genome size and gene content in parasitic plants predominantly results from loss of photosynthetic genes. Plants from the family Orobanchaceae are used as models for studying plastid genome evolution in the transition from an autotrophic to parasitic lifestyle. Diphelypaea is a poorly studied genus of the Orobanchaceae, comprising two species of non-photosynthetic root holoparasites. In this study, we sequenced the plastid genome of Diphelypaea coccinea and compared it with other Orobanchaceae, to elucidate patterns of plastid genome evolution. In addition, we used plastid genome data to define the phylogenetic position of Diphelypaea spp.MethodsThe complete nucleotide sequence of the plastid genome of D. coccinea was obtained from total plant DNA, using pyrosequencing technology.ResultsThe D. coccinea plastome is only 66,616 bp in length, and is highly rearranged; however, it retains a quadripartite structure. It contains only four rRNA genes, 25 tRNA genes and 25 protein-coding genes, being one of the most highly reduced plastomes among the parasitic Orobanchaceae. All genes related to photosynthesis, including the ATP synthase genes, had been lost, whereas most housekeeping genes remain intact. The plastome contains two divergent, but probably intact clpP genes. Intron loss had occurred in some protein-coding and tRNA genes. Phylogenetic analysis yielded a fully resolved tree for the Orobanchaceae, with Diphelypaea being a sister group to Orobanche sect. Orobanche.

Comparative plastid genomics of Pinus species: Insights into sequence variations and phylogenetic relationships

AbstractPinus L. is the largest genus of conifers and provides a classical model for studying species divergence and phylogenetic evolution by gymnosperms. However, our poor understanding of sequence divergence in the whole plastid genomes of Pinus species severely hinders studies of their evolution and phylogeny. Thus, we analyzed the sequences of 97 Pinus plastid genomes, including four newly sequenced genomes and 93 previously published plastomes, to explore the evolution and phylogenetic relationships in the genus Pinus. The complete chloroplast genomes of Pinus species ranged in size from 109 640 bp (P. cembra L.) to 121 976 bp (P. glabra Walter), and these genomes comprised circular DNA molecules in a similar manner to those of most gymnosperms. We identified 9108 repeats where most of the repeats comprised the dispersed type with 3983 (44%), followed by tandem repeats with 2999 (33%), and then palindromic repeats with 2126 (23%). Sixteen divergence hotspot regions were identified in Pinus plastid genomes, which could be useful molecular markers for future population genetics studies. Phylogenetic analysis showed that Pinus species could be divided into two diverged clades comprising the subgenera Strobus (single needle section) and Pinus (double needles section). Molecular dating suggested that the genus Pinus originated approximately 130.38 Mya during the late Cretaceous. The two subgenera subsequently split 85.86 Mya, which was largely consistent with the other molecular results based on partial DNA markers. These findings provide important insights into the sequence variations and phylogenetic evolution of Pinus plastid genomes.

Celebrating Research Devoted to Seed‐Free Land Plants

Celebrating Research Devoted to Seed‐Free Land Plants

Evolution of plastid genomes of Holcoglossum (Orchidaceae) with recent radiation

BackgroundThe plastid is a semiautonomous organelle with its own genome. Plastid genomes have been widely used as models for studying phylogeny, speciation and adaptive evolution. However, most studies focus on comparisons of plastid genome evolution at high taxonomic levels, and comparative studies of the process of plastome evolution at the infrageneric or intraspecific level remain elusive. Holcoglossum is a small genus of Orchidaceae, consisting of approximately 20 species of recent radiation. This made it an ideal group to explore the plastome mutation mode at the infrageneric or intraspecific level.ResultsIn this paper, we reported 15 complete plastid genomes from 12 species of Holcoglossum and 1 species of Vanda. The plastid genomes of Holcoglossum have a total length range between 145 kb and 148 kb, encoding a set of 102 genes. The whole set of ndh-gene families in Holcoglossum have been truncated or pseudogenized. Hairpin inversion in the coding region of the plastid gene ycf2 has been found.ConclusionsUsing a comprehensive comparative plastome analysis, we found that all the indels between different individuals of the same species resulted from the copy number variation of the short repeat sequence, which may be caused by replication slippage. Annotation of tandem repeats shows that the variation introduced by tandem repeats is widespread in plastid genomes. The hairpin inversion found in the plastid gene ycf2 occurred randomly in the Orchidaceae.

Plastid Genome Degradation in the Endangered, Mycoheterotrophic, North American Orchid Hexalectris warnockii.

Heterotrophic plants provide evolutionarily independent, natural experiments in the genomic consequences of radically altered nutritional regimes. Here, we have sequenced and annotated the plastid genome of the endangered mycoheterotrophic orchid Hexalectris warnockii. This orchid bears a plastid genome that is ∼80% the total length of the leafy, photosynthetic Phalaenopsis, and contains just over half the number of putatively functional genes of the latter. The plastid genome of H. warnockii bears pseudogenes and has experienced losses of genes encoding proteins directly (e.g., psa/psb, rbcL) and indirectly involved in photosynthesis (atp genes), suggesting it has progressed beyond the initial stages of plastome degradation, based on previous models of plastid genome evolution. Several dispersed and tandem repeats were detected, that are potentially useful as conservation genetic markers. In addition, a 29-kb inversion and a significant contraction of the inverted repeat boundaries are observed in this plastome. The Hexalectris warnockii plastid genome adds to a growing body of data useful in refining evolutionary models in parasites, and provides a resource for conservation studies in these endangered orchids.

Dense infraspecific sampling reveals rapid and independent trajectories of plastome degradation in a heterotrophic orchid complex.

Heterotrophic plants provide excellent opportunities to study the effects of altered selective regimes on genome evolution. Plastid genome (plastome) studies in heterotrophic plants are often based on one or a few highly divergent species or sequences as representatives of an entire lineage, thus missing important evolutionary-transitory events. Here, we present the first infraspecific analysis of plastome evolution in any heterotrophic plant. By combining genome skimming and targeted sequence capture, we address hypotheses on the degree and rate of plastome degradation in a complex of leafless orchids (Corallorhiza striata) across its geographic range. Plastomes provide strong support for relationships and evidence of reciprocal monophyly between C.involuta and the endangered C.bentleyi. Plastome degradation is extensive, occurring rapidly over a few million years, with evidence of differing rates of genomic change among the two principal clades of the complex. Genome skimming and targeted sequence capture differ widely in coverage depth overall, with depth in targeted sequence capture datasets varying immensely across the plastome as a function of GC content. These findings will help to fill a knowledge gap in models of heterotrophic plastid genome evolution, and have implications for future studies in heterotrophs.

Evolution of Geosiris (Iridaceae): historical biogeography and plastid-genome evolution in a genus of non-photosynthetic tropical rainforest herbs disjunct across the Indian Ocean

Mycoheterotrophs, i.e. plants that acquire carbon from root-associated soil fungi, often have highly degraded plastomes, reflecting relaxed selective constraints on plastid genes following the loss of photosynthesis. Geosiris Baill. is the only mycoheterotrophic genus in Iridaceae and comprises two species in Madagascar and nearby islands, and a third recently discovered species in north-eastern Australia. Here, we characterise the plastomes of the Australian and one Madagascan species to compare patterns of plastome degradation in relation to autotrophic and other mycoheterotrophic taxa and investigate the evolutionary and biogeographical history of the genus in Iridaceae. Both examined species have lost approximately half their plastid-encoded genes and a small but significant reduction in purifying selection in retained non-photosynthetic genes was observed. Geosiris is confirmed as monophyletic, with initial divergence of the genus occurring c. 53 million years ago, and subsequent diversification occurring c. 30 million years ago. Africa (including Madagascar) is reconstructed as the most likely ancestral area of the genus, implying a major range-expansion event of one lineage to Australia after its divergence in the Oligocene. Our study has highlighted the dynamic evolutionary history of Geosiris, contributed to the characterisation of mycoheterotrophic plastomes, and furthered our understanding of plastome structure and function.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

Plastome Evolution in the Sole Hemiparasitic Genus Laurel Dodder (Cassytha) and Insights into the Plastid Phylogenomics of Lauraceae.

To date, little is known about the evolution of plastid genomes (plastomes) in Lauraceae. As one of the top five largest families in tropical forests, the Lauraceae contain many species that are important ecologically and economically. Lauraceous species also provide wonderful materials to study the evolutionary trajectory in response to parasitism because they contain both nonparasitic and parasitic species. This study compared the plastomes of nine Lauraceous species, including the sole hemiparasitic and herbaceous genus Cassytha (laurel dodder; here represented by Cassytha filiformis). We found differential contractions of the canonical inverted repeat (IR), resulting in two IR types present in Lauraceae. These two IR types reinforce Cryptocaryeae and Neocinnamomum—Perseeae–Laureae as two separate clades. Our data reveal several traits unique to Cas. filiformis, including loss of IRs, loss or pseudogenization of 11 ndh and rpl23 genes, richness of repeats, and accelerated rates of nucleotide substitutions in protein-coding genes. Although Cas. filiformis is low in chlorophyll content, our analysis based on dN/dS ratios suggests that both its plastid house-keeping and photosynthetic genes are under strong selective constraints. Hence, we propose that short generation time and herbaceous lifestyle rather than reduced photosynthetic ability drive the accelerated rates of nucleotide substitutions in Cas. filiformis.

AccD nuclear transfer of Platycodon grandiflorum and the plastid of early Campanulaceae

BackgroundCampanulaceae species are known to have highly rearranged plastid genomes lacking the acetyl-CoA carboxylase (ACC) subunit D gene (accD), and instead have a nuclear (nr)-accD. Plastid genome information has been thought to depend on studies concerning Trachelium caeruleum and genome announcements for Adenophora remotiflora, Campanula takesimana, and Hanabusaya asiatica. RNA editing information for plastid genes is currently unavailable for Campanulaceae. To understand plastid genome evolution in Campanulaceae, we have sequenced and characterized the chloroplast (cp) genome and nr-accD of Platycodon grandiflorum, a basal member of Campanulaceae.ResultsWe sequenced the 171,818 bp cp genome containing a 79,061 bp large single-copy (LSC) region, a 42,433 bp inverted repeat (IR) and a 7840 bp small single-copy (SSC) region, which represents the cp genome with the largest IR among species of Campanulaceae. The genome contains 110 genes and 18 introns, comprising 77 protein-coding genes, four RNA genes, 29 tRNA genes, 17 group II introns, and one group I intron. RNA editing of genes was detected in 18 sites of 14 protein-coding genes. Platycodon has an IR containing a 3′ rps12 operon, which occurs in the middle of the LSC region in four other species of Campanulaceae (T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica), but lacks accD, clpP, infA, and rpl23, as has been found in these four species. Platycodon nr-accD contains about 3.2 kb intron between nr-accD.e1 and nr-accD.e2 at the same insertion point as in other Campanulaceae. The phylogenies of the plastid genomes and accD show that Platycodon is basal in the Campanulaceae clade, indicating that IR disruption in Campanulaceae occurred after the loss of accD, clpP, infA, and rpl23 in the cp genome, which occurred during plastid evolution in Campanulaceae.ConclusionsThe plastid genome of P. grandiflorum lacks the rearrangement of the IR found in T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica. The absence of accD, clpP, infA, and rpl23 in the plastid genome is a synapomorphic characteristic of Campanulaceae. The chloroplast genome phylogeny supports the hypothesis that chloroplast genomic arrangement occurred after accD nuclear transfer and loss of the four genes in the plastid of early Campanulaceae as a lineage of asterids.

Contrasting Patterns of Nucleotide Substitution Rates Provide Insight into Dynamic Evolution of Plastid and Mitochondrial Genomes of Geranium.

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid.

Positive Selection in Rapidly Evolving Plastid-Nuclear Enzyme Complexes.

Rates of sequence evolution in plastid genomes are generally low, but numerous angiosperm lineages exhibit accelerated evolutionary rates in similar subsets of plastid genes. These genes include clpP1 and accD, which encode components of the caseinolytic protease (CLP) and acetyl-coA carboxylase (ACCase) complexes, respectively. Whether these extreme and repeated accelerations in rates of plastid genome evolution result from adaptive change in proteins (i.e., positive selection) or simply a loss of functional constraint (i.e., relaxed purifying selection) is a source of ongoing controversy. To address this, we have taken advantage of the multiple independent accelerations that have occurred within the genus Silene (Caryophyllaceae) by examining phylogenetic and population genetic variation in the nuclear genes that encode subunits of the CLP and ACCase complexes. We found that, in species with accelerated plastid genome evolution, the nuclear-encoded subunits in the CLP and ACCase complexes are also evolving rapidly, especially those involved in direct physical interactions with plastid-encoded proteins. A massive excess of nonsynonymous substitutions between species relative to levels of intraspecific polymorphism indicated a history of strong positive selection (particularly in CLP genes). Interestingly, however, some species are likely undergoing loss of the native (heteromeric) plastid ACCase and putative functional replacement by a duplicated cytosolic (homomeric) ACCase. Overall, the patterns of molecular evolution in these plastid-nuclear complexes are unusual for anciently conserved enzymes. They instead resemble cases of antagonistic coevolution between pathogens and host immune genes. We discuss a possible role of plastid-nuclear conflict as a novel cause of accelerated evolution.

Coevolution between Nuclear-Encoded DNA Replication, Recombination, and Repair Genes and Plastid Genome Complexity.

Disruption of DNA replication, recombination, and repair (DNA-RRR) systems has been hypothesized to cause highly elevated nucleotide substitution rates and genome rearrangements in the plastids of angiosperms, but this theory remains untested. To investigate nuclear–plastid genome (plastome) coevolution in Geraniaceae, four different measures of plastome complexity (rearrangements, repeats, nucleotide insertions/deletions, and substitution rates) were evaluated along with substitution rates of 12 nuclear-encoded, plastid-targeted DNA-RRR genes from 27 Geraniales species. Significant correlations were detected for nonsynonymous (dN) but not synonymous (dS) substitution rates for three DNA-RRR genes (uvrB/C, why1, and gyrA) supporting a role for these genes in accelerated plastid genome evolution in Geraniaceae. Furthermore, correlation between dN of uvrB/C and plastome complexity suggests the presence of nucleotide excision repair system in plastids. Significant correlations were also detected between plastome complexity and 13 of the 90 nuclear-encoded organelle-targeted genes investigated. Comparisons revealed significant acceleration of dN in plastid-targeted genes of Geraniales relative to Brassicales suggesting this correlation may be an artifact of elevated rates in this gene set in Geraniaceae. Correlation between dN of plastid-targeted DNA-RRR genes and plastome complexity supports the hypothesis that the aberrant patterns in angiosperm plastome evolution could be caused by dysfunction in DNA-RRR systems.

Diversity of transcripts and transcript processing forms in plastids of the dinoflagellate alga Karenia mikimotoi.

Plastids produce a vast diversity of transcripts. These include mature transcripts containing coding sequences, and their processing precursors, as well as transcripts that lack direct coding functions, such as antisense transcripts. Although plastid transcriptomes have been characterised for many plant species, less is known about the transcripts produced in other plastid lineages. We characterised the transcripts produced in the fucoxanthin-containing plastids of the dinoflagellate alga Karenia mikimotoi. This plastid lineage, acquired through tertiary endosymbiosis, utilises transcript processing pathways that are very different from those found in plants and green algae, including 3′ poly(U) tail addition, and extensive substitutional editing of transcript sequences. We have sequenced the plastid transcriptome of K. mikimotoi, and have detected evidence for divergent evolution of fucoxanthin plastid genomes. We have additionally characterised polycistronic and monocistronic transcripts from two plastid loci, psbD-tRNAMet-ycf4 and rpl36-rps13-rps11. We find evidence for a range of transcripts produced from each locus that differ in terms of editing state, 5′ end cleavage position, and poly(U) tail addition. Finally, we identify antisense transcripts in K. mikimotoi, which appear to undergo different processing events from the corresponding sense transcripts. Overall, our study provides insights into the diversity of transcripts and processing intermediates found in plastid lineages across the eukaryotes.Electronic supplementary materialThe online version of this article (doi:10.1007/s11103-015-0408-9) contains supplementary material, which is available to authorized users.

Direct Estimate of the Spontaneous Mutation Rate Uncovers the Effects of Drift and Recombination in the Chlamydomonas reinhardtii Plastid Genome.

Plastids perform crucial cellular functions, including photosynthesis, across a wide variety of eukaryotes. Since endosymbiosis, plastids have maintained independent genomes that now display a wide diversity of gene content, genome structure, gene regulation mechanisms, and transmission modes. The evolution of plastid genomes depends on an input of de novo mutation, but our knowledge of mutation in the plastid is limited to indirect inference from patterns of DNA divergence between species. Here, we use a mutation accumulation experiment, where selection acting on mutations is rendered ineffective, combined with whole-plastid genome sequencing to directly characterize de novo mutation in Chlamydomonas reinhardtii. We show that the mutation rates of the plastid and nuclear genomes are similar, but that the base spectra of mutations differ significantly. We integrate our measure of the mutation rate with a population genomic data set of 20 individuals, and show that the plastid genome is subject to substantially stronger genetic drift than the nuclear genome. We also show that high levels of linkage disequilibrium in the plastid genome are not due to restricted recombination, but are instead a consequence of increased genetic drift. One likely explanation for increased drift in the plastid genome is that there are stronger effects of genetic hitchhiking. The presence of recombination in the plastid is consistent with laboratory studies in C. reinhardtii and demonstrates that although the plastid genome is thought to be uniparentally inherited, it recombines in nature at a rate similar to the nuclear genome.

Complete plastome sequence of Thalictrum coreanum (Ranunculaceae) and transfer of the rpl32 gene to the nucleus in the ancestor of the subfamily Thalictroideae.

BackgroundPlastids 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.ResultsThe 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.ConclusionsThe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0432-6) contains supplementary material, which is available to authorized users.

Correlation between sequence divergence and polymorphism reveals similar evolutionary mechanisms acting across multiple timescales in a rapidly evolving plastid genome.

BackgroundAlthough the plastid genome is highly conserved across most angiosperms, multiple lineages have increased rates of structural rearrangement and nucleotide substitution. These lineages exhibit an excess of nonsynonymous substitutions (i.e., elevated dN/dS ratios) in similar subsets of plastid genes, suggesting that similar mechanisms may be leading to relaxed and/or positive selection on these genes. However, little is known regarding whether these mechanisms continue to shape sequence diversity at the intraspecific level.ResultsWe examined patterns of interspecific divergence and intraspecific polymorphism in the plastid genome of Campanulastrum americanum, and across plastid genes found a significant correlation between dN/dS and pN/pS (i.e., the within-species equivalent of dN/dS). A number of genes including ycf1, ycf2, clpP, and ribosomal protein genes exhibited high dN/dS ratios. McDonald-Kreitman tests detected little evidence for positive selection acting on these genes, likely due to the presence of substantial intraspecific divergence.ConclusionsThese results suggest that mechanisms leading to increased nucleotide substitution rates in the plastid genome are continuing to act at the intraspecific level. Accelerated plastid genome evolution may increase the likelihood of intraspecific cytonuclear genetic incompatibilities, and thereby contribute to the early stages of the speciation process.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0268-y) contains supplementary material, which is available to authorized users.

The Largest Plastid Genome of Monocots: a Novel Genome Type Containing AT Residue Repeats in the Slipper Orchid Cypripedium japonicum

Plastid genome studies on Orchidaceae, one of the largest families of flowering plants, have been limited to the subfamily Epidendroideae. To increase understanding of the plastid genome evolution in Orchidaceae, we performed a complete plastid genome analysis of the slipper orchid Cypripedium japonicum positioned in the basal clade of the family. No significant gene deletions were observed in the plastid genome of the slipper orchid composed of 85 coding genes including members of the ndh gene family, which are mostly deleted or show pseudogenisation in orchids, although matK was identified as a pseudogene due to a frameshift mutation. Results also revealed that C. japonicum contains the largest plastid genome (174,417 bp) within monocots and the third largest one in Magnoliophyta. This is a new type of plastid genome extended due to abnormally frequent AT residue repeats within noncoding regions without inverted repeat (IR) expansion. In addition, we detected 25 plastid microsatellites and compared them among seven populations from Korea and Japan. These microsatellites may be applicable to population genetics and conservation biology studies on slipper orchids, which have become rare and are nearing extinction. AT-rich regions in the introns may also play a role in effective splicing of the genome after matK lost its function. As AT-rich regions are difficult to amplify and sequence using normal sequencing technologies, we propose a revised methodology for sequencing these regions.