Population genetic structure of Arabidopsis lyrata in Europe

Cloning and Selfing Affect Population Genetic Variation in Simulations of Outcrossing, Sexual Sea Stars

AbstractMany sea stars are well known for facultative or obligate asexual reproduction in both the adult and larval life-cycle stages. Some species and lineages are also capable of facultative or obligate hermaphroditic reproduction with self-fertilization. However, models of population genetic variation and empirical analyses of genetic data typically assume only sexual reproduction and outcrossing. A recent reanalysis of previously published empirical data (microsatellite genotypes) from two studies of one of the most well-known sea star species (the crown-of-thorns sea star; Acanthaster sp.) concluded that cloning and self-fertilization in that species are rare and contribute little to patterns of population genetic variation. Here we reconsider that conclusion by simulating the contribution of cloning and selfing to genetic variation in a series of models of sea star demography. Simulated variation in two simple models (analogous to previous analyses of empirical data) was consistent with high rates of cloning or selfing or both. More realistic scenarios that characterize population flux in sea stars of ecological significance, including outbreaks of crown-of-thorns sea stars that devastate coral reefs, invasions by Asterias amurensis, and epizootics of sea star wasting disease that kill Pisaster ochraceus, also showed significant but smaller effects of cloning and selfing on variation within subpopulations and differentiation between subpopulations. Future models or analyses of genetic variation in similar study systems might benefit from simulation modeling to characterize possible contributions of cloning or selfing to genetic variation in population samples or to understand the limits on inferring the effects of cloning or selfing in nature.

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.

Host preference and invasiveness of commensal bacteria in the Lotus and Arabidopsis root microbiota

Roots of different plant species are colonized by bacterial communities, that are distinct even when hosts share the same habitat. It remains unclear to what extent the host actively selects these communities and whether commensals are adapted to a specific plant species. To address this question, we assembled a sequence-indexed bacterial culture collection from roots and nodules of Lotus japonicus that contains representatives of most species previously identified using metagenomics. We analysed taxonomically paired synthetic communities from L. japonicus and Arabidopsis thaliana in a multi-species gnotobiotic system and detected signatures of host preference among commensal bacteria in a community context, but not in mono-associations. Sequential inoculation experiments revealed priority effects during root microbiota assembly, where established communities are resilient to invasion by latecomers, and that host preference of commensal bacteria confers a competitive advantage in their cognate host. Our findings show that host preference in commensal bacteria from diverse taxonomic groups is associated with their invasiveness into standing root-associated communities.

Interspecific introgression mediates adaptation to whole genome duplication

Adaptive gene flow is a consequential phenomenon across all kingdoms. Although recognition is increasing, there is no study showing that bidirectional gene flow mediates adaptation at loci that manage core processes. We previously discovered concerted molecular changes among interacting members of the meiotic machinery controlling crossover number upon adaptation to whole-genome duplication (WGD) in Arabidopsis arenosa. Here we conduct a population genomic study to test the hypothesis that adaptation to WGD has been mediated by adaptive gene flow between A. arenosa and A. lyrata. We find that A. lyrata underwent WGD more recently than A. arenosa, suggesting that pre-adapted alleles have rescued nascent A. lyrata, but we also detect gene flow in the opposite direction at functionally interacting loci under the most extreme levels of selection. These data indicate that bidirectional gene flow allowed for survival after WGD, and that the merger of these species is greater than the sum of their parts.

Whole genome duplication (WGD) presents new challenges to the establishment of optimal allelic combinations and to the meiotic machinery. Here, the authors show that adaptive gene flow from Arabidopsis arenosa could rescue the nascent A. lyrata from extinction following WGD.

Phenotypic Trait Variation as a Response to Altitude-Related Constraints in Arabidopsis Populations

Natural variations help in identifying genetic mechanisms of morphologically and developmentally complex traits. Mountainous habitats provide an altitudinal gradient where one species encounters different abiotic conditions. We report the study of 341 individuals of Arabidopsis thaliana derived from 30 natural populations not belonging to the 1001 genomes, collected at increasing altitudes, between 200 and 1800 m in the Pyrenees. Class III peroxidases and ribosomal RNA sequences were used as markers to determine the putative genetic relationships among these populations along their altitudinal gradient. Using Bayesian-based statistics and phylogenetic analyses, these Pyrenean populations appear with significant divergence from the other regional accessions from 1001 genome (i.e., from north Spain or south France). Individuals of these populations exhibited varying phenotypic changes, when grown at sub-optimal temperature (22 vs. 15°C). These phenotypic variations under controlled conditions reflected intraspecific morphological variations. This study could bring new information regarding the west European population structure of A. thaliana and its phenotypic variations at different temperatures. The integrative analysis combining genetic, phenotypic variation and environmental datasets is used to analyze the acclimation of population in response to temperature changes. Regarding their geographical proximity and environmental diversity, these populations represent a tool of choice for studying plant response to temperature variation.

