Demographic history of european populations of Arabidopsis thaliana

Ecological trade-offs drive phenotypic and genetic differentiation of Arabidopsis thaliana in Europe

Plant diversity is shaped by trade-offs between traits related to competitive ability, propagule dispersal, and stress resistance. However, we still lack a clear understanding of how these trade-offs influence species distribution and population dynamics. In Arabidopsis thaliana, recent genetic analyses revealed a group of cosmopolitan genotypes that successfully recolonized Europe from its center after the last glaciation, excluding older (relict) lineages from the distribution except for their north and south margins. Here, we tested the hypothesis that cosmopolitans expanded due to higher colonization ability, while relicts persisted at the margins due to higher tolerance to competition and/or stress. We compared the phenotypic and genetic differentiation between 71 European genotypes originating from the center, and the south and north margins. We showed that a trade-off between plant fecundity and seed mass shapes the differentiation of A. thaliana in Europe, suggesting that the success of the cosmopolitan groups could be explained by their high dispersal ability. However, at both north and south margins, we found evidence of selection for alleles conferring low dispersal but highly competitive and stress-resistance abilities. This study sheds light on the role of ecological trade-offs as evolutionary drivers of the distribution and dynamics of plant populations.

The Roles of Drift and Selection on Short Stamen Loss in Arabidopsis thaliana along an Elevational Gradient in the Spanish Pyrenees

Traits that have lost function sometimes persist through evolutionary time. These traits may be maintained by a lack of standing genetic variation for the trait, if selection against the trait is weak relative to drift, or if they have a residual function. To determine the evolutionary processes shaping whether nonfunctional traits are retained or lost, we investigated short stamens in 16 populations of Arabidopsis thaliana along an elevational cline in the Spanish Pyrenees. We found a cline in short stamen number from retention of short stamens in high elevation populations to incomplete loss in low elevation populations. We did not find evidence that limited genetic variation constrains the loss of short stamens at high elevations nor evidence for divergent selection on short stamens between high and low elevations. Finally, we identified loci associated with short stamens in the Spanish Pyrenees that are different from loci associated with variation in short stamen number across latitudes from a previous study. Overall, we did not identify the evolutionary mechanisms maintaining an elevational cline in short stamen number but did identify different genetic loci underlying the variation in short stamen along similar phenotypic clines.

Teaser text

The evolutionary mechanisms underlying loss or retention of traits that have lost function are poorly understood. Short stamens in Arabidopsis thaliana provide a compelling system to investigate the roles of genetic drift and selection in trait loss across latitudinal and elevational clines. This study investigates the role of drift and selection in short stamen loss in 16 populations of A. thaliana along an elevational gradient in the Spanish Pyrenees. An investigation of the genetic loci underlying variation in short stamen number suggests mutations in different genes may cause trait loss in similar phenotypic clines within a species.

Genetic adaptations in the population history of Arabidopsis thaliana

A population encounters a variety of environmental stresses, so the full source of its resilience can only be captured by collecting all the signatures of adaptation to the selection of the local environment in its population history. Based on the multiomic data of Arabidopsis thaliana, we constructed a database of phenotypic adaptations (p-adaptations) and gene expression (e-adaptations) adaptations in the population. Through the enrichment analysis of the identified adaptations, we inferred a likely scenario of adaptation that is consistent with the biological evidence from experimental work. We analyzed the dynamics of the allele frequencies at the 23,880 QTLs of 174 traits and 8,618 eQTLs of 1,829 genes with respect to the total SNPs in the genomes and identified 650 p-adaptations and 3,925 e-adaptations [false discovery rate (FDR) = 0.05]. The population underwent large-scale p-adaptations and e-adaptations along 4 lineages. Extremely cold winters and short summers prolonged seed dormancy and expanded the root system architecture. Low temperatures prolonged the growing season, and low light intensity required the increased chloroplast activity. The subtropical and humid environment enhanced phytohormone signaling pathways in response to the biotic and abiotic stresses. Exposure to heavy metals selected alleles for lower heavy metal uptake from soil, lower growth rate, lower resistance to bacteria, and higher expression of photosynthetic genes were selected. The p-adaptations are directly interpretable, while the coadapted gene expressions reflect the physiological requirements for the adaptation. The integration of this information characterizes when and where the population has experienced environmental stress and how the population responded at the molecular level.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

The genome sequence of thale cress, Arabidopsis thaliana (Heynh., 1842)

We present a genome assembly of an Arabidopsis thaliana specimen (thale cress; Streptophyta; Magnoliopsida; Brassicales; Brassicaceae). The genome sequence spans 138 megabases. Most of the assembly (98.76%) is scaffolded into five chromosomal pseudomolecules. The mitochondrial and plastid genomes were also assembled and are 368.8 and 154.5 kilobases in length respectively.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

An Evolutionary Approach to the History of Barley (Hordeum vulgare) Cultivation in the Canary Islands

