Electron microscopic localization of replication origins in Oenothera chloroplast DNA

Summary The origins of chloroplast DNA (cpDNA) replication were mapped in two plastome types of Oenothera in order to determine whether variation in the origin of cpDNA replication could account for the different transmission abilities associated with these plastomes. Two pairs of displacement loop (D-loop) initiation sites were observed on closed circular cpDNA molecules by electron microscopy. Each pair of D-loops was mapped to the inverted repeats of the Oenothera cpDNA by the analysis of restriction fragments. The starting points of the two adjacent D-loops are approximately 4 kb apart, bracketing the 16S rRNA gene. Although there are small DNA length variations near one of the D-loop initiation sites, no apparent differences in the number and the location of replication origins were observed between plastomes with the highest (type I) and lowest (type IV) transmission efficiencies.     

Sequence variation in the putative replication origins of the five genetically distinct basicEuoenothera plastid chromosomes (plastomes)

Approximately 4200 nucleotides of the 16S/23S rDNA spacer and the 5 region flanking therrn operon of the plastid chromosomes representing the five basic, phylogenetically relatedEuoenothera plastomes were sequenced and compared. The sequences that harbor the putative replication origins are almost identical except for a 785-bp intercistronic segment between the genes for the 16S rRNA andtrnI. Differences are mainly caused by insertions/deletions and duplications; the predicted potential for formation of quite extensive secondary structure differs among the plastomes. Unexpected intraplastome variation has also been noted. Furthermore, the sequence-based and published genetically deduced plastome pedigrees differ significantly.     

Sequence variation in the putative replication origins of the five genetically distinct basicEuoenothera plastid chromosomes (plastomes)

Approximately 4200 nucleotides of the 16S/23S rDNA spacer and the 5′ region flanking therrn operon of the plastid chromosomes representing the five basic, phylogenetically relatedEuoenothera plastomes were sequenced and compared. The sequences that harbor the putative replication origins are almost identical except for a 785-bp intercistronic segment between the genes for the 16S rRNA andtrnI. Differences are mainly caused by insertions/deletions and duplications; the predicted potential for formation of quite extensive secondary structure differs among the plastomes. Unexpected intraplastome variation has also been noted. Furthermore, the sequence-based and published genetically deduced plastome pedigrees differ significantly.     

Functional analysis of plastid DNA replication origins in tobacco by targeted inactivation

Sequences described as chloroplast DNA replication origins were analysed in vivo by creating deletion and insertion mutants via plastid transformation in tobacco. Deletion of the described oriA sequence, which is located within the intron of the trnI gene, resulted in heteroplastomic transformants, when the selection marker was inserted within the intron. Removal of the complete intron sequence together with the oriA sequence, however, yielded homoplastomic transformants of normal phenotype, in which wild-type signals were no longer detectable through Southern analysis, thus bringing the role of the described oriA sequence for plastome replication into question. Similarly, deletion of sequence elements upstream of trnI, which have a possible ori function in Oenothera, did not show any effect in tobacco. The two copies of oriB, which are located at the very end of the plastome Inverted Repeats, were targeted with two different transformation vectors in a cotransformation approach. While in initial transformants integration of the selection marker could be detected at both sites, the transgene was found exclusively at one site or the other after additional rounds of regeneration. Whereas the copy of oriB in Inverted Repeat B could be completely deleted, targeting of the copy in Inverted Repeat A resulted in heteroplastomic lines, as the essential ycf1 gene was also affected. Due to the strong selection against cotransformants we conclude that at least one copy of the oriB sequence is essential for plastome replication, whereas replication appears possible without oriA elements.     

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.     

Study of direct repeats in micro evolution of plant mitochondria and plastids based on protein clustering

The study focuses on insertions of perfect direct repeats of words of arbitrary length in plastomes and mitochondriomes. The approach is exemplified using seed plants. Plastomes of close species were analyzed to further develop and refine published evidence for the evolution of non-coding DNA. The results suggest that perfect repeats are common elementary events resulting from replication errors—duplication of DNA. The role of such duplications in the evolution of the plastome is discussed.     

Plastome characteristics of Cannabaceae

