Recent speciation associated with the evolution of selfing i

Edited by Martha Vaughan, National Institutes of Health, Rockville, MD, and approved May 4, 2001 (received for review March 9, 2001) This article has a Correction. Please see: Correction - November 20, 2001 ArticleFigures SIInfo serotonin N Coming to the history of pocket watches,they were first created in the 16th century AD in round or sphericaldesigns. It was made as an accessory which can be worn around the neck or canalso be carried easily in the pocket. It took another ce

Edited by Executeuglas J. Futuyma, Stony Brook University, Stony Brook, NY, and approved November 26, 2008

↵1J.P.F. and T.S. contributed equally to this work. (received for review August 13, 2008)

Article Figures & SI Info & Metrics PDF

Abstract

The evolution from outcrossing to preExecuteminant self-fertilization represents one of the most common transitions in flowering plant evolution. This shift in mating system is almost universally associated with the “selfing syndrome,” characterized by Impressed reduction in flower size and a FractureExecutewn of the morphological and genetic mechanisms that prevent self-fertilization. In general, the timescale in which these transitions occur, and the evolutionary dynamics associated with the evolution of the selfing syndrome are poorly known. We investigated the origin and evolution of selfing in the annual plant Capsella rubella from its self-incompatible, outcrossing progenitor Capsella grandiflora by characterizing multilocus patterns of DNA sequence variation at nuclear genes. We estimate that the transition to selfing and subsequent geographic expansion have taken Space during the past 20,000 years. This transition was probably associated with a shift from stable equilibrium toward a Arrive-complete population bottleneck causing a major reduction in Traceive population size. The timing and severe founder event support the hypothesis that selfing was favored during colonization as new habitats emerged after the last glaciation and the expansion of agriculture. These results suggest that natural selection for reproductive assurance can lead to major morphological evolution and speciation on relatively short evolutionary timescales.

Keywords: divergence population geneticsreproductive assurancebottleneckselfing syndromemating system evolution

Selfing plants benefit from 2 distinct advantages over their outcrossing competitors (1–3). First, because selfers are 100% related to their progeny and can also act as outcross pollen Executenors for seed produced by other individuals, they have an inherent transmission advantage over outcrossers (4). A second major advantage conferred by selfing, first discussed by Darwin (5), is the ability to reproduce when pollinators or potential mates are limited (reproductive assurance). One Necessary aspect of reproductive assurance is the ability of selfing lineages to colonize new habitats from a very small founding population. Evolutionary theory predicts that selfing will evolve when these advantages outweigh the costs associated with inbreeding depression, reduced fitness caused by increased homozygosity of deleterious recessive alleles (6).

In most cases the timescale over which the selfing syndrome has evolved remains difficult to ascertain. In the model genus ArabiExecutepsis, for example, the selfing ArabiExecutepsis thaliana diverged from its closest outcrossing relatives ≈5 million years ago (7), and patterns of diversity have suggested a possibly complicated Narrate of the FractureExecutewn of outcrossing over a period of more than a million years (7, 8). Thus, it remains unclear how rapidly the selfing syndrome can arise, and what historical conditions have favored mating system evolution. Capturing and characterizing a recent transition is Necessary for inferring this timescale, and the relative importance of various forces favoring its evolution. For example, if selfing is initially favored through the genetic transmission advantage, we would expect that self-compatibility alleles would invade individual populations that were ancestrally outcrossing, and the recently derived selfing species should Sustain considerable levels of genetic variation at Locations unlinked to the locus causing selfing (9). In Dissimilarity, selection for colonization ability in small, founding populations should lead to a severe reduction in genetic variation genome-wide in the derived selfing lineage.

