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

Post-glacial phylogeography and evolution of a wide-ranging highly-exploited keystone forest tree, eastern white pine (Pinus strobus) in North America: single refugium, multiple routes

Background

Knowledge of the historical distribution and postglacial phylogeography and evolution of a species is important to better understand its current distribution and population structure and potential fate in the future, especially under climate change conditions, and conservation of its genetic resources. We have addressed this issue in a wide-ranging and heavily exploited keystone forest tree species of eastern North America, eastern white pine (Pinus strobus). We examined the range-wide population genetic structure, tested various hypothetical population history and evolutionary scenarios and inferred the location of glacial refugium and post-glacial recolonization routes. Our hypothesis was that eastern white pine survived in a single glacial refugium and expanded through multiple post-glacial recolonization routes.

Results

We studied the range-wide genetic diversity and population structure of 33 eastern white pine populations using 12 nuclear and 3 chloroplast microsatellite DNA markers. We used Approximate Bayesian Computation approach to test various evolutionary scenarios. We observed high levels of genetic diversity, and significant genetic differentiation (F_ST = 0.104) and population structure among eastern white pine populations across its range. A south to north trend of declining genetic diversity existed, consistent with repeated founder effects during post-glaciation migration northwards. We observed broad consensus from nuclear and chloroplast genetic markers supporting the presence of two main post-glacial recolonization routes that originated from a single southern refugium in the mid-Atlantic plain. One route gave rise to populations at the western margin of the species’ range in Minnesota and western Ontario. The second route gave rise to central-eastern populations, which branched into two subgroups: central and eastern. We observed minimal sharing of chloroplast haplotypes between recolonization routes but there was evidence of admixture between the western and west-central populations.

Conclusions

Our study reveals a single southern refugium, two recolonization routes and three genetically distinguishable lineages in eastern white pine that we suggest to be treated as separate Evolutionarily Significant Units. Like many wide-ranging North American species, eastern white pine retains the genetic signatures of post-glacial recolonization and evolution, and its contemporary population genetic structure reflects not just the modern distribution and effects of heavy exploitation but also routes northward from its glacial refugium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0624-1) contains supplementary material, which is available to authorized users.

Phylogeographic breaks within Asian butternuts indicate the existence of a phytogeographic divide in East Asia

East Asia has been hypothesized to be subdivided into two distinct northern and southern areas, separated by a band of dry climate that was far more severe in the early Tertiary but still exists today. However, this biogeographic hypothesis has rarely been tested using a molecular phylogeographic approach. We genotyped 70 populations throughout the distributional range of Asian butternuts (Juglans section Cardiocaryon) using eight chloroplast DNA regions, one single-copy nuclear gene, and 17 nuclear microsatellite loci, supplemented with paleodistribution modeling of the major genetic clades. The genetic data consistently identified two clades, one northern, comprising Juglans mandshurica and Juglans ailantifolia, and one southern, comprising Juglans cathayensis. The two clades diverged through climate-induced vicariance of an ancestral northern range during the mid-Miocene and remained mostly separate thereafter, with geographical isolation of the Japanese Islands and refugial isolation or secondary contacts in the late Pleistocene producing further subdivision within the northern clade. But beyond all that, we also discovered a role of environmental adaptation in maintaining and/or reinforcing the north-south divergence. Asian butternuts offer a strong case for the existence of a biogeographic divide between the northern and southern parts of East Asia during the Neogene and into the Pleistocene.

Low genetic diversity in pygmy blue whales is due to climate-induced diversification rather than anthropogenic impacts

Unusually low genetic diversity can be a warning of an urgent need to mitigate causative anthropogenic activities. However, current low levels of genetic diversity in a population could also be due to natural historical events, including recent evolutionary divergence, or long-term persistence at a small population size. Here, we determine whether the relatively low genetic diversity of pygmy blue whales (Balaenoptera musculus brevicauda) in Australia is due to natural causes or overexploitation. We apply recently developed analytical approaches in the largest genetic dataset ever compiled to study blue whales (297 samples collected after whaling and representing lineages from Australia, Antarctica and Chile). We find that low levels of genetic diversity in Australia are due to a natural founder event from Antarctic blue whales (Balaenoptera musculus intermedia) that occurred around the Last Glacial Maximum, followed by evolutionary divergence. Historical climate change has therefore driven the evolution of blue whales into genetically, phenotypically and behaviourally distinct lineages that will likely be influenced by future climate change.

The eastern migratory caribou: the role of genetic introgression in ecotype evolution

