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.

A Reference Genome Assembly of the Bobcat, Lynx rufus.

The bobcat (Lynx rufus) is a medium-sized carnivore well adapted to various environments and an indicator species for landscape connectivity. It is one of the 4 species within the extant Lynx genus in the family Felidae. Because of its broad geographic distribution and central role in food webs, the bobcat is important for conservation. Here we present a high-quality de novo genome assembly of a male bobcat located in Mendocino County, CA, as part of the California Conservation Genomics Project (CCGP). The assembly was generated using the standard CCGP pipeline from a combination of Omni-C and HiFi technologies. The primary assembly comprises 76 scaffolds spanning 2.4 Gb, represented by a scaffold N50 of 142 Mb, a contig N50 of 66.2 Mb, and a BUSCO completeness score of 95.90%. The bobcat genome will be an important resource for the effective management and conservation of this species and comparative genomics exploration.

Spatial and environmental influences on selection in a clock gene coding trinucleotide repeat in Canada lynx (Lynx canadensis).

Clock genes exhibit substantial control over gene expression and ultimately life-histories using external cues such as photoperiod, and are thus likely to be critical for adaptation to shifting seasonal conditions and novel environments as species redistribute their ranges under climate change. Coding trinucleotide repeats (cTNRs) are found within several clock genes, and may be interesting targets of selection due to their containment within exonic regions and elevated mutation rates. Here, we conduct inter-specific characterization of the NR1D1 cTNR between Canada lynx and bobcat, and intra-specific spatial and environmental association analyses of neutral microsatellites and our functional cTNR marker, to investigate the role of selection on this locus in Canada lynx. We report signatures of divergent selection between lynx and bobcat, with the potential for hybrid-mediated gene flow in the area of range overlap. We also provide evidence that this locus is under selection across Canada lynx in eastern Canada, with both spatial and environmental variables significantly contributing to the explained variation, after controlling for neutral population structure. These results suggest that cTNRs may play an important role in the generation of functional diversity within some mammal species, and allow for contemporary rates of adaptation in wild populations in response to environmental change. We encourage continued investment into the study of cTNR markers to better understand their broader relevance to the evolution and adaptation of mammals.

Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts.

Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.

A novel understanding of global DNA methylation in bobcat (Lynx rufus).

Epigenetic mechanisms may provide a novel prospective of bobcat (Lynx rufus) adaptation to habitat loss/fragmentation. Previous research has focused on bobcat behavior and genetics, but epigenetics has not been studied in bobcat. The aim of this study was to determine the quantity of global DNA methylation in the liver of 30 bobcats. DNA was extracted from liver samples obtained from the Vermont Fish and Wildlife Department. The percent of global DNA methylation was quantified and calculated using the MethylFlashTM Methylated DNA 5-mC Quantification Kit from Epigentek (Farmingdale, NY, USA). Age, sex, and carcass weight data were collected at sampling and analyzed with percent of global DNA methylation. Global DNA methylation was found to range from 0.46% to 2.76%. Age ranged from <1 to 12 years old and weight ranged from 3.18 to 13.61 kg. Further analysis of differential methylation may provide insight into novel means of bobcat conservation within different regions of Vermont. These results reinforce the need for genome-wide epigenetic studies in conservation biology.

Genome-scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx.

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high-throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20-33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non-model systems.

Urbanization reduces genetic connectivity in bobcats (Lynx rufus) at both intra- and interpopulation spatial scales.

Urbanization is a major factor driving habitat fragmentation and connectivity loss in wildlife. However, the impacts of urbanization on connectivity can vary among species and even populations due to differences in local landscape characteristics, and our ability to detect these relationships may depend on the spatial scale at which they are measured. Bobcats (Lynx rufus) are relatively sensitive to urbanization and the status of bobcat populations is an important indicator of connectivity in urban coastal southern California. We genotyped 271 bobcats at 13,520 SNP loci to conduct a replicated landscape resistance analysis in five genetically distinct populations. We tested urban and natural factors potentially influencing individual connectivity in each population separately, as well as study-wide. Overall, landscape genomic effects were most frequently detected at the study-wide spatial scale, with urban land cover (measured as impervious surface) having negative effects and topographic roughness having positive effects on gene flow. The negative effect of urban land cover on connectivity was also evident when populations were analyzed separately despite varying substantially in spatial area and the proportion of urban development, confirming a pervasive impact of urbanization largely independent of spatial scale. The effect of urban development was strongest in one population where stream habitat had been lost to development, suggesting that riparian corridors may help mitigate reduced connectivity in urbanizing areas. Our results demonstrate the importance of replicating landscape genetic analyses across populations and considering how landscape genetic effects may vary with spatial scale and local landscape structure.

Genetic constraints of population expansion of the Carpathian lynx at the western edge of its native distribution range in Central Europe.

