Spatial inference of admixture proportions and secondary contact zones

The genetic signature left by the range expansion of red foxes in Australia is detectable after more than 80 years of population stability

Reconstructing biological invasions from historical sources can provide insights into how they occur but are difficult to do when invasions are poorly documented. Genetic signatures left by invaders can also offer insights into invasion routes, points of origin and general biology but often present conclusions that are contradictory to expectations. Here, we test the ability of continental-wide microsatellite genotype data from 29 loci and 3122 samples to reconstruct the well-documented invasion of red foxes Vulpes vulpes from the United Kingdom into Australia over 150 years ago, an invasion that has led to the extinction of many native species. Our analysis reveals several key signals of invasion evident in Australian foxes. They display lower levels of diversity than foxes sampled from the UK, exhibit clines in diversity from the point of introduction (south-east Australia) to the edge of their range, and show strong evidence of allele surfing in westerly and north-easterly directions. These characteristics are consistent with a single point of origin followed by rapid expansion in westerly and north-easterly directions as suggested by historical records. We also find little genetic structure in foxes across Australia with only the vast Nullarbor Plains and Great Victoria Desert region presenting a detectable barrier to their dispersal. As such, no mainland region within the current range of foxes can be considered genetically isolated and therefore appropriate for localised eradication efforts. Overall, our analyses demonstrate the ability of comprehensive population genetic studies to reconstruct invasion histories even after more than 80 years since colonisation was stabilised.

Population Genetic Structure and Diversity of Cryptic Species of the Plant Genus Macrocarpaea (Gentianaceae) from the Tropical Andes

The Pleistocene climatic oscillations (PCO) that provoked several cycles of glacial-interglacial periods are thought to have profoundly affected species distribution, richness and diversity around the world. While the effect of the PCO on population dynamics at temperate latitudes is well known, considerable questions remain about its impact on the biodiversity of neotropical mountains. Here, we use amplified fragment length polymorphism molecular markers (AFLPs) to investigate the phylogeography and genetic structure of 13 plant species belonging to the gentian genus Macrocarpaea (Gentianaceae) in the tropical Andes. These woody herbs, shrubs or small trees show complex and potentially reticulated relationships, including cryptic species. We show that populations of M. xerantifulva in the dry system of the Rio Marañón in northern Peru have lower levels of genetic diversity compared to other sampled species. We suggest that this is due to a recent demographic bottleneck resulting from the contraction of the montane wet forests into refugia because of the expansion of the dry system into the valley during the glacial cycles of the PCO. This may imply that the ecosystems of different valleys of the Andes might have responded differently to the PCO.

Genetic Variability and Admixture Zones in the Italian Populations of Turkey Oak (Quercus cerris L.)

The Turkey oak (Quercus cerris L.) is widely distributed in Italy, where it is the ecologically dominant oak on sandy and acidic soil. In this work, we analysed 23 natural populations by means of eight SSR (microsatellite) markers, to obtain the first synthetic map of genetic variability for this species and to study its dispersion during the Holocene, due to the possibility that at least one refugium during the Last Glacial Maximum was in Italy. The analyses showed a good amount of genetic variability together with fair differentiation between populations, as indicated by F_ST = 0.059. A Bayesian analysis of the amount of admixture among populations revealed the presence of four putative gene pools of origin and a rough subdivision of the populations according to their geographic location, as confirmed by the spatial analysis. No evidence for the existence of putative refugial populations was found; however, this study paves the way for the planning of conservation strategies also with regard to the relationship between Turkey oak and other oak species in Italy.

Genetic divergence at species boundaries of the dolphinfish (Coryphaena hippurus) in the Tropical Eastern Pacific

Background

Marine species constitute commercially important resources, and knowledge about mechanisms that shape phylogeographic patterns and genetic structure provides valuable information for conservation. The dolphinfish, Coryphaena hippurus, is one of the most important species caught in the Tropical Eastern Pacific (TEP). However, the lack of consensus about the existence of genetically differentiated populations in the area has hindered the adoption of management strategies to ensure its viability.

Methods

We assessed genetic variation and phylogeographic structure using two mitochondrial genes and 14 nuclear DNA microsatellite loci. Population genetic tools were used to characterize the spatial distribution of genetic variation of C. hippurus in the TEP, evaluate the extent of connectivity between dolphinfish populations, infer potential barriers to gene flow, and test for signals of contemporary and historical demographic expansions.

Results

Mitochondrial DNA sequences showed genetic homogeneity across locations in the TEP, as well as a strong signal of population expansion dated to the late Pleistocene. In contrast, nuclear microsatellite markers resolved four genetically distinct groups with a remarked genetic differentiation between the most distant locations, at the northern and southern boundaries of the species' range. High mean genetic diversity was found at all localities (Hs = 0.66-0.81). Notwithstanding, positive F _IS and low effective population size (Ne = 77.9-496.4) were also recorded.

Conclusions

