Factors influencing terrestriality in primates of the Americas and Madagascar.

Among mammals, the order Primates is exceptional in having a high taxonomic richness in which the taxa are arboreal, semiterrestrial, or terrestrial. Although habitual terrestriality is pervasive among the apes and African and Asian monkeys (catarrhines), it is largely absent among monkeys of the Americas (platyrrhines), as well as galagos, lemurs, and lorises (strepsirrhines), which are mostly arboreal. Numerous ecological drivers and species-specific factors are suggested to set the conditions for an evolutionary shift from arboreality to terrestriality, and current environmental conditions may provide analogous scenarios to those transitional periods. Therefore, we investigated predominantly arboreal, diurnal primate genera from the Americas and Madagascar that lack fully terrestrial taxa, to determine whether ecological drivers (habitat canopy cover, predation risk, maximum temperature, precipitation, primate species richness, human population density, and distance to roads) or species-specific traits (body mass, group size, and degree of frugivory) associate with increased terrestriality. We collated 150,961 observation hours across 2,227 months from 47 species at 20 sites in Madagascar and 48 sites in the Americas. Multiple factors were associated with ground use in these otherwise arboreal species, including increased temperature, a decrease in canopy cover, a dietary shift away from frugivory, and larger group size. These factors mostly explain intraspecific differences in terrestriality. As humanity modifies habitats and causes climate change, our results suggest that species already inhabiting hot, sparsely canopied sites, and exhibiting more generalized diets, are more likely to shift toward greater ground use.

Assessing the risk of climate maladaptation for Canadian polar bears.

The Arctic is warming four times faster than the rest of the world, threatening the persistence of many Arctic species. It is uncertain if Arctic wildlife will have sufficient time to adapt to such rapidly warming environments. We used genetic forecasting to measure the risk of maladaptation to warming temperatures and sea ice loss in polar bears (Ursus maritimus) sampled across the Canadian Arctic. We found evidence for local adaptation to sea ice conditions and temperature. Forecasting of genome-environment mismatches for predicted climate scenarios suggested that polar bears in the Canadian high Arctic had the greatest risk of becoming maladapted to climate warming. While Canadian high Arctic bears may be the most likely to become maladapted, all polar bears face potentially negative outcomes to climate change. Given the importance of the sea ice habitat to polar bears, we expect that maladaptation to future warming is already widespread across Canada.

The Space of Equidistant Phylogenetic Cactuses.

An equidistant X-cactus is a type of rooted, arc-weighted, directed acyclic graph with leaf set X, that is used in biology to represent the evolutionary history of a set X of species. In this paper, we introduce and investigate the space of equidistant X-cactuses. This space contains, as a subset, the space of ultrametric trees on X that was introduced by Gavryushkin and Drummond. We show that equidistant-cactus space is a CAT(0)-metric space which implies, for example, that there are unique geodesic paths between points. As a key step to proving this, we present a combinatorial result concerning ranked rooted X-cactuses. In particular, we show that such graphs can be encoded in terms of a pairwise compatibility condition arising from a poset of collections of pairs of subsets of X that satisfy certain set-theoretic properties. As a corollary, we also obtain an encoding of ranked, rooted X-trees in terms of partitions of X, which provides an alternative proof that the space of ultrametric trees on X is CAT(0). We expect that our results will provide the basis for novel ways to perform statistical analyses on collections of equidistant X-cactuses, as well as new directions for defining and understanding spaces of more general, arc-weighted phylogenetic networks.

Survival improvements of marine mammals in zoological institutions mirror historical advances in human longevity.

An intense public debate has fuelled governmental bans on marine mammals held in zoological institutions. The debate rests on the assumption that survival in zoological institutions has been and remains lower than in the wild, albeit the scientific evidence in support of this notion is equivocal. Here, we used statistical methods previously applied to assess historical improvements in human lifespan and data on 8864 individuals of four marine mammal species (harbour seal, Phoca vitulina; California sea lion, Zalophus californianus; polar bear, Ursus maritimus; common bottlenose dolphin, Tursiops truncatus) held in zoos from 1829 to 2020. We found that life expectancy increased up to 3.40 times, and first-year mortality declined up to 31%, during the last century in zoos. Moreover, the life expectancy of animals in zoos is currently 1.65-3.55 times longer than their wild counterparts. Like humans, these improvements have occurred concurrently with advances in management practices, crucial for population welfare. Science-based decisions will help effective legislative changes and ensure better implementation of animal care.