The Canary Islands are an archipelago that lies about 100 km west of North Africa. Barley (Hordeum vulgare) has been continuously cultivated since the colonization of the islands. To investigate the agricultural history of the islands, the DNA from multiple individuals of six extant landraces of barley was sequenced, and the resulting data were analyzed with ABC modeling. Estimates of separation times of barley populations on the different islands and the mainland were congruent with archaeological dating of the earliest settlements on the islands. The results of the genetic analyses were consistent with the continuous cultivation of barley on Lanzarote island since it was first colonized, but suggested cultivation was carried out at a smaller scale than on Gran Canaria and Tenerife. Contrary to archaeological evidence and early written historical sources, the genetic analyses suggest that barley was cultivated on a larger scale on Tenerife than on Gran Canaria. The genetic analysis of contemporary barley added support to the dating of the colonization of the islands and pointed to the need for more archaeological data concerning barley cultivation on Tenerife.

SPLATCHE3: simulation of serial genetic data under spatially explicit evolutionary scenarios including long-distance dispersal

SUMMARY

SPLATCHE3 simulates genetic data under a variety of spatially explicit evolutionary scenarios, extending previous versions of the framework. The new capabilities include long-distance migration, spatially and temporally heterogeneous short-scale migrations, alternative hybridization models, simulation of serial samples of genetic data and a large variety of DNA mutation models. These implementations have been applied independently to various studies, but grouped together in the current version.

AVAILABILITY AND IMPLEMENTATION

SPLATCHE3 is written in C++ and is freely available for non-commercial use from the website http://www.splatche.com/splatche3. It includes console versions for Linux, MacOs and Windows and a user-friendly GUI for Windows, as well as detailed documentation and ready-to-use examples.

Occurrence of apomictic conspecifics and ecological preferences rather than colonization history govern the geographic distribution of sexual Potentilla puberula

The geographic distribution of sexual-apomictic taxa (i.e., comprising individuals usually reproducing either sexually or asexually via seeds) is traditionally thought to be driven by their ecological preferences and colonization histories. Where sexuals and apomicts get into contact with each other, competitive and reproductive interactions can interfere with these factors, an aspect which hitherto received little attention in biogeographic studies. We disentangled and quantified the relative effects of the three factors on the distribution of tetraploid sexuals in Potentilla puberula in a latitudinal transect through the Eastern European Alps, in which they are codistributed with penta-, hepta-, and octoploid apomictic conspecifics. Effects were explored by means of binomial generalized linear regression models combining a single with a multiple predictor approach. Postglacial colonization history was inferred from population genetic variation (AFLPs and cpDNA) and quantified using a cost distance metric. The study was based on 235 populations, which were purely sexual, purely apomictic, or of mixed reproductive mode. The occurrence of apomicts explained most of the variation in the distribution of sexuals (31%). Specifically, the presence of sexual tetraploids was negatively related to the presence of each of the three apomictic cytotypes. Effects of ecological preferences were substantial too (7% and 12% of the total variation explained by ecological preferences alone, or jointly with apomicts' occurrence, respectively). In contrast, colonization history had negligible effects on the occurrence of sexuals. Taken together, our results highlight the potentially high impact of reproductive interactions on the geographic distribution of sexual and apomictic conspecifics and that resultant mutual exclusion interrelates to ecological differentiation, a situation potentially promoting their local coexistence.

Remember where you came from: ABA insensitivity is epigenetically inherited in mesophyll, but not seeds

Plants transmit their experiences of environmental conditions to their progeny through epigenetic inheritance, improving their progeny's fitness under prevailing conditions. Though ABA is known to regulate epigenetic-modification genes, no strong phenotypic link between those genes and intergenerational "memory" has been shown. Previously, we demonstrated that mesophyll insensitivity to ABA (FBPase::abi1-1{fa} transgenic plants) results in a range of developmental phenotypes, including early growth vigor and early flowering (i.e., stress-escape behavior). Here, we show that null plants, used as controls (segregates of FBPase::abi1 that are homozygote descendants of a heterozygous transgenic plant, but do not contain the transformed abi1-1 gene) phenotypically resembled their FBPase::abi1-1 parents. However, in germination and early seedling development assays, null segregants resembled WT plants. These FBPase::abi1-1 null segregants mesophyll-related phenotypes were reproducible and stable for at least three generations. These results suggest that the heritability of stress response is linked to ABA's epigenetic regulatory effect through ABI1 and mesophyll-related traits. The discrepancy between the epigenetic heritability of seed and mesophyll-related traits is an example of the complexity of epigenetic regulation, which is both gene and process-specific, and may be attributed to the fine-tuning of tradeoffs between flowering time, growth rate and levels of risk that allow annual plants to optimize their fitness in uncertain environments.

Variation among 532 genomes unveils the origin and evolutionary history of a global insect herbivore

The diamondback moth, Plutella xylostella is a cosmopolitan pest that has evolved resistance to all classes of insecticide, and costs the world economy an estimated US $4-5 billion annually. We analyse patterns of variation among 532 P. xylostella genomes, representing a worldwide sample of 114 populations. We find evidence that suggests South America is the geographical area of origin of this species, challenging earlier hypotheses of an Old-World origin. Our analysis indicates that Plutella xylostella has experienced three major expansions across the world, mainly facilitated by European colonization and global trade. We identify genomic signatures of selection in genes related to metabolic and signaling pathways that could be evidence of environmental adaptation. This evolutionary history of P. xylostella provides insights into transoceanic movements that have enabled it to become a worldwide pest.

Global gene flow releases invasive plants from environmental constraints on genetic diversity

When plants establish outside their native range, their ability to adapt to the new environment is influenced by both demography and dispersal. However, the relative importance of these two factors is poorly understood. To quantify the influence of demography and dispersal on patterns of genetic diversity underlying adaptation, we used data from a globally distributed demographic research network comprising 35 native and 18 nonnative populations of Plantago lanceolata Species-specific simulation experiments showed that dispersal would dilute demographic influences on genetic diversity at local scales. Populations in the native European range had strong spatial genetic structure associated with geographic distance and precipitation seasonality. In contrast, nonnative populations had weaker spatial genetic structure that was not associated with environmental gradients but with higher within-population genetic diversity. Our findings show that dispersal caused by repeated, long-distance, human-mediated introductions has allowed invasive plant populations to overcome environmental constraints on genetic diversity, even without strong demographic changes. The impact of invasive plants may, therefore, increase with repeated introductions, highlighting the need to constrain future introductions of species even if they already exist in an area.

Functional variants of DOG1 control seed chilling responses and variation in seasonal life-history strategies in Arabidopsis thaliana

The seasonal timing of seed germination is critical for plant fitness in different climates. To germinate at the right time of year, seeds respond to seasonal environmental cues, such as cold temperatures. We characterized genetic variation in seed dormancy responses to cold across the geographic range of a widespread annual plant. Induction of secondary seed dormancy during winter conditions (which restricts germination to autumn) was positively correlated with flowering time, constructing winter and spring seasonal life-history strategies. Variation in seed chilling responses was strongly associated with functional variants of a known dormancy gene. These variants showed evidence of ancient diversification associated with Pleistocene glacial cycles, and were associated with climate gradients across the species’ geographical range.

The seasonal timing of seed germination determines a plant’s realized environmental niche, and is important for adaptation to climate. The timing of seasonal germination depends on patterns of seed dormancy release or induction by cold and interacts with flowering-time variation to construct different seasonal life histories. To characterize the genetic basis and climatic associations of natural variation in seed chilling responses and associated life-history syndromes, we selected 559 fully sequenced accessions of the model annual species Arabidopsis thaliana from across a wide climate range and scored each for seed germination across a range of 13 cold stratification treatments, as well as the timing of flowering and senescence. Germination strategies varied continuously along 2 major axes: 1) Overall germination fraction and 2) induction vs. release of dormancy by cold. Natural variation in seed responses to chilling was correlated with flowering time and senescence to create a range of seasonal life-history syndromes. Genome-wide association identified several loci associated with natural variation in seed chilling responses, including a known functional polymorphism in the self-binding domain of the candidate gene DOG1. A phylogeny of DOG1 haplotypes revealed ancient divergence of these functional variants associated with periods of Pleistocene climate change, and Gradient Forest analysis showed that allele turnover of candidate SNPs was significantly associated with climate gradients. These results provide evidence that A. thaliana’s germination niche and correlated life-history syndromes are shaped by past climate cycles, as well as local adaptation to contemporary climate.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

