Search

An accurate assessment of species diversity is a cornerstone of biology and conservation. The lynx spiders (Araneae: Oxyopidae) represent one of the most diverse and widespread cursorial spider groups, however their species richness in Asia is highly underestimated. In this study, we revised species diversity with extensive taxon sampling in Taiwan and explored species boundaries based on morphological traits and genetic data using a two-step approach of molecular species delimitation. Firstly, we employed a single COI dataset and applied two genetic distance-based methods: ABGD and ASAP, and two topology-based methods: GMYC and bPTP. Secondly, we further analyzed the lineages that were not consistently delimited, and incorporated H3 to the dataset for a coalescent-based analysis using BPP. A total of eight morphological species were recognized, including five new species, Hamataliwa cordivulvasp. nov., Hamat. leporaurissp. nov., Tapponia auriolasp. nov., T. parvasp. nov. and T. rarobulbussp. nov., and three newly recorded species, Hamadruas hieroglyphica (Thorell, 1887), Hamat. foveata Tang & Li, 2012 and Peucetia latikae Tikader, 1970. All eight morphological species exhibited reciprocally monophyletic lineages. The results of molecular-based delimitation analyses suggested a variety of species hypotheses that did not fully correspond to the eight morphological species. We found that Hamat. cordivulvasp. nov. and Hamat. foveata showed shallow genetic differentiation in the COI, but they were unequivocally distinguishable according to their genitalia. In contrast, T. parvasp. nov. represented a deep divergent lineage, while differences of genitalia were not detected. This study highlights the need to comprehensively employ multiple evidence and methods to delineate species boundaries and the values of diagnostic morphological characters for taxonomic studies in lynx spiders. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

An empirical question of long-standing interest is how price promotions affect a brand’s sale shares in the fast-moving consumer-goods market. We investigated this question with concurrent promotions and sales records of specialty beer brands pooled over Tesco stores in the UK. Most brands were continuously promoted, rendering infeasible a conventional approach of establishing impact against an off-promotion sales baseline, and arguing in favor of a dynamics approach. Moreover, promotion/sales records were volatile without easily-discernable regularity. Past work conventionally attributed volatility to the impact of exogenous random shocks on stable markets, and reasoned that promotions have only an ephemeral impact on sales shares in stationary mean-reverting stochastic markets, or a persistent freely-wandering impact in nonstationary markets. We applied new empirical methods from the applied sciences to uncover an overlooked alternative: ‘systematic persistence’ in which promotional impacts evolve systematically in an endogenously-unstable market governed by deterministic-nonlinear dynamics. We reconstructed real-world market dynamics from the Tesco dataset, and detected deterministic-nonlinear market dynamics. We used reconstructed market dynamics to identify a complex network of systematic interactions between promotions and sales shares among competing brands, and quantified/characterized the dynamics of these interactions. For the majority of weeks in the study, we found that: (1) A brand’s promotions drove down own sales shares (a possibility recognized in the literature), but ‘cannibalized’ sales shares of competing brands (perhaps explaining why brands were promoted despite a negative marginal impact on own sales shares); and (2) Competitive interactions between brands owned by the same multinational brewery differed from those with outside brands. In particular, brands owned by the same brewery enjoyed a ‘mutually-beneficial’ relationship in which an incremental increase in the sales share of one marginally increased the sales share of the other. Alternatively, the sales shares of brands owned by different breweries preyed on each other’s market shares. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Coyotes (Canis latrans) colonized the eastern United States over the last century and formed a 3-species predator guild with bobcats (Lynx rufus) and gray foxes (Urocyon cinereoargenteus) across much of the southeastern United States. Diets among the three species vary along with respective impacts on game species such as white-tailed deer (Odocoileus virginianus) and wild turkeys (Meleagris gallopavo). To determine predation impacts on vertebrate prey and dietary overlap in consumption of prey items, we assessed diets of coyote, bobcat, and gray fox during spring, coinciding with white-tailed deer fawning and wild turkey nesting and brood rearing. We sampled across three sites along the Savannah River in South Carolina from mid-May through mid-June of 2020-2021. We collected 180 scat samples along 295.9 kilometers (71.1-122.4 km/site) of unpaved secondary roads and used DNA metabarcoding to determine vertebrate diet items. We identified predator species of scat using DNA metabarcoding and species-specific mtDNA fragment analysis (153 were coyote, 20 bobcat, and seven gray fox). Overall, we found evidence that two species, coyote and bobcat, consumed deer while all three consumed turkeys. Frequency of deer in the diet varied across sites for coyotes from 62-86% and wild turkey was present with a frequency of occurrence of 9% for coyotes, 5% for bobcats, and 14% for gray fox. Vertebrate diet specialization was evident across predator species with high frequency of deer in coyote diets, rabbits and small mammals in bobcat diets, and herpetofauna in gray fox diets. During deer fawning and wild turkey nesting and brood rearing, dietary overlap appears to be mediated by disparate selection of prey items, which reduced competition among coyotes, bobcats, and gray foxes. Use of DNA metabarcoding may augment our understanding of dietary preferences within this predator guild by providing increased resolution of diet composition among important game species., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Information regarding species’ status at a regional scale is instrumental for effective conservation planning. Some regions of southwestern Europe, such as Portugal, albeit included in the Mediterranean biodiversity hotspot, lack a detailed assessment of the distribution patterns of several taxonomic groups, such as carnivores. Moreover, information is scattered, often unreliable and biased towards some species or regions. This study aimed at reviewing the existing knowledge on mammalian terrestrial carnivores in Portugal, to analyse research trends, update the species checklist and assess their historical and current distribution patterns. We conducted a comprehensive review of 755 scientific studies to analyse several publication metrics and compiled 20,189 presence records of all mammalian terrestrial carnivores occurring in Portugal since historical times to evaluate their distribution patterns. Carnivore research in Portugal began in the 18th century, with a recent boost in the mid-1990s, and has been biased towards certain research topics and regionally threatened species. There are 15 extant species in Portugal, with nine occurring across the country, six showing a more limited range, as well as one additional species currently locally extinct (Ursus arctos). Over the last decades, the distribution ranges of seven species apparently remained stable, two expanded, two contracted, and three showed unclear trends. The presence of a new invasive carnivore, the raccoon (Procyon lotor), is also documented here. This study illustrates the relevance of a comprehensive analysis of non-systematic data to assess the historic and current status of mammalian terrestrial carnivores at a national level, and to identify knowledge gaps and research priorities. [ABSTRACT FROM AUTHOR]

