Body mass is associated with hibernation length, body temperature, and heart rate in free-ranging brown bears.

BACKGROUND: Despite centuries of research, debate remains on the scaling of metabolic rate to mass especially for intraspecific cases. The high variation of body mass within brown bears presents a unique opportunity to study the intraspecific effects of body mass on physiological variables. The amplitude of metabolic rate reduction in hibernators is dependent on body mass of the species. Small hibernators have high metabolic rates when euthermic but experience a drastic decrease in body temperature during torpor, which is necessary to reach a very low metabolic rate. Conversely, large hibernators, such as the brown bear (Ursus arctos), show a moderate decrease in temperature during hibernation, thought to be related to the bear's large size. We studied body mass, abdominal body temperature, heart rate, and accelerometer-derived activity from 63 free-ranging brown bears (1-15 years old, 15-233 kg). We tested for relationships between body mass and body temperature, heart rate, and hibernation duration. RESULTS: The smallest individuals maintained lower body temperatures during hibernation, hibernated longer, and ended hibernation later than large bears. Unlike body temperature, winter heart rates were not associated with body mass. In summer, the opposite pattern was found, with smaller individuals having higher body temperature and daytime heart rates. Body mass was associated with body temperature in the winter hypometabolic state, even in a large hibernating mammal. Smaller bears, which are known to have higher thermal conductance, reached lower body temperatures during hibernation. During summer, smaller bears had higher body temperatures and daytime heart rates, a phenomenon not previously documented within a single mammalian species. CONCLUSION: We conclude that the smallest bears hibernated more deeply and longer than large bears, likely from a combined effect of basic thermodynamics, the higher need for energy savings, and a lower cost of warming up a smaller body.

Large carnivores and naturalness affect forest recreational value.

Recreation is a crucial contribution of nature to people, relevant for forest ecosystems. Large carnivores (LCs) are important components of forests, however, their contribution to forest recreational value has not yet been evaluated. Given the current expansion of LC populations, the ongoing forest conservation debate, and the increasing use of nature for recreational purposes, this is a timely study. We used discrete choice experiments and willingness-to-travel to determine people' preferences for both forest structural characteristics and presence of four LC species in Poland (N = 1097 respondents) and Norway (N = 1005). In both countries, two-thirds of the respondents (termed 'wildness-positive') perceived LCs as contributing positively to forest recreational value and preferred to visit old forests with trees of different species and ages and presence of dead wood (i.e. natural forests). Respondents with negative preferences towards LCs preferred more intensively managed forest ('wildness-negative'); their preferences were stronger than in wildness-positive respondents and in Norway. Preferences towards wild nature were highly polarized and there were hardly neutral people. Our results showed a strong link between preferences for LC presence and forest structure, and reflected the dualism of human-nature relationships. This study highlights the need to consider the contribution of forests and LCs to human recreation services in ecosystem management policies.

The interplay between hunting rate, hunting selectivity, and reproductive strategies shapes population dynamics of a large carnivore.

Harvest, through its intensity and regulation, often results in selection on female reproductive traits. Changes in female traits can have demographic consequences, as they are fundamental in shaping population dynamics. It is thus imperative to understand and quantify the demographic consequences of changes in female reproductive traits to better understand and anticipate population trajectories under different harvest intensities and regulations. Here, using a dynamic, frequency-dependent, population model of the intensively hunted brown bear (Ursus arctos) population in Sweden, we quantify and compare population responses to changes in four reproductive traits susceptible to harvest-induced selection: litter size, weaning age, age at first reproduction, and annual probability to reproduce. We did so for different hunting quotas and under four possible hunting regulations: (i) no individuals are protected, (ii) mothers but not dependent offspring are protected, (iii) mothers and dependent offspring of the year (cubs) are protected, and (iv) entire family groups are protected (i.e., mothers and dependent offspring of any age). We found that population growth rate declines sharply with increasing hunting quotas. Increases in litter size and the probability to reproduce have the greatest potential to affect population growth rate. Population growth rate increases the most when mothers are protected. Adding protection on offspring (of any age), however, reduces the availability of bears for hunting, which feeds back to increase hunting pressure on the nonprotected categories of individuals, leading to reduced population growth. Finally, we found that changes in reproductive traits can dampen population declines at very high hunting quotas, but only when protecting mothers. Our results illustrate that changes in female reproductive traits may have context-dependent consequences for demography. Thus, to predict population consequences of harvest-induced selection in wild populations, it is critical to integrate both hunting intensity and regulation, especially if hunting selectivity targets female reproductive strategies.