Network-based hierarchical population structure analysis for large genomic data sets

Analysis of population structure in natural populations using genetic data is a common practice in ecological and evolutionary studies. With large genomic data sets of populations now appearing more frequently across the taxonomic spectrum, it is becoming increasingly possible to reveal many hierarchical levels of structure, including fine-scale genetic clusters. To analyze these data sets, methods need to be appropriately suited to the challenges of extracting multilevel structure from whole-genome data. Here, we present a network-based approach for constructing population structure representations from genetic data. The use of community-detection algorithms from network theory generates a natural hierarchical perspective on the representation that the method produces. The method is computationally efficient, and it requires relatively few assumptions regarding the biological processes that underlie the data. We show the approach by analyzing population structure in the model plant species Arabidopsis thaliana and in human populations. These examples illustrate how network-based approaches for population structure analysis are well-suited to extracting valuable ecological and evolutionary information in the era of large genomic data sets.

Biogeographic Patterns of Ectomycorrhizal Fungal Communities Associated With Castanopsis sieboldii Across the Japanese Archipelago

Biogeographic patterns in ectomycorrhizal (ECM) fungal communities and their drivers have been elucidated, including effects of host tree species and abiotic (climatic and edaphic) conditions. At these geographic scales, genotypic diversity and composition of single host tree species change with spatial and environmental gradients, reflecting their historical dispersal events. However, whether the host genotypes can be associated with the biogeographic patterns of ECM communities remains unclear. We investigated the biogeographic pattern of ECM fungal community associated with the single host species Castanopsis sieboldii (Fagaceae), whose genotypic diversity and composition across the Japanese archipelago has already been evaluated. ECM communities were investigated in 12 mature Castanopsis-dominated forests covering almost the entire distribution range of C. sieboldii, and we quantified the effect of host genotypes on the biogeographic pattern of ECM fungal communities. Richness and community composition of ECM fungi changed with latitude and longitude; these biogeographic changes of ECM community were significantly correlated with host genotypic variables. Quantitative analyses showed a higher relative explanatory power of climatic and spatial variables than that of host genotypic variables for the biogeographic patterns in the ECM community. Our results suggest historical events of host dispersal can affect the biogeographic patterns of the ECM fungal community, while their explanation power was lower than that for climatic filtering and/or fungal dispersal.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

