Phylogeographic history of grey wolves in Europe

Assessment of genetic diversity, population structure and wolf-dog hybridisation in the Eastern Romanian Carpathian wolf population

The Carpathian Mountains have been constantly inhabited by grey wolves and present one of the largest distribution areas in Europe, comprising between 2300 and 2700 individuals in Romania. To date, however, relatively little is known about the Romanian wolf population. We aimed to provide a first assessment of genetic diversity, population structure and wolf-dog hybridisation based on 444 mostly non-invasively collected samples in the Eastern Romanian Carpathians. Pack reconstruction and analysis of population genetic parameters were performed with mitochondrial DNA control-region sequencing and microsatellite genotyping. We found relatively high levels of genetic diversity, which is similar to values found in previous studies on Carpathian wolves from Poland and Slovakia, as well as to the long-lasting Dinaric-Balkan wolf population. We found no significant population structure in our study region, suggesting effective dispersal and admixture. Analysis of wolf-dog hybridisation using a Single Nucleotide Polymorphism panel optimised for hybrid detection revealed low rates of admixture between wolves and domestic dogs. Our results provide evidence for the existence of a genetically viable wolf population in the Romanian Carpathians. The genetic data obtained in this study may serve as valuable baseline information for the elaboration of monitoring standards and management plans for wolves in Romania.

Being a Dog: A Review of the Domestication Process

The process of canine domestication represents certainly one of the most interesting questions that evolutionary biology aims to address. A "multiphase" view of this process is now accepted, with a first phase during which different groups of wolves were attracted by the anthropogenic niche and a second phase characterized by the gradual establishment of mutual relationships between wolves and humans. Here, we provide a review of dog (Canis familiaris) domestication, highlighting the ecological differences between dogs and wolves, analyzing the molecular mechanisms which seem to have influenced the affiliative behaviors first observed in Belyaev's foxes, and describing the genetics of ancient European dogs. Then, we focus on three Mediterranean peninsulas (Balkan, Iberian and Italian), which together represent the main geographic area for studying canine domestication dynamics, as it has shaped the current genetic variability of dog populations, and where a well-defined European genetic structure was pinpointed through the analysis of uniparental genetic markers and their phylogeny.

Evaluating the functional, sexual and seasonal variation in the chemical constituents from feces of adult Iberian wolves (Canis lupus signatus)

Chemical signals deposited in feces play an important role in intraspecific and interspecific communication of many mammals. We collected fresh feces of adult wolves from wild breeding groups. All samples visually identified as belonging to wolves were subsequently identified to species level by sequencing a small fragment of mtDNA and sexed typing DBX6 and DBY7 sex markers. Using gas chromatography-mass spectrometry (GC-MS), we identified 56 lipophilic compounds in the feces, mainly heterocyclic aromatic organic compounds, such as indole or phenol, but also steroids, such as cholesterol, carboxylic acids and their esters between n-C₄ and n-C₁₈, aldehydes, alcohols and significant quantities of squalene and α-tocopherol, which would increase the chemical stability of feces on humid substrates. There was variability in the number and proportions of compounds between sexes, which could be indicative of their function as chemical signals. We also found variability in different reproductive states, especially in odorous compounds, steroids and α-tocopherol. Feces with a presumed marking function had higher proportions of α-tocopherol and steroids than feces with non-marking function. These compounds could be involved in intragroup and intergroup communication of wolves and their levels in feces could be directly related with the wolf's sex and physiological and reproductive status.

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

Genomic analyses reveal range-wide devastation of sea otter populations

The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species' range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 Mb of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76×) and three historical Californian samples from ~1500 and ~200 years ago to low coverage (1.5-3.5×). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

The origin and population genetics of wolves in the north Hungarian mountains

