Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates.

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17ß-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Archaic Introgression Shaped Human Circadian Traits.

When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.

Sequence-based machine learning reveals 3D genome differences between bonobos and chimpanzees.

Phenotypic divergence between closely related species, including bonobos and chimpanzees (genus Pan), is largely driven by variation in gene regulation. The 3D structure of the genome mediates gene expression; however, genome folding differences in Pan are not well understood. Here, we apply machine learning to predict genome-wide 3D genome contact maps from DNA sequence for 56 bonobos and chimpanzees, encompassing all five extant lineages. We use a pairwise approach to estimate 3D divergence between individuals from the resulting contact maps in 4,420 1 Mb genomic windows. While most pairs were similar, ∼17% were predicted to be substantially divergent in genome folding. The most dissimilar maps were largely driven by single individuals with rare variants that produce unique 3D genome folding in a region. We also identified 89 genomic windows where bonobo and chimpanzee contact maps substantially diverged, including several windows harboring genes associated with traits implicated in Pan phenotypic divergence. We used in silico mutagenesis to identify 51 3D-modifying variants in these bonobo-chimpanzee divergent windows, finding that 34 or 66.67% induce genome folding changes via CTCF binding motif disruption. Our results reveal 3D genome variation at the population-level and identify genomic regions where changes in 3D folding may contribute to phenotypic differences in our closest living relatives.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates.

SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17ß-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning.

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Machine learning reveals the diversity of human 3D chromatin contact patterns.

Understanding variation in chromatin contact patterns across human populations is critical for interpreting non-coding variants and their ultimate effects on gene expression and phenotypes. However, experimental determination of chromatin contacts at a population-scale is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the diversity 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We then apply this approach to thousands of diverse modern humans and the inferred human-archaic hominin ancestral genome. While patterns of 3D contact divergence genome-wide are qualitatively similar to patterns of sequence divergence, we find that 3D divergence in local 1-megabase genomic windows does not follow sequence divergence. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, 26% of genomic windows have rare 3D contact variation observed in a small number of individuals. Using in silico mutagenesis we find that most sequence changes to do not result in changes to 3D chromatin contacts. However in windows with substantial 3D divergence, just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals diverse contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

A comparison of faecal glucocorticoid metabolite concentration and gut microbiota diversity in bonobos (Pan paniscus).

Sex, age, diet, stress and social environment have all been shown to influence the gut microbiota. In several mammals, including humans, increased stress is related to decreasing gut microbial diversity and may differentially impact specific taxa. Recent evidence from gorillas shows faecal glucocorticoid metabolite concentration (FGMC) did not significantly explain gut microbial diversity, but it was significantly associated with the abundance of the family Anaerolineaceae. These patterns have yet to be examined in other primates, like bonobos (Pan paniscus). We compared FGMC to 16S rRNA amplicons for 202 bonobo faecal samples collected across 5 months to evaluate the impact of stress, measured with FGMC, on the gut microbiota. Alpha diversity measures (Chao's and Shannon's indexes) were not significantly related to FGMC. FGMC explained 0.80 % of the variation in beta diversity for Jensen-Shannon and 1.2% for weighted UniFrac but was not significant for unweighted UniFrac. We found that genus SHD-231, a member of the family Anaerolinaceae had a significant positive relationship with FGMC. These results suggest that bonobos are relatively similar to gorillas in alpha diversity and family Anaerolinaceae responses to FGMC, but different from gorillas in beta diversity. Members of the family Anaerolinaceae may be differentially affected by FGMC across great apes. FGMC appears to be context dependent and may be species-specific for alpha and beta diversity but this study provides an example of consistent change in two African apes. Thus, the relationship between physiological stress and the gut microbiome may be difficult to predict, even among closely related species.

Predicting Archaic Hominin Phenotypes from Genomic Data.

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes.

Estimating bonobo (Panpaniscus) and chimpanzee (Pantroglodytes) evolutionary history from nucleotide site patterns.

