Searching for genetic evidence of demographic decline in an arctic seabird: beware of overlapping generations.

Genetic data are useful for detecting sudden population declines in species that are difficult to study in the field. Yet this indirect approach has its own drawbacks, including population structure, mutation patterns, and generation overlap. The ivory gull (Pagophila eburnea), a long-lived Arctic seabird, is currently suffering from rapid alteration of its primary habitat (i.e., sea ice), and dramatic climatic events affecting reproduction and recruitment. However, ivory gulls live in remote areas, and it is difficult to assess the population trend of the species across its distribution. Here we present complementary microsatellite- and SNP-based genetic analyses to test a recent bottleneck genetic signal in ivory gulls over a large portion of their distribution. With attention to the potential effects of population structure, mutation patterns, and sample size, we found no significant signatures of population decline worldwide. At a finer scale, we found a significant bottleneck signal at one location in Canada. These results were compared with predictions from simulations showing how generation time and generation overlap can delay and reduce the bottleneck microsatellite heterozygosity excess signal. The consistency of the results obtained with independent methods strongly indicates that the species shows no genetic evidence of an overall decline in population size. However, drawing conclusions related to the species' population trends will require a better understanding of the effect of age structure in long-lived species. In addition, estimates of the effective global population size of ivory gulls were surprisingly low (~1000 ind.), suggesting that the evolutionary potential of the species is not assured.

Agricultural intensification alters marbled newt genetic diversity and gene flow through density and dispersal reduction.

Recent agricultural intensification threatens global biodiversity with amphibians being one of the most impacted groups. Because of their biphasic life cycle, amphibians are particularly vulnerable to habitat loss and fragmentation that often result in small, isolated populations and loss of genetic diversity. Here, we studied how landscape heterogeneity affects genetic diversity, gene flow and demographic parameters in the marbled newt, Triturus marmoratus, over a hedgerow network landscape in Western France. While the northern part of the study area consists of preserved hedged farmland, the southern part was more profoundly converted for intensive arable crops production after WWII. Based on 67 sampled ponds and 10 microsatellite loci, we characterized regional population genetic structure and evaluated the correlation between landscape variables and (i) local genetic diversity using mixed models and (ii) genetic distance using multiple regression methods and commonality analysis. We identified a single genetic population characterized by a spatially heterogeneous isolation-by-distance pattern. Pond density in the surrounding landscape positively affected local genetic diversity while arable crop land cover negatively affected gene flow and connectivity. We used demographic inferences to quantitatively assess differences in effective population density and dispersal between the contrasted landscapes characterizing the northern and southern parts of the study area. Altogether, results suggest recent land conversion affected T. marmoratus through reduction in both effective population density and dispersal due to habitat loss and reduced connectivity.

GSpace: an exact coalescence simulator of recombining genomes under isolation by distance.

