Genetic diversity and landscape genetic structure of otter (Lutra lutra) populations in Europe

Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods

Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.

Faecal DNA-based genetic survey of a relict Eurasian otter (Lutra lutra) population (Sila Massif, S Italy)

Faecal DNA-based genetic analysis is a suitable tool for assessing both population size and genetic diversity of threatened and elusive species. We applied microsatellite analysis and mtDNA sequencing for investigating the southernmost Italian (Sila Massif, Calabria Region) population of the Eurasian otter (Lutra lutra). This relict population, filed as extinct in the mid-1980s, is currently expanding but still quite isolated. On the two main rivers hosting otters permanently since 2014, we collected 47 spraints, out of which 24 (51.1%) were successfully genotyped (on average 2.0 alleles per locus). Thirteen individuals were identified: seven females and three males (sex identification success of 76.9%). Population size was assessed as 16 individuals (13–22), corresponding to a density of 0.15 (0.13–0.21) ind/km. Successfully amplified mtDNA samples (N = 16) confirmed the occurrence of a haplotype—H10—which had been previously reported only for Southern Italy, bringing new evidence of the unicity of the Italian otter population. Although density values complied with those reported for the core area of otter Italian range, the small size and genetic isolation of this population require special attention. To assist the ongoing re-colonisation of the Sila Massif, habitat management should aim to enhance fish availability and connectivity with the core area.

Health and Mortality Monitoring in Threatened Mammals: A First Post Mortem Study of Otters (Lutra lutra L.) in Italy

Dead specimens provide valuable data for the conservation of threatened species, allowing investigations of mortality, health conditions, and demographic parameters. The Eurasian otter (Lutra lutra) is a semiaquatic carnivore listed as endangered in Italy. In 2009, we started the first post mortem (PM) study of otters in Italy, through collaborative research between mammal ecologists and veterinary pathologists, using standardized protocols. Twenty-eight otters, mostly collected between 2009 and 2017, were examined. Most otters were males (67%), between 1 and 3 years old (64%), and predominantly in good nutritional condition. Adult males were significantly larger than adult females (p < 0.02), as expected for the species, although both sexes appeared to be smaller than otters examined in Central−northern Europe. The youngest sexually mature female was 3 years old. Road traffic collisions were the major cause of death, especially in young individuals, and mainly occurred in autumn−winter, particularly for females. Investigations of the scene of death contributed to revealing factors forcing otters to travel out of the water and move over the road, suggesting appropriate measures to reduce vehicle collision risk. Other causes of death included blunt chest trauma of uncertain origin, dog and conspecific attacks, or diseases of infectious or non-infectious origin, such as ulcerative gastritis, pleuropneumonia and peritonitis. Other diagnosed diseases included lymphoma. Ecto- and endoparasites were rarely detected, although we report the first documentation of heartworm and Ixodes hexagonus infestation in Italian otters. It is important to continue comprehensive, standardized PM investigations of otters in Italy to define baseline health, biometric and demographic parameters, collect biological samples for comparative analyses, and to reduce road-kill mortality. The present study suggests that the timely collection of carcasses and collaborative and coordinated research efforts are essential for obtaining useful data for the conservation of otters.

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 and viability assessment of a reintroduced Eurasian otter Lutra lutra population on the River Ticino, Italy

On the River Ticino in northern Italy, a small number of captive Eurasian otters Lutra lutra, belonging to the European breeding programme for self-sustaining captive populations, were reintroduced in 1997, after the species had been declared locally extinct in the 1980s. We surveyed for otter signs in 2008, 2010, 2016–2017 and 2018, confirming the presence of what is probably a small population. To assess the abundance and viability of the population, we genotyped fresh spraints collected during the last two surveys, using 11 microsatellite markers, and modelled the population trend using Vortex. A minimum of six individuals were identified from 25 faecal samples. The analysis of mitochondrial DNA determined that the reintroduced otters share a transversion that is characteristic of the Asiatic subspecies Lutra lutra barang, confirming the contribution of the Asiatic subspecies to the genetic pool of the captive-bred founder population. Population size was consistent with the release of three pairs of otters and all models implied that the number of founders was too small to ensure the long-term survival of the population. Stochastic factors are therefore likely to threaten the success of this reintroduction.