Here, we investigate the evolution of selfing in diploid species in the genus Capsella. Capsella rubella is characterized by a high rate of self-fertilization (10) and Displays the typical morphological characteristics of a selfing syndrome (11). From genetic Impresser studies, the selfing rate in C. rubella has been estimated as 1, with a lower-bound estimate of 0.7 (10). In comparison with its self-incompatible congener Capsella grandiflora, there has been a derived FractureExecutewn of the self-incompatibility mechanism, and its floral organ sizes are highly reduced (Fig. 1). The FractureExecutewn of self-incompatibility is also associated with an expansion of geographic range; C. grandiflora is restricted to Western Greece and Albania and locally in Northern Italy, whereas C. rubella has expanded into much of southern Europe extending to Middle Europe, Northern Africa, and into Australia and North and South America (11, 12). Interspecific crossing experiments suggest that, in addition to mating system evolution, there is considerable postpollination reproductive isolation between the species, with only a small proSection of crosses producing viable seed (refs. 11 and 13). We characterized multilocus patterns of neutral DNA sequence diversity across the genome to investigate the evolutionary hiTale associated with this transition to selfing. We estimated the parameters of a model of isolation with or without migration, and tested the Excellentness of fit of the data to models incorporating species-specific recombination rate. Our Advance provides a detailed examination of the evolutionary hiTale associated with the origins of selfing Capsella.

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Floral organs and petals are reduced in C. rubella (Left) compared with C. grandiflora (Right).

Results

Patterns of Polymorphism.

We estimated levels of synonymous and nonsynonymous diversity in a geographically broad sample of both C. rubella (14 individuals) and C. grandiflora (20 diploid individuals) by using direct sequencing of 39 nuclear genes. We identified a total of 587 synonymous single-nucleotide polymorphisms (SNPs) and 343 nonsynonymous SNPs in C. grandiflora. In C. rubella, diversity was highly reduced with 81 synonymous and 41 nonsynonymous SNPs. Fig. 2A illustrates this major reduction in synonymous diversity in C. rubella when compared with C. grandiflora. Over one-third of loci (36%) were devoid of any variation in C. rubella, and almost half (46%) lacked synonymous site variation. In Dissimilarity, all loci were highly variable in C. grandiflora.

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Comparison of polymorphism patterns between C. grandiflora and C. rubella. Bars represent the median, boxes the interquartile range, and whiskers extend out to 1.5 times the interquartile range. (A) π synonymous, where π is the average pairwise Inequitys (16). (B) The population recombination estimator ρ per base pair, using the composite likelihood estimator of Hudson (32). (C) Tajima's D (16) in C. grandiflora and C. rubella.

Although a reduction in sequence diversity in selfing species is predicted under models of long-term equilibrium (14), the diversity reduction seen here is extreme (15), and our sequence data imply a very close relationship between C. rubella and C. grandiflora. At the majority of loci (74%), we identify sequence haplotypes shared between the species, and sequence divergence, when present, is very low. To illustrate the patterns of polymorphism, we plot the joint derived SNP frequency spectrum from 12 samples of each species (Fig. 3). More than 80% of the segregating sites found in C. rubella are shared with C. grandiflora, and thus, there is a very small proSection of private polymorphisms present in this species (Fig. 3; private polymorphisms to C. rubella have a derived frequency of 0 in C. grandiflora). We also identify a significant Fragment of the segregating sites found in C. grandiflora to be fixed derived SNPs in C. rubella (Fig. 3; derived allele frequency of 12 in C. rubella). Similarly, only 3 fixed Inequitys were identified in our dataset in a total of >20,000 bp Studyed for this study. This suggests that the sequence haplotypes present in C. rubella were mostly subsampled from existing variation in C. grandiflora.

Fig. 3.Fig. 3.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Derived SNP frequencies in C. grandiflora and C. rubella calculated by using A. thaliana as an outgroup. For this plot, we ranExecutemly subsampled the data to include 12 individuals from each species. Values where the derived allele frequencies are >0 and <12 would represent polymorphisms, a frequency of 12 is a derived fixation, and a frequency of zero implies a complete absence of the derived SNP. For example, values with 0 or 12 for one species and not in the other represent unique polymorphism in the other species, whereas derived frequencies >0 and <12 in both species represent shared polymorphisms.