Helitron distribution in Brassicaceae and whole Genome Helitron density as a character for distinguishing plant species

BACKGROUND

Helitron is a rolling-circle DNA transposon; it plays an important role in plant evolution. However, Helitron distribution and contribution to evolution at the family level have not been previously investigated.

RESULTS

We developed the software easy-to-annotate Helitron (EAHelitron), a Unix-like command line, and used it to identify Helitrons in a wide range of 53 plant genomes (including 13 Brassicaceae species). We determined Helitron density (abundance/Mb) and visualized and examined Helitron distribution patterns. We identified more than 104,653 Helitrons, including many new Helitrons not predicted by other software. Whole genome Helitron density is independent from genome size and shows stability at the species level. Using linear discriminant analysis, de novo genomes (next-generation sequencing) were successfully classified into Arabidopsis thaliana groups. For most Brassicaceae species, Helitron density negatively correlated with gene density, and Helitron distribution patterns were similar to those of A. thaliana. They preferentially inserted into sequence around the centromere and intergenic region. We also associated 13 Helitron polymorphism loci with flowering-time phenotypes in 18 A. thaliana ecotypes.

CONCLUSION

EAHelitron is a fast and efficient tool to identify new Helitrons. Whole genome Helitron density can be an informative character for plant classification. Helitron insertion polymorphism could be used in association analysis.

On the postglacial spread of human commensal Arabidopsis thaliana: journey to the East

With more sequenced genomes, our understanding of the demographic history of Arabidopsis thaliana is rapidly expanding. However, no-one has yet compiled previous data to investigate patterns of genetic variation across Eurasia. While sub-Saharan accessions have been reported to be the most divergent group, in the nuclear genome we found accessions from Yunnan, China to be genetically closest to the sub-Saharan group. In chloroplast, several deeply diverged haplogroups exist only in Eurasia, and African populations have lower variation in many haplogroups that they share with the Eurasian populations. These patterns cannot be easily explained by a single out-of-Africa event suggested previously. For more recent demographic history, we dated the nonrelict expansion to 10 ka. In the Chinese Yangtze nonrelicts, we found clear traces of gene flow with local relicts, and genes under strong selection were enriched for traces of relict introgression, especially those related to biotic and immune responses. The results suggest the ability of nonrelicts to obtain locally adaptive alleles through admixture with relicts is important for the expansion across environmental gradients of Eurasia. Our re-analyses provide another model for the early history as well as elucidating factors contributing to the recent demographic turnover event of this species.

SNP genotyping elucidates the genetic diversity of Magna Graecia grapevine germplasm and its historical origin and dissemination

BACKGROUND

Magna Graecia is the ancient name for the modern geopolitical region of South Italy extensively populated by Greek colonizers, shown by archeological and historical evidence to be the oldest wine growing region of Italy, crucial for the spread of specialized viticulture around Mediterranean shores. Here, the genetic diversity of Magna Graecia grape germplasm was assessed and its role in grapevine propagation around the Mediterranean basin was underlined.

RESULTS

A large collection of grapevines from Magna Graecia was compared with germplasm from Georgia to the Iberian Peninsula using the 18 K SNP array. A high level of genetic diversity of the analyzed germplasm was determined; clustering, structure analysis and DAPC (Discriminant Analysis of Principal Components) highlighted the genetic relationships among genotypes from South Italy and the Eastern Mediterranean (Greece). Gene flow from east (Georgia) to west (Iberian Peninsula) was identified throughout the large number of detected admixed samples. Pedigree analysis showed a complex and well-structured network of first degree relationships, where the cultivars from Magna Graecia were mainly involved.

CONCLUSIONS

This study provided evidence that Magna Graecia germplasm was shaped by historical events that occurred in the area due to the robust link between South Italian and Greek genotypes, as well as, by the availability of different thermal resources for cultivars growing in such different winegrowing areas. The uniqueness of this ampelographic platform was mainly an outcome of complex natural or human-driven crosses involving elite cultivars.

Demographic inferences after a range expansion can be biased: the test case of the blacktip reef shark (Carcharhinus melanopterus)

The evolutionary history of species is a dynamic process as they modify, expand, and contract their spatial distributions over time. Range expansions (REs) occur through a series of founder events that are followed by migration among neighboring demes. The process usually results in structured metapopulations and leaves a distinct signature in the genetic variability of species. Explicitly modeling the consequences of complex demographic events such as REs is computationally very intensive. Here we propose an an alternative approach that requires less computational effort than a comprehensive RE model, but that can recover the demography of species undergoing a RE, by combining spatially explicit modelling with simplified but realistic metapopulation models. We examine the demographic and colonization history of Carcharhinus melanopterus, an abundant reef-associated shark, as a test case. We first used a population genomics approach to statistically confirm the occurrence of a RE in C. melanopterus, and identify its origin in the Indo-Australian Archipelago. Spatial genetic modelling identified two waves of stepping-stone colonization: an eastward wave moving through the Pacific and a westward one moving through the Indian Ocean. We show that metapopulation models best describe the demographic history of this species and that not accounting for this may lead to incorrectly interpreting the observed genetic variation as signals of widespread population bottlenecks. Our study highlights insights that can be gained about demography by coupling metapopulation models with spatial modeling and underscores the need for cautious interpretation of population genetic data when advancing conservation priorities.