Fifty Years of Research on European Mink Mustela lutreola L., 1761 Genetics: Where Are We Now in Studies on One of the Most Endangered Mammals?

The purpose of this review is to present the current state of knowledge about the genetics of European mink Mustela lutreola L., 1761, which is one of the most endangered mammalian species in the world. This article provides a comprehensive description of the studies undertaken over the last 50 years in terms of cytogenetics, molecular genetics, genomics (including mitogenomics), population genetics of wild populations and captive stocks, phylogenetics, phylogeography, and applied genetics (including identification by genetic methods, molecular ecology, and conservation genetics). An extensive and up-to-date review and critical analysis of the available specialist literature on the topic is provided, with special reference to conservation genetics. Unresolved issues are also described, such as the standard karyotype, systematic position, and whole-genome sequencing, and hotly debated issues are addressed, like the origin of the Southwestern population of the European mink and management approaches of the most distinct populations of the species. Finally, the most urgent directions of future research, based on the research questions arising from completed studies and the implementation of conservation measures to save and restore M. lutreola populations, are outlined. The importance of the popularization of research topics related to European mink genetics among scientists is highlighted.

Initial Investigations of the Cranial Size and Shape of Adult Eurasian Otters (Lutra lutra) in Great Britain

Three-dimensional (3D) surface scans were carried out in order to determine the shapes of the upper sections of (skeletal) crania of adult Eurasian otters (Lutra lutra) from Great Britain. Landmark points were placed on these shapes using a graphical user interface (GUI) and distance measurements (i.e., the length, height, and width of the crania) were found by using the landmark points. Male otters had significantly larger skulls than females (P < 0.001). Differences in size also occurred by geographical area in Great Britain (P < 0.05). Multilevel Principal Components Analysis (mPCA) indicated that sex and geographical area explained 31.1% and 9.6% of shape variation in "unscaled" shape data and that they explained 17.2% and 9.7% of variation in "scaled" data. The first mode of variation at level 1 (sex) correctly reflected size changes between males and females for "unscaled" shape data. Modes at level 2 (geographical area) also showed possible changes in size and shape. Clustering by sex and geographical area was observed in standardized component scores. Such clustering in a cranial shape by geographical area might reflect genetic differences in otter populations in Great Britain, although other potentially confounding factors (e.g., population age-structure, diet, etc.) might also drive regional differences. This work provides a successful first test of the effectiveness of 3D surface scans and multivariate methods, such as mPCA, to study the cranial morphology of otters.

First telemetry data on wild individuals from the threatened, isolated Italian otter (Lutra lutra) population

The Eurasian otter (Lutra lutra) is among the most endangered mammals in Italy, its population being isolated and genetically distinct from other European populations. Yet, its ecology and behavior are largely unknown in Italy, preventing its much-needed effective conservation. We radio-tracked one male and one female otter daily in Southern Italy, for 1 year and 8 months, respectively. The male and female used nearly 30 and 20 km of river, respectively. The extension of the ranges varied among seasons. The male was more conservative in its spatial utilization, being mostly located around a trout farm. These first direct data on space use by wild Italian otters provided preliminary information potentially useful for a sound management of this endangered population.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Genetic population structure of fishers (Pekania pennanti) in the Great Lakes region: remnants and reintroductions