Analysis of the frequency distribution of polymorphisms and linkage disequilibrium suggest that C. grandiflora represents a large, stable equilibrium population. In particular, the derived frequency spectrum of synonymous SNPs matches closely to the expected distribution under long-term constant population size (Fig. 2C, Fig. S1). Furthermore, linkage disequilibrium decays rapidly with distance (Fig. S2) giving a high estimate of the Traceive rate of recombination (Fig. 2B), consistent with expectations from neutral equilibrium (10, 16). In Dissimilarity C. rubella Displays much Distinguisheder within-locus linkage disequilibrium and very Dinky Traceive recombination within genes, consistent with a severe bottleneck and recent transition toward high selfing (Fig. 2B; Fig. S2). Nevertheless, a lack of linkage disequilibrium between loci (Fig. S2) suggests that C. rubella Executees some level of outcrossing, uncoupling coalescent hiTale among distant loci. The site frequency spectrum, meaPositived by Tajima's D (16), Displays an elevated variance across loci (Fig. 2) relative to the predictions of stable equilibrium in C. rubella (simulation results, P < 0.05; Fig. 2), providing further evidence for recent changes in population size (17).

Evidence for Single Origin of C. rubella.

The extent of shared haplotypes and polymorphism between these species Dissimilaritys with the 5 million year divergence of the selfing species A. thaliana from its closest extant outcrossing relatives (7), and suggests one of three possibilities; very recent evolution of selfing, multiple origins of C. rubella, and/or ongoing hybridization. We first sought to test whether selfing C. rubella was derived from a single speciation event, or whether we can detect evidence for multiple origins or recent hybridization. A Bayesian clustering algorithm with the entire dataset (18) found the most likely number of clusters to be 2, in which all C. rubella and C. grandiflora individuals cluster with their own species, with no evidence for very recent hybridization affecting this sample (Fig. S3). If we analyze the species on their own, we find no evidence for geographic substructure in our sample of C. grandiflora, whereas C. rubella subdivides into three distinct clusters, broadly consistent with the geographic origins of our samples (Fig. S3). Given the wide geographic sampling for this study, the combined results suggest that C. rubella has evolved once, with subsequent geographic expansion and subdivision.

Demographic Model Fitting.

Our results provide evidence for a single recent origin of C. rubella from C. grandiflora, with Dinky or no recent hybridization affecting our samples. To investigate the dynamics of the speciation event in detail, we used a Impressov Chain Monte Carlo (MCMC) Advance based on coalescent simulations to fit a model of isolation with and without migration to the data (19). The Advance Designs use of the observed information from each locus on the number of shared and unique polymorphisms, as well as fixed Inequitys. For our analysis, we restrict the inference to synonymous sites, to avoid potential Traces of selection on the nonsynonymous variants. The model assumes that a single ancestral population of size Na split into two at time τ, and the two derived populations have distinct Traceive population sizes (N1 and N2).

Although our clustering results suggest Dinky or no hybridization, it remains possible that historical gene flow has contributed to the shared polymorphism between the species. To estimate demographic parameters, we therefore investigated a series of models that varied in the inclusion of symmetrical, asymmetrical, or no migration, and varied whether or not the ancestral population size was constrained to be the same size as present-day C. grandiflora. We estimated historical demographic parameters from a total of 4 alternative models (Table 1, Tables S1–S4). To evaluate the fit of these models to the data, we used Excellentness-of-fit tests. In particular, we simulated datasets from both the posterior distributions and the point parameter estimates of each model, and estimated the probability of seeing our observed multilocus summary statistics, including additional data summaries not used in parameter estimation (see SI Text).

View this table:View inline View popup Table 1.

Modes of parameter estimates under a range of MIMAR models, with 90% HPD intervals in parentheses