The distribution of genetic variation could be related to expansion-contraction cycles following seasonal temperature changes at transitional areas, promoting population subdivisions. However, we cannot rule out the effect of oceanographic dynamics to the observed patterns. Although this marine species remains highly abundant despite commercial exploitation, the low Ne values are of conservation concern and must be considered in fishery management plans.

Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann’s tortoise Testudo hermanni

Background

Hotspots of intraspecific genetic diversity represent invaluable resources for species to cope with environmental changes, and their identification is increasingly recognized as a major goal of conservation ecology research. However, even for iconic and endangered species, conservation strategies are often planned without thorough information on the geographic patterns of genetic variation. Here, we investigated the spatial patterns of genetic variation of the endangered Hermann’s tortoise Testudo hermanni in the Italian Peninsula by genotyping 174 individuals at 7 microsatellite loci, with the aim to contribute to planning effective conservation strategies.

Results

Ordination-based and Bayesian clustering analyses consistently identified three main genetic clusters, one spread in the central and northern part of the peninsula, and two restricted to southern Italy and Sicily, respectively. The highest levels of genetic diversity were found in populations of the southern cluster and, in particular, at the northern edges of its distribution (He > 0.6, Ar > 2.8 ), that correspond to areas of putative secondary contact and admixture between distinct lineages. Our results clearly identify a hotspot of genetic diversity for the Hermann’s tortoise in southern Italy.

Conclusion

We inferred the evolutionary history and the spatial patterns of genetic variation of the Hermann’s tortoise in the Italian Peninsula. We identified three main genetic clusters along the peninsula and a hotspot of intraspecific diversity in southern Italy. Our results underline the urgent need for conservation actions to warrant the long-term persistence of viable tortoise populations in this area. Furthrmore, these data add further evidence to the role of southern Italy as a biodiversity hotspot for temperate fauna, claiming for higher consideration of this area in large scale conservation programs.

Contrasting genetic trajectories of endangered and expanding red fox populations in the western U.S

As anthropogenic disturbances continue to drive habitat loss and range contractions, the maintenance of evolutionary processes will increasingly require targeting measures to the population level, even for common and widespread species. Doing so requires detailed knowledge of population genetic structure, both to identify populations of conservation need and value, as well as to evaluate suitability of potential donor populations. We conducted a range-wide analysis of the genetic structure of red foxes in the contiguous western U.S., including a federally endangered distinct population segment of the Sierra Nevada subspecies, with the objectives of contextualizing field observations of relative scarcity in the Pacific mountains and increasing abundance in the cold desert basins of the Intermountain West. Using 31 autosomal microsatellites, along with mitochondrial and Y-chromosome markers, we found that populations of the Pacific mountains were isolated from one another and genetically depauperate (e.g., estimated Ne range = 3–9). In contrast, red foxes in the Intermountain regions showed relatively high connectivity and genetic diversity. Although most Intermountain red foxes carried indigenous western matrilines (78%) and patrilines (85%), the presence of nonindigenous haplotypes at lower elevations indicated admixture with fur-farm foxes and possibly expanding midcontinent populations as well. Our findings suggest that some Pacific mountain populations could likely benefit from increased connectivity (i.e., genetic rescue) but that nonnative admixture makes expanding populations in the Intermountain basins a non-ideal source. However, our results also suggest contact between Pacific mountain and Intermountain basin populations is likely to increase regardless, warranting consideration of risks and benefits of proactive measures to mitigate against unwanted effects of Intermountain gene flow.

Genomics and ecological modelling clarify species integrity in a confusing group of butterflies

Recent advances in both genomics and ecological modelling present new, multidisciplinary opportunities for resolving species boundaries and understanding the mechanisms that maintain their integrity in regions of contact. Here, we use a combination of high-throughput DNA sequencing and ecological niche modelling to resolve species boundaries and niche divergence within the Speyeria atlantis-hesperis (Lepidoptera: Nymphalidae) complex, a confusing group of North American butterflies. This complex is notorious for its muddled species delimitations, morphological ambiguity, and extensive mito-nuclear discordance. Our admixture and multispecies coalescent-based analyses of single nucleotide polymorphisms identified substantial divergences between S. atlantis and S. hesperis in areas of contact, as well as between distinct northern and southern lineages within S. hesperis. Our results also provide evidence of past introgression relating to another species, S. zerene, which previous work has shown to be more distantly related to the S. atlantis-hesperis complex. We then used ecological models to predict habitat suitability for each of the three recovered genomic lineages in the S. atlantis-hesperis complex and assess their pairwise niche divergence. These analyses resolved that these three lineages are significantly diverged in their respective niches and are not separated by discontinuities in suitable habitat that might present barriers to gene flow. We therefore infer that ecologically-mediated selection resulting in disparate habitat associations is a principal mechanism reinforcing their genomic integrity. Overall, our results unambiguously support significant evolutionary and ecological divergence between the northern and southern lineages of S. hesperis, sufficient to recognize the southern evolutionary lineage as a distinct species, called S. nausicaa based on taxonomic priority.

