Analysis of 100 high-coverage genomes from a pedigreed captive baboon colony

Comparative lifespan and healthspan of nonhuman primate species common to biomedical research

There is a critical need to generate age- and sex-specific survival curves to characterize chronological aging consistently across nonhuman primates (NHP) used in biomedical research. Accurate measures of chronological aging are essential for inferences into genetic, demographic, and physiological variables driving differences in NHP lifespan within and between species. Understanding NHP lifespans is relevant to public health because unraveling the demographic, molecular, and clinical bases of health across the life course in translationally relevant NHP species is fundamentally important to the study of human aging. Data from more than 110,000 captive individual NHP were contributed by 15 major research institutions to generate sex-specific Kaplan-Meier survival curves using uniform methods in 12 translational aging models: Callithrix jacchus (common marmoset), Chlorocebus aethiops sabaeus (vervet/African green), Macaca fascicularis (cynomolgus macaque), M. fuscata (Japanese macaque), M. mulatta (rhesus macaque), M. nemestrina (pigtail macaque), M. radiata (bonnet macaque), Pan troglodytes spp. (chimpanzee), Papio hamadryas spp. (baboon), Plecturocebus cupreus (coppery titi monkey), Saguinus oedipus (cotton-top tamarin), and Saimiri spp. (squirrel monkey). After employing strict inclusion criteria, primary analysis results are based on 12,269 NHP that survived to adulthood and died of natural/health-related causes. A secondary analysis was completed for 32,616 NHP that died of any cause. For the primary analyses, we report ages of 25th, 50th, 75th, and 85th percentiles of survival, maximum observed ages, rates of survivorship, and sex-based differences captured by quantile regression models and Kolmogorov-Smirnov tests. Our findings show a pattern of reduced male survival among catarrhines (African and Asian primates), especially macaques, but not platyrrhines (Central and South American primates). For many species, median lifespans were lower than previously reported. An important consideration is that these analyses may offer a better reflection of healthspan than lifespan. Captive NHP used in research are typically euthanized for humane welfare reasons before their natural end of life, often after diagnosis of their first major disease requiring long-term treatment with reduced quality of life (e.g., endometriosis, cancer, osteoarthritis). Supporting the idea that these data are capturing healthspan, for several species typical age at onset of chronic disease is similar to the median lifespan estimates. This data resource represents the most comprehensive characterization of sex-specific lifespan and age-at-death distributions for 12 biomedically relevant species, to date. The results clarify the relationships among NHP ages and will provide a valuable resource for the aging research community, improving human-NHP age equivalencies, informing investigators of the expected survival rates of NHP assigned to studies, providing a metric for comparisons in future studies, and contributing to our understanding of the factors that drive lifespan differences within and among species.

Global and local ancestry estimation in a captive baboon colony

The last couple of decades have highlighted the importance of studying hybridization, particularly among primate species, as it allows us to better understand our own evolutionary trajectory. Here, we report on genetic ancestry estimates using dense, full genome data from 881 olive (Papio anubus), yellow (Papio cynocephalus), or olive-yellow crossed captive baboons from the Southwest National Primate Research Center. We calculated global and local ancestry information, imputed low coverage genomes (n = 830) to improve marker quality, and updated the genetic resources of baboons available to assist future studies. We found evidence of historical admixture in some putatively purebred animals and identified errors within the Southwest National Primate Research Center pedigree. We also compared the outputs between two different phasing and imputation pipelines along with two different global ancestry estimation software. There was good agreement between the global ancestry estimation software, with R2 > 0.88, while evidence of phase switch errors increased depending on what phasing and imputation pipeline was used. We also generated updated genetic maps and created a concise set of ancestry informative markers (n = 1,747) to accurately obtain global ancestry estimates.

Genetic regulatory effects in response to a high-cholesterol, high-fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWASs), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using non-human primate models.

Genetic regulatory effects in response to a high cholesterol, high fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWAS), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using nonhuman primate models.

Landscape genomics analysis provides insights into future climate change-driven risk in rhesus macaque

Climate change significantly affects the suitability of wildlife habitats. Thus, understanding how animals adapt ecologically and genetically to climate change is important for targeted species protection. Rhesus macaques (Macaca mulatta) are widely distributed and multi-climatically adapted primates. This study explored how rhesus macaques adapt to climate change by integrating ecological and genetic methods and applying species distribution models (SDMs) and a gradient forest (GF) model. The findings suggested that temperature seasonality primarily affects habitat suitability and indicated that climate change will have a dramatic impact on macaque populations in the future. We also applied genotype-environment association (GEA) analyses and selection signature analyses to identify genes associated with climate change and provide possible explanations for the adaptation of rhesus macaques to climatic environments. The population genomics analyses suggested that the Taihang population has the highest genomic vulnerability with inbreeding and low heterozygosity. Combined with the higher ecological vulnerability, additional conservation strategies are required for this population under higher risk of climate change. Our work measured the impact of climate change and enabled the identification of populations that exhibit high vulnerability to severe climate change. Such information is useful for selecting populations of rhesus macaques as subject of long-term monitoring or evolutionary rescue under future climate change.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Combinatorial targeting of epigenome-modifying enzymes with decitabine and RN-1 synergistically increases HbF