Archaic lineages broaden our view on the history of Arabidopsis thaliana

Contents Summary 1194 I. Introduction 1194 II. Origin of the A. thaliana species 1194 III. The classic model of the history of A. thaliana 1195 IV. New genomic data from outside Eurasia challenge our view of A. thaliana history 1195 V. Conclusions 1197 Acknowledgements 1197 References 1197 SUMMARY: Natural variation in Arabidopsis thaliana has contributed to discoveries in diverse areas of plant biology. While A. thaliana has typically been considered a weed associated primarily with human-mediated environments, including agricultural and urban sites and railways, it has recently been shown that it is also native in remote natural areas, including high altitude sites in Eurasia and Africa, from the Atlas mountains in Morocco to the afro-alpine regions in Eastern and South Africa to Yunnan in China, the Himalayas and the Tibetan Plateau. This finding suggests that while A. thaliana has been extensively studied in Europe and Western Asia there are still many open questions about its population history, genotype-phenotype relationships and mechanisms of adaptation.

How Do Cold-Adapted Plants Respond to Climatic Cycles? Interglacial Expansion Explains Current Distribution and Genomic Diversity in Primula farinosa L

Understanding the effects of past climatic fluctuations on the distribution and population-size dynamics of cold-adapted species is essential for predicting their responses to ongoing global climate change. In spite of the heterogeneity of cold-adapted species, two main contrasting hypotheses have been proposed to explain their responses to Late Quaternary glacial cycles, namely, the interglacial contraction versus the interglacial expansion hypotheses. Here, we use the cold-adapted plant Primula farinosa to test two demographic models under each of the two alternative hypotheses and a fifth, null model. We first approximate the time and extent of demographic contractions and expansions during the Late Quaternary by projecting species distribution models across the last 72 ka. We also generate genome-wide sequence data using a Reduced Representation Library approach to reconstruct the spatial structure, genetic diversity, and phylogenetic relationships of lineages within P. farinosa. Finally, by integrating the results of climatic and genomic analyses in an Approximate Bayesian Computation framework, we propose the most likely model for the extent and direction of population-size changes in $P$. farinosa through the Late Quaternary. Our results support the interglacial expansion of $P$. farinosa, differing from the prevailing paradigm that the observed distribution of cold-adapted species currently fragmented in high altitude and latitude regions reflects the consequences of postglacial contraction processes.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500–6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east–west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

The Scandinavian peninsula was the last part of Europe to be colonized after the Last Glacial Maximum. The migration routes, cultural networks, and the genetic makeup of the first Scandinavians remain elusive and several hypotheses exist based on archaeology, climate modeling, and genetics. By analyzing the genomes of early Scandinavian hunter-gatherers, we show that their migrations followed two routes: one from the south and another from the northeast along the ice-free Norwegian Atlantic coast. These groups met and mixed in Scandinavia, creating a population more diverse than contemporaneous central and western European hunter-gatherers. As northern Europe is associated with cold and low light conditions, we investigated genomic patterns of adaptation to these conditions and genes known to be involved in skin pigmentation. We demonstrate that Mesolithic Scandinavians had higher levels of light pigmentation variants compared to the respective source populations of the migrations, suggesting adaptation to low light levels and a surprising signal of genetic continuity in TMEM131, a gene that may be involved in long-term adaptation to the cold.

High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition

BACKGROUND

Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition.

RESULTS

Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one.

CONCLUSIONS

Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.

Chlorosis caused by two recessively interacting genes reveals a role of RNA helicase in hybrid breakdown in Arabidopsis thaliana

Hybrids often differ in fitness from their parents. They may be superior, translating into hybrid vigour or heterosis, but they may also be markedly inferior, because of hybrid weakness or incompatibility. The underlying genetic causes for the latter can often be traced back to genes that evolve rapidly because of sexual or host-pathogen conflicts. Hybrid weakness may manifest itself only in later generations, in a phenomenon called hybrid breakdown. We have characterized a case of hybrid breakdown among two Arabidopsis thaliana accessions, Shahdara (Sha, Tajikistan) and Lövvik-5 (Lov-5, Northern Sweden). In addition to chlorosis, a fraction of the F₂ plants have defects in leaf and embryo development, and reduced photosynthetic efficiency. Hybrid chlorosis is due to two major-effect loci, of which one, originating from Lov-5, appears to encode an RNA helicase (AtRH18). To examine the role of the chlorosis allele in the Lövvik area, in addition to eight accessions collected in 2009, we collected another 240 accessions from 15 collections sites, including Lövvik, from Northern Sweden in 2015. Genotyping revealed that Lövvik collection site is separated from the rest. Crosses between 109 accessions from this area and Sha revealed 85 cases of hybrid chlorosis, indicating that the chlorosis-causing allele is common in this area. These results suggest that hybrid breakdown alleles not only occur at rapidly evolving loci, but also at genes that code for conserved processes.