Demographic history and divergence of sibling grouse species inferred from whole genome sequencing reveal past effects of climate change.

BACKGROUND: The boreal forest is one of the largest biomes on earth, supporting thousands of species. The global climate fluctuations in the Quaternary, especially the ice ages, had a significant influence on the distribution of boreal forest, as well as the divergence and evolution of species inhabiting this biome. To understand the possible effects of on-going and future climate change it would be useful to reconstruct past population size changes and relate such to climatic events in the past. We sequenced the genomes of 32 individuals from two forest inhabiting bird species, Hazel Grouse (Tetrastes bonasia) and Chinese Grouse (T. sewerzowi) and three representatives of two outgroup species from Europe and China. RESULTS: We estimated the divergence time of Chinese Grouse and Hazel Grouse to 1.76 (0.46-3.37) MYA. The demographic history of different populations in these two sibling species was reconstructed, and showed that peaks and bottlenecks of effective population size occurred at different times for the two species. The northern Qilian population of Chinese Grouse became separated from the rest of the species residing in the south approximately 250,000 years ago and have since then showed consistently lower effective population size than the southern population. The Chinese Hazel Grouse population had a higher effective population size at the peak of the Last Glacial Period (approx. 300,000 years ago) than the European population. Both species have decreased recently and now have low effective population sizes. CONCLUSIONS: Combined with the uplift history and reconstructed climate change during the Quaternary, our results support that cold-adapted grouse species diverged in response to changes in the distribution of palaeo-boreal forest and the formation of the Loess Plateau. The combined effects of climate change and an increased human pressure impose major threats to the survival and conservation of both species.

Hibernating brown bears are protected against atherogenic dyslipidemia.

To investigate mechanisms by which hibernators avoid atherogenic hyperlipidemia during hibernation, we assessed lipoprotein and cholesterol metabolisms of free-ranging Scandinavian brown bears (Ursus arctos). In winter- and summer-captured bears, we measured lipoprotein sizes and sub-classes, triglyceride-related plasma-enzyme activities, and muscle lipid composition along with plasma-levels of antioxidant capacities and inflammatory markers. Although hibernating bears increased nearly all lipid levels, a 36%-higher cholesteryl-ester transfer-protein activity allowed to stabilize lipid composition of high-density lipoproteins (HDL). Levels of inflammatory metabolites, i.e., 7-ketocholesterol and 11ß-prostaglandin F2α, declined in winter and correlated inversely with cardioprotective HDL2b-proportions and HDL-sizes that increased during hibernation. Lower muscle-cholesterol concentrations and lecithin-cholesterol acyltransferase activity in winter suggest that hibernating bears tightly controlled peripheral-cholesterol synthesis and/or release. Finally, greater plasma-antioxidant capacities prevented excessive lipid-specific oxidative damages in plasma and muscles of hibernating bears. Hence, the brown bear manages large lipid fluxes during hibernation, without developing adverse atherogenic effects that occur in humans and non-hibernators.

Cardiomyocyte Protection by Hibernating Brown Bear Serum: Toward the Identification of New Protective Molecules Against Myocardial Infarction.

Ischemic heart disease remains one of the leading causes of death worldwide. Despite intensive research on the treatment of acute myocardial infarction, no effective therapy has shown clinical success. Therefore, novel therapeutic strategies are required to protect the heart from reperfusion injury. Interestingly, despite physical inactivity during hibernation, brown bears (Ursus arctos) cope with cardiovascular physiological conditions that would be detrimental to humans. We hypothesized that bear serum might contain circulating factors that could provide protection against cell injury. In this study, we sought to determine whether addition of bear serum might improve cardiomyocyte survival following hypoxia-reoxygenation. Isolated mouse cardiomyocytes underwent 45 min of hypoxia followed by reoxygenation. At the onset of reoxygenation, cells received fetal bovine serum (FBS; positive control), summer (SBS) or winter bear serum (WBS), or adult serums of other species, as indicated. After 2 h of reoxygenation, propidium iodide staining was used to evaluate cell viability by flow cytometry. Whereas, 0.5% SBS tended to decrease reperfusion injury, 0.5% WBS significantly reduced cell death, averaging 74.04 ± 7.06% vs. 79.20 ± 6.53% in the FBS group. This cardioprotective effect was lost at 0.1%, became toxic above 5%, and was specific to the bear. Our results showed that bear serum exerts a therapeutic effect with an efficacy threshold, an optimal dose, and a toxic effect on cardiomyocyte viability after hypoxia-reoxygenation. Therefore, the bear serum may be a potential source for identifying new therapeutic molecules to fight against myocardial reperfusion injury and cell death in general.