HIGHLIGHTS

-Studying the natural diversity of A. thaliana in the Pyrenees mountains helps to understand European population structure and to evaluate the phenotypic trait variation in response to climate change.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

Natural variation in tolerance to sub-zero temperatures among populations of Arabidopsis lyrata ssp. petraea

BACKGROUND

Temperature is one of the most important abiotic factors limiting plant growth and productivity. Many plants exhibit cold acclimation to prepare for the likelihood of freezing as temperatures decrease towards 0 °C. The physiological mechanisms associated with enabling increased tolerance to sub-zero temperatures vary between species and genotypes. Geographically and climatically diverse populations of Arabidopsis lyrata ssp. petraea were examined for their ability to survive, maintain functional photosynthetic parameters and cellular electrolyte leakage integrity after being exposed to sub-zero temperatures. The duration of cold acclimation prior to sub-zero temperatures was also manipulated (2 and 14 days).

RESULTS

We found that there was significant natural variation in tolerances to sub-zero temperatures among populations of A. petraea. The origin of the population affected the acclimation response and survival after exposure to sub-zero temperatures. Cold acclimation of plants prior to sub-zero temperatures affected the maximum quantum efficiency of photosystem II (PSII) (Fv/Fm) in that plants that were cold acclimated for longer periods had higher values of Fv/Fm as a result of sub-zero temperatures. The inner immature leaves were better able to recover Fv/Fm from sub-zero temperatures than mature outer leaves. The Irish population (Leitrim) acclimated faster, in terms of survival and electrolyte leakage than the Norwegian population (Helin).

CONCLUSION

The ability to survive, recover photosynthetic processes and cellular electrolyte leakage after exposure to sub-zero temperatures is highly dependent on the duration of cold acclimation.

How much genetic variation is stored in the endangered and fragmented shrub Tetraena mongolica Maxim?

Tetraena mongolica Maxim (Zygophyllaceae) is an endangered species endemic to western Inner Mongolia and China, and is currently threatened by habitat loss and human over-exploitation. We explored the genetic background, its genetic diversity, population structure, and demographic history, based on 12 polymorphic nuclear microsatellite loci. Our results indicated high genetic diversity in extant populations, but no distinguishable gene cluster corresponding with a specific biogeography. Population demography analysis using a MSVAR indicated a strong, recent population decline approximately 5,455 years ago. These results suggest that the Yellow River and Zhuozi Mountain range may not prevent pollination between populations. Finally, we surmised that the population demography of T. mongolica was likely to have been affected by early mankind activities.

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4806-7) contains supplementary material, which is available to authorized users.

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.

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.

R-gene variation across Arabidopsis lyrata subspecies: effects of population structure, selection and mating system

Background

Examining allelic variation of R-genes in closely related perennial species of Arabidopsis thaliana is critical to understanding how population structure and ecology interact with selection to shape the evolution of innate immunity in plants. We finely sampled natural populations of Arabidopsis lyrata from the Great Lakes region of North America (A. l. lyrata) and broadly sampled six European countries (A. l. petraea) to investigate allelic variation of two R-genes (RPM1 and WRR4) and neutral genetic markers (Restriction Associated DNA sequences and microsatellites) in relation to mating system, phylogeographic structure and subspecies divergence.

Results

Fine-scale sampling of populations revealed strong effects of mating system and population structure on patterns of polymorphism for both neutral loci and R-genes, with no strong evidence for selection. Broad geographic sampling revealed evidence of balancing selection maintaining polymorphism in R-genes, with elevated heterozygosity and diversity compared to neutral expectations and sharing of alleles among diverged subspecies. Codon-based tests detected both positive and purifying selection for both R-genes, as commonly found for animal immune genes.

Conclusions

Our results highlight that combining fine and broad-scale sampling strategies can reveal the multiple factors influencing polymorphism and divergence at potentially adaptive genes such as R-genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0665-5) contains supplementary material, which is available to authorized users.

Phylogeography of Arabidopsis halleri (Brassicaceae) in mountain regions of Central Europe inferred from cpDNA variation and ecological niche modelling

The present study aimed to investigate phylogeographical patterns present within A. halleri in Central Europe. 1,281 accessions sampled from 52 populations within the investigated area were used in the study of genetic variation based on chloroplast DNA. Over 500 high-quality species occurrence records were used in ecological niche modelling experiments. We evidenced the presence of a clear phylogeographic structure within A. halleri in Central Europe. Our results showed that two genetically different groups of populations are present in western and eastern part of the Carpathians. The hypothesis of the existence of a glacial refugium in the Western Carpathians adn the Bohemian Forest cannot be rejected from our data. It seems, however, that the evidence collected during the present study is not conclusive. The area of Sudetes was colonised after LGM probably by migrants from the Bohemian Forest.