African genomes illuminate the early history and transition to selfing in Arabidopsis thaliana

Over the past 20 y, many studies have examined the history of the plant ecological and molecular model, Arabidopsis thaliana, in Europe and North America. Although these studies informed us about the recent history of the species, the early history has remained elusive. In a large-scale genomic analysis of African A. thaliana, we sequenced the genomes of 78 modern and herbarium samples from Africa and analyzed these together with over 1,000 previously sequenced Eurasian samples. In striking contrast to expectations, we find that all African individuals sampled are native to this continent, including those from sub-Saharan Africa. Moreover, we show that Africa harbors the greatest variation and represents the deepest history in the A. thaliana lineage. Our results also reveal evidence that selfing, a major defining characteristic of the species, evolved in a single geographic region, best represented today within Africa. Demographic inference supports a model in which the ancestral A. thaliana population began to split by 120-90 kya, during the last interglacial and Abbassia pluvial, and Eurasian populations subsequently separated from one another at around 40 kya. This bears striking similarities to the patterns observed for diverse species, including humans, implying a key role for climatic events during interglacial and pluvial periods in shaping the histories and current distributions of a wide range of species.

On the post-glacial spread of human commensal Arabidopsis thaliana

Recent work has shown that Arabidopsis thaliana contains genetic groups originating from different ice age refugia, with one particular group comprising over 95% of the current worldwide population. In Europe, relicts of other groups can be found in local populations along the Mediterranean Sea. Here we provide evidence that these 'relicts' occupied post-glacial Eurasia first and were later replaced by the invading 'non-relicts', which expanded through the east-west axis of Eurasia, leaving traces of admixture in the north and south of the species range. The non-relict expansion was likely associated with human activity and led to a demographic replacement similar to what occurred in humans. Introgressed genomic regions from relicts are associated with flowering time and enriched for genes associated with environmental conditions, such as root cap development or metal ion trans-membrane transport, which suggest that admixture with locally adapted relicts helped the non-relicts colonize new habitats.

Demographic history of the trace metal hyperaccumulator Noccaea caerulescens (J. Presl and C. Presl) F. K. Mey. in Western Europe

Noccaea caerulescens (Brassicaceae) is a major pseudometallophyte model for the investigation of the genetics and evolution of metal hyperaccumulation in plants. We studied the population genetics and demographic history of this species to advance the understanding of among-population differences in metal hyperaccumulation and tolerance abilities. Sampling of seven to 30 plants was carried out in 62 sites in Western Europe. Genotyping was carried out using a combination of new chloroplast and nuclear neutral markers. A strong genetic structure was detected, allowing the definition of three genetic subunits. Subunits showed a good geographic coherence. Accordingly, distant metallicolous populations generally belonged to distinct subunits. Approximate Bayesian computation analysis of demographic scenarios among subunits further supported a primary isolation of populations from the southern Massif Central prior to last glacial maximum, whereas northern populations may have derived during postglacial recolonization events. Estimated divergence times among subunits were rather recent in comparison with the species history, but certainly before the establishment of anthropogenic metalliferous sites. Our results suggest that the large-scale genetic structure of N. caerulescens populations pre-existed to the local adaptation to metalliferous sites. The population structure of quantitative variation for metal-related adaptive traits must have established independently in isolated gene pools. However, features of the most divergent genetic unit (e.g. extreme levels of Cd accumulation observed in previous studies) question the putative relationships between adaptive evolution of metal-related traits and subunits isolation. Finally, admixture signals among distant metallicolous populations suggest a putative role of human activities in facilitating long-distance genetic exchanges.

Population genomics of C. melanopterus using target gene capture data: demographic inferences and conservation perspectives

Population genetics studies on non-model organisms typically involve sampling few markers from multiple individuals. Next-generation sequencing approaches open up the possibility of sampling many more markers from fewer individuals to address the same questions. Here, we applied a target gene capture method to deep sequence ~1000 independent autosomal regions of a non-model organism, the blacktip reef shark (Carcharhinus melanopterus). We devised a sampling scheme based on the predictions of theoretical studies of metapopulations to show that sampling few individuals, but many loci, can be extremely informative to reconstruct the evolutionary history of species. We collected data from a single deme (SID) from Northern Australia and from a scattered sampling representing various locations throughout the Indian Ocean (SCD). We explored the genealogical signature of population dynamics detected from both sampling schemes using an ABC algorithm. We then contrasted these results with those obtained by fitting the data to a non-equilibrium finite island model. Both approaches supported an Nm value ~40, consistent with philopatry in this species. Finally, we demonstrate through simulation that metapopulations exhibit greater resilience to recent changes in effective size compared to unstructured populations. We propose an empirical approach to detect recent bottlenecks based on our sampling scheme.

On the role played by the carrying capacity and the ancestral population size during a range expansion

