Genetics at the verge of extinction: insights from the Iberian lynx

Serological and Molecular Survey of Leishmania infantum in a Population of Iberian Lynxes (Lynx pardinus)

Leishmania infantum, the sand fly-transmitted protozoan parasite responsible for leishmaniasis in humans, dogs, and cats, is endemic in the Iberian Peninsula. However, the impact of L. infantum infection on the conservation of the endangered Iberian lynx (Lynx pardinus) is unknown. Herein, we describe for the first time the occurrence of L. infantum infection among a population of reintroduced and wild-born L. pardinus living in the Portuguese Guadiana Valley Park. The presence of infection was addressed by molecular detection of Leishmania kinetoplast DNA (kDNA) in 35 lynxes, with further confirmation of L. infantum species performed by an internally transcribed spacer (ITS)-1 sequencing. Eight blood samples were positive for kDNA, and ITS-1 sequencing confirmed the presence of L. infantum in two of those samples. Exposure to Leishmania was screened in a group of 36 lynxes using an immunofluorescence antibody test (IFAT) and a multi-antigen enzyme-linked immunosorbent assay (ELISA), using SPLA, rK39, and CPX as Leishmania-specific antigens. Four animals presented a positive IFAT at a dilution of 1:40. Eight samples were considered seropositive to all ELISA Leishmania-specific antigens. Agreement between PCR, IFAT, and all ELISA antigens was found for 1 in 27 samples. These results highlight the susceptibility of autochthonous L. pardinus to L. infantum infection. Further investigation is required to assess the impact of L. infantum infection on this wild species conservation.

The first assessment of the genetic diversity and structure of the endangered West Indian manatee in Cuba

The coastal waters of Cuba are home to a small, endangered population of West Indian manatee, which would benefit from a comprehensive characterization of the population's genetic variation. We conducted the first genetic assessment of Cuban manatees to determine the extent of the population's genetic structure and characterize the neutral genetic diversity among regions within the archipelago. We genotyped 49 manatees at 18 microsatellite loci, a subset of 27 samples on 1703 single nucleotide polymorphisms (SNPs), and sequenced 59 manatees at the mitochondrial control region. The Cuba manatee population had low nuclear (microsatellites H_E = 0.44, and SNP H_E = 0.29) and mitochondrial genetic diversity (h = 0.068 and π = 0.00025), and displayed moderate departures from random mating (microsatellite F_IS = 0.12, SNP F_IS = 0.10). Our results suggest that the western portion of the archipelago undergoes periodic exchange of alleles based on the evidence of shared ancestry and low but significant differentiation. The southeast Guantanamo Bay region and the western portion of the archipelago were more differentiated than southwest and northwest manatees. The genetic distinctiveness observed in the southeast supports its recognition as a demographically independent unit for natural resource management regardless of whether it is due to historical isolation or isolation by distance. Estimates of the regional effective population sizes, with the microsatellite and SNP datasets, were small (all N_e < 60). Subsequent analyses using additional samples could better examine how the observed structure is masking simple isolation by distance patterns or whether ecological or biogeographic forces shape genetic patterns.

Purging of deleterious burden in the endangered Iberian lynx

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.

New developments in the field of genomic technologies and their relevance to conservation management

Recent technological advances in the field of genomics offer conservation managers and practitioners new tools to explore for conservation applications. Many of these tools are well developed and used by other life science fields, while others are still in development. Considering these technological possibilities, choosing the right tool(s) from the toolbox is crucial and can pose a challenging task. With this in mind, we strive to inspire, inform and illuminate managers and practitioners on how conservation efforts can benefit from the current genomic and biotechnological revolution. With inspirational case studies we show how new technologies can help resolve some of the main conservation challenges, while also informing how implementable the different technologies are. We here focus specifically on small population management, highlight the potential for genetic rescue, and discuss the opportunities in the field of gene editing to help with adaptation to changing environments. In addition, we delineate potential applications of gene drives for controlling invasive species. We illuminate that the genomic toolbox offers added benefit to conservation efforts, but also comes with limitations for the use of these novel emerging techniques.

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