Concurrent BMP Signaling Maintenance and TGF-ß Signaling Inhibition Is a Hallmark of Natural Resistance to Muscle Atrophy in the Hibernating Bear.

Muscle atrophy arises from a multiplicity of physio-pathological situations and has very detrimental consequences for the whole body. Although knowledge of muscle atrophy mechanisms keeps growing, there is still no proven treatment to date. This study aimed at identifying new drivers for muscle atrophy resistance. We selected an innovative approach that compares muscle transcriptome between an original model of natural resistance to muscle atrophy, the hibernating brown bear, and a classical model of induced atrophy, the unloaded mouse. Using RNA sequencing, we identified 4415 differentially expressed genes, including 1746 up- and 2369 down-regulated genes, in bear muscles between the active versus hibernating period. We focused on the Transforming Growth Factor (TGF)-ß and the Bone Morphogenetic Protein (BMP) pathways, respectively, involved in muscle mass loss and maintenance. TGF-ß- and BMP-related genes were overall down- and up-regulated in the non-atrophied muscles of the hibernating bear, respectively, and the opposite occurred for the atrophied muscles of the unloaded mouse. This was further substantiated at the protein level. Our data suggest TGF-ß/BMP balance is crucial for muscle mass maintenance during long-term physical inactivity in the hibernating bear. Thus, concurrent activation of the BMP pathway may potentiate TGF-ß inhibiting therapies already targeted to prevent muscle atrophy.

Correcting for enzyme immunoassay changes in long term monitoring studies.

Enzyme immunoassays (EIAs) are a common tool for measuring steroid hormones in wildlife due to their low cost, commercial availability, and rapid results. Testing technologies improve continuously, sometimes requiring changes in protocols or crucial assay components. Antibody replacement between EIA kits can cause differences in EIA sensitivity, which can hinder monitoring hormone concentration over time. The antibody in a common cortisol EIA kit used for long-term monitoring of stress in wildlife was replaced in 2014, causing differences in cross reactivity and standard curve concentrations. Therefore, the objective of this study was to develop a method to standardize results following changes in EIA sensitivity. We validated this method using cortisol concentrations measured in the hair of brown bears (Ursus arctos).•We used a simple linear regression to model the relationship between cortisol concentrations using kit 1 and kit 2.•We found a linear relationship between the two kits (R2 = 0.85) and used the regression equation (kit2 = (0.98 × kit1) + 1.65) to predict cortisol concentrations in re-measured samples.•Mean predicted percent error was 16% and 72% of samples had a predicted percent error <20%, suggesting that this method is well-suited for correcting changes in EIA sensitivity.

Harvest is associated with the disruption of social and fine-scale genetic structure among matrilines of a solitary large carnivore.

Harvest can disrupt wildlife populations by removing adults with naturally high survival. This can reshape sociospatial structure, genetic composition, fitness, and potentially affect evolution. Genetic tools can detect changes in local, fine-scale genetic structure (FGS) and assess the interplay between harvest-caused social and FGS in populations. We used data on 1614 brown bears, Ursus arctos, genotyped with 16 microsatellites, to investigate whether harvest intensity (mean low: 0.13 from 1990 to 2005, mean high: 0.28 from 2006 to 2011) caused changes in FGS among matrilines (8 matrilines; 109 females ≥4 years of age), sex-specific survival and putative dispersal distances, female spatial genetic autocorrelation, matriline persistence, and male mating patterns. Increased harvest decreased FGS of matrilines. Female dispersal distances decreased, and male reproductive success was redistributed more evenly. Adult males had lower survival during high harvest, suggesting that higher male turnover caused this redistribution and helped explain decreased structure among matrilines, despite shorter female dispersal distances. Adult female survival and survival probability of both mother and daughter were lower during high harvest, indicating that matriline persistence was also lower. Our findings indicate a crucial role of regulated harvest in shaping populations, decreasing differences among "groups," even for solitary-living species, and potentially altering the evolutionary trajectory of wild populations.