The grey wolf (Canis lupus) is one of the most challenging species to conserve in our modern and crowded world. Due to various factors, most European wolf populations are currently growing. In Hungary, numbers have increased since the 2000s. Although spontaneous recolonisation from Slovakia is considered to be the most likely mechanism by the majority of experts, some stakeholders claim that hand-reared individuals have been released. To determine the origin of wolves in northern Hungary, we analysed samples of free-ranging wolves collected in Slovakia and Hungary as well as samples from wolves in private enclosures in the region. We also included reference samples from domestic dogs. All samples were genotyped at 14 canine autosomal tetranucleotide microsatellite loci (STR) and analysed using multivariate, Bayesian methods. Hungarian wolf samples were also analysed using kinship methods. In the free-ranging wolf samples, all loci were polymorphic with 3–12 alleles. The overall observed (Ho) and unbiased expected (uHE) heterozygosities were 0.60–0.66 and 0.69–0.71, respectively. Parental and sibling relationships were also found among Hungarian individuals: three generations of a pack in the Bükk Mountains were identified. Samples from free-ranging wolves clustered separately from those of captive wolves and dogs. However, genetic similarities were found between Slovakian and Hungarian wolf samples. Our analyses indicate a Slovakian origin of the sampled Hungarian wolves, and we found no evidence that individuals originating in captivity have played any role in the recolonisation process. Kinship relationships and moderate genetic diversity suggest that there is ongoing gene flow across the Slovakian–Hungarian border.

The Diversity in the Genus Canis Challenges Conservation Biology: A Review of Available Data on Asian Wolves

Taxa belonging to the Genus Canis can challenge taxonomists because species boundaries and distribution ranges are often gradual. Species delineation within Canis is currently not based on consistent criteria, and is hampered by geographical bias and lack of taxonomic research. But a consistent taxonomy is critical, given its importance for assigning legal protection, conservation priorities, and financial resources. We carried out a qualitative review of the major wolf lineages so far identified from Asia from historical to contemporary time and considered relevant morphological, ecological, and genetic evidence. We present full mitochondrial phylogenies and genetic distances between these lineages. This review aims to summarize the available data on contemporary Asian wolf lineages within the context of the larger phylogenetic Canis group and to work toward a taxonomy that is consistent within the Canidae. We found support for the presence and taxon eligibility of Holarctic gray, Himalayan/Tibetan, Indian, and Arabian wolves in Asia and recommend their recognition at the taxonomic levels consistent within the group.

A reduced SNP panel to trace gene flow across southern European wolf populations and detect hybridization with other Canis taxa

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ Q_i ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Q_i ≤ 0.827 to parental groups). An assay design and test with invasive and non-invasive canid samples performed successfully for 178 SNPs. By separating natural population admixture from inter-specific hybridization, our reduced panel can help advance evolutionary research, monitoring, and timely conservation management.

Population genetic structure of wolves in the northwestern Dinaric-Balkan region

The Balkan Peninsula and the Dinaric Mountains possess extraordinary biodiversity and support one of the largest and most diverse wolf (Canis lupus) populations in Europe. Results obtained with diverse genetic markers show west-east substructure, also seen in various other species, despite the absence of obvious barriers to movement. However, the spatial extent of the genetic clusters remains unresolved, and our aim was to combine fine-scale sampling with population and spatial genetic analyses to improve resolution of wolf genetic clusters. We analyzed 16 autosomal microsatellites from 255 wolves sampled in Slovenia, Croatia, Bosnia and Herzegovina (BIH), and Serbia and documented three genetic clusters. These comprised (1) Slovenia and the regions of Gorski kotar and Lika in Croatia, (2) the region of Dalmatia in southern Croatia and BIH, and (3) Serbia. When we mapped the clusters geographically, we observed west-east genetic structure across the study area, together with some specific structure in BIH-Dalmatia. We observed that cluster 1 had a smaller effective population size, consistent with earlier reports of population recovery since the 1980s. Our results provide foundation for future genomic studies that would further resolve the observed west-east population structure and its evolutionary history in wolves and other taxa in the region and identify focal areas for habitat conservation. They also have immediate importance for conservation planning for the wolves in one of the most important parts of the species' European range.

Disentangling the admixed trails of grey wolf evolution