Admixture appears increasingly ubiquitous in the evolutionary history of various taxa, including humans. Such gene flow likely also occurred among our closest living relatives: bonobos (Pan paniscus) and chimpanzees (Pan troglodytes). However, our understanding of their evolutionary history has been limited by studies that do not consider all Pan lineages or do not analyze all lineages simultaneously, resulting in conflicting demographic models. Here, we investigate this gap in knowledge using nucleotide site patterns calculated from whole-genome sequences from the autosomes of 71 bonobos and chimpanzees, representing all five extant Pan lineages. We estimated demographic parameters and compared all previously proposed demographic models for this clade. We further considered sex bias in Pan evolutionary history by analyzing the site patterns from the X chromosome. We show that 1) 21% of autosomal DNA in eastern chimpanzees derives from western chimpanzee introgression and that 2) all four chimpanzee lineages share a common ancestor about 987,000 y ago, much earlier than previous estimates. In addition, we suggest that 3) there was male reproductive skew throughout Pan evolutionary history and find evidence of 4) male-biased dispersal from western to eastern chimpanzees. Collectively, these results offer insight into bonobo and chimpanzee evolutionary history and suggest considerable differences between current and historic chimpanzee biogeography.

A Test of Foraging Models Using Dietary Diversity Indices for the Lomako Forest Bonobos.

