Spontaneous Chloroplast Mutants Mostly Occur by Replication Slippage and Show a Biased Pattern in the Plastome of Oenothera

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Leaf variegation caused by plastome structural variation: an example from Dianella tasmanica

Variegated plants often exhibit plastomic heteroplasmy due to single-nucleotide mutations or small insertions/deletions in their albino sectors. Here, however, we identified a plastome structural variation in albino sectors of the variegated plant Dianella tasmanica (Asphodelaceae), a perennial herbaceous plant widely cultivated as an ornamental in tropical Asia. This structural variation, caused by intermolecular recombination mediated by an 11-bp inverted repeat flanking a 92-bp segment in the large single-copy region (LSC), generates a giant plastome (228 878 bp) with the largest inverted repeat of 105 226 bp and the smallest LSC of 92 bp known in land plants. It also generates an ~7-kb deletion on the boundary of the LSC, which eliminates three protein coding genes (psbA, matK, and rps16) and one tRNA gene (trnK). Albino sectors exhibit dramatic changes in expression of many plastid genes, including negligible expression of psbA, matK, and rps16, reduced expression of photosynthesis-related genes, and increased expression of genes related to the translational apparatus. Microscopic and ultrastructure observations showed that albino tissues were present in both green and albino sectors of the variegated individuals, and chloroplasts were poorly developed in the mesophyll cells of the albino tissues of the variegated individuals. These poorly developed chloroplasts likely carry the large and rearranged plastome, which is likely responsible for the loss of photosynthesis and albinism in the leaf margins. Considering that short repeats are relatively common in plant plastomes and that photosynthesis is not necessary for albino sectors, structural variation of this kind may not be rare in the plastomes of variegated plants.

Crops under past diversification and ongoing climate change: more than just selection of nuclear genes for flowering

Diversification and breeding following domestication and under current climate change across the globe are the two most significant evolutionary events experienced by major crops. Diversification of crops from their wild ancestors has favored dramatic changes in the sensitivity of the plants to the environment, particularly significantly in transducing light inputs to the circadian clock, which has allowed the growth of major crops in the relatively short growing season experienced in the Northern Hemisphere. Historically, mutants and the mapping of quantitative trait loci (QTL) have facilitated the identification and the cloning of genes that underlie major changes of the clock and the regulation of flowering. Recent studies have suggested that the thermal plasticity of the circadian clock output, and not just the core genes that follow temperature compensation, has also been under selection during diversification and breeding. Wild alleles that accelerate output rhythmicity could be beneficial for crop resilience. Furthermore, wild alleles with beneficial and flowering-independent effects under stress indicate their possible role in maintaining a balanced source-sink relationship, thereby allowing productivity under climatic change. Because the chloroplast genome also regulates the plasticity of the clock output, mapping populations including cytonuclear interactions should be utilized within an integrated field and clock phenomics framework. In this review, we highlight the need to integrate physiological and developmental approaches (physio-devo) to gain a better understanding when re-domesticating wild gene alleles into modern cultivars to increase their robustness under abiotic heat and drought stresses.

High-throughput discovery of plastid genes causing albino phenotypes in ornamental chimeric plants

Chimeric plants composed of green and albino tissues have great ornamental value. To unveil the functional genes responsible for albino phenotypes in chimeric plants, we inspected the complete plastid genomes (plastomes) in green and albino leaf tissues from 23 ornamental chimeric plants belonging to 20 species, including monocots, dicots, and gymnosperms. In nine chimeric plants, plastomes were identical between green and albino tissues. Meanwhile, another 14 chimeric plants were heteroplasmic, showing a mutation between green and albino tissues. We identified 14 different point mutations in eight functional plastid genes related to plastid-encoded RNA polymerase (rpo) or photosystems which caused albinism in the chimeric plants. Among them, 12 were deleterious mutations in the target genes, in which early termination appeared due to small deletion-mediated frameshift or single nucleotide substitution. Another was single nucleotide substitution in an intron of the ycf3 and the other was a missense mutation in coding region of the rpoC2 gene. We inspected chlorophyll structure, protein functional model of the rpoC2, and expression levels of the related genes in green and albino tissues of Reynoutria japonica. A single amino acid change, histidine-to-proline substitution, in the rpoC2 protein may destabilize the peripheral helix of plastid-encoded RNA polymerase, impairing the biosynthesis of the photosynthesis system in the albino tissue of R. japonica chimera plant.

Long-read sequencing characterizes mitochondrial and plastid genome variants in Arabidopsis msh1 mutants