The predominant phylogenetic patterns within a genome do not always reflect correctly the history of evolutionary divergence and speciation, and the true phylogenetic signal tends to be concentrated within low-recombination regions of the genome. In this issue of Molecular Ecology, Hennelly et al. show that this is also the case for intraspecific relationships that are characterized by considerable gene flow between lineages. The study reconstructs the phylogenetic relationships of Indian and Tibetan wolves with other grey wolf (Canis lupus) populations worldwide, and demonstrates that these two populations represent phylogenetically distinct lineages. This inference was supported by using low-recombination regions of autosomal chromosomes and the X chromosome, which proved to be essential for correct inference of the lineage splitting order. Their study illustrates the power of analytical approaches that implement knowledge of genome evolution patterns to reconstruct complex intraspecific evolutionary relationships. The study also provides a compelling example of the application of modern phylogenomic approaches in the identification of evolutionarily significant units for the purpose of species conservation.

How the west was won: genetic reconstruction of rapid wolf recolonization into Germany's anthropogenic landscapes

Following massive persecution and eradication, strict legal protection facilitated a successful reestablishment of wolf packs in Germany, which has been ongoing since 2000. Here, we describe this recolonization process by mitochondrial DNA control-region sequencing, microsatellite genotyping and sex identification based on 1341 mostly non-invasively collected samples. We reconstructed the genealogy of German wolf packs between 2005 and 2015 to provide information on trends in genetic diversity, dispersal patterns and pack dynamics during the early expansion process. Our results indicate signs of a founder effect at the start of the recolonization. Genetic diversity in German wolves is moderate compared to other European wolf populations. Although dispersal among packs is male-biased in the sense that females are more philopatric, dispersal distances are similar between males and females once only dispersers are accounted for. Breeding with close relatives is regular and none of the six male wolves originating from the Italian/Alpine population reproduced. However, moderate genetic diversity and inbreeding levels of the recolonizing population are preserved by high sociality, dispersal among packs and several immigration events. Our results demonstrate an ongoing, rapid and natural wolf population expansion in an intensively used cultural landscape in Central Europe.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples

Background

Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples.

Results

Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations.

Conclusions

We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07761-5.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Ancient DNA from the Asiatic Wild Dog (Cuon alpinus) from Europe

The Asiatic wild dog (Cuon alpinus), restricted today largely to South and Southeast Asia, was widespread throughout Eurasia and even reached North America during the Pleistocene. Like many other species, it suffered from a huge range loss towards the end of the Pleistocene and went extinct in most of its former distribution. The fossil record of the dhole is scattered and the identification of fossils can be complicated by an overlap in size and a high morphological similarity between dholes and other canid species. We generated almost complete mitochondrial genomes for six putative dhole fossils from Europe. By using three lines of evidence, i.e., the number of reads mapping to various canid mitochondrial genomes, the evaluation and quantification of the mapping evenness along the reference genomes and phylogenetic analysis, we were able to identify two out of six samples as dhole, whereas four samples represent wolf fossils. This highlights the contribution genetic data can make when trying to identify the species affiliation of fossil specimens. The ancient dhole sequences are highly divergent when compared to modern dhole sequences, but the scarcity of dhole data for comparison impedes a more extensive analysis.

First wolves in Luxembourg since 1893, originating from the Alpine and Central European populations

Considering the spread of the grey wolf (Canis lupus) in Europe over the past 30 years, Luxembourg took some measures to prepare for the return of this apex predator, including the establishment of a management plan that notably addresses the issue of wolf depredation on livestock. Here we present the results of genetic analyses of putative wolf saliva, hair and scat samples collected from or near prey carcasses between 2015 and 2020. In two cases, the wolf was confirmed via DNA analysis: in July 2017 near Garnich and in April 2020 near Niederanven, both assigned to category C1 (hard evidence). A third case was classified as C2 (confirmed observation) based on prey carcass characteristics, while genetic analysis yielded no result. These are the first confirmed records of wolves in Luxemburg since 1893. Moreover, the two C1-cases originated from the Alpine (Garnich) and Central European (Niederanven) populations. Given similar developments in the neighboring countries and regions, we conclude that the area including the Benelux countries as well as Rhineland-Palatinate, Saarland and Northern France may well become a melting pot for wolves of the two aforementioned populations in the coming years and decades.