Bottleneck-associated changes in the genomic landscape of genetic diversity in wild lynx populations

Demographic bottlenecks generally reduce genetic diversity through more intense genetic drift, but their net effect may vary along the genome due to the random nature of genetic drift and to local effects of recombination, mutation, and selection. Here, we analyzed the changes in genetic diversity following a bottleneck by comparing whole-genome diversity patterns in populations with and without severe recent documented declines of Iberian (Lynx pardinus, n = 31) and Eurasian lynx (Lynx lynx, n = 29). As expected, overall genomic diversity correlated negatively with bottleneck intensity and/or duration. Correlations of genetic diversity with divergence, chromosome size, gene or functional site content, GC content, or recombination were observed in nonbottlenecked populations, but were weaker in bottlenecked populations. Also, functional features under intense purifying selection and the X chromosome showed an increase in the observed density of variants, even resulting in higher θ W diversity than in nonbottlenecked populations. Increased diversity seems to be related to both a higher mutational input in those regions creating a large collection of low-frequency variants, a few of which increase in frequency during the bottleneck to the point they become detectable with our limited sample, and the reduced efficacy of purifying selection, which affects not only protein structure and function but also the regulation of gene expression. The results of this study alert to the possible reduction of fitness and adaptive potential associated with the genomic erosion in regulatory elements. Further, the detection of a gain of diversity in ultra-conserved elements can be used as a sensitive and easy-to-apply signature of genetic erosion in wild populations.

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDI) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

The boon and bane of boldness: movement syndrome as saviour and sink for population genetic diversity

BACKGROUND

Many felid species are of high conservation concern, and with increasing human disturbance the situation is worsening. Small isolated populations are at risk of genetic impoverishment decreasing within-species biodiversity. Movement is known to be a key behavioural trait that shapes both demographic and genetic dynamics and affects population survival. However, we have limited knowledge on how different manifestations of movement behaviour translate to population processes. In this study, we aimed to 1) understand the potential effects of movement behaviour on the genetic diversity of small felid populations in heterogeneous landscapes, while 2) presenting a simulation tool that can help inform conservation practitioners following, or considering, population management actions targeting the risk of genetic impoverishment.

METHODS

We developed a spatially explicit individual-based population model including neutral genetic markers for felids and applied this to the example of Eurasian lynx. Using a neutral landscape approach, we simulated reintroductions into a three-patch system, comprising two breeding patches separated by a larger patch of differing landscape heterogeneity, and tested for the effects of various behavioural movement syndromes and founder population sizes. We explored a range of movement syndromes by simulating populations with various movement model parametrisations that range from 'shy' to 'bold' movement behaviour.

RESULTS

We find that movement syndromes can lead to a higher loss of genetic diversity and an increase in between population genetic structure for both "bold" and "shy" movement behaviours, depending on landscape conditions, with larger decreases in genetic diversity and larger increases in genetic differentiation associated with bold movement syndromes, where the first colonisers quickly reproduce and subsequently dominate the gene pool. In addition, we underline the fact that a larger founder population can offset the genetic losses associated with subpopulation isolation and gene pool dominance.

CONCLUSIONS

We identified a movement syndrome trade-off for population genetic variation, whereby bold-explorers could be saviours - by connecting populations and promoting panmixia, or sinks - by increasing genetic losses via a 'founder takes all' effect, whereas shy-stayers maintain a more gradual genetic drift due to their more cautious behaviour. Simulations should incorporate movement behaviour to provide better projections of long-term population viability and within-species biodiversity, which includes genetic diversity. Simulations incorporating demographics and genetics have great potential for informing conservation management actions, such as population reintroductions or reinforcements. Here, we present such a simulation tool for solitary felids.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

