Managing the Cayo Santiago rhesus macaque population: The role of density

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.

Differential effects of early life adversity on male and female rhesus macaque lifespan

Early life adversity predicts shorter adult lifespan in several animal taxa. Yet, work on long-lived primate populations suggests the evolution of mechanisms that contribute to resiliency and long lives despite early life insults. Here, we tested associations between individual and cumulative early life adversity and lifespan on rhesus macaques at the Cayo Santiago Biological Field Station using 50 years of demographic data. We performed sex-specific survival analyses at different life stages to contrast short-term effects of adversity (i.e., infant survival) with long-term effects (i.e., adult survival). Female infants showed vulnerability to multiple adversities at birth, but affected females who survived to adulthood experienced a reduced risk later in life. In contrast, male infants showed vulnerability to a lower number of adversities at birth, but those who survived to adulthood were negatively affected by both early life individual and cumulative adversity. Our study shows profound immediate effects of insults  on female infant cohorts and suggests that affected female adults are more robust. In contrast, adult males who experienced harsh conditions early in life showed an increased mortality risk at older ages as expected from hypotheses within the life course perspective. Our analysis suggests sex-specific selection pressures on life histories and highlights the need for studies addressing the effects of early life adversity across multiple life stages.

The biology of aging in a social world: Insights from free-ranging rhesus macaques

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.

The provisioned primate: patterns of obesity across lemurs, monkeys, apes and humans

Non-human primates are potentially informative but underutilized species for investigating obesity. I examined patterns of obesity across the Primate order, calculating the ratio of body mass in captivity to that in the wild. This index, relative body mass, for n = 40 non-human primates (mean ± s.d.: females: 1.28 ± 0.30, range 0.67-1.78, males: 1.24 ± 0.28, range 0.70-1.97) overlapped with a reference value for humans (women: 1.52, men: 1.44). Among non-human primates, relative body mass was unrelated to dietary niche, and was marginally greater among female cohorts of terrestrial species. Males and females had similar relative body masses, but species with greater sexual size dimorphism (male/female mass) in wild populations had comparatively larger female body mass in captivity. Provisioned populations in wild and free-ranging settings had similar relative body mass to those in research facilities and zoos. Compared to the wild, captive diets are unlikely to be low in protein or fat, or high in carbohydrate, suggesting these macronutrients are not driving overeating in captive populations. Several primate species, including chimpanzees, a sister-species to humans, had relative body masses similar to humans. Humans are not unique in the propensity to overweight and obesity. This article is part of a discussion meeting issue 'Causes of obesity: theories, conjectures and evidence (Part II)'.

A single-cell multi-omic atlas spanning the adult rhesus macaque brain

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.

Hurricanes affect diversification among individual life courses of a primate population

Extreme climatic events may influence individual-level variability in phenotypes, survival and reproduction, and thereby drive the pace of evolution. Climate models predict increases in the frequency of intense hurricanes, but no study has measured their impact on individual life courses within animal populations. We used 45 years of demographic data of rhesus macaques to quantify the influence of major hurricanes on reproductive life courses using multiple metrics of dynamic heterogeneity accounting for life course variability and life-history trait variances. To reduce intraspecific competition, individuals may explore new reproductive stages during years of major hurricanes, resulting in higher temporal variation in reproductive trajectories. Alternatively, individuals may opt for a single optimal life-history strategy due to trade-offs between survival and reproduction. Our results show that heterogeneity in reproductive life courses increased by 4% during years of major hurricanes, despite a 2% reduction in the asymptotic growth rate due to an average decrease in mean fertility and survival by that is, shortened life courses and reduced reproductive output. In agreement with this, the population is expected to achieve stable population dynamics faster after being perturbed by a hurricane (ρ = 1.512 ; 95% CI: 1.488, 1.538), relative to ordinary yearsρ = 1.482 ; 1.475 , 1.490 . Our work suggests that natural disasters force individuals into new demographic roles to potentially reduce competition during unfavourable environments where mean reproduction and survival are compromised. Variance in lifetime reproductive success and longevity are differently affected by hurricanes, and such variability is mostly driven by survival.