Reintroduction programs have been pivotal in augmenting populations of fishers (Pekania pennanti (Erxleben, 1777)) and re-establishing them to their former range in North America. The majority of reintroduction efforts in fishers have been considered demographically successful, but reintroductions can alter genetic population structure and success has rarely been evaluated in fishers from a genetic standpoint. We used microsatellite data (n = 169) to examine genetic population structure of fishers in the Great Lakes region and comment on the success of past reintroductions at two different spatial scales. We found significant genetic population structure among source and reintroduced populations within the Great Lakes region and large-scale genetic structure between fisher populations located in two geographically distant regions (Great Lakes and Northeast) in the eastern United States. Reintroductions associated with the Great Lakes produced results that were largely consistent with other studies of fisher reintroductions in the Northeast. However, our data are the first to support a measurable impact on genetic population structure in Pekania pennanti pennanti (Erxleben, 1777) from a reintroduction using geographically distant source and reintroduced populations. When feasible, we strongly recommend that reintroduction programs include an investigation of the underlying genetic structure to better define intended goals and supplement measures of demographic success.

The effect of a major drainage divide on the gene flow of a semiaquatic carnivore, the Eurasian otter

Major drainage divides (separating the water flow draining to different seas) may significantly affect the dispersal, distribution, and genetic structure of semiaquatic animals. Assessing this effect is important for the proper management of both endangered and invasive riparian species. Here, I determined the fine-scale impact of the drainage divide between the Baltic Sea and the Black Sea drainage basins on the genetic structure and dispersal of the Eurasian otter (Lutra lutra) in the Bieszczady Mountains (Poland and Slovakia). I investigated the genetic structure of the otter population and assessed 3 alternative dispersal models using 6 landscape genetics methods. The analyses were based on 48 individual genotypes obtained from 622 fecal and gland secretion samples collected in 2008–2011. Results indicate that the major drainage divide is not a barrier to gene flow for this population. This was established by analyses of population genetic structure and confirmed by analysis of the spatial distribution of samples originating from closely related individuals. In line with these findings, the best-supported dispersal model assumed that otters migrate through mountain passes, away from streams, thus revealing that they are able to cross a drainage divide. The genetic structure of the population studied exhibits an isolation-by-distance pattern; however, the locations of several repeatedly recorded or closely related individuals revealed the occurrence of long-distance movements. Confirmation of the high mobility of otters and their ability to cross a major drainage divide ridge suggest that landscape obstacles and discontinuity of river network are unlikely to stop otter dispersal, gene flow, and recolonization of new areas.

Assessing Genetic Structure in Common but Ecologically Distinct Carnivores: The Stone Marten and Red Fox

The identification of populations and spatial genetic patterns is important for ecological and conservation research, and spatially explicit individual-based methods have been recognised as powerful tools in this context. Mammalian carnivores are intrinsically vulnerable to habitat fragmentation but not much is known about the genetic consequences of fragmentation in common species. Stone martens (Martes foina) and red foxes (Vulpes vulpes) share a widespread Palearctic distribution and are considered habitat generalists, but in the Iberian Peninsula stone martens tend to occur in higher quality habitats. We compared their genetic structure in Portugal to see if they are consistent with their differences in ecological plasticity, and also to illustrate an approach to explicitly delineate the spatial boundaries of consistently identified genetic units. We analysed microsatellite data using spatial Bayesian clustering methods (implemented in the software BAPS, GENELAND and TESS), a progressive partitioning approach and a multivariate technique (Spatial Principal Components Analysis-sPCA). Three consensus Bayesian clusters were identified for the stone marten. No consensus was achieved for the red fox, but one cluster was the most probable clustering solution. Progressive partitioning and sPCA suggested additional clusters in the stone marten but they were not consistent among methods and were geographically incoherent. The contrasting results between the two species are consistent with the literature reporting stricter ecological requirements of the stone marten in the Iberian Peninsula. The observed genetic structure in the stone marten may have been influenced by landscape features, particularly rivers, and fragmentation. We suggest that an approach based on a consensus clustering solution of multiple different algorithms may provide an objective and effective means to delineate potential boundaries of inferred subpopulations. sPCA and progressive partitioning offer further verification of possible population structure and may be useful for revealing cryptic spatial genetic patterns worth further investigation.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

Patterns of genetic variation in the endangered European mink (Mustela lutreola L., 1761)

BACKGROUND