Planning the peninsula-wide recovery of the Iberian lynx: identification of favourable habitat areas

The creation of new populations through reintroductions in their former range is still necessary for the long-term conservation of the Iberian lynx (Lynx pardinus). A large-scale assessment of habitat suitability is a necessary prerequisite for evaluating the possibility of carrying out lynx reintroductions. We modelled habitat suitability for the Iberian lynx in the Iberian Peninsula using MaxEnt. Lynx presence data for 2010–2013 from Andujar-Cardeña and Doñana populations were used. Habitat variables were quantified using the CORINE Land Cover. This habitat suitability model predicted an environmentally suitable area of 87,747 km2, which represents 14.08% of the whole Iberian Peninsula. Of the total suitable habitat identified, 45% is located in a large continuous area in the south-western quadrant coinciding with the historical range of the species in the mid-20th century. However, there are also relevant patches mainly in the north-eastern quadrant. About 55% of the potential area suitable for the lynx falls outside protected areas. If Iberian lynx populations are to attain levels that will facilitate their long-term survival, it will be necessary to consider potential habitat for the species throughout the whole of the Iberian Peninsula.

Anthropogenic Effects on Natural Mammalian Populations: Correlation Between Telomere Length and Coal Exposure

The Candiota coal mine in Rio Grande do Sul (RS) is one of the largest in Brazil. Coal is a fossil fuel that causes environmental impacts from its extraction to combustion due to the release of different agents, such as polycyclic aromatic hydrocarbons (PAH) and heavy metals. Ctenomys torquatus are herbivorous and subterranean rodents that dig tunnels with their paws and teeth and can be exposed to coal through contaminated food. Exposure to pollutants can cause DNA damage and affect different tissues, inducing alterations in the population structure and genetic diversity. Our study aimed to evaluate the effect of exposure to coal and its derivatives on the C. torquatus population and to examine the relationship of coal exposure with variations in absolute telomere length (aTL), global DNA methylation and genotoxicity. Our study showed an inverse correlation between telomere length and coal exposure in addition to an increase in DNA damage. The results indicate that coal and its byproducts can contribute to the alteration of the C. torquatus population structure, as evidenced by a reduction in the number of adults.

Genetic evaluation of the Iberian lynx ex situ conservation programme

Ex situ programmes have become critical for improving the conservation of many threatened species, as they establish backup populations and provide individuals for reintroduction and reinforcement of wild populations. The Iberian lynx was considered the most threatened felid species in the world in the wake of a dramatic decline during the second half of the 20th century that reduced its numbers to around only 100 individuals. An ex situ conservation programme was established in 2003 with individuals from the two well-differentiated, remnant populations, with great success from a demographic point of view. Here, we evaluate the genetic status of the Iberian lynx captive population based on molecular data from 36 microsatellites, including patterns of relatedness and representativeness of the two remnant genetic backgrounds among founders, the evolution of diversity and inbreeding over the years, and genetic differentiation among breeding facilities. In general terms, the ex situ population harbours most of the genetic variability found in the two wild populations and has been able to maintain reasonably low levels of inbreeding and high diversity, thus validating the applied management measures and potentially representing a model for other species in need of conservation.

Genomic effects of population collapse in a critically endangered ironwood tree Ostrya rehderiana

Increased human activity and climate change are driving numerous tree species to endangered status, and in the worst cases extinction. Here we examine the genomic signatures of the critically endangered ironwood tree Ostrya rehderiana and its widespread congener O. chinensis. Both species have similar demographic histories prior to the Last Glacial Maximum (LGM); however, the effective population size of O. rehderiana continued to decrease through the last 10,000 years, whereas O. chinensis recovered to Pre-LGM numbers. O. rehderiana accumulated more deleterious mutations, but purged more severely deleterious recessive variations than in O. chinensis. This purging and the gradually reduced inbreeding depression together may have mitigated extinction and contributed to the possible future survival of the outcrossing O. rehderiana. Our findings provide critical insights into the evolutionary history of population collapse and the potential for future recovery of the endangered trees.