The biology of aging in a social world:insights from free-ranging rhesus macaques

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 Island, which is home to the world's longest continuously monitored free-ranging population of rhesus macaques. 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.

Multiregion transcriptomic profiling of the primate brain reveals signatures of aging and the social environment

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.

Variation and heritability of retinal cone ratios in a free-ranging population of rhesus macaques

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.

Social connections predict brain structure in a multidimensional free-ranging primate society

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.

Rhesus macaques compensate for reproductive delay following ecological adversity early in life

Adversity early in life can shape the reproductive potential of individuals through negative effects on health and life span. However, long-lived populations with multiple reproductive events may present alternative life history strategies to optimize reproductive schedules and compensate for shorter life spans. Here, we quantify the effects of major hurricanes and density dependence as sources of early-life ecological adversity on Cayo Santiago rhesus macaque female reproduction and decompose their effects onto the mean age-specific fertility, reproductive pace, and lifetime reproductive success (LRS). Females experiencing major hurricanes exhibit a delayed reproductive debut but maintain the pace of reproduction past debut and show a higher mean fertility during prime reproductive ages, relative to unaffected females. Increasing density at birth is associated to a decrease in mean fertility and reproductive pace, but such association is absent at intermediate densities. When combined, our study reveals that hurricanes early in life predict a delay-overshoot pattern in mean age-specific fertility that supports the maintenance of LRS. In contrast to predictive adaptive response models of accelerated reproduction, this long-lived population presents a novel reproductive strategy where females who experience major natural disasters early in life ultimately overcome their initial reproductive penalty with no major negative fitness outcomes. Density presents a more complex relation with reproduction that suggests females experiencing a population regulated at intermediate densities early in life will escape density dependence and show optimized reproductive schedules. Our results support hypotheses about life history trade-offs in which adversity-affected females ensure their future reproductive potential by allocating more energy to growth or maintenance processes at younger adult ages.

Population Numbers and Reproductive Health

A recent study published in The Lancet predicts a remarkable drop in population numbers following a peak that will be reached by 2064. A unique feature of the upcoming population drop is that it will be almost exclusively caused by decreased reproduction, rather than factors that increase rates of mortality. The reasons for decreased reproduction are also unique, as, unlike previous centuries, limited reproduction today is hardly due to a shortage in resources. In other words, the predicted population drop is almost exclusively due to changes in reproductive behavior and reproductive physiology. Today, global changes in reproductive behavior are mostly explained by social sciences in a framework of demographic transition hypotheses, while changes in reproductive physiology are usually attributed to effects of endocrine-disrupting pollutants. This review outlines a complementary/alternative hypothesis, which connects reproductive trends with population densities. Numerous wildlife and experimental studies of a broad range of animal species have demonstrated that reproductive behavior and reproductive physiology are negatively controlled via endocrine and neural signaling in response to increasing population densities. The causal chain of this control system, although not fully understood, includes suppression of every level of hypothalamic-pituitary-gonadal cascade by hypothalamic-pituitary-adrenal axis, activated in response to increasing stress of social interactions. This paper discusses evidence in support of a hypothesis that current trends in reproductive physiology and behavior may be partly explained by increasing population densities. Better understanding of the causal chain involved in reproduction suppression by population density-related factors may help in developing interventions to treat infertility and other reproductive conditions.

Sociality predicts individual variation in the immunity of free-ranging rhesus macaques

