BackgroundGenomic 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.ResultsWe 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.ConclusionsThe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1090-1) contains supplementary material, which is available to authorized users.

Population viability might become compromised by the loss of genetic diversity and the accumulation of inbreeding resulting from population decline and fragmentation. The Iberian lynx (Lynx pardinus) provides a paradigmatic example of a species at the verge of extinction, and because of the well-documented and different demographic histories of the two remaining populations (Doñana and Andújar), it provides the opportunity to evaluate the performance of analytical methods commonly applied to recently declined populations. We used mitochondrial sequences and 36 microsatellite markers to evaluate the current genetic status of the species and to assess the genetic signatures of its past history. Mitochondrial diversity was extremely low with only two haplotypes, alternatively fixed in each population. Both remnant populations have low levels of genetic diversity at microsatellite markers, particularly the population from Doñana, and genetic differentiation between the two populations is high. Bayesian coalescent-based methods suggest an earlier decline starting hundreds of years ago, while heterozygosity excess and M-ratio tests did not provide conclusive and consistent evidence for recent bottlenecks. Also, a model of gene flow received overwhelming support over a model of pure drift. Results that are in conflict with the known recent demography of the species call for caution in the use of these methods, especially when no information on previous demographic history is available. Overall, our results suggest that current genetic patterns in the Iberian lynx are mainly the result of its recent decline and fragmentation and alerts on possible genetic risks for its persistence. Conservation strategies should explicitly consider this threat and incorporate an integrated genetic management of wild, captive and re-introduced populations, including genetic restoration through translocations.

Theory suggests that demographic and genetic traits deteriorate (i.e., fitness and genetic diversity decrease) when populations become small, and that such deterioration could precipitate positive feedback loops called extinction vortices. We examined whether demographic attributes and genetic traits have changed over time in one of the 2 remaining small populations of the highly endangered Iberian lynx (Lynx pardinus) in Doñana, Spain. From 1983 to 2008, we recorded nontraumatic mortality rates, litter size, offspring survival, age at territory acquisition, and sex ratio. We combined these demographic attributes with measures of inbreeding and genetic diversity at neutral loci (microsatellites) and genes subjected to selection (major histocompatibility complex). Data on demographic traits were obtained through capture and radio tracking, checking dens during breeding, track surveys, and camera trapping. For genetic analyses, we obtained blood or tissue samples from captured or necropsied individuals or from museum specimens. Over time a female-biased sex ratio developed, age of territory acquisition decreased, mean litter size decreased, and rates of nontraumatic mortality increased, but there were no significant changes in overall mortality rates, standardized individual heterozygosity declined steadily, and allelic diversity of exon 2 of class II major histocompatibility complex DRB genes remained constant (2 allelic variants present in all individuals analyzed). Changes in sex ratio and age of territory acquisition may have resulted from demographic stochasticity, whereas changes in litter size and nontraumatic mortality may be related to observed increases in inbreeding. Concomitant deterioration of both demographic attributes and genetic traits is consistent with an extinction vortex. The co-occurrence, with or without interaction, of demographic and genetic deterioration may explain the lack of success of conservation efforts with the Doñana population of Iberian lynx.

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.