Spatiotemporal Dynamics of Genetic Variation in the Iberian Lynx along Its Path to Extinction Reconstructed with Ancient DNA

There is the tendency to assume that endangered species have been both genetically and demographically healthier in the past, so that any genetic erosion observed today was caused by their recent decline. The Iberian lynx (Lynx pardinus) suffered a dramatic and continuous decline during the 20th century, and now shows extremely low genome- and species-wide genetic diversity among other signs of genomic erosion. We analyze ancient (N = 10), historical (N = 245), and contemporary (N = 172) samples with microsatellite and mitogenome data to reconstruct the species' demography and investigate patterns of genetic variation across space and time. Iberian lynx populations transitioned from low but significantly higher genetic diversity than today and shallow geographical differentiation millennia ago, through a structured metapopulation with varying levels of diversity during the last centuries, to two extremely genetically depauperate and differentiated remnant populations by 2002. The historical subpopulations show varying extents of genetic drift in relation to their recent size and time in isolation, but these do not predict whether the populations persisted or went finally extinct. In conclusion, current genetic patterns were mainly shaped by genetic drift, supporting the current admixture of the two genetic pools and calling for a comprehensive genetic management of the ongoing conservation program. This study illustrates how a retrospective analysis of demographic and genetic patterns of endangered species can shed light onto their evolutionary history and this, in turn, can inform conservation actions.

The return of the Iberian lynx to Portugal: local voices

BACKGROUND

Ethnographic research can help to establish dialog between conservationists and local people in reintroduction areas. Considering that predator reintroductions may cause local resistance, we assessed attitudes of different key actor profiles to the return of the Iberian lynx (Lynx pardinus) to Portugal before reintroduction started in 2015. We aimed to characterize a social context from an ethnoecological perspective, including factors such as local knowledge, perceptions, emotions, and opinions.

METHODS

We conducted semi-structured interviews (n = 131) in three different protected areas and observed practices and public meetings in order to describe reintroduction contestation, emotional involvement with the species, and local perceptions about conservation. Detailed content data analysis was undertaken and an open-ended codification of citations was performed with the support of ATLAS.ti. Besides the qualitative analyses, we further explored statistic associations between knowledge and opinions and compared different geographical areas and hunters with non-hunters among key actors.

RESULTS

Local ecological knowledge encompassed the lynx but was not shared by the whole community. Both similarities and differences between local and scientific knowledge about the lynx were found. The discrepancies with scientific findings were not necessarily a predictor of negative attitudes towards reintroduction. Contestation issues around reintroduction differ between geographical areas but did not hinder an emotional attachment to the species and its identification as a territory emblem. Among local voices, financial compensation was significantly associated to hunters and nature tourism was cited the most frequent advantage of lynx presence. Materialistic discourses existed in parallel with non-economic factors and the existence of moral agreement with its protection. The considerable criticism and reference to restrictions by local actors concerning protected areas and conservation projects indicated the experience of an imposed model of nature conservation. Opinions about participation in the reintroduction process highlighted the need for a closer dialog between all actors and administration.

CONCLUSIONS

Local voices analyzed through an ethnoecological perspective provide several views on reintroduction and nature conservation. They follow two main global trends of environmental discourse: (1) nature becomes a commodified object to exploit while contestation about wildlife is centered on financial return and (2) emblematic wild species create an emotional attachment, become symbolic, and gather moral agreement for nature protection. Lynx reintroduction has been not only just a nature protection theme but also a negotiation process with administration. Western rural communities are not the "noble savages" and nature protectors as are other traditional groups, and actors tend to claim for benefits in a situation of reintroduction. Both parties comprehend a similar version of appropriated nature. Understanding complexity and diverse interests in local communities are useful in not oversimplifying local positions towards predator conservation. We recommend that professional conservation teams rethink their image among local populations and increase proximity with different types of key actors.