Even though populations of many large carnivores are expanding throughout Europe, the Eurasian lynx population in the Western Carpathians seems unable to spread beyond the western boundaries of its current distributional range. Many factors, both extrinsic and intrinsic, can influence the potential for range expansion: landscape fragmentation, natal philopatry, low natural fecundity and high mortality, and low and sex-biased dispersal rates. In this study we used non-invasive genetic sampling to determine population size fluctuation, sub-structuring and social organisation of the peripheral lynx population at the Czech-Slovak border. Even though the population size has been relatively stable over the period studied (2010-2016), the individual inbreeding coefficients of residents at the end of the study were much higher than those of founders at the beginning of the study. While non-resident individuals (predominantly males) occurred regularly in the study population, only resident individuals with well-established home ranges participated in breeding and produced offspring. Almost half the offspring detected in the study (predominantly females) settled in or near the natal area. Subsequent incestuous mating resulted in production of inbred individuals, reduction of effective population size of the population, and sub-structuring of the population through formation of two distinct family lineages. Our study illustrates how social constraints, such as territoriality, breeding of residents and natal philopatry of females, lead to incestuous mating in small-sized populations, especially at the periphery of their distribution. This threat should be taken into account in planning of conservation and population recovery of species with similar social structure.

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.

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.

Genome-wide expression reveals multiple systemic effects associated with detection of anticoagulant poisons in bobcats (Lynx rufus).

Anticoagulant rodenticides (ARs) are indiscriminate toxicants that threaten nontarget predatory and scavenger species through secondary poisoning. Accumulating evidence suggests that AR exposure may have disruptive sublethal consequences on individuals that can affect fitness. We evaluated AR-related effects on genome-wide expression patterns in a population of bobcats in southern California. We identify differential expression of genes involved in xenobiotic metabolism, endoplasmic reticulum stress response, epithelial integrity and both adaptive and innate immune function. Further, we find that differential expression of immune-related genes may be attributable to AR-related effects on leucocyte differentiation. Collectively, our results provide an unprecedented understanding of the sublethal effects of AR exposure on a wild carnivore. These findings highlight potential detrimental effects of ARs on a wide variety of species worldwide that may consume poisoned rodents and indicate the need to investigate gene expression effects of other toxicants added to natural environments by humans.

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.

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.

A test of somatic mosaicism in the androgen receptor gene of Canada lynx (Lynx canadensis).

BACKGROUND: The androgen receptor, an X-linked gene, has been widely studied in human populations because it contains highly polymorphic trinucleotide repeat motifs that have been associated with a number of adverse human health and behavioral effects. A previous study on the androgen receptor gene in carnivores reported somatic mosaicism in the tissues of a number of species including Eurasian lynx (Lynx lynx). We investigated this claim in a closely related species, Canada lynx (Lynx canadensis). The presence of somatic mosaicism in lynx tissues could have implications for the future study of exonic trinucleotide repeats in landscape genomic studies, in which the accurate reporting of genotypes would be highly problematic. METHODS: To determine whether mosaicism occurs in Canada lynx, two lynx individuals were sampled for a variety of tissue types (lynx 1) and tissue locations (lynx 1 and 2), and 1,672 individuals of known sex were genotyped to further rule out mosaicism. RESULTS: We found no evidence of mosaicism in tissues from the two necropsied individuals, or any of our genotyped samples. CONCLUSIONS: Our results indicate that mosaicism does not manifest in Canada lynx. Therefore, the use of hide samples for further work involving trinucleotide repeat polymorphisms in Canada lynx is warranted.

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.

Bobcats (Lynx rufus) as a Model Organism to Investigate the Effects of Roads on Wide-Ranging Carnivores.

We are using bobcats (Lynx rufus) as a model organism to examine how roads affect the abundance, distribution, and genetic structure of a wide-ranging carnivore. First, we compared the distribution of bobcat-vehicle collisions to road density and then estimated collision probabilities for specific landscapes using a moving window with road-specific traffic volume. Next, we obtained incidental observations of bobcats from the public, camera-trap detections, and locations of bobcats equipped with GPS collars to examine habitat selection. These data were used to generate a cost-surface map to investigate potential barrier effects of roads. Finally, we have begun an examination of genetic structure of bobcat populations in relation to major road networks. Distribution of vehicle-killed bobcats was correlated with road density, especially state and interstate highways. Collision models suggested that some regions may function as demographic sinks. Simulated movements in the context of the cost-surface map indicated that some major roads may be barriers. These patterns were supported by the genetic structure of bobcats. The sharpest divisions among genetically distinct demes occurred along natural barriers (mountains and large lakes) and in road-dense regions. In conclusion, our study has demonstrated the utility of using bobcats as a model organism to understand the variety of threats that roads pose to a wide-ranging species. Bobcats may also be useful as one of a group of focal species while developing approaches to maintain existing connectivity or mitigate the negative effects of roads.