Timing of shoot development transitions affects degree of perenniality in Arabidopsis lyrata (Brassicaceae)

BACKGROUND

Perenniality is best understood in quantitative terms, involving the relationship between production vs. turnover of meristems, biomass, or energy reserves. Previous quantitative trait locus (QTL) studies using divergent populations of the perennial rock cress Arabidopsis lyrata have shown that trade-offs in vegetative growth vs. reproduction are due to cascading effects of differences in early vegetative development, which contribute to local adaptation. However, details of the developmental differences and how they affect perenniality remained unclear. In this study, we investigated in detail the developmental differences in perenniality between populations. A. lyrata from Norway and North Carolina populations, representing contrasting environments and degrees of perenniality, were grown under controlled conditions, and data were collected on plant phenology and shoot-level development. We tested hypotheses that differences in perenniality involve strict allocation of lateral meristems to vegetative vs. reproductive fates, or alternatively quantitative effects of pre-reproductive vegetative development.

RESULTS

The two populations showed large differences in the degree of vegetative development on individual shoots prior to reproductive transitions. The number of leaves produced on shoots prior to bolting, and not strict meristem allocation or variation in apical dominance, was able to explain variation in the number of inflorescences on individual plants. These results suggested that allocation of time to shoot vegetative vs. reproductive development could be a major factor in resource allocation differences between the populations.

CONCLUSIONS

Based on these results and those of previous QTL studies, we propose a model in which the degree of shoot vegetative development shapes the developmental context for reproduction and subsequent vegetative growth in different environments. Climate-specific effects of shoot development patterns on reproductive output and survival may result in divergent evolutionary trajectories along a perenniality continuum, which may have broader relevance for plant life history evolution.

Taming the wild: resolving the gene pools of non-model Arabidopsis lineages

BACKGROUND

Wild relatives in the genus Arabidopsis are recognized as useful model systems to study traits and evolutionary processes in outcrossing species, which are often difficult or even impossible to investigate in the selfing and annual Arabidopsis thaliana. However, Arabidopsis as a genus is littered with sub-species and ecotypes which make realizing the potential of these non-model Arabidopsis lineages problematic. There are relatively few evolutionary studies which comprehensively characterize the gene pools across all of the Arabidopsis supra-groups and hypothesized evolutionary lineages and none include sampling at a world-wide scale. Here we explore the gene pools of these various taxa using various molecular markers and cytological analyses.

RESULTS

Based on ITS, microsatellite, chloroplast and nuclear DNA content data we demonstrate the presence of three major evolutionary groups broadly characterized as A. lyrata group, A. halleri group and A. arenosa group. All are composed of further species and sub-species forming larger aggregates. Depending on the resolution of the marker, a few closely related taxa such as A. pedemontana, A. cebennensis and A. croatica are also clearly distinct evolutionary lineages. ITS sequences and a population-based screen based on microsatellites were highly concordant. The major gene pools identified by ITS sequences were also significantly differentiated by their homoploid nuclear DNA content estimated by flow cytometry. The chloroplast genome provided less resolution than the nuclear data, and it remains unclear whether the extensive haplotype sharing apparent between taxa results from gene flow or incomplete lineage sorting in this relatively young group of species with Pleistocene origins.

CONCLUSIONS

Our study provides a comprehensive overview of the genetic variation within and among the various taxa of the genus Arabidopsis. The resolved gene pools and evolutionary lineages will set the framework for future comparative studies on genetic diversity. Extensive population-based phylogeographic studies will also be required, however, in particular for A. arenosa and their affiliated taxa and cytotypes.

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.

Genetic changes in flowering and morphology in response to adaptation to a high-latitude environment in Arabidopsis lyrata

BACKGROUND AND AIMS

The adaptive plastic reactions of plant populations to changing climatic factors, such as winter temperatures and photoperiod, have changed during range shifts after the last glaciation. Timing of flowering is an adaptive trait regulated by environmental cues. Its genetics has been intensively studied in annual plants, but in perennials it is currently not well characterized. This study examined the genetic basis of differentiation in flowering time, morphology, and their plastic responses to vernalization in two locally adapted populations of the perennial Arabidopsis lyrata: (1) to determine whether the two populations differ in their vernalization responses for flowering phenology and morphology; and (2) to determine the genomic areas governing differentiation and vernalization responses.

METHODS