Understanding the evolutionary history of contemporary animal groups is essential for conservation and management of endangered species like caribou (Rangifer tarandus). In central Canada, the ranges of two caribou subspecies (barren-ground/woodland caribou) and two woodland caribou ecotypes (boreal/eastern migratory) overlap. Our objectives were to reconstruct the evolutionary history of the eastern migratory ecotype and to assess the potential role of introgression in ecotype evolution. STRUCTURE analyses identified five higher order groups (i.e. three boreal caribou populations, eastern migratory ecotype and barren-ground). The evolutionary history of the eastern migratory ecotype was best explained by an early genetic introgression from barren-ground into a woodland caribou lineage during the Late Pleistocene and subsequent divergence of the eastern migratory ecotype during the Holocene. These results are consistent with the retreat of the Laurentide ice sheet and the colonization of the Hudson Bay coastal areas subsequent to the establishment of forest tundra vegetation approximately 7000 years ago. This historical reconstruction of the eastern migratory ecotype further supports its current classification as a conservation unit, specifically a Designatable Unit, under Canada's Species at Risk Act. These findings have implications for other sub-specific contact zones for caribou and other North American species in conservation unit delineation.

The Story of a Hitchhiker: Population Genetic Patterns in the Invasive Barnacle Balanus(Amphibalanus) improvisus Darwin 1854

Understanding the ecological and evolutionary forces that determine the genetic structure and spread of invasive species is a key component of invasion biology. The bay barnacle, Balanus improvisus (= Amphibalanus improvisus), is one of the most successful aquatic invaders worldwide, and is characterised by broad environmental tolerance. Although the species can spread through natural larval dispersal, human-mediated transport through (primarily) shipping has almost certainly contributed to the current global distribution of this species. Despite its worldwide distribution, little is known about the phylogeography of this species. Here, we characterize the population genetic structure and model dispersal dynamics of the barnacle B. improvisus, and describe how human-mediated spreading via shipping as well as natural larval dispersal may have contributed to observed genetic variation. We used both mitochondrial DNA (cytochrome c oxidase subunit I: COI) and nuclear microsatellites to characterize the genetic structure in 14 populations of B. improvisus on a global and regional scale (Baltic Sea). Genetic diversity was high in most populations, and many haplotypes were shared among populations on a global scale, indicating that long-distance dispersal (presumably through shipping and other anthropogenic activities) has played an important role in shaping the population genetic structure of this cosmopolitan species. We could not clearly confirm prior claims that B. improvisus originates from the western margins of the Atlantic coasts; although there were indications that Argentina could be part of a native region. In addition to dispersal via shipping, we show that natural larval dispersal may play an important role for further colonisation following initial introduction.

Molecular and iridescent feather reflectance data reveal recent genetic diversification and phenotypic differentiation in a cloud forest hummingbird

The present day distribution and spatial genetic diversity of Mesoamerican biota reflects a long history of responses to habitat change. The hummingbird Lampornis amethystinus is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling across the species range using mitochondrial DNA (mtDNA) sequences and nuclear microsatellites jointly analysed with phenotypic and climatic data, we (1) test whether the fragmented distribution is correlated with main evolutionary lineages, (2) assess body size and plumage color differentiation of populations in geographic isolation, and (3) evaluate a set of divergence scenarios and demographic patterns of the hummingbird populations. Analysis of genetic variation revealed four main groups: blue‐throated populations (Sierra Madre del Sur); two groups of amethyst‐throated populations (Trans‐Mexican Volcanic Belt and Sierra Madre Oriental); and populations east of the Isthmus of Tehuantepec (IT) with males showing an amethyst throat. The most basal split is estimated to have originated in the Pleistocene, 2.39–0.57 million years ago (MYA), and corresponded to groups of populations separated by the IT. However, the estimated recent divergence time between blue‐ and amethyst‐throated populations does not correspond to the 2‐MY needed to be in isolation for substantial plumage divergence, likely because structurally iridescent colors are more malleable than others. Results of species distribution modeling and Approximate Bayesian Computation analysis fit a model of lineage divergence west of the Isthmus after the Last Glacial Maximum (LGM), and that the species’ suitable habitat was disjunct during past and current conditions. These results challenge the generality of the contraction/expansion glacial model to cloud forest‐interior species and urges management of cloud forest, a highly vulnerable ecosystem to climate change and currently facing destruction, to prevent further loss of genetic diversity or extinction.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Mimetic Divergence and the Speciation Continuum in the Mimic Poison Frog Ranitomeya imitator