Most species are structured and various population genetics models have been proposed to investigate their history. For mathematical tractability, most of these models make the simplifying assumption of equilibrium. Here we focus on the properties of a nonequilibrium spatial explicit model, range expansions (REs). Despite their abundance, many details of their genetic consequences need yet to be fully investigated. The model we studied is characterized by four main parameters: the effective population size of each deme (N), the migration rate per generation per deme (m), the time of the expansion (T_exp) and the effective size of the deme from which the expansion started (N_anc). By means of extensive coalescent simulations, we focused on two aspects of range expansions for fixed Nm: (1) the separate influence of N and m and (2) the role of N_anc. We compared our results with an equilibrium stepping stone model and found two main features typical of REs: an excess of rare variants for larger N and a complex interaction between N, T_exp and N_anc in shaping the degree of population differentiation (which depends only on Nm in the stepping stone model). Finally, we developed an approximate Bayesian computation approach to jointly estimate N and m and to infer N_anc. When applied to pseudo-observed data sets, we could correctly recover both N and m (but not N_anc), provided a large number of demes were sampled. These findings highlight how it will be possible to estimate the dispersal rate in nonequilibrium metapopulations by using population genetics approaches.

Investigating past range dynamics for a weed of cultivation, Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.

1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana

Arabidopsis thaliana serves as a model organism for the study of fundamental physiological, cellular, and molecular processes. It has also greatly advanced our understanding of intraspecific genome variation. We present a detailed map of variation in 1,135 high-quality re-sequenced natural inbred lines representing the native Eurasian and North African range and recently colonized North America. We identify relict populations that continue to inhabit ancestral habitats, primarily in the Iberian Peninsula. They have mixed with a lineage that has spread to northern latitudes from an unknown glacial refugium and is now found in a much broader spectrum of habitats. Insights into the history of the species and the fine-scale distribution of genetic diversity provide the basis for full exploitation of A. thaliana natural variation through integration of genomes and epigenomes with molecular and non-molecular phenotypes.

The Arabidopsis thaliana mobilome and its impact at the species level

Transposable elements (TEs) are powerful motors of genome evolution yet a comprehensive assessment of recent transposition activity at the species level is lacking for most organisms. Here, using genome sequencing data for 211 Arabidopsis thaliana accessions taken from across the globe, we identify thousands of recent transposition events involving half of the 326 TE families annotated in this plant species. We further show that the composition and activity of the 'mobilome' vary extensively between accessions in relation to climate and genetic factors. Moreover, TEs insert equally throughout the genome and are rapidly purged by natural selection from gene-rich regions because they frequently affect genes, in multiple ways. Remarkably, loci controlling adaptive responses to the environment are the most frequent transposition targets observed. These findings demonstrate the pervasive, species-wide impact that a rich mobilome can have and the importance of transposition as a recurrent generator of large-effect alleles.

Genomic signature of successful colonization of Eurasia by the allopolyploid shepherd's purse (Capsella bursa-pastoris)

Polyploidization is a dominant feature of flowering plant evolution. However, detailed genomic analyses of the interpopulation diversification of polyploids following genome duplication are still in their infancy, mainly because of methodological limits, both in terms of sequencing and computational analyses. The shepherd's purse (Capsella bursa-pastoris) is one of the most common weed species in the world. It is highly self-fertilizing, and recent genomic data indicate that it is an allopolyploid, resulting from hybridization between the ancestors of the diploid species Capsella grandiflora and Capsella orientalis. Here, we investigated the genomic diversity of C. bursa-pastoris, its population structure and demographic history, following allopolyploidization in Eurasia. To that end, we genotyped 261 C. bursa-pastoris accessions spread across Europe, the Middle East and Asia, using genotyping-by-sequencing, leading to a total of 4274 SNPs after quality control. Bayesian clustering analyses revealed three distinct genetic clusters in Eurasia: one cluster grouping samples from Western Europe and Southeastern Siberia, the second one centred on Eastern Asia and the third one in the Middle East. Approximate Bayesian computation (ABC) supported the hypothesis that C. bursa-pastoris underwent a typical colonization history involving low gene flow among colonizing populations, likely starting from the Middle East towards Europe and followed by successive human-mediated expansions into Eastern Asia. Altogether, these findings bring new insights into the recent multistage colonization history of the allotetraploid C. bursa-pastoris and highlight ABC and genotyping-by-sequencing data as promising but still challenging tools to infer demographic histories of selfing allopolyploids.

Genome Wide Association Mapping for the Tolerance to the Polyamine Oxidase Inhibitor Guazatine in Arabidopsis thaliana

