Spatial inference of admixture proportions and secondary contact zones

Guns, germs and dogs: On the origin of Leishmania chagasi

The evolutionary history of Leishmania chagasi, the aetiological agent of visceral leishmaniasis in South America has been widely debated. This study addresses the problem of the origin of L. chagasi, its timing and demography with fast evolving genetic markers, a suite of Bayesian clustering algorithms and coalescent modelling. Here, using 14 microsatellite markers, 450 strains from the Leishmania donovani complex, we show that the vast majority of the Central and South American L. chagasi were nested within the Portuguese Leishmania infantum clade. Moreover, L. chagasi allelic richness was half that of their Old World counterparts. The bottleneck signature was estimated to be about 500 years old and the settlement of L. chagasi in the New World, probably via infected dogs, was accompanied by a thousand-fold population decrease. Visceral leishmaniasis, lethal if untreated, is therefore one more disease that the Conquistadores brought to the New World.

Lineage admixture during postglacial range expansion is responsible for the increased gene diversity of Kalopanax septemlobus in a recently colonised territory

We aimed to reveal the effects of range expansion and subsequent lineage admixture from separated glacial refugia on genetic diversity of Kalopanax septemlobus in Japan, by combining nuclear microsatellite data and ecological niche modelling. Allelic richness and gene diversity were compared at the population and regional level. We also statistically examined these indices as a function of population accessibility to the last glacial maximum (LGM) palaeodistribution reconstructed by ecological niche modelling to test a simple range expansion scenario from glacial refugia. Genetic diversity was highest in the populations of southern Japan and gradually decreased towards the north. However, an additional centre of genetic diversity, when measured as gene diversity, was found in northern Honshu Island, where distinct lineages were shown to be in contact. Positive effects of population accessibility to the LGM range were detected in both diversity indices at different spatial scales. The combined data support independent postglacial range expansions towards the north from the edge populations on the exposed coastal shelf of Pacific and Sea of Japan in northern Honshu during the LGM, which subsequently resulted in markedly low genetic diversity in the northernmost extant range, Hokkaido. The regional increase in gene diversity in northern Honshu is likely to be the result of postglacial lineage admixture. Relative difference in the spatial scales best relating population genetic diversity with the LGM distribution can be explained by a higher rate of allelic richness diversity loss during range expansions and stronger effects of lineage admixture on gene diversity.

Genetic signature of population fragmentation varies with mobility in seven bird species of a fragmented Kenyan cloud forest

Habitat fragmentation can restrict geneflow, reduce neighbourhood effective population size, and increase genetic drift and inbreeding in small, isolated habitat remnants. The extent to which habitat fragmentation leads to population fragmentation, however, differs among landscapes and taxa. Commonly, researchers use information on the current status of a species to predict population effects of habitat fragmentation. Such methods, however, do not convey information on species-specific responses to fragmentation. Here, we compare levels of past population differentiation, estimated from microsatellite genotypes, with contemporary dispersal rates, estimated from multi-strata capture-recapture models, to infer changes in mobility over time in seven sympatric, forest-dependent bird species of a Kenyan cloud forest archipelago. Overall, populations of sedentary species were more strongly differentiated and clustered compared to those of vagile ones, while geographical patterning suggested an important role of landscape structure in shaping genetic variation. However, five of seven species with broadly similar levels of genetic differentiation nevertheless differed substantially in their current dispersal rates. We conclude that post-fragmentation levels of vagility, without reference to past population connectivity, may not be the best predictor of how forest fragmentation affects the life history of forest-dependent species. As effective conservation strategies often hinge on accurate prediction of shifts in ecological and genetic relationships among populations, conservation practices based solely upon current population abundances or movements may, in the long term, prove to be inadequate.

Predictions of native American population structure using linguistic covariates in a hidden regression framework

BACKGROUND

The mainland of the Americas is home to a remarkable diversity of languages, and the relationships between genes and languages have attracted considerable attention in the past. Here we investigate to which extent geography and languages can predict the genetic structure of Native American populations.

METHODOLOGY/PRINCIPAL FINDINGS

Our approach is based on a Bayesian latent cluster regression model in which cluster membership is explained by geographic and linguistic covariates. After correcting for geographic effects, we find that the inclusion of linguistic information improves the prediction of individual membership to genetic clusters. We further compare the predictive power of Greenberg's and The Ethnologue classifications of Amerindian languages. We report that The Ethnologue classification provides a better genetic proxy than Greenberg's classification at the stock and at the group levels. Although high predictive values can be achieved from The Ethnologue classification, we nevertheless emphasize that Choco, Chibchan and Tupi linguistic families do not exhibit a univocal correspondence with genetic clusters.

CONCLUSIONS/SIGNIFICANCE

The Bayesian latent class regression model described here is efficient at predicting population genetic structure using geographic and linguistic information in Native American populations.