Two A. lyrata populations, from central Europe and Scandinavia, were grown in growth-chamber conditions with and without cold treatment. A QTL analysis was performed to find genomic regions that interact with vernalization.

KEY RESULTS

The population from central Europe flowered more rapidly and invested more in inflorescence growth than the population from alpine Scandinavia, especially after vernalization. The alpine population had consistently a low number of inflorescences and few flowers, suggesting strong constraints due to a short growing season, but instead had longer leaves and higher leaf rosettes. QTL mapping in the F2 population revealed genomic regions governing differentiation in flowering time and morphology and, in some cases, the allelic effects from the two populations on a trait were influenced by vernalization (QTL × vernalization interactions).

CONCLUSIONS

The results indicate that many potentially adaptive genetic changes have occurred during colonization; the two populations have diverged in their plastic responses to vernalization in traits closely connected to fitness through changes in many genomic areas.

Adaptation at range margins: common garden trials and the performance of Arabidopsis lyrata across its northwestern European range

Widely distributed species, such as the perennial plant Arabidopsis lyrata, face a range of environmental conditions across space, creating selective pressures for local evolutionary adaptation. The species' fragmented distribution may reduce gene flow, which could either reduce or increase adaptive potential. The substantial variation in phenotypic traits observed across this species' northwestern European range may reflect a combination of plastic responses to environmental conditions, evolutionary adaptation and nonadaptive genetic differentiation. We conducted multi-site common garden experiments to study differences in plant performance in core and marginal areas. Plants from eight source populations representing the species' full geographic and altitudinal range in northwestern Europe were planted out in Iceland, Sweden, Scotland and Wales. We found evidence of both strong plastic responses and apparently adaptive differentiation in performance. Most evidence for local adaptation was found at range margins, with the strongest effects on reproductive output. Both biotic and abiotic factors affected performance, especially at range margins. Performance of most plants was best in the Scottish and Swedish common garden sites, in the core of the species' distribution. Despite adaptations at range margins, the performance of the species declines at distributional limits, with extreme southern populations looking particularly vulnerable.

Role of vernalization and of duplicated FLOWERING LOCUS C in the perennial Arabidopsis lyrata

FLOWERING LOCUS C (FLC) is one of the main genes influencing the vernalization requirement and natural flowering time variation in the annual Arabidopsis thaliana. Here we studied the effects of vernalization on flowering and its genetic basis in the perennial Arabidopsis lyrata. Two tandemly duplicated FLC genes (FLC1 and FLC2) were compared with respect to expression and DNA sequence. The effect of vernalization on flowering and on the expression of FLC1 was studied in three European populations. The genetic basis of the FLC1 expression difference between two of the populations was further studied by expression quantitative trait locus (eQTL) mapping and sequence analysis. FLC1 was shown to have a likely role in the vernalization requirement for flowering in A. lyrata. Vernalization decreased its expression and the northern study populations showed higher FLC1 expression than the southern one. eQTL mapping between two of the populations revealed one eQTL affecting FLC1 expression in the genomic region containing the FLC genes. Most FLC1 sequence differences between the study populations were found in the promoter region and in the first intron. Variation in the FLC1 sequence may cause differences in FLC1 expression between late- and early-flowering A. lyrata populations.

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

PREMISE OF THE STUDY

Plant populations that face new environments adapt and diverge simultaneously, and both processes leave footprints in their genetic diversity. Arabidopsis lyrata is an excellent species for studying these processes. Pairs of populations and subspecies of A. lyrata represent different stages of divergence. These populations are also known to be locally adapted and display various stages of emerging reproductive isolation.

METHODS

We used nucleotide diversity data from 19 loci to estimate divergence times and levels of diversity among nine A. lyrata populations. Traditional distance-based methods and model-based clustering analysis were used to supplement pairwise coalescence-based analysis of divergence.

KEY RESULTS

Estimated divergence times varied from 130,000 generations between North American and European subspecies to 39,000 generations between central European and Scandinavian populations. In concordance with previous studies, the highest level of diversity was found in Central Europe and the lowest in North America and a diverged Russian Karhumäki population. Local adaptation among Northern and central European populations has emerged during the last 39,000 generations. Populations that are reproductively isolated by prezygotic mechanisms have been separated for a longer time period of ∼70,000 generations but still have shared nucleotide polymorphism.

CONCLUSIONS

In A. lyrata, reproductively isolated populations started to diverge ∼70,000 generations ago and more closely related, locally adapted populations have been separate lineages for ∼39,000 generations. However, based on the posterior distribution of divergence times, the processes leading to reproductive isolation and local adaptation are likely to temporally coincide.

Nuclear and chloroplast DNA phylogeography reveals vicariance among European populations of the model species for the study of metal tolerance, Arabidopsis halleri (Brassicaceae)

