Selection at work in self-incompatible Arabidopsis lyrata. II. Spatial distribution of S haplotypes in Iceland

Breakdown of self-incompatibility due to genetic interaction between a specific S-allele and an unlinked modifier

Breakdown of self-incompatibility has frequently been attributed to loss-of-function mutations of alleles at the locus responsible for recognition of self-pollen (i.e. the S-locus). However, other potential causes have rarely been tested. Here, we show that self-compatibility of S1S1-homozygotes in selfing populations of the otherwise self-incompatible Arabidopsis lyrata is not due to S-locus mutation. Between-breeding-system cross-progeny are self-compatible if they combine S1 from the self-compatible cross-partner with recessive S1 from the self-incompatible cross-partner, but self-incompatible with dominant S-alleles. Because S1S1 homozygotes in outcrossing populations are self-incompatible, mutation of S1 cannot explain self-compatibility in S1S1 cross-progeny. This supports the hypothesis that an S1-specific modifier unlinked to the S-locus causes self-compatibility by functionally disrupting S1. Self-compatibility in S19S19 homozygotes may also be caused by an S19-specific modifier, but we cannot rule out a loss-of-function mutation of S19. Taken together, our findings indicate that breakdown of self-incompatibility is possible without disruptive mutations at the S-locus.

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.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.

Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation

Local adaptation is assumed to occur when populations differ in a phenotypic trait or a set of traits, and such variation has a genetic basis. Here, we introduce Arabidopsis halleri and its life history as a perennial model system to study population differentiation and local adaptation. Studies on altitudinal adaptation have been conducted in two regions: Mt. Ibuki in Japan and the European Alps. Several studies have demonstrated altitudinal adaptation in ultraviolet-B (UV-B) tolerance, leaf water repellency against spring frost and anti-herbivore defences. Studies on population differentiation in A. halleri have also focused on metal hyperaccumulation and tolerance to heavy metal contamination. In these study systems, genome scans to identify candidate genes under selection have been applied. Lastly, we briefly discuss how RNA-Seq can broaden phenotypic space and serve as a link to underlying mechanisms. In conclusion, A. halleri provides us with opportunities to study population differentiation and local adaptation, and relate these to the genetic systems underlying target functional traits.

The impact of self-incompatibility systems on the prevention of biparental inbreeding

Inbreeding in hermaphroditic plants can occur through two different mechanisms: biparental inbreeding, when a plant mates with a related individual, or self-fertilization, when a plant mates with itself. To avoid inbreeding, many hermaphroditic plants have evolved self-incompatibility (SI) systems which prevent or limit self-fertilization. One particular SI system-homomorphic SI-can also reduce biparental inbreeding. Homomorphic SI is found in many angiosperm species, and it is often assumed that the additional benefit of reduced biparental inbreeding may be a factor in the success of this SI system. To test this assumption, we developed a spatially-explicit, individual-based simulation of plant populations that displayed three different types of homomorphic SI. We measured the total level of inbreeding avoidance by comparing each population to a self-compatible population (NSI), and we measured biparental inbreeding avoidance by comparing to a population of self-incompatible plants that were free to mate with any other individual (PSI). Because biparental inbreeding is more common when offspring dispersal is limited, we examined the levels of biparental inbreeding over a range of dispersal distances. We also tested whether the introduction of inbreeding depression affected the level of biparental inbreeding avoidance. We found that there was a statistically significant decrease in autozygosity in each of the homomorphic SI populations compared to the PSI population and, as expected, this was more pronounced when seed and pollen dispersal was limited. However, levels of homozygosity and inbreeding depression were not reduced. At low dispersal, homomorphic SI populations also suffered reduced female fecundity and had smaller census population sizes. Overall, our simulations showed that the homomorphic SI systems had little impact on the amount of biparental inbreeding in the population especially when compared to the overall reduction in inbreeding compared to the NSI population. With further study, this observation may have important consequences for research into the origin and evolution of homomorphic self-incompatibility systems.