Cannabaceae is an economically important family that includes ten genera and ca.117 accepted species. To explore the structure and size variation of their plastomes,we sequenced ten plastomes representing all ten genera of Cannabaceae.Each plastome possessed the typical angiosperm quadripartite structure and contained a total of 128 genes.The Inverted Repeat (IR) regions in five plastomes had experienced small expansions (330-983 bp) into the Large Single-Copy (LSC) region.The plastome of Chaetachme aristata has experienced a 942-bp IR contraction and lost rpl22 and rps19 in its IRs.The substitution rates of rps19 and rpl22 decreased after they shifted from the LSC to IR.A 270-bp inversion was detected in the Parasponia rugosa plastome,which might have been mediated by 18-bp inverted repeats.Repeat sequences,simple sequence repeats,and nucleotide substitution rates varied among these plastomes. Molecular markers with more than 13% variable sites and 5% parsimony-informative sites were identified,which may be useful for further phylogenetic analysis and species identification.Our results show strong support for a sister relationship between Gironniera and Lozanell (BS=100).Celtis,Cannabis-Humulus,Chaetachme-Pteroceltis,and Trema-Parasponia formed a strongly supported clade,and their relationships were well resolved with strong support (BS=100).The availability of these ten plastomes provides valuable genetic information for accurately identifying species,clarifying taxonomy and reconstructing the intergeneric phylogeny of Cannabaceae.     

Mutational meltdown or controlled chain reaction: The dynamics of rapid plastome evolution in the hyperdiversity of Poaceae

The study of genomic structural evolution associated with accelerated evolutionary rates that result in avoidance of meltdown and increase biodiversity is becoming ever more possible as the number of available plastomes increases. To more comprehensively analyze rate heterogeneity among monocots and within Poaceae, we sequenced plastomes from four Poaceae species, combined them with publicly available data from ~200 plastomes, and conducted comparative analyses to quantify the pattern of rate heterogeneity between different lineages, functional groups, and periods of evolutionary time. We compared structural differences across the Poaceae to quantify how changes in plastome size correspond to different genomic subunits and the evolution of IR–SC junction boundaries. The substitution rates among ancestral Poaceae were inferred to be exceptionally rapid compared to other monocots but slowed after divergence into extant lineages, which could not be sufficiently explained by positive selection. As such, rapid rates in the ancestral lineage leading to Poaceae might be more closely linked to large-scale structural changes like the loss of ycf1 and ycf2. The total increase in plastome size across Poaceae was positively correlated with the total length of intergenic spacers, tandem repeats, and dispersed repeats as well as large single copy, and inverted repeats (IRs). The continuous evolution of IR–SC junction boundaries was asynchronous with sizes of total genome and subunits across Poaceae. Future work is needed to better understand what factors in ancestral Poaceae evolved to harness such rapid rates of plastome evolution, avoid a mutational meltdown, and escape the stagnation of strong purifying selection as well as if these factors could be utilized to synthetically control rates.     

Passiflora plastome sequencing reveals widespread genomic rearrangements

Although past studies have included Passiflora among angiosperm lineages with highly rearranged plastid genomes (plastomes), knowledge about plastome organization in the genus is limited. So far only one draft and one complete plastome have been published. Expanded sampling of Passiflora plastomes is needed to understand the extent of the genomic rearrangement in the genus, which is also unusual in having biparental plastid inheritance and plastome‐genome incompatibility. We sequenced 15 Passiflora plastomes using either Illumina paired‐end or shotgun cloning and Sanger sequencing approaches. Assembled plastomes were annotated using Dual Organellar GenoMe Annotator (DOGMA) and tRNAscan‐SE. The Populus trichocarpa plastome was used as a reference to estimate genomic rearrangements in Passiflora by performing whole genome alignment in progressiveMauve. The phylogenetic distribution of rearrangements was plotted on the maximum likelihood tree generated from 64 plastid encoded protein genes. Inverted repeat (IR) expansion/contraction and loss of the two largest hypothetical open reading frames, ycf1 and ycf2, account for most plastome size variation, which ranges from 139 262 base pairs (bp) in P. biflora to 161 494 bp in P. pittieri. Passiflora plastomes have experienced numerous inversions, gene and intron losses along with multiple independent IR expansions and contractions resulting in a distinct organization in each of the three subgenera examined. Each Passiflora subgenus has a unique plastome structure in terms of gene content, order and size. The phylogenetic distribution of rearrangements shows that Passiflora has experienced widespread genomic changes, suggesting that such events may not be reliable phylogenetic markers.     

The first complete plastome sequence of the basal asterid family Styracaceae (Ericales) reveals a large inversion