Arabidopsis halleri is a pseudometallophyte involved in numerous molecular studies of the adaptation to anthropogenic metal stress. In order to test the representativeness of genetic accessions commonly used in these studies, we investigated the A. halleri population genetic structure in Europe. Microsatellite and nucleotide polymorphisms from the nuclear and chloroplast genomes, respectively, were used to genotype 65 populations scattered over Europe. The large-scale population structure was characterized by a significant phylogeographic signal between two major genetic units. The localization of the phylogeographic break was assumed to result from vicariance between large populations isolated in southern and central Europe, on either side of ice sheets covering the Alps during the Quaternary ice ages. Genetic isolation was shown to be maintained in western Europe by the high summits of the Alps, whereas admixture was detected in the Carpathians. Considering the phylogeographic literature, our results suggest a distinct phylogeographic pattern for European species occurring in both mountain and lowland habitats. Considering the evolution of metal adaptation in A. halleri, it appears that recent adaptations to anthropogenic metal stress that have occurred within either phylogeographic unit should be regarded as independent events that potentially have involved the evolution of a variety of genetic mechanisms.

Does speciation between Arabidopsis halleri and Arabidopsis lyrata coincide with major changes in a molecular target of adaptation?

Ever since Darwin proposed natural selection as the driving force for the origin of species, the role of adaptive processes in speciation has remained controversial. In particular, a largely unsolved issue is whether key divergent ecological adaptations are associated with speciation events or evolve secondarily within sister species after the split. The plant Arabidopsis halleri is one of the few species able to colonize soils highly enriched in zinc and cadmium. Recent advances in the molecular genetics of adaptation show that the physiology of this derived ecological trait involves copy number expansions of the AhHMA4 gene, for which orthologs are found in single copy in the closely related A. lyrata and the outgroup A. thaliana. To gain insight into the speciation process, we ask whether adaptive molecular changes at this candidate gene were contemporary with important stages of the speciation process. We first inferred the scenario and timescale of speciation by comparing patterns of variation across the genomic backgrounds of A. halleri and A. lyrata. Then, we estimated the timing of the first duplication of AhHMA4 in A. halleri. Our analysis suggests that the historical split between the two species closely coincides with major changes in this molecular target of adaptation in the A. halleri lineage. These results clearly indicate that these changes evolved in A. halleri well before industrial activities fostered the spread of Zn- and Cd-polluted areas, and suggest that adaptive processes related to heavy-metal homeostasis played a major role in the speciation process.

The relative importance of factors determining genetic drift: mating system, spatial genetic structure, habitat and census size in Arabidopsis lyrata

• The mating system, dispersal and census size are predicted to determine the magnitude of genetic drift, but little is known about their relative importance in nature. • We estimated the contributions of several population-level features to genetic drift in 18 populations of Arabidopsis lyrata. The factors were outcrossing rate, within-population spatial genetic structure, census size and substrate type. The expected heterozygosity (H(E)) at 10 microsatellite loci was taken to reflect the effective population size (N(e)) and the strength of genetic drift. • The mating system explained most of the variation in H(E) (60%), followed by substrate (10%), genetic structure (9%) and census size (6%). The most outcrossing population had a +0.32 higher predicted H(E) than the most selfing population; the estimated N(e) of selfing populations was less than half that of outcrossing populations. Rocky outcrops supported populations with a +0.14 higher H(E) than did sandy substrates. The most structured population had a +0.24 higher H(E) than the least structured population, and the largest population had a +0.18 higher H(E) than the smallest population. • This study illustrates the importance of outcrossing, genetic structure and the physical environment--together with census size--in maintaining H(E), and suggests that multiple population-level characteristics influence N(e) and the action of genetic drift.

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.

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.

Influence of mutation rate on estimators of genetic differentiation--lessons from Arabidopsis thaliana

BACKGROUND

The analysis of molecular variation within and between populations is crucial to establish strategies for conservation as well as to detect the footprint of spatially heterogeneous selection. The traditional estimator of genetic differentiation (F(ST)) has been shown to be misleading if genetic diversity is high. Alternative estimators of F(ST) have been proposed, but their robustness to variation in mutation rate is not clearly established. We first investigated the effect of mutation and migration rate using computer simulations and examined their joint influence on Q(ST), a measure of genetic differentiation for quantitative traits. We further used experimental data in natural populations of Arabidopsis thaliana to characterize the effect of mutation rate on various estimates of population differentiation. Since natural species exhibit various degrees of self-fertilisation, we also investigated the effect of mating system on the different estimators.

RESULTS