Background: One of the classic approaches in environmental economics is the environmental Kuznets curve, which predicts that when a national economy grows from low to medium levels, threats to biodiversity conservation increase, but they decrease when the economy moves from medium to high. We evaluated this approach by examining how population densities of the brown bear (Ursus arctos), gray wolf (Canis lupus), and Eurasian lynx (Lynx lynx) were related to the national economy in 24 European countries. Methodology/Principal findings: We used forest proportions, the existence of a compensation system, and country group (former socialist countries, Nordic countries, other countries) as covariates in a linear model with the first- and the second-order polynomial terms of per capita gross domestic product (GDP). Country group was treated as a random factor, but remained insignificant and was ignored. All models concerning brown bear and wolf provided evidence that population densities decreased with increasing GDP, but densities of lynx were virtually independent of GDP. Models for the wolf explained >80% of the variation in densities, without a difference between the models with all independent variables and the model with only GDP. For the bear, the model with GDP alone accounted for 10%, and all three variables 33%, of the variation in densities. Conclusions: Wolves exhibit a higher capacity for dispersal and reproduction than bear or lynx, but still exists at the lowest densities in wealthy European countries. We are aware that several other factors, not available for our models, influenced large carnivore densities. Based on the pronounced differences among large carnivore species in their countrywide relationships between densities and GDP, and a strikingly high relationship for the gray wolf, we suggest that our results reflected differences in political history and public acceptance of these species among countries. The compensation paid for the damages caused by the carnivores is not a key to higher carnivore densities, but might be necessity for the presence of large carnivores to be accepted in countries with high GDP. [ABSTRACT FROM AUTHOR]

Landscape connectivity is integral to the persistence of metapopulations of wide ranging carnivores and other terrestrial species. The objectives of this research were to investigate the landscape characteristics essential to use of areas by lynx and bobcats in northern New England, map a habitat availability model for each species, and explore connectivity across areas of the region likely to experience future development pressure. A Mahalanobis distance analysis was conducted on location data collected between 2005 and 2010 from 16 bobcats in western Vermont and 31 lynx in northern Maine to determine which variables were most consistent across all locations for each species using three scales based on average 1) local (15 minute) movement, 2) linear distance between daily locations, and 3) female home range size. The bobcat model providing the widest separation between used locations and random study area locations suggests that they cue into landscape features such as edge, availability of cover, and development density at different scales. The lynx model with the widest separation between random and used locations contained five variables including natural habitat, cover, and elevation—all at different scales. Shrub scrub habitat—where lynx’s preferred prey is most abundant—was represented at the daily distance moved scale. Cross validation indicated that outliers had little effect on models for either species. A habitat suitability value was calculated for each 30 m2 pixel across Vermont, New Hampshire, and Maine for each species and used to map connectivity between conserved lands within selected areas across the region. Projections of future landscape change illustrated potential impacts of anthropogenic development on areas lynx and bobcat may use, and indicated where connectivity for bobcats and lynx may be lost. These projections provided a guide for conservation of landscape permeability for lynx, bobcat, and species relying on similar habitats in the region. [ABSTRACT FROM AUTHOR]

Ideal free distribution theory predicts that increased conspecific density redistributes individuals to low-density, suboptimal habitat. However, possible lags in response to population density remain poorly documented. Snowshoe hares (Lepus americanus) may exhibit density-dependent habitat selection due to its marked variation in population density. Based on 11 years (2004–2014) of snow tracking in Quebec (Canada), we investigated snowshoe hares’ short-term and delayed habitat selection responses to population density. We predicted that at high densities, hare distribution expands into low-density habitat, thus weakening the association between hares and high-density habitat. We surveyed hare tracks along 95 km of transects on average each year and georeferenced 14,240 tracks. We used Generalized Estimating Equations for track count per 100 m transect segment as a function of the proportion of different forest age classes (0–20 y, 20–40 y and 40–80 y) within 50 m of the segments. We used model coefficients for each age class as a measure of habitat preference, and modeled those coefficients as a function of a population density index in current and previous winters. Coefficients for 20- to 40-y-old forests were positive each year, indicating that this habitat was preferred. The association between track counts and 20- to 40-y-old forest significantly declined with density during the previous winter, suggesting that hare spread from preferred forest with a lagged response to density. To our knowledge, no previous empirical studies have documented a lagged habitat selection response to population density. Time lags offer possible explanation for documented deviations from ideal free distribution models. [ABSTRACT FROM AUTHOR]