Evidence for Adaptive Introgression of Disease Resistance Genes Among Closely Related Arabidopsis Species

The generation and maintenance of functional variation in the pathogen defense system of plants is central to the constant evolutionary battle between hosts and parasites. If a species is susceptible to a given pathogen, hybridization and subsequent introgression of a resistance allele from a related species can potentially be an important source of new immunity and is therefore expected to be selected for in a process referred to as adaptive introgression. Here, we survey sequence variation in 10 resistance (R-) genes and compare them with 37 reference genes in natural populations of the two closely related and interfertile species: Arabidopsis lyrata and A. halleri. The R-genes are highly polymorphic in both species and show clear signs of trans-species polymorphisms. We show that A. lyrata and A. halleri have had a history of limited introgression for the reference genes. For the R-genes, the introgression rate has been significantly higher than for the reference genes, resulting in fewer fixed differences between species and a higher sharing of identical haplotypes. We conclude that R-genes likely cross the species boundaries at a higher rate than reference genes and therefore also that some of the increased diversity and trans-specific polymorphisms in R-genes is due to adaptive introgression.

Tissue Culture as a Source of Replicates in Nonmodel Plants: Variation in Cold Response in Arabidopsis lyrata ssp. petraea

While genotype-environment interaction is increasingly receiving attention by ecologists and evolutionary biologists, such studies need genetically homogeneous replicates-a challenging hurdle in outcrossing plants. This could be potentially overcome by using tissue culture techniques. However, plants regenerated from tissue culture may show aberrant phenotypes and "somaclonal" variation. Here, we examined somaclonal variation due to tissue culturing using the response to cold treatment of photosynthetic efficiency (chlorophyll fluorescence measurements for Fv/Fm, Fv'/Fm', and ΦPSII, representing maximum efficiency of photosynthesis for dark- and light-adapted leaves, and the actual electron transport operating efficiency, respectively, which are reliable indicators of photoinhibition and damage to the photosynthetic electron transport system). We compared this to variation among half-sibling seedlings from three different families of Arabidopsis lyrata ssp. petraea Somaclonal variation was limited, and we could detect within-family variation in change in chlorophyll fluorescence due to cold shock successfully with the help of tissue-culture derived replicates. Icelandic and Norwegian families exhibited higher chlorophyll fluorescence, suggesting higher performance after cold shock, than a Swedish family. Although the main effect of tissue culture on Fv/Fm, Fv'/Fm', and ΦPSII was small, there were significant interactions between tissue culture and family, suggesting that the effect of tissue culture is genotype-specific. Tissue-cultured plantlets were less affected by cold treatment than seedlings, but to a different extent in each family. These interactive effects, however, were comparable to, or much smaller than the single effect of family. These results suggest that tissue culture is a useful method for obtaining genetically homogenous replicates for studying genotype-environment interaction related to adaptively-relevant phenotypes, such as cold response, in nonmodel outcrossing plants.

What causes mating system shifts in plants? Arabidopsis lyrata as a case study

The genetic breakdown of self-incompatibility (SI) and subsequent mating system shifts to inbreeding has intrigued evolutionary geneticists for decades. Most of our knowledge is derived from interspecific comparisons between inbreeding species and their outcrossing relatives, where inferences may be confounded by secondary mutations that arose after the initial loss of SI. Here, we study an intraspecific breakdown of SI and its consequences in North American Arabidopsis lyrata to test whether: (1) particular S-locus haplotypes are associated with the loss of SI and/or the shift to inbreeding; (2) a population bottleneck may have played a role in driving the transition to inbreeding; and (3) the mutation(s) underlying the loss of SI are likely to have occurred at the S-locus. Combining multiple approaches for genotyping, we found that outcrossing populations on average harbour 5 to 9 S-locus receptor kinase (SRK) alleles, but only two, S1 and S19, are shared by most inbreeding populations. Self-compatibility (SC) behaved genetically as a recessive trait, as expected from a loss-of-function mutation. Bulked segregant analysis in SC × SI F2 individuals using deep sequencing confirmed that all SC plants were S1 homozygotes but not all S1 homozygotes were SC. This was also revealed in population surveys, where only a few S1 homozygotes were SC. Together with crossing data, this suggests that there is a recessive factor that causes SC that is physically unlinked to the S-locus. Overall, our results emphasise the value of combining classical genetics with advanced sequencing approaches to resolve long outstanding questions in evolutionary biology.

