A quantitative genetic analysis of nuclear-cytoplasmic male sterility in structured populations of Silene vulgaris

Chromosome-Level Genome Assembly for the Angiosperm Silene conica

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.

Chromosome-level genome assembly for the angiosperm Silene conica

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on PacBio HiFi, Hi-C and Bionano technologies. The assembly produced 10 scaffolds (one per chromosome) with a total length of 862 Mb and only ~1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.

Significance

Whole-genome sequences have been largely lacking for species in the genus Silene even though these flowering plants have been used for studying ecology, evolution, and genetics for over a century. Here, we address this gap by providing a high-quality nuclear genome assembly for S. conica, a species known to have greatly accelerated rates of sequence and structural divergence in its mitochondrial and plastid genomes. This resource will be valuable in understanding the coevolutionary interactions between nuclear and cytoplasmic genomes and in comparative analyses across this highly diverse genus.

Relative rates of evolution of male-beneficial nuclear compensatory mutations and male-harming Mother's Curse mitochondrial alleles

Mother's Curse alleles represent a significant source of potential male fitness defects. The maternal inheritance of mutations with the pattern of sex-specific fitness effects, s♀>0>s♂, allows Mother's Curse alleles to spread through a population even though they reduce male fitness. Although the mitochondrial genomes of animals contain only a handful of protein-coding genes, mutations in many of these genes have been shown to have a direct effect on male fertility. The evolutionary process of nuclear compensation is hypothesized to counteract the male-limited mitochondrial defects that spread via Mother's Curse. Here we use population genetic models to investigate the evolution of compensatory autosomal nuclear mutations that act to restore the loss of fitness caused by mitochondrial mutation pressures. We derive the rate of male fitness deterioration by Mother's Curse and the rate of restoration by nuclear compensatory evolution. We find that the rate of nuclear gene compensation is many times slower than that of its deterioration by cytoplasmic mutation pressure, resulting in a significant lag in the recovery of male fitness. Thus, the numbers of nuclear genes capable of restoring male mitochondrial fitness defects must be large in order to sustain male fitness in the face of mutation pressures.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Scale dependence of sex ratio in wild plant populations: implications for social selection

Social context refers to the composition of an individual's social interactants, including potential mates. In spatially structured populations, social context can vary among individuals within populations, generating the opportunity for social selection to drive differences in fitness functions among individuals at a fine spatial scale. In sexually polymorphic plants, the local sex ratio varies at a fine scale and thus has the potential to generate this opportunity. We measured the spatial distribution of two wild populations of the gynodioecious plant Silene vulgaris and show that there is fine‐scale heterogeneity in the local distribution of the sexes within these populations. We demonstrate that the largest variance in sex ratio is among nearest neighbors. This variance is greatly reduced as the spatial scale of social interactions increases. These patterns suggest the sex of neighbors has the potential to generate fine‐scale differences in selection differentials among individuals. One of the most important determinants of social interactions in plants is the behavior of pollinators. These results suggest that the potential for selection arising from sex ratio will be greatest when pollen is shared among nearest neighbors. Future studies incorporating the movement of pollinators may reveal whether and how this fine‐scale variance in sex ratio affects the fitness of individuals in these populations.

Molecular characterization and application of a novel cytoplasmic male sterility-associated mitochondrial sequence in rice

Background

Cytoplasmic male sterility (CMS) is a maternally inherited inability to produce functional pollen found in numerous flowering plant species. CMS is associated with mitochondrial DNA mutation, novel chimeric open reading frames (ORFs), and rearrangement of coding and noncoding regions of the mitochondrial genome.

Results

BLAST (Basic Local Alignment Search Tool) analysis indicated that L-sp1, a new sequence-characterized amplified region, is non-homologous to atp6-orfH79 (or atp6-orf79) and WA352 cloned CMS-associated genes. L-sp1 was found in 11 of 102 wild rice accessions belonging to four AA genome species: Oryza rufipogon, Oryza nivara, Oryza glumaepatula, and Oryza meridionalis. Using L-sp1, two new CMS lines were developed, from either low natural fertility plants or sterile plants, by backcrossing BC₁F₁ with Yuetai B. Northern blot and RT-PCR revealed that L-sp1 was only expressed in the anthers of w1/YTB, w2/YTB, w1/YTB//YTB, and w2/YTB//YTB when in the same cytoplasm background.

Conclusions

L-sp1 is a single-copy chimeric CMS-associated gene found in the mitochondrial genome. It can be expressed in anthers with the same specific cytoplasm background, and will be a useful molecular marker for the development and marker-assisted selection of new CMS lines.

Electronic supplementary material

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

Regular bottlenecks and restrictions to somatic fusion prevent the accumulation of mitochondrial defects in Neurospora