Increased fetal hemoglobin (HbF) levels reduce the symptoms of sickle cell disease (SCD) and increase the lifespan of patients. Because curative strategies for bone marrow transplantation and gene therapy technologies remain unavailable to a large number of patients, the development of a safe and effective pharmacological therapy that increases HbF offers the greatest potential for disease intervention. Although hydroxyurea increases HbF, a substantial proportion of patients fail to demonstrate an adequate response. Pharmacological inhibitors of DNA methyltransferase (DNMT1) and lysine-specific demethylase 1A (LSD1), 2 epigenome-modifying enzymes associated with the multiprotein corepressor complex recruited to the repressed γ-globin gene, are powerful in vivo inducers of HbF. The hematological side effects of these inhibitors limit feasible clinical exposures. We evaluated whether administering these drugs in combination could reduce the dose and/or time of exposure to any single agent to minimize adverse effects, while achieving additive or synergistic increases in HbF. The DNMT1 inhibitor decitabine (0.5 mg/kg per day) and the LSD1 inhibitor RN-1 (0.25 mg/kg per day) administered in combination 2 days per week produced synergistic increases in F-cells, F-reticulocytes, and γ-globin messenger RNA in healthy baboons. Large increases in HbF and F-cells were observed in healthy, nonanemic, and anemic (phlebotomized) baboons. Combinatorial therapy targeting epigenome-modifying enzymes could thus be a useful strategy for producing larger increases in HbF to modify the clinical course of SCD.

Age- and sex-related differences in baboon (Papio anubis) gray matter covariation

Age-related changes in cognition, brain morphology, and behavior are exhibited in several primate species. Baboons, like humans, naturally develop Alzheimer's disease-like pathology and cognitive declines with age and are an underutilized model for studies of aging. To determine age-related differences in gray matter covariation of 89 olive baboons (Papio anubis), we used source-based morphometry (SBM) to analyze data from magnetic resonance images. We hypothesized that we would find significant age effects in one or more SBM components, particularly those which include regions influenced by age in humans and other nonhuman primates (NHPs). A multivariate analysis of variance revealed that individual weighted gray matter covariation scores differed across the age classes. Elderly baboons contributed significantly less to gray matter covariation components including the brainstem, superior parietal cortex, thalamus, and pallidum compared to juveniles, and middle and superior frontal cortex compared to juveniles and young adults (p < 0.05). Future studies should examine the relationship between the changes in gray matter covariation reported here and age-related cognitive decline.

Ecology and age, but not genetic ancestry, predict fetal loss in a wild baboon hybrid zone

OBJECTIVES

Pregnancy failure represents a major fitness cost for any mammal, particularly those with slow life histories such as primates. Here, we quantified the risk of fetal loss in wild hybrid baboons, including genetic, ecological, and demographic sources of variance. We were particularly interested in testing the hypothesis that hybridization increases fetal loss rates. Such an effect would help explain how baboons may maintain genetic and phenotypic integrity despite interspecific gene flow.

MATERIALS AND METHODS

We analyzed outcomes for 1020 pregnancies observed over 46 years in a natural yellow baboon-anubis baboon hybrid zone. Fetal losses and live births were scored based on records of female reproductive state and the appearance of live neonates. We modeled the probability of fetal loss as a function of a female's genetic ancestry (the proportion of her genome estimated to be descended from anubis [vs. yellow] ancestors), age, number of previous fetal losses, dominance rank, group size, climate, and habitat quality using binomial mixed effects models.

RESULTS

Female genetic ancestry did not predict fetal loss. Instead, the risk of fetal loss is elevated for very young and very old females. Fetal loss is most robustly predicted by ecological factors, including poor habitat quality prior to a home range shift and extreme heat during pregnancy.

DISCUSSION

Our results suggest that gene flow between yellow and anubis baboons is not impeded by an increased risk of fetal loss for hybrid females. Instead, ecological conditions and female age are key determinants of this component of female reproductive success.

Selection against admixture and gene regulatory divergence in a long-term primate field study