Contrasting patterns of genetic structuring in natural populations of Arabidopsis lyrata Subsp. petraea across different regions in northern Europe

Level and partitioning of genetic diversity is expected to vary between contrasting habitats, reflecting differences in strength of ecological and evolutionary processes. Therefore, it is necessary to consider processes acting on different time scales when trying to explain diversity patterns in different parts of species' distributions. To explore how historical and contemporary factors jointly may influence patterns of genetic diversity and population differentiation, we compared genetic composition in the perennial herb Arabidopsis lyrata ssp. petraea from the northernmost parts of its distribution range on Iceland to that previously documented in Scandinavia. Leaf tissue and soil were sampled from ten Icelandic populations of A. lyrata. Seedlings were grown from soil samples, and tissue from above-ground and seed bank individuals were genotyped with 21 microsatellite markers. Seed bank density in Icelandic populations was low but not significantly different from that observed in Norwegian populations. While within-population genetic diversity was relatively high on Iceland (H(E) = 0.35), among-population differentiation was low (F(ST) = 0.10) compared to Norwegian and Swedish populations. Population differentiation was positively associated with geographical distance in both Iceland and Scandinavia, but the strength of this relationship varied between regions. Although topography and a larger distribution range may explain the higher differentiation between mountainous Norwegian relative to lowland populations in Sweden, these factors cannot explain the lower differentiation in Icelandic compared to Swedish populations. We propose that low genetic differentiation among Icelandic populations is not caused by differences in connectivity, but is rather due to large historical effective population sizes. Thus, rather than contemporary processes, historical factors such as survival of Icelandic lineages in northern refugia during the last glacial period may have contributed to the observed pattern.

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

Gradual molecular evolution of a sex determination switch through incomplete penetrance of femaleness

Some genes regulate phenotypes that are either present or absent. They are often important regulators of developmental switches and are involved in morphological evolution. We have little understanding of the molecular mechanisms by which these absence/presence gene functions have evolved, because the phenotype and fitness of molecular intermediate forms are unknown. Here, we studied the sex-determining switch of 14 natural sequence variants of the csd gene among 76 genotypes of the honeybee (Apis mellifera). Heterozygous genotypes (different specificities) of the csd gene determine femaleness, while hemizygous genotypes (single specificity) determine maleness. Homozygous genotypes of the csd gene (same specificity) are lethal. We found that at least five amino acid differences and length variation between Csd specificities in the specifying domain (PSD) were sufficient to regularly induce femaleness. We estimated that, on average, six pairwise amino acid differences evolved under positive selection. We also identified a natural evolutionary intermediate that showed only three amino acid length differences in the PSD relative to its parental allele. This genotype showed an intermediate fitness because it implemented lethality regularly and induced femaleness infrequently (i.e., incomplete penetrance). We suggest incomplete penetrance as a mechanism through which new molecular switches can gradually and adaptively evolve.

The population genetics of sporophytic self-incompatibility in three hybridizing senecio (asteraceae) species with contrasting population histories