If mutation rate is high and migration rate low, classical measures of genetic differentiation are misleading. Only Phi(ST), an estimator that takes the mutational distances between alleles into account, is independent of mutation rate, for all migration rates. However, the performance of Phi(ST) depends on the underlying mutation model and departures from this model cause its performance to degrade. We further show that Q(ST) has the same bias. We provide evidence that, in A. thaliana, microsatellite variation correlates with mutation rate. We thereby demonstrate that our results on estimators of genetic differentiation have important implications, even for species that are well established models in population genetics and molecular biology.

CONCLUSIONS

We find that alternative measures of differentiation like F'(ST) and D are not suitable for estimating effective migration rate and should not be used in studies of local adaptation. Genetic differentiation should instead be measured using an estimator that takes mutation rate into account, such as Phi(ST). Furthermore, in systems where migration between populations is low, such as A. thaliana, Q(ST) < F(ST) cannot be taken as evidence for homogenising selection as has been traditionally thought.

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.

Population structure and historical biogeography of European Arabidopsis lyrata

Understanding the natural history of model organisms is important for the effective use of their genomic resources. Arabidopsis lyrata has emerged as a useful plant for studying ecological and evolutionary genetics, based on its extensive natural variation, sequenced genome and close relationship to A. thaliana. We studied genetic diversity across the entire range of European Arabidopsis lyrata ssp. petraea, in order to explore how population history has influenced population structure. We sampled multiple populations from each region, using nuclear and chloroplast genome markers, and combined population genetic and phylogeographic approaches. Within-population diversity is substantial for nuclear allozyme markers (mean P=0.610, A(e)=1.580, H(e)=0.277) and significantly partitioned among populations (F(ST)=0.271). The Northern populations have modestly increased inbreeding (F(IS)=0.163 verses F(IS)=0.093), but retain comparable diversity to central European populations. Bottlenecks are common among central and northern Europe populations, indicating recent demographic history as a dominant factor in structuring the European diversity. Although the genetic structure was detected at all geographic scales, two clear differentiated units covering northern and central European areas (F(CT) =0.155) were identified by Bayesian analysis and supported by regional pairwise F(CT) calculations. A highly similar geographic pattern was observed from the distribution of chloroplast haplotypes, with the dominant northern haplotypes absent from central Europe. We conclude A. l. petraea's cold-tolerance and preference for disturbed habitats enabled glacial survival between the alpine and Nordic glaciers in central Europe and an additional cryptic refugium. While German populations are probable peri-glacial leftovers, Eastern Austrian populations have diversity patterns possibly compatible with longer-term survival.

The evolution of transposable elements in natural populations of self-fertilizing Arabidopsis thaliana and its outcrossing relative Arabidopsis lyrata

Background

Transposable Elements (TEs) make up the majority of plant genomes, and thus understanding TE evolutionary dynamics is key to understanding plant genome evolution. Plant reproductive systems are diverse and mating type variation is one factor among many hypothesized to influence TE evolutionary dynamics. Here, we collected a large TE-display data set in self-fertilizing Arabidopsis thaliana, and compared it to data gathered in outcrossing Arabidopsis lyrata. We analyzed seven TE families in four natural populations of each species to tease apart the effects of mating system, demography, transposition, and selection in determining patterns of TE diversity.

Results

Measures of TE band differentiation were largely consistent across TE families. However, patterns of diversity in A. thaliana Ac elements differed significantly from that other TEs, perhaps signaling a lack of recent transposition. Across TE families, we estimated higher allele frequencies and lower selection coefficients on A. thaliana TE insertions relative to A. lyrata TE insertions.

Conclusions

The differences in TE distributions between the two Arabidopsis species represents a synthesis of evolutionary forces that include the transposition dynamics of individual TE families and the demographic histories of populations. There are also species-specific differences that could be attributed to the effects of mating system, including higher overall allele frequencies in the selfing lineage and a greater proportion of among population TE diversity in the outcrossing lineage.

Ecogeographic genetic epidemiology

Complex diseases such as cancer and heart disease result from interactions between an individual's genetics and environment, i.e. their human ecology. Rates of complex diseases have consistently demonstrated geographic patterns of incidence, or spatial "clusters" of increased incidence relative to the general population. Likewise, genetic subpopulations and environmental influences are not evenly distributed across space. Merging appropriate methods from genetic epidemiology, ecology and geography will provide a more complete understanding of the spatial interactions between genetics and environment that result in spatial patterning of disease rates. Geographic information systems (GIS), which are tools designed specifically for dealing with geographic data and performing spatial analyses to determine their relationship, are key to this kind of data integration. Here the authors introduce a new interdisciplinary paradigm, ecogeographic genetic epidemiology, which uses GIS and spatial statistical analyses to layer genetic subpopulation and environmental data with disease rates and thereby discern the complex gene-environment interactions which result in spatial patterns of incidence.

