Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae)

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

Demographic Processes Linked to Genetic Diversity and Positive Selection across a Species' Range

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range.

Impact of demography on linked selection in two outcrossing Brassicaceae species

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.

Inflorescence shoot elongation, but not flower primordia formation, is photoperiodically regulated in Arabidopsis lyrata

BACKGROUND AND AIMS

Photoperiod contains information about the progress of seasons. Plants use the changing photoperiod as a cue for the correct timing of important life history events, including flowering. Here the effect of photoperiod on flowering in four Arabidopsis lyrata populations originating from different latitudes was studied, as well as expression levels of candidate genes for governing the between-population differences.

METHODS

Flowering of plants from four A. lyrata populations was studied in three different photoperiods after vernalization. Flowering development was separated into three steps: flower primordia formation, inflorescence shoot elongation and opening of the first flower. Circadian expression rhythms of the A. lyrata homologues of GIGANTEA (GI), FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (FKF1), CONSTANS (CO) and FLOWERING LOCUS T (FT) were studied in three of the populations in the intermediate (14 h) photoperiod treatment.

KEY RESULTS

Most plants in all populations formed visible flower primordia during vernalization. Further inflorescence development after vernalization was strongly inhibited by short days in the northern European population (latitude 61°N), only slightly in the central European population (49°N) and not at all in the North American populations (36°N and 42°N). In the 14 h daylength, where all plants from the three southernmost populations but only 60 % of the northernmost population flowered, the circadian expression rhythm of the A. lyrata FT was only detected in the southern populations, suggesting differentiation in the critical daylength for activation of the long-day pathway. However, circadian expression rhythms of A. lyrata GI, FKF1 and CO were similar between populations.

CONCLUSIONS

The results indicate that in A. lyrata, transition to flowering can occur through pathways independent of long days, but elongation of inflorescences is photoperiodically regulated.

Local adaptation and ecological differentiation under selection, migration, and drift in Arabidopsis lyrata

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole-genome-based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4Ne m < 11) and have small effective population sizes (Ne <  6000), suggesting that both migration and drift can counteract local selection. However, the three-year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.

Genome-Wide Analysis of Colonization History and Concomitant Selection in Arabidopsis lyrata

The high climatic variability in the past hundred thousand years has affected the demographic and adaptive processes in many species, especially in boreal and temperate regions undergoing glacial cycles. This has also influenced the patterns of genome-wide nucleotide variation, but the details of these effects are largely unknown. Here we study the patterns of genome-wide variation to infer colonization history and patterns of selection of the perennial herb species Arabidopsis lyrata, in locally adapted populations from different parts of its distribution range (Germany, UK, Norway, Sweden, and USA) representing different environmental conditions. Using site frequency spectra based demographic modeling, we found strong reduction in the effective population size of the species in general within the past 100,000 years, with more pronounced effects in the colonizing populations. We further found that the northwestern European A. lyrata populations (UK and Scandinavian) are more closely related to each other than with the Central European populations, and coalescent based population split modeling suggests that western European and Scandinavian populations became isolated relatively recently after the glacial retreat. We also highlighted loci showing evidence for local selection associated with the Scandinavian colonization. The results presented here give new insights into postglacial Scandinavian colonization history and its genome-wide effects.

An Arabidopsis introgression zone studied at high spatio-temporal resolution: interglacial and multiple genetic contact exemplified using whole nuclear and plastid genomes

Background

Gene flow between species, across ploidal levels, and even between evolutionary lineages is a common phenomenon in the genus Arabidopsis. However, apart from two genetically fully stabilized allotetraploid species that have been investigated in detail, the extent and temporal dynamics of hybridization are not well understood. An introgression zone, with tetraploid A. arenosa introgressing into A. lyrata subsp. petraea in the Eastern Austrian Forealps and subsequent expansion towards pannonical lowlands, was described previously based on morphological observations as well as molecular data using microsatellite and plastid DNA markers. Here we investigate the spatio-temporal context of this suture zone, making use of the potential of next-generation sequencing and whole-genome data. By utilizing a combination of nuclear and plastid genomic data, the extent, direction and temporal dynamics of gene flow are elucidated in detail and Late Pleistocene evolutionary processes are resolved.

Results

Analysis of nuclear genomic data significantly recognizes the clinal structure of the introgression zone, but also reveals that hybridization and introgression is more common and substantial than previously thought. Also tetraploid A. lyrata and A. arenosa subsp. borbasii from outside the previously defined suture zone show genomic signals of past introgression. A. lyrata is shown to serve usually as the maternal parent in these hybridizations, but one exception is identified from plastome-based phylogenetic reconstruction. Using plastid phylogenomics with secondary time calibration, the origin of A. lyrata and A. arenosa lineages is pre-dating the last three glaciation complexes (approx. 550,000 years ago). Hybridization and introgression followed during the last two glacial-interglacial periods (since approx. 300,000 years ago) with later secondary contact at the northern and southern border of the introgression zone during the Holocene.

Conclusions

Footprints of adaptive introgression in the Northeastern Forealps are older than expected and predate the Last Glaciation Maximum. This correlates well with high genetic diversity found within areas that served as refuge area multiple times. Our data also provide some first hints that early introgressed and presumably preadapted populations account for successful and rapid postglacial re-colonization and range expansion.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4220-6) contains supplementary material, which is available to authorized users.

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson-Dobzhansky-Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.

