Evolutionary rate variation at multiple levels of biological organization in plant mitochondrial DNA

Genome skimming and microsatellite analysis reveal contrasting patterns of genetic diversity in a rare sandhill endemic (Erysimum teretifolium, Brassicaceae)

Barriers between islands often inhibit gene flow creating patterns of isolation by distance. In island species, the majority of genetic diversity should be distributed among isolated populations. However, a self-incompatible mating system leads to higher genetic variation within populations and very little between-population subdivision. We examine these two contrasting predictions in Erysimum teretifolium, a rare self-incompatible plant endemic to island-like sandhill habitats in Santa Cruz County, California. We used genome skimming and nuclear microsatellites to assess the distribution of genetic diversity within and among eight of the 13 remaining populations. Phylogenetic analyses of the chloroplast genomes revealed a deep separation of three of the eight populations. The nuclear ribosomal DNA cistron showed no genetic subdivision. Nuclear microsatellites suggest 83% of genetic variation resides within populations. Despite this, 18 of 28 between-population comparisons exhibited significant population structure (mean F_ST = 0.153). No isolation by distance existed among all populations, however when one outlier population was removed from the analysis due to uncertain provenance, significant isolation by distance emerged (r² = 0.5611, p = 0.005). Population census size did not correlate with allelic richness as predicted on islands. Bayesian population assignment detected six genetic groupings with substantial admixture. Unique genetic clusters were concentrated at the periphery of the species’ range. Since the overall distribution of nuclear genetic diversity reflects E. tereifolium’s self-incompatible mating system, the vast majority of genetic variation could be sampled within any individual population. Yet, the chloroplast genome results suggest a deep split and some of the nuclear microsatellite analyses indicate some island-like patterns of genetic diversity. Restoration efforts intending to maximize genetic variation should include representatives from both lineages of the chloroplast genome and, for maximum nuclear genetic diversity, should include representatives of the smaller, peripheral populations.

The dual-targeted RNA editing factor AEF1 is universally conserved among angiosperms and reveals only minor adaptations upon loss of its chloroplast or its mitochondrial target

Upon loss of either its chloroplast or mitochondrial target, a uniquely dual-targeted factor for C-to-U RNA editing in angiosperms reveals low evidence for improved molecular adaptation to its remaining target. RNA-binding pentatricopeptide repeat (PPR) proteins specifically recognize target sites for C-to-U RNA editing in the transcriptomes of plant chloroplasts and mitochondria. Among more than 80 PPR-type editing factors that have meantime been characterized, AEF1 (or MPR25) is a special case given its dual targeting to both organelles and addressing an essential mitochondrial (nad5eU1580SL) and an essential chloroplast (atpFeU92SL) RNA editing site in parallel in Arabidopsis. Here, we explored the angiosperm-wide conservation of AEF1 and its two organelle targets. Despite numerous independent losses of the chloroplast editing site by C-to-T conversion and at least four such conversions at the mitochondrial target site in other taxa, AEF1 remains consistently conserved in more than 120 sampled angiosperm genomes. Not a single case of simultaneous loss of the chloroplast and mitochondrial editing target or of AEF1 disintegration or loss could be identified, contrasting previous findings for editing factors targeted to only one organelle. Like in most RNA editing factors, the PPR array of AEF1 reveals potential for conceptually "improved fits" to its targets according to the current PPR-RNA binding code. Surprisingly, we observe only minor evidence for adaptation to the mitochondrial target also after deep losses of the chloroplast target among Asterales, Caryophyllales and Poales or, vice versa, for the remaining chloroplast target after a deep loss of the mitochondrial target among Malvales. The evolutionary observations support the notion that PPR-RNA mismatches may be essential for proper function of editing factors.

Whole Plastome Sequencing Within Silene Section Psammophilae Reveals Mainland Hybridization and Divergence With the Balearic Island Populations

Reconstructing the phylogenetic relationships within Caryophyllaceae tribe Sileneae has been obscured by hybridization and incomplete lineage sorting. Silene is the largest genus in the Caryophyllaceae, and unraveling its evolutionary history has been particularly challenging. In order to infer the phylogenetic relationships among the five species in Silene section Psammophilae, we have performed a genome skimming approach to acquire the complete plastid genome (cpDNA), nuclear ribosomal cistron (nrDNA), and partial mitochondrial genome (mtDNA). We have included 26 populations, representing the range of each species' distribution. This section includes five morphologically similar species endemic to the Iberian Peninsula and Balearic Islands (Ibiza and Formentera), yet some of them occupy distinct edaphic habitats (e.g. maritime sands, calcareous sandstones). In addition to phylogeographic analyses, genetic structuring using the chloroplast data set was inferred with Discriminant Analysis of Principal Components (DAPC), analyses of molecular variance (AMOVA), and a partial Mantel test. Reference-guided assembly of 50 bp single-end and 250 bp paired-end Illumina reads produced the nearly complete cpDNA genome (154 kbp), partial mtDNA genome (from 81 to 114 kbp), and the nrDNA cistron (6.4 kbp). Selected variable regions of the cpDNA and mtDNA assemblies were confirmed by Sanger sequencing. Phylogenetic analyses of the mainland populations reveal incongruence among the three genomes. None of the three data sets produced relationships consistent with taxonomy or geography. In contrast, Silene cambessedesii, present in the Balearic Islands, is the only species that forms a strongly supported monophyletic clade in the cpDNA genome and is strongly differentiated with respect to the remaining taxa of the Iberian Peninsula. These results contrast with those obtained for mainland populations. Across the entire analysis, only one well-supported mainland clade of Silene littorea and Silene stockenii emerges from the southern region of the Iberian Peninsula. DAPC and AMOVA results suggest the absence of genetic structure among mainland populations of Silene section Psammophilae, whereas partial Mantel test discarded spatial correlation of genetic differentiation. The widespread incongruence between morphology-based taxonomic boundaries and phylogeography suggests a history of interspecific hybridization, in which only a substantial geographic barrier, like isolation by the Mediterranean Sea, was sufficient to create and maintain species boundaries in Silene section Psammophilae.

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

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