Assessing the effectiveness of protected areas for panda conservation under future climate and land use change scenarios.

While the relatively stable land use and land cover (LULC) patterns is an important feature of protected areas (PAs), the influence of this feature on future species distribution and the effectiveness of the PAs has rarely been explored. Here, we assessed the role of land use patterns within PAs on the projected range of the giant panda (Ailuropoda melanoleuca) by comparing projections inside and outside of PAs for four model configurations: (1) only climate covariates, (2) climate and dynamic land use covariates, (3) climate and static land use covariates and (4) climate and hybrid dynamic-static land use covariates. Our objectives were twofold: to understand the role of protected status on projected panda habitat suitability and evaluate the relative efficacy of different climate modeling approaches. The climate and land use change scenarios used in the models include two shared socio-economic pathways (SSPs) scenarios: SSP126 [an optimistic scenario] and SSP585 [a pessimistic scenario]. We found that models including land-use covariates performed significantly better than climate-only models and that these projected more suitable habitat than climate-only models. Static land-use models projected more suitable habitat than both the dynamic and hybrid models under SSP126, while these models did not differ under SSP585. China's panda reserve system was projected to effectively maintain suitable habitat inside PAs. Panda dispersal ability also significantly impacted outcomes, with most models assuming unlimited dispersal forecasting range expansion and models assuming zero dispersal consistently forecasting range contraction. Our findings highlight that policies targeting improved land-use practices should be an effective means for offsetting some of the negative effects of climate change on pandas. As the effectiveness of PAs is projected to be maintained, we recommend the judicious management and expansion of the PA system to ensure the resilience of panda populations into the future.

Long-term increases in pathogen seroprevalence in polar bears (Ursus maritimus) influenced by climate change.

The influence of climate change on wildlife disease dynamics is a burgeoning conservation and human health issue, but few long-term studies empirically link climate to pathogen prevalence. Polar bears (Ursus maritimus) are vulnerable to the negative impacts of sea ice loss as a result of accelerated Arctic warming. While studies have associated changes in polar bear body condition, reproductive output, survival, and abundance to reductions in sea ice, no long-term studies have documented the impact of climate change on pathogen exposure. We examined 425 serum samples from 381 adult polar bears, collected in western Hudson Bay (WH), Canada, for antibodies to selected pathogens across three time periods: 1986-1989 (n = 157), 1995-1998 (n = 159) and 2015-2017 (n = 109). We ran serological assays for antibodies to seven pathogens: Toxoplasma gondii, Neospora caninum, Trichinella spp., Francisella tularensis, Bordetella bronchiseptica, canine morbillivirus (CDV) and canine parvovirus (CPV). Seroprevalence of zoonotic parasites (T. gondii, Trichinella spp.) and bacterial pathogens (F. tularensis, B. bronchiseptica) increased significantly between 1986-1989 and 1995-1998, ranging from +6.2% to +20.8%, with T. gondii continuing to increase into 2015-2017 (+25.8% overall). Seroprevalence of viral pathogens (CDV, CPV) and N. caninum did not change with time. Toxoplasma gondii seroprevalence was higher following wetter summers, while seroprevalences of Trichinella spp. and B. bronchiseptica were positively correlated with hotter summers. Seroprevalence of antibodies to F. tularensis increased following years polar bears spent more days on land, and polar bears previously captured in human settlements were more likely to be seropositive for Trichinella spp. As the Arctic has warmed due to climate change, zoonotic pathogen exposure in WH polar bears has increased, driven by numerous altered ecosystem pathways.

Mean and Variance of Phylogenetic Trees.

We describe the use of the Fréchet mean and variance in the Billera-Holmes-Vogtmann (BHV) treespace to summarize and explore the diversity of a set of phylogenetic trees. We show that the Fréchet mean is comparable to other summary methods, and, despite its stickiness property, is more likely to be binary than the majority-rule consensus tree. We show that the Fréchet variance is faster and more precise than commonly used variance measures. The Fréchet mean and variance are more theoretically justified, and more robust, than previous estimates of this type and can be estimated reasonably efficiently, providing a foundation for building more advanced statistical methods and leading to applications such as mean hypothesis testing and outlier detection.