The replication and segregation of multi-copy mitochondrial DNA (mtDNA) are not under strict control of the nuclear DNA. Within-cell selection may thus favour variants with an intracellular selective advantage but a detrimental effect on cell fitness. High relatedness among the mtDNA variants of an individual is predicted to disfavour such deleterious selfish genetic elements, but experimental evidence for this hypothesis is scarce. We studied the effect of mtDNA relatedness on the opportunities for suppressive mtDNA variants in the fungus Neurospora carrying the mitochondrial mutator plasmid pKALILO. During growth, this plasmid integrates into the mitochondrial genome, generating suppressive mtDNA variants. These mtDNA variants gradually replace the wild-type mtDNA, ultimately culminating in growth arrest and death. We show that regular sequestration of mtDNA variation is required for effective selection against suppressive mtDNA variants. First, bottlenecks in the number of mtDNA copies from which a 'Kalilo' culture started significantly increased the maximum lifespan and variation in lifespan among cultures. Second, restrictions to somatic fusion among fungal individuals, either by using anastomosis-deficient mutants or by generating allotype diversity, prevented the accumulation of suppressive mtDNA variants. We discuss the implications of these results for the somatic accumulation of mitochondrial defects during ageing.

The roles of female and hermaphroditic flowers in the gynodioecious-gynomonoecious Silene littorea: insights into the phenology of sex expression

Some gynodioecious species have intermediate individuals that bear both female and hermaphroditic flowers. This phenomenon is known as a gynodioecious-gynomonoecious sexual system. Gender expression in such species has received little attention in the past, and the phenologies of male and female functions have also yet to be explored. In this study, we examined variations in gender patterns, their effects on female reproductive success and sex expression in depth throughout the flowering period in two populations. The studied populations of Silene littorea contained mostly gynomonoecious plants and the number of pure females was very low. The gynomonoecious plants showed high variability in the total proportion of female flowers. In addition, the proportion of female flowers in each plant varied widely across the flowering season. Although there was a trend towards maleness, our measures of functional gender suggested that most plants transmit their genes via both pollen and ovules. Fruit set and seed set were not significantly different among populations; in contrast, flower production significantly varied between the two populations - and among plants - with consequent variation in total seed production. Conversely, gender and sex expression were similar in both populations. Plants with higher phenotypic femaleness did not have higher fruit set, seed set or total female fecundity. The mating environment fluctuated little across the flowering period, but fluctuations were higher in the population with low flower production. We therefore conclude that the high proportion of gynomonoecious individuals in our studied populations of S. littorea may be advantageous for the species, providing the benefits of both hermaphroditic and female flowers.

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.

Genomic diversity in two related plant species with and without sex chromosomes--Silene latifolia and S. vulgaris

Background

Genome size evolution is a complex process influenced by polyploidization, satellite DNA accumulation, and expansion of retroelements. How this process could be affected by different reproductive strategies is still poorly understood.

Methodology/Principal Findings

We analyzed differences in the number and distribution of major repetitive DNA elements in two closely related species, Silene latifolia and S. vulgaris. Both species are diploid and possess the same chromosome number (2n = 24), but differ in their genome size and mode of reproduction. The dioecious S. latifolia (1C = 2.70 pg DNA) possesses sex chromosomes and its genome is 2.5× larger than that of the gynodioecious S. vulgaris (1C = 1.13 pg DNA), which does not possess sex chromosomes. We discovered that the genome of S. latifolia is larger mainly due to the expansion of Ogre retrotransposons. Surprisingly, the centromeric STAR-C and TR1 tandem repeats were found to be more abundant in S. vulgaris, the species with the smaller genome. We further examined the distribution of major repetitive sequences in related species in the Caryophyllaceae family. The results of FISH (fluorescence in situ hybridization) on mitotic chromosomes with the Retand element indicate that large rearrangements occurred during the evolution of the Caryophyllaceae family.

Conclusions/Significance

Our data demonstrate that the evolution of genome size in the genus Silene is accompanied by the expansion of different repetitive elements with specific patterns in the dioecious species possessing the sex chromosomes.

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.

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.

Genetic determination of male sterility in gynodioecious Silene nutans

Gynodioecy, the coexistence of female and hermaphrodite plants within a species, is often under nuclear-cytoplasmic sex determination, involving cytoplasmic male sterility (CMS) genes and nuclear restorers. A good knowledge of CMS and restorer polymorphism is essential for understanding the evolution and maintenance of gynodioecy, but reciprocal crossing studies remain scarce. Although mitochondrial diversity has been studied in a few gynodioecious species, the relationship between mitotype diversity and CMS status is poorly known. From a French sample of Silene nutans, a gynodioecious species whose sex determination remains unknown, we chose the four most divergent mitotypes that we had sampled at the cytochrome b gene and tested by reciprocal crosses whether they carry distinct CMS genes. We show that gynodioecy in S. nutans is under nuclear-cytoplasmic control, with at least two different CMSs and up to four restorers with epistatic interactions. Female occurrence and frequency were highly dependent on the mitotype, suggesting that the level of restoration varies greatly among CMSs. Two of the mitotypes, which have broad geographic distributions, represent different CMSs and are very unequally restored. We discuss the dynamics of gynodioecy at the large-scale meta-population level.