Disentangling direct and indirect determinants of the duration of maternal care in brown bears: Environmental context matters.

The duration of maternal care, an important life-history trait affecting population dynamics, varies greatly within species. Yet, our understanding of its predictors is limited, mostly correlative and subject to misinterpretations, due to difficulties to disentangle the role of maternal- and offspring-related characteristics. We conducted path analysis on a dataset including 217 brown bear litters captured over a 29-year period in two populations in Sweden ('North' and 'South') facing contrasting environmental conditions to identify and quantify the causes of variation in the duration of maternal care (1.5 or 2.5 years). We showed that the causal determinants of the duration of maternal care were context-dependent. Contrary to their expected central role in the determination of the duration of maternal care, yearling mass and its direct determinants (i.e. litter size and maternal mass) were only important in the North population, where environmental conditions are harsher and the cost of extended maternal care presumably higher. In the South, the duration of maternal care was not caused by yearling mass nor any maternal or litter characteristics. Extension of maternal care may thus result from factors independent from maternal and offspring condition in the South, such as an artificial hunting-induced selection for longer maternal care through the legal protection of family groups. Our results provide an important contribution to our very limited knowledge of the direct and indirect determinants of the duration of maternal care and highlight the importance of accounting for the environmental context when assessing maternal reproductive tactics.

Specific shifts in the endocannabinoid system in hibernating brown bears.

In small hibernators, global downregulation of the endocannabinoid system (ECS), which is involved in modulating neuronal signaling, feeding behavior, energy metabolism, and circannual rhythms, has been reported to possibly drive physiological adaptation to the hibernating state. In hibernating brown bears (Ursus arctos), we hypothesized that beyond an overall suppression of the ECS, seasonal shift in endocannabinoids compounds could be linked to bear's peculiar features that include hibernation without arousal episodes and capacity to react to external disturbance. We explored circulating lipids in serum and the ECS in plasma and metabolically active tissues in free-ranging subadult Scandinavian brown bears when both active and hibernating. In winter bear serum, in addition to a 2-fold increase in total fatty acid concentration, we found significant changes in relative proportions of circulating fatty acids, such as a 2-fold increase in docosahexaenoic acid C22:6 n-3 and a decrease in arachidonic acid C20:4 n-6. In adipose and muscle tissues of hibernating bears, we found significant lower concentrations of 2-arachidonoylglycerol (2-AG), a major ligand of cannabinoid receptors 1 (CB1) and 2 (CB2). Lower mRNA level for genes encoding CB1 and CB2 were also found in winter muscle and adipose tissue, respectively. The observed reduction in ECS tone may promote fatty acid mobilization from body fat stores, and favor carbohydrate metabolism in skeletal muscle of hibernating bears. Additionally, high circulating level of the endocannabinoid-like compound N-oleoylethanolamide (OEA) in winter could favor lipolysis and fatty acid oxidation in peripheral tissues. We also speculated on a role of OEA in the conservation of an anorexigenic signal and in the maintenance of torpor during hibernation, while sustaining the capacity of bears to sense stimuli from the environment.

Room without a view-Den excavation in relation to body size in brown bears.

Hibernation is an adaptive strategy to survive harsh winter conditions and food shortage. The use of well-insulated winter dens helps animals minimize energy loss during hibernation. Brown bears (Ursus arctos) commonly use excavated dens for hibernation. Physical attributes of excavated dens are expected to impact the bear's heat retention and energy conservation. The objective of this study was to examine the determinants of cavity size of excavated dens and the impact of physical attributes of excavated dens on energy conservation in hibernating bears, hypothesizing that bears excavate dens in a way to minimize heat loss and optimize energy conservation during hibernation. We predicted that den cavity size would be determined by the bear's body size and that older bears would excavate better-fitting cavities to minimize heat loss, due to their previous experience. We further predicted that physical attributes of excavated dens would affect the bears' posthibernation body condition. Our results revealed that bears excavated a den cavity in relation to their body size, regardless of sex, and that older bears tended to excavate better-fitting den cavities compared to young bears, as we expected. Older bears excavated better-fitting den cavities, suggesting a potentially experience-based shift with age in den-excavation behavior and an optimum cavity size relative to a bear's body size. Our key finding is that insulation of excavated dens provided by wall/rood thickness and bedding materials had a significant positive effect on bears' posthibernation body condition. We believe that our study provides new insight into how not only the quality of denning habitat, but also the quality of dens may affect hibernating animals, by presenting a potential adaptive aspect of den preparation (age effect on efficiency in den excavation) and effect of den attributes on the posthibernation body condition of brown bears.

