Spatially explicit summary statistics for historical population genetic inference

Landscape Genetics for Brazilian Equines

Optimization of DNA collection for National gene bank and conservation programs requires information on spatial and genetic distribution of animals countrywide. The relationship between genetic and geographic distances were examined in 8 Brazilian horse breeds (Baixadeiro, Crioulo, Campeiro, Lavradeiro, Marajoara, Mangalarga Marchador, Pantaneiro and Puruca) using Single Nucleotide Polymorphism markers and collection point locations. Mantel correlations, Genetic Landscape Shape Interpolation, Allelic Aggregation Index Analyses and Spatial autocorrelation tests indicated a nonrandom distribution of horses throughout the country. Minimum collection distances for the national Gene Bank should be 530km, with clear divisions seen in genetic structure of horse populations in both North/South and East/West directions. Comparing Pantaneiro and North/Northeastern breeds, physical distance is not necessarily the defining factor for genetic differentiation. This should be considered when sampling these local breeds. These data can help optimise GenBank collection routines and conservation strategies for these breeds.

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

Can secondary contact following range expansion be distinguished from barriers to gene flow?

Secondary contact is the reestablishment of gene flow between sister populations that have diverged. For instance, at the end of the Quaternary glaciations in Europe, secondary contact occurred during the northward expansion of the populations which had found refugia in the southern peninsulas. With the advent of multi-locus markers, secondary contact can be investigated using various molecular signatures including gradients of allele frequency, admixture clines, and local increase of genetic differentiation. We use coalescent simulations to investigate if molecular data provide enough information to distinguish between secondary contact following range expansion and an alternative evolutionary scenario consisting of a barrier to gene flow in an isolation-by-distance model. We find that an excess of linkage disequilibrium and of genetic diversity at the suture zone is a unique signature of secondary contact. We also find that the directionality index ψ, which was proposed to study range expansion, is informative to distinguish between the two hypotheses. However, although evidence for secondary contact is usually conveyed by statistics related to admixture coefficients, we find that they can be confounded by isolation-by-distance. We recommend to account for the spatial repartition of individuals when investigating secondary contact in order to better reflect the complex spatio-temporal evolution of populations and species.

Inferring the geographic origin of a range expansion: Latitudinal and longitudinal coordinates inferred from genomic data in an ABC framework with the program x-origin

Climatic or environmental change is not only driving distributional shifts in species today, but it has also caused distributions to expand and contract in the past. Inferences about the geographic locations of past populations especially regions that served as refugia (i.e., source populations) and migratory routes are a challenging endeavour. Refugial areas may be evidenced from fossil records or regions of temporal stability inferred from ecological niche models. Genomic data offer an alternative and broadly applicable source of information about the locality of refugial areas, especially relative to fossil data, which are either unavailable or incomplete for most species. Here, we present a pipeline we developed (called x-origin) for statistically inferring the geographic origin of range expansion using a spatially explicit coalescent model and an approximate Bayesian computation testing framework. In addition to assessing the probability of specific latitudinal and longitudinal coordinates of refugial or source populations, such inferences can also be made accounting for the effects of temporal and spatial environmental heterogeneity, which may impact migration routes. We demonstrate x-origin with an analysis of genomic data collected in the Collared pika that underwent postglacial expansion across Alaska, as well as present an assessment of its accuracy under a known model of expansion to validate the approach.