While divergent ecological adaptation can drive speciation, understanding the factors that facilitate or constrain this process remains a major goal in speciation research. Here, we study two mimetic transition zones in the poison frog Ranitomeya imitator, a species that has undergone a Müllerian mimetic radiation to establish four morphs in Peru. We find that mimetic morphs are strongly phenotypically differentiated, producing geographic clines with varying widths. However, distinct morphs show little neutral genetic divergence, and landscape genetic analyses implicate isolation by distance as the primary determinant of among-population genetic differentiation. Mate choice experiments suggest random mating at the transition zones, although certain allopatric populations show a preference for their own morph. We present evidence that this preference may be mediated by color pattern specifically. These results contrast with an earlier study of a third transition zone, in which a mimetic shift was associated with reproductive isolation. Overall, our results suggest that the three known mimetic transition zones in R. imitator reflect a speciation continuum, which we have characterized at the geographic, phenotypic, behavioral, and genetic levels. We discuss possible explanations for variable progress toward speciation, suggesting that multifarious selection on both mimetic color pattern and body size may be responsible for generating reproductive isolation.

Historical population size of the threatened New Zealand sea lion Phocarctos hookeri

Marine mammal species were exploited worldwide during periods of commercial sealing in the 18th and 19th centuries. For many of these species, an estimate of the pre-exploitation abundance of the species is lacking, as historical catch records are generally scarce and inaccurate. Genetic estimates of long-term effective population size provide a means to estimate the pre-exploitation abundance. Here, we apply genetic methods to estimate the long-term effective population size of the subantarctic lineage of the New Zealand sea lion (NZ sea lion), Phocarctos hookeri . This species is predominantly restricted to the subantarctic islands, south of mainland New Zealand, following commercial sealing in the 19th century. Today, the population consists of ~9,880 animals and population growth is slow. Auckland Island breeding colonies of NZ sea lion are currently impacted by commercial trawl fisheries via regular sea lion deaths as bycatch. In order to estimate sustainable levels of bycatch, an estimate of the population’s carrying capacity ( K ) is required. We apply the genetically estimated long-term effective population size of NZ sea lions as a proxy for the estimated historical carrying capacity of the subantarctic population. The historical abundance of subantarctic NZ sea lions was significantly higher than the target values of K employed by the contemporary management. The current management strategy may allow unsustainable bycatch levels, thereby limiting the recovery of the NZ sea lion population toward historical carrying capacity.

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter- and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.

Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations

Unravelling variation among taxonomic orders regarding the rate of evolution in microsatellites is crucial for evolutionary biology and population genetics research. The mean mutation rate of microsatellites tends to be lower in arthropods than in vertebrates, but data are scarce and mostly concern accumulation of mutations in model species. Based on parent-offspring segregations and a hierarchical Bayesian model, the mean rate of mutation in the orthopteran insect Schistocerca gregaria was estimated at 2.1e(-4) per generation per untranscribed dinucleotide locus. This is close to vertebrate estimates and one order of magnitude higher than estimates from species of other arthropod orders, such as Drosophila melanogaster and Daphnia pulex. We also found evidence of a directional bias towards expansions even for long alleles and exceptionally large ranges of allele sizes. Finally, at transcribed microsatellites, the mean rate of mutation was half the rate found at untranscribed loci and the mutational model deviated from that usually considered, with most mutations involving multistep changes that avoid disrupting the reading frame. Our direct estimates of mutation rate were discussed in the light of peculiar biological and genomic features of S. gregaria, including specificities in mismatch repair and the dependence of its activity to allele length. Shedding new light on the mutational dynamics of grasshopper microsatellites is of critical importance for a number of research fields. As an illustration, we showed how our findings improve microsatellite application in population genetics, by obtaining a more precise estimation of S. gregaria effective population size from a published data set based on the same microsatellites.

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep-sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep-sea ecosystems is scarce, especially for deep-sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep-sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497-bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ(ST) = -0.0091, F(ST) = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ(ST) = 0.0744, P < 0.01; F(ST) = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ(ST) = 0.3808, F(ST) = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.

Reliable ABC model choice via random forests

MOTIVATION

Approximate Bayesian computation (ABC) methods provide an elaborate approach to Bayesian inference on complex models, including model choice. Both theoretical arguments and simulation experiments indicate, however, that model posterior probabilities may be poorly evaluated by standard ABC techniques.

RESULTS

We propose a novel approach based on a machine learning tool named random forests (RF) to conduct selection among the highly complex models covered by ABC algorithms. We thus modify the way Bayesian model selection is both understood and operated, in that we rephrase the inferential goal as a classification problem, first predicting the model that best fits the data with RF and postponing the approximation of the posterior probability of the selected model for a second stage also relying on RF. Compared with earlier implementations of ABC model choice, the ABC RF approach offers several potential improvements: (i) it often has a larger discriminative power among the competing models, (ii) it is more robust against the number and choice of statistics summarizing the data, (iii) the computing effort is drastically reduced (with a gain in computation efficiency of at least 50) and (iv) it includes an approximation of the posterior probability of the selected model. The call to RF will undoubtedly extend the range of size of datasets and complexity of models that ABC can handle. We illustrate the power of this novel methodology by analyzing controlled experiments as well as genuine population genetics datasets.

