High gene flow in the American badger overrides habitat preferences and limits broadscale genetic structure

Disease-driven top predator decline affects mesopredator population genomic structure

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

Absence of genetic isolation across highly fragmented landscape in the ant Temnothorax nigriceps

Background

Human activities, including changes in agricultural landscapes, often impact biodiversity through habitat fragmentation. This potentially reduces genetic exchange between previously connected populations. Using a combination of nuclear and mitochondrial markers, we investigated (i) genetic diversity and population structure at multiple spatial scales and (ii) colony genetic structure and queen mating frequency in the ant species Temnothorax nigriceps in a highly anthropized environment.

Results

Although the results highlighted genetic structure on a European spatial scale, they did not reveal an impact of fragmentation on a regional scale, and we did not observe any genetic population structure on a regional scale. Across all populations, regardless of their geographical location, colony structure suggested monogyny (a single queen per colony) and monandry (single mating). However, nestmates were more related than expected, indicating that large-scale dispersal does not fully prevent genetic isolation.

Conclusions

Despite living in fragmented patches of habitat, populations of Temnothorax nigriceps are apparently genetically not isolated at a regional scale. However, large-scale dispersal alone does not prevent genetic isolation. The ecological requirements of T. nigriceps may explain their resilience to habitat fragmentation by allowing them to survive in very small patches of suitable habitat. The deeper investigation of the diversity of functional habitats for this species should allow to appreciate better the mechanisms permitting this species to overcome the negative impacts of fragmentation.

First Genome Sequence of the Gunnison's Prairie Dog (Cynomys gunnisoni), a Keystone Species and Player in the Transmission of Sylvatic Plague

Prairie dogs (genus Cynomys) are a charismatic symbol of the American West. Their large social aggregations and complex vocalizations have been the subject of scientific and popular interest for decades. A large body of literature has documented their role as keystone species of western North America's grasslands: They generate habitat for other vertebrates, increase nutrient availability for plants, and act as a food source for mammalian, squamate, and avian predators. An additional keystone role lies in their extreme susceptibility to sylvatic plague (caused by Yersinia pestis), which results in periodic population extinctions, thereby generating spatiotemporal heterogeneity in both biotic communities and ecological processes. Here, we report the first Cynomys genome for a Gunnison's prairie dog (C. gunnisoni gunnisoni) from Telluride, Colorado (USA). The genome was constructed using a hybrid assembly of PacBio and Illumina reads and assembled with MaSuRCA and PBJelly, which resulted in a scaffold N50 of 824 kb. Total genome size was 2.67 Gb, with 32.46% of the bases occurring in repeat regions. We recovered 94.9% (91% complete) of the single copy orthologs using the mammalian Benchmarking Universal Single-Copy Orthologs database and detected 49,377 gene models (332,141 coding regions). Pairwise Sequentially Markovian Coalescent showed support for long-term stable population size followed by a steady decline beginning near the end of the Pleistocene, as well as a recent population reduction. The genome will aid in studies of mammalian evolution, disease resistance, and the genomic basis of life history traits in ground squirrels.

Multiscale patterns of isolation by ecology and fine-scale population structure in Texas bobcats