Genetic admixture is central to primate evolution. We combined 50 years of field observations of immigration and group demography with genomic data from ~9 generations of hybrid baboons to investigate the consequences of admixture in the wild. Despite no obvious fitness costs to hybrids, we found signatures of selection against admixture similar to those described for archaic hominins. These patterns were concentrated near genes where ancestry is strongly associated with gene expression. Our analyses also show that introgression is partially predictable across the genome. This study demonstrates the value of integrating genomic and field data for revealing how "genomic signatures of selection" (e.g., reduced introgression in low-recombination regions) manifest in nature; moreover, it underscores the importance of other primates as living models for human evolution.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species

Background

Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes.

Results

We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the “northern” Papio clade, and signal the presence of population structure within P. ursinus.

Conclusions

The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01999-7.

Mechanisms of inbreeding avoidance in a wild primate

Inbreeding often imposes net fitness costs,1-5 leading to the expectation that animals will engage in inbreeding avoidance when the costs of doing so are not prohibitive.4-9 However, one recent meta-analysis indicates that animals of many species do not avoid mating with kin in experimental settings,6 and another reports that behavioral inbreeding avoidance generally evolves only when kin regularly encounter each other and inbreeding costs are high.9 These results raise questions about the processes that separate kin, how these processes depend on kin class and context, and whether kin classes differ in how effectively they avoid inbreeding via mate choice-in turn, demanding detailed demographic and behavioral data within individual populations. Here, we address these questions in a wild mammal population, the baboons of the Amboseli ecosystem in Kenya. We find that death and dispersal are very effective at separating opposite-sex pairs of close adult kin. Nonetheless, adult kin pairs do sometimes co-reside, and we find strong evidence for inbreeding avoidance via mate choice in kin classes with relatedness ≥0.25. Notably, maternal kin avoid inbreeding more effectively than paternal kin despite having identical coefficients of relatedness, pointing to kin discrimination as a potential constraint on effective inbreeding avoidance. Overall, demographic and behavioral processes ensure that inbred offspring are rare in undisturbed social groups (1% of offspring). However, in an anthropogenically disturbed social group with reduced male dispersal, we find inbreeding rates 10× higher. Our study reinforces the importance of demographic and behavioral contexts for understanding the evolution of inbreeding avoidance.9.

High-Resolution Estimates of Crossover and Noncrossover Recombination from a Captive Baboon Colony

Homologous recombination has been extensively studied in humans and a handful of model organisms. Much less is known about recombination in other species, including nonhuman primates. Here, we present a study of crossovers (COs) and noncrossover (NCO) recombination in olive baboons (Papio anubis) from two pedigrees containing a total of 20 paternal and 17 maternal meioses, and compare these results to linkage disequilibrium (LD) based recombination estimates from 36 unrelated olive baboons. We demonstrate how COs, combined with LD-based recombination estimates, can be used to identify genome assembly errors. We also quantify sex-specific differences in recombination rates, including elevated male CO and reduced female CO rates near telomeres. Finally, we add to the increasing body of evidence suggesting that while most NCO recombination tracts in mammals are short (e.g., <500 bp), there is a non-negligible fraction of longer (e.g., >1 kb) NCO tracts. For NCO tracts shorter than 10 kb, we fit a mixture of two (truncated) geometric distributions model to the NCO tract length distribution and estimate that >99% of all NCO tracts are very short (mean 24 bp), but the remaining tracts can be quite long (mean 4.3 kb). A single geometric distribution model for NCO tract lengths is incompatible with the data, suggesting that LD-based methods for estimating NCO recombination rates that make this assumption may need to be modified.

Genetic ancestry predicts male-female affiliation in a natural baboon hybrid zone

Opposite-sex social relationships are important predictors of fitness in many animals, including several group-living mammals. Consequently, understanding sources of variance in the tendency to form opposite-sex relationships is important for understanding social evolution. Genetic contributions are of particular interest due to their importance in long-term evolutionary change, but little is known about genetic effects on male-female relationships in social mammals, especially outside of the mating context. Here, we investigate the effects of genetic ancestry on male-female affiliative behaviour in a hybrid zone between the yellow baboon, Papio cynocephalus, and the anubis baboon, Papio anubis, in a population in which male-female social bonds are known predictors of life span. We place our analysis within the context of other social and demographic predictors of affiliative behaviour in baboons. Genetic ancestry was the most consistent predictor of opposite-sex affiliative behaviour we observed, with the exception of strong effects of dominance rank. Our results show that increased anubis genetic ancestry is associated with a subtle, but significantly higher, probability of opposite-sex affiliative behaviour, in both males and females. Additionally, pairs of anubis-like males and anubis-like females were the most likely to socially affiliate, resulting in moderate assortativity in grooming and proximity behaviour as a function of genetic ancestry. Our findings indicate that opposite-sex affiliative behaviour partially diverged during baboon evolution to differentiate yellow and anubis baboons, despite overall similarities in their social structures and mating systems. Furthermore, they suggest that affiliative behaviour may simultaneously promote and constrain baboon admixture, through additive and assortative effects of ancestry, respectively.

