Evaluating the use of non-invasive hair sampling and ddRAD to characterize populations of endangered species: Application to a peripheral population of the European mink.

The application of next-generation sequencing (NGS) to non-invasive samples is one of the most promising methods in conservation genomics, but these types of samples present significant challenges for NGS. The European mink (Mustela lutreola) is critically endangered throughout its range. However, important aspects such as census size and inbreeding remain still unknown in many populations, so it is crucial to develop new methods to monitor this species. In this work, we placed hair tubes along riverbanks in a border area of the Iberian population, which allowed the genetic identification of 76 European mink hair samples. We then applied a reduced representation genomic sequencing (ddRAD) technique to a subset of these samples to test whether we could extract sufficient genomic information from them. We show that several problems with the DNA, including contamination, fragmentation, oxidation, and possibly sample mixing, affected the samples. Using various bioinformatic techniques to reduce these problems, we were able to unambiguously genotype 19 hair samples belonging to six individuals. This small number of individuals showed that the demographic status of the species in this peripheral population is worse than expected. The data obtained also allowed us to perform preliminary analyses of relatedness and inbreeding. Although further improvements in sampling and analysis are needed, the application of the ddRAD technique to non-invasively obtained hairs represents a significant advance in the genomic study of endangered species.

Methodological challenges in the genomic analysis of an endangered mammal population with low genetic diversity.

Recently, populations of various species with very low genetic diversity have been discovered. Some of these persist in the long term, but others could face extinction due to accelerated loss of fitness. In this work, we characterize 45 individuals of one of these populations, belonging to the Iberian desman (Galemys pyrenaicus). For this, we used the ddRADseq technique, which generated 1421 SNPs. The heterozygosity values of the analyzed individuals were among the lowest recorded for mammals, ranging from 26 to 91 SNPs/Mb. Furthermore, the individuals from one of the localities, highly isolated due to strong barriers, presented extremely high inbreeding coefficients, with values above 0.7. Under this scenario of low genetic diversity and elevated inbreeding levels, some individuals appeared to be almost genetically identical. We used different methods and simulations to determine if genetic identification and parentage analysis were possible in this population. Only one of the methods, which does not assume population homogeneity, was able to identify all individuals correctly. Therefore, genetically impoverished populations pose a great methodological challenge for their genetic study. However, these populations are of primary scientific and conservation interest, so it is essential to characterize them genetically and improve genomic methodologies for their research.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species.

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.

Publisher Correction: Size increase without genetic divergence in the Eurasian water shrew Neomys fodiens.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Size increase without genetic divergence in the Eurasian water shrew Neomys fodiens.

When a population shows a marked morphological change, it is important to know whether that population is genetically distinct; if it is not, the novel trait could correspond to an adaptation that might be of great ecological interest. Here, we studied a subspecies of water shrew, Neomys fodiens niethammeri, which is found in a narrow strip of the northern Iberian Peninsula. This subspecies presents an abrupt increase in skull size when compared to the rest of the Eurasian population, which has led to the suggestion that it is actually a different species. Skulls obtained from owl pellets collected over the last 50 years allowed us to perform a morphometric analysis in addition to an extensive multilocus analysis based on short intron fragments successfully amplified from these degraded samples. Interestingly, no genetic divergence was detected using either mitochondrial or nuclear data. Additionally, an allele frequency analysis revealed no significant genetic differentiation. The absence of genetic divergence and differentiation revealed here indicate that the large form of N. fodiens does not correspond to a different species and instead represents an extreme case of size increase, of possible adaptive value, which deserves further investigation.

Quantitative analysis of connectivity in populations of a semi-aquatic mammal using kinship categories and network assortativity.

Analysing the impact of anthropogenic and natural river barriers on the dispersal of aquatic and semi-aquatic species may be critical for their conservation. Knowledge of kinship relationships between individuals and reconstructions of pedigrees obtained using genomic data can be extremely useful, not only for studying the social organization of animals, but also inferring contemporary dispersal and quantifying the effect of specific barriers on current connectivity. In this study, we used kinship data to analyse connectivity patterns in a small semi-aquatic mammal, the Pyrenean desman (Galemys pyrenaicus), in an area comprising two river systems with close headwaters and dams of various heights and types. Using a large SNP dataset from 70 specimens, we obtained kinship categories and reconstructed pedigrees. To quantify the barrier effect of specific obstacles, we built kinship networks and devised a method based on the assortativity coefficient, which measures the proportion between observed and expected kinship relationships across a barrier. The estimation of this parameter enabled us to infer that the most important barrier in the area was the watershed divide between the rivers, followed by a dam on one of the rivers. Other barriers did not significantly reduce the expected number of kinship relationships across them. This strategy and the information obtained with it may be crucial in determining the most important connectivity problems in an area and help develop conservation plans aimed at improving genetic exchange between populations of threatened species.

Using relatedness networks to infer contemporary dispersal: Application to the endangered mammal Galemys pyrenaicus.

Information about the degree of contemporary dispersal is important when trying to understand how populations interchange individuals and identify the specific barriers that prevent these movements. In the case of endangered species, this can represent crucial information when designing appropriate conservation strategies. Here we analyse relatedness between individuals from different localities and use these data to infer whether dispersal occurred in recent generations. We applied this approach to the Pyrenean desman (Galemys pyrenaicus), a semiaquatic and endangered species endemic to the Iberian Peninsula. We studied this species in four primary rivers in the Iberian Range, where two ancient mitochondrial lineages are separated by a strict contact zone, suggesting the existence of complex dispersal patterns. Using next-generation sequencing, we obtained 912 SNPs from each specimen and estimated relatedness values between them. While relatedness networks were dense within each river, we found surprisingly few relationships between individuals from different rivers despite their close proximity in some cases, indicating much lower dispersal between rivers compared to dispersal within a single river. In agreement with this result, the degree of inbreeding was exceedingly high in most individuals. These data show that relatedness information can be crucial to understand the contemporary dispersal patterns and conservation status of specific populations of endangered species.