The abundant repeats in plant mitochondrial genomes can cause rapid genome rearrangements and are also a major obstacle in short-read sequencing studies. Nuclear-encoded proteins such as MSH1 are known to suppress the generation of repeat-associated mitochondrial genome variants, but our understanding of these mechanisms has been constrained by the limitations of short-read technologies. Here, we used highly accurate long-read sequencing (PacBio HiFi) to characterize mitochondrial and plastid genome variants in Arabidopsis thaliana msh1 mutant individuals. The HiFi reads provided a global view of recombination dynamics with detailed quantification of parental and crossover recombination products for both large and small repeats. We found that recombination breakpoints were distributed relatively evenly across the length of repeated sequences and detected widespread internal exchanges of sequence variants between pairs of imperfect repeats in the mitochondrial genome of msh1 mutants. Long-read assemblies of mitochondrial genomes from seven other A. thaliana wild-type accessions differed by repeat-mediated structural rearrangements similar to those observed in msh1 mutants, but they were all in a simple low-heteroplasmy state. The Arabidopsis plastid genome generally lacks small repeats and exhibited a very different pattern of variant accumulation in msh1 mutants compared with the mitochondrial genome. Our data illustrate the power of HiFi technology in studying repeat-mediated recombination in plant organellar genomes and improved the sequence resolution for recombinational processes suppressed by MSH1. Plant organellar genomes can undergo rapid rearrangements. Long-read sequencing provides a detailed and quantitative view of mitochondrial and plastid genome variants normally suppressed by MSH1, advancing our understanding of plant organellar genome dynamics.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 (MSH1), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis. We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes (N) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria (N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

The evolution of multi-gene families and metabolic pathways in the evening primroses (Oenothera: Onagraceae): A comparative transcriptomics approach

The plant genus Oenothera has played an important role in the study of plant evolution of genomes and plant defense and reproduction. Here, we build on the 1kp transcriptomic dataset by creating 44 new transcriptomes and analyzing a total of 63 transcriptomes to present a large-scale comparative study across 29 Oenothera species. Our dataset included 30.4 million reads per individual and 2.3 million transcripts on average. We used this transcriptome resource to examine genome-wide evolutionary patterns and functional diversification by searching for orthologous genes and performed gene family evolution analysis. We found wide heterogeneity in gene family evolution across the genus, with section Oenothera exhibiting the most pronounced evolutionary changes. Overall, more significant gene family expansions occurred than contractions. We also analyzed the molecular evolution of phenolic metabolism by retrieving proteins annotated for phenolic enzymatic complexes. We identified 1,568 phenolic genes arranged into 83 multigene families that varied widely across the genus. All taxa experienced rapid phenolic evolution (fast rate of genomic turnover) involving 33 gene families, which exhibited large expansions, gaining about 2-fold more genes than they lost. Upstream enzymes phenylalanine ammonia-lyase (PAL) and 4-coumaroyl: CoA ligase (4CL) accounted for most of the significant expansions and contractions. Our results suggest that adaptive and neutral evolutionary processes have contributed to Oenothera diversification and rapid gene family evolution.

Graph-based models of the Oenothera mitochondrial genome capture the enormous complexity of higher plant mitochondrial DNA organization

Plant mitochondrial genomes display an enormous structural complexity, as recombining repeat-pairs lead to the generation of various sub-genomic molecules, rendering these genomes extremely challenging to assemble. We present a novel bioinformatic data-processing pipeline called SAGBAC (Semi-Automated Graph-Based Assembly Curator) that identifies recombinogenic repeat-pairs and reconstructs plant mitochondrial genomes. SAGBAC processes assembly outputs and applies our novel ISEIS (Iterative Sequence Ends Identity Search) algorithm to obtain a graph-based visualization. We applied this approach to three mitochondrial genomes of evening primrose (Oenothera), a plant genus used for cytoplasmic genetics studies. All identified repeat pairs were found to be flanked by two alternative and unique sequence-contigs defining so-called 'double forks', resulting in four possible contig-repeat-contig combinations for each repeat pair. Based on the inferred structural models, the stoichiometry of the different contig-repeat-contig combinations was analyzed using Illumina mate-pair and PacBio RSII data. This uncovered a remarkable structural diversity of the three closely related mitochondrial genomes, as well as substantial phylogenetic variation of the underlying repeats. Our model allows predicting all recombination events and, thus, all possible sub-genomes. In future work, the proposed methodology may prove useful for the investigation of the sub-genome organization and dynamics in different tissues and at various developmental stages.

Repairing TALEN-mediated double-strand break by microhomology-mediated recombination in tobacco plastids generates abundant subgenomic DNA

Transcription activator-like effector nuclease (TALEN) technology has been widely used to edit nuclear genomes in plants but rarely for editing organellar genomes. In addition, ciprofloxacin, commonly used to cause the double-strand break of organellar DNA for studying the repair mechanism in plants, confers no organellar selectivity and site-specificity. To demonstrate the feasibility of TALEN-mediated chloroplast DNA editing and to use it for studying the repair mechanism in plastids, we developed a TALEN-mediated editing technology fused with chloroplast transit peptide (cpTALEN) to site-specifically edit the rpoB gene via Agrobacteria-mediated transformation of tobacco leaf. Transgenic plants showed various degrees of chlorotic phenotype. Repairing damaged plastid DNA resulted in point mutation, large deletion and small inversion surrounding the rpoB gene by homologous recombination and/or microhomology-mediated recombination. In an albino line, microhomology-mediated recombination via a pair of 12-bp direct repeats between rpoC2 and ycf2 genes generated the chimeric ycf2-rpoC2 subgenome, with the level about 3- to 5-fold higher for subgenomic DNA than ycf2. Additionally, the expression of chimeric ycf2-rpoC2 transcripts versus ycf2 mRNA agreed well with the level of corresponding DNA. The ycf2-rpoC2 subgenomic DNA might independently and preferentially replicate in plastids.

Comparative Analysis of Complete Chloroplast Genomes of 13 Species in Epilobium, Circaea, and Chamaenerion and Insights Into Phylogenetic Relationships of Onagraceae

The evening primrose family, Onagraceae, is a well defined family of the order Myrtales, comprising 22 genera widely distributed from boreal to tropical areas. In this study, we report and characterize the complete chloroplast genome sequences of 13 species in Circaea, Chamaenerion, and Epilobium using a next-generation sequencing method. We also retrieved chloroplast sequences from two other Onagraceae genera to characterize the chloroplast genome of the family. The complete chloroplast genomes of Onagraceae encoded an identical set of 112 genes (with exclusion of duplication), including 78 protein-coding genes, 30 transfer RNAs, and four ribosomal RNAs. The chloroplast genomes are basically conserved in gene arrangement across the family. However, a large segment of inversion was detected in the large single copy region of all the samples of Oenothera subsect. Oenothera. Two kinds of inverted repeat (IR) region expansion were found in Oenothera, Chamaenerion, and Epilobium samples. We also compared chloroplast genomes across the Onagraceae samples in some features, including nucleotide content, codon usage, RNA editing sites, and simple sequence repeats (SSRs). Phylogeny was inferred by the chloroplast genome data using maximum-likelihood (ML) and Bayesian inference methods. The generic relationship of Onagraceae was well resolved by the complete chloroplast genome sequences, showing potential value in inferring phylogeny within the family. Phylogenetic relationship in Oenothera was better resolved than other densely sampled genera, such as Circaea and Epilobium. Chloroplast genomes of Oenothera subsect. Oenothera, which are biparental inheritated, share a syndrome of characteristics that deviate from primitive pattern of the family, including slightly expanded inverted repeat region, intron loss in clpP, and presence of the inversion.

A photosynthesis operon in the chloroplast genome drives speciation in evening primroses

Genetic incompatibility between the cytoplasm and the nucleus is thought to be a major factor in species formation, but mechanistic understanding of this process is poor. In evening primroses (Oenothera spp.), a model plant for organelle genetics and population biology, hybrid offspring regularly display chloroplast-nuclear incompatibility. This usually manifests in bleached plants, more rarely in hybrid sterility or embryonic lethality. Hence, most of these incompatibilities affect photosynthetic capability, a trait that is under selection in changing environments. Here we show that light-dependent misregulation of the plastid psbB operon, which encodes core subunits of photosystem II and the cytochrome b6f complex, can lead to hybrid incompatibility, and this ultimately drives speciation. This misregulation causes an impaired light acclimation response in incompatible plants. Moreover, as a result of their different chloroplast genotypes, the parental lines differ in photosynthesis performance upon exposure to different light conditions. Significantly, the incompatible chloroplast genome is naturally found in xeric habitats with high light intensities, whereas the compatible one is limited to mesic habitats. Consequently, our data raise the possibility that the hybridization barrier evolved as a result of adaptation to specific climatic conditions.

Correction of frameshift mutations in the atpB gene by translational recoding in chloroplasts of Oenothera and tobacco

Translational recoding, also known as ribosomal frameshifting, is a process that causes ribosome slippage along the messenger RNA, thereby changing the amino acid sequence of the synthesized protein. Whether the chloroplast employs recoding is unknown. I-iota, a plastome mutant of Oenothera (evening primrose), carries a single adenine insertion in an oligoA stretch [11A] of the atpB coding region (encoding the β-subunit of the ATP synthase). The mutation is expected to cause synthesis of a truncated, nonfunctional protein. We report that a full-length AtpB protein is detectable in I-iota leaves, suggesting operation of a recoding mechanism. To characterize the phenomenon, we generated transplastomic tobacco lines in which the atpB reading frame was altered by insertions or deletions in the oligoA motif. We observed that insertion of two adenines was more efficiently corrected than insertion of a single adenine, or deletion of one or two adenines. We further show that homopolymeric composition of the oligoA stretch is essential for recoding, as an additional replacement of AAA lysine codon by AAG resulted in an albino phenotype. Our work provides evidence for the operation of translational recoding in chloroplasts. Recoding enables correction of frameshift mutations and can restore photoautotrophic growth in the presence of a mutation that otherwise would be lethal.

Epimutations and mutations, nurturing phenotypic diversity

Epimutations and mutations are two dissimilar mechanisms that have contributed to the phenotypic diversities in organisms. Though dissimilar, many previous studies have revealed that the consequences of epimutations and mutations are not mutually exclusive. DNA rich in epigenetic modifications can be prone to mutations and vice versa. In order to get a better insight into the molecular evolution in organisms, it is important to consider the information of both genetic and epigenetic changes in their genomes. Understanding the similarities and differences between the consequences of epimutations and mutations is required for a better interpretation of phenotypic diversities in organisms. Factors contributing to epigenetic changes such as paramutations and mutation hotspots and, the correlation of the interdependence of mutations and epigenetic changes in DNA are important aspects that need to be considered for molecular evolutionary studies. Thus, this review explains what epimutations are, their causes, how they are similar/different from mutations, and the influence of epigenetic changes and mutations on each other, further emphasizing how molecular evolution involving both mutations and epimutations can lead to speciation. Considering this approach will aid in reorganizing taxonomic classifications, importantly, solving disparities in species identification.

Ultra-deep sequencing reveals dramatic alteration of organellar genomes in Physcomitrella patens due to biased asymmetric recombination

Destabilization of organelle genomes causes organelle dysfunction that appears as abnormal growth in plants and diseases in human. In plants, loss of the bacterial-type homologous recombination repair (HRR) factors RECA and RECG induces organelle genome instability. In this study, we show the landscape of organelle genome instability in Physcomitrella patens HRR knockout mutants by deep sequencing in combination with informatics approaches. Genome-wide maps of rearrangement positions in the organelle genomes, which exhibited prominent mutant-specific patterns, were highly biased in terms of direction and location and often associated with dramatic variation in read depth. The rearrangements were location-dependent and mostly derived from the asymmetric products of microhomology-mediated recombination. Our results provide an overall picture of organelle-specific gross genomic rearrangements in the HRR mutants, and suggest that chloroplasts and mitochondria share common mechanisms for replication-related rearrangements.

Masaki Odahara and Kensuke Nakamura et al. use deep paired-end sequencing to examine organellar genome recombination when homologous recombination repair genes are individually knocked out in the moss, Physcomitrella patens. Their results suggest that chloroplasts and mitochondria share a common mechanism for replication-related rearrangements.

Plastid genome evolution in Amazonian açaí palm (Euterpe oleracea Mart.) and Atlantic forest açaí palm (Euterpe edulis Mart.)

The plastomes of E. edulis and E. oleracea revealed several molecular markers useful for genetic studies in natural populations and indicate specific evolutionary features determined by vicariant speciation. Arecaceae is a large and diverse family occurring in tropical and subtropical ecosystems worldwide. E. oleracea is a hyperdominant species of the Amazon forest, while E. edulis is a keystone species of the Atlantic forest. It has reported that E. edulis arose from vicariant speciation after the emergence of the belt barrier of dry environment (Cerrado and Caatinga biomes) between Amazon and Atlantic forests, isolating the E. edulis in the Atlantic forest. We sequenced the complete plastomes of E. edulis and E. oleracea and compared them concerning plastome structure, SSRs, tandem repeats, SNPs, indels, hotspots of nucleotide polymorphism, codon Ka/Ks ratios and RNA editing sites aiming to investigate evolutionary traits possibly affected by distinct environments. Our analyses revealed 303 SNPs, 91 indels, and 82 polymorphic SSRs among both species. Curiously, the narrow correlation among localization of repetitive sequences and indels strongly suggests that replication slippage is involved in plastid DNA mutations in Euterpe. Moreover, most non-synonymous substitutions represent amino acid variants in E. edulis that evolved specifically or in a convergent manner across the palm phylogeny. Amino acid variants observed in several plastid proteins in E. edulis were also identified as positive signatures across palm phylogeny. The higher incidence of specific amino acid changes in plastid genes of E. edulis in comparison with E. oleracea probably configures adaptive genetic variations determined by vicariant speciation. Our data indicate that the environment generates a selective pressure on the plastome making it more adapted to specific conditions.

Complete chloroplast genome sequencing and comparative analysis reveals changes to the chloroplast genome after allopolyploidization in Cucumis

Allopolyploids undergo "genomic shock" leading to significant genetic and epigenetic modifications. Previous studies have mainly focused on nuclear changes, while little is known about the inheritance and changes of organelle genome in allopolyploidization. The synthetic allotetraploid Cucumis ×hytivus, which is generated via hybridization between C. hystrix and C. sativus, is a useful model system for studying cytonuclear variation. Here, we report the chloroplast genome of allotetraploid C. ×hytivus and its diploid parents via sequencing and comparative analysis. The size of the obtained chloroplast genomes ranged from 154 673 to 155 760 bp, while their gene contents, gene orders, and GC contents were similar to each other. Comparative genome analysis supports chloroplast maternal inheritance. However, we identified 51 indels and 292 SNP genetic variants in the chloroplast genome of the allopolyploid C. ×hytivus relative to its female parent C. hystrix. Nine intergenic regions with rich variation were identified through comparative analysis of the chloroplast genomes within the subgenus Cucumis. The phylogenetic network based on the chloroplast genome sequences clarified the evolution and taxonomic position of the synthetic allotetraploid C. ×hytivus. The results of this study provide us with an insight into the changes of organelle genome after allopolyploidization, and a new understanding of the cytonuclear evolution.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

Under the rug: Abandoning persistent misconceptions that obfuscate organelle evolution

The advent and advance of next generation sequencing over the past two decades made it possible to accumulate large quantities of sequence reads that could be used to assemble complete or nearly complete organelle genomes (plastome or mitogenome). The result has been an explosive increase in the availability of organelle genome sequences with over 4000 different species of green plants currently available on GenBank. During the same time period, plant molecular biologists greatly enhanced the understanding of the structure, repair, replication, recombination, transcription and translation, and inheritance of organelle DNA. Unfortunately many plant evolutionary biologists are unaware of or have overlooked this knowledge, resulting in misrepresentation of several phenomena that are critical for phylogenetic and evolutionary studies using organelle genomes. We believe that confronting these misconceptions about organelle genome organization, composition, and inheritance will improve our understanding of the evolutionary processes that underly organelle evolution. Here we discuss four misconceptions that can limit evolutionary biology studies and lead to inaccurate phylogenies and incorrect structure of the organellar DNA used to infer organelle evolution.

A point mutation in the photosystem I P700 chlorophyll a apoprotein A1 gene confers variegation in Helianthus annuus L

Even a point mutation in the psaA gene mediates chlorophyll deficiency. The role of the plastid signal may perform the redox state of the compounds on the acceptor-side of PSI. Two extranuclear variegated mutants of sunflower, Var1 and Var33, were investigated. The yellow sectors of both mutants were characterized by an extremely low chlorophyll and carotenoid content, as well as poorly developed, unstacked thylakoid membranes. A full-genome sequencing of the cpDNA revealed mutations in the psaA gene in both Var1 and Var33. The cpDNA from the yellow sectors of Var1 differs from those in the wild type by only a single, non-synonymous substitution (Gly734Glu) in the psaA gene, which encodes a subunit of photosystem (PS) I. In the cpDNA from the yellow sectors of Var33, the single-nucleotide insertion in the psaA gene was revealed, leading to frameshift at the 580 amino acid position. Analysis of the photosynthetic electron transport demonstrated an inhibition of the PSI and PSII activities in the yellow tissues of the mutant plants. It has been suggested that mutations in the psaA gene of both Var1 and Var33 led to the disruption of PSI. Due to the non-functional PSI, photosynthetic electron transport is blocked, which, in turn, leads to photodamage of PSII. These data are confirmed by immunoblotting analysis, which showed a significant reduction in PsbA in the yellow leaf sectors, but not PsaA. The expression of chloroplast and nuclear genes encoding the PSI subunits (psaA, psaB, and PSAN), the PSII subunits (psbA, psbB, and PSBW), the antenna proteins (LHCA1, LHCB1, and LHCB4), the ribulose 1.5-bisphosphate carboxylase subunits (rbcL and RbcS), and enzymes of chlorophyll biosynthesis were down-regulated in the yellow leaf tissue. The extremely reduced transcriptional activity of the two protochlorophyllide oxidoreductase (POR) genes involved in chlorophyll biosynthesis is noteworthy. The disruption of NADPH synthesis, due to the non-functional PSI, probably led to a significant reduction in NADPH-protochlorophyllide oxidoreductase in the yellow sectors of Var1 and Var33. A dramatic decrease in chlorophyllide was shown in the yellow sectors. A reduction in NADPH-protochlorophyllide oxidoreductase, along with photodegradation, has been suggested as a result of chlorophyll deficiency.

MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

Mitochondrial plastid DNA can cause DNA barcoding paradox in plants

The transfer of ancestral plastid genomes into mitochondrial genomes to generate mitochondrial plastid DNA (MTPT) is known to occur in plants, but its impacts on mitochondrial genome complexity and the potential for causing a false-positive DNA barcoding paradox have been underestimated. Here, we assembled the organelle genomes of Cynanchum wilfordii and C. auriculatum, which are indigenous medicinal herbs in Korea and China, respectively. In both species, it is estimated that 35% of the ancestral plastid genomes were transferred to mitochondrial genomes over the past 10 million years and remain conserved in these genomes. Some plastid barcoding markers co-amplified the conserved MTPTs and caused a barcoding paradox, resulting in mis-authentication of botanical ingredients and/or taxonomic mis-positioning. We identified dynamic and lineage-specific MTPTs that have contributed to mitochondrial genome complexity and might cause a putative barcoding paradox across 81 plant species. We suggest that a DNA barcoding guidelines should be developed involving the use of multiple markers to help regulate economically motivated adulteration.

Comparative analyses of chloroplast genomes from 13 Lagerstroemia (Lythraceae) species: identification of highly divergent regions and inference of phylogenetic relationships

Seven divergence hotspots as plastid markers for DNA barcoding was selected, and the phylogeny of 13 Lagerstroemia species based on the cp genome data was reconstructed within Myrtales. The Lagerstroemia species used in this study originated in China and have high economic and ecological value. The shared interspecific morphological characteristics and intraspecific morphological variation resulting from hybridization among Lagerstroemia taxa have made resolving their classification problems and phylogenetic relationships difficult. Systematic comparative genomic analysis has been shown to resolve phylogenetic relationships. We sequenced and annotated 6 Lagerstroemia cp genomes (Lagerstroemia excelsa, Lagerstroemia limii, Lagerstroemia siamica, Lagerstroemia tomentosa, Lagerstroemia venusta, and Lagerstroemia calyculata) for the first time and combined them with previously published genomes for Lagerstroemia species. Bioinformatics was used to analyse the 13 cp genomes in terms of gene structure and organization, codon usage, contraction and expansion of inverted repeat regions, repeat structure, divergence hotspots, species pairwise Ka/Ks ratios and phylogenetic relationships. The length varied between 152,049 bp in Lagerstroemia subcostata and 152,521 bp in L. venusta. We selected seven divergence hotspots in the cp genomes that had the potential to act as plastid markers to distinguish Lagerstroemia species. The phylogenetic relationships within Myrtales inferred from the cp genomes of 13 Lagerstroemia species and 27 other Myrtales species were highly supported, which illustrated several novel relationships within Myrtales. Taken together, our results provide comprehensive chloroplast genomic resources, which can be used further for species identification and molecular breeding of Lagerstroemia species.

Recurrent gene duplication in the angiosperm tribe Delphinieae (Ranunculaceae) inferred from intracellular gene transfer events and heteroplasmic mutations in the plastid matK gene

The study of intracellular gene transfer may allow for the detection of interesting evolutionary processes such as ancient polyploidization. We compared 24 plastid genomes (plastomes) from tribe Delphinieae, one from tribe Nigelleae and one from tribe Ranunculeae, including five newly sequenced genomes. The functional transfers of the plastids rpl32 and rps16 to the nucleus in tribe Delphinieae were identified. Unexpectedly, we discovered multiple divergent copies of the nuclear-encoded plastid rpl32 in the genus Aconitum. Phylogenetic and synonymous substitution rate analyses revealed that the nuclear-encoded plastid rpl32 underwent two major duplication events. These ancient gene duplication events probably occurred via multiple polyploidization events in Aconitum between 11.9 and 24.7 Mya. Furthermore, our sequence rate analysis indicated that the eight plastid-encoded rpl subunits in Aconitum had a significantly accelerated evolutionary rate compared to those in other genera, suggesting that highly divergent paralogs targeted to the plastid may contribute to an elevated rate of evolution in plastid rpl genes. In addition, heteroplasmy of the plastid matK from two Aconitum species suggested the existence of potentially functional plastid maturases in its plastome. Our results provide insight into the evolutionary history of the tribe Delphinieae.

DNA Repair and the Stability of the Plant Mitochondrial Genome

The mitochondrion stands at the center of cell energy metabolism. It contains its own genome, the mtDNA, that is a relic of its prokaryotic symbiotic ancestor. In plants, the mitochondrial genetic information influences important agronomic traits including fertility, plant vigor, chloroplast function, and cross-compatibility. Plant mtDNA has remarkable characteristics: It is much larger than the mtDNA of other eukaryotes and evolves very rapidly in structure. This is because of recombination activities that generate alternative mtDNA configurations, an important reservoir of genetic diversity that promotes rapid mtDNA evolution. On the other hand, the high incidence of ectopic recombination leads to mtDNA instability and the expression of gene chimeras, with potential deleterious effects. In contrast to the structural plasticity of the genome, in most plant species the mtDNA coding sequences evolve very slowly, even if the organization of the genome is highly variable. Repair mechanisms are probably responsible for such low mutation rates, in particular repair by homologous recombination. Herein we review some of the characteristics of plant organellar genomes and of the repair pathways found in plant mitochondria. We further discuss how homologous recombination is involved in the evolution of the plant mtDNA.

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

Chloroplast competition is controlled by lipid biosynthesis in evening primroses

In most eukaryotes, organellar genomes are transmitted preferentially by the mother, but molecular mechanisms and evolutionary forces underlying this fundamental biological principle are far from understood. It is believed that biparental inheritance promotes competition between the cytoplasmic organelles and allows the spread of so-called selfish cytoplasmic elements. Those can be, for example, fast-replicating or aggressive chloroplasts (plastids) that are incompatible with the hybrid nuclear genome and therefore maladaptive. Here we show that the ability of plastids to compete against each other is a metabolic phenotype determined by extremely rapidly evolving genes in the plastid genome of the evening primrose Oenothera Repeats in the regulatory region of accD (the plastid-encoded subunit of the acetyl-CoA carboxylase, which catalyzes the first and rate-limiting step of lipid biosynthesis), as well as in ycf2 (a giant reading frame of still unknown function), are responsible for the differences in competitive behavior of plastid genotypes. Polymorphisms in these genes influence lipid synthesis and most likely profiles of the plastid envelope membrane. These in turn determine plastid division and/or turnover rates and hence competitiveness. This work uncovers cytoplasmic drive loci controlling the outcome of biparental chloroplast transmission. Here, they define the mode of chloroplast inheritance, as plastid competitiveness can result in uniparental inheritance (through elimination of the "weak" plastid) or biparental inheritance (when two similarly "strong" plastids are transmitted).

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

Background

The 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.

Results

In 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.

Conclusions

Using 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.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1384-5) contains supplementary material, which is available to authorized users.

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