AVAILABILITY AND IMPLEMENTATION

The proposed methodology is implemented in the R package abcrf available on the CRAN.

CONTACT

jean-michel.marin@umontpellier.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Population genetic structure of the predatory, social wasp Vespula pensylvanica in its native and invasive range

Invasive species cause extensive damage to their introduced ranges. Ocean archipelagos are particularly vulnerable to invasive taxa. In this study, we used polymorphic microsatellite markers to investigate the genetic structure of the social wasp Vespula pensylvanica in its native range of North America and its introduced range in the archipelago of Hawaii. Our goal was to gain a better understanding of the invasion dynamics of social species and the processes affecting biological invasions. We found that V. pensylvanica showed no significant genetic isolation by distance and little genetic structure over a span of 2000 km in its native range. This result suggests that V. pensylvanica can successfully disperse across large distances either through natural- or human-mediated mechanisms. In contrast to the genetic patterns observed in the native range, we found substantial genetic structure in the invasive V. pensylvanica range in Hawaii. The strong patterns of genetic differentiation within and between the Hawaiian Islands may reflect the effects of geographic barriers and invasion history on gene flow. We also found some evidence for gene flow between the different islands of Hawaii which was likely mediated through human activity. Overall, this study provides insight on how geographic barriers, invasion history, and human activity can shape population genetic structure of invasive species.

The genomics of ecological vicariance in threespine stickleback fish

Populations occurring in similar habitats and displaying similar phenotypes are increasingly used to explore parallel evolution at the molecular level. This generally ignores the possibility that parallel evolution can be mimicked by the fragmentation of an ancestral population followed by genetic exchange with ecologically different populations. Here we demonstrate such an ecological vicariance scenario in multiple stream populations of threespine stickleback fish divergent from a single adjacent lake population. On the basis of demographic and population genomic analyses, we infer the initial spread of a stream-adapted ancestor followed by the emergence of a lake-adapted population, that selective sweeps have occurred mainly in the lake population, that adaptive lake–stream divergence is maintained in the face of gene flow from the lake into the streams, and that this divergence involves major inversion polymorphisms also important to marine-freshwater stickleback divergence. Overall, our study highlights the need for a robust understanding of the demographic and selective history in evolutionary investigations.

Threespine stickleback fish are adapted to lake and stream habitats in Central Europe. Here, the authors show colonization of a lake basin by a stream-adapted ancestor, followed by the emergence of a lake-adapted population in the face of gene flow across lake–stream boundaries.

Reconstructing the invasion history of Heracleum persicum (Apiaceae) into Europe

Sparse, incomplete and inappropriate historical records of invasive species often hamper invasive species management interventions. Population genetic analyses of invaders might provide a suitable context for the identification of their source populations and possible introduction routes. Here, we describe the population genetics of Heracleum persicum Desf. ex Fisch and trace its route of introduction into Europe. Microsatellite markers revealed a significantly higher genetic diversity of H. persicum in its native range, and the loss of diversity in the introduced range may be attributed to a recent genetic bottleneck. Bayesian cluster analysis on regional levels identified three and two genetic clusters in the native and the introduced ranges, respectively. A global structure analysis revealed two worldwide distinct genetic groups: one primarily in Iran and Denmark, the other primarily in Norway. There were also varying degrees of admixture in England, Sweden, Finland and Latvia. Approximate Bayesian computation indicated two independent introductions of H. persicum from Iran to Europe: the first one in Denmark and the second one in England. Finland was subsequently colonized by English populations. In contrast to the contemporary hypothesis of English origin of Norwegian populations, we found Finland to be a more likely source for Norwegian populations, a scenario supported by higher estimated historical migration from Finland to Norway. Genetic diversity per se is not a primary determinant of invasiveness in H. persicum. Our results indicate that, due to either pre-adaptations or rapid local adaptations, introduced populations may have acquired invasiveness after subsequent introductions, once a suitable environment was encountered.

Ecological opportunities and specializations shaped genetic divergence in a highly mobile marine top predator

Environmental conditions can shape genetic and morphological divergence. Release of new habitats during historical environmental changes was a major driver of evolutionary diversification. Here, forces shaping population structure and ecotype differentiation ('pelagic' and 'coastal') of bottlenose dolphins in the North-east Atlantic were investigated using complementary evolutionary and ecological approaches. Inference of population demographic history using approximate Bayesian computation indicated that coastal populations were likely founded by the Atlantic pelagic population after the Last Glacial Maxima probably as a result of newly available coastal ecological niches. Pelagic dolphins from the Atlantic and the Mediterranean Sea likely diverged during a period of high productivity in the Mediterranean Sea. Genetic differentiation between coastal and pelagic ecotypes may be maintained by niche specializations, as indicated by stable isotope and stomach content analyses, and social behaviour. The two ecotypes were only weakly morphologically segregated in contrast to other parts of the World Ocean. This may be linked to weak contrasts between coastal and pelagic habitats and/or a relatively recent divergence. We suggest that ecological opportunity to specialize is a major driver of genetic and morphological divergence. Combining genetic, ecological and morphological approaches is essential to understanding the population structure of mobile and cryptic species.

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species' evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.