Bordetella bronchiseptica-reactive antibodies in Canadian polar bears.

Bordetella bronchiseptica is a promiscuous bacterium that infects a variety of species but has not been reported in free-ranging polar bears (Ursus maritimus). Sera from 385 polar bears from the western Hudson Bay region, 1986 to 2017, were tested for reactivity to B. bronchiseptica with enzyme-linked immunosorbent assays using anti-canine IgG and Streptococcus protein G as secondary reagents. Sera from bears had variable reactivity to B. bronchiseptica antigens, and there was no difference among bears that had a history of coming near the town of Churchill, Manitoba, and bears that did not. Although the sources of exposure were not determined, equivalent results in both groups suggest that potential exposure to humans (aside from handling during sampling) and their animals (dogs) was not an important co-factor in sero-positivity to B. bronchiseptica.

Anticorps réactifs à Bordetella bronchiseptica chez les ours polaires du Canada. Bordetella bronchiseptica est une bactérie qui infecte une variété d'espèces mais qui n'a pas été signalée chez les ours polaires en liberté (Ursus maritimus). Les sérums de 385 ours polaires de la région ouest de la baie d'Hudson, de 1986 à 2017, ont été testés pour leur réactivité à B. bronchiseptica par une épreuve ELISA utilisant des anticorps anti-IgG canines et de la protéine G de Streptococcus comme réactifs secondaires. Les sérums d'ours avaient une réactivité variable aux antigènes de B. bronchiseptica, et il n'y avait aucune différence entre les ours qui avaient l'habitude de s'approcher de la ville de Churhill, au Manitoba, et les ours qui n'en avaient pas. Bien que les sources d'exposition n'aient pas été déterminées, des résultats équivalents dans les deux groupes suggèrent que l'exposition potentielle des humains (en dehors de la manipulation pendant l'échantillonnage) et de leurs animaux (chiens) n'était pas un cofacteur important de la séropositivité à B. bronchiseptica.(Traduit par Dr Serge Messier).

Climate change and landscape-use patterns influence recent past distribution of giant pandas.

Climate change is one of the most pervasive threats to biodiversity globally, yet the influence of climate relative to other drivers of species depletion and range contraction remain difficult to disentangle. Here, we examine climatic and non-climatic correlates of giant panda (Ailuropoda melanoleuca) distribution using a large-scale 30 year dataset to evaluate whether a changing climate has already influenced panda distribution. We document several climatic patterns, including increasing temperatures, and alterations to seasonal temperature and precipitation. We found that while climatic factors were the most influential predictors of panda distribution, their importance diminished over time, while landscape variables have become relatively more influential. We conclude that the panda's distribution has been influenced by changing climate, but conservation intervention to manage habitat is working to increasingly offset these negative consequences.

Comprehensive Breeding Techniques for the Giant Panda.

The dramatic growth of the captive giant panda (Ailuropoda melanoleuca) population exemplifies how the application of scientific findings to animal care and reproductive management can improve conservation breeding outcomes. Detailed behavioral studies of giant panda estrus, pregnancy and cub rearing have demonstrated the importance of husbandry management that supports natural reproductive behavior to enhance breeding success. Natural breeding has been valuably augmented by the development of assisted reproductive techniques founded through detailed studies of the reproductive physiology of the giant panda and outlining fundamental information about reproductive seasonality, male fertility and characterization of the estrous cycle. The resultant holistic understanding of giant panda reproduction has improved reproductive success in the captive population to such an extent that it is now self-sustaining and provides surplus animals for reintroduction. Despite these significant advances, there are knowledge gaps and remaining challenges to be addressed. Pregnancy detection remains the single biggest challenge when determining if natural mating or assisted breeding have been successful. Because pregnancy can only be determined in the few weeks prior to parturition, there are gaps in understanding and detecting delayed implantation and early embryonic loss. Additionally, dynamic management practices and standard of care for reproductive assistance needs to be developed. Only large breeding centers in China have the ability to promote normal reproductive behaviors and allow mate choice for the giant panda. These challenges need to be addressed in the near future in order to maintain a self-sustaining, genetically diverse and behaviorally competent captive population. This chapter documents the development of successful giant panda managed breeding programs by focusing on three key areas, (1) the development of science-driven reproductive techniques to improve fecundity in a species where the mating system was poorly understood, (2) how targeted research and adaptive management of social settings surrounding estrus and breeding improved reproductive success, and (3) insights and solutions to challenges faced across the program's history with future directions for research.