BACKGROUND

The 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.

RESULTS

In 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.

CONCLUSIONS

Using 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.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

The first next-generation sequencing approach to the mitochondrial phylogeny of African monogenean parasites (Platyhelminthes: Gyrodactylidae and Dactylogyridae)

BACKGROUND

Monogenean flatworms are the main ectoparasites of fishes. Representatives of the species-rich families Gyrodactylidae and Dactylogyridae, especially those infecting cichlid fishes and clariid catfishes, are important parasites in African aquaculture, even more so due to the massive anthropogenic translocation of their hosts worldwide. Several questions on their evolution, such as the phylogenetic position of Macrogyrodactylus and the highly speciose Gyrodactylus, remain unresolved with available molecular markers. Also, diagnostics and population-level research would benefit from the development of higher-resolution genetic markers. We aim to offer genetic resources for work on African monogeneans by providing mitogenomic data of four species (two belonging to Gyrodactylidae, two to Dactylogyridae), and analysing their gene sequences and gene order from a phylogenetic perspective.

RESULTS

Using Illumina technology, the first four mitochondrial genomes of African monogeneans were assembled and annotated for the cichlid parasites Gyrodactylus nyanzae, Cichlidogyrus halli, Cichlidogyrus mbirizei (near-complete mitogenome) and the catfish parasite Macrogyrodactylus karibae (near-complete mitogenome). Complete nuclear ribosomal operons were also retrieved, as molecular vouchers. The start codon TTG is new for Gyrodactylus and for Dactylogyridae, as is the incomplete stop codon TA for Dactylogyridae. Especially the nad2 gene is promising for primer development. Gene order was identical for protein-coding genes and differed between the African representatives of these families only in a tRNA gene transposition. A mitochondrial phylogeny based on an alignment of nearly 12,500 bp including 12 protein-coding and two ribosomal RNA genes confirms that the Neotropical oviparous Aglaiogyrodactylus forficulatus takes a sister group position with respect to the other gyrodactylids, instead of the supposedly 'primitive' African Macrogyrodactylus. Inclusion of the African Gyrodactylus nyanzae confirms the paraphyly of Gyrodactylus. The position of the African dactylogyrid Cichlidogyrus is unresolved, although gene order suggests it is closely related to marine ancyrocephalines.