Mitogenomic Meta-Analysis Identifies Two Phases of Migration in the History of Eastern Eurasian Sheep

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637-0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8-11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20-60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5-6.8 ka (lineages A and B: ∼6.4-6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West-East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a "transportation hub" in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Challenging the inbreeding hypothesis in a eusocial mammal: population genetics of the naked mole-rat, Heterocephalus glaber

The role of genetic relatedness in the evolution of eusociality has been the topic of much debate, especially when contrasting eusocial insects with vertebrates displaying reproductive altruism. The naked mole-rat, Heterocephalus glaber, was the first described eusocial mammal. Although this discovery was based on an ecological constraints model of eusocial evolution, early genetic studies reported high levels of relatedness in naked mole-rats, providing a compelling argument that low dispersal rates and consanguineous mating (inbreeding as a mating system) are the driving forces for the evolution of this eusocial species. One caveat to accepting this long-held view is that the original genetic studies were based on limited sampling from the species' geographic distribution. A growing body of evidence supports a contrary view, with the original samples not representative of the species-rather reflecting a single founder event, establishing a small population south of the Athi River. Our study is the first to address these competing hypotheses by examining patterns of molecular variation in colonies sampled from north and south of the Athi and Tana rivers, which based on our results, serve to isolate genetically distinct populations of naked mole-rats. Although colonies south of the Athi River share a single mtDNA haplotype and are fixed at most microsatellite loci, populations north of the Athi River are considerably more variable. Our findings support the position that the low variation observed in naked mole-rat populations south of the Athi River reflects a founder event, rather than a consequence of this species' unusual mating system.

Genetic landscape with sharp discontinuities shaped by complex demographic history in moose ( Alces alces )

The moose ( Alces alces ) is the most intensely managed game species in Fennoscandia; approximately one-third of the population, ca. 160,000 animals, is harvested annually. Despite the species’ biological and socioeconomic importance, there are knowledge gaps with respect to its intraspecific diversity and genetic structure. Recent studies of moose in neighboring countries report 2 genetic groups in Finland, 3 in Norway with one of them suggested to be of ancient origin, and no indications of bottlenecks. To delineate the spatial genetic landscape of the Swedish moose, we used allozyme variability from over 20,000 georeferenced moose collected all over Sweden in combination with 12 microsatellites ( n &gt; 1,200) and mitochondrial DNA (mtDNA) sequences ( n = 44). We combined individual-based and traditional statistical approaches with coalescence-based simulations. The results indicate a complex history with bottlenecks and recent expansions that is consistent with historical records. Swedish moose are separated into 2 major genetic groups, a northern and a southern one, where the southern group is further divided into 3 subgroups. The 2 main subpopulations are moderately differentiated ( FST = 0.1; RST = 0.07) and separated by sharp genetic discontinuities occurring over a relatively narrow transition zone in central Sweden that coincides with a similar, previously reported transition zone in Norway. This differentiation is not reflected in mtDNA variation, where no significant divergence was observed. Together with the FST and RST similarities, this suggests that the 2 major subpopulations in Sweden reflect divergence shaped after the postglacial recolonization of Scandinavia. Neighborhood size assessments indicate that gene flow is relatively restricted with an estimated average dispersal distance of 3.5–11.1 km, and spatial autocorrelograms suggest that genetic similarity decreases almost linearly over space resulting in continuous genetic clines within major subgroups. Management areas largely coincide with genetic clusters, simplifying the integration of genetic information into management.

Higher genetic diversity in recolonized areas than in refugia of Alnus glutinosa triggered by continent-wide lineage admixture

Genetic admixture is supposed to be an important trigger of species expansions because it can create the potential for selection of genotypes suitable for new climatic conditions. Up until now, however, no continent-wide population genetic study has performed a detailed reconstruction of admixture events during natural species expansions. To fill this gap, we analysed the postglacial history of Alnus glutinosa, a keystone species of European swamp habitats, across its entire distribution range using two molecular markers, cpDNA and nuclear microsatellites. CpDNA revealed multiple southern refugia located in the Iberian, Apennine, Balkan and Anatolian Peninsulas, Corsica and North Africa. Analysis of microsatellites variation revealed three main directions of postglacial expansion: (i) from the northern part of the Iberian Peninsula to Western and Central Europe and subsequently to the British Isles, (ii) from the Apennine Peninsula to the Alps and (iii) from the eastern part of the Balkan Peninsula to the Carpathians followed by expansion towards the Northern European plains. This challenges the classical paradigm that most European populations originated from refugial areas in the Carpathians. It has been shown that colonizing lineages have met several times and formed secondary contact zones with unexpectedly high population genetic diversity in Central Europe and Scandinavia. On the contrary, limited genetic admixture in southern refugial areas of A. glutinosa renders rear-edge populations in the Mediterranean region more vulnerable to extinction due to climate change.

Population Genomics of the Euryhaline Teleost Poecilia latipinna

Global climate change and increases in sea levels will affect coastal marine communities. The conservation of these ecologically important areas will be a challenge because of their wide geographic distribution, ecological diversity and species richness. To address this problem, we need to better understand how the genetic variation of the species in these communities is distributed within local populations, among populations and between distant regions. In this study we apply genotyping by sequencing (GBS) and examine 955 SNPs to determine Sailfin molly (Poecilia latipinna) genetic diversity among three geographically close mangrove salt marsh flats in the Florida Keys compared to populations in southern and northern Florida. The questions we are asking are whether there is sufficient genetic variation among isolated estuarine fish within populations and whether there are significant divergences among populations. Additionally, we want to know if GBS approaches agree with previous studies using more traditional molecular approaches. We are able to identify large genetic diversity within each saltmarsh community (π ≈ 36%). Additionally, among the Florida Key populations and the mainland or between southern and northern Florida regions, there are significant differences in allele frequencies seen in population structure and evolutionary relationships among individuals. Surprisingly, even though the cumulative F_ST value using all 955 SNPs within the three Florida Key populations is small, there are 29 loci with significant F_ST values, and 11 of these were outliers suggestive of adaptive divergence. These data suggest that among the salt marsh flats surveyed here, there is significant genetic diversity within each population and small but significant differences among populations. Much of the genetic variation within and among populations found here with GBS is very similar to previous studies using allozymes and microsatellites. However, the meaningful difference between GBS and these previous measures of genetic diversity is the number of loci examined, which allows more precise delineations of population structure as well as facilitates identifying loci with excessive F_ST values that could indicate adaptive divergence.

Low genetic diversity despite multiple introductions of the invasive plant species Impatiens glandulifera in Europe

Background

Invasive species can be a major threat to native biodiversity and the number of invasive plant species is increasing across the globe. Population genetic studies of invasive species can provide key insights into their invasion history and ensuing evolution, but also for their control. Here we genetically characterise populations of Impatiens glandulifera, an invasive plant in Europe that can have a major impact on native plant communities. We compared populations from the species’ native range in Kashmir, India, to those in its invaded range, along a latitudinal gradient in Europe. For comparison, the results from 39 other studies of genetic diversity in invasive species were collated.

Results

Our results suggest that I. glandulifera was established in the wild in Europe at least twice, from an area outside of our Kashmir study area. Our results further revealed that the genetic diversity in invasive populations of I. glandulifera is unusually low compared to native populations, in particular when compared to other invasive species. Genetic drift rather than mutation seems to have played a role in differentiating populations in Europe. We find evidence of limitations to local gene flow after introduction to Europe, but somewhat less restrictions in the native range. I. glandulifera populations with significant inbreeding were only found in the species’ native range and invasive species in general showed no increase in inbreeding upon leaving their native ranges. In Europe we detect cases of migration between distantly located populations. Human activities therefore seem to, at least partially, have facilitated not only introductions, but also further spread of I. glandulifera across Europe.

Conclusions

Although multiple introductions will facilitate the retention of genetic diversity in invasive ranges, widespread invasive species can remain genetically relatively invariant also after multiple introductions. Phenotypic plasticity may therefore be an important component of the successful spread of Impatiens glandulifera across Europe.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0242-8) contains supplementary material, which is available to authorized users.