Population Dynamics in Italian Canids between the Late Pleistocene and Bronze Age

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a ~24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

Development of D-Loop mitochondrial markers for amplification of prey DNA from wolf scat

Analysis of wolves dietary is a currently important theme because of the discussion about wolves preying on livestock as sheep or goats. We developed molecular markers to especially amplify the DNA of the prey out of wolf scat. For this purpose, we used the mitochondrial D-Loop using public available sequences for wolf and seven potential prey species (even-toed ungulates). We developed special primers amplifying either the wolves DNA or the prey DNA. In a fragment of 223-225 basepairs (bp) length we identified 21 SNPs, two 1-bp indels and one 3-bp indel, and three microsatellites to separate seven prey species from each other. Validation of the markers was performed by sequencing the PCR products of 12 fresh prey tissues and 20 wolf scat samples using the different primer pairs.

Genomic evidence for the Old divergence of Southern European wolf populations

The grey wolf (Canis lupus) is one of the most widely distributed mammals in which a variety of distinct populations have been described. However, given their currently fragmented distribution and recent history of human-induced population decline, little is known about the events that led to their differentiation. Based on the analysis of whole canid genomes, we examined the divergence times between Southern European wolf populations and their ancient demographic history. We found that all present-day Eurasian wolves share a common ancestor ca 36 000 years ago, supporting the hypothesis that all extant wolves derive from a single population that subsequently expanded after the Last Glacial Maximum. We also estimated that the currently isolated European populations of the Iberian Peninsula, Italy and the Dinarics-Balkans diverged very closely in time, ca 10 500 years ago, and maintained negligible gene flow ever since. This indicates that the current genetic and morphological distinctiveness of Iberian and Italian wolves can be attributed to their isolation dating back to the end of the Pleistocene, predating the recent human-induced extinction of wolves in Central Europe by several millennia.

Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts

Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.

Dynamic range expansion leads to establishment of a new, genetically distinct wolf population in Central Europe

Local extinction and recolonization events can shape genetic structure of subdivided animal populations. The gray wolf (Canis lupus) was extirpated from most of Europe, but recently recolonized big part of its historical range. An exceptionally dynamic expansion of wolf population is observed in the western part of the Great European Plain. Nonetheless, genetic consequences of this process have not yet been fully understood. We aimed to assess genetic diversity of this recently established wolf population in Western Poland (WPL), determine its origin and provide novel data regarding the population genetic structure of the grey wolf in Central Europe. We utilized both spatially explicit and non-explicit Bayesian clustering approaches, as well as a model-independent, multivariate method DAPC, to infer genetic structure in large dataset (881 identified individuals) of wolf microsatellite genotypes. To put the patterns observed in studied population into a broader biogeographic context we also analyzed a mtDNA control region fragment widely used in previous studies. In comparison to a source population, we found slightly reduced allelic richness and heterozygosity in the newly recolonized areas west of the Vistula river. We discovered relatively strong west-east structuring in lowland wolves, probably reflecting founder-flush and allele surfing during range expansion, resulting in clear distinction of WPL, eastern lowland and Carpathian genetic groups. Interestingly, wolves from recently recolonized mountainous areas (Sudetes Mts, SW Poland) clustered together with lowland, but not Carpathian wolf populations. We also identified an area in Central Poland that seems to be a melting pot of western, lowland eastern and Carpathian wolves. We conclude that the process of dynamic recolonization of Central European lowlands lead to the formation of a new, genetically distinct wolf population. Together with the settlement and establishment of packs in mountains by lowland wolves and vice versa, it suggests that demographic dynamics and possibly anthropogenic barriers rather than ecological factors (e.g. natal habitat-biased dispersal patterns) shape the current wolf genetic structure in Central Europe.