In the Southern Appalachian region of the United States, harvest data has indicated the occurrence of low deer densities while exposing a trend of declining white-tailed deer (Odocoileus virginianus) populations over the past several decades in northern Georgia. A triumvirate of increasing fawn predator populations reside in the Southern Appalachian Mountains including coyotes (Canis latrans), black bears (Ursus americanus) and bobcats (Lynx rufus). This region is also characterized by a homogenous landscape composed of mature forests and sparse understory vegetation, likely lacking adequate cover to offer fawns refugia from predators. Our objectives were to estimate survival and cause-specific mortality rates of fawns while assessing a possible link between mortality risk, intrinsic fawn characteristics (i.e., birth mass, Julian birth date, sibling status), and landscape features within fawn usage areas. During 2018-2020, we radio-collared 71 fawns within the Chattahoochee National Forest of northern Georgia, USA and monitored survival to 12 weeks of age. We observed low fawn survival (cumulative = 0.157, 95% CI = 0.091-0.273; vaginal implant transmitter = 0.196, 95% CI = 0.096-0.403) with predation as the leading cause of all known mortalities (45 of 55 mortalities; 82%) due primarily to coyotes (n = 22), black bears (n = 12), and bobcats (n = 7). Relationships between landscape features and fawn predation risk were minimal with only one informative covariate. Increasing amounts of early successional land cover within fawn usage areas decreased fawn mortality risk within the first 20 days of life, but elevated mortality risk thereafter. All fawns with any amount of early successional land cover in their usage areas died of predation (n = 13) at various time intervals, suggesting limited areas of potential fawning cover may be targeted by predators. However, fawn predation risk seemed to be high regardless of landscape covariates due to the limited number of surviving fawns. Coyote-caused mortality occurred over a longer period at a consistently higher magnitude than all other forms of mortality, indicating possible delayed prey-switching behavior and coyote predation as an important factor of fawn survival. The low recruitment of fawns influenced by high predation rates and homogenous habitat conditions is likely the cause of deer population declines in the region., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Survival and predation of snowshoe hares (Lepus americanus) has been widely studied, yet there has been little quantification of the changes in vulnerability of hares to specific predators that may result from seasonal changes in vegetation and cover. We investigated survival and causes of mortalities of snowshoe hares during the late increase, peak, and decline of a population in interior Alaska. From June 2008 to May 2012, we radio-tagged 288 adult and older juvenile hares in early successional and black spruce (Picea mariana) forests and, using known-fate methods in program MARK, evaluated 85 survival models that included variables for sex, age, and body condition of hares, as well as trapping site, month, season, year, snowfall, snow depth, and air temperature. We compared the models using Akaike’s information criterion with correction for small sample size. Model results indicated that month, capture site, and body condition were the most important variables in explaining survival rates. Survival was highest in July, and more generally during summer, when alternative prey was available to predators of hares. Low survival rates coincided with molting periods, breeding activity in the spring, and the introduction of juveniles to the sample population in the fall. We identified predation as the cause of mortality in 86% of hare deaths. When the source of predation could be determined, hares were killed more often by goshawks (Accipiter gentilis) than other predators in early successional forest (30%), and more often by lynx (Lynx canadensis) than other predators in black spruce forest (31%). Great horned owls (Bubo virginianus) and coyotes (Canis latrans) represented smaller proportions of hare predation, and non-predatory causes were a minor source (3%) of mortality. Because hares rely on vegetative cover for concealment from predators, we measured cover in predation sites and habitats that the hares occupied and concluded that habitat type had a greater influence on the sources of predation than the amount of cover in any given location within a habitat. Our observations illustrate the vulnerability of hares to predators in even the densest coniferous habitat available in the boreal forest, and indicate strong seasonal changes in the rates and sources of predation. [ABSTRACT FROM AUTHOR]

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

In Central Europe, protected areas are too small to ensure survival of populations of large carnivores. In the surrounding areas, these species are often persecuted due to competition with game hunters. Therefore, understanding how predation intensity varies spatio-temporally across areas with different levels of protection is fundamental. We investigated the predation patterns of Eurasian lynx (Lynx lynx) on roe deer (Capreolus capreolus) and red deer (Cervus elaphus) in both protected areas and multi-use landscapes of the Bohemian Forest Ecosystem. Based on 359 roe and red deer killed by 10 GPS-collared lynx, we calculated the species-specific annual kill rates and tested for effects of season and lynx age, sex and reproductive status. Because roe and red deer in the study area concentrate in unprotected lowlands during winter, we modeled spatial distribution of kills separately for summer and winter and calculated-the probability of a deer killed by lynx and-the expected number of kills for areas with different levels of protection. Significantly more roe deer (46.05–74.71/year/individual lynx) were killed than red deer (1.57–9.63/year/individual lynx), more deer were killed in winter than in summer, and lynx family groups had higher annual kill rates than adult male, single adult female and subadult female lynx. In winter the probability of a deer killed and the expected number of kills were higher outside the most protected part of the study area than inside; in summer, this probability did not differ between areas, and the expected number of kills was slightly larger inside than outside the most protected part of the study area. This indicates that the intensity of lynx predation in the unprotected part of the Bohemian Forest Ecosystem increases in winter, thus mitigation of conflicts in these areas should be included as a priority in the lynx conservation strategy. [ABSTRACT FROM AUTHOR]