Comparison of Bayesian clustering and edge detection methods for inferring boundaries in landscape genetics

Recently, techniques available for identifying clusters of individuals or boundaries between clusters using genetic data from natural populations have expanded rapidly. Consequently, there is a need to evaluate these different techniques. We used spatially-explicit simulation models to compare three spatial Bayesian clustering programs and two edge detection methods. Spatially-structured populations were simulated where a continuous population was subdivided by barriers. We evaluated the ability of each method to correctly identify boundary locations while varying: (i) time after divergence, (ii) strength of isolation by distance, (iii) level of genetic diversity, and (iv) amount of gene flow across barriers. To further evaluate the methods' effectiveness to detect genetic clusters in natural populations, we used previously published data on North American pumas and a European shrub. Our results show that with simulated and empirical data, the Bayesian spatial clustering algorithms outperformed direct edge detection methods. All methods incorrectly detected boundaries in the presence of strong patterns of isolation by distance. Based on this finding, we support the application of Bayesian spatial clustering algorithms for boundary detection in empirical datasets, with necessary tests for the influence of isolation by distance.

Glacial refugia in pathogens: European genetic structure of anther smut pathogens on Silene latifolia and Silene dioica

Climate warming is predicted to increase the frequency of invasions by pathogens and to cause the large-scale redistribution of native host species, with dramatic consequences on the health of domesticated and wild populations of plants and animals. The study of historic range shifts in response to climate change, such as during interglacial cycles, can help in the prediction of the routes and dynamics of infectious diseases during the impending ecosystem changes. Here we studied the population structure in Europe of two Microbotryum species causing anther smut disease on the plants Silene latifolia and Silene dioica. Clustering analyses revealed the existence of genetically distinct groups for the pathogen on S. latifolia, providing a clear-cut example of European phylogeography reflecting recolonization from southern refugia after glaciation. The pathogen genetic structure was congruent with the genetic structure of its host species S. latifolia, suggesting dependence of the migration pathway of the anther smut fungus on its host. The fungus, however, appeared to have persisted in more numerous and smaller refugia than its host and to have experienced fewer events of large-scale dispersal. The anther smut pathogen on S. dioica also showed a strong phylogeographic structure that might be related to more northern glacial refugia. Differences in host ecology probably played a role in these differences in the pathogen population structure. Very high selfing rates were inferred in both fungal species, explaining the low levels of admixture between the genetic clusters. The systems studied here indicate that migration patterns caused by climate change can be expected to include pathogen invasions that follow the redistribution of their host species at continental scales, but also that the recolonization by pathogens is not simply a mirror of their hosts, even for obligate biotrophs, and that the ecology of hosts and pathogen mating systems likely affects recolonization patterns.

Regional differentiation and post-glacial expansion of the Atlantic silverside, Menidia menidia, an annual fish with high dispersal potential

The coastal marine environment of the Northwest Atlantic contains strong environmental gradients that create distinct marine biogeographic provinces by limiting dispersal, recruitment, and survival. This region has also been subjected to numerous Pleistocene glacial cycles, resulting in repeated extirpations and recolonizations in northern populations of marine organisms. In this study, we examined patterns of genetic structure and historical demography in the Atlantic silverside, Menidia menidia, an annual marine fish with high dispersal potential but with well-documented patterns of clinal phenotypic adaptation along the environmental gradients of the Northwest Atlantic. Contrary to previous studies indicating genetic homogeneity that should preclude regional adaptation, results demonstrate subtle but significant (FST = 0.07; P < 0.0001) genetic structure among three phylogeographic regions that partially correspond with biogeographic provinces, suggesting regional limits to gene flow. Tests for non-equilibrium population dynamics and latitudinal patterns in genetic diversity indicate northward population expansion from a single southern refugium following the last glacial maximum, suggesting that phylogeographic and phenotypic patterns have relatively recent origins. The recovery of phylogeographic structure and the partial correspondence of these regions to recognized biogeographic provinces suggest that the environmental gradients that shape biogeographic patterns in the Northwest Atlantic may also limit gene flow in M. menidia, creating phylogeographic structure and contributing to the creation of latitudinal phenotypic clines in this species.

The use of genetics for the management of a recovering population: temporal assessment of migratory peregrine falcons in North America

BACKGROUND

Our ability to monitor populations or species that were once threatened or endangered and in the process of recovery is enhanced by using genetic methods to assess overall population stability and size over time. This can be accomplished most directly by obtaining genetic measures from temporally-spaced samples that reflect the overall stability of the population as given by changes in genetic diversity levels (allelic richness and heterozygosity), degree of population differentiation (F(ST) and D(EST)), and effective population size (N(e)). The primary goal of any recovery effort is to produce a long-term self-sustaining population, and these genetic measures provide a metric by which we can gauge our progress and help make important management decisions.