Comparing the population structure of the specialist Butler’s Gartersnake (Thamnophis butleri) and the generalist Eastern Gartersnake (Thamnophis sirtalis sirtalis) in Ontario (Canada) and Michigan (USA)

Species differing in life-history attributes vary in their responses to features within a shared landscape. We evaluated genetic structure of sympatric gartersnake species in southwestern Ontario (Canada) and southeastern Michigan (USA), where habitat fragmentation is high due to agriculture and urbanization. We surveyed genetic structure of a habitat specialist, Butler’s Gartersnake (Thamnophis butleri (Cope, 1889)), and a habitat generalist, Eastern Gartersnake (Thamnophis sirtalis sirtalis (Linnaeus, 1758)), using DNA microsatellites. Bayesian clustering, discriminant analysis of principal components, and pairwise population comparisons revealed genetic differentiation among three major regional clusters of Butler’s Gartersnake with evidence of further division within one. Genetic clustering of Butler’s Gartersnake suggest that inhospitable habitat limits dispersal. Eastern Gartersnakes showed less structure, with assignment tests implying a single genetic cluster. We found positive significant Mantel’s r for both species in the smallest distance class (<15 km), but significant isolation by distance for Butler’s Gartersnake only. These findings together imply that connectivity for Eastern Gartersnakes is less impacted by habitat loss and fragmentation or that we were less able to detect their effects. Our study shows the value of multispecies comparisons in studies seeking to understand the underlying causes of genetic structure in natural populations.

Landscape and climatic features drive genetic differentiation processes in a South American coastal plant

BACKGROUND

Historical and ecological processes shape patterns of genetic diversity in plant species. Colonization to new environments and geographical landscape features determine, amongst other factors, genetic diversity within- and differentiation between-populations. We analyse the genetic diversity and population structure of Calibrachoa heterophylla to infer the influence of abiotic landscape features on the level of gene flow in this coastal species of the South Atlantic Coastal Plain.

RESULTS

The C. heterophylla populations located on early-deposited coastal plain regions show higher genetic diversity than those closer to the sea. The genetic differentiation follows a pattern of isolation-by-distance. Landscape features, such as water bodies and wind corridors, and geographical distances equally explain the observed genetic differentiation, whereas the precipitation seasonality exhibits a strong signal for isolation-by-environment in marginal populations. The estimated levels of gene flow suggest that marginal populations had restricted immigration rates enhancing differentiation.

CONCLUSIONS

Topographical features related to coastal plain deposition history influence population differentiation in C. heterophylla. Gene flow is mainly restricted to nearby populations and facilitated by wind fields, albeit without any apparent influence of large water bodies. Furthermore, differential rainfall regimes in marginal populations seem to promote genetic differentiation.

Impact of Limited Dispersion Capacity and Natural Barriers on the Population Structure of the Grasshopper Ommexecha virens (Orthoptera: Ommexechidae)

The grasshopper Ommexecha virens Serville has low dispersion capacity, and it is regarded as a specialist, only being found in sandy, dry environments with high incidence of sunlight. Considering these aspects, we evaluated the diversity and genetic structure of O. virens natural populations using ISSR (Inter Simple Sequence Repeat) markers. The data pointed to low expected heterozygosity for some populations (H_E = 0.06-0.09), probably a consequence of positive inbreeding, which is typical of species showing low or null dispersion indices. Moreover, significant genetic differentiation was observed (F_ST = 0.50 and G_ST = 0.51), as well as low number of migrants (N_m = 0.47), indicating that the populations are genetically differentiated. This is likely related to the limitation in dispersing and fragmentation of suitable environment localities colonized by O. virens. The populations of O. virens were structured in three genetic groups associated to different landscapes, revealing the presence of a secondary contact zone, possibly arisen from isolation followed by genetic divergence among populations and subsequent gene flow of divergent individuals of O. virens. At last, we found positive isolation by distance (IBD; r: 0.427; P: 0.025) which is an important factor, since it may be adding to the emergence of reproductive barriers among individuals of O. virens that have been experiencing isolation.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F_ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F_ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F_ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species.

Integrative assessment of intraspecific diversification in Loggerhead Shrike (Lanius ludovicianus) provides insight on the geographic pattern of phenotypic divergence and process of speciation

Integrated studies of the geographical, ecological, and historical factors that shape intraspecific phenotypic and genetic variation can help us to decipher the processes leading to geographic patterns of population divergence and speciation. We quantify and compare morphological and genetic variation in the Loggerhead Shrike (Lanius ludovicianus Linnaeus, 1766), a broadly distributed passerine in North America with both migratory and non-migratory populations that occupy a diversity of habitats and topographies. The geographic distributions and patterns of differentiation among subspecies suggest that migration has strongly impacted population divergence, including the habit of migrating itself, but also dispersal. Patterns of mitochondrial and nuclear genetic differentiation can be attributed to female-biased dispersal and to increased dispersal rates and distances in migratory populations. Weak phenotypic differentiation among migratory versus migratory and non-migratory populations suggest that migration may more strongly affect morphology than adaptation to local habitats. Our results generally support previous subspecific designations with two notable exceptions. We found little genetic differentiation between two subspecies (Lanius ludovicianus gambeli Ridgway, 1887 and Lanius ludovicianus mexicanus C.L. Brehm, 1854), but identify a new, distinct subspecies, which we refer to as Lanius ludovicianus centralis ssp. nov.