The Andean bear alopecia syndrome may be caused by social housing.

The Andean bear alopecia syndrome is a progressive and chronic condition documented in ex situ populations. Recent advances focus on treating symptoms, not preventing future cases. We therefore explored the epidemiology of this syndrome through an analysis of husbandry and veterinary conditions of 63 Andean bears (26M:37F) housed in North and South American zoos and other ex situ circumstances. We had the most complete information for the North American population and found that 29% of females (n = 24) were affected. No males (n = 26) were affected. An analysis of generalized linear models indicated that three models were competitive in describing the occurrence of the condition (i.e., ΔAICc ≤ 2): the model including only the individual's sex (χ2 = 13.41, df = 1, p < .001), the model including both individual sex and social housing status (χ2 = 1.36, df = 2, p < .001), and the model including both individual sex and the expression of stereotypical behaviors (χ2 = 13.82, df = 2, p = .001). Stereotypical behaviors were common among both males (50%, n = 26) and females (51.9%, n = 27) whether or not they were affected, but the syndrome was seen only in females who had been socially housed. Therefore, we suggest that the Andean bear alopecia syndrome is a symptomatic response to the long-term social housing of bears that would otherwise not live socially. To prevent new cases, we recommend that female Andean bears be housed with adult conspecifics only when females choose to cohabitate.

Energetic costs of locomotion in bears: is plantigrade locomotion energetically economical?

Ursids are the largest mammals to retain a plantigrade posture. This primitive posture has been proposed to result in reduced locomotor speed and economy relative to digitigrade and unguligrade species, particularly at high speeds. Previous energetics research on polar bears (Ursus maritimus) found locomotor costs were more than double predictions for similarly sized quadrupedal mammals, which could be a result of their plantigrade posture or due to adaptations to their Arctic marine existence. To evaluate whether polar bears are representative of terrestrial ursids or distinctly uneconomical walkers, this study measured the mass-specific metabolism, overall dynamic body acceleration, and gait kinematics of polar bears and grizzly bears (Ursus arctos) trained to rest and walk on a treadmill. At routine walking speeds, we found polar bears and grizzly bears exhibited similar costs of locomotion and gait kinematics, but differing measures of overall dynamic body acceleration. Minimum cost of transport while walking in the two species (2.21 J kg-1 m-1) was comparable to predictions for similarly sized quadrupedal mammals, but these costs doubled (4.42 J kg-1 m-1) at speeds ≥5.4 km h-1 Similar to humans, another large plantigrade mammal, bears appear to exhibit a greater economy while moving at slow speeds.

Approaches to studying behavior in captive sloth bears through animal keeper feedback.

Animal keepers at zoos and wildlife rescue centers often possess in-depth knowledge of the health and behavior of the individuals under their care. While it is often not feasible for keepers to regularly collect behavior data through formal scientific methods, efforts should be made to find alternative means to capture this knowledge. We investigated the use of keeper feedback to study the behavior of sloth bears at the Agra Bear Rescue Facility (ABRF; Agra, India). We prepared a survey with 5 questions focused on behaviors indicative of playfulness, boldness, aggressiveness, and the tendency to express self-directed behaviors (SDB). We asked keepers to rate bears on a Likert scale from 1 (least likely to exhibit a behavior) to 5 (most likely) for 44 adult female bears (5-21 years of age). We validated this method by comparing keeper ratings of SDB with formal behavior observations, finding that time of day had an influence on the accuracy of keeper assessments. We found a significant negative correlation between housing bears in larger groups (>15) and SDB. In addition, we correlated ratings given by keepers for all study behaviors. Social play had significant negative correlation with aggression toward people. There was no correlation between social play and aggression toward other bears, possibly due to the existence of cohesive social groups in group housing or high dimensionality of the data. We found that keeper feedback is an efficient tool to gather behavior data on captive sloth bears and recommend its use in future studies.

A systematic survey of the integration of animal behavior into conservation.