Plastome sequences are rich sources of information for resolving difficult phylogenetic relationships and provide genomic data for conservation studies. Here, the complete plastome sequence of Alniphyllum eberhardtii Guillaumin is reported, representing the first plastome of the basal asterid family Styracaceae (Ericales). The plastome is 155,384 bp in length and contains 79 protein-coding genes, 30 tRNA genes and 4 rRNA genes, totaling 113 unique genes with 19 genes in the inverted repeat region. Unusual features of the plastome include the presence a large 20-kb inversion in the Large Single-Copy region, the pseudogenization of the accD gene, and the loss of the second intron from clpP. The 20-kb inversion includes 14 genes and has not been previously reported in other Ericales plastomes. Thirty-nine plastid simple sequence repeats (SSRs) that may provide genetic resources for the conservation of this economically import timber plant are characterized. Phylogenetic results inferred from ML and MP analyses of 66 plastid genes and 26 taxa reveal that the Styracaceae are sister to a clade including Actinidiaceae and Ericaceae and suggest that complete plastomes are likely to be very helpful in resolving the basal relationships among Ericales families, which have resisted resolution in smaller phylogenetic data sets.     

Plastomes of Mimosoideae: structural and size variation,sequence divergence,and phylogenetic implication

Plastomes of Fabaceae show both significant structural and size variation; however, most published plastomes are from subfamily Papilionoideae and only a few are from the other two subfamilies. In order to address the plastome structural and size variation of subfamily Mimosoideae, we integrated 11 newly sequenced plastomes from representing genera with three previously published ones. Each mimosoid plastome presented a typical quadripartite structure and contained 111 unique genes. Their inverted repeats (IRs) experienced multiple expansion/contraction; a ca. 13-kb IR expansion into small single copy (SSC) was detected in plastomes of a clade formed by tribe Ingeae and Acacia sensu stricto (s.s.), and a ca. 1.7-kb IR expansion into and a ca. 1.9-kb contraction out of large single copy (LSC) were found in Pithecellobium flexicaule and Acacia dealbata, respectively. Linear regression analysis showed decreased synonymous substitution rates of genes relocating from SSC into IR. A loss of both introns of clpP occurred in A. dealbata and Faidherbia albida, and a duplicated clpP copy was detected in A. dealbata. Furthermore, a 421-bp inversion that containing rps18 was found in A. dealbata. The size of mimosoid plastomes was found significantly affected by a IR-SC boundary shift, and also associated with repeat content. Plastome coding and noncoding regions with variable sequence divergence may supply valuable markers for molecular evolutionary and phylogenetic studies at different taxonomic levels. Plastid phylogenomics well resolved relationships among sampled mimosoid species.     

Plastome mutator-induced alterations arise in Oenothera chloroplast DNA through template slippage.

The plastome mutator of Oenothera hookeri strain Johansen causes deletions and duplications at target sites defined by direct repeats in the plastid genome. Previous studies characterized the mutations long after they had occurred and could not discriminate between the possibilities that the plastome mutator acted through unequal homologous recombination or template slippage. From the known hotspots, the rRNA spacer in the large inverted repeat was chosen for this study because it contains both direct and indirect repeats. Identical deletions were recovered from independently derived plants; the altered regions were always flanked by direct repeats. The regions in which the deletions occurred have the potential to form secondary structures that would stabilize the intervening sequence. Of the two affected regions, the one with the stronger potential secondary structure was altered more frequently. Because no duplication products or inversions were recovered, it is proposed that the plastome mutator acts through template slippage rather than through a recombination mechanism.     

Plastid Genome Assembly Using Long-read data

Although plastid genome (plastome) structure is highly conserved across most seed plants, investigations during the past two decades have revealed several disparately related lineages that experienced substantial rearrangements. Most plastomes contain a large inverted repeat and two single-copy regions, and a few dispersed repeats; however, the plastomes of some taxa harbour long repeat sequences (>300 bp). These long repeats make it challenging to assemble complete plastomes using short-read data, leading to misassemblies and consensus sequences with spurious rearrangements. Single-molecule, long-read sequencing has the potential to overcome these challenges, yet there is no consensus on the most effective method for accurately assembling plastomes using long-read data. We generated a pipeline, plastid Genome Assembly Using Long-read data (ptGAUL), to address the problem of plastome assembly using long-read data from Oxford Nanopore Technologies (ONT) or Pacific Biosciences platforms. We demonstrated the efficacy of the ptGAUL pipeline using 16 published long-read data sets. We showed that ptGAUL quickly produces accurate and unbiased assemblies using only ~50× coverage of plastome data. Additionally, we deployed ptGAUL to assemble four new Juncus (Juncaceae) plastomes using ONT long reads. Our results revealed many long repeats and rearrangements in Juncus plastomes compared with basal lineages of Poales. The ptGAUL pipeline is available on GitHub: https://github.com/Bean061/ptgaul .     