CONCLUSIONS

The amount of mitogenomic data available for gyrodactylids and dactylogyrids is increased by roughly one-third. Our study underscores the potential of mitochondrial genes and gene order in flatworm phylogenetics, and of next-generation sequencing for marker development for these non-model helminths for which few primers are available.

Authentication Markers for Five Major Panax Species Developed via Comparative Analysis of Complete Chloroplast Genome Sequences

Ginseng represents a set of high-value medicinal plants of different species: Panax ginseng (Asian ginseng), Panax quinquefolius (American ginseng), Panax notoginseng (Chinese ginseng), Panax japonicus (Bamboo ginseng), and Panax vietnamensis (Vietnamese ginseng). Each species is pharmacologically and economically important, with differences in efficacy and price. Accordingly, an authentication system is needed to combat economically motivated adulteration of Panax products. We conducted comparative analysis of the chloroplast genome sequences of these five species, identifying 34-124 InDels and 141-560 SNPs. Fourteen InDel markers were developed to authenticate the Panax species. Among these, eight were species-unique markers that successfully differentiated one species from the others. We generated at least one species-unique marker for each of the five species, and any of the species can be authenticated by selection among these markers. The markers are reliable, easily detectable, and valuable for applications in the ginseng industry as well as in related research.

Localization of a defensive volatile 4-hydroxy-4-methylpentan-2-one in the capitate glandular trichomes of Oenothera glazioviana