The effects of predation on ungulate populations depend on several factors. One of the most important factors is the proportion of predation that is additive or compensatory respectively to other mortality in the prey, i.e., the relative effect of top-down and bottom-up processes. We estimated Eurasian lynx (Lynx lynx) kill rate on roe deer (Capreolus capreolus) using radio-collared lynx. Kill rate was strongly affected by lynx social status. For males it was 4.85 ± 1.30 S.E. roe deer per 30 days, for females with kittens 6.23 ± 0.83 S.E. and for solitary females 2.71 ± 0.47 S.E. We found very weak support for effects of prey density (both for Type I (linear) and Type II (non-linear) functional responses) and of season (winter, summer) on lynx kill rate. Additionally, we analysed the growth rate in a roe deer population from 1985 to 2005 in an area, which lynx naturally re-colonized in 1996. The annual roe deer growth rate was lower after lynx re-colonized the study area, but it was also negatively influenced by roe deer density. Before lynx colonized the area roe deer growth rate was λ = 1.079 (± 0.061 S.E.), while after lynx re-colonization it was λ = 0.94 (± 0.051 S.E.). Thus, the growth rate in the roe deer population decreased by Δλ = 0.14 (± 0.080 S.E.) after lynx re-colonized the study area, which corresponded to the estimated lynx predation rate on roe deer (0.11 ± 0.042 S.E.), suggesting that lynx predation was mainly additive to other mortality in roe deer. To conclude, this study suggests that lynx predation together with density dependent factors both influence the roe deer population dynamics. Thus, both top-down and bottom-up processes operated at the same time in this predator-prey system. [ABSTRACT FROM AUTHOR]

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 Białowieża 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Over the last decades the phylogeography and genetic structure of a multitude of species inhabiting Europe and North America have been described. The flora and fauna of the vast landmasses of north-eastern Eurasia are still largely unexplored in this respect. The Eurasian lynx is a large felid that is relatively abundant over much of the Russian sub-continent and the adjoining countries. Analyzing 148 museum specimens collected throughout its range over the last 150 years we have described the large-scale genetic structuring in this highly mobile species. We have investigated the spatial genetic patterns using mitochondrial DNA sequences (D-loop and cytochrome b) and 11 microsatellite loci, and describe three phylogenetic clades and a clear structuring along an east-west gradient. The most likely scenario is that the contemporary Eurasian lynx populations originated in central Asia and that parts of Europe were inhabited by lynx during the Pleistocene. After the Last Glacial Maximum (LGM) range expansions lead to colonization of north-western Siberia and Scandinavia from the Caucasus and north-eastern Siberia from a refugium further east. No evidence of a Berinigan refugium could be detected in our data. We observed restricted gene flow and suggest that future studies of the Eurasian lynx explore to what extent the contemporary population structure may be explained by ecological variables. [ABSTRACT FROM AUTHOR]

Over the last decades the phylogeography and genetic structure of a multitude of species inhabiting Europe and North America have been described. The flora and fauna of the vast landmasses of north-eastern Eurasia are still largely unexplored in this respect. The Eurasian lynx is a large felid that is relatively abundant over much of the Russian sub-continent and the adjoining countries. Analyzing 148 museum specimens collected throughout its range over the last 150 years we have described the large-scale genetic structuring in this highly mobile species. We have investigated the spatial genetic patterns using mitochondrial DNA sequences (D-loop and cytochrome b) and 11 microsatellite loci, and describe three phylogenetic clades and a clear structuring along an east-west gradient. The most likely scenario is that the contemporary Eurasian lynx populations originated in central Asia and that parts of Europe were inhabited by lynx during the Pleistocene. After the Last Glacial Maximum (LGM) range expansions lead to colonization of north-western Siberia and Scandinavia from the Caucasus and north-eastern Siberia from a refugium further east. No evidence of a Berinigan refugium could be detected in our data. We observed restricted gene flow and suggest that future studies of the Eurasian lynx explore to what extent the contemporary population structure may be explained by ecological variables. [ABSTRACT FROM AUTHOR]