Plastid genome of <Emphasis Type="Italic">Seseli montanum</Emphasis>: Complete sequence and comparison with plastomes of other members of the Apiaceae family

This work reports the complete plastid (pt) DNA sequence of Seseli montanum L. of the Apiaceae family, determined using next-generation sequencing technology. The complete genome sequence has been deposited in GenBank with accession No. KM035851. The S. montanum plastome is 147,823 bp in length. The plastid genome has a typical structure for angiosperms and contains a large single-copy region (LSC) of 92,620 bp and a small single-copy region (SSC) of 17,481 bp separated by a pair of 18,861 bp inverted repeats (IRa and IRb). The composition, gene order, and AT-content in the S. montanum plastome are similar to that of a typical flowering plant pt DNA. One hundred fourteen unique genes have been identified, including 30 tRNA genes, four rRNA genes, and 80 protein genes. Of 18 intron-containing genes found, 16 genes have one intron, and two genes (ycf3, clpP) have two introns. Comparative analysis of Apiaceae plastomes reveals in the S. montanum plastome a LSC/IRb junction shift, so that the part of the ycf2 (4980 bp) gene is located in the LSC, but the other part of ycf2 (1301 bp) is within the inverted repeat. Thus, structural rearrangements in the plastid genome of S. montanum result in an enlargement of the LSC region by means of capture of a large part of ycf2, in contrast to eight Apiaceae plastomes where the complete ycf2 gene sequence is located in the inverted repeat.     

Structural features of a wheat plastome as revealed by complete sequencing of chloroplast DNA

Structural features of the wheat plastome were clarified by comparison of the complete sequence of wheat chloroplast DNA with those of rice and maize chloroplast genomes. The wheat plastome consists of a 134,545-bp circular molecule with 20,703-bp inverted repeats and the same gene content as the rice and maize plastomes. However, some structural divergence was found even in the coding regions of genes. These alterations are due to illegitimate recombination between two short direct repeats and/or replication slippage. Overall comparison of chloroplast DNAs among the three cereals indicated the presence of some hot-spot regions for length mutations. Whereas the region with clustered tRNA genes and that downstream of rbcL showed divergence in a species-specific manner, the deletion patterns of ORFs in the inverted-repeat regions and the borders between the inverted repeats and the small single-copy region support the notion that wheat and rice are related more closely to each other than to maize.     

Complete sequencing and comparative analyses of the pepper (Capsicum annuum L.) plastome revealed high frequency of tandem repeats and large insertion/deletions on pepper plastome

Plants in the family Solanaceae are used as model systems in comparative and evolutionary genomics. The complete chloroplast genomes of seven solanaceous species have been sequenced, including tobacco, potato and tomato, but not peppers. We analyzed the complete chloroplast genome sequence of the hot pepper, Capsicum annuum. The pepper chloroplast genome was 156,781 bp in length, including a pair of inverted repeats (IR) of 25,783 bp. The content and the order of 133 genes in the pepper chloroplast genome were identical to those of other solanaceous plastomes. To characterize pepper plastome sequence, we performed comparative analysis using complete plastome sequences of pepper and seven solanaceous plastomes. Frequency and contents of large indels and tandem repeat sequences and distribution pattern of genome-wide sequence variations were investigated. In addition, a phylogenetic analysis using concatenated alignments of coding sequences was performed to determine evolutionary position of pepper in Solanaceae. Our results revealed two distinct features of pepper plastome compared to other solanaceous plastomes. Firstly, large indels, including insertions on accD and rpl20 gene sequences, were predominantly detected in the pepper plastome compared to other solanaceous plastomes. Secondly, tandem repeat sequences were particularly frequent in the pepper plastome. Taken together, our study represents unique features of evolution of pepper plastome among solanaceous plastomes.     

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine proteincoding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.     

Do nonasterid holoparasitic flowering plants have plastid genomes?