Global Phylogeographic and Admixture Patterns in Grey Wolves and Genetic Legacy of An Ancient Siberian Lineage

The evolutionary relationships between extinct and extant lineages provide important insight into species' response to environmental change. The grey wolf is among the few Holarctic large carnivores that survived the Late Pleistocene megafaunal extinctions, responding to that period's profound environmental changes with loss of distinct lineages and phylogeographic shifts, and undergoing domestication. We reconstructed global genome-wide phylogeographic patterns in modern wolves, including previously underrepresented Siberian wolves, and assessed their evolutionary relationships with a previously genotyped wolf from Taimyr, Siberia, dated at 35 Kya. The inferred phylogeographic structure was affected by admixture with dogs, coyotes and golden jackals, stressing the importance of accounting for this process in phylogeographic studies. The Taimyr lineage was distinct from modern Siberian wolves and constituted a sister lineage of modern Eurasian wolves and domestic dogs, with an ambiguous position relative to North American wolves. We detected gene flow from the Taimyr lineage to Arctic dog breeds, but population clustering methods indicated closer similarity of the Taimyr wolf to modern wolves than dogs, implying complex post-divergence relationships among these lineages. Our study shows that introgression from ecologically diverse con-specific and con-generic populations was common in wolves' evolutionary history, and could have facilitated their adaptation to environmental change.

The Sicilian Wolf: Genetic Identity of a Recently Extinct Insular Population

Historically, many local grey wolf (Canis lupus) populations have undergone substantial reductions in size or become extinct. Among these, the wolf population once living in Sicily, the largest island in the Mediterranean Sea, was completely eradicated by human activity in the early decades of the 20th century. To gain a better understanding of the genetic identity of the Sicilian wolf, we used techniques for the study of ancient DNA to analyze the mitochondrial (mt) variability of six specimens stored in Italian museums. We were able to amplify a diagnostic mtDNA fragment of the control region (CR) in four of the samples. Two of the samples shared the same haplotype, differing by two substitutions from the currently most diffused Italian wolf haplotype (W14) and one substitution from the only other Italian haplotype (W16). The third sample showed a previously unreported wolf-like haplotype, and the fourth a haplotype commonly found in dogs. All of the wolf haplotypes analyzed in this study belonged to the mitochondrial haplogroup that includes haplotypes detected in all the known European Pleistocene wolves and in several modern southern European populations. Unfortunately, this endemic island population, which exhibited unique mtDNA variability, was definitively lost before it was possible to understand its taxonomic uniqueness and conservational value.

The Apennines as a cryptic Pleistocene refugium of the bark beetle Pityogenes chalcographus (Coleoptera: Curculionidae)

The Apennine Mountains in Italy are an important biogeographical region and of particular interest in phylogeographical research, because they have been a refugium during Pleistocene glaciation events for numerous European species. We performed a genetic study on the Eurasian bark beetle Pityogenes chalcographus (Linnaeus, 1760), focusing on two Apennine (Italian) and two Central European (Austrian) locations to assess the influence of the Apennines in the evolutionary history of the beetle, particularly during the Pleistocene. We analysed a part of the mitochondrial COI gene and a set of 5470 informative genome-wide markers to understand its biogeography. We found 75 distinct mitochondrial haplotypes, which are structured in three main clades. In general, the Apennine locations harbour a higher number of mitochondrial clades than Central European sites, with one specific clade exclusively detected in the Apennines. Analysis of our genome-wide, multi-locus dataset reveals a clustering of P. chalcographus by geography, with Italian individuals clearly separated from Austrian samples. Our data highlight the significance of the Apennines for the genetic diversity of P. chalcographus and support the hypothesis that this area was an important refugium during unfavourable conditions in the Pleistocene. We discuss additional life-history traits and processes that shaped the evolution of this widespread beetle.

Winter temperature correlates with mtDNA genetic structure of yellow-necked mouse population in NE Poland