Population-level assessment of genetic diversity and habitat fragmentation in critically endangered Grauer's gorillas

OBJECTIVES

The critically endangered Grauer's gorilla (Gorilla beringei graueri) has experienced an estimated 77% population decline within a single generation. Although crucial for informed conservation decisions, there is no clear understanding about population structure and distribution of genetic diversity across the species' highly fragmented range. We fill this gap by studying several core and peripheral Grauer's gorilla populations throughout their distribution range.

MATERIALS AND METHODS

We generated genetic profiles for a sampling of an unstudied population of Grauer's gorillas from within the species' core range at 13 autosomal microsatellite loci and combined them with previously published and newly generated data from four other Grauer's gorilla populations, two mountain gorilla populations, and one western lowland gorilla population.

RESULTS

In agreement with previous studies, the genetic diversity of Grauer's gorillas is intermediate, falling between western lowland and mountain gorillas. Among Grauer's gorilla populations, we observe lower genetic diversity and high differentiation in peripheral compared with central populations, indicating a strong effect of genetic drift and limited gene flow among small, isolated forest fragments.

DISCUSSION

Although genetically less diverse, peripheral populations are frequently essential for the long-term persistence of a species and migration between peripheral and core populations may significantly enrich the overall species genetic diversity. Thus, in addition to central Grauer's gorilla populations from the core of the distribution range that clearly deserve conservation attention, we argue that conservation strategies aiming to ensure long-term species viability should include preserving peripheral populations and enhancing habitat connectivity.

Novel efficient genome-wide SNP panels for the conservation of the highly endangered Iberian lynx

BACKGROUND

The Iberian lynx (Lynx pardinus) has been acknowledged as the most endangered felid species in the world. An intense contraction and fragmentation during the twentieth century left less than 100 individuals split in two isolated and genetically eroded populations by 2002. Genetic monitoring and management so far have been based on 36 STRs, but their limited variability and the more complex situation of current populations demand more efficient molecular markers. The recent characterization of the Iberian lynx genome identified more than 1.6 million SNPs, of which 1536 were selected and genotyped in an extended Iberian lynx sample.

METHODS

We validated 1492 SNPs and analysed their heterozygosity, Hardy-Weinberg equilibrium, and linkage disequilibrium. We then selected a panel of 343 minimally linked autosomal SNPs from which we extracted subsets optimized for four different typical tasks in conservation applications: individual identification, parentage assignment, relatedness estimation, and admixture classification, and compared their power to currently used STR panels.

RESULTS

We ascribed 21 SNPs to chromosome X based on their segregation patterns, and identified one additional marker that showed significant differentiation between sexes. For all applications considered, panels of autosomal SNPs showed higher power than the currently used STR set with only a very modest increase in the number of markers.

CONCLUSIONS

These novel panels of highly informative genome-wide SNPs provide more powerful, efficient, and flexible tools for the genetic management and non-invasive monitoring of Iberian lynx populations. This example highlights an important outcome of whole-genome studies in genetically threatened species.

Retention of functional variation despite extreme genomic erosion: MHC allelic repertoires in the Lynx genus

BACKGROUND

Demographic bottlenecks erode genetic diversity and may increase endangered species' extinction risk via decreased fitness and adaptive potential. The genetic status of species is generally assessed using neutral markers, whose dynamic can differ from that of functional variation due to selection. The MHC is a multigene family described as the most important genetic component of the mammalian immune system, with broad implications in ecology and evolution. The genus Lynx includes four species differing immensely in demographic history and population size, which provides a suitable model to study the genetic consequences of demographic declines: the Iberian lynx being an extremely bottlenecked species and the three remaining ones representing common and widely distributed species. We compared variation in the most variable exon of the MHCI and MHCII-DRB loci among the four species of the Lynx genus.

RESULTS

