SPLATCHE2: a spatially explicit simulation framework for complex demography, genetic admixture and recombination

SUMMARY: SPLATCHE2 is a program to simulate the demography of populations and the resulting molecular diversity for a wide range of evolutionary scenarios. The spatially explicit simulation framework can account for environmental heterogeneity and fluctuations, and it can manage multiple population sources. A coalescent-based approach is used to generate genetic markers mostly used in population genetics studies (DNA sequences, SNPs, STRs or RFLPs). Various combinations of independent, fully or partially linked genetic markers can be produced under a recombination model based on the ancestral recombination graph. Competition between two populations (or species) can also be simulated with user-defined levels of admixture between the two populations. SPLATCHE2 may be used to generate the expected genetic diversity under complex demographic scenarios and can thus serve to test null hypotheses. For model parameter estimation, SPLATCHE2 can easily be integrated into an Approximate Bayesian Computation (ABC) framework. Availability and implementation: SPLATCHE2 is a C++ program compiled for Windows and Linux platforms. It is freely available at www.splatche.com, together with its related documentation and example data. CONTACT: mathias.currat@unige.ch     

cdpop: A spatially explicit cost distance population genetics program

Spatially explicit simulation of gene flow in complex landscapes is essential to explain observed population responses and provide a foundation for landscape genetics. To address this need, we wrote a spatially explicit, individual-based population genetics model (cdpop). The model implements individual-based population modelling with Mendelian inheritance and k-allele mutation on a resistant landscape. The model simulates changes in population and genotypes through time as functions of individual based movement, reproduction, mortality and dispersal on a continuous cost surface. This model will be a valuable tool for the study of landscape genetics by increasing our understanding about the effects of life history, vagility and differential models of landscape resistance on the genetic structure of populations in complex landscapes.     

Marked genetic structuring and extreme dispersal limitation in the Pyrenean brook newt Calotriton asper (Amphibia: Salamandridae) revealed by genome-wide AFLP but not mtDNA

Direct estimation of dispersal rates at large geographic scales can be technically and logistically challenging, especially in small animals of low vagility like amphibians. The use of molecular markers to reveal patterns of genetic structure provides an indirect way to infer dispersal rates and patterns of recent and historical gene flow among populations. Here, we use mitochondrial DNA (mtDNA) sequence data and genome-wide amplified fragment length polymorphism markers to examine population structure in the Pyrenean brook newt ( Calotriton asper ) across four main drainages in the French Pyrenees. mtDNA sequence data (2040 bp) revealed three phylogroups shallowly differentiated and with low genetic diversity. In sharp contrast, variation in 382 amplified fragment length polymorphism loci was high and revealed a clear pattern of isolation by distance consistent with long-term restriction of gene flow at three spatial scales: (i) among all four main drainages, (ii) between sites within drainages, and (iii) even between adjacent populations separated by less than 4 km. The high pairwise F _ST values between localities across numerous loci, together with the high frequency of fixed alleles in several populations, suggests a combination of marked geographic isolation, small population sizes and very limited dispersal in C. asper . The contrasting lack of variation detected in mtDNA sequence data is intriguing and underscores the importance of multilocus approaches to detect true patterns of gene flow in natural populations of amphibians.     

Mate finding, dispersal distances and population growth in invading species: a spatially explicit model

The expansion of populations into new areas is dependent upon dispersal distances and the ability of colonists to find mates. These factors interact through the spatial distribution of individuals. We develop a mechanistic, spatially explicit model to investigate the interaction between dispersal distances and mate finding in expanding populations. At high dispersal distances and low mate finding abilities, population growth was constrained by the inverse density dependent inability of adults to find mates (an Allee effect). In contrast, at low dispersal distances and high mate finding abilities, growth was constrained by the density dependent inability of dispersers to find vacant territories. Population growth was highest in between these extremes. We suggest that these spatial interactions play an important role in the expansion of populations into new areas and that this methodology provides a useful tool for investigating them.     

Coevolution across landscapes: a spatially explicit model of parasitoid-host coevolution