Genetical and molecular analysis reveals a cooperating relationship between cytoplasmic male sterility- and fertility restoration-related genes in Oryza species

Although the characterization of genes associated with cytoplasmic male sterility (CMS) and fertility restoration (Rf) has been well documented, the evolutionary relationship between nuclear Rf and CMS factors in mitochondria in Oryza species is still less understood. Here, 41 accessions from 7 Oryza species with AA genome were employed for analyzing the evolutionary relationships between the CMS factors and Rf candidates on chromosome 10. The phylogenetic tree based on restriction fragment length polymorphism patterns of CMS-associated mitochondrial genes showed that these 41 Oryza accessions fell into 3 distinct groups. Another phylogenetic tree based on PCR profiles of the nuclear Rf candidates on chromosome 10 was also established, and three groups were distinctively grouped. The accessions in each subgroup/group of the two phylogenetic trees are well parallel to each other. Furthermore, the 41 investigated accessions were test-crossed with Honglian (gametophytic type) and Wild-abortive (sporophytic type) CMS, and 5 groups were classified according to their restoring ability. The accessions in the same subgroup of the two phylogenetic trees shared similar fertility restoring pattern. Therefore, we conclude that the CMS-associated mitotypes are compatible to the Rf candidate-related nucleotypes, CMS and Rf have a parallel evolutionary relation in the Oryza species.

Experimental evidence for frequency dependent self-fertilization in the gynodioecious plant, Silene vulgaris

After over a half century of empirical and theoretical research regarding the evolution and maintenance of gynodioecy in plants, unexplored factors influencing the relative fitnesses of females and hermaphrodites remain. Theoretical studies suggest that hermaphrodite self-fertilization (selfing) rate influences the maintenance of gynodioecy and we hypothesized that population sex ratio may influence hermaphrodite selfing rate. An experimental test for frequency-dependent self-fertilization was conducted using replicated populations constructed with different sex ratios of the gynodioecious plant Silene vulgaris. We found that hermaphrodite selfing increased with decreased hermaphrodite frequency, whereas evidence for increased inbreeding depression was equivocal. We argue that incorporation of context dependent inbreeding into future models of the evolution of gynodioecy is likely to yield novel insights into sex ratio evolution.

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.

The ecology and genetics of a host shift: microbotryum as a model system

The need to prevent and cure emerging diseases often precludes their continuing study in situ. We present studies on the process of disease emergence by host shifts using the model system of anther-smut disease (Microbotryum violaceum) on the plant genus Silene (Caryophyllaceae). This system has little direct social impact, and it is readily amenable to experimental manipulation. Our microevolutionary studies have focused on the host shift of Microbotryum from Silene alba (=latifolia; white campion) onto Silene vulgaris (bladder campion) in a population in Virginia. Karyotypic variation shows that the host shift is recent and originates from the disease on sympatric S. alba. Analysis of the spatial pattern of disease shows that the host shift has been contingent on the co-occurrence of the two species at a local scale. Cross-inoculation studies show that families of the new host differ greatly in their susceptibility to the pathogen, indicating the potential for rapid evolution of resistance. Disease expression on the new host is frequently abnormal, suggesting that the pathogen is imperfectly adapted to its new host. In experimental populations, disease transmission within populations of the old host is greater than within populations of the new host. However, there is also a high transmission rate of the disease from the new host back to the old host, suggesting a feedback effect that increases disease prevalence in the community as a whole. Continuing studies of these populations are designed to determine whether this new host-pathogen system is likely to be self-sustaining and to quantify evolutionary changes in both the host and the pathogen.

Disentangling the causes of intrainflorescence variation in floral traits and fecundity in the hermaphrodite Silene acutifolia

Inflorescence architecture directly determines variations in floral traits and fecundity. Disentangling these patterns of variation is crucial to understanding intraplant variation, which sometimes is directly attributed to competition for resources with developing fruits. The dichasial cymes of Silene acutifolia were experimentally manipulated in the field to analyze whether the declines in petal size, ovule number, fruit set, and seed/ovule ratio along the inflorescence are constrained by ontogenetic development or are phenotypically plastic in response to environmental changes. At the same time, the level of pollen deficit was measured on different positions of the dichasia. The results showed clearly that all measured variables were more influenced by architecture than by resource competition with developing fruits; the removal of central (basal) and primary lateral flowers in the dichasia did not increase either the measures of floral characters or fecundity. On the other hand, although most of the decline in fecundity was due to architectural effects, there was also a pollen limitation, dependent to some degree on inflorescence position, which was probably due to lower pollen availability in the population when secondary flowers are in the female phase.

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.

Recombination and linkage disequilibrium among mitochondrial genes in structured populations of the gynodioecious plant Silene vulgaris