Interglacial refugia preserved high genetic diversity of the Chinese mole shrew in the mountains of southwest China

The mountains of southwest China (MSC) harbor extremely high species diversity; however, the mechanism behind this diversity is unknown. We investigated to what degree the topography and climate change shaped the genetic diversity and diversification in these mountains, and we also sought to identify the locations of microrefugia areas in these mountains. For these purposes, we sampled extensively to estimate the intraspecific phylogenetic pattern of the Chinese mole shrew (Anourosorex squamipes) in southwest China throughout its range of distribution. Two mitochondrial genes, namely, cytochrome b (CYT B) and NADH dehydrogenase subunit 2 (ND2), from 383 archived specimens from 43 localities were determined for phylogeographic and demographic analyses. We used the continuous-diffusion phylogeographic model, extensive Bayesian skyline plot species distribution modeling (SDM) and approximate Bayesian computation (ABC) to explore the changes in population size and distribution through time of the species. Two phylogenetic clades were identified, and significantly higher genetic diversity was preserved in the southern subregion of the mountains. The results of the SDM, continuous-diffusion phylogeographic model, extensive Bayesian skyline plot and ABC analyses were congruent and supported that the Last Interglacial Maximum (LIG) was an unfavorable period for the mole shrews because of a high degree of seasonality; A. squamipes survived in isolated interglacial refugia mainly located in the southern subregion during the LIG and rapidly expanded during the last glacial period. These results furnished the first evidence for major Pleistocene interglacial refugia and a latitudinal effect in southwest China, and the results shedding light on the higher level of species richness in the southern subregion.

