Coalescent-based analysis distinguishes between allo- and autopolyploid origin in Shepherd's Purse (Capsella bursa-pastoris).

Polyploidization plays an important role in plant speciation. The most recent estimates report that up to 15% of angiosperm speciation events and 31% in ferns are accompanied by changes in ploidy level. Polyploids can arise either through autopolyploidy, when the sets of chromosomes originate from a single species, or through allopolyploidy, when they originate from different species. In this study, we used two different coalescent-based methods to determine the date and mode of the polyploidization event that led to the tetraploid cosmopolitan weed, Capsella bursa-pastoris. We sampled 78 C. bursa-pastoris accessions, and 53 and 43 accessions from the only two other members of this genus, C. grandiflora and C. rubella, respectively, and sequenced these accessions at 14 unlinked nuclear loci with locus-specific primers in order to be able to distinguish the two homeologues in the tetraploid. A large fraction of fixed differences between homeologous genes in C. bursa-pastoris are segregating as polymorphisms in C. grandiflora, consistent with an autopolyploid origin followed by disomic inheritance. To test this, we first estimated the demographic parameters of an isolation-with-migration model in a pairwise fashion between C. grandiflora and both genomes of C. bursa-pastoris and used these parameters in coalescent simulations to test the mode of origin of C. bursa-pastoris. Second, we used Approximate Bayesian Computation to compare an allopolyploid and an autopolyploid model. Both analyses led to the conclusion that C. bursa-pastoris originated less than 1 Ma by doubling of the C. grandiflora genome.

Reconstructing origins of loss of self-incompatibility and selfing in North American Arabidopsis lyrata: a population genetic context.

Theoretical and empirical comparisons of molecular diversity in selfing and outcrossing plants have primarily focused on long-term consequences of differences in mating system (between species). However, improving our understanding of the causes of mating system evolution requires ecological and genetic studies of the early stages of mating system transition. Here, we examine nuclear and chloroplast DNA sequences and microsatellite variation in a large sample of populations of Arabidopsis lyrata from the Great Lakes region of Eastern North American that show intra- and interpopulation variation in the degree of self-incompatibility and realized outcrossing rates. Populations show strong geographic clustering irrespective of mating system, suggesting that selfing either evolved multiple times or has spread to multiple genetic backgrounds. Diversity is reduced in selfing populations, but not to the extent of the severe loss of variation expected if selfing evolved due to selection for reproductive assurance in connection with strong founder events. The spread of self-compatibility in this region may have been favored as colonization bottlenecks following glaciation or migration from Europe reduced standing levels of inbreeding depression. However, our results do not suggest a single transition to selfing in this system, as has been suggested for some other species in the Brassicaceae.

Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size.

Recent studies comparing genome-wide polymorphism and divergence in Drosophila have found evidence for a surprisingly high proportion of adaptive amino acid fixations, but results for other taxa are mixed. In particular, few studies have found convincing evidence for adaptive amino acid substitution in plants. To assess the generality of this finding, we have sequenced 257 loci in the outcrossing crucifer Capsella grandiflora, which has a large effective population size and low population structure. Using a new method that jointly infers selective and demographic effects, we estimate that 40% of amino acid substitutions were fixed by positive selection in this species, and we also infer a low proportion of slightly deleterious amino acid mutations. We contrast these estimates with those for a similar data set from the closely related Arabidopsis thaliana and find significantly higher rates of adaptive evolution and fewer nearly neutral mutations in C. grandiflora. In agreement with results for other taxa, genes involved in reproduction show the strongest evidence for positive selection in C. grandiflora. Taken together, these results imply that both positive and purifying selection are more effective in C. grandiflora than in A. thaliana, consistent with the contrasting demographic history and effective population sizes of these species.

Signature of diversifying selection on members of the pentatricopeptide repeat protein family in Arabidopsis lyrata.