MOTIVATION: Simulation-based inference can bypass the limitations of statistical methods based on analytical approximations, but software allowing simulation of structured population genetic data without the classical n-coalescent approximations (such as those following from assuming large population size) are scarce or slow. RESULTS: We present GSpace, a simulator for genomic data, based on a generation-by-generation coalescence algorithm taking into account small population size, recombination and isolation by distance. AVAILABILITY AND IMPLEMENTATION: Freely available at site web INRAe (http://www1.montpellier.inra.fr/CBGP/software/gspace/download.html).

Sampling schemes and drift can bias admixture proportions inferred by structure.

The interbreeding of individuals coming from genetically differentiated but incompletely isolated populations can lead to the formation of admixed populations, having important implications in ecology and evolution. In this simulation study, we evaluate how individual admixture proportions estimated by the software structure are quantitatively affected by different factors. Using various scenarios of admixture between two diverging populations, we found that unbalanced sampling from parental populations may seriously bias the inferred admixture proportions; moreover, proportionally large samples from the admixed population can also decrease the accuracy and precision of the inferences. As expected, weak differentiation between parental populations and drift after the admixture event strongly increase the biases caused by uneven sampling. We also show that admixture proportions are generally more biased when parental populations unequally contributed to the admixed population. Finally, with few exceptions, using a large number of markers reduces those biases, but using alternative priors for individual ancestry or the uncorrelated allele model only marginally affect the inference of admixture in most situations. We conclude that unbalanced sampling may cause important biases in the admixture proportions estimated by structure, especially when a small number of markers are used, and those biases can be worsened by the effect of drift and unequal genetic contribution of parental populations. Empirical studies should thus be careful with their sampling design and consider historical characteristics when using this software to estimate the ancestry of individuals from admixed populations.

Morphology and multilocus phylogeny of the Spiny-footed Lizard (Acanthodactylus erythrurus) complex reveal two new mountain species from the Moroccan Atlas.

We generated an extensive morphological and multilocus molecular dataset to investigate the taxonomy of Acanthodactylus erythrurus, a widespread species across the Mediterranean and semiarid habitats of the Iberian Peninsula and the Maghreb. Our integrated analyses revealed the existence of at least five basal lineages: (i) an Ibero-Moroccan clade widespread across Morocco and the Iberian Peninsula, from sea-level up to a maximal known altitude of 1,930 m, (ii) an Algero-Tunisian clade, distributed in coastal and inland areas of eastern Algeria and Tunisia, (iii) a Central Algerian clade, formed by two inland populations located in central Algeria (1,000-1,500 m a.s.l.), (iv) a western High Atlas clade including two montane populations from Jbel Siroua and Tizi n'Tichka (at 2,320 m and 2,176 m a.s.l., respectively) and (v) an eastern High Atlas clade, including at least two montane populations from Isli and Tislit (both localities around 2,275 m a.s.l.). An integrated species delimitation approach combining molecular and multivariate morphological analyses demonstrated complete reproductive isolation and hence speciation between the Ibero-Moroccan clade and the eastern High Atlas clade in their contact zone. The divergence between all five clades is broadly similar, supporting the existence of at least five species in the Acanthodactylus erythrurus complex. In the present work we describe the two well-differentiated endemic species from the Moroccan Atlas for which no name is available: Acanthodactylus lacrymae sp. nov. from Isli and Tislit and A. montanus sp. nov. from Jbel Siroua and Tizi n'Tichka. Further work will be needed to fully resolve the taxonomy of this species complex.

Integrating population genetics to define conservation units from the core to the edge of Rhinolophus ferrumequinum western range.

The greater horseshoe bat (Rhinolophus ferrumequinum) is among the most widespread bat species in Europe but it has experienced severe declines, especially in Northern Europe. This species is listed Near Threatened in the European IUCN Red List of Threatened Animals, and it is considered to be highly sensitive to human activities and particularly to habitat fragmentation. Therefore, understanding the population boundaries and demographic history of populations of this species is of primary importance to assess relevant conservation strategies. In this study, we used 17 microsatellite markers to assess the genetic diversity, the genetic structure, and the demographic history of R. ferrumequinum colonies in the western part of its distribution. We identified one large population showing high levels of genetic diversity and large population size. Lower estimates were found in England and northern France. Analyses of clustering and isolation by distance suggested that the Channel and the Mediterranean seas could impede R. ferrumequinum gene flow. These results provide important information to improve the delineation of R. ferrumequinum management units. We suggest that a large management unit corresponding to the population ranging from Spanish Basque Country to northern France must be considered. Particular attention should be given to mating territories as they seem to play a key role in maintaining high levels of genetic mixing between colonies. Smaller management units corresponding to English and northern France colonies must also be implemented. These insular or peripheral colonies could be at higher risk of extinction in the near future.

Population genetics and historical demographic inferences of the blue crab Callinectes sapidus in the US based on microsatellites.

The native range of the blue crab Callinectes sapidus spans Nova Scotia to northern Argentina. In the US, it constitutes a keystone species in estuarine habitats of the Atlantic coast and Gulf of Mexico (GOM), serving as both predator and prey to other species, and also has historically represented a multi-billion dollar fishery. Knowledge relevant to effective management and monitoring of this ecologically and economically important species, such as levels of population genetic differentiation and genetic diversity, is necessary. Although several population genetics studies have attempted to address these questions in one or more parts of its distribution, conflicting results and potential problems with the markers used, as well as other issues, have obscured our understanding on them. In this study, we examined large-scale genetic connectivity of the blue crab in the US, using 16 microsatellites, and genotyped individuals from Chesapeake Bay, in the US Atlantic, and from nine localities along the US GOM coast. Consistent with the high long-distance dispersal potential of this species, very low levels of genetic differentiation were detected for the blue crab among the ten US localities examined, suggesting it constitutes a large panmictic population within this region. Estimations of genetic diversity for the blue crab appear to be high in the US, and provide a baseline for monitoring temporal changes in this species. Demographic analyses indicate a recent range expansion of the US population, probably during the Holocene. In addition, capitalizing on published microsatellite data from southern Brazil, our analyses detected high genetic differentiation between localities in the US and Brazil. These results point to the need for examination of genetic diversity and differentiation along the area spanning the US to southern Brazil.

Limited genetic structure and demographic expansion of the Brassicogethes aeneus populations in France and in Europe.

BACKGROUND: The pollen beetle, Brassicogethes aeneus (Fabricius, 1775), is one of the most significant pests of oilseed rape. To shed light on past and current pollen beetle demography (dispersal, population size), 12 microsatellite markers were developed, and population genetic diversity and structure were analysed at different spatial scales in France and in Europe from 433 individuals collected in 18 winter oilseed rape fields. RESULTS: Genetic differentiation among the population samples was low but was significant between the Estonian sample and the rest of Europe. Isolation by distance was significant only at the European scale. Genetic variability was similar among the 18 population samples. Demographic inferences suggested a recent expansion of B. aeneus population size over Europe, possibly corresponding to an increase in oilseed rape crop area during past decades. CONCLUSION: Current population size and dispersal are not straightforward to estimate from the distribution of genetic variability in B. aeneus over Europe because of the complexity of the demographic history of this pest. Nevertheless, because gene flow was important enough to prevent strong genetic differentiation at large geographical scales, the management of pollen beetle populations should likely be thought of at a continental Europe level. © 2018 Society of Chemical Industry.

Measuring Genetic Differentiation from Pool-seq Data.

The advent of high throughput sequencing and genotyping technologies enables the comparison of patterns of polymorphisms at a very large number of markers. While the characterization of genetic structure from individual sequencing data remains expensive for many nonmodel species, it has been shown that sequencing pools of individual DNAs (Pool-seq) represents an attractive and cost-effective alternative. However, analyzing sequence read counts from a DNA pool instead of individual genotypes raises statistical challenges in deriving correct estimates of genetic differentiation. In this article, we provide a method-of-moments estimator of [Formula: see text] for Pool-seq data, based on an analysis-of-variance framework. We show, by means of simulations, that this new estimator is unbiased and outperforms previously proposed estimators. We evaluate the robustness of our estimator to model misspecification, such as sequencing errors and uneven contributions of individual DNAs to the pools. Finally, by reanalyzing published Pool-seq data of different ecotypes of the prickly sculpin Cottus asper, we show how the use of an unbiased [Formula: see text] estimator may question the interpretation of population structure inferred from previous analyses.

Genetic diversity and structure related to expansion history and habitat isolation: stone marten populating rural-urban habitats.

BACKGROUND: Population genetic diversity and structure are determined by past and current evolutionary processes, among which spatially limited dispersal, genetic drift, and shifts in species distribution boundaries have major effects. In most wildlife species, environmental modifications by humans often lead to contraction of species' ranges and/or limit their dispersal by acting as environmental barriers. However, in species well adapted to anthropogenic habitat or open landscapes, human induced environmental changes may facilitate dispersal and range expansions. In this study, we analysed whether isolation by distance and deforestation, among other environmental features, promotes or restricts dispersal and expansion in stone marten (Martes foina) populations. RESULTS: We genotyped 298 martens from eight sites at twenty-two microsatellite loci to characterize the genetic variability, population structure and demographic history of stone martens in Poland. At the landscape scale, limited genetic differentiation between sites in a mosaic of urban, rural and forest habitats was mostly influenced by isolation by distance. Statistical clustering and multivariate analyses showed weak genetic structuring with two to four clusters and a high rate of gene flow between them. Stronger genetic differentiation was detected for one stone marten population (NE1) located inside a large forest complex. Genetic differentiation between this site and all others was 20% higher than between other sites separated by similar distances. The genetic uniqueness index of NE1 was also twofold higher than in other sites. Past demographic history analyses showed recent expansion of this species in north-eastern Poland. A decrease in genetic diversity from south to north, and MIGRAINE analyses indicated the direction of expansion of stone marten. CONCLUSIONS: Our results showed that two processes, changes in species distribution boundaries and limited dispersal associated with landscape barriers, affect genetic diversity and structure in stone marten. Analysis of local barriers that reduced dispersal and large scale analyses of genetic structure and demographic history highlight the importance of isolation by distance and forest cover for the past colonization of central Europe by stone marten. This confirmed the hypothesis that human-landscape changes (deforestation) accelerated stone marten expansion, to which climate warming probably has also been contributing over the last few decades.

Demographic inference through approximate-Bayesian-computation skyline plots.

The skyline plot is a graphical representation of historical effective population sizes as a function of time. Past population sizes for these plots are estimated from genetic data, without a priori assumptions on the mathematical function defining the shape of the demographic trajectory. Because of this flexibility in shape, skyline plots can, in principle, provide realistic descriptions of the complex demographic scenarios that occur in natural populations. Currently, demographic estimates needed for skyline plots are estimated using coalescent samplers or a composite likelihood approach. Here, we provide a way to estimate historical effective population sizes using an Approximate Bayesian Computation (ABC) framework. We assess its performance using simulated and actual microsatellite datasets. Our method correctly retrieves the signal of contracting, constant and expanding populations, although the graphical shape of the plot is not always an accurate representation of the true demographic trajectory, particularly for recent changes in size and contracting populations. Because of the flexibility of ABC, similar approaches can be extended to other types of data, to multiple populations, or to other parameters that can change through time, such as the migration rate.

A reassessment of explanations for discordant introgressions of mitochondrial and nuclear genomes.

Hybridization is increasingly recognized as a significant evolutionary process, in particular because it can lead to introgression of genes from one species to another. A striking pattern of discordance in the amount of introgression between mitochondrial and nuclear markers exists such that substantial mitochondrial introgression is often found in combination with no or little nuclear introgression. Multiple mechanisms have been proposed to explain this discordance, including positive selection for introgressing mitochondrial variants, several types of sex-biases, drift, negative selection against introgression in the nuclear genome, and spatial expansion. Most of these hypotheses are verbal, and have not been quantitatively evaluated so far. We use individual-based, multilocus, computer simulations of secondary contact under a wide range of demographic and genetic scenarios to evaluate the ability of the different mechanisms to produce discordant introgression. Sex-biases and spatial expansions fail to produce substantial mito-nuclear discordance. Drift and nuclear selection can produce strong discordance, but only under a limited range of conditions. In contrast, selection on the mitochondrial genome produces strong discordance, particularly when dispersal rates are low. However, commonly used statistical tests have little power to detect this selection. Altogether, these results dismiss several popular hypotheses, and provide support for adaptive mitochondrial introgression.

Spatial heterogeneity in landscape structure influences dispersal and genetic structure: empirical evidence from a grasshopper in an agricultural landscape.

Dispersal may be strongly influenced by landscape and habitat characteristics that could either enhance or restrict movements of organisms. Therefore, spatial heterogeneity in landscape structure could influence gene flow and the spatial structure of populations. In the past decades, agricultural intensification has led to the reduction in grassland surfaces, their fragmentation and intensification. As these changes are not homogeneously distributed in landscapes, they have resulted in spatial heterogeneity with generally less intensified hedged farmland areas remaining alongside streams and rivers. In this study, we assessed spatial pattern of abundance and population genetic structure of a flightless grasshopper species, Pezotettix giornae, based on the surveys of 363 grasslands in a 430-km² agricultural landscape of western France. Data were analysed using geostatistics and landscape genetics based on microsatellites markers and computer simulations. Results suggested that small-scale intense dispersal allows this species to survive in intensive agricultural landscapes. A complex spatial genetic structure related to landscape and habitat characteristics was also detected. Two P. giornae genetic clusters bisected by a linear hedged farmland were inferred from clustering analyses. This linear hedged farmland was characterized by high hedgerow and grassland density as well as higher grassland temporal stability that were suspected to slow down dispersal. Computer simulations demonstrated that a linear-shaped landscape feature limiting dispersal could be detected as a barrier to gene flow and generate the observed genetic pattern. This study illustrates the relevance of using computer simulations to test hypotheses in landscape genetics studies.

Congruent signals of population history but radically different patterns of genetic diversity between mitochondrial and nuclear markers in a mountain lizard.

Historical factors, current population size, population connectivity and selective processes at linked loci contribute to shaping contemporary patterns of neutral genetic diversity. It is now widely acknowledged that nuclear and mitochondrial markers react differently to current demography as well as to past history, so the use of both types of markers is often advocated to gain insight on both historical and contemporary processes. We used 12 microsatellite loci genotyped in 13 populations of a mountain lizard (Iberolacerta bonnali) to test whether the historical scenario favoured by a previous mitochondrial study was also supported by nuclear markers and thereby evaluated the consequences of postglacial range movements on nuclear diversity. Congruent signals of recent history were revealed by nuclear and mitochondrial markers using an Approximate Bayesian computation approach, but contemporary patterns of mtDNA and nuclear DNA diversity were radically different. Although dispersal in this species is probably highly restricted at all spatial scales, colonization abilities have been historically good, suggesting capability for reestablishment of locally extinct populations except in fully disconnected habitats.

Blacktip reef sharks, Carcharhinus melanopterus, have high genetic structure and varying demographic histories in their Indo-Pacific range.

For free-swimming marine species like sharks, only population genetics and demographic history analyses can be used to assess population health/status as baseline population numbers are usually unknown. We investigated the population genetics of blacktip reef sharks, Carcharhinus melanopterus; one of the most abundant reef-associated sharks and the apex predator of many shallow water reefs of the Indian and Pacific Oceans. Our sampling includes 4 widely separated locations in the Indo-Pacific and 11 islands in French Polynesia with different levels of coastal development. Four-teen microsatellite loci were analysed for samples from all locations and two mitochondrial DNA fragments, the control region and cytochrome b, were examined for 10 locations. For microsatellites, genetic diversity is higher for the locations in the large open systems of the Red Sea and Australia than for the fragmented habitat of the smaller islands of French Polynesia. Strong significant structure was found for distant locations with FST values as high as ~0.3, and a smaller but still significant structure is found within French Polynesia. Both mitochondrial genes show only a few mutations across the sequences with a dominant shared haplotype in French Polynesia and New Caledonia suggesting a common lineage different to that of East Australia. Demographic history analyses indicate population expansions in the Red Sea and Australia that may coincide with sea level changes after climatic events. Expansions and flat signals are indicated for French Polynesia as well as a significant recent bottleneck for Moorea, the most human-impacted lagoon of the locations in French Polynesia.

Maximum-likelihood inference of population size contractions from microsatellite data.

Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method (MSVAR) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package.

Genetic structure of populations of whale sharks among ocean basins and evidence for their historic rise and recent decline.

This study presents genetic evidence that whale sharks, Rhincodon typus, are comprised of at least two populations that rarely mix and is the first to document a population expansion. Relatively high genetic structure is found when comparing sharks from the Gulf of Mexico with sharks from the Indo-Pacific. If mixing occurs between the Indian and Atlantic Oceans, it is not sufficient to counter genetic drift. This suggests whale sharks are not all part of a single global metapopulation. The significant population expansion we found was indicated by both microsatellite and mitochondrial DNA. The expansion may have happened during the Holocene, when tropical species could expand their range due to sea-level rise, eliminating dispersal barriers and increasing plankton productivity. However, the historic trend of population increase may have reversed recently. Declines in genetic diversity are found for 6 consecutive years at Ningaloo Reef in Australia. The declines in genetic diversity being seen now in Australia may be due to commercial-scale harvesting of whale sharks and collision with boats in past decades in other countries in the Indo-Pacific. The study findings have implications for models of population connectivity for whale sharks and advocate for continued focus on effective protection of the world's largest fish at multiple spatial scales.

DIYABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data.

MOTIVATION: DIYABC is a software package for a comprehensive analysis of population history using approximate Bayesian computation on DNA polymorphism data. Version 2.0 implements a number of new features and analytical methods. It allows (i) the analysis of single nucleotide polymorphism data at large number of loci, apart from microsatellite and DNA sequence data, (ii) efficient Bayesian model choice using linear discriminant analysis on summary statistics and (iii) the serial launching of multiple post-processing analyses. DIYABC v2.0 also includes a user-friendly graphical interface with various new options. It can be run on three operating systems: GNU/Linux, Microsoft Windows and Apple Os X. AVAILABILITY: Freely available with a detailed notice document and example projects to academic users at http://www1.montpellier.inra.fr/CBGP/diyabc CONTACT: estoup@supagro.inra.fr Supplementary information: Supplementary data are available at Bioinformatics online.

The impact of human conflict on the genetics of Mastomys natalensis and Lassa virus in West Africa.

Environmental changes have been shown to play an important role in the emergence of new human diseases of zoonotic origin. The contribution of social factors to their spread, especially conflicts followed by mass movement of populations, has not been extensively investigated. Here we reveal the effects of civil war on the phylogeography of a zoonotic emerging infectious disease by concomitantly studying the population structure, evolution and demography of Lassa virus and its natural reservoir, the rodent Mastomys natalensis, in Guinea, West Africa. Analysis of nucleoprotein gene sequences enabled us to reconstruct the evolutionary history of Lassa virus, which appeared 750 to 900 years ago in Nigeria and only recently spread across western Africa (170 years ago). Bayesian demographic inferences revealed that both the host and the virus populations have gone recently through severe genetic bottlenecks. The timing of these events matches civil war-related mass movements of refugees and accompanying environmental degradation. Forest and habitat destruction and human predation of the natural reservoir are likely explanations for the sharp decline observed in the rodent populations, the consequent virus population decline, and the coincident increased incidence of Lassa fever in these regions. Interestingly, we were also able to detect a similar pattern in Nigeria coinciding with the Biafra war. Our findings show that anthropogenic factors may profoundly impact the population genetics of a virus and its reservoir within the context of an emerging infectious disease.