The Iberian lynx was characterised by lower number of MHC alleles than its sister species (the Eurasian lynx). However, it maintained most of the functional genetic variation at MHC loci present in the remaining and genetically healthier lynx species at all nucleotide, amino acid, and supertype levels.

CONCLUSIONS

Species-wide functional genetic diversity can be maintained even in the face of severe population bottlenecks, which caused devastating whole genome genetic erosion. This could be the consequence of divergent alleles being retained across paralogous loci, an outcome that, in the face of frequent gene conversion, may have been favoured by balancing selection.

Territoriality ensures paternity in a solitary carnivore mammal

In solitary carnivorous mammals, territoriality is assumed to benefit male fitness by ensuring the exclusivity of matings within territories via mate guarding and female defence. However, this hypothesis remains empirically untested. Here, we examined this hypothesis for solitary territorial carnivores using the Iberian lynx (Lynx pardinus) as a case study. We expected that territorial males sire all litters born within their territories, translating into the absence of multi-paternity cases within the same litter. We analysed parentage in 43 kittens, belonging to 20 different litters. For 42 kittens, a father could be assigned using microsatellites and always coincided with the individual holding the territory. For 16 kittens from 10 litters for which we also had information on SNPs, paternity assignments coincided with microsatellites, except for a litter (two kittens) from the same litter for which a different male was assigned, but the territorial male could not be excluded. Our results indicated that multi-paternity in the Iberian lynx must be a rare event, and that territorial males sire all litters born from the females with which they share territories. We propose that both the low number of mature individuals in the lynx population and the fact that female oestrus is induced by male presence may explain results.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

PCR Strategies for Complete Allele Calling in Multigene Families Using High-Throughput Sequencing Approaches

The characterization of multigene families with high copy number variation is often approached through PCR amplification with highly degenerate primers to account for all expected variants flanking the region of interest. Such an approach often introduces PCR biases that result in an unbalanced representation of targets in high-throughput sequencing libraries that eventually results in incomplete detection of the targeted alleles. Here we confirm this result and propose two different amplification strategies to alleviate this problem. The first strategy (called pooled-PCRs) targets different subsets of alleles in multiple independent PCRs using different moderately degenerate primer pairs, whereas the second approach (called pooled-primers) uses a custom-made pool of non-degenerate primers in a single PCR. We compare their performance to the common use of a single PCR with highly degenerate primers using the MHC class I of the Iberian lynx as a model. We found both novel approaches to work similarly well and better than the conventional approach. They significantly scored more alleles per individual (11.33 ± 1.38 and 11.72 ± 0.89 vs 7.94 ± 1.95), yielded more complete allelic profiles (96.28 ± 8.46 and 99.50 ± 2.12 vs 63.76 ± 15.43), and revealed more alleles at a population level (13 vs 12). Finally, we could link each allele's amplification efficiency with the primer-mismatches in its flanking sequences and show that ultra-deep coverage offered by high-throughput technologies does not fully compensate for such biases, especially as real alleles may reach lower coverage than artefacts. Adopting either of the proposed amplification methods provides the opportunity to attain more complete allelic profiles at lower coverages, improving confidence over the downstream analyses and subsequent applications.

Apoptosis-Related Factors in the Luteal Phase of the Domestic Cat and Their Involvement in the Persistence of Corpora Lutea in Lynx