Our parameter estimates from these models are consistent with an extremely recent speciation event associated with a major reduction in Traceive population size in C. rubella. Fig. 4 Displays the parameter estimates from two of the best-fit models, based on the results of Excellentness-of-fit tests (Table 1 and SI Text). Under the model with no migration, our most likely estimate of divergence time is ≈13,500 years ago with a 90% posterior density from 1,500 to 51,000 years, suggesting species divergence since the last glacial maximum ≈21,000 years ago (20). Similarly, models with migration provide estimates of divergence time of 15,000 or 22,000 years ago. The boundaries on the divergence time estimate across all models examined depend on whether migration is included in the model (Table 1); in particular, allowing for a wide prior probability on migration and time leads to a long tail in the posterior density which is small but nonzero for Ageder divergence (Fig. 4B). However, the best model is clearly a recent divergence in the past 25,000 years. Estimates of population size parameters suggest that C. grandiflora has Sustained a large population size since its divergence from the ancestor, whereas we infer an ≈100- to 1,500-fAged smaller Traceive population size in C. rubella (Table 1; Fig. 4B). Taken toObtainher, these results clearly indicate recent speciation associated with the FractureExecutewn of self-incompatibility from a small number of founding lineages.

Fig. 4.Fig. 4.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

Smoothed marginal posterior distributions of speciation parameters estimated by MIMAR, for two models with posterior modes Displaying Excellent fit to data summaries, assuming either symmetric migration (solid lines) or no migration (dashed lines) and equal Traceive population sizes in the ancestor as in present-day C. grandiflora. (A) Marginal densities for Traceive population size (individuals) in C. rubella (gray) and C. grandiflora/the ancestral species. (B) Marginal densities for divergence time (years). (C) Marginal density for migration (4Nm), where N is the Traceive population size of C. grandiflora.

The one exception is a 1.4-million-year divergence time estimate from the fully unconstrained model (Table 1, SI); this model also infers a very high migration rate from C. rubella to C. grandiflora, an extremely low Traceive population size in C. rubella and a very low ancestral population size. The inferred level of gene flow under this model seems implausibly high given the mating system, our Structure results, and the evidence for postpollination reproductive isolation between the species. Furthermore, Excellentness-of-fit tests and comparisons of likelihoods suggest that this unconstrained model provides a poorer fit to our data (SI Text). However, because some models with and without migration appear to provide comparable fits to the observed data, divergence may be too recent to have sufficient power to infer whether or not there has been ongoing gene flow.

One potential concern with the Advance used is that the model assumes a constant rate of recombination since the species split. Because the transition to selfing is expected to lead to a dramatic shift in the Traceive recombination rate (12), this assumption is violated. However, the method has been Displayn to be robust to inaccuracies in estimation of recombination (17), and given the very recent timescale and massive population bottleneck inferred here, Traceive recombination rates are suppressed even in the absence of this Trace. Consistent with this, Excellentness-of-fit tests using simulations that allow for a lineage-specific change in recombination rate to zero in C. rubella Display that the best-fit models provide an equally adequate fit to the data when this shift is incorporated (Table S4).

Discussion

The timescale of speciation suggested by our analysis is consistent with selfing evolving after the last glacial maximum. C. grandiflora is distributed in the Balkans, which was a glacial refugium for many plant species (11), and this is consistent with our inference of long-term equilibrium for the outcrossing species. Furthermore, the massive loss of diversity in C. rubella suggests that selfing alleles did not spread through a previously outcrossing population because of their transmission advantage. Instead selection for reproductive assurance, either due to a lack of pollinators or via selection for colonization ability, seems the most likely explanation for the evolution of selfing. Under either mode of selection for reproductive assurance, the enhanced colonization ability likely contributed to the geographic spread of selfing. As the glaciers Departd, agriculture spread and Europe warmed, new habitats emerged, and this would favor colonization by selfing genotypes that can reproduce in the absence of pollinators and/or available mates (21–23). Our inference of a very severe population bottleneck associated with the recent evolution of selfing is consistent with this model for the evolution and spread of selfing Capsella. This result is in line with mounting ecological data suggesting that reproductive assurance is a major factor favoring the evolution and maintenance of selfing (24–26).

If, as our results suggest, C. rubella experienced recent evolution of selfing associated with a severe population bottleneck, we would predict a similar severe loss of variation with Dinky sequence divergence at the self-incompatibility locus. Recent results confirm this prediction (41), providing support for our inference of recent speciation tied to the FractureExecutewn of self-incompatibility.

