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

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus)

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

Application of real-time PCR for the identification of the endangered species Galemys pyrenaicus through faecal samples

BACKGROUND

Currently, many micromammals are important targets for study. The endangered Galemys pyrenaicus is an outstanding example. Globally, their populations have suffered a substantial decline in last 20 years. In the surveyed area, the capture of desman is legally forbidden due to the high conservation concerns. Reason by non-invasive sampling through faeces is proposed for its monitoring. Furthermore, the confusion between faeces from desman and Mediterranean water shrews must be considered. Thus, the aim of this study was focused on developing RT-PCR assays to determine the presence of Galemys pyrenaicus and N. a. anomalus from non-invasive samples.

METHODS AND RESULTS

The study was conducted in the mountains of the System Central of Extremadura (Spain). A total of 186 samples were collected from 2018 to 2021 by experts where historically reported and/or our previous studies confirmed their presence. RT-PCR assays using hydrolysis probes were designed to detect genetic material from both desman and Mediterranean water shrews and its specificity was confirmed. The reliability of the method was further assessed by PCR sequencing of mitochondrial Cyb and d-loop, resulting fully compatible with the RT-PCR approach. Intraspecific phylogenetic relationship was reported to improve knowledge about mtDNA variability in the desman from the Central System.

CONCLUSIONS

We demonstrated that RT-PCR gives a gold opportunity to further map the species using faeces which minimizes disturbance and reports both population status and individual presence. Cost-effective RT-PCR combined with field-collected faeces allows us to better investigate the full range of occurrence of the species.

Non-Invasive Wildlife Disease Surveillance Using Real Time PCR Assays: The Case of the Endangered Galemys pyrenaicus Populations from the Central System Mountains (Extremadura, Spain)

The Iberian desman (Galemys pyrenaicus) is a small semi-aquatic mammal that inhabits mountainous areas from the centre to the north of the Iberian Peninsula and the Pyrenees and is listed as endangered because it has suffered a serious decline. Since 1960, only three species of digeneans (Omphalometra flexuosa, Maritrema pyrenaica and Mathovius galemydis) and two nematodes (Aonchotheca galemydis and Paracuaria hispanica) have been reported from the desman, but no further information on health status and no data from Extremadura has been available. The aim of our study was to characterise the diversity and distribution of parasites and microbiomes of desmans in different areas of the Central System of Extremadura. Between 2019 and 2021 we collected 238 fecal samples and one tissue (intestine) sample that was obtained from a dead desman. DNA templates were processed by commercial or customised real-time PCR using TaqMan probes. Representative data were obtained for Cryptosporidium spp., Omphalometra spp., Eimeria spp., Salmonella spp., Staphylococcus spp. and Leptospira spp. Omphalometra spp. was studied using a newly developed PCR test. The screening of the dead desman allowed us to obtain, for the first time, a partial sequence of the 18SrDNA. This study is the most complete study of the desman, allowing us to identify parasites and the microbiome in populations of G. pyrenaicus using non-invasive sampling.

Koala Genome Survey: An Open Data Resource to Improve Conservation Planning

Genome sequencing is a powerful tool that can inform the management of threatened species. Koalas (Phascolarctos cinereus) are a globally recognized species that captured the hearts and minds of the world during the 2019/2020 Australian megafires. In 2022, koalas were listed as ‘Endangered’ in Queensland, New South Wales, and the Australian Capital Territory. Populations have declined because of various threats such as land clearing, habitat fragmentation, and disease, all of which are exacerbated by climate change. Here, we present the Koala Genome Survey, an open data resource that was developed after the Australian megafires. A systematic review conducted in 2020 demonstrated that our understanding of genomic diversity within koala populations was scant, with only a handful of SNP studies conducted. Interrogating data showed that only 6 of 49 New South Wales areas of regional koala significance had meaningful genome-wide data, with only 7 locations in Queensland with SNP data and 4 locations in Victoria. In 2021, we launched the Koala Genome Survey to generate resequenced genomes across the Australian east coast. We have publicly released 430 koala genomes (average coverage: 32.25X, range: 11.3–66.8X) on the Amazon Web Services Open Data platform to accelerate research that can inform current and future conservation planning.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Bactrocera dorsalis in the Indian Ocean: A tale of two invasions