Optimal diet and functional response models are used to understand the evolution of primate foraging strategies. The predictions of these models can be tested by examining the geographic and seasonal variation in dietary diversity. Dietary diversity is a useful tool that allows dietary comparisons across differing sampling locations and time periods. Bonobos (Pan paniscus) are considered primarily frugivorous and consume fruits, leaves, insects, vertebrates, terrestrial herbaceous vegetation, and flowers. Frugivores, like bonobos, are valuable for examining dietary diversity and testing foraging models because they eat a variety of species and are subject to seasonal shifts in fruit availability. Frugivorous primate species thus allow for tests of how variation in dietary diversity is correlated with variation in ecological factors. We investigated measures of dietary diversity in bonobos at two research camps across field seasons within the same protected area (N'dele and Iyema) in Lomako Forest, Democratic Republic of the Congo. We compared the results of behavioral observation (1984/1985, 1991, 1995, 2014, and 2017) and fecal washing analysis (2007 and 2009) between seasons and study period using three diversity indices (Shannon's, Simpson's, and SW evenness). The average yearly dietary diversity indices at N'dele were Shannon's H' = 2.04, Simpson's D = 0.82, and SW evenness = 0.88 while at Iyema, the indices were Shannon's H' = 2.02, Simpson's D = 0.82, and SW evenness = 0.88. Behavioral observation data sets yielded significantly higher dietary diversity indices than fecal washing data sets. We found that food item (fruit, leaf, and flower) consumption was not associated with seasonal food availability for the 2017 behavioral observation data set. Shannon's index was lower during periods when fewer bonobo dietary items were available to consume and higher when fruit was abundant. Finally, we found that optimal diet models best-explained patterns of seasonal food availability and dietary diversity. Dietary diversity is an essential factor to consider when understanding primate diets and can be a tool in understanding variation in primate diets, particularly among frugivores. Dietary diversity varies across populations of the same species and across time, and it is critical in establishing a complete understanding of how primate diets change over time.

Age, sex, and inflammatory markers predict chronic conditions, cardiac disease, and mortality among captive western lowland gorillas (Gorilla gorilla gorilla).

In humans, inflammatory markers predict health risks. As great apes experience many similar conditions, measuring inflammation may provide valuable health information. We examined four serum inflammatory markers in zoo-housed gorillas (n = 48): albumin, CRP, IL-6, and TNF-α. We first analyzed age- and sex-associated patterns, then used multimodel inference to evaluate models with age, sex, and inflammatory markers as predictors of all-cause morbidity, cardiac disease, and mortality. Older gorillas had lower albumin and higher IL-6, and males had higher albumin, lower CRP, and lower TNF-α. All-cause morbidity was best predicted by age, sex, and TNF-α, but the second model containing only age and sex was equivalent. Cardiac disease was best predicted by TNF-α alongside age and sex, with lower levels associated with increased risk. When outliers were removed, the model with TNF-α was second to the model containing only age and sex. Finally, mortality risk was best predicted by the model with only age and sex. Other models containing individual inflammatory markers were within top model sets for each health outcome. Our results indicate that age and sex are robust for predicting all-cause morbidity and mortality risk in gorillas; while models which include individual inflammatory markers also predict risk, they may not improve predictions over age and sex alone. However, given the prevalence of cardiac disease in great apes, these results suggest that TNF-α warrants further investigation. With their potential to provide valuable health information, data on inflammatory markers may contribute to the care and management of gorillas in human care.

Sex differences in bonobo (Pan paniscus) terrestriality: implications for human evolution.

Recent finds in hominin fossil environments place the transition to terrestriality in a wooded or forested habitat. Therefore, forest-dwelling apes can aid in understanding this important evolutionary transition. Sex differences in ape locomotion have been previously attributed to sexual dimorphism or ecological niche differences between males and females. This study examined the hypothesis that differential advantages of terrestrial travel may impact mating success in male bonobos. We examined whether males are more terrestrial when there are mating benefits for fast travel. We analyzed behavioral data on wild bonobos over a ten-month period in the Lomako Forest, DRC and examined the proportion of time spent at lower heights compared to higher heights between adult females and males relative to their location to feeding contexts with high mating frequencies. We found a significant interaction between sex and height class away from food patches (F=4.65, df =1, p <0.05) such that females were primarily arboreal whereas there was no difference between males across height classes. However, there was also a significant interaction between sex and height class (F =29.35, df =1, p <0.0001) for adults traveling near or entering a food patch. Males often arrived at food patches terrestrially and females arrived almost exclusively arboreally. We found a significant difference between the expected and observed distribution of matings by food patch context (G =114.36, df =4, p <0.0001) such that most mating occurred near or in a food patch. These results suggest that males may travel terrestrially to arrive at food patches before cohesive parties of females arrive arboreally, in order to compete with other males for mating access to these females. Such intrasexual selection for sex differences in locomotion may be important in considerations of the evolution of locomotion strategies in hominins in a forested environment.

Variation in Adult Male-Juvenile Affiliative Behavior in Japanese Macaques (Macaca fuscata).

Adult males of some primate species are known to positively interact with juveniles. In cases where paternal certainty is high, these behaviors have been largely attributed to the paternal investment hypothesis. Males have also been observed to interact with nonkin juveniles, which has often been explained in terms of mating effort. Here, we examined variation in adult male-juvenile affiliation in semi-free ranging Japanese macaques (Macaca fuscata) at the Oregon National Primate Research Center against possible influencing factors such as age, dominance rank, and female affiliation, and we also tested for fitness benefits. We conducted 154 h of focal observations of 14 adult males from June to September 2018. Males differed significantly in their rate of juvenile-directed affiliation, but not in their fitness in terms of number of offspring. There was a significant positive correlation between rank and age in the group, indicating that, in this group, rank does not conform to the classic inverted-U pattern observed elsewhere in this species. Although there was a significant positive correlation between rank and juvenile-directed affiliation, the highest-ranking male had few offspring and exhibited little juvenile-directed affiliation. These results suggest little to no preliminary support for either the paternal investment or mating effort hypotheses as explanations for juvenile-directed affiliation. This study suggests that there are multiple behavioral strategies for older males that may influence reproductive success.

Abundance, density, and social structure of African forest elephants (Loxodonta cyclotis) in a human-modified landscape in southwestern Gabon.

Population monitoring is critical to effective conservation, but forest living taxa can be difficult to directly observe. This has been true of African forest elephants (Loxodonta cyclotis), for which we have limited information regarding population size and social behavior despite their threatened conservation status. In this study, we estimated demographic parameters using genetic capture-recapture of forest elephants in the southern Industrial Corridor of the Gamba Complex of Protected Areas in southwestern Gabon, which is considered a global stronghold for forest elephants. Additionally, we examined social networks, predicting that we would find matrilineal structure seen in both savanna and forest elephants. Given 95% confidence intervals, we estimate population size in the sampled area to be between 754 and 1,502 individuals and our best density estimate ranges from 0.47 to 0.80 elephants per km2. When extrapolated across the entire Industrial Corridor, this estimate suggests an elephant population size of 3,033 to 6,043 based on abundance or 1,684 to 2,832 based on density, approximately 40-80% smaller than previously suggested. Our social network analysis revealed approximately half of network components included females with different mitochondrial haplotypes suggesting a wider range of variation in forest elephant sociality than previously thought. This study emphasizes the threatened status of forest elephants and demonstrates the need to further refine baseline estimates of population size and knowledge on social behavior in this taxon, both of which will aid in determining how population dynamics in this keystone species may be changing through time in relation to increasing conservation threats.

Social hair plucking is a grooming convention in a group of captive bonobos (Pan paniscus).

Hair plucking is observed in many captive primate species and is often characterized as an abnormal behavior. However, this behavior may be both self-directed and social and may have different etiologies. Early research in captive macaques (Macaca mulatta) described the aggressive nature of social hair plucking while more recent observations did not find an association with aggression or grooming, but the behavior was initiated most frequently by individuals with more secure dominance rank. Here, we investigate patterns of social hair plucking in a group of captive bonobos at the Columbus Zoo. We tested the hypothesis that social plucking reflects the dominance hierarchy by examining the association between social plucking and grooming, dominance, and kinship. We collected 128 h of grooming data on 16 captive bonobos using all-occurrence sampling. We ran three Mantel tests between a directed grooming matrix and (1) a plucking matrix, (2) a matrix reflecting dominance, and (3) matrix of relatedness. Grooming and hair plucking were significantly correlated (r = 0.25, p < 0.01), however, there was no association between plucking and dominance (r = - 0.04, p = 0.67), or plucking and relatedness (r = 0.07, p = 0.24). These results support the hypothesis that social plucking in bonobos is a grooming convention and is unrelated to dominance.

New Observations of Meat Eating and Sharing in Wild Bonobos (Pan paniscus) at Iyema, Lomako Forest Reserve, Democratic Republic of the Congo.

Bonobos (Pan paniscus) consume a variety of vertebrates, although direct observations remain relatively rare compared to chimpanzees (Pan troglodytes). We report the first direct observations of meat eating and sharing among bonobos at Iyema, Lomako Forest, Democratic Republic of Congo. We collected meat consumption data ad libitum from June to November 2017 over 176.5 observation hours and conducted monthly censuses to measure the abundance of potential prey species. We observed 3 occasions of duiker consumption and found indirect evidence of meat consumption twice (n = 5). We identified the prey species as Weyn's duiker (Cephalophus weynsi) in all 4 cases that we saw the carcass. This species was the most abundant duiker species at Iyema, but other potential prey species were also available. Meat sharing was observed or inferred during all 3 observations. However, the individual controlling the carcass frequently resisted sharing, and aggressive attempts to take the carcass were observed. This report contributes to a growing body of data suggesting that wild bonobos consume meat at higher rates than previously thought, female control of carcasses is frequent but not exclusive, and meat sharing in bonobos is primarily passive but not without aggression.

Prevalence and characteristics of hair plucking in captive bonobos (Pan paniscus) in North American zoos.

When primates exhibit hair loss and are observed to engage in self or social hair plucking (a rapid jerking away of the hair shaft and follicle by the hand or mouth, often accompanied by inspection, and consumption) the altered appearance, and behavior patterns are thought to reflect individual physiological, and psychological well-being. Hair loss and hair plucking occur in many captive primate species, including all of the great apes. We present the first survey of this behavior among captive bonobos (N = 88; 50 females and 38 males) in seven zoos in the United States. We found that 43% of the population engaged in this behavior pattern and discounting youngsters (who are not observed to hair pluck until the age of five), 58% of individuals hair plucked. Of the individuals who hair plucked, 97% engaged in social plucking, whereas 46% engaged in self-directed plucking. We regressed the occurrence of hair plucking with multiple predictor variables using binary logistic regression and multimodel inference to determine which predictors best explained the prevalence of self-directed and social plucking. We also analyzed publicly available data on hair plucking in captive chimpanzees. We found that the occurrence of another abnormal behavior, age, origin, and pelage condition best explained self-directed plucking in bonobos. Social plucking was explained by age, origin, pelage, and sex. Our analysis of chimpanzee hair plucking revealed that plucking is strongly influenced by rearing and sex. This study demonstrates that hair plucking is more prevalent in captive bonobos compared to captive chimpanzees and gorillas, however, the covariates associated with hair plucking for each species are different. Our data suggest a potential link between self-directed plucking and well-being. However, the higher prevalence and etiology of social hair plucking is more difficult to explain.

Laterality of Grooming and Tool Use in a Group of Captive Bonobos (Pan paniscus).

Humans exhibit population level handedness for the right hand; however, the evolution of this behavioral phenotype is poorly understood. Here, we compared the laterality of a simple task (grooming) and a complex task (tool use) to investigate whether increasing task difficulty elicited individual hand preference among a group of captive bonobos (Pan paniscus). Subjects were 17 bonobos housed at the Columbus Zoo and Aquarium. Laterality of grooming was recorded using group scans; tool use was recorded using all-occurrence sampling. Grooming was characterized as unimanual or bimanual, and both tasks were scored as right-handed or left-handed. Most individuals did not exhibit significant hand preference for unimanual or bimanual (asymmetrical hand use) grooming, although 1 individual was lateralized for each. For the 8 subjects who engaged in termite fishing enough for statistical testing, 7 individuals exhibited significant laterality and strong individual hand preference. Four subjects preferred their left hand, 3 preferred their right, and 1 had no preference. Grooming, a simple behavior, was not lateralized in this group, yet a more complex behavior revealed a strong individual hand preference, and these results are congruent with other recent findings that demonstrate complex tasks elicit hand preference in bonobos.

Hair plucking, stress, and urinary cortisol among captive bonobos (Pan paniscus).

Hair plucking has been observed in many captive primate species, including the great apes; however, the etiology of this behavioral pattern is poorly understood. While this behavior has not been reported in wild apes, an ethologically identical behavior in humans, known as trichotillomania, is linked to chronic psychosocial stress and is a predominantly female disorder. This study examines hair plucking (defined here as a rapid jerking away of the hair shaft and follicle by the hand or mouth, often accompanied by inspection and consumption of the hair shaft and follicle) in a captive group of bonobos (N = 13) at the Columbus Zoo and Aquarium in Columbus, Ohio. Plucking data were collected using behavior and all-occurrence sampling; 1,450 social and self-directed grooming bouts were recorded during 128 hr of observation. Twenty-one percent of all grooming bouts involved at least one instance of plucking. Urine samples (N = 55) were collected and analyzed for the stress hormone cortisol. Analyses of urinary cortisol levels showed a significant positive correlation between mean cortisol and self-directed plucking for females (r = 0.88, P < 0.05) but not for males (r = -0.73, P = 0.09). These results demonstrate an association between relative self-directed hair plucking and cortisol among female bonobos. This is the first study to investigate the relationship between hair plucking and cortisol among apes. Overall, these data add to our knowledge of a contemporary issue in captive ape management. Zoo Biol. 35:415-422, 2016. © Wiley Periodicals, Inc.