Maternal transmission, sex ratio distortion, and mitochondria

In virtually all multicellular eukaryotes, mitochondria are transmitted exclusively through one parent, usually the mother. In this short review, we discuss some of the major consequences of uniparental transmission of mitochondria, including deleterious effects in males and selection for increased transmission through females. Many of these consequences, particularly sex ratio distortion, have well-studied parallels in other maternally transmitted genetic elements, such as bacterial endosymbionts of arthropods. We also discuss the consequences of linkage between mitochondria and other maternally transmitted genetic elements, including the role of cytonuclear incompatibilities in maintaining polymorphism. Finally, as a case study, we discuss a recently discovered maternally transmitted sex ratio distortion in an insect that is associated with extraordinarily divergent mitochondria.

Nucleotide substitution analyses of the glaucophyte Cyanophora suggest an ancestrally lower mutation rate in plastid vs mitochondrial DNA for the Archaeplastida

A lot is known about the evolution and architecture of plastid, mitochondrial, and nuclear genomes, but surprisingly little is known about their relative rates of mutation. Most available relative-rate data come from seed plants, which, with few exceptions, have a mitochondrial mutation rate that is lower than those of the plastid and nucleus. But new findings from diverse plastid-bearing lineages have shown that for some eukaryotes the mitochondrial mutation rate is an order of magnitude greater than those of the plastid and nucleus. Here, we explore for the first time relative rates of mutation within the Glaucophyta-one of three main lineages that make up the Archaeplastida (or Plantae sensu lato). Nucleotide substitution analyses from distinct isolates of the unicellular glaucophyte Cyanophora paradoxa reveal 4-5-fold lower rates of mutation in the plastid and nucleus than the mitochondrion, which is similar to the mutational pattern observed in red algae and haptophytes, but opposite to that of seed plants. These data, together with data from previous reports, suggest that for much of the known photosynthetic eukaryotic diversity, plastid DNA mutations occur less frequently than those in mitochondrial DNA.

Comparative analyses of two Geraniaceae transcriptomes using next-generation sequencing

BACKGROUND

Organelle genomes of Geraniaceae exhibit several unusual evolutionary phenomena compared to other angiosperm families including accelerated nucleotide substitution rates, widespread gene loss, reduced RNA editing, and extensive genomic rearrangements. Since most organelle-encoded proteins function in multi-subunit complexes that also contain nuclear-encoded proteins, it is likely that the atypical organellar phenomena affect the evolution of nuclear genes encoding organellar proteins. To begin to unravel the complex co-evolutionary interplay between organellar and nuclear genomes in this family, we sequenced nuclear transcriptomes of two species, Geranium maderense and Pelargonium x hortorum.

RESULTS

Normalized cDNA libraries of G. maderense and P. x hortorum were used for transcriptome sequencing. Five assemblers (MIRA, Newbler, SOAPdenovo, SOAPdenovo-trans [SOAPtrans], Trinity) and two next-generation technologies (454 and Illumina) were compared to determine the optimal transcriptome sequencing approach. Trinity provided the highest quality assembly of Illumina data with the deepest transcriptome coverage. An analysis to determine the amount of sequencing needed for de novo assembly revealed diminishing returns of coverage and quality with data sets larger than sixty million Illumina paired end reads for both species. The G. maderense and P. x hortorum transcriptomes contained fewer transcripts encoding the PLS subclass of PPR proteins relative to other angiosperms, consistent with reduced mitochondrial RNA editing activity in Geraniaceae. In addition, transcripts for all six plastid targeted sigma factors were identified in both transcriptomes, suggesting that one of the highly divergent rpoA-like ORFs in the P. x hortorum plastid genome is functional.

CONCLUSIONS

The findings support the use of the Illumina platform and assemblers optimized for transcriptome assembly, such as Trinity or SOAPtrans, to generate high-quality de novo transcriptomes with broad coverage. In addition, results indicated no major improvements in breadth of coverage with data sets larger than six billion nucleotides or when sampling RNA from four tissue types rather than from a single tissue. Finally, this work demonstrates the power of cross-compartmental genomic analyses to deepen our understanding of the correlated evolution of the nuclear, plastid, and mitochondrial genomes in plants.

Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes

Plant mitochondrial genomes are usually transmitted to the progeny from the maternal parent. However, cases of paternal transmission are known and are perhaps more common than once thought. This review will consider recent evidence, both direct and indirect, of paternal transmission (leakage) of the mitochondrial genome of seed plants, especially in natural populations, and how this can result in offspring that carry a mixture of maternally and paternally derived copies of the genome; a type of heteroplasmy. It will further consider how this heteroplasmy facilitates recombination between genetically distinct partners; a process that can enhance mitochondrial genotypic diversity. This will then form the basis for a discussion of five evolutionary questions that arise from these observations. Questions include how plant mitochondrial genome evolution can be placed on a sexual to asexual continuum, whether cytoplasmic male sterility (CMS) facilitates the evolution of paternal leakage, whether paternal leakage is more likely in populations undergoing admixture, how leakage influences patterns of gene flow, and whether heteroplasmy occurs in natural populations at a frequency greater than predicted by crossing experiments. It is proposed that each of these questions offers fertile ground for future research on a diversity of plant species.

Disentangling the effects of mating systems and mutation rates on cytoplasmic [correction of cytoplamic] diversity in gynodioecious Silene nutans and dioecious Silene otites

Many flowering plant species exhibit a variety of distinct sexual morphs, the two most common cases being the co-occurrence of females and males (dioecy) or the co-occurrence of hermaphrodites and females (gynodioecy). In this study, we compared DNA sequence variability of the three genomes (nuclear, mitochondrial and chloroplastic) of a gynodioecious species, Silene nutans, with that of a closely related dioecious species, Silene otites. In the light of theoretical models, we expect cytoplasmic diversity to differ between the two species due to the selective dynamics that acts on cytoplasmic genomes in gynodioecious species: under an epidemic scenario, the gynodioecious species is expected to exhibit lower cytoplasmic diversity than the dioecious species, while the opposite is expected in the case of balancing selection maintaining sterility cytoplasms in the gynodioecious species. We found no difference between the species for nuclear gene diversity, but, for the cytoplasmic loci, the gynodioecious S. nutans had more haplotypes, and higher nucleotide diversity, than the dioecious relative, S. otites, even though the latter has a relatively high rate of mitochondrial synonymous substitutions, and therefore presumably a higher mutation rate. Therefore, as the mitochondrial mutation rate cannot account for the higher cytoplasmic diversity found in S. nutans, our findings support the hypothesis that gynodioecy in S. nutans has been maintained by balancing selection rather than by epidemic-like dynamics.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

Intraspecific variation in mitochondrial genome sequence, structure, and gene content in Silene vulgaris, an angiosperm with pervasive cytoplasmic male sterility

In angiosperms, mitochondrial-encoded genes can cause cytoplasmic male sterility (CMS), resulting in the coexistence of female and hermaphroditic individuals (gynodioecy). We compared four complete mitochondrial genomes from the gynodioecious species Silene vulgaris and found unprecedented amounts of intraspecific diversity for plant mitochondrial DNA (mtDNA). Remarkably, only about half of overall sequence content is shared between any pair of genomes. The four mtDNAs range in size from 361 to 429 kb and differ in gene complement, with rpl5 and rps13 being intact in some genomes but absent or pseudogenized in others. The genomes exhibit essentially no conservation of synteny and are highly repetitive, with evidence of reciprocal recombination occurring even across short repeats (< 250 bp). Some mitochondrial genes exhibit atypically high degrees of nucleotide polymorphism, while others are invariant. The genomes also contain a variable number of small autonomously mapping chromosomes, which have only recently been identified in angiosperm mtDNA. Southern blot analysis of one of these chromosomes indicated a complex in vivo structure consisting of both monomeric circles and multimeric forms. We conclude that S. vulgaris harbors an unusually large degree of variation in mtDNA sequence and structure and discuss the extent to which this variation might be related to CMS.

Relative rates of evolution among the three genetic compartments of the red alga Porphyra differ from those of green plants and do not correlate with genome architecture

In photosynthetic eukaryotes, relative silent-site nucleotide substitution rates (which can be used to approximate relative mutation rates) among mitochondrial, plastid, and nuclear genomes (mtDNAs, ptDNAs, and nucDNAs) are estimated to be 1:3:10 respectively for seed plants and roughly equal for green algae. These estimates correlate with certain genome characteristics, such as size and coding density, and have therefore been taken to support a relationship between mutation rate and genome architecture. Plants and green algae, however, represent a small fraction of the major eukaryotic plastid-bearing lineages. Here, we investigate relative rates of mutation within the model red algal genus Porphyra. In contrast to plants, we find that the levels of silent-site divergence between the Porphyra purpurea and Porphyra umbilicalis mtDNAs are three times that of their ptDNAs and five times that of their nucDNAs. Moreover, relative mutation rates do not correlate with genome architecture: despite an estimated three-fold difference in their mutation rate, the mitochondrial and plastid genome coding densities are equivalent - an observation that extends to organisms with secondary red algal plastids. These findings are supported by within-species silent-site polymorphism data from P. purpurea.

Mosaic origins of a complex chimeric mitochondrial gene in Silene vulgaris