Wolf habitat selection when sympatric or allopatric with brown bears in Scandinavia.

Habitat selection of animals depends on factors such as food availability, landscape features, and intra- and interspecific interactions. Individuals can show several behavioral responses to reduce competition for habitat, yet the mechanisms that drive them are poorly understood. This is particularly true for large carnivores, whose fine-scale monitoring is logistically complex and expensive. In Scandinavia, the home-range establishment and kill rates of gray wolves (Canis lupus) are affected by the coexistence with brown bears (Ursus arctos). Here, we applied resource selection functions and a multivariate approach to compare wolf habitat selection within home ranges of wolves that were either sympatric or allopatric with bears. Wolves selected for lower altitudes in winter, particularly in the area where bears and wolves are sympatric, where altitude is generally higher than where they are allopatric. Wolves may follow the winter migration of their staple prey, moose (Alces alces), to lower altitudes. Otherwise, we did not find any effect of bear presence on wolf habitat selection, in contrast with our previous studies. Our new results indicate that the manifestation of a specific driver of habitat selection, namely interspecific competition, can vary at different spatial-temporal scales. This is important to understand the structure of ecological communities and the varying mechanisms underlying interspecific interactions.

Correction: Diet of the brown bear in Himalaya: Combining classical and molecular genetic techniques.

[This corrects the article DOI: 10.1371/journal.pone.0225698.].

Do follicles matter? Testing the effect of follicles on hair cortisol levels.

Cortisol concentrations in hair are used increasingly as a biomarker of long-term stress in free-ranging wildlife. Cortisol is believed to be integrated into hair primarily during its active growth phase, typically occurring over weeks to months or longer periods, depending on latitude. Cortisol concentrations in hair thus reflect the activity of the hypothalamic-pituitary-adrenal axis over this time. However, local, independent cortisol secretion within the skin, which includes hair follicles, may also contribute to cortisol levels in growing hair. Methodological differences between studies include the measurement of cortisol in only the hair shaft (i.e. follicle absent, as with shaved hair) versus the whole hair (i.e. follicle present, as with plucked hair). If the concentration of cortisol in the follicle is high enough to influence the overall hair cortisol concentration (HCC), this could confound comparisons between studies using different types of hair samples (hair shafts vs. whole hair) and collection methods. Here, we test the hypothesis that cortisol present in follicles influences HCC. We compared HCC in paired subsamples of hair with and without follicles from 30 free-ranging Scandinavian brown bears (Ursus arctos) and observed significantly greater HCC in samples with follicles present. The effect of follicles remained significant also with sex and age of sampled bears taken into account in a linear mixed model. Finally, we provide an overview of collection methods and types of hair samples used for HCC analysis in 77 studies dealing with stress in wild mammal species. Our findings highlight the need to unify methods of hair collection and preparation to allow for valid comparisons, and to optimize labour input in ecophysiological studies.

MicroRNAs facilitate skeletal muscle maintenance and metabolic suppression in hibernating brown bears.

Hibernating brown bears, Ursus arctos, undergo extended periods of inactivity and yet these large hibernators are resilient to muscle disuse atrophy. Physiological characteristics associated with atrophy resistance in bear muscle have been examined (e.g., muscle mechanics, neural activity) but roles for molecular signaling/regulatory mechanisms in the resistance to muscle wasting in bears still require investigation. Using quantitative reverse transcription PCR (RT-qPCR), the present study characterized the responses of 36 microRNAs linked with development, metabolism, and regeneration of skeletal muscle, in the vastus lateralis of brown bears comparing winter hibernating and summer active animals. Relative levels of mRNA of selected genes (mef2a, pax7, id2, prkaa1, and mstn) implicated upstream and downstream of the microRNAs were examined. Results indicated that hibernation elicited a myogenic microRNA, or "myomiR", response via MEF2A-mediated signaling. Upregulation of MEF2A-controlled miR-1 and miR-206 and respective downregulation of pax7 and id2 mRNA are suggestive of responses that promote skeletal muscle maintenance. Increased levels of metabolic microRNAs, such as miR-27, miR-29, and miR-33, may facilitate metabolic suppression during hibernation via mechanisms that decrease glucose uptake and fatty acid oxidation. This study identified myomiR-mediated mechanisms for the promotion of muscle regeneration, suppression of ubiquitin ligases, and resistance to muscle atrophy during hibernation mediated by observed increases in miR-206, miR-221, miR-31, miR-23a, and miR-29b. This was further supported by the downregulation of myomiRs associated with a muscle injury and inflammation (miR-199a and miR-223) during hibernation. The present study provides evidence of myomiR-mediated signaling pathways that are activated during hibernation to maintain skeletal muscle functionality in brown bears.

