Revisiting the phylogeography and demography of European badgers (Meles meles) based on broad sampling, multiple markers and simulations

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Genetic evidence further elucidates the history and extent of badger introductions from Great Britain into Ireland

The colonization of Ireland by mammals has been the subject of extensive study using genetic methods and forms a central problem in understanding the phylogeography of European mammals after the Last Glacial Maximum. Ireland exhibits a depauperate mammal fauna relative to Great Britain and continental Europe, and a range of natural and anthropogenic processes have given rise to its modern fauna. Previous Europe-wide surveys of the European badger (Meles meles) have found conflicting microsatellite and mitochondrial DNA evidence in Irish populations, suggesting Irish badgers have arisen from admixture between human imported British and Scandinavian animals. The extent and history of contact between British and Irish badger populations remains unclear. We use comprehensive genetic data from Great Britain and Ireland to demonstrate that badgers in Ireland's northeastern and southeastern counties are genetically similar to contemporary British populations. Simulation analyses suggest this admixed population arose in Ireland 600-700 (CI 100-2600) years before present most likely through introduction of British badgers by people. These findings add to our knowledge of the complex colonization history of Ireland by mammals and the central role of humans in facilitating it.

Holocene changes in the distributions of Asian and European badgers (Carnivora: Mustelidae: Meles) inferred from ancient DNA analysis

Although the present-day distributional boundary between the European badger (Meles meles) and the Asian badger (Meles leucurus) is around the Volga River, studies of ancient bone remains have indicated changes in the distribution of M. meles and M. leucurus in the Urals–Volga region during the Holocene. To examine past changes in distribution using genetic data, changes in genetic diversity, and the relationships of Holocene to modern populations, we sequenced ~150 bp of the mitochondrial DNA control region from the 44 ancient badger remains excavated from European Russian, Ural and Western Siberian sites, and we detected 12 haplotypes. Our study revealed Holocene changes in the distributional boundary between these badger species. Meles meles inhabited the Ural Mountains east of the Volga River in the Early Holocene, whereas M. leucurus expanded its distribution westwards, starting ~2500 years ago. Thereafter, M. leucurus rapidly replaced M. meles in the region between the Urals and the Volga, resulting in the present-day boundary in the Volga–Kama region. Among the 12 haplotypes detected, three for M. leucurus and four for M. meles were identical to partial sequences of haplotypes detected in modern populations, indicating considerable genetic continuity between Holocene and modern populations.

Genetic analyses reveal population structure and recent decline in leopards (Panthera pardus fusca) across the Indian subcontinent

Background

Large carnivores maintain the stability and functioning of ecosystems. Currently, many carnivore species face declining population sizes due to natural and anthropogenic pressures. The leopard, Panthera pardus, is probably the most widely distributed and highly adaptable large felid globally, still persisting in most of its historic range. However, we lack subspecies-level data on country or regional scale on population trends, as ecological monitoring approaches are difficult to apply on such wide-ranging species. We used genetic data from leopards sampled across the Indian subcontinent to investigate population structure and patterns of demographic decline.

Methods

We collected faecal samples from the Terai-Arc landscape of northern India and identified 56 unique individuals using a panel of 13 microsatellite markers. We merged this data with already available 143 leopard individuals and assessed genetic structure at country scale. Subsequently, we investigated the demographic history of each identified subpopulations and compared genetic decline analyses with countrywide local extinction probabilities.

Results

Our genetic analyses revealed four distinct subpopulations corresponding to Western Ghats, Deccan Plateau-Semi Arid, Shivalik and Terai region of the north Indian landscape, each with high genetic variation. Coalescent simulations with microsatellite loci revealed a possibly human-induced 75–90% population decline between ∼120–200 years ago across India. Population-specific estimates of genetic decline are in concordance with ecological estimates of local extinction probabilities in these subpopulations obtained from occupancy modeling of the historic and current distribution of leopards in India.

Conclusions

Our results confirm the population decline of a widely distributed, adaptable large carnivore. We re-iterate the relevance of indirect genetic methods for such species in conjunction with occupancy assessment and recommend that detailed, landscape-level ecological studies on leopard populations are critical to future conservation efforts. Our approaches and inference are relevant to other widely distributed, seemingly unaffected carnivores such as the leopard.

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.

Investigation into the genetic diversity in toll-like receptors 2 and 4 in the European badger Meles meles