Glandular trichomes of plants produce a wide variety of secondary metabolites which are considered as major defensive chemicals. The capitate glandular trichomes of Oenothera glazioviana (Onagraceae) were collected with laser microdissection and analyzed by gas chromatography-mass spectrometry. The volatile compound 4-hydroxy-4-methylpentan-2-one (1) was identified. We found that compound 1 displays antimicrobial, insecticidal, and phytotoxic activities. These results suggest that compound 1 might function as a defensive compound in the capitate glandular trichomes of O. glazioviana against pathogens, insect herbivores, and presumably competitive plants as well.

A Dynamic Tandem Repeat in Monocotyledons Inferred from a Comparative Analysis of Chloroplast Genomes in Melanthiaceae

Chloroplast genomes (cpDNA) are highly valuable resources for evolutionary studies of angiosperms, since they are highly conserved, are small in size, and play critical roles in plants. Slipped-strand mispairing (SSM) was assumed to be a mechanism for generating repeat units in cpDNA. However, research on the employment of different small repeated sequences through SSM events, which may induce the accumulation of distinct types of repeats within the same region in cpDNA, has not been documented. Here, we sequenced two chloroplast genomes from the endemic species Heloniopsis tubiflora (Korea) and Xerophyllum tenax (USA) to cover the gap between molecular data and explore "hot spots" for genomic events in Melanthiaceae. Comparative analysis of 23 complete cpDNA sequences revealed that there were different stages of deletion in the rps16 region across the Melanthiaceae. Based on the partial or complete loss of rps16 gene in cpDNA, we have firstly reported potential molecular markers for recognizing two sections (Veratrum and Fuscoveratrum) of Veratrum. Melathiaceae exhibits a significant change in the junction between large single copy and inverted repeat regions, ranging from trnH_GUG to a part of rps3. Our results show an accumulation of tandem repeats in the rpl23-ycf2 regions of cpDNAs. Small conserved sequences exist and flank tandem repeats in further observation of this region across most of the examined taxa of Liliales. Therefore, we propose three scenarios in which different small repeated sequences were used during SSM events to generate newly distinct types of repeats. Occasionally, prior to the SSM process, point mutation event and double strand break repair occurred and induced the formation of initial repeat units which are indispensable in the SSM process. SSM may have likely occurred more frequently for short repeats than for long repeat sequences in tribe Parideae (Melanthiaceae, Liliales). Collectively, these findings add new evidence of dynamic results from SSM in chloroplast genomes which can be useful for further evolutionary studies in angiosperms. Additionally, genomics events in cpDNA are potential resources for mining molecular markers in Liliales.