We analysed a fragment (247 bp) of cytochrome b of mitochondrial DNA sequenced using 353 samples of yellow-necked mice Apodemus flavicollis trapped in seven forests and along three woodlot transects in north-eastern Poland. Our aims were to identify the phylogeographic pattern and mtDNA structure of the population and to evaluate the role of environmental conditions in shaping the spatial pattern of mtDNA diversity. We found out that three European haplogroups occurred sympatrically in north-eastern Poland. Inferences based on mtDNA haplotype distribution and frequency defined five subpopulations. The mtDNA-based structure of mice significantly correlated with winter temperature: frequency of Haplogroup 1 was positively, and that of Haplogroup 3 negatively correlated to mean temperature of January in the year of trapping. Synthesis of the published pan-European data on the species phylogeography also showed that the possibly 'thermophilous' Haplogroup 1 has the westernmost occurrence, whereas the more 'cold-resistant' Haplogroup 3 occurs much further to north-east than the other haplogroups. The observed patter may be a byproduct of the tight coevolution with nuclear genes, as we have earlier found that - in mice population in NE Poland - the spatial pattern of nuclear DNA was best explained by January temperature. Alternatively, the observed association of mitochondrial genetic variation with temperature is possible to be adaptive as cytochrome b is involved in the process of ATP production via oxidative phosphorylation.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1420-5) contains supplementary material, which is available to authorized users.

The evolutionary history of grey wolf Y chromosomes

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Old wild wolves: ancient DNA survey unveils population dynamics in Late Pleistocene and Holocene Italian remains

Background

The contemporary Italian wolf (Canis lupus italicus) represents a case of morphological and genetic uniqueness. Today, Italian wolves are also the only documented population to fall exclusively within the mitochondrial haplogroup 2, which was the most diffused across Eurasian and North American wolves during the Late Pleistocene. However, the dynamics leading to such distinctiveness are still debated.

Methods

In order to shed light on the ancient genetic variability of this wolf population and on the origin of its current diversity, we collected 19 Late Pleistocene-Holocene samples from northern Italy, which we analyzed at a short portion of the hypervariable region 1 of the mitochondrial DNA, highly informative for wolf and dog phylogenetic analyses.

Results

Four out of the six detected haplotypes matched the ones found in ancient wolves from northern Europe and Beringia, or in modern European and Chinese wolves, and appeared closely related to the two haplotypes currently found in Italian wolves. The haplotype of two Late Pleistocene samples matched with primitive and contemporary dog sequences from the canine mitochondrial clade A. All these haplotypes belonged to haplogroup 2. The only exception was a Holocene sample dated 3,250 years ago, affiliated to haplogroup 1.

Discussion

In this study we describe the genetic variability of the most ancient wolf specimens from Italy analyzed so far, providing a preliminary overview of the genetic make-up of the population that inhabited this area from the last glacial maximum to the Middle Age period. Our results endorsed that the genetic diversity carried by the Pleistocene wolves here analyzed showed a strong continuity with other northern Eurasian wolf specimens from the same chronological period. Contrarily, the Holocene samples showed a greater similarity only with modern sequences from Europe and Asia, and the occurrence of an haplogroup 1 haplotype allowed to date back previous finding about its presence in this area. Moreover, the unexpected discovery of a 24,700-year-old sample carrying a haplotype that, from the fragment here obtained, falls within the canine clade A, could represent the oldest evidence in Europe of such dog-rich clade. All these findings suggest complex population dynamics that deserve to be further investigated based on mitochondrial or whole genome sequencing.

Two decades of non-invasive genetic monitoring of the grey wolves recolonizing the Alps support very limited dog introgression