On the young age of intraspecific herbaceous taxa

Dated phylogenies rarely include the divergence times of sister intraspecific taxa, and when they do little is said about this subject. We show that over 90% of the intraspecific plant taxa found in a literature search are estimated to be 5 million yr old or younger, with only 4% of taxa estimated to be over 10 million yr old or older. A Bayesian analysis of intraspecific taxon ages indicates that indeed these taxa are expected to be < 10 million yr old. This result for the young age of intraspecific taxa is consistent with the earlier observation that post-pollination reproductive barriers develop between 5 and 10 million yr after lineage splitting, thus leading to species formation. If lineages have not graduated to the species level of divergence by 10 million yr or so, they are likely to have gone extinct by that time as a result of narrow geographical distributions, narrow niche breadths, and relatively small numbers across populations.

Selection for population-specific adaptation shaped patterns of variation in the photoperiod pathway genes in Arabidopsis lyrata during post-glacial colonization

Spatially varying selection can lead to population-specific adaptation, which is often recognized at the phenotypic level; however, the genetic evidence is weaker in many groups of organisms. In plants, environmental shifts that occur due to colonization of a novel environment may require adaptive changes in the timing of growth and flowering, which are often governed by location-specific environmental cues such as day length. We studied locally varying selection in 19 flowering time loci in nine populations of the perennial herb Arabidopsis lyrata, which has a wide but patchy distribution in temperate and boreal regions of the northern hemisphere. The populations differ in their recent population demographic and colonization histories and current environmental conditions, especially in the growing season length. We searched for population-specific molecular signatures of directional selection by comparing a set of candidate flowering time loci with a genomic reference set within each population using multiple approaches and contrasted the patterns of different populations. The candidate loci possessed approximately 20% of the diversity of the reference loci. On average the flowering time loci had more rare alleles (a smaller Tajima's D) and an excess of highly differentiated sites relative to the reference, suggesting positive selection. The strongest signal of selection was detected in photoperiodic pathway loci in the colonizing populations of Northwestern Europe, whereas no evidence of positive selection was detected in the Central European populations. These findings emphasized the population-specific nature of selection and suggested that photoperiodic adaptation was important during postglacial colonization of the species.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

Cytoplasmic male sterility contributes to hybrid incompatibility between subspecies of Arabidopsis lyrata

In crosses between evolutionarily diverged populations, genomic incompatibilities may result in sterile hybrids, indicating evolution of reproductive isolation. In several plant families, crosses within a population can also lead to male sterile progeny because of conflict between the maternally and biparentally inherited genomes. We examined hybrid fertility between subspecies of the perennial outcrossing self-incompatible Lyrate rockcress (Arabidopsis lyrata) in large reciprocal F2 progenies and three generations of backcrosses. In one of the reciprocal F2 progenies, almost one-fourth of the plants were male-sterile. Correspondingly, almost one-half of the plants in one of the four reciprocal backcross progenies expressed male sterility. In an additional four independent F2 and backcross families, three segregated male sterility. The observed asymmetrical hybrid incompatibility is attributable to male sterility factors in one cytoplasm, for which the other population lacks effective fertility restorers. Genotyping of 96 molecular markers and quantitative trait locus mapping revealed that only 60% of the plants having the male sterile cytoplasm and lacking the corresponding restorers were phenotypically male-sterile. Genotyping data showed that there is only one restorer locus, which mapped to a 600-kb interval at the top of chromosome 2 in a region containing a cluster of pentatricopeptide repeat genes. Male fertility showed no trade-off with seed production. We discuss the role of cytoplasm and genomic conflict in incipient speciation and conclude that cytoplasmic male sterility–lowering hybrid fitness is a transient effect with limited potential to form permanent reproductive barriers between diverged populations of hermaphrodite self-incompatible species.

Investigating incipient speciation in Arabidopsis lyrata from patterns of transmission ratio distortion

Our understanding of the development of intrinsic reproductive isolation is still largely based on theoretical models and thorough empirical studies on a small number of species. Theory suggests that reproductive isolation develops through accumulation of epistatic genic incompatibilities, also known as Bateson-Dobzhansky-Muller (BDM) incompatibilities. We can detect these from marker transmission ratio distortion (TRD) in hybrid progenies of crosses between species or populations, where TRD is expected to result from selection against heterospecific allele combinations in hybrids. TRD may also manifest itself because of intragenomic conflicts or competition between gametes or zygotes. We studied early stage speciation in Arabidopsis lyrata by investigating patterns of TRD across the genome in F2 progenies of three reciprocal crosses between four natural populations. We found that the degree of TRD increases with genetic distance between crossed populations, but also that reciprocal progenies may differ substantially in their degree of TRD. Chromosomes AL6 and especially AL1 appear to be involved in many single- and two-locus distortions, but the location and source of TRD vary between crosses and between reciprocal progenies. We also found that the majority of single- and two-locus TRD appears to have a gametic, as opposed to zygotic, origin. Thus, while theory on BDM incompatibilities is typically illustrated with derived nuclear alleles proving incompatible in hybrid zygotes, our results suggest a prominent role for distortions emerging before zygote formation.