A naturally occurring conditional albino mutant in rice caused by defects in the plastid-localized OsABCI8 transporter

A wide range of molecules are transported across membranes by the ATP binding cassette (ABC) transporters. Plants possess a collection of ABC proteins bearing similarities to the components of prokaryotic multi subunit ABC transporters, designed as ABC group I. However the functions of most of them are not well understood. Here, we characterized a naturally occurring rice mutant that exhibited albino phenotype under continuous rainy days in the field, but gradually recovered to normal green after the rainy season. Molecular and genetic analyses revealed that the phenotypes were caused by a mutation in the OsABCI8 that encoded a member of the ABCI family. Subcellular localization demonstrated that OsABCI8 is a chloroplast ABC transporter. Expression of OsABCI8 is significantly enhanced in rainy days compared to sunny days. Besides defects in chloroplast development and chlorophyll biosynthesis, the mutant phenotype is accompanied by a higher accumulation of iron, suggesting that OsABCI8 is involved in iron transportation and/or homeostasis in rice. Our results demonstrate that OsABCI8 represents a conserved ABCI protein involved in transition metals transportation and/or homeostasis and suggest an important role of the plastid-localized OsABCI8 for chloroplast development.

What can we do with 1000 plastid genomes?

The plastid genome of plants is the smallest and most gene-rich of the three genomes in each cell and the one generally present in the highest copy number. As a result, obtaining plastid DNA sequence is a particularly cost-effective way of discovering genetic information about a plant. Until recently, the sequence information gathered in this way was generally limited to small portions of the genome amplified by polymerase chain reaction, but recent advances in sequencing technology have stimulated a substantial rate of increase in the sequencing of complete plastid genomes. Within the last year, the number of complete plastid genomes accessible in public sequence repositories has exceeded 1000. This sudden flood of data raises numerous challenges in data analysis and interpretation, but also offers the keys to potential insights across large swathes of plant biology. We examine what has been learnt so far, what more could be learnt if we look at the data in the right way, and what we might gain from the tens of thousands more genome sequences that will surely arrive in the next few years. The most exciting new discoveries are likely to be made at the interdisciplinary interfaces between molecular biology and ecology.