Ancestral origins and invasion pathways in a globally invasive bird correlate with climate and influences from bird trade

Invasive species present a major threat to global biodiversity. Understanding genetic patterns and evolutionary processes that reinforce successful establishment is paramount for elucidating mechanisms underlying biological invasions. Among birds, the ring-necked parakeet (Psittacula krameri) is one of the most successful invasive species, established in over 35 countries. However, little is known about the evolutionary genetic origins of this species and what population genetic signatures tell us about patterns of invasion. We reveal the ancestral origins of populations across the invasive range and explore the potential influence of climate and propagule pressure from the pet trade on observed genetic patterns. Ring-necked parakeet samples representing the ancestral native range (n = 96) were collected from museum specimens, and modern samples from the invasive range (n = 855) were gathered from across Europe, Mauritius and Seychelles, and sequenced for two mitochondrial DNA markers comprising 868 bp of cytochrome b and control region, and genotyped at 10 microsatellite loci. Invasive populations comprise birds that originate predominantly from Pakistan and northern areas of India. Haplotypes associated with more northerly distribution limits in the ancestral native range were more prevalent in invasive populations in Europe, and the predominance of Asian haplotypes in Europe is consistent with the higher number of Asian birds transported by the pet trade outside the native range. Successful establishment of invasive species is likely to be underpinned by a combination of environmental and anthropogenic influences.

The Adaptive Change of HLA-DRB1 Allele Frequencies Caused by Natural Selection in a Mongolian Population That Migrated to the South of China

Pathogen-driven balancing selection determines the richness of human leukocyte antigen (HLA) alleles. Changes in the pathogen spectrum may cause corresponding changes in HLA loci. Approximately 700 years ago, a Mongolian population moved from the north of China to the Yunnan region in the south of China. The pathogen spectrum in the south of China differs from that in the north. In this study, changes in the HLA genes in the Yunnan Mongolian population, as well as the underlying mechanism, were investigated. A sequence-based typing method (SBT) was used to genotype HLA-DRB1 in 470 individuals from two Mongolian populations and another five ethnic groups. Meanwhile, 10 autosomal short tandem repeats (STRs) were genotyped to assess the influence of genetic background on HLA-DRB1 frequencies. The frequencies of certain alleles changed significantly in the Mongolian population that migrated to Yunnan. For example, DRB1*12:02:01 increased from 6.1% to 35.4%. STR analysis excluded the possibility of a recent bottleneck and indicated that 50% of the genetic consistency between northern and southern Mongolians; Tajima's D value for HLA-DRB1 exon2 and dN/dS analysis showed that the HLA-DRB1 genes in both Mongolian populations were under balancing selection. However, the sites under natural selection changed. We proposed that the dramatically change of HLA frequencies in southern Mongolian was caused by a combination of inter-population gene flow and natural selection. Certain diseases specific to the south of China, such as malaria, may be the driving force behind the enhanced DRB1*12:02:01 frequency.

Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

Interpreting the genetic structure of a metapopulation as the outcome of gene flow over a variety of timescales is essential for the proper understanding of how changes in landscape affect biological connectivity. Here we contrast historical and contemporary connectivity in two metapopulations of the freshwater fish Galaxias platei in northern and southernmost Patagonia where paleolakes existed during the Holocene and Pleistocene, respectively. Contemporary gene flow was mostly high and asymmetrical in the northern system while extremely reduced in the southernmost system. Historical migration patterns were high and symmetric in the northern system and high and largely asymmetric in the southern system. Both systems showed a moderate structure with a clear pattern of isolation by distance (IBD). Effective population sizes were smaller in populations with low contemporary gene flow. An approximate Bayesian computation (ABC) approach suggests a late Holocene colonization of the lakes in the northern system and recent divergence of the populations from refugial populations from east and west of the Andes. For the southern system, the ABC approach reveals that some of the extant G. platei populations most likely derive from an ancestral population inhabiting a large Pleistocene paleolake while the rest derive from a higher-altitude lake. Our results suggest that neither historical nor contemporary processes individually fully explain the observed structure and geneflow patterns and both are necessary for a proper understanding of the factors that affect diversity and its distribution. Our study highlights the importance of a temporal perspective on connectivity to analyse the diversity of spatially complex metapopulations.

