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.     

The chloroplast genome of Anomochloa marantoidea (Anomochlooideae; Poaceae) comprises a mixture of grass-like and unique features

Features in the complete plastome of Anomochloa marantoidea (Poaceae) were investigated. This species is one of four of Anomochlooideae, the crown node of which diverged before those of any other grass subfamily. The plastome was sequenced from overlapping amplicons using previously designed primers. The plastome of A. marantoidea is 138412 bp long with a typical gene content for Poaceae. Five regions were examined in detail because of prior surveys that identified structural alterations among graminoid Poales. Anomochloa marantoidea was found to have an intron in rpoC1, unlike other Poaceae. The insertion region of rpoC2 is unusually short in A. marantoidea compared with those of other grasses, but with atypically long subrepeats. Both ycf1 and ycf2 are nonfunctional as is typical in grasses, but A. marantoidea has a uniquely long ψycf1. Finally, the rbcL-psaI spacer in A. marantoidea is atypically short with no evidence of the ψrpl23 locus found in all other Poaceae. Some of these features are of noteworthy dissimilarity between A. marantoidea and those crown grasses for which entire plastomes have been sequenced. Complete plastome sequences of other Anomochlooideae and outgroups will further advance our understanding of the evolutionary events in the plastome that accompanied graminoid diversification.     

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 .     

Implications of the Plastid Genome Sequence of Typha (Typhaceae, Poales) for Understanding Genome Evolution in Poaceae

Plastid genomes of the grasses (Poaceae) are unusual in their organization and rates of sequence evolution. There has been a recent surge in the availability of grass plastid genome sequences, but a comprehensive comparative analysis of genome evolution has not been performed that includes any related families in the Poales. We report on the plastid genome of Typha latifolia, the first non-grass Poales sequenced to date, and we present comparisons of genome organization and sequence evolution within Poales. Our results confirm that grass plastid genomes exhibit acceleration in both genomic rearrangements and nucleotide substitutions. Poaceae have multiple structural rearrangements, including three inversions, three genes losses (accD, ycf1, ycf2), intron losses in two genes (clpP, rpoC1), and expansion of the inverted repeat (IR) into both large and small single-copy regions. These rearrangements are restricted to the Poaceae, and IR expansion into the small single-copy region correlates with the phylogeny of the family. Comparisons of 73 protein-coding genes for 47 angiosperms including nine Poaceae genera confirm that the branch leading to Poaceae has significantly accelerated rates of change relative to other monocots and angiosperms. Furthermore, rates of sequence evolution within grasses are lower, indicating a deceleration during diversification of the family. Overall there is a strong correlation between accelerated rates of genomic rearrangements and nucleotide substitutions in Poaceae, a phenomenon that has been noted recently throughout angiosperms. The cause of the correlation is unknown, but faulty DNA repair has been suggested in other systems including bacterial and animal mitochondrial genomes.     

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 evolution of chloroplast genes and genomes in ferns

Most of the publicly available data on chloroplast (plastid) genes and genomes come from seed plants, with relatively little information from their sister group, the ferns. Here we describe several broad evolutionary patterns and processes in fern plastid genomes (plastomes), and we include some new plastome sequence data. We review what we know about the evolutionary history of plastome structure across the fern phylogeny and we compare plastome organization and patterns of evolution in ferns to those in seed plants. A large clade of ferns is characterized by a plastome that has been reorganized with respect to the ancestral gene order (a similar order that is ancestral in seed plants). We review the sequence of inversions that gave rise to this organization. We also explore global nucleotide substitution patterns in ferns versus those found in seed plants across plastid genes, and we review the high levels of RNA editing observed in fern plastomes.     

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.     

Exploring the plastid genome disparity of liverworts

Sequencing the plastid genomes of land plants provides crucial improvements to our understanding of the plastome evolution of land plants. Although the number of available complete plastid genome sequences has rapidly increased in the recent years, only a few sequences have been yet released for the three bryophyte lineages, namely hornworts, liverworts, and mosses. Here, we explore the disparity of the plastome structure of liverworts by increasing the number of sequenced liverwort plastomes from five to 18. The expanded sampling included representatives of all major lineages of liverworts including the genus Haplomitrium. The disparity of the liverwort genomes was compared with other 2386 land plant plastomes with emphasis on genome size and GC‐content. We found evidence for structural conservatism of the plastid genomes in liverworts and a trend towards reduced plastome sequence length in liverworts and derived mosses compared to other land plants, including hornworts and basal lineages of mosses. Furthermore, Aneura and Haplomitrium were distinct from other liverworts by an increased GC content, with the one found in Haplomitrium only second to the lycophyte Selaginella. The results suggest the hypothesis that liverworts and other land plants inherited and conserved the plastome structure of their most recent algal ancestors.     

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.     

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.     

