Insights into subspecies classification and conservation priorities of Central Asian lynx populations revealed by morphometric and genetic analyses.

The Eurasian lynx (Lynx lynx) exhibits geographic variability and phylogenetic intraspecific relationships. Previous morphological studies have suggested the existence of multiple lynx subspecies, but recent genetic research has questioned this classification, particularly in Central Asia. In this study, we aimed to analyse the geographic and genetic variation in Central Asian lynx populations, particularly the Turkestan lynx and Altai lynx populations, using morphometric data and mtDNA sequences to contribute to their taxonomic classification. The comparative analysis of morphometric data revealed limited clinal variability between lynx samples from the Altai and Tien Shan regions. By examining mtDNA fragments (control region and cytochrome b) obtained from Kazakhstani lynx populations, two subspecies were identified: L. l. isabellinus (represented by a unique haplotype of the South clade, H46) and L. l. wrangeli (represented by haplotypes H36, H45, and H47 of the East clade). L. l. isabellinus was recognized only in Tien Shan Mountain, while Altai lynx was likely identical to L. l. wrangeli and found in northern Kazakhstan, Altai Mountain, Saur and Tarbagatai Mountains, and Tien Shan Mountain. The morphological and mtDNA evidence presented in this study, although limited in sample size and number of genetic markers, renders the differentiation of the two subspecies challenging. Further sampling and compilation of whole-genome sequencing data are necessary to confirm whether the proposed subspecies warrant taxonomic standing.

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.

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.

Individual genotypes from environmental DNA: Fingerprinting snow tracks of three large carnivore species.

Continued advancements in environmental DNA (eDNA) research have made it possible to access intraspecific variation from eDNA samples, opening new opportunities to expand non-invasive genetic studies of wildlife populations. However, the use of eDNA samples for individual genotyping, as typically performed in non-invasive genetics, still remains elusive. We present successful individual genotyping of eDNA obtained from snow tracks of three large carnivores: brown bear (Ursus arctos), European lynx (Lynx lynx) and wolf (Canis lupus). DNA was extracted using a protocol for isolating water eDNA and genotyped using amplicon sequencing of short tandem repeats (STR), and for brown bear a sex marker, on a high-throughput sequencing platform. Individual genotypes were obtained for all species, but genotyping performance differed among samples and species. The proportion of samples genotyped to individuals was higher for brown bear (5/7) and wolf (7/10) than for lynx (4/9), and locus genotyping success was greater for brown bear (0.88). The sex marker was typed in six out of seven brown bear samples. Results for three species show that reliable individual genotyping, including sex identification, is now possible from eDNA in snow tracks, underlining its vast potential to complement the non-invasive genetic methods used for wildlife. To fully leverage the application of snow track eDNA, improved understanding of the ideal species- and site-specific sampling conditions, as well as laboratory methods promoting genotyping success, is needed. This will also inform efforts to retrieve and type nuclear DNA from other eDNA samples, thereby advancing eDNA-based individual and population-level studies.

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.

Beating the heat: ecology of desert bobcats.

BACKGROUND: Relative to temperate regions, little is known about bobcats (Lynx rufus) in the Sonoran Desert portion of their range, in part due to the difficulty of sampling an elusive carnivore in harsh desert environments. Here, we quantify habitat selection and evaluate diet of bobcats at Kofa National Wildlife Refuge, Arizona, USA, using multiple sampling techniques including GPS telemetry, camera traps, and DNA metabarcoding. RESULTS: Home ranges during the hot season were smaller than during the cool season. Camera trapping failed to yield a high enough detection rate to identify habitat occupancy trends but third-order resource selection from GPS-collar data showed a preference for higher elevations and rugged terrain at lower elevations. Diet composition consisted of a diverse range of available small prey items, including a higher frequency of avian prey than previously observed in bobcats. CONCLUSIONS: Desert bobcats in our study maintained smaller home ranges and primarily consumed smaller prey than their more northern relatives. This study illustrates the benefit of employing multiple, complementary sampling methods to understand the ecology of elusive species.

Phylogenetic analyses of Eurasian lynx (Lynx lynx Linnaeus, 1758) including new mitochondrial DNA sequences from Iran.