Past work involving the plastid genome (plastome) of holoparasitic plants has been confined to Scrophulariaceae (or Orobanchaceae) which have truncated plastomes owing to loss of photosynthetic and other genes. Nonasterid holoparasites from Balanophoraceae (Corynaea), Hydnoraceae (Hydnora) and Cytinaceae (Cytinus) were tested for the presence of plastid genes and a plastome. Using PCR, plastid 16S rDNA was successfully amplified and sequenced from the above three holoparasites. The sequence of Cytinus showed 121 single base substitutions relative to Nicotiana (8% of the molecule) whereas higher sequence divergence was observed in Hydnora and Corynaea (287 and 513 changes, respectively). Secondary structural models for these 16S rRNAs show that most changes are compensatory, thus suggesting they are functional. Probes constructed for 16S rDNA and for four plastid-encoded ribosomal protein genes (rps2, rps4, rps7 and rpl16) were used in Southern blots of digested genomic DNA from the three holoparasites. Positive hybridizations were obtained using each of the five probes only for Cytinus. For SmaI digests, all plastid gene probes hybridized to a common fragment ca. 20 kb in length in this species. Taken together, these data provide preliminary evidence suggestive of the retention of highly diverged and truncated plastid genome in Cytinus. The greater sequence divergence for 16S rDNA and the negative hybridization results for Hydnora and Corynaea suggests two possibilities: the loss of typically conserved elements of their plastomes or the complete absence of a plastome.     

Plastome organization and evolution of chloroplast genes in Cardamine species adapted to contrasting habitats

Background

Plastid genomes, also known as plastomes, are shaped by the selective forces acting on the fundamental cellular functions they code for and thus they are expected to preserve signatures of the adaptive path undertaken by different plant species during evolution. To identify molecular signatures of positive selection associated to adaptation to contrasting ecological niches, we sequenced with Solexa technology the plastomes of two congeneric Brassicaceae species with different habitat preference, Cardamine resedifolia and Cardamine impatiens.

Results

Following in-depth characterization of plastome organization, repeat patterns and gene space, the comparison of the newly sequenced plastomes between each other and with 15 fully sequenced Brassicaceae plastomes publically available in GenBank uncovered dynamic variation of the IR boundaries in the Cardamine lineage. We further detected signatures of positive selection in ten of the 75 protein-coding genes of the examined plastomes, identifying a range of chloroplast functions putatively involved in adaptive processes within the family. For instance, the three residues found to be under positive selection in RUBISCO could possibly be involved in the modulation of RUBISCO aggregation/activation and enzymatic specificty in Brassicaceae. In addition, our results points to differential evolutionary rates in Cardamine plastomes.

Conclusions

Overall our results support the existence of wider signatures of positive selection in the plastome of C. resedifolia, possibly as a consequence of adaptation to high altitude environments. We further provide a first characterization of the selective patterns shaping the Brassicaceae plastomes, which could help elucidate the driving forces underlying adaptation and evolution in this important plant family.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1498-0) contains supplementary material, which is available to authorized users.     

The Complete Chloroplast Genome of <Emphasis Type="Italic">Coix lacryma-jobi</Emphasis> and a Comparative Molecular Evolutionary Analysis of Plastomes in Cereals

Graminoid molecular evolution was investigated by chloroplast genome (plastome) scale analyses. A complete plastome from Coix lacryma-jobi (Poaceae) and a draft plastome from Joinvillea plicata (Joinvilleaceae) were sequenced and analyzed. The draft plastome included conserved protein-coding loci routinely analyzed in previous studies plus one additional locus of demonstrated phylogenetic utility. The methodological approach was to directly sequence overlapping amplicons from known plastome regions. Over 100 pairs of amplification and sequencing primers were designed and positioned to flank overlapping 1,200-base pair fragments around the entire plastome. Newly determined sequences were analyzed with published plastomes from representatives of Panicoideae, Ehrhartoideae, and Pooideae. Considerable variation was found for studies within the family and even within Andropogoneae. Readily interpreted mutation patterns were observed, such as small inversions in hairpin-loop regions and indels, which were common in intergenic spacers. Maximum or near-maximum bootstrap support was observed in all analyses resolving relationships between subfamilies. However, the addition of characters from noncoding regions increased the number of parsimony-informative characters and lengthened short internal branches (Andropogoneae), better defining intergeneric relationships. Thus, characters in complete plastomes can be used over a wide scope of phylogenetic studies.