Potential hybridization between wolves and dogs has fueled the sensitive conservation and political debate underlying the recovery of the grey wolf throughout Europe. Here we provide the first genetic analysis of wolf-dog admixture in an area entirely recolonized, the northwestern Alps. As part of a long-term monitoring program, we performed genetic screening of thousands of non-invasive samples collected in Switzerland and adjacent territories since the return of the wolf in the mid-1990s. We identified a total of 115 individuals, only 2 of them showing significant signs of admixture stemming from past interbreeding with dogs, followed by backcrossing. This low rate of introgression (<2% accounting for all wolves ever detected over 1998–2017) parallels those from other European populations, especially in Western Europe (<7%). Despite potential hybridization with stray dogs, few founders and strong anthropogenic pressures, the genetic integrity of the Alpine population has remained intact throughout the entire recolonization process. In a context of widespread misinformation, this finding should reduce conflicts among the different actors involved and facilitate wolf conservation. Real-time genetic monitoring will be necessary to identify potential hybrids and support an effective management of this emblematic population.

Howling from the past: historical phylogeography and diversity losses in European grey wolves

Genetic bottlenecks resulting from human-induced population declines make alarming symbols for the irreversible loss of our natural legacy worldwide. The grey wolf (Canis lupus) is an iconic example of extreme declines driven by anthropogenic factors. Here, we assessed the genetic signatures of 150 years of wolf persecution throughout the Western Palaearctic by high-throughput mitochondrial DNA sequencing of historical specimens in an unprecedented spatio-temporal framework. Despite Late Pleistocene bottlenecks, we show that historical genetic variation had remained high throughout Europe until the last several hundred years. In Western Europe, where wolves nearly got fully exterminated, diversity dramatically collapsed at the turn of the twentieth century and recolonization from few homogeneous relict populations induced drastic shifts of genetic composition. By contrast, little genetic displacement and steady levels of diversity were maintained in Eastern European regions, where human persecution had lesser effects on wolf demography. By comparing prehistoric, historic and modern patterns of genetic diversity, our study hence traces the timeframe and the active human role in the decline of the grey wolf, an emblematic yet controversial animal which symbolizes the complex relationship between human societies and nature conservation.

Ancient DNA from Giant Panda (Ailuropoda melanoleuca) of South-Western China Reveals Genetic Diversity Loss during the Holocene

The giant panda was widely distributed in China and south-eastern Asia during the middle to late Pleistocene, prior to its habitat becoming rapidly reduced in the Holocene. While conservation reserves have been established and population numbers of the giant panda have recently increased, the interpretation of its genetic diversity remains controversial. Previous analyses, surprisingly, have indicated relatively high levels of genetic diversity raising issues concerning the efficiency and usefulness of reintroducing individuals from captive populations. However, due to a lack of DNA data from fossil specimens, it is unknown whether genetic diversity was even higher prior to the most recent population decline. We amplified complete cytb and 12s rRNA, partial 16s rRNA and ND1, and control region sequences from the mitochondrial genomes of two Holocene panda specimens. We estimated genetic diversity and population demography by analyzing the ancient mitochondrial DNA sequences alongside those from modern giant pandas, as well as from other members of the bear family (Ursidae). Phylogenetic analyses show that one of the ancient haplotypes is sister to all sampled modern pandas and the second ancient individual is nested among the modern haplotypes, suggesting that genetic diversity may indeed have been higher earlier during the Holocene. Bayesian skyline plot analysis supports this view and indicates a slight decline in female effective population size starting around 6000 years B.P., followed by a recovery around 2000 years ago. Therefore, while the genetic diversity of the giant panda has been affected by recent habitat contraction, it still harbors substantial genetic diversity. Moreover, while its still low population numbers require continued conservation efforts, there seem to be no immediate threats from the perspective of genetic evolutionary potential.

New insights into the genetic composition and phylogenetic relationship of wolves and dogs in the Iberian Peninsula