Plastome comparison and evolution within the tribes of Plantaginaceae: Insights from an Asian gypsyweed

In spite of availability of several plastomes representing different tribes of Plantaginaceae, sparse attempts have been made to understand the plastome structure, evolution, and phylogenomics. In the present study, we have made an effort to understand the gene content and plastome evolution in the family Plantaginaceae using the newly generated plastome sequence of Veronica ovata subsp. kiusiana, a taxon native to SE Asia. In the first-ever attempt, plastomes of seven out of 10 tribes of Plantaginaceae have been compared to understand the evolution across the tribes of Plantaginaceae. The size of the plastome of V. ovata subsp. kiusiana is 152,249 bp, showing a typical quadripartite structure containing LSC, SSC, and two IRs with the sizes of 83,187, 17,704, and 25,679 respectively. The plastome comparison revealed the unique deletions in ycf2 and ndhF genes of members of different tribes, and also revealed high nucleotide variable hotspots. The study also revealed six highly variable genes and intergenic spacer viz. rps16, rps15-ycf1, ccsA-ndhD, ndhC-trnV, petN-psbM, and ycf1-trnN as potential DNA barcodes for the genus Veronica. The phylogenomic study revealed the sister relationship between V. ovata subsp. kiusiana and V. persica and also suggested the tentative placement of seven tribes in the family Plantaginaceae.     

Total duplication of the small single copy region in the angiosperm plastome: Rearrangement and inverted repeat instability in Asarum

Premise of the Study

As more plastomes are assembled, it is evident that rearrangements, losses, intergenic spacer expansion and contraction, and syntenic breaks within otherwise functioning plastids are more common than was thought previously, and such changes have developed independently in disparate lineages. However, to date, the magnoliids remain characterized by their highly conserved plastid genomes (plastomes).

Methods

Illumina HiSeq and MiSeq platforms were used to sequence the plastomes of Saruma henryi and those of representative species from each of the six taxonomic sections of Asarum. Sequenced plastomes were compared in a phylogenetic context provided by maximum likelihood and parsimony inferences made using an additional 18 publicly available plastomes from early‐diverging angiosperm lineages.

Key Results

In contrast to previously published magnoliid plastomes and the newly sequenced Saruma henryi plastome published here, Asarum plastomes have undergone extensive disruption and contain extremely lengthy AT‐repeat regions. The entirety of the small single copy region (SSC) of A. canadense and A. sieboldii var. sieboldii has been incorporated into the inverted repeat regions (IR), and the SSC of A. delavayi is only 14 bp long. All sampled Asarum plastomes share an inversion of a large portion of the large single copy region (LSC) such that trnE‐UUC is adjacent to the LSC‐IR boundary.

Conclusions

Plastome divergence in Asarum appears to be consistent with trends seen in highly rearranged plastomes of the monocots and eudicots. We propose that plastome instability in Asarum is due to repetitive motifs that serve as recombinatory substrates and reduce genome stability.     

Characterization of chloroplast DNA microsatellites from Saccharum spp and related species

Microsatellites, or simple sequence repeats (SSRs), and their flanking regions in chloroplast genomes (plastomes) of some species of the family Poaceae were analyzed in silico to look for DNA sequence variations. Comparison of the complete chloroplast DNA sequences (cpDNAs) of sugarcane (Saccharum hybrid cv. SP-80-3280 and S. officinarum cv. NCo310) and related species, Agrostis stolonifera, Brachypodium distachyon, Hordeum vulgare subsp vulgare, Lolium perenne, Oryza nivara, O. sativa subsp indica, O. sativa subsp japonica, Sorghum bicolor, Triticum aestivum, Zea mays, and Z. mays cv. B73, allowed us to examine the organization of chloroplast SSRs (cpSSRs) in genic and intergenic regions. We identified 204 cpSSRs in the sugarcane cpDNA; 22.5% were in genic regions. The ndh, rps, trn, and rpl gene clusters of the chloroplasts had the most repeats. Mononucleotide repeats were the most abundant cpSSRs in these species; however, di-, tri-, tetra-, penta-, and hexanucleotide repeats were also identified. Many base substitutions and deletions/insertions were identified in the cpSSR loci and their flanking regions. Multiple alignments of all cpSSR sequences of Poaceae species made identification of nucleotide variability possible; repeat motifs are not uniformly distributed across the Poaceae plastomes, but are mostly confined to intergenic regions. Phylogeny was determined by maximum parsimony and neighbor-joining inference methods. The cpSSRs of these species were found to be polymorphic. It appears that individual cpSSRs in the Poaceae are stable, at least over short periods of evolutionary time. We conclude that the plastome database can be exploited for phylogenetic analysis and biotechnological development.     