Felids generally follow a poly-estrous reproductive strategy. Eurasian lynx (Lynx lynx) display a different pattern of reproductive cyclicity where physiologically persistent corpora lutea (CLs) induce a mono-estrous condition which results in highly seasonal reproduction. The present study was based around a sono-morphological and endocrine study of captive Eurasian lynx, and a control-study on free-ranging lynx. We verified that CLs persist after pregnancy and pseudo-pregnancy for at least a two-year period. We could show that lynx are able to enter estrus in the following year, while CLs from the previous years persisted in structure and only temporarily reduced their function for the period of estrus onset or birth, which is unique among felids. The almost constant luteal progesterone secretion (average of 5 ng/ml serum) seems to prevent folliculogenesis outside the breeding season and has converted a poly-estrous general felid cycle into a mono-estrous cycle specific for lynx. The hormonal regulation mechanism which causes lynx to have the longest CL lifespan amongst mammals remains unclear. The described non-felid like ovarian physiology appears to be a remarkably non-plastic system. The lynx's reproductive ability to adapt to environmental and anthropogenic changes needs further investigation. [ABSTRACT FROM AUTHOR]

Understanding the factors shaping the dynamics of carnivore–livestock conflicts is vital to facilitate large carnivore conservation in multi-use landscapes. We investigated how the density of their main wild prey, roe deer Capreolus capreolus, modulates individual Eurasian lynx Lynx lynx kill rates on free-ranging domestic sheep Ovis aries across a range of sheep and roe deer densities. Lynx kill rates on free-ranging domestic sheep were collected in south-eastern Norway from 1995 to 2011 along a gradient of different livestock and wild prey densities using VHF and GPS telemetry. We used zero-inflated negative binomial (ZINB) models including lynx sex, sheep density and an index of roe deer density as explanatory variables to model observed kill rates on sheep, and ranked the models based on their AICc values. The model including the effects of lynx sex and sheep density in the zero-inflation model and the effect of lynx sex and roe deer density in the negative binomial part received most support. Irrespective of sheep density and sex, we found the lowest sheep kill rates in areas with high densities of roe deer. As roe deer density decreased, males killed sheep at higher rates, and this pattern held for both high and low sheep densities. Similarly, females killed sheep at higher rates in areas with high densities of sheep and low densities of roe deer. However, when sheep densities were low females rarely killed sheep irrespective of roe deer density. Our quantification of depredation rates can be the first step towards establishing fairer compensation systems based on more accurate and area specific estimation of losses. This study demonstrates how we can use ecological theory to predict where losses of sheep will be greatest, and can be used to identify areas where mitigation measures are most likely to be needed. [ABSTRACT FROM AUTHOR]

Understanding habitat selection of top predators is critical to predict their impacts on ecological communities and interactions with humans, particularly in recovering populations. We analyzed habitat selection in a recovering population of bobcats (Lynx rufus) in south-central Indiana using a Random Forest model. We predicted that bobcats would select forest habitat and forest edges but avoid agriculture to maximize encounters with prey species. We also predicted that bobcats would avoid developed areas and roads to minimize potential antagonistic interactions with humans. Results partially supported our predictions and were consistent with bobcats in the early stages of population expansion. Bobcats exhibited elevated use near forest edges, thresholds of avoidance near agriculture, and thresholds of selection for low and intermediate habitat heterogeneity. Bobcats exhibited peak probability of use 1-3 km from major roads, >800 m from minor roads, and <1km from developed areas, suggesting tradeoffs in reward for high-quality hunting areas and mortality risk. Our Random Forest model highlighted complex non-linear patterns and revealed that most shifts in habitat use occurred within 1 km of the edge of each habitat type. These results largely supported previous studies in the Midwest and across North America but also produced refinements of bobcat habitat use in our system, particularly at habitat boundaries. Refined models of habitat selection by carnivores enable improved prediction of the most suitable habitat for recovering populations and provides useful information for conservation., Competing Interests: The authors have declared that no competing interests exist.

Recent studies have indicated the existence of an extensive trans-genomic trans-mural co-metabolism between gut microbes and animal hosts that is diet-, host phylogeny- and provenance-influenced. Here, we analyzed the biodiversity at the level of small subunit rRNA gene sequence and the metabolic composition of 18 Mbp of consensus metagenome sequences and activity characteristics of bacterial intra-cellular extracts, in wild Iberian lynx (Lynx pardinus) fecal samples. Bacterial signatures (14.43% of all of the Firmicutes reads and 6.36% of total reads) related to the uncultured anaerobic commensals Anaeroplasma spp., which are typically found in ovine and bovine rumen, were first identified. The lynx gut was further characterized by an over-representation of 'presumptive' aquaporin aqpZ genes and genes encoding 'active' lysosomal-like digestive enzymes that are possibly needed to acquire glycerol, sugars and amino acids from glycoproteins, glyco(amino)lipids, glyco(amino)glycans and nucleoside diphosphate sugars. Lynx gut was highly enriched (28% of the total glycosidases) in genes encoding α-amylase and related enzymes, although it exhibited low rate of enzymatic activity indicative of starch degradation. The preponderance of β-xylosidase activity in protein extracts further suggests lynx gut microbes being most active for the metabolism of β-xylose containing plant N-glycans, although β-xylosidases sequences constituted only 1.5% of total glycosidases. These collective and unique bacterial, genetic and enzymatic activity signatures suggest that the wild lynx gut microbiota not only harbors gene sets underpinning sugar uptake from primary animal tissues (with the monotypic dietary profile of the wild lynx consisting of 80-100% wild rabbits) but also for the hydrolysis of prey-derived plant biomass. Although, the present investigation corresponds to a single sample and some of the statements should be considered qualitative, the data most likely suggests a tighter, more coordinated and complex evolutionary and nutritional ecology scenario of carnivore gut microbial communities than has been previously assumed. [ABSTRACT FROM AUTHOR]