Diet of the brown bear in Himalaya: Combining classical and molecular genetic techniques.

The ecological requirements of brown bears are poorly known in the Himalaya region, which complicates conservation efforts. We documented the diet of the Himalayan brown bear (Ursus arctos isabellinus) by combining classical scat analysis and a newly developed molecular genetic technique (the trnL approach), in Deosai National Park, Pakistan. Brown bears consumed over 50 plant species, invertebrates, ungulates, and several rodents. Eight plant families; Poaceae, Polygonaceae, Cyperaceae, Apiaceae, Asteraceae, Caryophyllaceae, Lamiaceae, and Rubiaceae were commonly eaten with graminoids comprising the bulk of the diet. Golden marmots comprised the major mammalian biomass in the park, and were also the main meat source for bears. Animal matter, making 36% of dietary content, contributed half of the digestible energy, due to its higher nutritious value. We did not find a significant temporal pattern in diet, perhaps because the availability of the major diet (graminoids) did not change over the foraging period. Male brown bears were more carnivorous than females, probably because of their larger size, which requires higher energy and also makes them more efficient in capturing marmots. Frequencies of three plant species were also significantly higher in male brown bears; Bistorta affinis, Carex diluta, and Carex sp. Diet of the brown bear differed significantly between the park and surrounding valleys. In valleys, diet consisted predominantly of graminoids and crops, whereas the park provided more nutritious and diverse foodThe estimated digestible energy available to brown bears in Deosai was the lowest documented among brown bear populations, due to the lack of fruits and a relatively lower meat content. The low nutritious diet and high cost of metabolism in a high-altitude environment, probably explains the very low reproductive potential of this population.

Phylogeography, genetic diversity, and connectivity of brown bear populations in Central Asia.

Knowledge of genetic diversity and population structure is critical for conservation and management planning at the population level within a species' range. Many brown bear populations in Central Asia are small and geographically isolated, yet their phylogeographic relationships, genetic diversity, and contemporary connectivity are poorly understood. To address this knowledge gap, we collected brown bear samples from the Gobi Desert (n = 2360), Altai, Sayan, Khentii, and Ikh Khyangan mountains of Mongolia (n = 79), and Deosai National Park in the Himalayan Mountain Range of Pakistan (n = 5) and generated 927 base pairs of mitochondrial DNA (mtDNA) sequence data and genotypes at 13 nuclear DNA microsatellite loci. We documented high levels of mtDNA and nDNA diversity in the brown bear populations of northern Mongolia (Altai, Sayan, Buteeliin nuruu and Khentii), but substantially lower diversity in brown bear populations in the Gobi Desert and Himalayas of Pakistan. We detected 3 brown bear mtDNA phylogeographic groups among bears of the region, with clade 3a1 in Sayan, Khentii, and Buteeliin nuruu mountains, clade 3b in Altai, Sayan, Buteeliin nuruu, Khentii, and Ikh Khyangan, and clade 6 in Gobi and Pakistan. Our results also clarified the phylogenetic relationships and divergence times with other brown bear mtDNA clades around the world. The nDNA genetic structure analyses revealed distinctiveness of Gobi bears and different population subdivisions compared to mtDNA results. For example, genetic distance for nDNA microsatellite loci between the bears in Gobi and Altai (FST = 0.147) was less than that of the Gobi and Pakistan (FST = 0.308) suggesting more recent male-mediated nuclear gene flow between Gobi and Altai than between Gobi and the Pakistan bears. Our results provide valuable information for conservation and management of bears in this understudied region of Central Asia and highlight the need for special protection and additional research on Gobi brown bears.