The impact of intergenic recombination on the population genetics of plant mitochondrial genomes is unknown. In an effort to study this in the gynodioecious plant Silene vulgaris three-locus PCR/RFLP genotypes (based on the mitochondrial genes atpA, cox1, and cob) were determined for 239 individuals collected from 20 North American populations. Seventeen three-locus PCR/RFLP genotypes were found. Recombination was indicated by observation of each of the four two-locus genotypes possible when the two most common alleles are considered for each of two loci. Based on these common alleles the absolute values of standardized linkage disequilibrium |D'| between pairs of loci range from 0.17 to 0.78. This indicates modest disequilibrium, rather than the maximum value expected in the absence of recombination |D'=1|, or the linkage equilibrium expected if recombination is pervasive (D'=0). Values of D' did not depend on which pair of loci contributed alleles to the analysis. The direction of D' obtained for the common atpA and cox1 alleles was comparable in sign and magnitude to that obtained by examining similar information obtained in a prior study of European samples. All three loci indicated a high degree of population structure (average FST=0.63), which would limit the within-population genetic diversity required for intergenic recombination to create novel genotypes, if most mating is local. Thus, population structure acts as a constraint on the approach to linkage equilibrium.

Sex-ratio evolution in nuclear-cytoplasmic gynodioecy when restoration is a threshold trait

Gynodioecious plant species, which have populations consisting of female and hermaphrodite individuals, usually have complex sex determination involving cytoplasmic male sterility (CMS) alleles interacting with nuclear restorers of fertility. In response to recent evidence, we present a model of sex-ratio evolution in which restoration of male fertility is a threshold trait. We find that females are maintained at low frequencies for all biologically relevant parameter values. Furthermore, this model predicts periodically high female frequencies (>50%) under conditions of lower female seed fecundity advantages (compensation, x = 5%) and pleiotropic fitness effects associated with restorers of fertility (costs of restoration, y = 20%) than in other models. This model explains the maintenance of females in species that have previously experienced invasions of CMS alleles and the evolution of multiple restorers. Sensitivity of the model to small changes in cost and compensation values and to initial conditions may explain why populations of the same species vary widely for sex ratio.

A cytonuclear incompatibility causes anther sterility in Mimulus hybrids

Multilocus interactions (also known as Dobzhansky-Muller incompatibilities) are thought to be the major source of hybrid inviability and sterility. Because cytoplasmic and nuclear genomes have conflicting evolutionary interests and are often highly coevolved, cytonuclear incompatibilities may be among the first to develop in incipient species. Here, we report the discovery of cytoplasm-dependent anther sterility in hybrids between closely related Mimulus species, outcrossing M. guttatus and selfing M. nasutus. A novel pollenless anther phenotype was observed in F2 hybrids with the M. guttatus cytoplasm (F2G) but not in the reciprocal F2N hybrids, F1 hybrids or parental genotypes. The pattern of phenotypic segregation in the F2G hybrids and two backcross populations fit a Mendelian single-locus recessive model, allowing us to map the underlying nuclear locus to a small region on LG7 of the Mimulus linkage map. Anther sterility was associated with a 20% reduction in flower size in backcross hybrids and we mapped a major cytoplasm-dependent corolla width QTL with its peak at the anther sterility locus. We argue that the cytonuclear anther sterility seen in hybrids reflects the presence of a cryptic cytoplasmic male sterility (CMS) and restorer system within the hermaphroditic M. guttatus population and therefore name the anther sterility locus restorer-of-male-fertility (RMF). The genetic mapping of RMF is a first step toward testing hypotheses about the molecular basis, individual fitness consequences, and ecological context of CMS and restoration in a system without stable CMS-restorer polymorphism (i.e., gynodioecy). The discovery of cryptic CMS in a hermaphroditic wildflower further suggests that selfish cytoplasmic evolution may play an important, but often undetected, role in shaping patterns of hybrid incompatibility and interspecific introgression in plants.

Maternal sex and mate relatedness affect offspring quality in the gynodioecious Silene acaulis

In gynodioecious species, females sacrifice fitness by not producing pollen, and hence must have a fitness advantage over hermaphrodites. Because females are obligately outcrossed, they may derive a fitness advantage by avoiding selfing and inbreeding depression. However, both sexes are capable of biparental inbreeding, and there are currently few estimates of the independent effects of maternal sex and multiple levels of inbreeding on female advantage. To test these hypotheses, females and hermaphrodites from six Alaskan populations of Silene acaulis were crossed with pollen from self (hermaphrodites only), a sibling, a random plant within the same population, and a plant from a different population. Germination, survivorship and early growth revealed inbreeding depression for selfs and higher germination but reduced growth in sib-crosses, relative to outcrosses. Independent of mate relatedness, females germinated more seeds that grew faster than offspring from hermaphrodites. This indicates that inbreeding depression as well as maternal sex can influence breeding system evolution. The effect of maternal sex may be explained by higher performance of female genotypes and a greater abundance of female genotypes among the offspring of female mothers.

Fine scale spatial structuring of sex and mitochondria in Silene vulgaris

Fine scale spatial structure (FSSS) of cytoplasmic genes in plants is thought to be generated via founder events and can be amplified when seeds germinate close to their mother. In gynodioecious species these processes are expected to generate FSSS in sex ratio because maternally inherited cytoplasmic male sterility genes partially influence sex expression. Here we document a striking example of FSSS in both mitochondrial genetic markers and sex in roadside populations of Silene vulgaris. We show that in one population FSSS of sexes influences relative fruit production of females compared to hermaphrodites. Furthermore, FSSS in sex ratio is expected to persist into future generations because offspring sex ratios from females are female-biased whereas offspring sex ratios from hermaphrodites are hermaphrodite-biased. Earlier studies indicated that pollen limitation is the most likely mechanism underlying negative frequency dependent fitness of females. Our results support the theoretical predictions that FSSS in sex ratio can reduce female fitness by decreasing the frequency at which females experience hermaphrodites. We argue that the influence of FSSS on female fitness is complementary to the influence of larger scale population structure on female fitness, and that population structure at both scales will act to decrease female frequencies in gynodioecious species. Better comprehension of the spatial structure of genders and genes controlling sex expression at a local scale is required for future progress toward understanding sex ratio evolution in gynodioecious plants.

Variable populations within variable populations: quantifying mitochondrial heteroplasmy in natural populations of the gynodioecious plant Silene vulgaris

Populations of mitochondria reside within individuals. Among angiosperms, these populations are rarely considered as genetically variable entities and typically are not found to be heteroplasmic in nature, leading to the widespread assumption that plant mitochondrial populations are homoplasmic. However, empirical studies of mitochondrial variation in angiosperms are relatively uncommon due to a paucity of sequence variation. Recent greenhouse studies of Silene vulgaris suggested that heteroplasmy might occur in this species at a level that it is biologically relevant. Here, we use established qualitative methods and a novel quantitative PCR method to study the intraindividual population genetics of mitochondria across two generations in natural populations of S. vulgaris. We show incidences of heteroplasmy for mitochondrial atpA and patterns of inheritance that are suggestive of more widespread heteroplasmy at both atpA and cox1. Further, our results demonstrate that quantitative levels of mitochondrial variation within individuals are high, constituting 26% of the total in one population. These findings are most consistent with a biparental model of mitochondrial inheritance. However, selection within individuals may be instrumental in the maintenance of variation because S. vulgaris is gynodioecious. Male sterility is, in part, regulated by the mitochondrial genome, and strong selection pressures appear to influence the frequency of females in these populations.

Heteroplasmy as a common state of mitochondrial genetic information in plants and animals

Plant and animal mitochondrial genomes, although quite distinct in size, structure, expression and evolutionary dynamics both may exhibit the state of heteroplasmy--the presence of more than one type of mitochondrial genome in an organism. This review is focused on heteroplasmy in plants, but we also highlight the most striking similarities and differences between plant and animal heteroplasmy. First we summarize the information on heteroplasmy generation and methods of its detection. Then we describe examples of quantitative changes in heteroplasmic populations of mitochondrial DNA (mtDNA) and consequences of such events. We also summarize the current knowledge about transmission and somatic segregation of heteroplasmy in plants and animals. Finally, factors which influence the stoichiometry of heteroplasmic mtDNA variants are discussed. Despite the apparent differences between the plant and animal heteroplasmy, the observed similarities allow one to conclude that this condition must play an important role in the mitochondrial biology of living organisms.

Inbreeding effects on progeny sex ratio and gender variation in the gynodioecious Silene vulgaris (Caryophyllaceae)

In gynodioecious species, sex expression is generally determined through cytoplasmic male sterility genes interacting with nuclear restorers of the male function. With dominant restorers, there may be an excess of females in the progeny of self-fertilized compared with cross-fertilized hermaphrodites. Moreover, the effect of inbreeding on late stages of the life cycle remains poorly explored. Here, we used hermaphrodites of the gynodioecious Silene vulgaris originating from three populations located in different valleys in the Alps to investigate the effects of two generations of self- and cross-fertilization on sex ratio and gender variation. We detected an increase in females in the progeny of selfed compared with outcrossed hermaphrodites and inbreeding depression for female and male fertility. Male fertility correlated positively with sex ratio differences between outbred and inbred progeny, suggesting that dominant restorers are likely to influence male fertility qualitatively and quantitatively in S. vulgaris. We argue that the excess of females in the progeny of selfed compared with outcrossed hermaphrodites and inbreeding depression for gamete production may contribute to the maintenance of females in gynodioecious populations of S. vulgaris because purging of the genetic load is less likely to occur.

Evidence for paternal transmission and heteroplasmy in the mitochondrial genome of Silene vulgaris, a gynodioecious plant

Gynodioecy refers to the co-occurrence of females and hermaphrodites in the same population. In many gynodioecious plants, sex is determined by an epistatic interaction between mitochondrial and nuclear genes, resulting in intragenomic evolutionary conflict, should the mitochondrial genome be maternally inherited. While maternal inheritance of the mitochondrial genome is common in angiosperms, few gynodioecious species have been studied. Here, the inheritance of the mitochondrial genes atpA and coxI was studied in 318 Silene vulgaris individuals distributed among 23 crosses. While maternal inheritance was indicated in 96% of the individuals studied, one or more individuals from each of four sib groups displayed a genotype that was identical to the father, or that did not match either parent. Given evidence that inheritance is not strictly maternal, it was hypothesized that some individuals could carry a mixture of maternally and paternally derived copies of the mitochondrial genome, a condition known as heteroplasmy. Since heteroplasmy might be difficult to detect should multiple versions of the mitochondrial genome co-occur in highly unequal copy number, a method was devised to amplify low-copy number forms of atpA differentially. Evidence for heteroplasmy was found in 23 of the 99 individuals studied, including cases in which the otherwise cryptic form of atpA matched the paternal genotype. The distribution of shared nucleotide sequence polymorphism among atpA haplotypes and the results of a population survey of the joint distribution of atpA and coxI haplotypes across individuals supports the hypothesis that heteroplasmy facilitates formation of novel mitochondrial genotypes by recombination.

Paternal leakage sustains the cytoplasmic polymorphism underlying gynodioecy but remains invasible by nuclear restorers

Cytoplasmic male sterility (CMS) in plants often results in gynodioecious populations, composed of hermaphrodites and male-sterile females. All models of gynodioecy assume maternal inheritance of the cytoplasmic alleles and postulate a variety of negatively frequency-dependent mechanisms to maintain the cytoplasmic polymorphisms observed in many natural populations. However, in some plant species, mitochondria are transmitted at least occasionally by pollen, a process called paternal leakage. We show that even a small amount of paternal leakage is sufficient to sustain a permanent, stable cytoplasmic polymorphism. Because only hermaphrodites provide pollen in gynodioecious species, the effects of paternal leakage are biased and occur more often from the non-CMS male-fertile haplotype to the CMS male-sterile haplotype. We also show that a nuclear restorer disrupts the polymorphic cytoplasmic equilibrium, leading to fixation of both the CMS allele and the restorer. Although a dominant nuclear restorer fixes, it fixes much more slowly than in the standard CMS models. Although a stable cytonuclear polymorphism is possible with "matching alleles" nuclear restoration, oscillations to low frequencies present a risk of loss by drift. Paternal leakage enhances the stability of joint cytonuclear polymorphism by reducing the chance that a CMS allele is lost by drift.

Genetics and adaptation in structured populations: sex ratio evolution in Silene vulgaris

Theoretical models suggest that population structure can interact with frequency dependent selection to affect fitness in such a way that adaptation is dependent not only on the genotype of an individual and the genotypes with which it co-occurs within populations (demes), but also the distribution of genotypes among populations. A canonical example is the evolution of altruistic behavior, where the costs and benefits of cooperation depend on the local frequency of other altruists, and can vary from one population to another. Here we review research on sex ratio evolution that we have conducted over the past several years on the gynodioecious herb Silene vulgaris in which we combine studies of negative frequency dependent fitness on female phenotypes with studies of the population structure of cytoplasmic genes affecting sex expression. This is presented as a contrast to a hypothetical example of selection on similar genotypes and phenotypes, but in the absence of population structure. Sex ratio evolution in Silene vulgaris provides one of the clearest examples of how selection occurs at multiple levels and how population structure, per se, can influence adaptive evolution.

Comparison between mitochondrial and chloroplast DNA variation in the native range of Silene vulgaris

A detailed survey of mitochondrial and chloroplast diversity in eight populations of Silene vulgaris from Central Europe was conducted for comparison with previously published data on diversity from S. vulgaris populations in the introduced range. Mitochondrial DNA (mtDNA) variation around the coxI gene was assessed with Southern blotting/restriction fragment length polymorphism methods. Chloroplast variation was assessed by sequencing the intergenic spacer separating the trnH and psbA genes. Thirty mtDNA haplotypes and 24 chloroplast DNA (cpDNA) haplotypes were found within 86 individuals. The overall genetic diversity h (0.941 for mitochondrial, and 0.893 for chloroplast markers) and within-population diversity were higher than reported in previous population studies of S. vulgaris in the USA and Europe. The frequency of private alleles was surprisingly high - more than 90% for both kinds of markers. Most of our populations were large and located in relatively undisturbed meadows, whereas surveys in Virginia consisted of smaller roadside populations. The slow rate of population turnover in European populations is discussed as a factor responsible for the relatively high diversity of S. vulgaris in undisturbed areas of its native range. Association between mtDNA and cpDNA haplotypes was also demonstrated. Finally, gender and mtDNA haplotype were associated in the Alps populations, where females were very rare.

Variation in restorer genes and primary sexual investment in gynodioecious Plantago coronopus: the trade-off between male and female function

In many gynodioecious species the nuclear inheritance of male fertility is complex and involves multiple (restorer) genes. In addition to restoring plants from the female (male sterile) to the hermaphrodite (male fertile) state, these genes are also thought to play a role in the determination of the quantity of pollen produced by hermaphrodites. The more restorer alleles a hermaphroditic plant possesses, the higher the pollen production. To test this hypothesis I combined the results of crossing studies of the genetics of male sterility with phenotypic data on investment in stamens and ovules among the progeny of plants involved in these studies. The sex ratio (i.e. the frequency of hermaphrodites among the progeny), being a measure of the number of restorer alleles of the maternal plant, was positively related to the investment in pollen (male function), but negatively related to the investment in ovules (female function), in both field and greenhouse experiments. Consequently, a negative correlation between male and female function was observed (trade-off) and it is suggested that antagonistic pleiotropic effects of restorer genes might be the cause. Phenotypic gender, a measure combining investment in both pollen and ovules, was highly repeatable between field and greenhouse, indicating genetic determination of a more male- or female-biased allocation pattern among the studied plants.

The hierarchical spatial distribution of chloroplast DNA polymorphism across the introduced range ofSilene vulgaris

Silene vulgaris was introduced into North America sometime prior to 1800. In order to document the population structure that has developed since that time, collections were made from 56 local populations distributed among 9 geographical regions in eastern North America. Individual plants were characterized for chloroplast DNA (cpDNA) haplotype by restriction fragment size analysis of four noncoding regions of cpDNA amplified by polymerase chain reaction. A total of 19 cpDNA haplotypes were detected using this method. The overall gene diversity of 0.85 is quite similar to the diversity detected in these same regions of cpDNA in a previously published sample of S. vulgaris taken from across much of Europe. The spatial distribution of the North American cpDNA diversity was quantified by hierarchical F-statistics that partitioned the genetic variance into variation among local populations within regions, and variation among regions. The average FST among populations within regions was 0.66 and the FST among regions was 0.09. The among-region variation was due to both differences among regions in the frequency of two most common haplotypes, and to the presence of a number of region-specific haplotypes. In order to test for isolation by distance at the regional level, FST values were calculated for all possible pairs of regions, and regressed against the geographical distance between those regions. There was no evidence for isolation by distance. It is suggested that the local population structure is generated by recent extinction/colonization dynamics, and that the among-region structure reflects demographic events associated with range expansion following introduction to North America.

A synthetic review of the theory of gynodioecy

Cytoplasmic male sterility alleles (CMS) and corresponding nuclear restorer alleles (R) determine gender expression in gynodioecious populations. In this article, we combine cytonuclear epistasis theory, multilevel selection theory, and resource allocation theory to partition the many selective forces acting in gynodioecious species. This description of selective forces reveals that the notion that genomic conflict leads to the spread of restorers is erroneous. By elucidating the effect of each selective force alone and in combination with the others, our approach shows why change in a single parameter, such as the degree of restorer dominance or the cost of restoration, can result in a myriad of opposing selective effects, making it difficult to test theoretical predictions in experimental systems. In particular, our approach allows us to identify the harmonic mean resource allocation between pollen and ovules and the attendant Fisherian selection as playing a critical role in the evolution of the restorer allele and gender polymorphism. We use this conceptual framework to propose empirical methods that emphasize the role of Fisherian selection acting in gynodioecious populations. Because our approach increases the number and specificity of theoretical predictions, we argue that it provides a better framework for empirical testing.

Modeling gynodioecy: novel scenarios for maintaining polymorphism

Nuclear-cytoplasmic gynodioecy is a breeding system of plants in which females and hermaphrodites co-occur in populations, and gender is jointly determined by cytoplasmic male sterility (CMS) genes and nuclear restorers of male fertility. Persistent polymorphism at both CMS and nuclear-restorer loci is necessary to maintain this breeding system. Theoretical models have explained how nuclear-cytoplasmic gynodioecy can be stable for certain assumptions. However, recent advances in our understanding of the genetics, population biology, and molecular mechanisms of sex determination in nuclear-cytoplasmic gynodioecious species suggest the utility of new models with different underlying assumptions. In this article, we examine different negative pleiotropic fitness effects of nuclear restorers (costs of restoration) using genetic and population assumptions based on recent literature. Specifically, we model populations with two CMS types and separate nuclear restorer loci for each CMS type. Under these assumptions, both overdominance for fitness and frequency-dependent selection at nuclear-restorer loci can support nuclear-cytoplasmic gynodioecy. Costs of restoration can be either dependent or independent of the cytoplasmic background. Seed fitness costs are more vulnerable to fixation of CMS types than pollen costs. Survivorship costs are effective at maintaining polymorphism even when total reproductive effects are low. Overall, our models display differences in the stability of nuclear-cytoplasmic gynodioecy and predicted population sex ratios that should be informative to researchers studying gynodioecy in the wild.

Associations among cytoplasmic molecular markers, gender, and components of fitness in Silene vulgaris, a gynodioecious plant

It has been suggested that the dynamics of chloroplast DNA (cpDNA) or mitochondrial DNA (mtDNA) genetic markers used in studies of plant populations could be influenced by natural selection acting elsewhere in the genome. This could be particularly true in gynodioecious plants if cpDNA or mtDNA genetic markers are in gametic disequilibrium with genes responsible for sex expression. In order to investigate this possibility, a natural population of the gynodioecious plant Silene vulgaris was used to study associations among mtDNA haplotype, cpDNA haplotype, sex and some components of fitness through seed. Individuals were sampled for mtDNA and cpDNA haplotype as determined using restriction fragment length polymorphism (RFLP) methods, sex (female or hermaphrodite), fruit number, fruit set, seeds/fruit and seed germination. The sex of surviving germinating seeds was also noted. All individuals in the population fell into one of two cytoplasmic categories, designated haplotypes f and g by a unique electrophoretic signature in both the mtDNA and cpDNA. The subset of the population carrying haplotype g included a significantly higher proportion of females when compared with the sex ratio of the subset carrying the f haplotype. Haplotype g had a significantly higher fitness when measured by fruit number, fruit set and seeds/fruit, whereas haplotype f had significantly higher fitness when measured by seed germination. Offspring of individuals carrying haplotype g included a significantly greater proportion of females when compared with offspring of individuals carrying the f haplotype. Other studies of gynodioecious plants have shown that females generally have higher fitness through seed than hermaphrodites, but in this study not all fitness differences between haplotypes could be predicted from differences in haplotype-specific sex ratio alone. Rather, some differences in haplotype-specific fitness were due to differences in fitness between individuals of the same sex, but carrying different haplotypes. The results are discussed with regard to the potential for hitchhiking selection to influence the dynamics of the noncoding regions used to designate the cpDNA and mtDNA haplotypes.

Karyotypic similarity identifies multiple host-shifts of a pathogenic fungus in natural populations

The detection of incipient host-shifts is important to the study of emergent diseases because it allows the examination of ecological and genetic conditions that favor novel inter-species transmission. Mixed populations of Silene latifolia and Silene vulgaris were investigated for the putative occurrence of host-shifts by the fungal plant pathogen Microbotryum violaceum (the cause of anther-smut disease) between S. latifolia (a common host for the pathogen) and S. vulgaris (a rare host). Samples of the fungus from mixed and pure host populations were studied for variation in their electrophoretic karyotypes. A karyotype distance matrix showed that fungal samples clustered by locality, but not by host species. Fungal samples from S. vulgaris were indistinguishable from sympatric samples from S. latifolia in multiple cases. The results indicated at least two independent host-shifts, one in the US and perhaps two in Italy. The karyotype and ecological data indicate that the direction of the host-shifts is from S. latifolia to S. vulgaris.

Genealogical evidence for epidemics of selfish genes

Some genetic elements spread infectiously in populations by increasing their rate of genetic transmission at the expense of other genes in the genome. These so-called selfish genetic elements comprise a substantial portion of eukaryotic genomes and have long been viewed as a potent evolutionary force. Despite this view, little is known about the evolutionary history of selfish genetic elements in natural populations, or their genetic effects on other portions of the genome. Here we use nuclear and chloroplast gene genealogies in two species of Silene to show the historical pattern of selection on a well known selfish genetic element, cytoplasmic male sterility. We provide evidence that evolution of cytoplasmic male sterility has been characterized by frequent turnovers of mutations in natural populations, thus supporting an epidemic model for the evolution of selfish genes, where new mutations repeatedly arise and rapidly sweep through populations.

Transmission rates and phenotypic effects of mitochondrial plasmids and cytotypes in Silene vulgaris

We investigated the transmission properties and the phenotypic effects of two mitochondrial plasmids in a population of the bladder campion, Silene vulgaris. In reciprocal crosses between plasmid-free and plasmid-carrying plants, no cases of paternal transmission or loss during maternal transmission were recorded. Neither was any transmission via pollen observed when plasmid-carrying plants of S. vulgaris were used to pollinate plasmid-free plants of the closely related species Silene uniflora. The phenotypic effects of the plasmids were investigated by comparing germination rate, early growth properties, and the gender of plants grown from seeds with and without plasmids. A significant association between plasmid status, on the one hand, and germination propensity and offspring gender, on the other, was found. However, because all plants carrying plasmids in the experiment shared the same cytoplasmic background, the exact contribution of the plasmid to the phenotypic variation could not be determined. Taken together, our experiments show that in S. vulgaris the mt-plasmids are not currently involved in any strong genetic conflict, but that they evolve in close association with their mitochondrial host.