The subtle role of climate change on population genetic structure in Canada lynx.

Anthropogenically driven climatic change is expected to reshape global patterns of species distribution and abundance. Given recent links between genetic variation and environmental patterns, climate change may similarly impact genetic population structure, but we lack information on the spatial and mechanistic underpinnings of genetic-climate associations. Here, we show that current genetic variability of Canada lynx (Lynx canadensis) is strongly correlated with a winter climate gradient (i.e. increasing snow depth and winter precipitation from west-to-east) across the Pacific-North American (PNO) to North Atlantic Oscillation (NAO) climatic systems. This relationship was stronger than isolation by distance and not explained by landscape variables or changes in abundance. Thus, these patterns suggest that individuals restricted dispersal across the climate boundary, likely in the absence of changes in habitat quality. We propose habitat imprinting on snow conditions as one possible explanation for this unusual phenomenon. Coupling historical climate data with future projections, we also found increasingly diverging snow conditions between the two climate systems. Based on genetic simulations using projected climate data (2041-2070), we predicted that this divergence could lead to a threefold increase in genetic differentiation, potentially leading to isolated east-west populations of lynx in North America. Our results imply that subtle genetic structure can be governed by current climate and that substantive genetic differentiation and related ecological divergence may arise from changing climate patterns.

Modular real-time PCR screening assay for common European animal families.

A screening assay based on real-time PCR and melt curve analysis was developed to detect DNA from nine common European animal families/species and human. The assay consists of a 10-cycle universal pre-amplification followed by specific nested PCR and was designed to exploit the different melting temperatures (T m) of family/species-specific 12S ribosomal ribonucleic acid and cytochrome b fragments, which are amplified in duplex reactions. Case-related modular application is possible. Beyond determination of the animal family and discrimination from human DNA, evaluation of the melt curve in some cases additionally allows for species determination (e.g. cat vs. lynx). The method presents a quick, flexible and sample-saving approach to assess non-human DNA at low expenses, and it is especially useful in resolution of DNA mixtures.

Recent increase in species-wide diversity after interspecies introgression in the highly endangered Iberian lynx.

Genetic diversity is lost in small and isolated populations, affecting many globally declining species. Interspecific admixture events can increase genetic variation in the recipient species' gene pool, but empirical examples of species-wide restoration of genetic diversity by admixture are lacking. Here we present multi-fold coverage genomic data from three ancient Iberian lynx (Lynx pardinus) approximately 2,000-4,000 years old and show a continuous or recurrent process of interspecies admixture with the Eurasian lynx (Lynx lynx) that increased modern Iberian lynx genetic diversity above that occurring millennia ago despite its recent demographic decline. Our results add to the accumulating evidence for natural admixture and introgression among closely related species and show that this can result in an increase of species-wide genetic diversity in highly genetically eroded species. The strict avoidance of interspecific sources in current genetic restoration measures needs to be carefully reconsidered, particularly in cases where no conspecific source population exists.

Rehabilitation and release of orphaned Eurasian lynx (Lynx lynx) in Europe: Implications for management and conservation.

Rehabilitation of injured or immature individuals has become an increasingly used conservation and management tool. However, scientific evaluation of rehabilitations is rare, raising concern about post-release welfare as well as the cost-effectiveness of spending scarce financial resources. Over the past 20 years, events of juvenile Eurasian lynx presumably orphaned have been observed in many European lynx populations. To guide the management of orphaned lynx, we documented survival, rehabilitation and fate after the release and evaluated the potential relevance of lynx orphan rehabilitation for population management and conservation implications. Data on 320 orphaned lynx was collected from 1975 to 2022 from 13 countries and nine populations. The majority of orphaned lynx (55%) were taken to rehabilitation centres or other enclosures. A total of 66 orphans were released back to nature. The portion of rehabilitated lynx who survived at least one year after release was 0.66. Release location was the best predictor for their survival. Of the 66 released lynx, ten have reproduced at least once (8 females and 2 males). Conservation implications of rehabilitation programmes include managing genetic diversity in small, isolated populations and reintroducing species to historical habitats. The lynx is a perfect model species as most reintroduced populations in Central Europe show significantly lower observed heterozygosity than most of the autochthonous populations, indicating that reintroduction bottlenecks, isolation and post-release management have long-term consequences on the genetic composition of populations. The release of translocated orphans could be a valuable contribution to Eurasian lynx conservation in Europe. It is recommended to release orphans at the distribution edge or in the frame of reintroduction projects instead of a release in the core area of a population where it is not necessary from a demographic and genetic point of view. Rehabilitation programmes can have conservation implications that extend far beyond individual welfare benefits.