Floral and insect-induced volatile formation in Arabidopsis lyrata ssp. petraea, a perennial, outcrossing relative of A. thaliana

Volatile organic compounds have been reported to serve some important roles in plant communication with other organisms, but little is known about the biological functions of most of these substances. To gain insight into this problem, we have compared differences in floral and vegetative volatiles between two closely related plant species with different life histories. The self-pollinating annual, Arabidopsis thaliana, and its relative, the outcrossing perennial, Arabidopsis lyrata, have markedly divergent life cycles and breeding systems. We show that these differences are in part reflected in the formation of distinct volatile mixtures in flowers and foliage. Volatiles emitted from flowers of a German A. lyrata ssp. petraea population are dominated by benzenoid compounds in contrast to the previously described sesquiterpene-dominated emissions of A. thaliana flowers. Flowers of A. lyrata ssp. petraea release benzenoid volatiles in a diurnal rhythm with highest emission rates at midday coinciding with observed visitations of pollinating insects. Insect feeding on leaves of A. lyrata ssp. petraea causes a variable release of the volatiles methyl salicylate, C11- and C16-homoterpenes, nerolidol, plus the sesquiterpene (E)-beta-caryophyllene, which in A. thaliana is emitted exclusively from flowers. An insect-induced gene (AlCarS) with high sequence similarity to the florally expressed (E)-beta-caryophyllene synthase (AtTPS21) from A. thaliana was identified from individuals of a German A. lyrata ssp. petraea population. Recombinant AlCarS converts the sesquiterpene precursor, farnesyl diphosphate, into (E)-beta-caryophyllene with alpha-humulene and alpha-copaene as minor products indicating its close functional relationship to the A. thaliana AtTPS21. Differential regulation of these genes in flowers and foliage is consistent with the different functions of volatiles in the two Arabidopsis species.

RPW8 and resistance to powdery mildew pathogens in natural populations of Arabidopsis lyrata

It is not clear to what extent the orthologues of genes that are adaptively important in one species also contribute to adaptive variation in others. Here, we examine Arabidopsis lyrata to assess the functional and evolutionary significance of natural variation in an orthologue of the gene RPW8 known to be a major determinant of powdery mildew resistance in Arabidopsis thaliana. We assessed the sequence variation at RPW8 and the associated resistance reaction in populations of A. lyrata ssp. petraea. Neutrality tests were performed to understand the importance of local adaptation in maintaining variation at the locus. Highly truncated RPW8 proteins were frequent in all populations and were associated with an increased risk of susceptibility. Haplotypes encoding full-length proteins were highly significantly associated with resistance. There were no signatures of selection at the species-wide level, but some evidence for positive selection in two populations. RPW8 in A. lyrata appears to have a role in powdery mildew resistance, similar to its orthologue in A. thaliana. Unlike A. thaliana, A. lyrata contains a genetic component that can act independently of RPW8 to confer resistance to powdery mildew pathogens. Infrequent local selective sweeps may favour different alleles in different populations, and thereby contribute to the maintenance of species-wide variation at the locus.

Antagonism between local dispersal and self-incompatibility systems in a continuous plant population

Many self-incompatible plant species exist in continuous populations in which individuals disperse locally. Local dispersal of pollen and seeds facilitates inbreeding because pollen pools are likely to contain relatives. Self-incompatibility promotes outbreeding because relatives are likely to carry incompatible alleles. Therefore, populations can experience an antagonism between these forces. In this study, a novel computational model is used to explore the effects of this antagonism on gene flow, allelic diversity, neighbourhood sizes, and identity by descent. I confirm that this antagonism is sensitive to dispersal levels and linkage. However, the results suggest that there is little to no difference between the effects of gametophytic and sporophytic self-incompatibility systems (GSI and SSI) on unlinked loci. More importantly, both GSI and SSI affect unlinked loci in a manner similar to obligate outcrossing without mating types. This suggests that the primary evolutionary impact of self-incompatibility systems may be to prevent selfing, and prevention of biparental inbreeding might be a beneficial side-effect.

Variation at range margins across multiple spatial scales: environmental temperature, population genetics and metabolomic phenotype

Range margins are spatially complex, with environmental, genetic and phenotypic variations occurring across a range of spatial scales. We examine variation in temperature, genes and metabolomic profiles within and between populations of the subalpine perennial plant Arabidopsis lyrata ssp. petraea from across its northwest European range. Our surveys cover a gradient of fragmentation from largely continuous populations in Iceland, through more fragmented Scandinavian populations, to increasingly widely scattered populations at the range margin in Scotland, Wales and Ireland. Temperature regimes vary substantially within some populations, but within-population variation represents a larger fraction of genetic and especially metabolomic variances. Both physical distance and temperature differences between sites are found to be associated with genetic profiles, but not metabolomic profiles, and no relationship was found between genetic and metabolomic population structures in any region. Genetic similarity between plants within populations is the highest in the fragmented populations at the range margin, but differentiation across space is the highest there as well, suggesting that regional patterns of genetic diversity may be scale dependent.