The Toll-like receptor (TLR) genes are a conserved family of genes central to the innate immune response to pathogen infection. They encode receptor proteins, recognise pathogen associated molecular patterns (PAMPs) and trigger initial immune responses. In some host-pathogen systems, it is reported that genetic differences, such as single nucleotide polymorphisms (SNPs), associate with disease resistance or susceptibility. Little is known about TLR gene diversity in the European badger (Meles meles). We collected DNA from UK badgers, carried out PCR amplification of the badger TLR2 gene and exon 3 of TLR4 and determined DNA sequences for individual badgers for TLR2 (n = 61) and TLR4 exon 3 (n = 59). No polymorphism was observed in TLR4. Three TLR2 amino acid haplotype variants were found. Ninety five percent of badgers were homozygous for one common haplotype (H1), the remaining three badgers had genotypes H1/H3, H1/H2 and H2/H2. By broad comparison with other species, diversity in TLR genes in badgers seems low. This could be due to a relatively localised sampling or inherent low genetic diversity. Further studies are required to assess the generality of the low observed diversity and the relevance to the immunological status of badgers.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

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

Habitat associations are a function of habitat preferences and dispersal capabilities, both of which can influence how species responded to Quaternary climatic changes and contemporary habitat heterogeneity. Predicting resultant genetic structure is not always straightforward, especially in species where high dispersal potential and habitat preferences yield opposing predictions. The American badger has high dispersal capabilities that predict widespread panmixia, but avoids closed-canopy forests and clay soils, which could restrict gene flow and create ecologically based population genetic structure. We used mitochondrial sequence and microsatellite data sets to characterize how these opposing forces contribute to genetic structure in badgers at a continent-wide scale. Our data revealed an overall lack of ecologically based population genetic structure, suggesting that high dispersal capabilities were sufficiently realized to overcome most habitat-based genetic structure. At a broadscale, badger gene flow is limited only by geographic distance (isolation by distance) and large water barriers (Lake Michigan and the Mississippi River). The absence of genetic structure in a species with strong avoidance of unsuitable habitats advances our understanding of when and how genetic structure emerges in widespread, highly mobile species.

Genetic diversity and phylogeography of Siberian roe deer, Caproulus pygargus, in central and peripheral populations

Current understanding of phylogeographical structure and genetic diversity of Siberian roe deer remains limited mainly due to small sample size and/or low geographical coverage in previous studies. Published data suggest at least two phylogroups: western (Ural Mountains and Western Siberia) and eastern (east from lake Baikal, including the Korean peninsula), but their phylogenetic relationship remains unclear. Combined sequences of cytochrome b (1140 bp) and the mtDNA control region (963 bp) were analyzed from 219 Siberian roe deer from 12 locations in Russia, Mongolia, and South Korea, which cover a large part of its range, to assess genetic diversity and phylogeographical status. Special emphasis was placed on the demographic history and genetic features of central, peripheral, and isolated populations. Results of median-joining network and phylogenetic tree analyses indicate that Siberian roe deer from the Urals to the Pacific Ocean are genetically diverse and that geographical distribution and composition of haplogroups coincide with previously described ranges of the subspecies Capreolus pygargus pygargus and Capreolus pygargus tianschanicus. We found that peripheral populations in the northwestern parts of the species range (Urals), as well as the isolated population from Jeju Island, are genetically distinct from those in the core part of the range, both in terms of genetic diversity and quantitative composition of haplogroups. We also found that northwestern (Urals) and northern (Yakutia) peripheral populations share the same haplogroup and fall into the same phylogenetic clade with the isolated population from Jeju Island. This finding sheds light on the taxonomic status of the Jeju Island population and leads to hypotheses about the discordance of morphological and genetic evolution in isolated populations and specific genetic features of peripheral populations.

Genetic Structure and Effective Population Sizes in European Red Deer (Cervus elaphus) at a Continental Scale: Insights from Microsatellite DNA

We analyzed more than 600 red deer (Cervus elaphus) from large parts of its European distribution range at 13 microsatellite loci, presenting the first continent-wide study of this species using nuclear markers. Populations were clearly differentiated (overall F ST = 0.166, Jost's D est = 0.385), and the BAPS clustering algorithm yielded mainly geographically limited and adjacent genetic units. When forced into only 3 genetic clusters our data set produced a very similar geographic pattern as previously found in mtDNA phylogeographic studies: a western group from Iberia to central and parts of Eastern Europe, an eastern group from the Balkans to Eastern Europe, and a third group including the threatened relict populations from Sardinia and Mesola in Italy. This result was also confirmed by a multivariate approach to analyzing our data set, a discriminant analysis of principal components. Calculations of genetic diversity and effective population sizes (linkage disequilibrium approach) yielded the lowest results for Italian (Sardinia, Mesola; N e between 2 and 8) and Scandinavian red deer, in line with known bottlenecks in these populations. Our study is the first to present comparative nuclear genetic data in red deer across Europe and may serve as a baseline for future analyses of genetic diversity and structuring in this widespread ungulate.