Genome-environment association analyses reveal geographically restricted adaptive divergence across the range of the widespread Eurasian carnivore Lynx lynx (Linnaeus, 1758).

Local adaptations to the environment are an important aspect of the diversity of a species and their discovery, description and quantification has important implications for the fields of taxonomy, evolutionary and conservation biology. In this study, we scan genomes from several populations across the distributional range of the Eurasian lynx, with the objective of finding genomic windows under positive selection which may underlie local adaptations to different environments. A total of 394 genomic windows are found to be associated to local environmental conditions, and they are enriched for genes involved in metabolism, behaviour, synaptic organization and neural development. Adaptive genetic structure, reconstructed from SNPs in candidate windows, is considerably different than the neutral genetic structure of the species. A widespread adaptively homogeneous group is recovered occupying areas of harsher snow and temperature climatic conditions in the north-western, central and eastern parts of the distribution. Adaptively divergent populations are recovered in the westernmost part of the range, especially within the Baltic population, but also predicted for different patches in the western and southern part of the range, associated with different snow and temperature regimes. Adaptive differentiation driven by climate does not correlate much with the subspecies taxonomic delimitations, suggesting that subspecific divergences are mostly driven by neutral processes of genetic drift and gene flow. Our results will aid the selection of source populations for assisted gene flow or genetic rescue programs by identifying what climatic patterns to look for as predictors of pre-adaptation of individuals. Particularly, the Carpathian population is confirmed as the best source of individuals for the genetic rescue of the endangered, isolated and genetically eroded Balkan population. Additionally, reintroductions in central and western Europe, currently based mostly on Carpathian lynxes, could consider the Baltic population as an additional source to increase adaptive variation and likely improve adaptation to their milder climate.

Timing and synchrony of birth in Eurasian lynx across Europe.

The ecology and evolution of reproductive timing and synchrony have been a topic of great interest in evolutionary ecology for decades. Originally motivated by questions related to behavioral and reproductive adaptation to environmental conditions, the topic has acquired new relevance in the face of climate change. However, there has been relatively little research on reproductive phenology in mammalian carnivores. The Eurasian lynx (Lynx lynx) occurs across the Eurasian continent, covering three of the four main climate regions of the world. Thus, their distribution includes a large variation in climatic conditions, making it an ideal species to explore reproductive phenology. Here, we used data on multiple reproductive events from 169 lynx females across Europe. Mean birth date was May 28 (April 23 to July 1), but was ~10 days later in northern Europe than in central and southern Europe. Birth dates were relatively synchronized across Europe, but more so in the north than in the south. Timing of birth was delayed by colder May temperatures. Severe and cold weather may affect neonatal survival via hypothermia and avoiding inclement weather early in the season may select against early births, especially at northern latitudes. Overall, only about half of the kittens born survived until onset of winter but whether kittens were born relatively late or early did not affect kitten survival. Lynx are strict seasonal breeders but still show a degree of flexibility to adapt the timing of birth to surrounding environmental conditions. We argue that lynx give birth later when exposed to colder spring temperatures and have more synchronized births when the window of favorable conditions for raising kittens is shorter. This suggests that lynx are well adapted to different environmental conditions, from dry and warm climates to alpine, boreal, and arctic climates. This variation in reproductive timing will be favorable in times of climate change, as organisms with high plasticity are more likely to adjust to new environmental conditions.

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.

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.

Ungulate management in European national parks: Why a more integrated European policy is needed.