Chimeric genes are significant sources of evolutionary innovation that are normally created when portions of two or more protein coding regions fuse to form a new open reading frame. In plant mitochondria astonishingly high numbers of different novel chimeric genes have been reported, where they are generated through processes of rearrangement and recombination. Nonetheless, because most studies do not find or report nucleotide variation within the same chimeric gene, evolution after the origination of these chimeric genes remains unstudied. Here we identify two alleles of a complex chimera in Silene vulgaris that are divergent in nucleotide sequence, genomic position relative to other mitochondrial genes, and expression patterns. Structural patterns suggest a history partially influenced by gene conversion between the chimeric gene and functional copies of subunit 1 of the mitochondrial ATP synthase gene (atp1). We identified small repeat structures within the chimeras that are likely recombination sites allowing generation of the chimera. These results establish the potential for chimeric gene divergence in different plant mitochondrial lineages within the same species. This result contrasts with the absence of diversity within mitochondrial chimeras found in crop species.

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

BACKGROUND

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

METHODOLOGY/PRINCIPAL FINDINGS

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

CONCLUSIONS/SIGNIFICANCE

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.

Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates

A pair of species within the genus Silene have evolved the largest known mitochondrial genomes, coinciding with extreme changes in mutation rate, recombination activity, and genome structure.

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure.

A fundamental challenge in evolutionary biology is to explain why organisms exhibit dramatic variation in genome size and complexity. One hypothesis predicts that high rates of mutation in DNA sequence create selection against large and complex genomes, which are more susceptible to mutational disruption. Species of flowering plants in the genus Silene vary by approximately 100-fold in the rates of mutation in their mitochondrial DNA, providing an excellent opportunity to test the predicted effects of high mutation rates on genome evolution. Contrary to expectation, Silene species with elevated mutation rates have experienced dramatic expansions in mitochondrial genome size compared to their slowly evolving relatives, resulting in the largest known mitochondrial genomes. In addition to the increases in size and mutation rate, these genomes also reveal a history of rapid change in genome structure. They have been fragmented into dozens of chromosomes and appear to have experienced major reductions in recombination activity. All of these changes have occurred in just the past few million years. This mitochondrial genome diversity within the genus Silene provides a striking example of rapid genomic change and raises new hypotheses regarding the relationship between mutation rate and genome evolution.

Recent acceleration of plastid sequence and structural evolution coincides with extreme mitochondrial divergence in the angiosperm genus Silene

The angiosperm genus Silene exhibits some of the most extreme and rapid divergence ever identified in mitochondrial genome architecture and nucleotide substitution rates. These patterns have been considered mitochondrial specific based on the absence of correlated changes in the small number of available nuclear and plastid gene sequences. To better assess the relationship between mitochondrial and plastid evolution, we sequenced the plastid genomes from four Silene species with fully sequenced mitochondrial genomes. We found that two species with fast-evolving mitochondrial genomes, S. noctiflora and S. conica, also exhibit accelerated rates of sequence and structural evolution in their plastid genomes. The nature of these changes, however, is markedly different from those in the mitochondrial genome. For example, in contrast to the mitochondrial pattern, which appears to be genome wide and mutationally driven, the plastid substitution rate accelerations are restricted to a subset of genes and preferentially affect nonsynonymous sites, indicating that altered selection pressures are acting on specific plastid-encoded functions in these species. Indeed, some plastid genes in S. noctiflora and S. conica show strong evidence of positive selection. In contrast, two species with more slowly evolving mitochondrial genomes, S. latifolia and S. vulgaris, have correspondingly low rates of nucleotide substitution in plastid genes as well as a plastid genome structure that has remained essentially unchanged since the origin of angiosperms. These results raise the possibility that common evolutionary forces could be shaping the extreme but distinct patterns of divergence in both organelle genomes within this genus.

Twenty-fold difference in evolutionary rates between the mitochondrial and plastid genomes of species with secondary red plastids

Within plastid-bearing species, the relative rates of evolution between mitochondrial and plastid genomes are poorly studied, but for the few lineages in which they have been explored, including land plants and green algae, the mitochondrial DNA mutation rate is nearly always estimated to be lower than or equal to that of the plastid DNA. Here, we show that in protists from three distinct lineages with secondary, red algal-derived plastids, the opposite is true: their mitochondrial genomes are evolving 5-30 times faster than their plastid genomes, even when the plastid is nonphotosynthetic. These findings have implications for understanding the origins and evolution of organelle genome architecture and the genes they encode.

De novo transcriptome assembly and polymorphism detection in the flowering plant Silene vulgaris (Caryophyllaceae)