The geographic mosaic theory of coevolution suggests that population spatial structure may have a strong impact on coevolutionary dynamics. Therefore, coevolution must be studied across geographic scales, not just in single populations. To examine the impact of movement rate on coevolutionary dynamics, we developed a spatially explicit model of host–parasitoid coevolution. We described space as a coupled-map lattice and assumed that resistance (defined as the ability of a host to encapsulate a parasitoid egg) and virulence (defined as the successful parasitization of a host) traits were graded and costly. The model explicitly detailed population and evolutionary dynamics. When holding all parameters constant and varying only the movement rate of the host and parasitoid, profoundly different dynamics were observed. We found that fluctuations in the mean levels of resistance and virulence in the global population were greatest when the movement rate of the host and parasitoid was high. In addition, we found that the variation in resistance and virulence levels among neighboring patches was greatest when the movement rates of the host and parasitoid was low. However, as the distance among patches increased, so did the variation in resistance and virulence levels regardless of movement rate. These generalizations did not hold when spatial patterns in the distribution of resistance and virulence traits, such as spirals, were observed. Finally, we found that the evolution of resistance and virulence caused the abundance of hosts to increase and the abundance of parasitoids to decrease. As a result, the spatial distribution of hosts and parasitoids was influenced.     

Origin of spatial genetic structure in an expanding oak population

Spatial genetic structure (SGS) results from the interplay of several demographical processes that are difficult to tease apart. In this study, we explore the specific effects of seed and pollen dispersal and of early postdispersal mortality on the SGS of a seedling cohort (N = 786) recruiting within and around an expanding pedunculate oak (Quercus robur) stand. Using data on dispersal (derived from parentage analysis) and mortality (monitored in the field through two growing seasons), we decompose the overall SGS of the cohort into its components by contrasting the SGS of dispersed (i.e. growing away from their mother tree) vs. nondispersed (i.e. growing beneath their mother tree) and initial vs. surviving seedlings. Patterns differ strongly between nondispersed and dispersed seedlings. Nondispersed seedlings are largely responsible for the positive kinship values observed at short distances in the studied population, whereas dispersed seedlings determine the overall SGS at distances beyond c. 30 m. The paternal alleles of nondispersed seedlings show weak yet significantly positive kinships up to c. 15 m, indicating some limitations in pollen flow that should further promote pedigree structures at short distances. Seedling mortality does not alter SGS, except for a slight increase in the nondispersed group. Field data reveal that mortality in this group is negatively density‐dependent, probably because of small‐scale variation in light conditions. Finally, we observe a remarkable similarity between the SGS of the dispersed seedlings and that of the adults, which probably reflects dispersal processes during the initial expansion of the population. Overall, this study demonstrates that incorporating individual‐level complementary information into analyses can greatly improve the detail and confidence of ecological inferences drawn from SGS.     

Spatial genetic structure within a metallicolous population of Arabidopsis halleri, a clonal, self-incompatible and heavy-metal-tolerant species

Arabidopsis halleri, a close wild relative of A. thaliana, is a clonal, insect-pollinated herb tolerant to heavy metals (Zn, Pd, Cd) and a hyperaccumulator of Zn and Cd. It is of particular interest in the study of evolutionary processes and phytoremediation. However, little is known about its population gene flow patterns and the structure of its genetic diversity. We used five microsatellite loci to investigate the genetic structure at a fine spatial scale (10 cm to 500 m) in a metallicolous population of A. halleri. We also studied the contributions made by clonal propagation and sexual reproduction (seed and pollen dispersal) to the genetic patterns. Clonal diversity was high (D(G) > 0.9). Clonal spread occurs only at short distances (< 1 m). Both clonal spread and limited dispersal, associated with sexual reproduction, contribute to the significant spatial genetic structure revealed by spatial autocorrelation analysis. The shape of the autocorrelogram suggests that seed dispersal is restricted and pollen flow extensive, which may be related to intense activity by insect pollinators. Clonal spread was more extensive in the lowly polluted zone than in the highly polluted zone. This cannot be interpreted as a strategy for promoting the propagation of adapted genotypes under the harshest ecological constraints (highest heavy metal concentrations). The higher fine-scale spatial genetic structure found in the lowly polluted zone can be ascribed to plant densities that were lower than in the highly polluted zone. No evidence of genetic divergence due to spatial heavy metal heterogeneity was found between lowly and highly polluted zones.     

A population genomic approach to map recent positive selection in model species

Based on nearly complete genome sequences from a variety of organisms data on naturally occurring genetic variation on the scale of hundreds of loci to entire genomes have been collected in recent years. In parallel, new statistical tests have been developed to infer evidence of recent positive selection from these data and to localize the target regions of selection in the genome. These methods have now been successfully applied to Drosophila melanogaster , humans, mice and a few plant species. In genomic regions of normal recombination rates, the targets of positive selection have been mapped down to the level of individual genes.     

Evidence of recent and continuous speciation in a biodiversity hotspot: a population genetic approach in southern African gladioli (Gladiolus; Iridaceae)