The European mink (Mustela lutreola, L. 1761) is a critically endangered mustelid, which inhabits several main river drainages in Europe. Here, we assess the genetic variation of existing populations of this species, including new sampling sites and additional molecular markers (newly developed microsatellite loci specific to European mink) as compared to previous studies. Probabilistic analyses were used to examine genetic structure within and between existing populations, and to infer phylogeographic processes and past demography.

RESULTS

According to both mitochondrial and nuclear microsatellite markers, Northeastern (Russia, Estonia and Belarus) and Southeastern (Romania) European populations showed the highest intraspecific diversity. In contrast, Western European (France and Spain) populations were the least polymorphic, featuring a unique mitochondrial DNA haplotype. The high differentiation values detected between Eastern and Western European populations could be the result of genetic drift in the latter due to population isolation and reduction. Genetic differences among populations were further supported by Bayesian clustering and two main groups were confirmed (Eastern vs. Western Europe) along with two contained subgroups at a more local scale (Northeastern vs. Southeastern Europe; France vs. Spain).

CONCLUSIONS

Genetic data and performed analyses support a historical scenario of stable European mink populations, not affected by Quaternary climate oscillations in the Late Pleistocene, and posterior expansion events following river connections in both North- and Southeastern European populations. This suggests an eastern refuge during glacial maxima (as already proposed for boreal and continental species). In contrast, Western Europe was colonised more recently following either natural expansions or putative human introductions. Low levels of genetic diversity observed within each studied population suggest recent bottleneck events and stress the urgent need for conservation measures to counteract the demographic decline experienced by the European mink.

Impact of population expansion on genetic diversity and structure of river otters (Lontra canadensis) in Central North America

Populations of North American river otters (Lontra canadensis) declined throughout large portions of the continent during the early 1900s due to habitat degradation and unregulated trapping. River otters had been extirpated in North Dakota (ND), but the Red River Valley has since been recolonized, with potential source populations including the neighboring states of Minnesota or South Dakota, or the Canadian province of Manitoba (MB). We genotyped 9 microsatellite loci in 121 samples to determine the source population of river otters in the Red River Valley of ND, as well as to assess population structure and diversity of river otters in central North America. Overall, genetic diversity was high, with an average observed heterozygosity of 0.58. Genetic differentiation was low (F ST < 0.05) between river otters in ND and those of Minnesota, suggesting that eastern ND was recolonized by river otters from Minnesota. River otters from MB were genetically distinct from all other sampled populations. Low genetic differentiation (F ST = 0.044) between South Dakota and Louisiana (LA) suggested that reintroductions using LA stock were successful. The genetic distinctiveness of river otters from different geographic regions should be considered when deciding on source populations for future translocations.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus

Background

Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia.

Results

To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers.

Conclusions

The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Genetic diversity of the Eurasian Otter (Lutra lutra) population in Israel

The Israeli population of Eurasian otter (Lutra lutra) marks the Palearctic southern boundary of the species' distribution in the Levant. During the 20th century, the otter population in Israel experienced a dramatic decline due to anthropogenic habitat alterations. Currently, the otter population in Israel is estimated at about 100 individuals and defined as "Critically Endangered". The aim of this research was to characterize the Israeli otter population in order to determine its genetic diversity and fragmentation state for conservation purposes. Monitoring spraint sites during 2000-2011 along active and historic otter distribution regions indicate both stable and unstable otter subpopulations, mainly along the Jordan River. Four otter subpopulations, representing 57 individuals, were characterized by 12 microsatellites, previously used to characterize the European otter populations. The genetic results indicated three subpopulations correlating with three geographical regions: the Hula Valley, Sea of Galilee, and the Harod Valley. A moderate genetic diversity (F (st) = 0.087-0.123) was found among the subpopulations, suggesting sporadic interactions between individuals from distinct geographical locations along the Jordan Rift Valley. The Israeli otter population was found to be very small, demographically remote and genetically distinct, harboring unique alleles absent from the studied European populations. Therefore, immediate conservation actions are recommended to prevent the deterioration of the isolated, unique, and critically endangered otter population in Israel.