Pentatricopeptide repeat (PPR) proteins compose a family of nuclear-encoded transcriptional regulators of cytoplasmic genes. They have shown dramatic expansion in copy number in plants, and although the functional importance of many remains unclear, a subset has been repeatedly implicated as nuclear restorers for cytoplasmic male sterility. Here we investigate the molecular population genetics and molecular evolution of seven single-copy PPR genes in the outcrossing model plant Arabidopsis lyrata. In comparison with neutral reference loci, we find, on average, elevated levels of polymorphism and an excess of high-frequency variants at these PPR genes, suggesting that natural selection is maintaining polymorphism at some of these loci. This elevation in diversity persists when we control for divergence and generally decreases in the flanking regions, suggesting that these genes are themselves the targets of selection. Some of the PPR genes also demonstrate elevated population differentiation, which is consistent with spatially varying selection. In contrast, no comparable patterns are observed at these loci in A. thaliana, providing no evidence for the action of balancing selection in this selfing species. Taken together, these results suggest that a subset of PPR genes may be subject to balancing selection associated with ongoing cytonuclear coevolution in the outcrossing A. lyrata, which is possibly mediated either by intergenomic conflict or by compensatory evolution.

Recent speciation associated with the evolution of selfing in Capsella.

The evolution from outcrossing to predominant 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 marked reduction in flower size and a breakdown 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 place during the past 20,000 years. This transition was probably associated with a shift from stable equilibrium toward a near-complete population bottleneck causing a major reduction in effective 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.

Patterns of polymorphism and demographic history in natural populations of Arabidopsis lyrata.

BACKGROUND: Many of the processes affecting genetic diversity act on local populations. However, studies of plant nucleotide diversity have largely ignored local sampling, making it difficult to infer the demographic history of populations and to assess the importance of local adaptation. Arabidopsis lyrata, a self-incompatible, perennial species with a circumpolar distribution, is an excellent model system in which to study the roles of demographic history and local adaptation in patterning genetic variation. PRINCIPAL FINDINGS: We studied nucleotide diversity in six natural populations of Arabidopsis lyrata, using 77 loci sampled from 140 chromosomes. The six populations were highly differentiated, with a median FST of 0.52, and structure analysis revealed no evidence of admixed individuals. Average within-population diversity varied among populations, with the highest diversity found in a German population; this population harbors 3-fold higher levels of silent diversity than worldwide samples of A. thaliana. All A. lyrata populations also yielded positive values of Tajima's D. We estimated a demographic model for these populations, finding evidence of population divergence over the past 19,000 to 47,000 years involving non-equilibrium demographic events that reduced the effective size of most populations. Finally, we used the inferred demographic model to perform an initial test for local adaptation and identified several genes, including the flowering time gene FCA and a disease resistance locus, as candidates for local adaptation events. CONCLUSIONS: Our results underscore the importance of population-specific, non-equilibrium demographic processes in patterning diversity within A. lyrata. Moreover, our extensive dataset provides an important resource for future molecular population genetic studies of local adaptation in A. lyrata.

Selection on amino acid substitutions in Arabidopsis.

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.

Effective population size and tests of neutrality at cytoplasmic genes in Arabidopsis.

Cytoplasmic genomes typically lack recombination, implying that genetic hitch-hiking could be a predominant force structuring nucleotide polymorphism in the chloroplast and mitochondria. We test this hypothesis by analysing nucleotide polymorphism data at 28 loci across the chloroplast and mitochondria of the outcrossing plant Arabidopsis lyrata, and compare patterns with multiple nuclear loci, and the highly selfing Arabidopsis thaliana. The maximum likelihood estimate of the ratio of effective population size at cytoplasmic relative to nuclear genes in A. lyrata does not depart from the neutral expectation of 0.5. Similarly, the ratio of effective size in A. thaliana is close to unity, the neutral expectation for a highly selfing species. The results are thus consistent with neutral organelle polymorphism in these species or with comparable effects of hitch-hiking in both cytoplasmic and nuclear genes, in contrast to the results of recent studies on gynodioecious taxa. The four-gamete test and composite likelihood estimation provide evidence for very low levels of recombination in the organelles of A. lyrata, although permutation tests do not suggest that adjacent polymorphic sites are more closely linked than more distant sites across the two genomes, suggesting that mutation hotspots or very low rates of gene conversion could explain the data.

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata.

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.