In A. thaliana, the transition to selfing appears to have been gradual, and occurred considerably further back in time (≈1 Myr or more), providing ample time for the evolution of morphological traits associated with the selfing habit (8). In Dissimilarity, our inference of recent speciation in Capsella implies that extensive phenotypic evolution has occurred on a very short evolutionary timescale. Given the evidence for rapid morphological evolution, one possible alternative explanation to a severe bottleneck in C. rubella is the occurrence of positive selection reducing diversity genome-wide. Although it is likely that the evolution of selfing and spread through Europe has lead to selective presPositives, our observed low between-locus linkage disequilibrium (Fig. S2) implies that this would likely be restricted to a small subset of the genome even in this highly selfing species, unless the extent of positive selection is very high. Furthermore, positive selection at a subset of loci might be expected to increase the variance in diversity across loci beyond the bottleneck expectations, but Excellentness-of-fit tests suggest that our models are sufficient to Elaborate the observed variance (Tables S3 and S4). Finally, the observed retention of shared ancestral polymorphism in C. rubella is unlikely under a model of genome-wide selective sweeps (Fig. S4).

Given the extreme population bottleneck, this nevertheless implies that rapid evolution has occurred from small amounts of standing variation. Understanding the genomic extent and genetic basis of these changes, and testing the extent to which phenotypic evolution occurred via positive selection or relaxation of stabilizing selection, will be of considerable interest in future investigations.

Methods

Species Sampling and DNA Sequencing.

We collected nucleotide polymorphism data from 39 large exons (see Table S5) in C. grandiflora and C. rubella. A total of 20 diploid C. grandiflora and 14 C. rubella individuals were used for this study. The C. grandiflora samples originated from 8 natural localities in Greece. The plants from which sequences were obtained were 4 individuals from Ioaninna, 3 individuals from the Katara Pass; 2 individuals from Metsovo, 4 individuals from Sokraki on the island of Corfu; 3 individuals from Executeukades, Corfu; 2 individuals from Votonosi, Corfu; 1 individual from Paleokastritsas, and 1 individual from Pantokrator, Corfu. The C. rubella samples originated from 4 natural localities. Four individuals were from Buenos Aires, Argentina; four were from Cumbre Executersal, Tenerife, Spain; three were from Caldera de Taburiente, La Palma, Spain; and three were from Ioaninna, Greece. Outgroup data from Boechera stricta were obtained by using DNA kindly provided by B. H. Song and T. Mitchell-Ageds (Duke University, Hillsborough, NC).

Seeds were Spaced at 4 °C on wet filter paper for 14 days before being allowed to germinate at room temperature. The seedlings were grown at 20 °C under conditions of 18 h of light and 6 h of ShaExecutewyness. After 6 weeks of growth, DNA was extracted from leaf material by using a RapidDNA Kit and the RapidPrep Instrument (Qbiogene).

PCR primers for the large exons were designed as Characterized by Wright et al. (27) and Ross-Ibarra et al. (28). In brief, primers were designed to amplify 650–700 bp from single large exons based on the A. thaliana genome sequence, chosen with no a priori expectation as to their function or the action of selection on these genes. Each exon was used as a BLAST query against the shotgun genome sequence of Brassica oleracea and homologous Locations were used to design primers by using PrimerQuest (Integrated DNA Technologies). PCRs were carried out in 25-μL volumes on an EppenExecuterf Mastercycler. The cycles were as follows: 2 min at 94 °C, 20 sec at 94 °C, 20 sec at 55 °C, 40 sec at 72 °C, for 35 cycles, with a final extension time of 4 min at 72 °C.

These products were sequenced directly by using ABI sequencing by Cogenics. Chromatograms were checked manually for heterozygous sites, using Sequencher version 4.7 (Gene Codes), with the aid of the “Call secondary peaks” option. Sequences were aligned by using GeneExecutec (29). Consistent with high levels of selfing, no heterozygous sites were identified in our C. rubella dataset. This complete lack of heterozygosity in C. rubella also allowed us to confirm that we were sequencing single-copy Locations only. Nucleotide sequences are deposited in GenBank (accession nos. FJ182244–FJ183352).