Social integration and social status can substantially affect an individual’s health and survival. One route through which this occurs is by altering immune function, which can be highly sensitive to changes in the social environment. However, we currently have limited understanding of how sociality influences markers of immunity in naturalistic populations where social dynamics can be fully realized. To address this gap, we asked if social integration and social status in free-ranging rhesus macaques (Macaca mulatta) predict anatomical and physiological markers of immunity. We used data on agonistic interactions to determine social status, and social network analysis of grooming interactions to generate measures of individual variation in social integration. As measures of immunity, we included the size of two of the major organs involved in the immune response, the spleen and liver, and counts of three types of blood cells (red blood cells, platelets, and white blood cells). Controlling for body mass and age, we found that neither social status nor social integration predicted the size of anatomical markers of immunity. However, individuals that were more socially connected, i.e., with more grooming partners, had lower numbers of white blood cells than their socially isolated counterparts, indicating lower levels of inflammation with increasing levels of integration. These results build upon and extend our knowledge of the relationship between sociality and the immune system in humans and captive animals to free-ranging primates, demonstrating generalizability of the beneficial role of social integration on health.

Age and sex-associated variation in the multi-site microbiome of an entire social group of free-ranging rhesus macaques

BACKGROUND

An individual's microbiome changes over the course of its lifetime, especially during infancy, and again in old age. Confounding factors such as diet and healthcare make it difficult to disentangle the interactions between age, health, and microbial changes in humans. Animal models present an excellent opportunity to study age- and sex-linked variation in the microbiome, but captivity is known to influence animal microbial abundance and composition, while studies of free-ranging animals are typically limited to studies of the fecal microbiome using samples collected non-invasively. Here, we analyze a large dataset of oral, rectal, and genital swabs collected from 105 free-ranging rhesus macaques (Macaca mulatta, aged 1 month-26 years), comprising one entire social group, from the island of Cayo Santiago, Puerto Rico. We sequenced 16S V4 rRNA amplicons for all samples.

RESULTS

Infant gut microbial communities had significantly higher relative abundances of Bifidobacterium and Bacteroides and lower abundances of Ruminococcus, Fibrobacter, and Treponema compared to older age groups, consistent with a diet high in milk rather than solid foods. The genital microbiome varied widely between males and females in beta-diversity, taxonomic composition, and predicted functional profiles. Interestingly, only penile, but not vaginal, microbiomes exhibited distinct age-related changes in microbial beta-diversity, taxonomic composition, and predicted functions. Oral microbiome composition was associated with age, and was most distinctive between infants and other age classes.

CONCLUSIONS

Across all three body regions, with notable exceptions in the penile microbiome, while infants were distinctly different from other age groups, microbiomes of adults were relatively invariant, even in advanced age. While vaginal microbiomes were exceptionally stable, penile microbiomes were quite variable, especially at the onset of reproductive age. Relative invariance among adults, including elderly individuals, is contrary to findings in humans and mice. We discuss potential explanations for this observation, including that age-related microbiome variation seen in humans may be related to changes in diet and lifestyle. Video abstract.

Hurricane-induced demographic changes in a non-human primate population

Major disturbance events can have large impacts on the demography and dynamics of animal populations. Hurricanes are one example of an extreme climatic event, predicted to increase in frequency due to climate change, and thus expected to be a considerable threat to population viability. However, little is understood about the underlying demographic mechanisms shaping population response following these extreme disturbances. Here, we analyse 45 years of the most comprehensive free-ranging non-human primate demographic dataset to determine the effects of major hurricanes on the variability and maintenance of long-term population fitness. For this, we use individual-level data to build matrix population models and perform perturbation analyses. Despite reductions in population growth rate mediated through reduced fertility, our study reveals a demographic buffering during hurricane years. As long as survival does not decrease, our study shows that hurricanes do not result in detrimental effects at the population level, demonstrating the unbalanced contribution of survival and fertility to population fitness in long-lived animal populations.

Higher early life mortality with lower infant body mass in a free-ranging primate