Patterns of spatial genetic variation can be generated by a variety of ecological processes, including individual preferences based on habitat. These ecological processes act at multiple spatial and temporal scales, generating scale-dependent effects on gene flow. In this study, we focused on bobcats (Lynx rufus), a highly mobile, generalist felid that exhibits ecological and behavioral plasticity, high abundance, and broad connectivity across much of their range. However, bobcats also show genetic differentiation along habitat breaks, a pattern typically observed in cases of isolation-by-ecology (IBE). The IBE observed in bobcats is hypothesized to occur due to habitat-biased dispersal, but it is unknown if this occurs at other habitat breaks across their range or at what spatial scale IBE becomes most apparent. Thus, we used a multiscale approach to examine isolation by ecology (IBE) patterns in bobcats (Lynx rufus) at both fine and broad spatial scales in western Texas. We genotyped 102 individuals at nine microsatellite loci and used partial redundancy analysis (pRDA) to test if a suite of landscape variables influenced genetic variation in bobcats. Bobcats exhibited a latitudinal cline in population structure with a spatial signature of male-biased dispersal, and no clear barriers to gene flow. Our pRDA tests revealed high genetic similarity in similar habitats, and results differed by spatial scale. At the fine spatial scale, herbaceous rangeland was an important influence on gene flow whereas mixed rangeland and agriculture were significant at the broad spatial scale. Taken together, our results suggests that complex interactions between spatial-use behavior and landscape heterogeneity can create non-random gene flow in highly mobile species like bobcats. Furthermore, our results add to the growing body of data highlighting the importance of multiscale study designs when assessing spatial genetic structure.

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (∼0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

Canada lynx (Lynx canadensis) gene flow across a mountain transition zone in western North America

Mountain ecotones have the potential to cause multiple patterns in divergence, from simple barrier effects to more fundamental ecological divergence. Most work in mountain ecotones in North America has focused on reinforcement between refugial populations, making prediction of how mountains impact species that are not restricted to separate glacial refugia remains difficult. This study focused on the Canada lynx (Lynx canadensis Kerr, 1792), a highly mobile felid considered to be a habitat and dietary specialist. Specifically, we used 14 microsatellite loci and landscape genetic tools to investigate if the Rocky Mountains and associated climatic transitions influence lynx genetic differentiation in western North America. Although lynx exhibited high gene flow across the region, analyses detected structuring of neutral genetic variation across our study area. Gene flow for lynx most strongly related to temperature and elevation compared with other landscape variables (terrain roughness, percent forest cover, and habitat suitability index) and geographic distance alone. Overall, genetic structure in lynx is most consistent with barrier effects created by the Rocky Mountains rather than ecological divergence. Furthermore, warmer temperatures had a measurable impact on gene flow, which suggests connectivity may further decrease in peripheral or fragmented populations as climate warms.

Population genetic structures at multiple spatial scales: importance of social groups in European badgers

Population viability and metapopulation dynamics are strongly affected by gene flow. Identifying ecological correlates of genetic structure and gene flow in wild populations is therefore a major issue both in evolutionary ecology and species management. Studying the genetic structure of populations also enables identification of the spatial scale at which most gene flow occurs, hence the scale of the functional connectivity, which is of paramount importance for species ecology. In this study, we examined the genetic structure of a social, continuously distributed mammal, the European badger (Meles meles), both at large spatial scales (among populations) and fine (within populations) spatial scales. The study was carried out in 11 sites across France utilizing a noninvasive hair trapping protocol at 206 monitored setts. We identified 264 badgers genotyped at 24 microsatellite DNA loci. At the large scale, we observed high and significant genetic differentiation among populations (global Fst = 0.139; range of pairwise Fst [0.046–0.231]) that was not related to the geographic distance among sites, suggesting few large-scale dispersal events. Within populations, we detected a threshold value below which badgers were genetically close (< 400 m), highlighting that sociality is the major structuring process within badger populations at the fine scale.

Spatiotemporal analyses suggest the role of glacial history and the ice-free corridor in shaping American badger population genetic variation

Recurring glacial cycles through the Quaternary period drastically altered the size and distribution of natural populations of North American flora and fauna. The "southerly refugia model" has been the longstanding framework for testing the effects of glaciation on contemporary genetic patterns; however, insights from ancient DNA have contributed to the reconstruction of more complex histories for some species. The American badger, Taxidea taxus, provides an interesting species for exploring the genetic legacy of glacial history, having been hypothesized to have postglacially emerged from a single, southerly refugium to recolonize northern latitudes. However, previous studies have lacked genetic sampling from areas where distinct glacial refugia have been hypothesized, including the Pacific Northwest and American Far North (Yukon, Alaska). In order to further investigate the phylogeographic history of American badgers, we collected mitochondrial DNA sequence data from ancient subfossil material collected within the historical range (Alaska, Yukon) and combined them with new and previously published data from across the species' contemporary distribution (n = 1,207). We reconstructed a mostly unresolved phylogenetic tree and star-like haplotype network indicative of emergence from a largely panmictic glacial refugium and recent population expansion, the latter further punctuated by significantly negative Tajima's D and Fu's Fs values. Although directionality of migration cannot be unequivocally inferred, the moderate to high levels of genetic variation exhibited by American badgers, alongside the low frequency of haplotypes with indels in the Midwest, suggest a potential recolonization into central North America after the hypothesized ice-free corridor reopened ~13,000 years ago. Overall, the expanded reconstruction of phylogeographic history of American badgers offers a broader understanding of contemporary range-wide patterns and identifies unique genetic units that can likely be used to inform conservation of at-risk populations at the northern periphery.

Partitioning drivers of spatial genetic variation for a continuously distributed population of boreal caribou: Implications for management unit delineation

Isolation by distance (IBD) is a natural pattern not readily incorporated into theoretical models nor traditional metrics for differentiating populations, although clinal genetic differentiation can be characteristic of many wildlife species. Landscape features can also drive population structure additive to baseline IBD resulting in differentiation through isolation-by-resistance (IBR). We assessed the population genetic structure of boreal caribou across western Canada using nonspatial (STRUCTURE) and spatial (MEMGENE) clustering methods and investigated the relative contribution of IBD and IBR on genetic variation of 1,221 boreal caribou multilocus genotypes across western Canada. We further introduced a novel approach to compare the partitioning of individuals into management units (MU) and assessed levels of genetic connectivity under different MU scenarios. STRUCTURE delineated five genetic clusters while MEMGENE identified finer-scale differentiation across the study area. IBD was significant and did not differ for males and females both across and among detected genetic clusters. MEMGENE landscape analysis further quantified the proportion of genetic variation contributed by IBD and IBR patterns, allowing for the relative importance of spatial drivers, including roads, water bodies, and wildfires, to be assessed and incorporated into the characterization of population structure for the delineation of MUs. Local population units, as currently delineated in the boreal caribou recovery strategy, do not capture the genetic variation and connectivity of the ecotype across the study area. Here, we provide the tools to assess fine-scale spatial patterns of genetic variation, partition drivers of genetic variation, and evaluate the best management options for maintaining genetic connectivity. Our approach is highly relevant to vagile wildlife species that are of management and conservation concern and demonstrate varying degrees of IBD and IBR with clinal spatial genetic structure that challenges the delineation of discrete population boundaries.

The population and landscape genetics of the European badger (Meles meles) in Ireland

The population genetic structure of free-ranging species is expected to reflect landscape-level effects. Quantifying the role of these factors and their relative contribution often has important implications for wildlife management. The population genetics of the European badger (Meles meles) have received considerable attention, not least because the species acts as a potential wildlife reservoir for bovine tuberculosis (bTB) in Britain and Ireland. Herein, we detail the most comprehensive population and landscape genetic study of the badger in Ireland to date-comprised of 454 Irish badger samples, genotyped at 14 microsatellite loci. Bayesian and multivariate clustering methods demonstrated continuous clinal variation across the island, with potentially distinct differentiation observed in Northern Ireland. Landscape genetic analyses identified geographic distance and elevation as the primary drivers of genetic differentiation, in keeping with badgers exhibiting high levels of philopatry. Other factors hypothesized to affect gene flow, including earth worm habitat suitability, land cover type, and the River Shannon, had little to no detectable effect. By providing a more accurate picture of badger population structure and the factors effecting it, these data can guide current efforts to manage the species in Ireland and to better understand its role in bTB.