There has been much debate over the origin of species diversity in biodiversity hotspots, particularly the rate of speciation over extinction and the geographic mode of speciation. Here, we looked at speciation with varying degrees of sympatry in a biodiversity hotspot, focusing on a distinct morphological clade in the Cape Floristic Region in southern Africa, the Gladiolus carinatus species complex (Iridaceae). We investigate the mechanisms involved in population and species differentiation through a combination of ecological and genomic approaches. We estimated spatial and phenological overlap, differences in floral morphology, genetic isolation and genomic selection. A genetic coalescent analysis estimated that the time of divergence between lineages followed the establishment of available habitat in the Cape littoral plain where these species currently overlap geographically. Marked shifts in flowering time and morphology, which act as barriers to gene flow, have developed to varying degrees over the last 0.3–1.4 million years. An amplified fragment length polymorphism genome scan revealed signatures of divergent and balancing selection, although half of the loci consistently behaved neutrally. Divergent species outliers (1%) and floral morph outliers (3%) represent a small proportion of the genome, but these loci produced clear genetic clusters of species and significant associations with floral traits. These results indicate that the G. carinatus complex represents a continuum of recent speciation. We provide further evidence for ecological adaptation in the face of gene flow.     

Genetic differentiation,hybridization and adaptive divergence in two subspecies of the acorn barnacle Tetraclita japonica in the northwestern Pacific

Two acorn barnacles, Tetraclita japonica japonica and Tetraclita japonica formosana, have been recently reclassified as two subspecies, because they are morphologically similar and genetically indistinguishable in mitochondrial DNA sequences. The two barnacles are distinguishable by parietes colour and exhibit parapatric distributions, coexisting in Japan, where T. j. formosana is very low in abundance. Here we investigated the genetic differentiation between the subspecies using 209 polymorphic amplified fragment length polymorphism markers and 341 individuals from 12 locations. The subspecies are genetically highly differentiated (Φ_CT = 0.267). Bayesian analysis and principal component analysis indicate the presence of hybrids in T. j. formosana samples from Japan. Strong differentiation between the northern and southern populations of T. j. japonica was revealed, and a break between Taiwan and Okinawa was also found in T. j. formosana. The differentiation between the two taxa at individual loci does not deviate from neutral expectation, suggesting that the oceanographic pattern which restricts larval dispersal is a more important factor than divergent selection in maintaining genetic and phenotypic differentiation. The T. j. formosana in Japan are probably recent migrants from Okinawa, and their presence in Japan may represent a poleward range shift driven by global warming. This promotes hybridization and might lead to a breakdown of the boundary between the subspecies. However, both local adaptation and larval dispersal are crucial in determining the population structure within each subspecies. Our study provides new insights into the interplay of local adaptation and dispersal in determining the distribution and genetic structure of intertidal biota and the biogeography of the northwestern Pacific.     

Temporal change in genetic integrity suggests loss of local adaptation in a wild Atlantic salmon (Salmo salar) population following introgression by farmed escapees

In some wild Atlantic salmon populations, rapid declines in numbers of wild returning adults has been associated with an increase in the prevalence of farmed salmon. Studies of phenotypic variation have shown that interbreeding between farmed and wild salmon may lead to loss of local adaptation. Yet, few studies have attempted to assess the impact of interbreeding at the genome level, especially among North American populations. Here, we document temporal changes in the genetic makeup of the severely threatened Magaguadavic River salmon population (Bay of Fundy, Canada), a population that might have been impacted by interbreeding with farmed salmon for nearly 20 years. Wild and farmed individuals caught entering the river from 1980 to 2005 were genotyped at 112 single-nucleotide polymorphisms (SNPs), and/or eight microsatellite loci, to scan for potential shifts in adaptive genetic variation. No significant temporal change in microsatellite-based estimates of allele richness or gene diversity was detected in the wild population, despite its precipitous decline in numbers over the last two decades. This might reflect the effect of introgression from farmed salmon, which was corroborated by temporal change in linkage-disequilibrium. Moreover, SNP genome scans identified a temporal decrease in candidate loci potentially under directional selection. Of particular interest was a SNP previously shown to be strongly associated with an important quantitative trait locus for parr mark number, which retained its genetic distinctiveness between farmed and wild fish longer than other outliers. Overall, these results indicate that farmed escapees have introgressed with wild Magaguadavic salmon resulting in significant alteration of the genetic integrity of the native population, including possible loss of adaptation to wild conditions.