Plant Mitochondrial Genomes: Dynamics and Mechanisms of Mutation

The large mitochondrial genomes of angiosperms are unusually dynamic because of recombination activities involving repeated sequences. These activities generate subgenomic forms and extensive genomic variation even within the same species. Such changes in genome structure are responsible for the rapid evolution of plant mitochondrial DNA and for the variants associated with cytoplasmic male sterility and abnormal growth phenotypes. Nuclear genes modulate these processes, and over the past decade, several of these genes have been identified. They are involved mainly in pathways of DNA repair by homologous recombination and mismatch repair, which appear to be essential for the faithful replication of the mitogenome. Mutations leading to the loss of any of these activities release error-prone repair pathways, resulting in increased ectopic recombination, genome instability, and heteroplasmy. We review the present state of knowledge of the genes and pathways underlying mitochondrial genome stability.

Insights into the Existence of Isomeric Plastomes in Cupressoideae (Cupressaceae)

The cypress family (Cupressaceae) possesses highly rearranged plastomes that lack a pair of large inverted repeats typically found in land plants. A few cypress species have been reported to contain isomeric plastomes, but whether the existence of isomeric plastomes is ubiquitous in the family remains to be investigated with a broader taxon sampling. In this study, we sequenced the complete plastomes of ten species in Cupressoideae, the largest cypress subfamily. Cupressoideae showed relatively accelerated rates of substitutions at both nonsynonymous and synonymous sites as compared with other subfamilies of Cupressaceae. Our PCR and read mapping analyses together suggested the existence of isomeric plastomes in eight of the ten sequenced Cupressoideae species. The isomeric plastomes were also detected in 176 individuals from nine wild populations of four Cupressoideae species. Within Calocedrus macrolepis, we discovered a new type of isomeric plastomes that was likely derived from homologous recombination mediated by an 11-bp repeat. We conclude that isomeric plastomes are commonly present in Cupressoideae, thereby contributing to increased plastomic complexity.