Population genetic considerations for using biobanks as international resources in the pandemic era and beyond

The past years have seen the rise of genomic biobanks and mega-scale meta-analysis of genomic data, which promises to reveal the genetic underpinnings of health and disease. However, the over-representation of Europeans in genomic studies not only limits the global understanding of disease risk but also inhibits viable research into the genomic differences between carriers and patients. Whilst the community has agreed that more diverse samples are required, it is not enough to blindly increase diversity; the diversity must be quantified, compared and annotated to lead to insight. Genetic annotations from separate biobanks need to be comparable and computable and to operate without access to raw data due to privacy concerns. Comparability is key both for regular research and to allow international comparison in response to pandemics. Here, we evaluate the appropriateness of the most common genomic tools used to depict population structure in a standardized and comparable manner. The end goal is to reduce the effects of confounding and learn from genuine variation in genetic effects on phenotypes across populations, which will improve the value of biobanks (locally and internationally), increase the accuracy of association analyses and inform developmental efforts.

Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow-up study required to verify their geographic or evolutionary limits. Here, we follow-up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species' range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome-wide divergence (F ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species' range. Within lineages, gene flow was high, with low differentiation (mean F ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non-contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.

Analysis of Genetic Diversity and Population Structure of Orobanche foetida Populations From Tunisia Using RADseq

Orobanche foetida Poiret is a holoparasitic plant that lacks chlorophyll and totally depending on its host for its growth. Orobanche foetida parasitizes host plant roots and extract nutrient and water via a haustorium. Although O. foetida distributes in the Mediterranean region as a wild plant parasite, it parasitizes faba bean causing serious damages which may reach 90% yield losses in Tunisia. Analysis of genetic diversity of the parasite is important to better understand its evolution and spread, remained largely unknown. In this work, we present the first study on genetic diversity and population structure using the robust technique Restriction-site-Associated DNA sequencing (RADseq) for Orobanche spp. We collected 244 samples of O. foetida from 18 faba bean fields in the north of Tunisia including 17 populations from the north-west and one population form the north-east. To overcome the difficulty of SNP discovery in O. foetida genome as a non-model and tetraploid plant, we utilized three different informatics pipelines, namely UNEAK, pyRAD and Stacks. This study showed that genetic differentiation occurred in the Tunisian O. foetida emphasizing the isolation by distance effect. However, no strong population clustering was detected in this work basing on the three data sets and clustering methods used. The present study shed the light on the current distribution and the genetic variation situation of the fetid broomrape in Tunisia, highlighting the importance of understanding the evolution of this parasite and its genetic background. This will aid in developing efficient strategies to control this parasite and its expansion in Tunisia and worldwide.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Introduction to Population Genomics Methods

High-throughput sequencing technologies have provided an unprecedented opportunity to study the different evolutionary forces that have shaped present-day patterns of genetic diversity, with important implications for many directions in plant biology research. To manage such massive quantities of sequencing data, biologists, however, need new additional skills in informatics and statistics. In this chapter, our objective is to introduce population genomics methods to beginners following a learning-by-doing strategy in order to help the reader to analyze the sequencing data by themselves. Conducted analyses cover several main areas of evolutionary biology, such as an initial description of the evolutionary history of a given species or the identification of genes targeted by natural or artificial selection. In addition to the practical advices, we performed re-analyses of two cases studies with different kind of data: a domesticated cereal (African rice) and a non-domesticated tree species (sessile oak). All the code needed to replicate this work is publicly available on github ( https://github.com/ThibaultLeroyFr/Intro2PopGenomics/ ).

A clue to invasion success: genetic diversity quickly rebounds after introduction bottlenecks

One of the fundamental questions in invasion biology is to understand the genetic mechanisms behind success or failure during the establishment of a species. However, major limitations to understanding are usually a lack of spatiotemporal population data and information on the populations’ colonisation history. In a large-scale, detailed study on the bush-cricket Metrioptera roeselii 70 groups of founders were introduced in areas outside the species’ distribution range. We examined how (1) the number of founders (2–32 individuals), (2) the time since establishment (7 or 15 years after introduction) and (3) possible gene flow affected establishment success and temporal genetic changes of the introduced populations. We found higher establishment success in introductions with larger propagule sizes but genetic diversity indices were only partly correlated to propagule size. As expected, introduced populations were more similar to their founder population the larger the propagule size was. However, even if apparent at first, most of the differentiation in the small propagule introductions disappeared over time. Surprisingly, genetic variability was regained to a level comparable to the large and outbreeding founder population only 15 generations after severe demographic bottlenecks. We suggest that the establishment of these populations could be a result of several mechanisms acting in synergy. Here, rapid increase in genetic diversity of few introductions could potentially be attributed to limited gene flow from adjacent populations, behavioural adaptations and/or even increased mutation rate. We present unique insights into genetic processes that point towards traits that are important for understanding species’ invasiveness.

Limited divergent adaptation despite a substantial environmental cline in wild pea

Isolation by environment (IBE) is a widespread phenomenon in nature. It is commonly expected that the degree of difference among environments is proportional to the level of divergence between populations in their respective environments. It is therefore assumed that a species' genetic diversity displays a pattern of IBE in the presence of a strong environmental cline if gene flow does not mitigate isolation. We tested this common assumption by analysing the genetic diversity and demographic history of Pisum fulvum, which inhabits contrasting habitats in the southern Levant and is expected to display only minor migration rates between populations, making it an ideal test case. Ecogeographical and subpopulation structure were analysed and compared. The correlation of genetic with environmental distances was calculated to test the effect of isolation by distance and IBE and detect the main drivers of these effects. Historical effective population size was estimated using stairway plot. Limited overlap of ecogeographical and genetic clustering was observed, and correlation between genetic and environmental distances was statistically significant but small. We detected a sharp decline of effective population size during the last glacial period. The low degree of IBE may be the result of genetic drift due to a past bottleneck. Our findings contradict the expectation that strong environmental clines cause IBE in the absence of extensive gene flow.

The differential genetic signatures related to climatic landscapes for jaguars and pumas on a continental scale

Modern and paleoclimate changes may have altered species dynamics by shifting species' niche suitability over space and time. We analyze whether the current genetic structure and isolation of the two large American felids, jaguar (Panthera onca) and puma (Puma concolor), are mediated by changes in climatic suitability and connection routes over modern and paleoclimatic landscapes. We estimate species distribution under 5 climatic landscapes (modern, Holocene, last maximum glaciations [LMG], average suitability, and climatic instability) and correlate them with individuals' genetic isolation through causal modeling on a resemblance matrix. Both species exhibit genetic isolation patterns correlated with LMG climatic suitability, suggesting that these areas may have worked as "allele refuges." However, the jaguar showed higher vulnerability to climate changes, responding to modern climatic suitability and connection routes, whereas the puma showed a continuous and gradual transition of genetic variation. Despite differential responsiveness to climate change, both species are subjected to the climatic effects on genetic configuration, which may make them susceptible to future climatic changes, since these are progressing faster and with higher intensity than changes in the paleoclimate. Thus, the effects of climatic changes should be considered in the design of conservation strategies to ensure evolutionary and demographic processes mediated by gene flow for both species.

Spatial genetic structure in the rock hyrax (Procavia capensis) across the Namaqualand and western Fynbos areas of South Africa — a mitochondrial and microsatellite perspective

The interplay between biotic and abiotic environments is increasingly recognized as a major determinant of spatial genetic patterns. Among spatial genetic studies, saxicolous or rock-dwelling species remain underrepresented in spite of their strict dependence on landscape structure. Here we investigated patterns and processes operating at different spatial (fine and regional scales) and time scales (using mitochondrial and microsatellite markers) in the rock hyrax (Procavia capensis (Pallas, 1766)). Our focus was on the western seaboard of South Africa and included two recognized biodiversity hotspots (Cape Floristic Region and Succulent Karoo). At fine spatial scale, significant genetic structure was present between four rocky outcrops in an isolated population, likely driven by the social system of this species. At a broader spatial scale, ecological dependence on rocky habitat and population-level processes, in conjunction with landscape structure, appeared to be the main drivers of genetic diversity and structure. Large areas devoid of suitable rocky habitat (e.g., the Knersvlakte, Sandveld, and Cape Flats, South Africa) represent barriers to gene flow in the species, although genetic clusters closely follow climatic, geological, and phytogeographic regions, possibly indicating ecological specialization or adaptation as contributing factors enforcing isolation. Taken together, our study highlights the need to consider both intrinsic and extrinsic factors when investigating spatial genetic structures within species.

SPLATCHE3: simulation of serial genetic data under spatially explicit evolutionary scenarios including long-distance dispersal

SUMMARY

SPLATCHE3 simulates genetic data under a variety of spatially explicit evolutionary scenarios, extending previous versions of the framework. The new capabilities include long-distance migration, spatially and temporally heterogeneous short-scale migrations, alternative hybridization models, simulation of serial samples of genetic data and a large variety of DNA mutation models. These implementations have been applied independently to various studies, but grouped together in the current version.

AVAILABILITY AND IMPLEMENTATION

SPLATCHE3 is written in C++ and is freely available for non-commercial use from the website http://www.splatche.com/splatche3. It includes console versions for Linux, MacOs and Windows and a user-friendly GUI for Windows, as well as detailed documentation and ready-to-use examples.

Genetic Structure and Phylogeography of Tuber magnatum Populations

The ectomycorrhizal fungus Tuber magnatum produces the white truffle appreciated worldwide for its unique aroma. With respect to other Tuber spp. of economic interest, T. magnatum presents a narrower geographical range. This species has, in fact, long been considered endemic to Italy. However, over the last few decades several reports have documented the presence of white truffles in different Mediterranean countries and in particular in various areas of south-east Europe. In this study, samples from several Pannonian and Balkan countries such as Hungary, Serbia, Romania, Bulgaria and Greece have been collected and genotyped with microsatellite markers and the data merged with those available for Italian populations. Our objectives were to test whether Italian and south-east European populations are differentiated and to evaluate the genetic diversity of T. magnatum all over its distributional range. We show the genetic structure of T. magnatum populations with the differentiation of four main groups: northern Italy, central-northern Italy, southern Italy and the Balkan/Pannonian region. The present study allowed us to refine the evolutionary history of T. magnatum and track the possible post-glacial expansion route of this species. The assessment of T. magnatum’s genetic structure is not only of scientific relevance, but it is also important for the conservation and market traceability of this prestigious fungus.

The Multispecies Coalescent Over-Splits Species in the Case of Geographically Widespread Taxa

Many recent species delimitation studies rely exclusively on limited analyses of genetic data analyzed under the multispecies coalescent (MSC) model, and results from these studies often are regarded as conclusive support for taxonomic changes. However, most MSC-based species delimitation methods have well-known and often unmet assumptions. Uncritical application of these genetic-based approaches (without due consideration of sampling design, the effects of a priori group designations, isolation by distance, cytoplasmic-nuclear mismatch, and population structure) can lead to over-splitting of species. Here, we argue that in many common biological scenarios, researchers must be particularly cautious regarding these limitations, especially in cases of well-studied, geographically variable, and parapatrically distributed species complexes. We consider these points with respect to a historically controversial species group, the American milksnakes (Lampropeltis triangulum complex), using genetic data from a recent analysis (Ruane et al. 2014). We show that over-reliance on the program Bayesian Phylogenetics and Phylogeography, without adequate consideration of its assumptions and of sampling limitations, resulted in over-splitting of species in this study. Several of the hypothesized species of milksnakes instead appear to represent arbitrary slices of continuous geographic clines. We conclude that the best available evidence supports three, rather than seven, species within this complex. More generally, we recommend that coalescent-based species delimitation studies incorporate thorough analyses of geographic variation and carefully examine putative contact zones among delimited species before making taxonomic changes.

Fragmented and isolated: limited gene flow coupled with weak isolation by environment in the paleoendemic giant sequoia (Sequoiadendron giganteum)

PREMISE

Patterns of genetic structure across a species' range reflect the long-term interplay between genetic drift, gene flow, and selection. Given the importance of gene flow in preventing the loss of diversity through genetic drift among spatially isolated populations, understanding the dynamics of gene flow and the factors that influence connectivity across a species' range is a major goal for conservation of genetic diversity. Here we present a detailed look at gene flow dynamics of Sequoiadendron giganteum, a paleoendemic tree species that will likely face numerous threats due to climate change.

METHODS

We used microsatellite markers to examine nineteen populations of S. giganteum for patterns of genetic structure and to estimate admixture and rates of gene flow between eight population pairs. Also, we used Generalized Dissimilarity Models to elucidate landscape factors that shape genetic differentiation among populations.

RESULTS

We found minimal gene flow between adjacent groves in the northern disjunct range. In most of the southern portion of the range, groves showed a signal of connectivity which degrades to isolation in the extreme south. Geographic distance was the most important predictor of genetic dissimilarity across the range, with environmental conditions related to precipitation and temperature explaining a small, but significant, portion of the genetic variance.

CONCLUSIONS

Due to their isolation and unique genetic composition, northern populations of S. giganteum should be considered a high conservation priority. In this region, we suggest germplasm conservation as well as restoration planting to enhance genetic diversity.

Population genetic structure and habitat connectivity for jaguar (Panthera onca) conservation in Central Belize

BACKGROUND

Connectivity among jaguar (Panthera onca) populations will ensure natural gene flow and the long-term survival of the species throughout its range. Jaguar conservation efforts have focused primarily on connecting suitable habitat in a broad-scale. Accelerated habitat reduction, human-wildlife conflict, limited funding, and the complexity of jaguar behaviour have proven challenging to maintain connectivity between populations effectively. Here, we used non-invasive genetic sampling and individual-based conservation genetic analyses to assess genetic diversity and levels of genetic connectivity between individuals in the Cockscomb Basin Wildlife Sanctuary and the Maya Forest Corridor. We used expert knowledge and scientific literature to develop models of landscape permeability based on circuit theory with fine-scale landscape features as ecosystem types, distance to human settlements and roads to predict the most probable jaguar movement across central Belize.

RESULTS

We used 12 highly polymorphic microsatellite loci to identify 50 individual jaguars. We detected high levels of genetic diversity across loci (HE = 0.61, HO = 0.55, and NA = 9.33). Using Bayesian clustering and multivariate models to assess gene flow and genetic structure, we identified one single group of jaguars (K = 1). We identified critical areas for jaguar movement that fall outside the boundaries of current protected areas in central Belize. We detected two main areas of high landscape permeability in a stretch of approximately 18 km between Sittee River Forest Reserve and Manatee Forest Reserve that may increase functional connectivity and facilitate jaguar dispersal from and to Cockscomb Basin Wildlife Sanctuary. Our analysis provides important insights on fine-scale genetic and landscape connectivity of jaguars in central Belize, an area of conservation concern.

CONCLUSIONS

The results of our study demonstrate high levels of relatively recent gene flow for jaguars between two study sites in central Belize. Our landscape analysis detected corridors of expected jaguar movement between the Cockscomb Basin Wildlife Sanctuary and the Maya Forest Corridor. We highlight the importance of maintaining already established corridors and consolidating new areas that further promote jaguar movement across suitable habitat beyond the boundaries of currently protected areas. Continued conservation efforts within identified corridors will further maintain and increase genetic connectivity in central Belize.

Landscape genetics reveals unique and shared effects of urbanization for two sympatric pool-breeding amphibians

Metapopulation-structured species can be negatively affected when landscape fragmentation impairs connectivity. We investigated the effects of urbanization on genetic diversity and gene flow for two sympatric amphibian species, spotted salamanders (Ambystoma maculatum) and wood frogs (Lithobates sylvaticus), across a large (>35,000 km2) landscape in Maine, USA, containing numerous natural and anthropogenic gradients. Isolation-by-distance (IBD) patterns differed between the species. Spotted salamanders showed a linear and relatively high variance relationship between genetic and geographic distances (r = .057, p < .001), whereas wood frogs exhibited a strongly nonlinear and lower variance relationship (r = 0.429, p < .001). Scale dependence analysis of IBD found gene flow has its most predictable influence (strongest IBD correlations) at distances up to 9 km for spotted salamanders and up to 6 km for wood frogs. Estimated effective migration surfaces revealed contrasting patterns of high and low genetic diversity and gene flow between the two species. Population isolation, quantified as the mean IBD residuals for each population, was associated with local urbanization and less genetic diversity in both species. The influence of geographic proximity and urbanization on population connectivity was further supported by distance-based redundancy analysis and multiple matrix regression with randomization. Resistance surface modeling found interpopulation connectivity to be influenced by developed land cover, light roads, interstates, and topography for both species, plus secondary roads and rivers for wood frogs. Our results highlight the influence of anthropogenic landscape features within the context of natural features and broad spatial genetic patterns, in turn supporting the premise that while urbanization significantly restricts interpopulation connectivity for wood frogs and spotted salamanders, specific landscape elements have unique effects on these two sympatric species.

Multiple introductions and population structure during the rapid expansion of the invasive Sahara mustard (Brassica tournefortii)

The specific mechanisms that result in the success of any species invasion case are difficult to document. Reproductive strategies are often cited as a primary driver of invasive success, with human activities further facilitating invasions by, for example, acting as seed vectors for dispersal via road, train, air, and marine traffic, and by producing efficient corridors for movement including canals, drainages, and roadways. Sahara mustard (Brassica tournefortii) is a facultative autogamous annual native to Eurasia that has rapidly invaded the southwestern United States within the past century, displacing natives, and altering water-limited landscapes in the southwest. We used a genotyping-by-sequencing approach to study the population structure and spatial geography of Sahara mustard from 744 individuals from 52 sites across the range of the species' invasion. We also used herbaria records to model range expansion since its initial introduction in the 1920s. We found that Sahara mustard occurs as three populations in the United States unstructured by geography, identified three introduction sites, and combined herbaria records with genomic analyses to map the spread of the species. Low genetic diversity and linkage disequilibrium are consistent with self-fertilization, which likely promoted rapid invasive spread. Overall, we found that Sahara mustard experienced atypical expansion patterns, with a relatively constant rate of expansion and without the lag phase that is typical of many invasive species.

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R ST and F ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad-scale patterns of genetic structure and closely aligned with the R ST network. The F ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.

Killer whale genomes reveal a complex history of recurrent admixture and vicariance

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

The evolutionary history of the Cape hare (Lepus capensis sensu lato): insights for systematics and biogeography

Inferring the phylogeography of species with large distributions helps deciphering major diversification patterns that may occur in parallel across taxa. Here, we infer the evolutionary history of the Cape hare, Lepus capensis sensu lato, a species distributed from southern Africa to Asia, by analyzing variation at 18 microsatellites and 9 DNA (1 mitochondrial and 8 nuclear) sequenced loci, from field and museum-collected samples. Using a combination of assignment and coalescent-based methods, we show that the Cape hare is composed of five evolutionary lineages, distributed in distinct biogeographic regions-north-western Africa, eastern Africa, southern Africa, the Near East and the Arabian Peninsula. A deep phylogenetic break possibly dating to the Early Pleistocene was inferred between the African and Asian L. capensis groups, and the latter appear more closely related to other Eurasian hare species than to African Cape hares. The inferred phylogeographic structure is shared by numerous taxa distributed across the studied range, suggesting that environmental changes, such as the progressive aridification of the Saharo-Arabian desert and the fluctuations of savannah habitats in Sub-Saharan Africa, had comparable impacts across species. Fine-scale analyses of the western Sahara-Sahel populations showed rich fragmentation patterns for mitochondrial DNA but not for microsatellites, compatible with the environmental heterogeneity of the region and female philopatry. The complex evolutionary history of L. capensis sensu lato, which possibly includes interspecific gene flow, is not reflected by taxonomy. Integrating evolutionary inference contributes to an improved characterization of biodiversity, which is fundamental to foster the conservation of relevant evolutionary units.

Genes, geology and germs: gut microbiota across a primate hybrid zone are explained by site soil properties, not host species

Gut microbiota in geographically isolated host populations are often distinct. These differences have been attributed to between-population differences in host behaviours, environments, genetics and geographical distance. However, which factors are most important remains unknown. Here, we fill this gap for baboons by leveraging information on 13 environmental variables from 14 baboon populations spanning a natural hybrid zone. Sampling across a hybrid zone allowed us to additionally test whether phylosymbiosis (codiversification between hosts and their microbiota) is detectable in admixed, closely related primates. We found little evidence of genetic effects: none of host genetic ancestry, host genetic relatedness nor genetic distance between host populations were strong predictors of baboon gut microbiota. Instead, gut microbiota were best explained by the baboons' environments, especially the soil's geologic history and exchangeable sodium. Indeed, soil effects were 15 times stronger than those of host-population FST, perhaps because soil predicts which foods are present, or because baboons are terrestrial and consume soil microbes incidentally with their food. Our results support an emerging picture in which environmental variation is the dominant predictor of host-associated microbiomes. We are the first to show that such effects overshadow host species identity among members of the same primate genus.

Pleistocene-dated biogeographic barriers drove divergence within the Australo-Papuan region in a sex-specific manner: an example in a widespread Australian songbird

Understanding how environmental change has shaped species evolution can inform predictions of how future climate change might continue to do so. Research of widespread biological systems spanning multiple climates that have been subject to environmental change can yield generalizable inferences about the neutral and adaptive processes driving lineage divergence during periods of environmental change. We contribute to the growing body of multi-locus phylogeographic studies investigating the effect of Pleistocene climate change on species evolution by focusing on a widespread Australo-Papuan songbird with several mitochondrial lineages that diverged during the Pleistocene, the grey shrike-thrush (Colluricincla harmonica). We employed multi-locus phylogenetic, population genetic and coalescent analyses to (1) assess whether nuclear genetic diversity suggests a history congruent with that based on phenotypically defined subspecies ranges, mitochondrial clade boundaries and putative biogeographical barriers, (2) estimate genetic diversity within and genetic differentiation and gene flow among regional populations and (3) estimate population divergence times. The five currently recognized subspecies of grey shrike-thrush are genetically differentiated in nuclear and mitochondrial genomes, but connected by low levels of gene flow. Divergences among these populations are concordant with recognized historical biogeographical barriers and date to the Pleistocene. Discordance in the order of population divergence events based on mitochondrial and nuclear genomes suggests a history of sex-biased gene flow and/or mitochondrial introgression at secondary contacts. This study demonstrates that climate change can impact sexes with different dispersal biology in different ways. Incongruence between population and mitochondrial trees calls for a genome-wide investigation into dispersal, mitochondrial introgression and mitonuclear evolution.

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.

Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species' range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals.

Is population genetic structure of vascular plants shaped more by ecological or geographic factors? A study case on the Mediterranean endemic Centaurea filiformis (Asteraceae)

All known populations of the Sardinian endemic Centaurea filiformis Viv. (Asteraceae) were studied in order to understand the impact of both geographic and ecological factors on the genetic structuring of this species. Fourteen populations and 234 individuals were sampled. The demographic structure of the populations and the reproductive ecology were estimated in 28 plots. Population genetic analyses were based on SSR markers. Genetic structure was investigated by spatial Bayesian methods. Average densities of 0.51 individuals m^-2 were detected, with a prevalence of adults. Ten species of pollinators were identified; C. filiformis ability to self-pollinate and myrmecochory were demonstrated experimentally. The populations displayed an average heterozygosity value of H_e  = 0.576 and high genetic differentiation (overall F_ST  = 0.218). Bayesian analysis suggests that five is the most probable number of gene pools of origin. A strong correlation between geographic distances and genetic distances among populations was highlighted. The demographic population structure of C. filiformis is dominated by adults, suggesting that it is a stable-regressive or senile species, investing more in local persistence than colonisation ability. Despite the scattered distribution, the populations studied do not present evidence of genetic erosion. The analysis of genetic differentiation reveals very high differentiation levels among populations, thus indicating that effective barriers exist against gene flow. A general conclusion is that population distribution results in a clear genetic structure for the populations studied, and that geography and not ecology is shaping the present distribution of this species.