Members of the angiosperm genus Silene are widely used in studies of ecology and evolution, but available genomic and population genetic resources within Silene remain limited. Deep transcriptome (i.e. expressed sequence tag or EST) sequencing has proven to be a rapid and cost-effective means to characterize gene content and identify polymorphic markers in non-model organisms. In this study, we report the results of 454 GS-FLX Titanium sequencing of a polyA-selected and normalized cDNA library from Silene vulgaris. The library was generated from a single pool of transcripts, combining RNA from leaf, root and floral tissue from three genetically divergent European subpopulations of S. vulgaris. A single full-plate 454 run produced 959,520 reads totalling 363.6 Mb of sequence data with an average read length of 379.0 bp after quality trimming and removal of custom library adaptors. We assembled 832,251 (86.7%) of these reads into 40,964 contigs, which have a total length of 25.4 Mb and can be organized into 18,178 graph-based clusters or 'isogroups'. Assembled sequences were annotated based on homology to genes in multiple public databases. Analysis of sequence variants identified 13,432 putative single-nucleotide polymorphisms (SNPs) and 1320 simple sequence repeats (SSRs) that are candidates for microsatellite analysis. Estimates of nucleotide diversity from 1577 contigs were used to generate genome-wide distributions that revealed several outliers with high diversity. All of these resources are publicly available through NCBI and/or our website (http://silenegenomics.biology.virginia.edu) and should provide valuable genomic and population genetic tools for the Silene research community.

Structural and content diversity of mitochondrial genome in beet: a comparative genomic analysis

Despite their monophyletic origin, mitochondrial (mt) genomes of plants and animals have developed contrasted evolutionary paths over time. Animal mt genomes are generally small, compact, and exhibit high mutation rates, whereas plant mt genomes exhibit low mutation rates, little compactness, larger sizes, and highly rearranged structures. We present the (nearly) whole sequences of five new mt genomes in the Beta genus: four from Beta vulgaris and one from B. macrocarpa, a sister species belonging to the same Beta section. We pooled our results with two previously sequenced genomes of B. vulgaris and studied genome diversity at the species level with an emphasis on cytoplasmic male-sterilizing (CMS) genomes. We showed that, contrary to what was previously assumed, all three CMS genomes belong to a single sterile lineage. In addition, the CMSs seem to have undergone an acceleration of the rates of substitution and rearrangement. This study suggests that male sterility emergence might have been favored by faster rates of evolution, unless CMS itself caused faster evolution.

Species delimitation and biogeography of two fir species (Abies) in central China: cytoplasmic DNA variation

It remains unclear how speciation history might contribute to species-specific variation and affect species delimitation. We examined concordance between cytoplasmic genetic variation and morphological taxonomy in two fir species, Abies chensiensis and A. fargesii, with overlapping distributions in central China. Range-wide genetic variation was investigated using mitochondrial (mt) and plastid (pt) DNA sequences, which contrast in their rates of gene flow. Four mtDNA haplotypes were recovered and showed no obvious species' bias in terms of relative frequency. In contrast, a high level of ptDNA variation was recorded in both species with 3 common ptDNA haplotypes shared between them and 21 rare ptDNA haplotypes specific to one or other species. We argue that the lack of concordance between morphological and molecular variation between the two fir species most likely reflects extensive ancestral polymorphism sharing for both forms of cytoplasmic DNA variation. It is feasible that a relatively fast mutation rate for ptDNA contributed to the production of many species-specific ptDNA haplotypes, which remained rare due to insufficient time passing for their spread and fixation in either species, despite high levels of intraspecific ptDNA gene flow. Our phylogeographic analyses further suggest that polymorphisms in both organelle genomes most likely originated during and following glacial intervals preceding the last glacial maximum, when species distributions became fragmented into several refugia and then expanded in range across central China.

Extensive loss of translational genes in the structurally dynamic mitochondrial genome of the angiosperm Silene latifolia

Background

Mitochondrial gene loss and functional transfer to the nucleus is an ongoing process in many lineages of plants, resulting in substantial variation across species in mitochondrial gene content. The Caryophyllaceae represents one lineage that has experienced a particularly high rate of mitochondrial gene loss relative to other angiosperms.

Results

In this study, we report the first complete mitochondrial genome sequence from a member of this family, Silene latifolia. The genome can be mapped as a 253,413 bp circle, but its structure is complicated by a large repeated region that is present in 6 copies. Active recombination among these copies produces a suite of alternative genome configurations that appear to be at or near "recombinational equilibrium". The genome contains the fewest genes of any angiosperm mitochondrial genome sequenced to date, with intact copies of only 25 of the 41 protein genes inferred to be present in the common ancestor of angiosperms. As observed more broadly in angiosperms, ribosomal proteins have been especially prone to gene loss in the S. latifolia lineage. The genome has also experienced a major reduction in tRNA gene content, including loss of functional tRNAs of both native and chloroplast origin. Even assuming expanded wobble-pairing rules, the mitochondrial genome can support translation of only 17 of the 61 sense codons, which code for only 9 of the 20 amino acids. In addition, genes encoding 18S and, especially, 5S rRNA exhibit exceptional sequence divergence relative to other plants. Divergence in one region of 18S rRNA appears to be the result of a gene conversion event, in which recombination with a homologous gene of chloroplast origin led to the complete replacement of a helix in this ribosomal RNA.

Conclusions

These findings suggest a markedly expanded role for nuclear gene products in the translation of mitochondrial genes in S. latifolia and raise the possibility of altered selective constraints operating on the mitochondrial translational apparatus in this lineage.

Extensive loss of RNA editing sites in rapidly evolving Silene mitochondrial genomes: selection vs. retroprocessing as the driving force

Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have "fast-forwarded" a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.

Testing for selection on synonymous sites in plant mitochondrial DNA: the role of codon bias and RNA editing

Since plant mitochondrial genomes exhibit some of the slowest known synonymous substitution rates, it is generally believed that they experience exceptionally low mutation rates. However, the use of synonymous substitution rates to infer mutation rates depends on the implicit assumption that synonymous sites are evolving neutrally (or nearly so). To assess the validity of this assumption in plant mitochondrial genomes, we examined coding sequence for footprints of selection acting at synonymous sites. We found that synonymous sites exhibit an AT rich and pyrimidine skewed nucleotide composition compared to both non-synonymous sites and non-coding regions. We also found some evidence for selection associated with both biased codon usage and conservation of regulatory sequences involved in mRNA processing, although some of these findings are subject to alternative non-adaptive interpretations. Regardless, the inferred strength of selection appears too weak to account for the variation in substitution rates between the mitochondrial genomes of plants and other multicellular eukaryotes. Therefore, these results are consistent with the interpretation that plant mitochondrial genomes experience a substantially lower mutation rate rather than increased functional constraints acting on synonymous sites. Nevertheless, there are important nucleotide composition patterns (particularly the differences between synonymous sites and non-coding DNA) that remain largely unexplained.

Gene flow and species delimitation: a case study of two pine species with overlapping distributions in southeast China

Species delimitation detected by molecular markers is complicated by introgression and incomplete lineage sorting between species. Recent modeling suggests that fixed genetic differences between species are highly related to rates of intraspecific gene flow. However, it remains unclear whether such differences are due to high levels of intraspecific gene flow overriding the spread of introgressed alleles or favoring rapid lineage sorting between species. In pines, chloroplast (cp) and mitochondrial (mt) DNAs are normally paternally and maternally inherited, respectively, and thus their relative rates of intraspecific gene flow are expected to be high and low, respectively. In this study, we used two pine species with overlapping geographical distributions in southeast China, P. massoniana and P. hwangshanensis, as a model system to examine the association between organelle gene flow and variation within and between species. We found that cpDNA variation across these two pine species is more species specific than mtDNA variation and almost delimits taxonomic boundaries. The shared mt/cp DNA genetic variation between species shows no bias in regard to parapatric versus allopatric species' distributions. Our results therefore support the hypothesis that high intraspecific gene flow has accelerated cpDNA lineage sorting between these two pine species.

Heteroplasmy and stoichiometric complexity of plant mitochondrial genomes--though this be madness, yet there's method in't

Mitochondrial heteroplasmy is defined as the coexistence of divergent mitochondrial genotypes in a cell. The ratio of the alternative genomes may be variable, but in plants, the usually prevalent main genome is accompanied by sublimons--substoichiometric mitochondrial DNA (mtDNA) molecules. Plant mitochondrial heteroplasmy was originally viewed as being associated with pathological mutations or was found in non-natural plant populations. Currently, it is considered to be a common situation in plants. Recent years have changed the previous view on the role of homologous recombination, small-scale mutations, and paternal leakage of mtDNA in the generation of heteroplasmy. Newly developed sensitive techniques have allowed the precise estimation of mtDNA stoichiometry. Mechanisms of maintenance and transmission of heteroplasmic genomes, including DNA recombination and replication, as well as mitochondrial fusion and fission, have been studied. This review describes the high level of plant mitochondrial genome complication--the 'madness' resulting from the heteroplasmic state and explains the method hidden in this madness. Heteroplasmy is described as the evolutionary strategy of uniparentally inherited plant mitochondrial genomes which do not undergo sexual recombination. In order to compensate for this deficiency, alternative types of mtDNA are substoichiometrically accumulated as a reservoir of genetic variability and may undergo accelerated evolution. Occasionally, sublimons are selected and amplified in the process called substoichiometric shifting, to take over the role of the main genome. Alternative mitochondrial genomes may recombine, yielding new mtDNA variants, or segregate during plant growth resulting in plants with mosaic phenotypes. Two opposite roles of mitochondrial heteroplasmy with respect to acceleration or counteracting of mutation accumulation are also discussed. Finally, nuclear control of heteroplasmy and substoichiometric shifting is described.

Phylogenetic analysis of mitochondrial substitution rate variation in the angiosperm tribe Sileneae

BACKGROUND

Recent phylogenetic studies have revealed that the mitochondrial genome of the angiosperm Silene noctiflora (Caryophyllaceae) has experienced a massive mutation-driven acceleration in substitution rate, placing it among the fastest evolving eukaryotic genomes ever identified. To date, it appears that other species within Silene have maintained more typical substitution rates, suggesting that the acceleration in S. noctiflora is a recent and isolated evolutionary event. This assessment, however, is based on a very limited sampling of taxa within this diverse genus.

RESULTS

We analyzed the substitution rates in 4 mitochondrial genes (atp1, atp9, cox3 and nad9) across a broad sample of 74 species within Silene and related genera in the tribe Sileneae. We found that S. noctiflora shares its history of elevated mitochondrial substitution rate with the closely related species S. turkestanica. Another section of the genus (Conoimorpha) has experienced an acceleration of comparable magnitude. The phylogenetic data remain ambiguous as to whether the accelerations in these two clades represent independent evolutionary events or a single ancestral change. Rate variation among genes was equally dramatic. Most of the genus exhibited elevated rates for atp9 such that the average tree-wide substitution rate for this gene approached the values for the fastest evolving branches in the other three genes. In addition, some species exhibited major accelerations in atp1 and/or cox3 with no correlated change in other genes. Rates of non-synonymous substitution did not increase proportionally with synonymous rates but instead remained low and relatively invariant.

CONCLUSION

The patterns of phylogenetic divergence within Sileneae suggest enormous variability in plant mitochondrial mutation rates and reveal a complex interaction of gene and species effects. The variation in rates across genomic and phylogenetic scales raises questions about the mechanisms responsible for the evolution of mutation rates in plant mitochondrial genomes.

Recent advances in the study of gynodioecy: the interface of theory and empiricism

BACKGROUND

In this review we report on recent literature concerned with studies of gynodioecy, or the co-occurrence of female and hermaphrodite individuals in natural plant populations. Rather than review this literature in its entirety, our focus is on the interplay between theoretical and empirical approaches to the study of gynodioecy.

SCOPE

Five areas of active inquiry are considered. These are the cost of restoration, the influence of population structure on spatial sex-ratio variation, the influence of inbreeding on sex expression, the signature of cyto-nuclear coevolution on the mitochondrial genome, and the consequences of mitochondrial paternal leakage.

CONCLUSIONS

Recent advances in the study of gynodioecy have been made by considering both the ecology of female:hermaphrodite fitness differences and the genetics of sex expression. Indeed theory has guided empiricism and empiricism has guided theory. Future advances will require that some of the methods currently available only for model organisms be applied to a wider range of species.

Silene as a model system in ecology and evolution

The genus Silene, studied by Darwin, Mendel and other early scientists, is re-emerging as a system for studying interrelated questions in ecology, evolution and developmental biology. These questions include sex chromosome evolution, epigenetic control of sex expression, genomic conflict and speciation. Its well-studied interactions with the pathogen Microbotryum has made Silene a model for the evolution and dynamics of disease in natural systems, and its interactions with herbivores have increased our understanding of multi-trophic ecological processes and the evolution of invasiveness. Molecular tools are now providing new approaches to many of these classical yet unresolved problems, and new progress is being made through combining phylogenetic, genomic and molecular evolutionary studies with ecological and phenotypic data.

The effect of breeding system on polymorphism in mitochondrial genes of Silene

Gynodioecy is a breeding system characterized by the co-occurrence of hermaphrodite and female individuals, generally as the result of nuclear-cytoplasmic interactions. The question remains whether the genetic factors controlling gynodioecy are maintained in species over long evolutionary timescales by balancing selection or are continually arising and being replaced in epidemic sweeps. If balancing selection maintains these factors, then neutral cytoplasmic diversity should be greater in gynodioecious than hermaphroditic species. In contrast, epidemic sweeps of factors controlling gynodioecy should decrease cytoplasmic diversity in gynodioecious relative to hermaphroditic species. We took a comparative approach in which we sequenced two mitochondrial genes, cytochrome b (cob) and cytochrome oxidase (cox1), for multiple populations of several hermaphroditic, gynodioecious, and dioecious species in the genus Silene. Breeding system was predictive of polymorphism. Gynodioecious species harbor many old haplotypes while hermaphroditic and dioecious species have little to no nucleotide diversity. The genealogical structure of neither gene departed from neutral expectations. Taken together, our results suggest that balancing selection acts on cytoplasmic male-sterility factors in several gynodioecious species in the genus.

Mitochondrial heteroplasmy and paternal leakage in natural populations of Silene vulgaris, a gynodioecious plant

It is currently thought that most angiosperms transmit their mitochondrial genomes maternally. Maternal transmission limits opportunities for genetic heterogeneity (heteroplasmy) of the mitochondrial genome within individuals. Recent studies of the gynodioecious species Silene vulgaris and Silene acaulis, however, document both direct and indirect evidence of mitochondrial heteroplasmy, suggesting that the mitochondrial genome is at times transmitted via paternal leakage. This heteroplasmy allows the generation of multi-locus recombinants, as documented in recent studies of both species. A prior study that employed quantitative PCR (q-PCR) on a limited sample provided direct evidence of heteroplasmy in the mitochondrial gene atp1 in S. vulgaris. Here, we apply the q-PCR methods to a much larger sample and extend them to incorporate the study of an additional atp1 haplotype along with two other haplotypes of the mitochondrial gene cox1 to evaluate the origin, extent, and transmission of mitochondrial genome heteroplasmy in S. vulgaris. We first calibrate our q-PCR methods experimentally and then use them to quantify heteroplasmy in 408 S. vulgaris individuals sampled from 22 natural populations located in Virginia, New York, and Tennessee. Sixty-one individuals exhibit heteroplasmy, including five that exhibited the joint heteroplasmy at both loci that is a prerequisite for effective recombination. The heteroplasmic individuals were distributed among 18 of the populations studied, demonstrating that heteroplasmy is a widespread phenomenon in this species. Further, we compare mother and offspring from 71 families to determine the rate of heteroplasmy gained and lost via paternal leakage and vegetative sorting across generations. Of 17 sibships exhibiting cox1 heteroplasmy and 14 sibships exhibiting atp1 heteroplasmy, more than half of the observations of heteroplasmy are generated via paternal leakage at the time of fertilization, with the rest being inherited from a heteroplasmic mother. Moreover, we show that the average paternal contribution during paternal leakage is about 12%. These findings are surprising, given that the current understanding of gynodioecy assumes that mitochondrial cytoplasmic male sterility elements are strictly maternally inherited. Knowledge of the dynamics of mitochondrial populations within individuals plays an important role in understanding the evolution of gynodioecy, and we discuss our findings within this context.

Tracking the evolution of a cold stress associated gene family in cold tolerant grasses

BACKGROUND

Grasses are adapted to a wide range of climatic conditions. Species of the subfamily Pooideae, which includes wheat, barley and important forage grasses, have evolved extreme frost tolerance. A class of ice binding proteins that inhibit ice re-crystallisation, specific to the Pooideae subfamily lineage, have been identified in perennial ryegrass and wheat, and these proteins are thought to have evolved from a leucine-rich repeat phytosulfokine receptor kinase (LRR-PSR)-like ancestor gene. Even though the ice re-crystallisation inhibition function of these proteins has been studied extensively in vitro, little is known about the evolution of these genes on the molecular level.

RESULTS

We identified 15 putative novel ice re-crystallisation inhibition (IRI)-like protein coding genes in perennial ryegrass, barley, and wheat. Using synonymous divergence estimates we reconstructed the evolution of the IRI-like gene family. We also explored the hypothesis that the IRI-domain has evolved through repeated motif expansion and investigated the evolutionary relationship between a LRR-domain containing IRI coding gene in carrot and the Pooideae IRI-like genes. Our analysis showed that the main expansion of the IRI-gene family happened ~36 million years ago (Mya). In addition to IRI-like paralogs, wheat contained several sequences that likely were products of polyploidisation events (homoeologs). Through sequence analysis we identified two short motifs in the rice LRR-PSR gene highly similar to the repeat motifs of the IRI-domain in cold tolerant grasses. Finally we show that the LRR-domain of carrot and grass IRI proteins both share homology to an Arabidopsis thaliana LRR-trans membrane protein kinase (LRR-TPK).

CONCLUSION

The diverse IRI-like genes identified in this study tell a tale of a complex evolutionary history including birth of an ice binding domain, a burst of gene duplication events after cold tolerant grasses radiated from rice, protein domain structure differentiation between paralogs, and sub- and/or neofunctionalisation of IRI-like proteins. From our sequence analysis we provide evidence for IRI-domain evolution probably occurring through increased copy number of a repeated motif. Finally, we discuss the possibility of parallel evolution of LRR domain containing IRI proteins in carrot and grasses through two completely different molecular adaptations.

Reticulate or tree-like chloroplast DNA evolution in Sileneae (Caryophyllaceae)?

Despite sampling of up to 25kb of chloroplast DNA sequence from 24 species in Sileneae a number of nodes in the phylogeny remain poorly supported and it is not expected that additional sequence sampling will converge to a reliable phylogenetic hypothesis in these parts of the tree. The main reason for this is probably a combination of rapid radiation and substitution rate heterogeneity. Poor resolution among closely related species are often explained by low levels of variation in chloroplast data, but the problem with our data appear to be high levels of homoplasy. Tree-like cpDNA evolution cannot be rejected, but apparent incongruent patterns between different regions are evaluated with the possibility of ancient interspecific chloroplast recombination as explanatory model. However, several major phylogenetic relationships, previously not recognized, are confidently resolved, e.g. the grouping of the two SW Anatolian taxa S. cryptoneura and S. sordida strongly disagrees with previous studies on nuclear DNA sequence data, and indicate a possible case of homoploid hybrid origin. The closely related S. atocioides and S. aegyptiaca form a sister group to Lychnis and the rest of Silene, thus suggesting that Silene may be paraphyletic, despite recent revisions based on molecular data.

Do recent findings in plant mitochondrial molecular and population genetics have implications for the study of gynodioecy and cytonuclear conflict?

The coexistence of females and hermaphrodites in plant populations, or gynodioecy, is a puzzle recognized by Darwin. Correns identified cytoplasmic inheritance of one component of sex expression, now known as cytoplasmic male sterility (CMS). Lewis established cytonuclear inheritance of gynodioecy as an example of genetic conflict. Although biologists have since developed an understanding of the mechanisms allowing the joint maintenance of CMS and nuclear male fertility restorer genes, puzzles remain concerning the inheritance of sex expression and mechanisms governing the origination of CMS. Much of the theory of gynodioecy rests on the assumption of maternal inheritance of the mitochondrial genome. Here we review recent studies of the genetics of plant mitochondria, and their implications for the evolution and transmission of CMS. New studies of intragenomic recombination provide a plausible origin for the chimeric ORFs that characterize CMS. Moreover, evidence suggests that nonmaternal inheritance of mitochondria may be more common than once believed. These findings may have consequences for the maintenance of cytonuclear polymorphism, mitochondrial recombination, generation of gynomonoecious phenotypes, and interpretation of experimental crosses. Finally we point out that CMS can alter the nature of the cytonuclear conflict that may have originally selected for uniparental inheritance.