Genetic diversity of 1,845 rhesus macaques improves genetic variation interpretation and identifies disease models.

Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

Immune gene regulation is associated with age and environmental adversity in a nonhuman primate.

Phenotypic aging is ubiquitous across mammalian species, suggesting shared underlying mechanisms of aging. Aging is linked to molecular changes to DNA methylation and gene expression, and environmental factors, such as severe external challenges or adversities, can moderate these age-related changes. Yet, it remains unclear whether environmental adversities affect gene regulation via the same molecular pathways as chronological, or 'primary', aging. Investigating molecular aging in naturalistic animal populations can fill this gap by providing insight into shared molecular mechanisms of aging and the effects of a greater diversity of environmental adversities - particularly those that can be challenging to study in humans or laboratory organisms. Here, we characterised molecular aging - specifically, CpG methylation - in a sample of free-ranging rhesus macaques living off the coast of Puerto Rico (n samples = 571, n individuals = 499), which endured a major hurricane during our study. Age was associated with methylation at 78,661 sites (31% of all sites tested). Age-associated hypermethylation occurred more frequently in areas of active gene regulation, while hypomethylation was enriched in regions that show less activity in immune cells, suggesting these regions may become de-repressed in older individuals. Age-associated hypomethylation also co-occurred with increased chromatin accessibility while hypermethylation showed the opposite trend, hinting at a coordinated, multi-level loss of epigenetic stability during aging. We detected 32,048 CpG sites significantly associated with exposure to a hurricane, and these sites overlapped age-associated sites, most strongly in regulatory regions and most weakly in quiescent regions. Together, our results suggest that environmental adversity may contribute to aging-related molecular phenotypes in regions of active gene transcription, but that primary aging has specific signatures in non-regulatory regions.

Evolutionary and biomedical implications of sex differences in the primate brain transcriptome.

Humans exhibit sex differences in the prevalence of many neurodevelopmental disorders and neurodegenerative diseases. Here, we generated one of the largest multi-brain-region bulk transcriptional datasets for the rhesus macaque and characterized sex-biased gene expression patterns to investigate the translatability of this species for sex-biased neurological conditions. We identify patterns similar to those in humans, which are associated with overlapping regulatory mechanisms, biological processes, and genes implicated in sex-biased human disorders, including autism. We also show that sex-biased genes exhibit greater genetic variance for expression and more tissue-specific expression patterns, which may facilitate rapid evolution of sex-biased genes. Our findings provide insights into the biological mechanisms underlying sex-biased disease and support the rhesus macaque model for the translational study of these conditions.

Quantifying the relationship between bone and soft tissue measures within the rhesus macaques of Cayo Santiago.

OBJECTIVES: Interpretations of the primate and human fossil record often rely on the estimation of somatic dimensions from bony measures. Both somatic and skeletal variation have been used to assess how primates respond to environmental change. However, it is unclear how well skeletal variation matches and predicts soft tissue. Here, we empirically test the relationship between tissues by comparing somatic and skeletal measures using paired measures of pre- and post-mortem rhesus macaques from Cayo Santiago, Puerto Rico. MATERIALS AND METHODS: Somatic measurements were matched with skeletal dimensions from 105 rhesus macaque individuals to investigate paired signals of variation (i.e., coefficients of variation, sexual dimorphism) and bivariate codependence (reduced major axis regression) in measures of: (1) limb length; (2) joint breadth; and (3) limb circumference. Predictive models for the estimation of soft tissue dimensions from skeletons were built from Ordinary Least Squares regressions. RESULTS: Somatic and skeletal measurements showed statistically equivalent coefficients of variation and sexual dimorphism as well as high epiphyses-present ordinary least square (OLS) correlations in limb lengths (R2 >0.78, 0.82), joint breadths (R2 >0.74, 0.83) and, to a lesser extent, limb circumference (R2 >0.53, 0.68). CONCLUSION: Skeletal measurements are good substitutions for somatic values based on population signals of variation. OLS regressions indicate that skeletal correlates are highly predictive of somatic dimensions. The protocols and regression equations established here provide a basis for reliable reconstruction of somatic dimension from catarrhine fossils and validate our ability to compare or combine results of studies based on population data of either hard or soft tissue proxies.

Mechanical and morphometric approaches to body mass estimation in rhesus macaques: A test of skeletal variables.

OBJECTIVES: Estimation of body mass from skeletal metrics can reveal important insights into the paleobiology of archeological or fossil remains. The standard approach constructs predictive equations from postcrania, but studies have questioned the reliability of traditional measures. Here, we examine several skeletal features to assess their accuracy in predicting body mass. MATERIALS AND METHODS: Antemortem mass measurements were compared with common skeletal dimensions from the same animals postmortem, using 115 rhesus macaques (male: n = 43; female: n = 72). Individuals were divided into training (n = 58) and test samples (n = 57) to build and assess Ordinary Least Squares or multivariate regressions by residual sum of squares (RSS) and AIC weights. A leave-one-out approach was implemented to formulate the best fit multivariate models, which were compared against a univariate and a previously published catarrhine body-mass estimation model. RESULTS: Femur circumference represented the best univariate model. The best model overall was composed of four variables (femur, tibia and fibula circumference and humerus length). By RSS and AICw, models built from rhesus macaque data (RSS = 26.91, AIC = -20.66) better predicted body mass than did the catarrhine model (RSS = 65.47, AIC = 20.24). CONCLUSION: Body mass in rhesus macaques is best predicted by a 4-variable equation composed of humerus length and hind limb midshaft circumferences. Comparison of models built from the macaque versus the catarrhine data highlight the importance of taxonomic specificity in predicting body mass. This paper provides a valuable dataset of combined somatic and skeletal data in a primate, which can be used to build body mass equations for fragmentary fossil evidence.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques.

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.