Traits that reflect the amount of energy allocated to offspring by mothers, such as infant body mass, are predicted to have long-lasting effects on offspring fitness. In very long-lived species, such as anthropoid primates, where long-lasting and obligate parental care is required for successful recruitment of offspring, there are few studies on the fitness implications of low body mass among infants. Using body mass data collected from 253 free-ranging rhesus macaque Macaca mulatta infants on Cayo Santiago, Puerto Rico, we examined if lower infant body mass predicts lower chance of survival through to reproductive maturation (4th year of life). We also used data on inter-birth intervals and suckling behaviours to determine whether the duration of maternal care was adjusted to infant body mass. Rhesus macaque infants experienced on average 5% reduced hazard of death for an increase in body mass of 0.1 SD (~100 g) above the mean within their age-sex class. The positive association between body mass and early life survival was most pronounced in the 1st year of life. Infant body mass tended to be lower if mothers were young or old, but the link between infant body mass and early life survival remained after controlling for maternal age. This finding suggests that maternal effects on early life survival such as maternal age may act through their influence on infant body mass. Mothers of heavier infants were less likely to be delayed in subsequent reproduction, but the estimated association slightly overlapped with zero. The timing of the last week of suckling did not differ by infant body mass. Using infant body mass data that has been rarely available from free-ranging primates, our study provides comparative evidence to strengthen the existing body of literature on the fitness implications of variation in infant body mass.

Individuals in urban dwelling primate species face unequal benefits associated with living in an anthropogenic environment

In primates, living in an anthropogenic environment can significantly improve an individual's fitness, which is likely attributed to access to anthropogenic food resources. However, in non-professionally provisioned groups, few studies have examined whether individual attributes, such as dominance rank and sex, affect primates' ability to access anthropogenic food. Here, we investigated whether rank and sex explain individual differences in the proportion of anthropogenic food consumed by macaques. We observed 319 individuals living in nine urban groups across three macaque species. We used proportion of anthropogenic food in the diet as a proxy of access to those food resources. Males and high-ranking individuals in both sexes had significantly higher proportions of anthropogenic food in their diets than other individuals. We speculate that unequal access to anthropogenic food resources further increases within-group competition, and may limit fitness benefits in an anthropogenic environment to certain individuals.

Which male and female characteristics influence the probability of extragroup paternities in rhesus macaques, Macaca mulatta?

Extragroup paternity (EGP) is found across a wide range of species and may entail reproductive benefits, but may also entail costs to both sexes. While population and group parameters affecting the degree of EGPs are relatively well established, less is known about the individual characteristics that make males and females engage in alternative reproductive tactics such as EGP. Applying a combination of long-term demographic and genetic data from the rhesus macaque population of Cayo Santiago (Puerto Rico, U.S.A.), we investigate which male and female characteristics influence the probability of EGP to better understand the circumstances that shape the distribution and occurrence of EGP. Our results show that, against our expectations, higher-ranking females were more likely to produce EGP offspring than lower- ranking females. The probability of producing extragroup offspring was not significantly related to female or male age, male tenure or previous reproductive success. Furthermore, genetic relatedness between the parents did not affect the production of extragroup offspring, but extragroup offspring were more frequently produced early rather than late in a given mating season. Altogether, our analysis suggests that individual attributes and seasonal aspects create different opportunities and preferences for engaging in EGP as an alternative reproductive tactic. The observed patterns of EGP in rhesus macaques appear to be consistent with female mate choice for genetic benefits, which needs to be confirmed in future studies.

Drivers of Diversification in Individual Life Courses

Heterogeneity in life courses among individuals of a population influences the speed of adaptive evolutionary processes, but it is less clear how biotic and abiotic environmental fluctuations influence such heterogeneity. We investigate principal drivers of variability in sequence of stages during an individual's life in a stage-structured population. We quantify heterogeneity by measuring population entropy of a Markov chain, which computes the rate of diversification of individual life courses. Using individual data of a primate population, we show that density regulates the stage composition of the population but that its entropy and the generating moments of heterogeneity are independent of density. This lack of influence of density on heterogeneity is due to neither low year-to-year variation in entropy nor differences in survival among stages but is rather due to differences in stage transitions. Our analysis thus shows that well-known classical ecological selective forces, such as density regulation, are not linked to potential selective forces governing heterogeneity through underlying stage dynamics. Despite evolution acting heavily on individual variability in fitness components, our understanding is poor whether observed heterogeneity is adaptive and how it evolves and is maintained. Our analysis illustrates how entropy represents a more integrated measure of diversity compared to the population structural composition, giving us new insights about the underlying drivers of individual heterogeneity within populations and potential evolutionary mechanisms.

