Search

Genetic variation that impacts gene regulation, rather than protein function, can have strong effects on trait variation both within and between species. Epigenetic mechanisms, such as DNA methylation, are often an important intermediate link between genotype and phenotype, yet genetic effects on DNA methylation remain understudied in natural populations. To address this gap, we used reduced representation bisulfite sequencing to measure DNA methylation levels at 555,856 CpGs in peripheral whole blood of 573 samples collected from free-ranging rhesus macaques (Macaca mulatta) living on the island of Cayo Santiago, Puerto Rico. We used allele-specific methods to map cis-methylation quantitative trait loci (meQTL) and tested for effects of 243,389 single nucleotide polymorphisms (SNPs) on local DNA methylation levels. Of 776,092 tested SNP-CpG pairs, we identified 516,213 meQTL, with 69.12% of CpGs having at least one meQTL (FDR < 5%). On average, meQTL explained 21.2% of nearby methylation variance, significantly more than age or sex. meQTL were enriched in genomic compartments where methylation is likely to impact gene expression, for example, promoters, enhancers and binding sites for methylation-sensitive transcription factors. In support, using mRNA-seq data from 172 samples, we confirmed 332 meQTL as whole blood cis-expression QTL (eQTL) in the population, and found meQTL-eQTL genes were enriched for immune response functions, like antigen presentation and inflammation. Overall, our study takes an important step towards understanding the genetic architecture of DNA methylation in natural populations, and more generally points to the biological mechanisms driving phenotypic variation in our close relatives., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

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., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