The Eurasian lynx (Lynx lynx) is one of the widespread felids in Eurasia; however, relatively little is known about the Asian subspecies, and especially the Iranian populations, which comprise the most southwestern part of its range. The current study aimed to assess the phylogenetic status of Iranian populations relative to other populations of Eurasia, by sequencing a 613 bp fragment of the mitochondrial control region. In total, 44 haplotypes were recorded from 83 sequences throughout Eurasia, two of which were found in Iran. The haplotype (H1) is dominant in all Iranian lynx populations and identical to specimens from SW Russia and central China. The second haplotype (H2) is unique and was recorded only from Ghazvin Province in the central Alborz Mountains. Both haplotypes occur in Ghazvin Province. The phylogenetic tree and a median-joining network identified four clades (i.e., East, West 1, West 2, and South). These results are congruent with previous studies and suggest that Eurasian lynx was restricted to the southern part of its range during the glacial maxima and expanded from there to East Asia and to Europe during several independent re-colonization events. The Caucasus region most like plays an important role as a refugium during glacial cycles.

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.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx.

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species' evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies' current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification.

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.

Genetic analysis indicates spatial-dependent patterns of sex-biased dispersal in Eurasian lynx in Finland.

Conservation and management of large carnivores requires knowledge of female and male dispersal. Such information is crucial to evaluate the population's status and thus management actions. This knowledge is challenging to obtain, often incomplete and contradictory at times. The size of the target population and the methods applied can bias the results. Also, population history and biological or environmental influences can affect dispersal on different scales within a study area. We have genotyped Eurasian lynx (180 males and 102 females, collected 2003-2017) continuously distributed in southern Finland (~23,000 km2) using 21 short tandem repeats (STR) loci and compared statistical genetic tests to infer local and sex-specific dispersal patterns within and across genetic clusters as well as geographic regions. We tested for sex-specific substructure with individual-based Bayesian assignment tests and spatial autocorrelation analyses. Differences between the sexes in genetic differentiation, relatedness, inbreeding, and diversity were analysed using population-based AMOVA, F-statistics, and assignment indices. Our results showed two different genetic clusters that were spatially structured for females but admixed for males. Similarly, spatial autocorrelation and relatedness was significantly higher in females than males. However, we found weaker sex-specific patterns for the Eurasian lynx when the data were separated in three geographical regions than when divided in the two genetic clusters. Overall, our results suggest male-biased dispersal and female philopatry for the Eurasian lynx in Southern Finland. The female genetic structuring increased from west to east within our study area. In addition, detection of male-biased dispersal was dependent on analytical methods utilized, on whether subtle underlying genetic structuring was considered or not, and the choice of population delineation. Conclusively, we suggest using multiple genetic approaches to study sex-biased dispersal in a continuously distributed species in which population delineation is difficult.

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.

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.

Non-invasive faecal sampling reveals spatial organization and improves measures of genetic diversity for the conservation assessment of territorial species: Caucasian lynx as a case species.

The Caucasian lynx, Lynx lynx dinniki, has one of the southernmost distributions in the Eurasian lynx range, covering Anatolian Turkey, the Caucasus and Iran. Little is known about the biology and the genetic status of this subspecies. To collect baseline genetic, ecological and behavioural data and benefit future conservation of L. l. dinniki, we monitored 11 lynx territories (396 km2) in northwestern Anatolia. We assessed genetic diversity of this population by non-invasively collecting 171 faecal samples and trapped and sampled 12 lynx individuals using box traps. We observed high allelic variation at 11 nuclear microsatellite markers, and found no signs of inbreeding despite the potential isolation of this population. We obtained similar numbers of distinct genotypes from the two sampling sources. Our results indicated that first order female relatives occupy neighbouring territories (female philopatry) and that territorial male lynx were highly unrelated to each other and to female territorial lynx, suggesting long distance male dispersal. Particular male and female resident territorial lynx and their offspring (kittens and subadults) were more likely to be trapped than resident floaters or dispersing (unrelated) lynx. Conversely, we obtained more data for unrelated lynx and higher numbers of territorials using non-invasive sampling (faeces). When invasive and non-invasive samples were analysed separately, the spatial organisation of lynx (in terms of female philopatry and females and males occupying permanent ranges) affected measures of genetic diversity in such a way that estimates of genetic diversity were reduced if only invasive samples were considered. It appears that, at small spatial scales, invasive sampling using box traps may underestimate the genetic diversity in carnivores with permanent ranges and philopatry such as the Eurasian lynx. As non-invasive sampling can also provide additional data on diet and spatial organisation, we advocate the use of such samples for conservation genetic studies of vulnerable, endangered or data deficient territorial species.

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.

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.

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.