1. Primary objectives of national parks usually include both, the protection of natural processes and species conservation. When these objectives conflict, as occurs because of the cascading effects of large mammals (i.e., ungulates and large carnivores) on lower trophic levels, park managers have to decide upon the appropriate management while considering various local circumstances. 2. To analyse if ungulate management strategies are in accordance with the objectives defined for protected areas, we assessed the current status of ungulate management across European national parks using the naturalness concept and identified the variables that influence the management. 3. We collected data on ungulate management from 209 European national parks in 29 countries by means of a large-scale questionnaire survey. Ungulate management in the parks was compared by creating two naturalness scores. The first score reflects ungulate and large carnivore species compositions, and the second evaluates human intervention on ungulate populations. We then tested whether the two naturalness score categories are influenced by the management objectives, park size, years since establishment, percentage of government-owned land, and human impact on the environment (human influence index) using two generalized additive mixed models. 4. In 67.9% of the national parks, wildlife is regulated by culling (40.2%) or hunting (10.5%) or both (17.2%). Artificial feeding occurred in 81.3% of the national parks and only 28.5% of the national parks had a non-intervention zone covering at least 75% of the area. Furthermore, ungulate management differed greatly among the different countries, likely because of differences in hunting traditions and cultural and political backgrounds. Ungulate management was also influenced by park size, human impact on the landscape, and national park objectives, but after removing these variables from the full model the reduced models only showed a small change in the deviance explained. In areas with higher anthropogenic pressure, wildlife diversity tended to be lower and a higher number of domesticated species tended to be present. Human intervention (culling and artificial feeding) was lower in smaller national parks and when park objectives followed those set by the International Union for the Conservation of Nature (IUCN). 5. Our study shows that many European national parks do not fulfil the aims of protected area management as set by IUCN guidelines. In contrast to the USA and Canada, Europe currently has no common ungulate management policy within national parks. This lack of a common policy together with differences in species composition, hunting traditions, and cultural or political context has led to differences in ungulate management among European countries. To fulfil the aims and objectives of national parks and to develop ungulate management strategies further, we highlight the importance of creating a more integrated European ungulate management policy to meet the aims of national parks.

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.

Linking spatial patterns of terrestrial herbivore community structure to trophic interactions.

Large herbivores influence ecosystem functioning via their effects on vegetation at different spatial scales. It is often overlooked that the spatial distribution of large herbivores results from their responses to interacting top-down and bottom-up ecological gradients that create landscape-scale variation in the structure of the entire community. We studied the complexity of these cascading interactions using high-resolution camera trapping and remote sensing data in the best-preserved European lowland forest, Bialowieza Forest, Poland. We showed that the variation in spatial distribution of an entire community of large herbivores is explained by species-specific responses to both environmental bottom-up and biotic top-down factors in combination with human-induced (cascading) effects. We decomposed the spatial variation in herbivore community structure and identified functionally distinct landscape-scale herbivory regimes ('herbiscapes'), which are predicted to occur in a variety of ecosystems and could be an important mechanism creating spatial variation in herbivory maintaining vegetation heterogeneity.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species.

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Update of the helminth fauna in Eurasian lynx (Lynx lynx) in Poland.

The Eurasian lynx (Lynx lynx) is a strictly protected species of large carnivore in Poland. It inhabits forest complexes in north-eastern Poland and the Carpathian region in southern Poland. The status of the lynx within Poland requires special attention because its range decreased between 1980 and 2001 and has not yet recovered. One of the factors negatively affecting lynx populations is diseases, particularly parasites. The helminth fauna of the Eurasian lynx is not fully known in Poland. Previous coprological studies revealed that Polish lynx have been infected with seven species of nematodes, three species of cestodes, and one species of trematode. In this study, we present new data based on examination of opportunistically collected lynx carcasses. The aim of the study was to complement data on the helminth fauna of Eurasian lynx inhabiting Poland based on morphological and molecular analysis of parasites. Four species of cestodes-Taenia lynciscapreoli, Mesocestoides lineatus, Spirometra sp., and Taenia krabbei-were found for the first time in Eurasian lynx from Poland and three previously reported species of nematodes-Ancylostoma tubaeforme, Toxascaris leonina, and Toxocara cati-were confirmed. Larvae of Trichinella britovi were also detected in Eurasian lynx in Poland for the first time.

Estimating large carnivore populations at global scale based on spatial predictions of density and distribution - Application to the jaguar (Panthera onca).

Broad scale population estimates of declining species are desired for conservation efforts. However, for many secretive species including large carnivores, such estimates are often difficult. Based on published density estimates obtained through camera trapping, presence/absence data, and globally available predictive variables derived from satellite imagery, we modelled density and occurrence of a large carnivore, the jaguar, across the species' entire range. We then combined these models in a hierarchical framework to estimate the total population. Our models indicate that potential jaguar density is best predicted by measures of primary productivity, with the highest densities in the most productive tropical habitats and a clear declining gradient with distance from the equator. Jaguar distribution, in contrast, is determined by the combined effects of human impacts and environmental factors: probability of jaguar occurrence increased with forest cover, mean temperature, and annual precipitation and declined with increases in human foot print index and human density. Probability of occurrence was also significantly higher for protected areas than outside of them. We estimated the world's jaguar population at 173,000 (95% CI: 138,000-208,000) individuals, mostly concentrated in the Amazon Basin; elsewhere, populations tend to be small and fragmented. The high number of jaguars results from the large total area still occupied (almost 9 million km2) and low human densities (< 1 person/km2) coinciding with high primary productivity in the core area of jaguar range. Our results show the importance of protected areas for jaguar persistence. We conclude that combining modelling of density and distribution can reveal ecological patterns and processes at global scales, can provide robust estimates for use in species assessments, and can guide broad-scale conservation actions.

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.