Extensive population decline in the Tasmanian devil predates European settlement and devil facial tumour disease

The Tasmanian devil (Sarcophilus harrisii) was widespread in Australia during the Late Pleistocene but is now endemic to the island of Tasmania. Low genetic diversity combined with the spread of devil facial tumour disease have raised concerns for the species' long-term survival. Here, we investigate the origin of low genetic diversity by inferring the species' demographic history using temporal sampling with summary statistics, full-likelihood and approximate Bayesian computation methods. Our results show extensive population declines across Tasmania correlating with environmental changes around the last glacial maximum and following unstable climate related to increased 'El Niño-Southern Oscillation' activity.

Evaluating an icon of population persistence: the Devil's Hole pupfish

The Devil's Hole pupfish Cyprinodon diabolis has iconic status among conservation biologists because it is one of the World's most vulnerable species. Furthermore, C. diabolis is the most widely cited example of a persistent, small, isolated vertebrate population; a chronic exception to the rule that small populations do not persist long in isolation. It is widely asserted that this species has persisted in small numbers (less than 400 adults) for 10 000-20 000 years, but this assertion has never been evaluated. Here, we analyse the time series of count data for this species, and we estimate time to coalescence from microsatellite data to evaluate this hypothesis. We conclude that mean time to extinction is approximately 360-2900 years (median 410-1800), with less than a 2.1% probability of persisting 10 000 years. Median times to coalescence varied from 217 to 2530 years, but all five approximations had wide credible intervals. Our analyses suggest that Devil's Hole pupfish colonized this pool well after the Pleistocene Lakes receded, probably within the last few hundred to few thousand years; this could have occurred through human intervention.

Three-dimensional post-glacial expansion and diversification of an exploited oceanic fish

Vertical divergence in marine organisms is being increasingly documented, yet much remains to be carried out to understand the role of depth in the context of phylogeographic reconstruction and the identification of management units. An ideal study system to address this issue is the beaked redfish, Sebastes mentella – one of four species of ‘redfish’ occurring in the North Atlantic – which is known for a widely distributed ‘shallow‐pelagic’ oceanic type inhabiting waters between 250 and 550 m, and a more localized ‘deep‐pelagic’ population dwelling between 550 and 800 m, in the oceanic habitat of the Irminger Sea. Here, we investigate the extent of population structure in relation to both depth and geographic spread of oceanic beaked redfish throughout most of its distribution range. By sequencing the mitochondrial control region of 261 redfish collected over a decadal interval, and combining 160 rhodopsin coding nuclear sequences and previously genotyped microsatellite data, we map the existence of two strongly divergent evolutionary lineages with significantly different distribution patterns and historical demography, and whose genetic variance is mostly explained by depth. Combined genetic data, analysed via independent approaches, are consistent with a Late Pleistocene lineage split, where segregation by depth probably resulted from the interplay of climatic and oceanographic processes with life history and behavioural traits. The ongoing process of diversification in North Atlantic S. mentella may serve as an ‘hourglass’ to understand speciation and adaptive radiation in Sebastes and in other marine taxa distributed across a depth gradient.

Consequences of a demographic bottleneck on genetic structure and variation in the Scandinavian brown bear

The Scandinavian brown bear went through a major decline in population size approximately 100 years ago, due to intense hunting. After being protected, the population subsequently recovered and today numbers in the thousands. The genetic diversity in the contemporary population has been investigated in considerable detail, and it has been shown that the population consists of several subpopulations that display relatively high levels of genetic variation. However, previous studies have been unable to resolve the degree to which the demographic bottleneck impacted the contemporary genetic structure and diversity. In this study, we used mitochondrial and microsatellite DNA markers from pre- and postbottleneck Scandinavian brown bear samples to investigate the effect of the bottleneck. Simulation and multivariate analysis suggested the same genetic structure for the historical and modern samples, which are clustered into three subpopulations in southern, central and northern Scandinavia. However, the southern subpopulation appears to have gone through a marked change in allele frequencies. When comparing the mitochondrial DNA diversity in the whole population, we found a major decline in haplotype numbers across the bottleneck. However, the loss of autosomal genetic diversity was less pronounced, although a significant decline in allelic richness was observed in the southern subpopulation. Approximate Bayesian computations provided clear support for a decline in effective population size during the bottleneck, in both the southern and northern subpopulations. These results have implications for the future management of the Scandinavian brown bear because they indicate a recent loss in genetic diversity and also that the current genetic structure may have been caused by historical ecological processes rather than recent anthropogenic persecution.