The role of behavioral ecology in improving wildlife conservation and management has been the subject of much recent debate. We sought to answer 2 foundational questions about the current use of behavioral knowledge in conservation: To what extent is behavioral knowledge used in wildlife conservation and management, and how does the use of animal behavior differ among conservation fields in both frequency and types of use? We searched the literature for intersections between key fields of animal behavior and conservation and created a systematic heat map (i.e., graphical representation of data where values are represented as colors) to visualize relative efforts. Some behaviors, such as dispersal and foraging, were commonly considered (mean [SE] of 1147.38 [353.11] and 439.44 [108.85] papers per cell, respectively). In contrast, other behaviors, such as learning, social, and antipredatory behaviors were rarely considered (mean [SE] of 33.88 [7.62], 44.81 [10.65], and 22.69 [6.37] papers per cell, respectively). In many cases, awareness of the importance of behavior did not translate into applicable management tools. Our results challenge previous suggestions that there is little association between the fields of behavioral ecology and conservation and reveals tremendous variation in the use of different behaviors in conservation. We recommend that researchers focus on examining underutilized intersections of behavior and conservation themes for which preliminary work shows a potential for improving conservation and management, translating behavioral theory into applicable and testable predictions, and creating systematic reviews to summarize the behavioral evidence within the behavior-conservation intersections for which many studies exist.

On Determining if Tree-based Networks Contain Fixed Trees.

We address an open question of Francis and Steel about phylogenetic networks and trees. They give a polynomial time algorithm to decide if a phylogenetic network, N, is tree-based and pose the problem: given a fixed tree T and network N, is N based on T? We show that it is [Formula: see text]-hard to decide, by reduction from 3-Dimensional Matching (3DM) and further that the problem is fixed-parameter tractable.

The influence of ambient noise on maternal behavior in a Bornean sun bear (Helarctos malayanus euryspilus).

Anthropogenic noise has become a pervasive feature of both marine and terrestrial habitats worldwide. While a comprehensive understanding of the biologically significant impacts of noise on wildlife is lacking, concerns regarding its influence persist. Noise is also a common feature in the zoological setting, and much can be learned regarding the species-typical response to ambient noise by studying animals in captivity. Here we correlate behavioral and vocal patterns in a Bornean sun bear (Helarctos malayanus euryspilus) mother and cub with ambient noise levels during the 6-month post-partum period. We hypothesized that loud ambient noise would be correlated with changes in behavior, and predicted that noise would be negatively correlated with maternal care behavior, potentially masking cub vocalizations or providing a distraction to the mother. Contrary to expectation, we found that the mother spent significantly more time attending to her cub (P=0.03) on loud days. We also found that she tended to spend less time feeding (P=0.08); however her time spent resting was not affected. The cub was approximately twice as vocal on loud days, although these results were not statistically significant (humming: P=0.10; squawks/cries: P=0.14). Taken together, these results suggest that the behavioral response to ambient noise may have potential energetic costs, and as a result efforts should be made to reduce ambient noise exposure during the post-partum period.

Geodesics to characterize the phylogenetic landscape.

Phylogenetic trees are fundamental for understanding evolutionary history. However, finding maximum likelihood trees is challenging due to the complexity of the likelihood landscape and the size of tree space. Based on the Billera-Holmes-Vogtmann (BHV) distance between trees, we describe a method to generate intermediate trees on the shortest path between two trees, called pathtrees. These pathtrees give a structured way to generate and visualize part of treespace. They allow investigating intermediate regions between trees of interest, exploring locally optimal trees in topological clusters of treespace, and potentially finding trees of high likelihood unexplored by tree search algorithms. We compared our approach against other tree search tools (Paup*, RAxML, and RevBayes) using the highest likelihood trees and number of new topologies found, and validated the accuracy of the generated treespace. We assess our method using two datasets. The first consists of 23 primate species (CytB, 1141 bp), leading to well-resolved relationships. The second is a dataset of 182 milksnakes (CytB, 1117 bp), containing many similar sequences and complex relationships among individuals. Our method visualizes the treespace using log likelihood as a fitness function. It finds similarly optimal trees as heuristic methods and presents the likelihood landscape at different scales. It found relevant trees that were not found with MCMC methods. The validation measures indicated that our method performed well mapping treespace into lower dimensions. Our method complements heuristic search analyses, and the visualization allows the inspection of likelihood terraces and exploration of treespace areas not visited by heuristic searches.

Effects of administering exogenous bovine somatotropin to beef heifers during the first trimester on conceptus development as well as steroid- and eicosanoid-metabolizing enzymes.