The corpus luteum (CL) is a transient gland formed in the ovary after ovulation and is the major source of progesterone. In the Iberian and Eurasian lynx, CL physiologically persist after parturition and retain their capacity to produce progesterone, thus suppressing the ovarian activity. This unique reproductive characteristic has a big impact on the success of assisted reproduction techniques in the endangered Iberian lynx. The mechanisms behind CL persistence are not yet understood and require extensive studies on potential luteotropic and luteolytic factors in felids. Because the apoptosis system has been shown to be involved in structural regression of CL in many species, we aimed to investigate the capacity of perCL to undergo apoptosis. In addition, we performed initial studies on the apoptosis system in the luteal phase of the domestic cat. No previous research on this system has been made in this species. Our factors of interest included agents of the intrinsic apoptosis pathway, i.e., pro-survival B-cell CLL/lymphoma 2 (BCL2) and pro-apoptotic BCL2-associated X protein (BAX), the executioner caspase-3 (CASP3), as well as of the extrinsic pathway, i.e., pro-apoptotic receptor FAS, and tumor necrosis factor (TNF) and its receptors (pro-apoptotic TNFRSF1A and pro-survival TNFRSF1B). We analyzed the relative mRNA levels of these factors, as well as protein localization of CASP3 and TNF during stages of pregnancy and the non-pregnant luteal phase in CL of domestic cats. The same factors were investigated in freshly ovulated CL (frCL) and perCL of Iberian and Eurasian lynx, which were histologically analyzed. All factors were present in the CL tissue of both domestic cat and lynx throughout all analyzed stages. The presence of pro-apoptotic factors BAX, CASP3, FAS and TNFRSF1A in perCL of the Eurasian and Iberian lynx might indicate the potential sensitivity of perCL to apoptotic signals. The expression of pro-survival factors BCL2 and TNFRSF1B was significantly higher in perCL compared to frCL of studied Iberian lynx, suggesting the potential involvement of these factors in the structural integrity of perCL. In both Iberian lynx and pregnant and non-pregnant domestic cats, the expression of TNFRSF1A was significantly higher in forming CL compared to other stages, suggesting the conserved involvement of this factor in the tissue reorganization during formation of the feline CL. The mRNA levels of CASP3 and TNFRSF1B were highest during regression stages of domestic cat CL. The current study provides initial results on the possible involvement of the apoptosis system in the structure and function of the feline CL and in its physiological persistence.

Population and genetic outcomes 20 years after reintroducing bobcats (Lynx rufus) to Cumberland Island, Georgia USA

In 1988-1989, 32 bobcats Lynx rufus were reintroduced to Cumberland Island (CUIS), Georgia, USA, from which they had previously been extirpated. They were monitored intensively for 3 years immediately post-reintroduction, but no estimation of the size or genetic diversity of the population had been conducted in over 20 years since reintroduction. We returned to CUIS in 2012 to estimate abundance and effective population size of the present-day population, as well as to quantify genetic diversity and inbreeding. We amplified 12 nuclear microsatellite loci from DNA isolated from scats to establish genetic profiles to identify individuals. We used spatially explicit capture-recapture population estimation to estimate abundance. From nine unique genetic profiles, we estimate a population size of 14.4 (SE = 3.052) bobcats, with an effective population size (N e) of 5-8 breeding individuals. This is consistent with predictions of a population viability analysis conducted at the time of reintroduction, which estimated the population would average 12-13 bobcats after 10 years. We identified several pairs of related bobcats (parent-offspring and full siblings), but ~75% of the pairwise comparisons were typical of unrelated individuals, and only one individual appeared inbred. Despite the small population size and other indications that it has likely experienced a genetic bottleneck, levels of genetic diversity in the CUIS bobcat population remain high compared to other mammalian carnivores. The reintroduction of bobcats to CUIS provides an opportunity to study changes in genetic diversity in an insular population without risk to this common species. Opportunities for natural immigration to the island are limited; therefore, continued monitoring and supplemental bobcat reintroductions could be used to evaluate the effect of different management strategies to maintain genetic diversity and population viability. The successful reintroduction and maintenance of a bobcat population on CUIS illustrates the suitability of translocation as a management tool for re-establishing felid populations.

Genetic Signatures of Demographic Changes in an Avian Top Predator during the Last Century: Bottlenecks and Expansions of the Eurasian Eagle Owl in the Iberian Peninsula