The long-standing view in ecology is that disparity in overall resource selection is the basis for identifying niche breadth patterns, with species having narrow selection being classified ''specialists'' and those with broader selection being ''generalists''. The standard model of niche breadth characterizes generalists and specialists as having comparable levels of overall total resource exploitation, with specialists exploiting resources at a higher level of performance over a narrower range of conditions. This view has gone largely unchallenged. An alternate model predicts total resource use being lower for the specialized species with both peaking at a comparable level of performance over a particular resource gradient. To reconcile the niche breadth paradigm we contrasted both models by developing range-wide species distribution models for Canada lynx, Lynx canadensis, and bobcat, Lynx rufus. Using a suite of environmental factors to define each species' niche, we determined that Canada lynx demonstrated higher total performance over a restricted set of variables, specifically those related to snow and altitude, while bobcat had higher total performance across most variables. Unlike predictions generated by the standard model, bobcat level of exploitation was not compromised by the trade-off with peak performance, and Canada lynx were not restricted to exploiting a narrower range of conditions. Instead, the emergent pattern was that specialist species have a higher total resource utilization and peak performance value within a smaller number of resources or environmental axes than generalists. Our results also indicate that relative differences in niche breadth are strongly dependent on the variable under consideration, implying that the appropriate model describing niche breadth dynamics between specialists and generalists may be more complex than either the traditional heuristic or our modified version. Our results demonstrate a need to re-evaluate traditional, but largely untested, assumptions regarding resource utilization in species with broad and narrow niches. [ABSTRACT FROM AUTHOR]

Background: As a conservation tool, supplementary feeding programs may be directed to specific individuals or sectors of the target population whose productivity or survival is thought to be limited by food scarcity. However, the use of supplemental food by different sex and age classes has received little attention. We studied individual variation in the access of the endangered Iberian lynx (Lynx pardinus) to supplementary food. Methodology/Principal Findings: From 5349 pictures taken with automatic cameras placed in 25 feeding stations, we identified 28 individuals whose sex and age were known. All individuals known to live in areas subjected to supplementation regularly visited feeding stations. Food consumption was not proportional to expected variations in energy demand within sex and age classes. Food consumption by males was higher than by females, and increased with age, in agreement with a despotic distribution. Food consumption also increased with lynx body mass, and this pattern held for individuals sharing the same breeding territories. The access of inferior competitors increased with the number of feeding stations available within lynx territories. Conclusions/Significance: All lynx exposed to food supplementation made a regular use of extra food but individuals predicted to be competitively dominant visited stations more frequently than subordinates of the same breeding territory. Our results suggest that insufficient provision of supplementary food could restrict the access of juveniles, or even adult females, to feeding stations. Limited consumption by these target individuals may compromise the efficiency of the supplementary feeding programme at the population level, in endangered species that, as the Iberian lynx, exhibit marked sexual dimorphism in body size. [ABSTRACT FROM AUTHOR]

Doñana National Park (DNP) in southern Spain is a UNESCO Biosphere Reserve where commercial hunting and wildlife artificial feeding do not take place and traditional cattle husbandry still exists. Herein, we hypothesized that Mycobacterium bovis infection prevalence in wild ungulates will depend on host ecology and that variation in prevalence will reflect variation in the interaction between hosts and environmental risk factors. Cattle bTB reactor rates increased in DNP despite compulsory testing and culling of infected animals. In this study, 124 European wild boar, 95 red deer, and 97 fallow deer were sampled from April 2006 to April 2007 and analyzed for M. bovis infection. Modelling and GIS were used to identify risk factors and intra and inter-species relationships. Infection with M. bovis was confirmed in 65 (52.4%) wild boar, 26 (27.4%) red deer and 18 (18.5%) fallow deer. In the absence of cattle, wild boar M. bovis prevalence reached 92.3% in the northern third of DNP. Wild boar showed more than twice prevalence than that in deer (p<0.001). Modelling revealed that M. bovis prevalence decreased from North to South in wild boar (p<0.001) and red deer (p<0.01), whereas no spatial pattern was evidenced for fallow deer. Infection risk in wild boar was dependent on wild boar M. bovis prevalence in the buffer area containing interacting individuals (p<0.01). The prevalence recorded in this study is among the highest reported in wildlife. Remarkably, this high prevalence occurs in the absence of wildlife artificial feeding, suggesting that a feeding ban alone would have a limited effect on wildlife M. bovis prevalence. In DNP, M. bovis transmission may occur predominantly at the intra-species level due to ecological, behavioural and epidemiological factors. The results of this study allow inferring conclusions on epidemiological bTB risk factors in Mediterranean habitats that are not managed for hunting purposes. Our results support the need to consider wildlife species for the control of bTB in cattle and strongly suggest that bTB may affect animal welfare and conservation. [ABSTRACT FROM AUTHOR]