A Fine-Scale Genetic Map for Vervet Monkeys

Despite its important biological role, the evolution of recombination rates remains relatively poorly characterized. This owes, in part, to the lack of high-quality genomic resources to address this question across diverse species. Humans and our closest evolutionary relatives, anthropoid apes, have remained a major focus of large-scale sequencing efforts, and thus recombination rate variation has been comparatively well studied in this group-with earlier work revealing a conservation at the broad- but not the fine-scale. However, in order to better understand the nature of this variation, and the time scales on which substantial modifications occur, it is necessary to take a broader phylogenetic perspective. I here present the first fine-scale genetic map for vervet monkeys based on whole-genome population genetic data from ten individuals and perform a series of comparative analyses with the great apes. The results reveal a number of striking features. First, owing to strong positive correlations with diversity and weak negative correlations with divergence, analyses suggest a dominant role for purifying and background selection in shaping patterns of variation in this species. Second, results support a generally reduced broad-scale recombination rate compared with the great apes, as well as a narrower fraction of the genome in which the majority of recombination events are observed to occur. Taken together, this data set highlights the great necessity of future research to identify genomic features and quantify evolutionary processes that are driving these rate changes across primates.

Genome-wide ancestry and introgression in a Zambian baboon hybrid zone

Hybridization in nature offers unique insights into the process of natural selection in incipient species and their hybrids. In order to evaluate the patterns and targets of selection, we examine a recently discovered baboon hybrid zone in the Kafue River Valley of Zambia, where Kinda baboons (Papio kindae) and grey-footed chacma baboons (P. ursinus griseipes) coexist with hybridization. We genotyped baboons at 14,962 variable genome-wide autosomal markers using double-digest RADseq. We compared ancestry patterns from this genome-wide data set to previously reported ancestry from mitochondrial-DNA and Y-chromosome sources. We also fit a Bayesian genomic cline model to scan for genes with extreme patterns of introgression. We show that the Kinda baboon Y chromosome has penetrated the species boundary to a greater extent than either mitochondrial DNA or the autosomal chromosomes. We also find evidence for overall restricted introgression in the JAK/STAT signalling pathway. Echoing results in other species including humans, we find evidence for enhanced and/or directional introgression of immune-related genes or pathways including the toll-like receptor pathway, the blood coagulation pathway, and the LY96 gene. Finally we show enhanced introgression and excess chacma baboon ancestry in the sperm tail gene ODF2. Together, our results elucidate the dynamics of introgressive hybridization in a primate system while identifying genes and pathways possibly under selection.

Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data

BACKGROUND

Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses.

FINDINGS

Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome.

CONCLUSIONS

We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.

A comparison of humans and baboons suggests germline mutation rates do not track cell divisions

In humans, most germline mutations are inherited from the father. This observation has been widely interpreted as reflecting the replication errors that accrue during spermatogenesis. If so, the male bias in mutation should be substantially lower in a closely related species with similar rates of spermatogonial stem cell divisions but a shorter mean age of reproduction. To test this hypothesis, we resequenced two 3-4 generation nuclear families (totaling 29 individuals) of olive baboons (Papio anubis), who reproduce at approximately 10 years of age on average, and analyzed the data in parallel with three 3-generation human pedigrees (26 individuals). We estimated a mutation rate per generation in baboons of 0.57×10-8 per base pair, approximately half that of humans. Strikingly, however, the degree of male bias in germline mutations is approximately 4:1, similar to that of humans-indeed, a similar male bias is seen across mammals that reproduce months, years, or decades after birth. These results mirror the finding in humans that the male mutation bias is stable with parental ages and cast further doubt on the assumption that germline mutations track cell divisions. Our mutation rate estimates for baboons raise a further puzzle, suggesting a divergence time between apes and Old World monkeys of 65 million years, too old to be consistent with the fossil record; reconciling them now requires not only a slowdown of the mutation rate per generation in humans but also in baboons.

Insights into the evolution of social systems and species from baboon studies

Baboons, members of the genus Papio, comprise six closely related species distributed throughout sub-Saharan Africa and southwest Arabia. The species exhibit more ecological flexibility and a wider range of social systems than many other primates. This article summarizes our current knowledge of the natural history of baboons and highlights directions for future research. We suggest that baboons can serve as a valuable model for complex evolutionary processes, such as speciation and hybridization. The evolution of baboons has been heavily shaped by climatic changes and population expansion and fragmentation in the African savanna environment, similar to the processes that acted during human evolution. With accumulating long-term data, and new data from previously understudied species, baboons are ideally suited for investigating the links between sociality, health, longevity and reproductive success. To achieve these aims, we propose a closer integration of studies at the proximate level, including functional genomics, with behavioral and ecological studies.