The aim of this study was to evaluate the effect of bovine somatotropin (bST) on fetal and placental development during the first third of gestation in beef heifers. Angus heifers (n = 97) were randomly assigned to either receive a 500-mg injection of bST (BST) biweekly on days 0, 15, 29, 43, and 57 of gestation or not receive bST (CTL) throughout the experiment. Body weight (BW) was assessed on days -9, -3, 0, 15, 22, 29, 43, 50, 57, 64, and 77, while blood samples were collected on days 0, 22, 50, and 64. Pregnancy status was determined via transrectal ultrasonography on days 29 and 64. A subset of pregnant heifers (BST, n = 7; CTL, n = 5) were harvested on day 84, and complete gravid reproductive tracts and liver tissue were collected for analysis. Cytochrome P450 1A (CYP1A), 2C (CYP2C), 3A (CYP3A), and uridine 5'-diphospho-glucuronosyltransferase (UGT) activities were determined. Mean change in BW and average daily gain of heifers between fixed-time artificial insemination (day 0) and day 77 did not differ between treatments (P ≥ 0.05). Mean concentrations of insulin-like growth factor 1 (IGF-1) were greater (P < 0.001) in BST (347 ± 27.7 ng/mL) compared with CTL (135 ± 32.8 ng/mL) heifers. Mean placental weight, fetal membrane weight, uterine weight, and ovarian and corpus luteum (CL) weights, as well as fetal morphometric data, did not differ (P ≥ 0.05) between treatments. However, BST heifers had greater (P = 0.03) quantities of combined fetal fluid compared with CTL (521.6 ± 22.9 vs. 429.6 ± 27.14 g, respectively). Tendencies were observed for BST heifers to have reproductive tracts with fewer placentomes (P = 0.08) and fetuses with greater umbilical diameters (P = 0.09) compared with CTL. The activity of CYP1A did not differ (P ≥ 0.05) within the maternal and fetal liver, caruncle, cotyledon, or CL tissue samples between treatments. Furthermore, CYP3A activity was only observed in maternal liver samples and was not different between treatments (P ≥ 0.05). Interestingly, CYP2C activity was greater (P = 0.01) in the liver of BST vs. CTL heifers, and UGT activity was greater (P = 0.02) in the CL from BST heifers compared with CTL. In conclusion, the administration of bST during the first third of gestation increased plasma concentrations of IGF-1, which resulted in an increase in fetal fluid, decrease in placentome number, and greater umbilical diameter, but failed to alter fetal development.

Maximum Covering Subtrees for Phylogenetic Networks.

Tree-based phylogenetic networks, which may be roughly defined as leaf-labeled networks built by adding arcs only between the original tree edges, have elegant properties for modeling evolutionary histories. We answer an open question of Francis, Semple, and Steel about the complexity of determining how far a phylogenetic network is from being tree-based, including non-binary phylogenetic networks. We show that finding a phylogenetic tree covering the maximum number of nodes in a phylogenetic network can be computed in polynomial time via an encoding into a minimum-cost flow problem.

Stereotypic behaviour predicts reproductive performance and litter sex ratio in giant pandas.

Breeding and welfare problems confront many conservation breeding programs. Stereotypies-repetitive, unvarying, functionless behaviours -are common abnormal behaviours that often arise in suboptimal conditions. While the role of stereotypies in welfare assessment is well studied, few investigations address the relationship between stereotypic behaviour and reproduction. We examined the correlation between stereotypic behaviour and reproductive performance in 101 giant pandas (Ailuropoda melanoleuca). High stereotyping males copulated more and produced more cubs, suggesting that highly sexually motivated males were prone to stereotypy but also had high reproductive competence. Female stereotypies were negatively associated with all reproductive measures closely tied to behavioural competence: high stereotyping females were less likely to copulate, less likely to mother-rear cubs, and-probably a result of poor maternal care-had lower cub survival. However, females that exhibited stereotypies were more likely to produce a cub, suggesting stereotypies are tied to behavioural but not physiological competence. High stereotyping female pandas also displayed strong and consistent bias toward production of female offspring while paternal relationship to sex allocation was the reverse. These results are consistent with stress-mediated sex allocation theory. Our findings raise concern about differential reproductive success among high and low stereotyping pandas, and possible genetic adaptation to captivity.