Guazatine is a potent inhibitor of polyamine oxidase (PAO) activity. In agriculture, guazatine is used as non-systemic contact fungicide efficient in the protection of cereals and citrus fruits against disease. The composition of guazatine is complex, mainly constituted by a mixture of synthetic guanidated polyamines (polyaminoguanidines). Here, we have studied the effects from exposure to guazatine in the weed Arabidopsis thaliana. We report that micromolar concentrations of guazatine are sufficient to inhibit growth of Arabidopsis seedlings and induce chlorosis, whereas germination is barely affected. We observed the occurrence of quantitative variation in the response to guazatine between 107 randomly chosen Arabidopsis accessions. This enabled us to undertake genome-wide association (GWA) mapping that identified a locus on chromosome one associated with guazatine tolerance. CHLOROPHYLLASE 1 (CLH1) within this locus was studied as candidate gene, together with its paralog (CLH2). The analysis of independent clh1-2, clh1-3, clh2-3, clh2-2, and double clh1-2 clh2-3 mutant alleles indicated that CLH1 and/or CLH2 loss-of-function or expression down-regulation promote guazatine tolerance in Arabidopsis. We report a natural mechanism by which Arabidopsis populations can overcome toxicity by the fungicide guazatine.

A Polynucleotide Repeat Expansion Causing Temperature-Sensitivity Persists in Wild Irish Accessions of Arabidopsis thaliana

Triplet repeat expansions underlie several human genetic diseases such as Huntington's disease and Friedreich's ataxia. Although such mutations are primarily known from humans, a triplet expansion associated genetic defect has also been reported at the IIL1 locus in the Bur-0 accession of the model plant Arabidopsis thaliana. The IIL1 triplet expansion is an example of cryptic genetic variation as its phenotypic effects are seen only under genetic or environmental perturbation, with high temperatures resulting in a growth defect. Here we demonstrate that the IIL1 triplet expansion associated growth defect is not a general stress response and is specific to particular environmental perturbations. We also confirm and map genetic modifiers that suppress the effect of IIL1 triplet repeat expansion. By collecting and analyzing accessions from the island of Ireland, we recover the repeat expansion in wild populations suggesting that the repeat expansion has persisted at least 60 years in Ireland. Through genome-wide genotyping, we show that the repeat expansion is present in diverse Irish populations. Our findings indicate that even deleterious alleles can persist in populations if their effect is conditional. Our study demonstrates that analysis of groups of wild populations is a powerful tool for understanding the dynamics of cryptic genetic variation.

Adapting through glacial cycles: insights from a long-lived tree (Taxus baccata)

Despite the large body of research devoted to understanding the role of Quaternary glacial cycles in the genetic divergence of European trees, the differential contribution of geographic isolation and/or environmental adaptation in creating population genetic divergence remains unexplored. In this study, we used a long-lived tree (Taxus baccata) as a model species to investigate the impact of Quaternary climatic changes on genetic diversity via neutral (isolation-by-distance) and selective (isolation-by-adaptation) processes. We applied approximate Bayesian computation to genetic data to infer its demographic history, and combined this information with past and present climatic data to assess the role of environment and geography in the observed patterns of genetic structure. We found evidence that yew colonized Europe from the East, and that European samples diverged into two groups (Western, Eastern) at the beginning of the Quaternary glaciations, c. 2.2 Myr before present. Apart from the expected effects of geographical isolation during glacials, we discovered a significant role of environmental adaptation during interglacials at the origin of genetic divergence between both groups. This process may be common in other organisms, providing new research lines to explore the effect of Quaternary climatic factors on present-day patterns of genetic diversity.

TESS3: fast inference of spatial population structure and genome scans for selection

Geography and landscape are important determinants of genetic variation in natural populations, and several ancestry estimation methods have been proposed to investigate population structure using genetic and geographic data simultaneously. Those approaches are often based on computer-intensive stochastic simulations and do not scale with the dimensions of the data sets generated by high-throughput sequencing technologies. There is a growing demand for faster algorithms able to analyse genomewide patterns of population genetic variation in their geographic context. In this study, we present TESS3, a major update of the spatial ancestry estimation program TESS. By combining matrix factorization and spatial statistical methods, TESS3 provides estimates of ancestry coefficients with accuracy comparable to TESS and with run-times much faster than the Bayesian version. In addition, the TESS3 program can be used to perform genome scans for selection, and separate adaptive from nonadaptive genetic variation using ancestral allele frequency differentiation tests. The main features of TESS3 are illustrated using simulated data and analysing genomic data from European lines of the plant species Arabidopsis thaliana.

Serial founder effects and genetic differentiation during worldwide range expansion of monarch butterflies

Range expansions can result in founder effects, increasing genetic differentiation between expanding populations and reducing genetic diversity along the expansion front. However, few studies have addressed these effects in long-distance migratory species, for which high dispersal ability might counter the effects of genetic drift. Monarchs (Danaus plexippus) are best known for undertaking a long-distance annual migration in North America, but have also dispersed around the world to form populations that do not migrate or travel only short distances. Here, we used microsatellite markers to assess genetic differentiation among 18 monarch populations and to determine worldwide colonization routes. Our results indicate that North American monarch populations connected by land show limited differentiation, probably because of the monarch's ability to migrate long distances. Conversely, we found high genetic differentiation between populations separated by large bodies of water. Moreover, we show evidence for serial founder effects across the Pacific, suggesting stepwise dispersal from a North American origin. These findings demonstrate that genetic drift played a major role in shaping allele frequencies and created genetic differentiation among newly formed populations. Thus, range expansion can give rise to genetic differentiation and declines in genetic diversity, even in highly mobile species.