Population genetic structure of two Medicago species shaped by distinct life form, mating system and seed dispersal

BACKGROUND AND AIMS

Life form, mating system and seed dispersal are important adaptive traits of plants. In the first effort to characterize in detail the population genetic structure and dynamics of wild Medicago species in China, a population genetic study of two closely related Medicago species, M. lupulina and M. ruthenica, that are distinct in these traits, are reported. These species are valuable germplasm resources for the improvement of Medicago forage crops but are under threat of habitat destruction.

METHODS

Three hundred and twenty-eight individuals from 16 populations of the annual species, M. lupulina, and 447 individuals from 15 populations of the perennial species, M. ruthenica, were studied using 15 and 17 microsatellite loci, respectively. Conventional and Bayesian-clustering analyses were utilized to estimate population genetic structure, mating system and gene flow.

KEY RESULTS

Genetic diversity of M. lupulina (mean H(E)=0.246) was lower than that of M. ruthenica (mean H(E)=0.677). Populations of M. lupulina were more highly differentiated (F(ST)=0.535) than those of M. ruthenica (F(ST)=0.130). For M. lupulina, 55.5% of the genetic variation was partitioned among populations, whereas 76.6% of the variation existed within populations of M. ruthenica. Based on the genetic data, the selfing rates of M. lupulina and M. ruthenica were estimated at 95.8% and 29.5%, respectively. The genetic differentiation among populations of both species was positively correlated with geographical distance.

CONCLUSIONS

The mating system differentiation estimated from the genetic data is consistent with floral morphology and observed pollinator visitation. There was a much higher historical gene flow in M. ruthenica than in M. lupulina, despite more effective seed dispersal mechanisms in M. lupulina. The population genetic structure and geographical distribution of the two Medicago species have been shaped by life form, mating systems and seed dispersal mechanisms.

Multilocus patterns of nucleotide diversity, population structure and linkage disequilibrium in Boechera stricta, a wild relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in approximately 10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have approximately 45% higher recombination per kilobase than outcrossing A. lyrata.

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.

Demography and weak selection drive patterns of transposable element diversity in natural populations of Arabidopsis lyrata

Transposable elements (TEs) are the major component of most plant genomes, and characterizing their population dynamics is key to understanding plant genome complexity. Yet there have been few studies of TE population genetics in plant systems. To study the roles of selection, transposition, and demography in shaping TE population diversity, we generated a polymorphism dataset for six TE families in four populations of the flowering plant Arabidopsis lyrata. The TE data indicated significant differentiation among populations, and maximum likelihood procedures suggested weak selection. For strongly bottlenecked populations, the observed TE band-frequency spectra fit data simulated under neutral demographic models constructed from nucleotide polymorphism data. Overall, we propose that TEs are subjected to weak selection, the efficacy of which varies as a function of demographic factors. Thus, demographic effects could be a major factor driving distributions of TEs among plant lineages.

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.

High DNA sequence diversity in pericentromeric genes of the plant Arabidopsis lyrata

Differences in neutral diversity at different loci are predicted to arise due to differences in mutation rates and from the "hitchhiking" effects of natural selection. Consistent with hitchhiking models, Drosophila melanogaster chromosome regions with very low recombination have unusually low nucleotide diversity. We compared levels of diversity from five pericentromeric regions with regions of normal recombination in Arabidopsis lyrata, an outcrossing close relative of the highly selfing A. thaliana. In contrast with the accepted theoretical prediction, and the pattern in Drosophila, we found generally high diversity in pericentromeric genes, which is consistent with the observation in A. thaliana. Our data rule out balancing selection in the pericentromeric regions, suggesting that hitchhiking is more strongly reducing diversity in the chromosome arms than the pericentromere regions.

Demographic processes shaping genetic variation

Demographic processes modulate genome-wide levels and patterns of genetic variation via impacting effective population size independently of natural selection. Such processes include the perturbation of population distributions from external events shaping habitat landscape and internal factors shaping the probability of contemporaneous alleles in a population (coalescence). Several patterns have recently emerged: spatial and temporal heterogeneity in population structure have different influences on the persistence of new mutations and genetic variation, multi-locus analyses indicate that gene flow continues to occur during speciation and the incorporation of demographic processes into models of molecular evolution and association genetics approaches has improved statistical power to detect deviations from neutral-equilibrium expectations and decreased false positive rates.