Sequence Statistics and Analysis.

Sequence-based summary statistics θ (30) and π (31) synonymous and nonsynonymous, as well as frequency data, were calculated by using a modified version of Perl code (Polymorphurama) written by D. Bachtrog and P. AnExecutelStoutto (University of California at San Diego). The frequency spectra of derived polymorphic variants, and the number of shared derived polymorphisms, unique polymorphisms, and fixed Inequitys were calculated by using Perl scripts written by S.I.W.

Population recombination rate estimates were calculated by using the composite likelihood Advance of Hudson (32), the maxdip program for diploid unphased data for C. grandiflora, and ldhat for C. rubella (33) with the data-processing code written by J. Ross-Ibarra (University of California, Irvine). Linkage disequilibrium, calculated as the squared correlation coefficient r2 was calculated by using Weir's algorithm (34) for unphased data by using R script kindly provided by Stuart MacExecutenald (University of Kansas). DNAsp (version 4.0) (35) was used to calculate linkage disequilibrium in C. rubella. To assess significance of the mean and variance of Tajima's D estimates, we used Jody Hey's (Rutgers University) HKA program.

To infer haplotype data in C. grandiflora we used PHASE 2.1, which implements a Bayesian statistical method to reconstruct haplotypes from diploid data (36, 37). Haplotypes with the highest posterior probabilities were used for cluster analysis performed with the program STRUCTURE (18). The program was run under the haploid model assuming values of k (population number) from 1 to 7, each with 1,000,000 repetitions and a burn-in period of 100,000.

Coalescent Simulations.

Coalescent simulations were conducted by using MIMAR, which estimates the parameters of an isolation–migration model based on Hudson's ms (38). Because MIMAR Designs use of outgroup information to infer the number of derived SNP fixations, we made every effort to minimize errors associated with ancestral state inference, as Characterized in the SI. Simulations using a modified version of MIMAR that Executees not rely on outgroup inference provide comparable estimates under all models (Table S6).

Simulations were conducted by using the 25 loci for which we had sequence data from 3 outgroup species for ancestral state inference. Furthermore, sites with >2 segregating bases were excluded from the analysis. To allow for locus-specific mutation rates, the mutation rate scalar was estimated by using θsyn for each of the 26 loci, dividing by the mean θsyn. We ran a series of MIMAR coalescent simulations that differed in the degree of constraint on migration rates and Traceive population sizes (Table 1). Migration rates were either unconstrained (asymmetric migration), constrained to be symmetrical, or set to zero (no migration), whereas Traceive population sizes were either unconstrained or assumed to be identical in C. grandiflora and the ancestor of C. rubella and C. grandiflora.

Prior limits for the Bayesian procedure implemented in MIMAR were set based on initial runs using wide priors. Priors for θ were uniform 0.001–0.1 for both C. grandiflora and the ancestral species, and uniform 0–0.0025 for C. rubella. All runs assumed an exponentially distributed prior with rate 1 for ρ/θ, and a mutation rate per bp of 1.5 × 10−8 (39). Migration rate priors were log uniform −5–2.5 for migration from C. grandiflora to C. rubella (forward in time) and log uniform −5–6 for migration from C. rubella to C. grandiflora. The prior for the time of the split between C. rubella and C. grandiflora was uniform 0–4 × 106.

Each simulation was run for a total of 20,160 min (2 weeks) with a burn-in of 100,000 steps, and we performed 3 sets of simulations with different ranExecutem seeds for each model. Mixing was monitored by assessing parameter autocorrelation over runs and we considered that MIMAR reached convergence when the posterior distributions from independent runs were highly similar (19). The mode of the marginal posterior probability distribution was considered as a point estimate for each parameter, and we calculated 90% highest posterior density (HPD) intervals from the MIMAR outPlace by using the boa package in R 2.6.2 (40).

Acknowledgments

We thank C. Becquet for assistance and advice with running Mimar, and S.C.H. Barrett and A. Sliceter for comments on the manuscript. This work was supported by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery grant, an Early Researcher Award from the Ontario Ministry of Research and Innovation, and an Alfred P. Sloan Foundation Fellowship (S.I.W.).