Catch-per-unit-effort (CPUE) is often used to monitor wildlife populations and to develop statistical population models. Animals caught and released are often not included in CPUE metrics and their inclusion may create more accurate indices of abundance. We used 21 years of detailed harvest records for bobcat (Lynx rufus) in Wisconsin, U.S.A., to calculate CPUE and 'actual CPUE' (ACPUE; including animals caught and released) from bobcat hunters and trappers. We calibrated these metrics to an independent estimate of bobcat abundance and attempted to create simple but effective models to estimate CPUE and ACPUE using harvest success data (i.e., bobcats harvested/available permits). CPUE showed virtually no relationship with bobcat abundance across all years, but both CPUE and ACPUE had stronger, non-linear, and negative relationships with abundance during the periods when the population was decreasing. Annual harvest success strongly predicted composite ACPUE and CPUE from hunters and trappers and hunter ACPUE and CPUE but was a poorer predictor of trapper ACPUE and CPUE. The non-linear, and sometimes weak, relationships with bobcat abundance likely reflect the increasing selectivity of bobcat hunters for trophy animals. Studies calibrating per-unit-effort metrics against abundance should account for population trajectories and different harvest methods (e.g., hunting and trapping). Our results also highlight the potential for estimating per-unit-effort metrics from relatively simple and inexpensive data sources and we encourage additional research into the use of per-unit-effort metrics for population estimation., Competing Interests: The authors declare they have no competing financial or non-financial interests.

Mammalian carnivores can be particularly sensitive to human disturbance, even within protected areas (PAs). Our objective was to understand how human disturbance affects carnivore communities in southern Arizona, USA by studying habitat occupancy based on data collected using non-invasive methods in three PAs with different levels of human disturbance. Carnivore occupancy varied based on human disturbance variables (i.e., roads, trails, etc.). Common carnivore species (coyotes, gray foxes, and bobcats) had high occupancy probability in highly disturbed sites, while all other carnivore species had a higher probability of occupancy in low disturbance protected areas. Additionally, overall carnivore diversity was higher in PAs with low human disturbance. Edges of PAs appeared to negatively impact occupancy of nearly all carnivore species. We also found the presence of roads and trails, and not necessarily how much they are used, had a significant negative impact on the occupancy of most carnivore species. Furthermore, the overall level of disturbance within a PA influenced how sensitive carnivores were to human disturbance variables. Carnivores were more sensitive in PAs with higher levels of disturbance and were relatively unaffected by disturbance variables in a PA with low base levels of disturbance. Increased visitation to PAs, expected with the region's high level of population growth, is likely to cause shifts in the carnivore communities favoring species that are less sensitive to disturbance.

The use of camera traps as a tool for studying wildlife populations is commonplace. However, few have considered how the number of detections of wildlife differ depending upon the number of camera traps placed at cameras-sites, and how this impacts estimates of occupancy and community composition. During December 2015-February 2016, we deployed four camera traps per camera-site, separated into treatment groups of one, two, and four camera traps, in southern Illinois to compare whether estimates of wildlife community metrics and occupancy probabilities differed among survey methods. The overall number of species detected per camera-site was greatest with the four-camera survey method (P<0.0184). The four-camera survey method detected 1.25 additional species per camera-site than the one-camera survey method, and was the only survey method to completely detect the ground-dwelling silvicolous community. The four-camera survey method recorded individual species at 3.57 additional camera-sites (P = 0.003) and nearly doubled the number of camera-sites where white-tailed deer (Odocoileus virginianus) were detected compared to one- and two-camera survey methods. We also compared occupancy rates estimated by survey methods; as the number of cameras deployed per camera-site increased, occupancy estimates were closer to naïve estimates, detection probabilities increased, and standard errors of detection probabilities decreased. Additionally, each survey method resulted in differing top-ranked, species-specific occupancy models when habitat covariates were included. Underestimates of occurrence and misrepresented community metrics can have significant impacts on species of conservation concern, particularly in areas where habitat manipulation is likely. Having multiple camera traps per site revealed significant shortcomings with the common one-camera trap survey method. While we realize survey design is often constrained logistically, we suggest increasing effort to at least two camera traps facing opposite directions per camera-site in habitat association studies, and to utilize camera-trap arrays when restricted by equipment availability., Competing Interests: The authors have declared that no competing interests exist.

Risk of predation is an evolutionary force that affects behaviors of virtually all animals. In this study, we examined how habitat selection by roe deer was affected by risk of predation by Eurasian lynx – the main predator of roe deer in Scandinavia. Specifically, we compared how habitat selection by roe deer varied (1) before and after lynx re-established in the study area and (2) in relation to habitat-specific risk of predation by lynx. All analyses were conducted at the spatial and temporal scales of home ranges and seasons. We did not find any evidence that roe deer avoided habitats in which the risk of predation by lynx was greatest and information-theoretic model selection showed that re-colonization by lynx had limited impact on habitat selection by roe deer despite lynx predation causing 65% of known mortalities after lynx re-colonized the area. Instead we found that habitat selection decreased when habitat availability increased for 2 of 5 habitat types (a pattern referred to as functional response in habitat selection). Limited impact of re-colonization by lynx on habitat selection by roe deer in this study differs from elk in North America altering both daily and seasonal patterns in habitat selection at the spatial scales of habitat patches and home ranges when wolves were reintroduced to Yellowstone National Park. Our study thus provides further evidence of the complexity by which animals respond to risk of predation and suggest that it may vary between ecosystems and predator-prey constellations. [ABSTRACT FROM AUTHOR]