The study of the demographic history of species can help to understand the negative impact of recent population declines in organisms of conservation concern. Here, we use neutral molecular markers to explore the genetic consequences of the recent population decline and posterior recovery of the Eurasian eagle owl (Bubo bubo) in the Iberian Peninsula. During the last century, the species was the object of extermination programs, suffering direct persecution by hunters until the 70's. Moreover, during the last decades the eagle owl was severely impacted by increased mortality due to electrocution and the decline of its main prey species, the European rabbit (Oryctolagus cuniculus). In recent times, the decrease of direct persecution and the implementation of some conservation schemes have allowed the species' demographic recovery. Yet, it remains unknown to which extent the past population decline and the later expansion have influenced the current species' pattern of genetic diversity. We used eight microsatellite markers to genotype 235 eagle owls from ten Spanish subpopulations and analyse the presence of genetic signatures attributable to the recent population fluctuations experienced by the species. We found moderate levels of differentiation among the studied subpopulations and Bayesian analyses revealed the existence of three genetic clusters that grouped subpopulations from central, south-western and south-eastern Spain. The observed genetic structure could have resulted from recent human-induced population fragmentation, a patchy distribution of prey populations and/or the philopatric behaviour and habitat selection of the species. We detected an old population bottleneck, which occurred approximately 10,000 years ago, and significant signatures of recent demographic expansions. However, we did not find genetic signatures for a recent bottleneck, which may indicate that population declines were not severe enough to leave detectable signals on the species genetic makeup or that such signals have been eroded by the rapid demographic recovery experienced by the species in recent years.

Individual identification of endangered species using mosquito blood meals: a proof-of-concept study in Iberian lynx

Host identification from mosquito blood meals has been routinely used to identify the feeding preferences of insects in studies on transmission of vector-borne pathogens. Here, we identified for the first time the susceptibility of the endangered Iberian lynx (Lynx pardinus) to the attack of a wild mosquito female, the mosquito Anopheles atroparvus. Furthermore, we used 11 microsatellite markers to test for the utility of vertebrate DNA isolated from insect blood meals for individual identification of wildlife. Only the three smallest markers were successfully amplified; however, this genotype did not match with any of the previously genotyped individuals in southern Spain. These results support the use of DNA from mosquito blood meals as a non-invasive source of DNA and a powerful tool on epidemiological and conservation biology studies. However, as may be the case of other non-invasive sampling methods, the utility of this technique is probably limited by the quantity and quality of vertebrate DNA.

Causes and consequences of fine-scale population structure in a critically endangered freshwater seal

Background

Small, genetically uniform populations may face an elevated risk of extinction due to reduced environmental adaptability and individual fitness. Fragmentation can intensify these genetic adversities and, therefore, dispersal and gene flow among subpopulations within an isolated population is often essential for maintaining its viability. Using microsatellite and mtDNA data, we examined genetic diversity, spatial differentiation, interregional gene flow, and effective population sizes in the critically endangered Saimaa ringed seal (Phoca hispida saimensis), which is endemic to the large but highly fragmented Lake Saimaa in southeastern Finland.

Results

Microsatellite diversity within the subspecies (H_E = 0.36) ranks among the lowest thus far recorded within the order Pinnipedia, with signs of ongoing loss of individual heterozygosity, reflecting very low effective subpopulation sizes. Bayesian assignment analyses of the microsatellite data revealed clear genetic differentiation among the main breeding areas, but interregional structuring was substantially weaker in biparentally inherited microsatellites (F_ST = 0.107) than in maternally inherited mtDNA (F_ST = 0.444), indicating a sevenfold difference in the gene flow mediated by males versus females.

Conclusions

Genetic structuring in the population appears to arise from the joint effects of multiple factors, including small effective subpopulation sizes, a fragmented lacustrine habitat, and behavioural dispersal limitation. The fine-scale differentiation found in the landlocked Saimaa ringed seal is especially surprising when contrasted with marine ringed seals, which often exhibit near-panmixia among subpopulations separated by hundreds or even thousands of kilometres. Our results demonstrate that population structures of endangered animals cannot be predicted based on data on even closely related species or subspecies.