Drivers of Diversification in Individual Life Courses

Heterogeneity in life courses among individuals of a population influences the speed of adaptive evolutionary processes, but it is less clear how biotic and abiotic environmental fluctuations influence such heterogeneity. We investigate principal drivers of variability in sequence of stages during an individual's life in a stage-structured population. We quantify heterogeneity by measuring population entropy of a Markov chain, which computes the rate of diversification of individual life courses. Using individual data of a primate population, we show that density regulates the stage composition of the population but that its entropy and the generating moments of heterogeneity are independent of density. This lack of influence of density on heterogeneity is due to neither low year-to-year variation in entropy nor differences in survival among stages but is rather due to differences in stage transitions. Our analysis thus shows that well-known classical ecological selective forces, such as density regulation, are not linked to potential selective forces governing heterogeneity through underlying stage dynamics. Despite evolution acting heavily on individual variability in fitness components, our understanding is poor whether observed heterogeneity is adaptive and how it evolves and is maintained. Our analysis illustrates how entropy represents a more integrated measure of diversity compared to the population structural composition, giving us new insights about the underlying drivers of individual heterogeneity within populations and potential evolutionary mechanisms.

Population Genetic Structure of the Cayo Santiago Colony of Rhesus Macaques (Macaca mulatta)

The rhesus macaque population at Cayo Santiago increases annually and is in urgent need of control. In-depth assessments of the colony's population genetic and pedigree structures provide a starting point for improving the colony's long-term management program. We evaluated the degree of genetic variation and coefficients of inbreeding and kinship of the Cayo Santiago colony by using pedigree and short tandem repeat (STR) data from 4738 rhesus macaques, which represent 7 extant social groups and a group of migrant males. Information on each animal's parentage, sex, birth date, and date of death or removal from the island were used to generate estimates of mean kinship, kinship value, gene value, genome uniqueness (GU), founder equivalents (fe), and founder genome equivalents (fg). Pedigree and STR analyses revealed that the social groups have not differentiated genetically from each other due to male-mediated gene flow (that is, FST estimates were in the negative range) and exhibit sufficient genetic variation, with mean estimates of allele numbers and observed and expected heterozygosity of 6.57, 0.72, and 0.70, respectively. Estimates of GU, fe, and fg show that a high effective number of founders has affected the colony's current genetic structure in a positive manner. As demographic changes occur, genetic and pedigree matrices need to be monitored consistently to ensure the health and wellbeing of the Cayo Santiago colony.

The Population Genetic Composition of Conventional and SPF Colonies of Rhesus Macaques (Macaca mulatta) at the Caribbean Primate Research Center

The SPF breeding program at the Caribbean Primate Research Center supplies Indian-origin rhesus macaques of known genetic and virologic background for biomedical research. In this study, population genetic analyses using 14 short tandem-repeat sequences showed that the SPF colony has remained genetically homogenous over time, with sufficient amounts of heterozygosity and minimal stratification from its founders. Intergenerational studies indicated that an average of 7 alleles have been retained, inbreeding levels have remained low, and the degree of Indian ancestry is one of the highest among several national primate research centers. The relative low genetic diversity in the free-ranging population as well as in the captive SPF and conventional colonies when compared with that of other primate centers indicates that the free-ranging population, from which the captive-colony animals were derived, has experienced significant founder effects and genetic drift during the years after its establishment. This study supports the historical origin of the free-ranging population and confirms the high value of this resource for biomedical research. Current genetic diversity levels within the SPF colony can be ensured with the practice of colony management approaches such as equalizing male:female ratios in each SPF breeding group and increasing breeding group sizes. Introducing new Indian-origin macaques from other captive colonies might help to maximize the genetic diversity of the breeding stock. Furthermore, genetic estimates must be used to rank breeders according to their genetic value or their genome uniqueness to increase founder-genome representation and curb future genetic bottlenecks and allele loss.