Human-carnivore conflicts are complex and are influenced by: the spatial distribution of the conflict species; the organisation and intensity of management measures such as zoning; historical experience with wildlife; land use patterns; and local cultural traditions. We have used a geographically stratified sampling of social values and attitudes to provide a novel perspective to the human - wildlife conflict. We have focused on acceptance by and disagreements between residents (measured as Potential Conflict Index; PCI) towards illegal hunting of four species of large carnivores (bear, lynx, wolf, wolverine). The study is based on surveys of residents in every municipality in Sweden and Norway who were asked their opinion on illegal hunting. Our results show how certain social values are associated with acceptance of poaching, and how these values differ geographically independent of carnivore abundance. Our approach differs from traditional survey designs, which are often biased towards urban areas. Although these traditional designs intend to be representative of a region (i.e. a random sample from a country), they tend to receive relatively few respondents from rural areas that experience the majority of conflict with carnivores. Acceptance of poaching differed significantly between Norway (12.7-15.7% of respondents) and Sweden (3.3-4.1% of respondents). We found the highest acceptance of illegal hunting in rural areas with free-ranging sheep and strong hunting traditions. Disagreements between residents (as measured by PCI) were highest in areas with intermediate population density. There was no correlation between carnivore density and either acceptance of illegal hunting or PCI. A strong positive correlation between acceptance of illegal hunting and PCI showed that areas with high acceptance of illegal hunting are areas with high potential conflict between people. Our results show that spatially-stratified surveys are required to reveal the large scale patterns in social dynamics of human-wildlife conflicts.

Anthropogenic landscape change can lead to increased opportunities for pathogen transmission between domestic and non-domestic animals. Pumas, bobcats, and domestic cats are sympatric in many areas of North America and share many of the same pathogens, some of which are zoonotic. We analyzed bobcat, puma, and feral domestic cat samples collected from targeted geographic areas. We examined exposure to three pathogens that are taxonomically diverse (bacterial, protozoal, viral), that incorporate multiple transmission strategies (vector-borne, environmental exposure/ingestion, and direct contact), and that vary in species-specificity. Bartonella spp., Feline Immunodeficiency Virus (FIV), and Toxoplasma gondii IgG were detected in all three species with mean respective prevalence as follows: puma 16%, 41% and 75%; bobcat 31%, 22% and 43%; domestic cat 45%, 10% and 1%. Bartonella spp. were highly prevalent among domestic cats in Southern California compared to other cohort groups. Feline Immunodeficiency Virus exposure was primarily associated with species and age, and was not influenced by geographic location. Pumas were more likely to be infected with FIV than bobcats, with domestic cats having the lowest infection rate. Toxoplasma gondii seroprevalence was high in both pumas and bobcats across all sites; in contrast, few domestic cats were seropositive, despite the fact that feral, free ranging domestic cats were targeted in this study. Interestingly, a directly transmitted species-specific disease (FIV) was not associated with geographic location, while exposure to indirectly transmitted diseases--vector-borne for Bartonella spp. and ingestion of oocysts via infected prey or environmental exposure for T. gondii--varied significantly by site. Pathogens transmitted by direct contact may be more dependent upon individual behaviors and intra-specific encounters. Future studies will integrate host density, as well as landscape features, to better understand the mechanisms driving disease exposure and to predict zones of cross-species pathogen transmission among wild and domestic felids.

Background: Several felids are endangered and threatened by the illegal wildlife trade. Establishing geographic origin of tissues of endangered species is thus crucial for wildlife crime investigations and effective conservation strategies. As shown in other species, stable isotope analysis of hydrogen and oxygen in hair (δD(h), δ(18)O(h)) can be used as a tool for provenance determination. However, reliably predicting the spatial distribution of δD(h) and δ(18)O(h) requires confirmation from animal tissues of known origin and a detailed understanding of the isotopic routing of dietary nutrients into felid hair., Methodology/findings: We used coupled δD(h) and δ(18)O(h) measurements from the North American bobcat (Lynx rufus) and puma (Puma concolor) with precipitation-based assignment isoscapes to test the feasibility of isotopic geo-location of felidae. Hairs of felid and rabbit museum specimens from 75 sites across the United States and Canada were analyzed. Bobcat and puma lacked a significant correlation between H/O isotopes in hair and local waters, and also exhibited an isotopic decoupling of δ(18)O(h) and δD(h). Conversely, strong δD and δ(18)O coupling was found for key prey, eastern cottontail rabbit (Sylvilagus floridanus; hair) and white-tailed deer (Odocoileus virginianus; collagen, bone phosphate)., Conclusions/significance: Puma and bobcat hairs do not adhere to expected pattern of H and O isotopic variation predicted by precipitation isoscapes for North America. Thus, using bulk hair, felids cannot be placed on δ(18)O and δD isoscapes for use in forensic investigations. The effective application of isotopes to trace the provenance of feline carnivores is likely compromised by major controls of their diet, physiology and metabolism on hair δ(18)O and δD related to body water budgets. Controlled feeding experiments, combined with single amino acid isotope analysis of diets and hair, are needed to reveal mechanisms and physiological traits explaining why felid hair does not follow isotopic patterns demonstrated in many other taxa.