METHODOLOGY/PRINCIPAL FINDINGS

The peregrine falcon in North America (Falco peregrinus tundrius and anatum) was delisted in 1994 and 1999, respectively, and its abundance will be monitored by the species Recovery Team every three years until 2015. Although the United States Fish and Wildlife Service makes a distinction between tundrius and anatum subspecies, our genetic results based on eleven microsatellite loci suggest limited differentiation that can be attributed to an isolation by distance relationship and warrant no delineation of these two subspecies in its northern latitudinal distribution from Alaska through Canada into Greenland. Using temporal samples collected at Padre Island, Texas during migration (seven temporal time periods between 1985-2007), no significant differences in genetic diversity or significant population differentiation in allele frequencies between time periods were observed and were indistinguishable from those obtained from tundrius/anatum breeding locations throughout their northern distribution. Estimates of harmonic mean N(e) were variable and imprecise, but always greater than 500 when employing multiple temporal genetic methods.

CONCLUSIONS/SIGNIFICANCE

These results, including those from simulations to assess the power of each method to estimate N(e), suggest a stable or growing population, which is consistent with ongoing field-based monitoring surveys. Therefore, historic and continuing efforts to prevent the extinction of the peregrine falcon in North America appear successful with no indication of recent decline, at least from the northern latitude range-wide perspective. The results also further highlight the importance of archiving samples and their use for continual assessment of population recovery and long-term viability.

Habitat distribution influences dispersal and fine-scale genetic population structure of eastern foxsnakes (Mintonius gloydi) across a fragmented landscape

Dispersal is a fundamental attribute of species in nature and shapes population dynamics, evolutionary trajectories and genetic variation across spatial and temporal scales. It is increasingly clear that landscape features have large impacts on dispersal patterns. Thus, understanding how individuals and species move through landscapes is essential for predicting impacts of landscape alterations. Information on dispersal patterns, however, is lacking for many taxa, particularly reptiles. Eastern foxsnakes (Mintoinus gloydi) are marsh and prairie specialists that avoid agricultural fields, but they have persisted across a fragmented region in southwestern Ontario and northern Ohio. Here, we combined habitat suitability modelling with population genetic analyses to infer how foxsnakes disperse through a habitat mosaic of natural and altered landscape features. Boundary regions between the eight genetic clusters, identified through assignment tests, were comprised of low suitability habitat (e.g. agricultural fields). Island populations were grouped into a single genetic cluster, and comparatively low F(ST) values between island and mainland populations suggest open water presents less of a barrier than nonsuitable terrestrial habitat. Isolation by resistance and least-cost path analysis produced similar results with matrices of pairwise individual genetic distance significantly more correlated to matrices of resistance values derived from habitat suitability than models with an undifferentiated landscape. Spatial autocorrelation results matched better with assignment results when incorporating resistance values rather than straight-line distances. All analyses used in our study produced similar results suggesting that habitat degradation limits dispersal for foxsnakes, which has had a strong effect on the genetic population structure across this region.

Spatially explicit Bayesian clustering models in population genetics

This article reviews recent developments in Bayesian algorithms that explicitly include geographical information in the inference of population structure. Current models substantially differ in their prior distributions and background assumptions, falling into two broad categories: models with or without admixture. To aid users of this new generation of spatially explicit programs, we clarify the assumptions underlying the models, and we test these models in situations where their assumptions are not met. We show that models without admixture are not robust to the inclusion of admixed individuals in the sample, thus providing an incorrect assessment of population genetic structure in many cases. In contrast, admixture models are robust to an absence of admixture in the sample. We also give statistical and conceptual reasons why data should be explored using spatially explicit models that include admixture.

Principal component analysis under population genetic models of range expansion and admixture

In a series of highly influential publications, Cavalli-Sforza and colleagues used principal component (PC) analysis to produce maps depicting how human genetic diversity varies across geographic space. Within Europe, the first axis of variation (PC1) was interpreted as evidence for the demic diffusion model of agriculture, in which farmers expanded from the Near East approximately 10,000 years ago and replaced the resident hunter-gatherer populations with little or no interbreeding. These interpretations of the PC maps have been recently questioned as the original results can be reproduced under models of spatially covarying allele frequencies without any expansion. Here, we study PC maps for data simulated under models of range expansion and admixture. Our simulations include a spatially realistic model of Neolithic farmer expansion and assume various levels of interbreeding between farmer and resident hunter-gatherer populations. An important result is that under a broad range of conditions, the gradients in PC1 maps are oriented along a direction perpendicular to the axis of the expansion, rather than along the same axis as the expansion. We propose that this surprising pattern is an outcome of the "allele surfing" phenomenon, which creates sectors of high allele-frequency differentiation that align perpendicular to the direction of the expansion.