An increasing number of invasive fruit fly pests are colonizing new grounds. With this study, we aimed to uncover the invasion pathways of the oriental fruit fly, Bactrocera dorsalis into the islands of the Indian Ocean. By using genome-wide SNP data and a multipronged approach consisting of PCA, ancestry analysis, phylogenetic inference, and kinship networks, we were able to resolve two independent invasion pathways. A western invasion pathway involved the stepping-stone migration of B. dorsalis from the east African coast into the Comoros, along Mayotte and into Madagascar with a decreasing genetic diversity. The Mascarene islands (Reunion and Mauritius), on the contrary, were colonized directly from Asia and formed a distinct cluster. The low nucleotide diversity suggests that only a few genotypes invaded the Mascarenes. The presence of many long runs of homozygosity (ROH) in the introduced populations is indicative of population bottlenecks, with evidence of a more severe bottleneck for populations along the western migration pathway than on the Mascarene islands. More strict phytosanitary regulations are recommended in order to prevent the further spread of B. dorsalis.

Satellitome Analysis on Talpa aquitania Genome and Inferences about the satDNAs Evolution on Some Talpidae

In the genus Talpa a new species, named Talpa aquitania, has been recently described. Only cytogenetic data are available for the nuclear genome of this species. In this work, we characterize the satellitome of the T. aquitania genome that presents 16 different families, including telomeric sequences, and they represent 1.24% of the genome. The first satellite DNA family (TaquSat1-183) represents 0.558%, and six more abundant families, including TaquSat1-183, comprise 1.13%, while the remaining 11 sat-DNAs represent only 0.11%. The average A + T content of the SatDNA families was 50.43% and the median monomer length was 289.24 bp. The analysis of these SatDNAs indicated that they have different grades of clusterization, homogenization, and degeneration. Most of the satDNA families are present in the genomes of the other Talpa species analyzed, while in the genomes of other more distant species of Talpidae, only some of them are present, in accordance with the library hypothesis. Moreover, chromosomal localization by FISH revealed that some satDNAs are localized preferentially on centromeric and non-centromeric heterochromatin in T. aquitania and also in the sister species T. occidentalis karyotype. The differences observed between T. aquitania and the close relative T. occidentalis and T. europaea suggested that the satellitome is a very dynamic component of the genomes and that the satDNAs could be responsible for chromosomal differences between the species. Finally, in a broad context, these data contribute to the understanding of the evolution of satellitomes on mammals.

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.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Divergence time estimation using ddRAD data and an isolation-with-migration model applied to water vole populations of Arvicola

Molecular dating methods of population splits are crucial in evolutionary biology, but they present important difficulties due to the complexity of the genealogical relationships of genes and past migrations between populations. Using the double digest restriction-site associated DNA (ddRAD) technique and an isolation-with-migration (IM) model, we studied the evolutionary history of water vole populations of the genus Arvicola, a group of complex evolution with fossorial and semi-aquatic ecotypes. To do this, we first estimated mutation rates of ddRAD loci using a phylogenetic approach. An IM model was then used to estimate split times and other relevant demographic parameters. A set of 300 ddRAD loci that included 85 calibrated loci resulted in good mixing and model convergence. The results showed that the two populations of A. scherman present in the Iberian Peninsula split 34 thousand years ago, during the last glaciation. In addition, the much greater divergence from its sister species, A. amphibius, may help to clarify the controversial taxonomy of the genus. We conclude that this approach, based on ddRAD data and an IM model, is highly useful for analyzing the origin of populations and species.

Ancient demographics determine the effectiveness of genetic purging in endangered lizards

The purging of deleterious alleles has been hypothesized to mitigate inbreeding depression, but its effectiveness in endangered species remains debatable. To understand how deleterious alleles are purged during population contractions, we analyzed genomes of the endangered Chinese crocodile lizard (Shinisaurus crocodilurus), which is the only surviving species of its family and currently isolated into small populations. Population genomic analyses revealed four genetically distinct conservation units and sharp declines in both effective population size and genetic diversity. By comparing the relative genetic load across populations and conducting genomic simulations, we discovered that seriously deleterious alleles were effectively purged during population contractions in this relict species, although inbreeding generally enhanced the genetic burden. However, despite with the initial purging, our simulations also predicted that seriously deleterious alleles will gradually accumulate under prolonged bottlenecking. Therefore, we emphasize the importance of maintaining a minimum population capacity and increasing the functional genetic diversity in conservation efforts to preserve populations of the crocodile lizard and other endangered species.