Hybridization generates evolutionary novelty and spreads adaptive variation. By promoting outcrossing, plant self-incompatibility (SI) systems also favor interspecific hybridization because the S locus is under strong negative frequency-dependent balancing selection. This study investigates the SI mating systems of three hybridizing Senecio species with contrasting population histories. Senecio aethnensis and S. chrysanthemifolius native to Sicily, form a hybrid zone at intermediate altitudes on Mount Etna, and their neo-homoploid hybrid species, S. squalidus, has colonized disturbed urban habitats in the UK during the last 150 years. We show that all three species express sporophytic SI (SSI), where pollen incompatibility is controlled by the diploid parental genome, and that SSI is inherited and functions normally in hybrids. Large-scale crossing studies of wild sampled populations allowed direct comparison of SSI between species and found that the main impacts of colonization in S. squalidus compared to Sicilian Senecio was a reduced number of S alleles, increased S allele frequencies, and increased interpopulation S allele sharing. In general, many S alleles were shared between species and the S locus showed reduced intra- and interspecific population genetic structure compared to molecular genetic markers, indicative of enhanced effective gene flow due to balancing selection.

Recent loss of self-incompatibility by degradation of the male component in allotetraploid Arabidopsis kamchatica

The evolutionary transition from outcrossing to self-fertilization (selfing) through the loss of self-incompatibility (SI) is one of the most prevalent events in flowering plants, and its genetic basis has been a major focus in evolutionary biology. In the Brassicaceae, the SI system consists of male and female specificity genes at the S-locus and of genes involved in the female downstream signaling pathway. During recent decades, much attention has been paid in particular to clarifying the genes responsible for the loss of SI. Here, we investigated the pattern of polymorphism and functionality of the female specificity gene, the S-locus receptor kinase (SRK), in allotetraploid Arabidopsis kamchatica. While its parental species, A. lyrata and A. halleri, are reported to be diploid and mainly self-incompatible, A. kamchatica is self-compatible. We identified five highly diverged SRK haplogroups, found their disomic inheritance and, for the first time in a wild allotetraploid species, surveyed the geographic distribution of SRK at the two homeologous S-loci across the species range. We found intact full-length SRK sequences in many accessions. Through interspecific crosses with the self-incompatible and diploid congener A. halleri, we found that the female components of the SI system, including SRK and the female downstream signaling pathway, are still functional in these accessions. Given the tight linkage and very rare recombination of the male and female components on the S-locus, this result suggests that the degradation of male components was responsible for the loss of SI in A. kamchatica. Recent extensive studies in multiple Brassicaceae species demonstrate that the loss of SI is often derived from mutations in the male component in wild populations, in contrast to cultivated populations. This is consistent with theoretical predictions that mutations disabling male specificity are expected to be more strongly selected than mutations disabling female specificity, or the female downstream signaling pathway.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this genomic region.

Investigating the effects of topography and clonality on genetic structuring within a large Norwegian population of Arabidopsis lyrata

BACKGROUND AND AIMS

The gene flow through pollen or seeds governs the extent of spatial genetic structure in plant populations. Another factor that can contribute to this pattern is clonal growth. The perennial species Arabidopsis lyrata ssp. petraea (Brassicaceae) is a self-incompatible, clonal species found in disjunctive populations in central and northern Europe.

METHODS

Fourteen microsatellite markers were employed to study the level of kinship and clonality in a high-altitude mountain valley at Spiterstulen, Norway. The population has a continuous distribution along the banks of the River Visa for about 1.5 km. A total of 17 (10 m x 10 m) squares were laid out in a north-south transect following the river on both sides.

KEY RESULTS

It is shown that clonal growth is far more common than previously shown in this species, although the overall size of the genets is small (mean diameter = 6.4 cm). Across the whole population there is no indication of isolation by distance, and spatial genetic structure is only visible on fine spatial scales. In addition, no effect of the river on the spatial distribution of genotypes was found.

CONCLUSIONS

Unexpectedly, the data show that populations of small perennials like A. lyrata can behave like panmictic units across relatively large areas at local sites, as opposed to earlier findings in central Europe.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.

The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier

BACKGROUND

The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.

RESULTS

Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.

CONCLUSIONS

Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.