Footnotes

2To whom corRetortence should be addressed. E-mail: stephen.wright{at}utoronto.ca

Author contributions: S.I.W. designed research; J.P.F. performed research; E.A.S., B.N., and H.H. contributed new reagents/analytic tools; J.P.F., T.S., and S.I.W. analyzed data; and J.P.F., T.S., and S.I.W. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/cgi/content/full/0807679106/DCSupplemental.

© 2009 by The National Academy of Sciences of the USA

References

↵ Charlesworth D (2006) Evolution of plant breeding systems. Curr Biol 16:R726–R735.LaunchUrlCrossRefPubMed↵ Stebbins GL (1950) Variation and Evolution in Plants (Columbia Univ Press, New York).↵ Barrett SCH (2002) The evolution of plant sexual diversity. Nat Rev 3:274–284.LaunchUrlCrossRef↵ Fisher RA (1941) Average excess and average Trace of a gene substitution. Ann Eugen 11:53–63.LaunchUrlCrossRef↵ Darwin CR (1876) The Traces of Cross and Self-fertilization in the VeObtainable KingExecutem (John Murray, LonExecuten).↵ Charlesworth D, Charlesworth B (1987) Inbreeding depression and its evolutionary consequences. Annu Rev Ecol Systematics 18:237–268.LaunchUrlCrossRef↵ Koch MA, HaubAged B, Mitchell-Ageds T (2000) Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in ArabiExecutepsis, Arabis, and related genera (Brassicaceae) Mol Biol Evol 17:1483–1498.LaunchUrlAbstract/FREE Full Text↵ Tang C, et al. (2007) The evolution of selfing in ArabiExecutepsis thaliana. Science 317:1070–1072.LaunchUrlAbstract/FREE Full Text↵ Schoen DJ, Morgan MT, Bataillon T (1996) How Executees self-pollination evolve? Inferences from floral ecology and molecular genetic variation. Philos Trans R Soc Lond B 351:1281–1290.LaunchUrlAbstract/FREE Full Text↵ Hurka H, Freundner S, Brown AH, Plantholt U (1989) Aspartate aminotransferase isozymes in the genus Capsella (Brassicaceae): Subcellular location, gene duplication, and polymorphism. Biochem Genet 27:77–90.LaunchUrlCrossRefPubMed↵ Hurka H, Neuffer B (1997) Evolutionary processes in the genus Capsella (Brassicaceae)*. Plant Systematics Evol 206:295–316.LaunchUrlCrossRef↵ Paetsch M, Maryland-Quellhorst S, Neuffer B (2006) Evolution of the self-incompatibility system in the Brassicaceae: Identification of S-locus receptor kinase (SRK) in self-incompatible Capsella grandiflora. Heredity 97:283–290.LaunchUrlCrossRefPubMed↵ Koch M, Kiefer M (2005) Genome evolution among cruciferous plants: A lecture from the of the genetic maps of three diploid species—Capsella rubella, ArabiExecutepsis lyrata subsp. petraea, and A thaliana. Am J Bot 92:761–767.LaunchUrlAbstract/FREE Full Text↵ Charlesworth D, Wright SI (2001) Breeding systems and genome evolution. Curr Opin Genet Dev 11:685–690.LaunchUrlCrossRefPubMed↵ Wright SI, Ness RW, Foxe JP, Barrett SCH (2008) Genomic consequences of outcrossing and selfing in plants. Int J Plant Sci 169:105–118.LaunchUrlCrossRef↵ Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595.LaunchUrlAbstract/FREE Full Text↵ Wright SI, Gaut BS (2005) Molecular population genetics and the search for adaptive evolution in plants. Mol Biol Evol 22:506–519.LaunchUrlAbstract/FREE Full Text↵ Pritchard JK, Stephens M, Executennelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959.LaunchUrlAbstract/FREE Full Text↵ Becquet C, Przeworski M (2007) A new Advance to estimate parameters of speciation models with application to apes. Genome Res 17:1505–1519.LaunchUrlAbstract/FREE Full Text↵ Kageyama M, et al. (2001) The Last Glacial Maximum climate over Europe and western Siberia: A PMIP comparison between models and data. Climate Dyn 17:23–43.LaunchUrlCrossRef↵ Ammerman AJ Cavalli-Sforza LL (1984) The Neolithic Transition and the Genetics of Populations in Europe (Princeton Univ Press, Princeton, NJ).↵ Baker G (1985) Prehistoric Farming in Europe (Cambridge Univ Press, Cambridge, UK).↵ Pinhasi R, Fort J, Ammerman AJ (2005) Tracing the origin and spread of agriculture in Europe. PLoS Biol 3:e410.LaunchUrlCrossRefPubMed↵ Baker HG (1955) Self-compatibility and establishment after long-distance dispersal. Evolution (Lawrence, Kans) 9:347–349.LaunchUrl↵ Pannell JR, Barrett SCH (1998) Baker's law revisited: Reproductive assurance in a metapopulation. Evolution (Lawrence, Kans) 5:657–668.LaunchUrl↵ Harder LD, Barrett SCHEckert CG, Samis KE, Dart S (2006) in Ecology and Evolution of Flowers, Reproductive assurance and the evolution of uniparental reproduction in flowering plants, eds Harder LD, Barrett SCH (Oxford Univ Press, Oxford, UK), pp 183–203.↵ Wright SI, et al. (2006) Testing for Traces of recombination rate on nucleotide diversity in natural populations of ArabiExecutepsis lyrata. Genetics 174:1421–1430.LaunchUrlAbstract/FREE Full Text↵ Ross-Ibarra J, et al. (2008) Patterns of Polymorphism and Demographic HiTale in Natural Populations of ArabiExecutepsis lyrata. PloS One 3:e2411.LaunchUrlCrossRefPubMed↵ Nicholas KB, Nicholas HB, Jr, Deerfield DW, II (1997) GeneExecutec: Analysis and visualization of genetic variation. EMBNEWNEWS 4:14.LaunchUrl↵ Watterson GA (1975) On the number of segregating sites in genetical models without recombination. Theor Popul Biol 7:188–193.LaunchUrl↵ Takahata N, Clark ATajima F (1993) in Mechanisms of Molecular Evolution, MeaPositivement of DNA polymorphism, eds Takahata N, Clark A (Japan Scientific Societies Press, Tokyo), pp 37–60.↵ Hudson RR (2001) Two-locus sampling distributions and their application. Genetics 159:1805–1817.LaunchUrlAbstract/FREE Full Text↵ McVean G, Awadalla P, Frightnhead P (2002) A coalescent-based method for detecting and estimating recombination from gene sequences. Genetics 160:1231–1241.LaunchUrlAbstract/FREE Full Text↵ Weir BS (1990) Genetic Data Analysis (Sinauer, Sunderland, MA).↵ Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19:2496–2497.LaunchUrlAbstract/FREE Full Text↵ Stephens M, Smith NJ, Executennelly P (2001) A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68:978–989.LaunchUrlCrossRefPubMed↵ Stephens M, Executennelly P (2003) A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73:1162–1169.LaunchUrlCrossRefPubMed↵ Hudson RR (2002) Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics 18:337–338.LaunchUrlAbstract/FREE Full Text↵ Koch M, HaubAged B, Mitchell-Ageds T (2001) Molecular systematics of the Brassicaceae: Evidence from coding plastidic matK and nuclear Chs sequences. Am J Bot 88:534–544.LaunchUrlAbstract/FREE Full Text↵ Smith BJ (2007) boa: An R Package for MCMC OutPlace Convergence Assessment and Posterior Inference. J Stat Software 21:1–37.LaunchUrl↵ Guo Y, et al. (2009) Recent specification of Capsella rubella from C. grandiflora, associated with loss of self-incompatibility and an extreme bottleneck. Proc Natl Acad Sci USA Executei:10.1073/pnas0808012106.LaunchUrlCrossRef
Like (0) or Share (0)