Mesopredator behavioral response to olfactory signals of an apex predator.

Olfactory signals constitute an important mechanism in interspecific interactions, but little is known regarding their role in communication between predator species. We analyzed the behavioral responses of a mesopredator, the red fox (Vulpes vulpes), to an olfactory cue (scat) of an apex predator, the lynx (Lynx lynx) in Bialowieza Primeval Forest, Poland, using video camera traps. Red fox visited sites with scats more often than expected and the duration of their visits was longer at scat sites than at control sites (no scat added). Vigilant behavior, sniffing and scent marking (including over-marking) occurred more often at scat sites compared to control sites, where foxes mainly passed by. Vigilance was most pronounced during the first days of the recordings. Red fox behavior was also influenced by foxes previously visiting scat sites. They sniffed and scent marked (multiple over-marking) more frequently when the lynx scat had been over-marked previously by red fox. Fox visits to lynx scats may be seen as a trade-off between obtaining information on a potential food source (prey killed by lynx) and the potential risk of predation by an apex predator.

Does kinship affect spatial organization in a small and isolated population of a solitary felid: The Eurasian lynx?

Social organization in wild carnivores is mostly determined by patterns of family bonds, which may shape the degree of relatedness among individuals in the population. We studied kinship in a small and isolated population of a solitary carnivore, the Eurasian lynx (Lynx lynx) to evaluate its effect on spatial distribution of individuals. We investigated the relationship between spatial location and pair-wise kinship among 28 lynx individuals identified in 2004-2011 by telemetry, non-invasive sampling and genotyping with the use of 12 autosomal microsatellites in the Bialowieza Primeval Forest, Poland. The average relatedness of the lynx population was relatively low (Lynch and Ritland's R = 0.03). Females were significantly more related to each other than males with other males. The inferred pedigree showed that the population was dominated by only 2 familial groups. We did not find significant correlations between the relatedness and the extent of home range overlap or the straight-line distances between the home ranges' central points. These results suggest that the dynamics of kinship in this solitary felid may not differ from the random mating processes described in social carnivores. Although the chances of random mating could be limited to a few resident males and females, the presence of unrelated floaters may provide a "breeding buffer" that may prevent an increase of relatedness and likely inbreeding in the population. This system is likely to fail in preserving genetic diversity in small, highly isolated populations; therefore, restoring habitat connectivity is crucial to ensure sufficient immigration from neighboring populations.

Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range.

The activity patterns of most terrestrial animals are regarded as being primarily influenced by light, although other factors, such as sexual cycle and climatic conditions, can modify the underlying patterns. However, most activity studies have been limited to a single study area, which in turn limit the variability of light conditions and other factors. Here we considered a range of variables that might potentially influence the activity of a large carnivore, the Eurasian lynx, in a network of studies conducted with identical methodology in different areas spanning latitudes from 49°7'N in central Europe to 70°00'N in northern Scandinavia. The variables considered both light conditions, ranging from a day with a complete day-night cycle to polar night and polar day, as well as individual traits of the animals. We analysed activity data of 38 individual free-ranging lynx equipped with GPS-collars with acceleration sensors, covering more than 11,000 lynx days. Mixed linear additive models revealed that the lynx activity level was not influenced by the daily daylight duration and the activity pattern was bimodal, even during polar night and polar day. The duration of the active phase of the activity cycle varied with the widening and narrowing of the photoperiod. Activity varied significantly with moonlight. Among adults, males were more active than females, and subadult lynx were more active than adults. In polar regions, the amplitude of the lynx daily activity pattern was low, likely as a result of the polycyclic activity pattern of their main prey, reindeer. At lower latitudes, the basic lynx activity pattern peaked during twilight, corresponding to the crepuscular activity pattern of the main prey, roe deer. Our results indicated that the basic activity of lynx is independent of light conditions, but is modified by both individual traits and the activity pattern of the locally most important prey.

What cues do ungulates use to assess predation risk in dense temperate forests?

Anti-predator responses by ungulates can be based on habitat features or on the near-imminent threat of predators. In dense forest, cues that ungulates use to assess predation risk likely differ from half-open landscapes, as scent relative to sight is predicted to be more important. We studied, in the Bialowieza Primeval Forest (Poland), whether perceived predation risk in red deer (Cervus elaphus) and wild boar (Sus scrofa) is related to habitat visibility or olfactory cues of a predator. We used camera traps in two different set-ups to record undisturbed ungulate behavior and fresh wolf (Canis lupus) scats as olfactory cue. Habitat visibility at fixed locations in deciduous old growth forest affected neither vigilance levels nor visitation rate and cumulative visitation time of both ungulate species. However, red deer showed a more than two-fold increase of vigilance level from 22% of the time present on control plots to 46% on experimental plots containing one wolf scat. Higher vigilance came at the expense of time spent foraging, which decreased from 32% to 12% while exposed to the wolf scat. These behavioral changes were most pronounced during the first week of the experiment but continuous monitoring of the plots suggested that they might last for several weeks. Wild boar did not show behavioral responses indicating higher perceived predation risk. Visitation rate and cumulative visitation time were not affected by the presence of a wolf scat in both ungulate species. The current study showed that perceived predation risk in red deer and wild boar is not related to habitat visibility in a dense forest ecosystem. However, olfactory cues of wolves affected foraging behavior of their preferred prey species red deer. We showed that odor of wolves in an ecologically equivalent dose is sufficient to create fine-scale risk factors for red deer.

Long-range gene flow and the effects of climatic and ecological factors on genetic structuring in a large, solitary carnivore: the Eurasian lynx.

Due to their high mobility, large terrestrial predators are potentially capable of maintaining high connectivity, and therefore low genetic differentiation among populations. However, previous molecular studies have provided contradictory findings in relation to this. To elucidate patterns of genetic structure in large carnivores, we studied the genetic variability of the Eurasian lynx, Lynx lynx throughout north-eastern Europe using microsatellite, mitochondrial DNA control region and Y chromosome-linked markers. Using SAMOVA we found analogous patterns of genetic structure based on both mtDNA and microsatellites, which coincided with a relatively little evidence for male-biased dispersal. No polymorphism for the cytochrome b and ATP6 mtDNA genes and Y chromosome-linked markers were found. Lynx inhabiting a large area encompassing Finland, the Baltic countries and western Russia formed a single genetic unit, while some marginal populations were clearly divergent from others. The existence of a migration corridor was suggested to correspond with distribution of continuous forest cover. The lowest variability (in both markers) was found in lynx from Norway and Bialowieza Primeval Forest (BPF), which coincided with a recent demographic bottleneck (Norway) or high habitat fragmentation (BPF). The Carpathian population, being monomorphic for the control region, showed relatively high microsatellite diversity, suggesting the effect of a past bottleneck (e.g. during Last Glacial Maximum) on its present genetic composition. Genetic structuring for the mtDNA control region was best explained by latitude and snow cover depth. Microsatellite structuring correlated with the lynx's main prey, especially the proportion of red deer (Cervus elaphus) in its diet. Eurasian lynx are capable of maintaining panmictic populations across eastern Europe unless they are severely limited by habitat continuity or a reduction in numbers. Different correlations of mtDNA and microsatellite population divergence patterns with climatic and ecological factors may suggest separate selective pressures acting on males and females in this solitary carnivore.

The first record of Aelurostrongylus abstrusus (Angistrongylidae: Nematoda) in Eurasian lynx (Lynx lynx L.) from Poland based on fecal analysis.

MATERIAL: Thirty eight fecal samples of Eurasian Lynx (Lynx lynx L.) collected in Bialowieza Primeval Forest (E Poland) in years 2001-2004 were analysed. RESULTS: The presence of Aelurostrongylus abstrusus (L1) larvae was evidenced by use of decantation and flotation methods. The general prevalence of the infection recorded during the study was 21.1%, whereas mean intensity was 11,5 (1-33 larvae per sample). To our knowledge, this is the first case of Aelurostrongylus abstrusus recorded in Euroasian lynx from Poland.