Evolution of the rpoB-psbZ region in fern plastid genomes: notable structural rearrangements and highly variable intergenic spacers

Background  

The rpoB-psbZ (BZ) region of some fern plastid genomes (plastomes) has been noted to go through considerable genomic changes. Unraveling its evolutionary dynamics across all fern lineages will lead to clarify the fundamental process shaping fern plastome structure and organization.     

The evolution of chloroplast genome structure in ferns

The plastid genome (plastome) is a rich source of phylogenetic and other comparative data in plants. Most land plants possess a plastome of similar structure. However, in a major group of plants, the ferns, a unique plastome structure has evolved. The gene order in ferns has been explained by a series of genomic inversions relative to the plastome organization of seed plants. Here, we examine for the first time the structure of the plastome across fern phylogeny. We used a PCR-based strategy to map and partially sequence plastomes. We found that a pair of partially overlapping inversions in the region of the inverted repeat occurred in the common ancestor of most ferns. However, the ancestral (seed plant) structure is still found in early diverging branches leading to the osmundoid and filmy fern lineages. We found that a second pair of overlapping inversions occurred on a branch leading to the core leptosporangiates. We also found that the unique placement of the gene matK in ferns (lacking a flanking intron) is not a result of a large-scale inversion, as previously thought. This is because the intron loss maps to an earlier point on the phylogeny than the nearby inversion. We speculate on why inversions may occur in pairs and what this may mean for the dynamics of plastome evolution.     

The complete chloroplast genome of tall fescue (Lolium arundinaceum; Poaceae) and comparison of whole plastomes from the family Poaceae

In this paper, we describe the complete chloroplast genome of Lolium arundinaceum. This sequence is the culmination of a long-term project completed by >400 undergraduates who took general genetics at Middle Tennessee State University from 2004-2007. It was undertaken in an attempt to introduce these students to an open-ended experiential/exploratory lesson to produce and analyze novel data. The data they produced should provide the necessary information for both phylogenetic comparisons and plastome engineering of tall fescue. The fescue plastome (GenBank FJ466687) is 136048 bp with a typical quadripartite structure and a gene order similar to other grasses; 56% of the plastome is coding region comprised of 75 protein-coding genes, 29 tRNAs, four rRNAs, and one hypothetical coding region (ycf). Comparisons of Poaceae plastomes reveal size differences between the PACC (subfamilies Panicoideae, Arundinoideae, Centothecoideae, and Chloridoideae) and BOP (subfamilies Bambusoideae, Oryzoideae, and Pooideae) clades. Alignment analysis suggests that several potentially conserved large deletions in previously identified intergenic length polymorphic regions are responsible for the majority of the size discrepancy. Phylogenetic analysis using whole plastome data suggests that fescue closely aligns with Lolium perenne. Some unique features as well as phylogenetic branch length calculations, however, suggest that a number of changes have occurred since these species diverged.     

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.     

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.     

Phylogenomics,molecular evolution,and estimated ages of lineages from the deep phylogeny of Poaceae

The deeply diverging subfamilies of grasses: Anomochlooideae, Pharoideae, and Puelioideae, today inhabit tropical forest floors as sparsely distributed depauperate lineages. The BEP/PACMAD grasses, which make up the majority of the family, are the result of a more recent radiation. Species in the deeply diverging subfamilies were here investigated to better understand molecular evolutionary processes and ages of divergence. Complete chloroplast genomes (plastomes) of Pharus latifolius L., P. lappulaceus Aubl., and Puelia olyriformis (Franch.) Clayton were determined. Four plastome loci from seven species of the deep subfamilies were also sequenced. Phylogenetic and mutation analyses and divergence estimations were conducted on all sequences together with homologous sequences from other Poaceae. Mutation analyses surveyed insertion/deletion mutations across the plastomes, clarified a trend in the molecular evolution of the rpoC2 locus, and indicated unique pseudogenizations in the plastomes of Pharus and Puelia. Phylogenetic analyses largely confirmed earlier multi-gene phylogenies. Phylogenomic and divergence analyses produced estimated origins of the crown nodes of Anomochlooideae at 65–104 Ma, Pharoideae at 44–71 Ma, and Puelioideae at 62–96 Ma. The upper ends of our estimated ranges are in general agreement with previous estimates. However, the lower ends of our ranges are considerably older than previous estimates, reflecting the influence of the less commonly used oldest fossil calibration point. The deeply diverging subfamilies exhibited the accumulation of numerous substitution and indel mutations consistent with a long evolutionary history that predated the radiation of the BEP/PACMAD grasses. We hypothesize that relatively rapid warming and drying in Africa at 55–56.5 Ma may have acted as selective forces stimulating adaptive radiations of grasses from the African tropical forests into diverse habitats.     

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.