The origin of primates has long been associated with an increased emphasis on manual grasping and touch. Precision touch, facilitated by specialized mechanoreceptors in glabrous skin, provides critical sensory feedback for grasping-related tasks and perception of ecologically-relevant stimuli. Despite its importance, studies of mechanoreceptors in primate hands are limited, in part due to challenges of sample availability and histological methods. Dermatoglyphs have been proposed as alternative proxies of mechanoreceptor density. We investigated the relationships between mechanoreceptors (Meissner and Pacinian corpuscles), dermatoglyphs, and demography in the apical finger pads of 15 juvenile to adult rhesus macaques (Macaca mulatta) from a free-ranging population at Cayo Santiago Primate Field Station (Puerto Rico). Our results indicate substantial interindividual variation in mechanoreceptor density (Meissner corpuscles: 11.9-43.3 corpuscles/mm 2 ; Pacinian corpuscles: 0-4.5 corpuscles/mm 2 ). While sex and digit were generally not associated with variation, there was strong evidence of a developmental effect. Specifically, apical pad length, Meissner corpuscle size, and Pacinian corpuscle depth increased while mechanoreceptor densities decreased throughout juvenescence, suggesting that primate mechanoreceptors change as fingers grow during adolescence and then stabilize at physical maturity. We also found Meissner corpuscle density was significantly associated with dermatoglyph ridge width and spacing, such that density predicted by a dermatoglyph model was strongly correlated with observed values. Dermatoglyphs thus offer a useful proxy of relative Meissner corpuscle density in primates, which opens exciting avenues of noninvasive research. Finally, our results underscore the importance of considering demographic factors and methodology in comparative studies of primate touch., (© 2024 American Association for Anatomy.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2024. The Author(s).)

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 (R 2  >0.78, 0.82), joint breadths (R 2  >0.74, 0.83) and, to a lesser extent, limb circumference (R 2  >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., (© 2024 Wiley Periodicals LLC.)

Extreme weather events radically alter ecosystems. When ecological damage persists, selective pressures on individuals can change, leading to phenotypic adjustments. For group-living animals, social relationships may be a mechanism enabling adaptation to ecosystem disturbance. Yet whether such events alter selection on sociality and whether group-living animals can, as a result, adaptively change their social relationships remain untested. We leveraged 10 years of data collected on rhesus macaques before and after a category 4 hurricane caused persistent deforestation, exacerbating monkeys' exposure to intense heat. In response, macaques demonstrated persistently increased tolerance and decreased aggression toward other monkeys, facilitating access to scarce shade critical for thermoregulation. Social tolerance predicted individual survival after the hurricane, but not before it, revealing a shift in the adaptive function of sociality.

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., (© 2024 Wiley Periodicals LLC.)

Increasing age is associated with dysregulated immune function and increased inflammation-patterns that are also observed in individuals exposed to chronic social adversity. Yet we still know little about how social adversity impacts the immune system and how it might promote age-related diseases. Here, we investigated how immune cell diversity varied with age, sex and social adversity (operationalized as low social status) in free-ranging rhesus macaques. We found age-related signatures of immunosenescence, including lower proportions of CD20 + B cells, CD20 + /CD3 + ratio, and CD4 + /CD8 + T cell ratio - all signs of diminished antibody production. Age was associated with higher proportions of CD3 + /CD8 + Cytotoxic T cells, CD16 + /CD3- Natural Killer cells, CD3 + /CD4 + /CD25 + and CD3 + /CD8 + /CD25 + T cells, and CD14 + /CD16 + /HLA-DR + intermediate monocytes, and lower levels of CD14 + /CD16-/HLA-DR + classical monocytes, indicating greater amounts of inflammation and immune dysregulation. We also found a sex-dependent effect of exposure to social adversity (i.e., low social status). High-status males, relative to females, had higher CD20 + /CD3 + ratios and CD16 + /CD3 Natural Killer cell proportions, and lower proportions of CD8 + Cytotoxic T cells. Further, low-status females had higher proportions of cytotoxic T cells than high-status females, while the opposite was observed in males. High-status males had higher CD20 + /CD3 + ratios than low-status males. Together, our study identifies the strong age and sex-dependent effects of social adversity on immune cell proportions in a human-relevant primate model. Thus, these results provide novel insights into the combined effects of demography and social adversity on immunity and their potential contribution to age-related diseases in humans and other animals., (© 2023. The Author(s), under exclusive licence to American Aging Association.)

Social adversity can increase the age-associated risk of disease and death, yet the biological mechanisms that link social adversities to aging remain poorly understood. Long-term naturalistic studies of nonhuman animals are crucial for integrating observations of social behavior throughout an individual's life with detailed anatomical, physiological, and molecular measurements. Here, we synthesize the body of research from one such naturalistic study system, Cayo Santiago, which is home to the world's longest continuously monitored free-ranging population of rhesus macaques (Macaca mulatta). We review recent studies of age-related variation in morphology, gene regulation, microbiome composition, and immune function. We also discuss ecological and social modifiers of age-markers in this population. In particular, we summarize how a major natural disaster, Hurricane Maria, affected rhesus macaque physiology and social structure and highlight the context-dependent and domain-specific nature of aging modifiers. Finally, we conclude by providing directions for future study, on Cayo Santiago and elsewhere, that will further our understanding of aging across different domains and how social adversity modifies aging processes., Competing Interests: Competing interests The authors declare that they have no competing interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Cataloging the diverse cellular architecture of the primate brain is crucial for understanding cognition, behavior, and disease in humans. Here, we generated a brain-wide single-cell multimodal molecular atlas of the rhesus macaque brain. Together, we profiled 2.58 M transcriptomes and 1.59 M epigenomes from single nuclei sampled from 30 regions across the adult brain. Cell composition differed extensively across the brain, revealing cellular signatures of region-specific functions. We also identified 1.19 M candidate regulatory elements, many previously unidentified, allowing us to explore the landscape of cis-regulatory grammar and neurological disease risk in a cell type-specific manner. Altogether, this multi-omic atlas provides an open resource for investigating the evolution of the human brain and identifying novel targets for disease interventions.

Nonhuman primates (NHPs) are vital translational research models due to their high genetic, physiological, and anatomical homology with humans. The "golden" rhesus macaque (Macaca mulatta) phenotype is a naturally occurring, inherited trait with a visually distinct pigmentation pattern resulting in light blonde colored fur. Retinal imaging also reveals consistent hypopigmentation and occasional foveal hypoplasia. Here, we describe the use of genome-wide association in 2 distinct NHP populations to identify candidate variants in genes linked to the golden phenotype. Two missense variants were identified in the Tyrosinase-related protein 1 gene (Asp343Gly and Leu415Pro) that segregate with the phenotype. An additional and distinct association was also found with a Tyrosinase variant (His256Gln), indicating the light-colored fur phenotype can result from multiple genetic mechanisms. The implicated genes are related through their contribution to the melanogenesis pathway. Variants in these 2 genes are known to cause pigmentation phenotypes in other species and to be associated with oculocutaneous albinism in humans. The novel associations presented in this study will permit further investigations into the role these proteins and variants play in the melanogenesis pathway and model the effects of genetic hypopigmentation and altered melanogenesis in a naturally occurring nonhuman primate model., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

While skin microbes are known to mediate human health and disease, there has been minimal research on the interactions between skin microbiota, social behavior, and year-to-year effects in non-human primates-important animal models for translational biomedical research. To examine these relationships, we analyzed skin microbes from 78 rhesus macaques living on Cayo Santiago Island, Puerto Rico. We considered age, sex, and social group membership, and characterized social behavior by assessing dominance rank and patterns of grooming as compared to nonsocial behaviors. To measure the effects of a shifting environment, we sampled skin microbiota (based on sequence analysis of the 16S rRNA V4 region) and assessed weather across sampling periods between 2013 and 2015. We hypothesized that, first, monkeys with similar social behavior and/or in the same social group would possess similar skin microbial composition due, in part, to physical contact, and, second, microbial diversity would differ across sampling periods. We found significant phylum-level differences between social groups in the core microbiome as well as an association between total grooming rates and alpha diversity in the complete microbiome, but no association between microbial diversity and measures of rank or other nonsocial behaviors. We also identified alpha and beta diversity differences in microbiota and differential taxa abundance across two sampling periods. Our findings indicate that social dynamics interact with yearly environmental changes to shape the skin microbiota in rhesus macaques, with potential implications for understanding the factors affecting the microbiome in humans, which share many biological and social characteristics with these animals. IMPORTANCE Primate studies are valuable for translational and evolutionary insights into the human microbiome. The majority of primate microbiome studies focus on the gut, so less is known about the factors impacting the microbes on skin and how their links affect health and behavior. Here, we probe the impact of social interactions and the yearly environmental changes on food-provisioned, free-ranging monkeys living on a small island. We expected animals that lived together and groomed each other would have more similar microbes on their skin, but surprisingly found that the external environment was a stronger influence on skin microbiome composition. These findings have implications for our understanding of the human skin microbiome, including potential manipulations to improve health and treat disease.

Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species-specific genotyping methods (e.g. array-based genotyping) and must instead use sequencing-based approaches. Although costs are decreasing, high-coverage whole-genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low-coverage WGS combined with imputation using a high-confidence reference panel is a cost-effective alternative, but the accuracy of genotyping using low-coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low-coverage sequencing (0.1-10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (median r 2  ≥ 0.92). We show that low-coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI-for-non-model-species). Our results endorse accurate and effective genotype imputation from low-coverage sequencing, enabling the cost-effective generation of population-scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation., (© 2023 John Wiley & Sons Ltd.)

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Aging is accompanied by a host of social and biological changes that correlate with behavior, cognitive health and susceptibility to neurodegenerative disease. To understand trajectories of brain aging in a primate, we generated a multiregion bulk (N = 527 samples) and single-nucleus (N = 24 samples) brain transcriptional dataset encompassing 15 brain regions and both sexes in a unique population of free-ranging, behaviorally phenotyped rhesus macaques. We demonstrate that age-related changes in the level and variance of gene expression occur in genes associated with neural functions and neurological diseases, including Alzheimer's disease. Further, we show that higher social status in females is associated with younger relative transcriptional ages, providing a link between the social environment and aging in the brain. Our findings lend insight into biological mechanisms underlying brain aging in a nonhuman primate model of human behavior, cognition and health., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

A defining feature of catarrhine primates is uniform trichromacy-the ability to distinguish red (long; L), green (medium; M), and blue (short; S) wavelengths of light. Although the tuning of photoreceptors is conserved, the ratio of L:M cones in the retina is variable within and between species, with human cone ratios differing from other catarrhines. Yet, the sources and structure of variation in cone ratios are poorly understood, precluding a broader understanding of color vision variability. Here, we report a large-scale study of a pedigreed population of rhesus macaques (Macaca mulatta). We collected foveal RNA and analyzed opsin gene expression using cDNA and estimated additive genetic variance of cone ratios. The average L:M ratio and standard error was 1.03:1 ± 0.02. There was no age effect, and genetic contribution to variation was negligible. We found marginal sex effects with females having larger ratios than males. S cone ratios (0.143:1 ± 0.002) had significant genetic variance with a heritability estimate of 43% but did not differ between sexes or age groups. Our results contextualize the derived human condition of L-cone dominance and provide new information about the heritability of cone ratios and variation in primate color vision., (© 2022 The Authors. Evolution published by Wiley Periodicals LLC on behalf of The Society for the Study of Evolution.)

Reproduction and survival in most primate species reflects management of both competitive and cooperative relationships. Here, we investigated the links between neuroanatomy and sociality in free-ranging rhesus macaques. In adults, the number of social partners predicted the volume of the mid-superior temporal sulcus and ventral-dysgranular insula, implicated in social decision-making and empathy, respectively. We found no link between brain structure and other key social variables such as social status or indirect connectedness in adults, nor between maternal social networks or status and dependent infant brain structure. Our findings demonstrate that the size of specific brain structures varies with the number of direct affiliative social connections and suggest that this relationship may arise during development. These results reinforce proposed links between social network size, biological success, and the expansion of specific brain circuits.

Weather-related disasters are increasing in frequency and severity, leaving survivors to cope with ensuing mental, financial, and physical hardships. This adversity can exacerbate existing morbidities, trigger new ones, and increase the risk of mortality-features that are also characteristic of advanced age-inviting the hypothesis that extreme weather events may accelerate aging. To test this idea, we examined the impact of Hurricane Maria and its aftermath on immune cell gene expression in large, age-matched, cross-sectional samples from free-ranging rhesus macaques ( Macaca mulatta ) living on an isolated island. A cross section of macaques was sampled 1 to 4 y before ( n = 435) and 1 y after ( n = 108) the hurricane. Hurricane Maria was significantly associated with differential expression of 4% of immune-cell-expressed genes, and these effects were correlated with age-associated alterations in gene expression. We further found that individuals exposed to the hurricane had a gene expression profile that was, on average, 1.96 y older than individuals that were not-roughly equivalent to an increase in 7 to 8 y of a human life. Living through an intense hurricane and its aftermath was associated with expression of key immune genes, dysregulated proteostasis networks, and greater expression of inflammatory immune cell-specific marker genes. Together, our findings illuminate potential mechanisms through which the adversity unleashed by extreme weather and potentially other natural disasters might become biologically embedded, accelerate age-related molecular immune phenotypes, and ultimately contribute to earlier onset of disease and death., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

Objective: Reconstructing the social lives of extinct primates is possible only through an understanding of the interplay between morphology, sexual selection pressures, and social behavior in extant species. Somatic sexual dimorphism is an important variable in primate evolution, in part because of the clear relationship between the strength and mechanisms of sexual selection and the degree of dimorphism. Here, we examine body size dimorphism across ontogeny in male and female rhesus macaques to assess whether it is primarily achieved via bimaturism as predicted by a polygynandrous mating system, faster male growth indicating polygyny, or both., Methods: We measured body mass in a cross-sectional sample of 364 free-ranging rhesus macaques from Cayo Santiago, Puerto Rico to investigate size dimorphism: 1) across the lifespan; and 2) as an outcome of sex-specific growth strategies, including: a) age of maturation; b) growth rate; and c) total growth duration, using regression models fit to sex-specific developmental curves., Results: Significant body size dimorphism was observed by prime reproductive age with males 1.51 times the size of females. Larger male size resulted from a later age of maturation (males: 6.8-7.8 years versus females: 5.5-6.5 years; logistic model) and elevated growth velocity through the pre-prime period (LOESS model). Though males grew to larger sizes overall, females maintained adult size for longer before senescence (quadratic model)., Discussion: The ontogeny of size dimorphism in rhesus macaques is achieved by bimaturism and a faster male growth rate. Our results provide new data for understanding the development and complexities of primate dimorphism., Competing Interests: Conflict of Interest The authors declare that there is no conflict of interest regarding this publication.