This study investigates the gene pool of Portuguese autochthonous dog breeds and their wild counterpart, the Iberian wolf subspecies (Canis lupus signatus), using standard molecular markers. A combination of paternal and maternal molecular markers was used to investigate the genetic composition, genetic differentiation and genetic relationship of native Portuguese dogs and the Iberian wolf. A total of 196 unrelated dogs, including breed and village dogs from Portugal, and other dogs from Spain and North Africa, and 56 Iberian wolves (wild and captive) were analyzed for nuclear markers, namely Y chromosome SNPs, Y chromosome STR loci, autosomal STR loci, and a mitochondrial fragment of the control region I. Our data reveal new variants for the molecular markers and confirm significant genetic differentiation between Iberian wolf and native domestic dogs from Portugal. Based on our sampling, no signs of recent introgression between the two subspecies were detected. Y chromosome data do not reveal genetic differentiation among the analyzed dog breeds, suggesting they share the same patrilineal origin. Moreover, the genetic distinctiveness of the Iberian wolf from other wolf populations is further confirmed with the description of new mtDNA variants for this endemism. Our research also discloses new molecular markers for wolf and dog subspecies assignment, which might become particularly relevant in the case of forensic or noninvasive genetic studies. The Iberian wolf represents a relic of the once widespread wolf population in Europe and our study reveals that it is a reservoir of unique genetic diversity of the grey wolf, Canis lupus. These results stress the need for conservation plans that will guarantee the sustainability of this threatened top predator in Iberia.

Phylogenetic evidence for the ancient Himalayan wolf: towards a clarification of its taxonomic status based on genetic sampling from western Nepal

Wolves in the Himalayan region form a monophyletic lineage distinct from the present-day Holarctic grey wolf Canis lupus spp. (Linnaeus 1758) found across Eurasia and North America. Here, we analyse phylogenetic relationships and the geographic distribution of mitochondrial DNA haplotypes of the contemporary Himalayan wolf (proposed in previous studies as Canis himalayensis) found in Central Asia. We combine genetic data from a living Himalayan wolf population collected in northwestern Nepal in this study with already published genetic data, and confirm the Himalayan wolf lineage based on mitochondrial genomic data (508 bp cytochrome b and 242 bp D-loop), and X- and Y-linked zinc-finger protein gene (ZFX and ZFY) sequences. We then compare the genetic profile of the Himalayan wolf lineage found in northwestern Nepal with canid reference sequences from around the globe with maximum likelihood and Bayesian phylogeny building methods to demonstrate that the Himalayan wolf forms a distinct monophyletic clade supported by posterior probabilities/bootstrap for D-loop of greater than 0.92/85 and cytochrome b greater than 0.99/93. The Himalayan wolf shows a unique Y-chromosome (ZFY) haplotype, and shares an X-chromosome haplotype (ZFX) with the newly postulated African wolf. Our results imply that the Himalayan wolf distribution range extends from the Himalayan range north across the Tibetan Plateau up to the Qinghai Lakes region in Qinghai Province in the People's Republic of China. Based on its phylogenetic distinction and its older age of divergence relative to the Holarctic grey wolf, the Himalayan wolf merits formal classification as a distinct taxon of special conservation concern.

Combining phylogenetic and demographic inferences to assess the origin of the genetic diversity in an isolated wolf population

The survival of isolated small populations is threatened by both demographic and genetic factors. Large carnivores declined for centuries in most of Europe due to habitat changes, overhunting of their natural prey and direct persecution. However, the current rewilding trends are driving many carnivore populations to expand again, possibly reverting the erosion of their genetic diversity. In this study we reassessed the extent and origin of the genetic variation of the Italian wolf population, which is expanding after centuries of decline and isolation. We genotyped wolves from Italy and other nine populations at four mtDNA regions (control-region, ATP6, COIII and ND4) and 39 autosomal microsatellites. Results of phylogenetic analyses and assignment procedures confirmed in the Italian wolves a second private mtDNA haplotype, which belongs to a haplogroup distributed mostly in southern Europe. Coalescent analyses showed that the unique mtDNA haplotypes in the Italian wolves likely originated during the late Pleistocene. ABC simulations concordantly showed that the extant wolf populations in Italy and in south-western Europe started to be isolated and declined right after the last glacial maximum. Thus, the standing genetic variation in the Italian wolves principally results from the historical isolation south of the Alps.