Natural disaster and immunological aging in a nonhuman primate.

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.

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.

Determination of major histocompatibility class I and class II genetic composition of the Caribbean Primate Center specific pathogen-free rhesus macaque (Macaca mulatta) colony based on massively parallel sequencing.

BACKGROUND: Knowledge of major histocompatibility complex (MHC) composition and distribution in rhesus macaque colonies is critical for management strategies that maximize the utility of this model for biomedical research. METHODS: Variation within the Mamu-A and Mamu-B (class I) and DRB, DQA/B, and DPA/B (class II) regions of 379 animals from the Caribbean Primate Research Center's (CPRC) specific pathogen free (SPF) colony was examined using massively parallel sequencing. RESULTS: Analyses of the 7 MHC loci revealed a background of Indian origin with high levels of variation despite past genetic bottlenecks. All loci exhibited mutual linkage disequilibria while conforming to Hardy-Weinberg expectations suggesting the achievement of mutation-selection balance. CONCLUSION: The CPRC's SPF colony is a significant resource for research on AIDS and other infectious agents. Characterizing colony-wide MHC variability facilitates the breeding and selection of animals bearing desired haplotypes and increases the investigator's ability to understand the immune responses mounted by these animals.

An alphavirus vector-based tetravalent dengue vaccine induces a rapid and protective immune response in macaques that differs qualitatively from immunity induced by live virus infection.

Despite many years of research, a dengue vaccine is not available, and the more advanced live attenuated vaccine candidate in clinical trials requires multiple immunizations with long interdose periods and provides low protective efficacy. Here, we report important contributions to the development of a second-generation dengue vaccine. First, we demonstrate that a nonpropagating vaccine vector based on Venezuelan equine encephalitis virus replicon particles (VRP) expressing two configurations of dengue virus E antigen (subviral particles [prME] and soluble E dimers [E85]) successfully immunized and protected macaques against dengue virus, while antivector antibodies did not interfere with a booster immunization. Second, compared to prME-VRP, E85-VRP induced neutralizing antibodies faster, to higher titers, and with improved protective efficacy. Third, this study is the first to map antigenic domains and specificities targeted by vaccination versus natural infection, revealing that, unlike prME-VRP and live virus, E85-VRP induced only serotype-specific antibodies, which predominantly targeted EDIII, suggesting a protective mechanism different from that induced by live virus and possibly live attenuated vaccines. Fourth, a tetravalent E85-VRP dengue vaccine induced a simultaneous and protective response to all 4 serotypes after 2 doses given 6 weeks apart. Balanced responses and protection in macaques provided further support for exploring the immunogenicity and safety of this vaccine candidate in humans.

Lamina Cribrosa Microstructure in Nonhuman Primates With Naturally Occurring Peripapillary Retinal Nerve Fiber Layer Thinning.

Purpose: The lamina cribrosa (LC) is hypothesized to be the site of initial axonal damage in glaucoma with the circumpapillary retinal nerve fiber layer thickness (RNFL-T) widely used as a standard metric for quantifying the glaucomatous damage. The purpose of this study was to determine in vivo, 3-dimensional (3D) differences in the microstructure of the LC in eyes of nonhuman primates (NHPs) with naturally occurring glaucoma. Methods: Spectral-domain optical coherence tomography (OCT) scans (Leica, Chicago, IL, USA) of the optic nerve head were acquired from a colony of 50 adult rhesus monkeys suspected of having high prevalence of glaucoma. The RNFL-T was analyzed globally and in quadrants using a semi-automated segmentation software. From a set of 100 eyes, 18 eyes with the thinnest global RNFL-T were selected as the study group and 18 eyes with RNFL-T values around the 50th percentile were used as controls. A previously described automated segmentation algorithm was used for LC microstructure analysis. Parameters included beam thickness, pore diameter and their ratio (beam-to-pore ratio [BPR]), pore area and shape parameters, beam and pore volume, and connective tissue volume fraction (CTVF; beam volume/total volume). The LC microstructure was analyzed globally and in the following volumetric sectors: quadrants, central and peripheral lamina, and three depth slabs (anterior, middle, and posterior). Results: Although no significant difference was detected between groups for age, weight, or disc size, the study group had significantly thinner RNFL than the control group (P < 0.01). The study group had significantly smaller global and sectoral pore diameter and larger BPR compared with the control group. Across eyes, the global RNFL-T was associated positively with pore diameter globally. BPR and CTVF were significantly and negatively associated with the corresponding RNFL-T in the superior quadrant. Conclusions: Global and sectoral microstructural differences were detected when comparing thin and normal RNFL-T eyes. Whether these LC differences are the cause of RNFL damage or the result of remodeling of the LC requires further investigation. Translational Relevance: Our findings indicate structural alterations in the LC of NHP exhibiting natural thinning of the RNFL, a common characteristic of glaucomatous damage.

Early life adversity has sex-dependent effects on survival across the lifespan in rhesus macaques.

Exposure to early life adversity is linked to detrimental fitness outcomes across taxa. Owing to the challenges of collecting longitudinal data, direct evidence for long-term fitness effects of early life adversity from long-lived species remains relatively scarce. Here, we test the effects of early life adversity on male and female longevity in a free-ranging population of rhesus macaques (Macaca mulatta) on Cayo Santiago, Puerto Rico. We leveraged six decades of data to quantify the relative importance of 10 forms of early life adversity for 6599 macaques. Individuals that experienced more early life adversity died earlier than those that experienced less adversity. Mortality risk was highest during early life, defined as birth to 4 years old, but heightened mortality risk was also present in macaques that survived to adulthood. Females and males were affected differently by some forms of adversity, and these differences might be driven by varying energetic demands and dispersal patterns. Our results show that the fitness consequences of early life adversity are not uniform across individuals but vary as a function of the type of adversity, timing and social context, and thus contribute to our limited but growing understanding of the evolution of early life sensitivities.This article is part of the discussion meeting issue 'Understanding age and society using natural populations'.

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.

Immune cell composition varies by age, sex and exposure to social adversity in free-ranging Rhesus Macaques.

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.

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.

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.

Age-Related Differences in Ocular Features of a Naturalistic Free-Ranging Population of Rhesus Macaques.

Purpose: Rhesus macaques (Macaca mulatta) are the premier nonhuman primate model for studying human health and disease. We investigated if age was associated with clinically relevant ocular features in a large cohort of free-ranging rhesus macaques from Cayo Santiago, Puerto Rico. Methods: We evaluated 120 rhesus macaques (73 males, 47 females) from 0 to 29 years old (mean ± SD: 12.6 ± 6.4) from September to December 2021. The ophthalmic evaluation included intraocular pressure (IOP) assessment, corneal pachymetry, biomicroscopy, A-scan biometry, automated refraction, and fundus photography after pupil dilation. The associations of age with the outcomes were investigated through multilevel mixed-effects models adjusted for sex and weight. Results: On average, IOP, pachymetry, axial length, and automated refraction spherical equivalent were 18.37 ± 4.68 mmHg, 474.43 ± 32.21 µm, 19.49 ± 1.24 mm, and 0.30 ± 1.70 diopters (D), respectively. Age was significantly associated with pachymetry (ß coefficient = -1.20; 95% confidence interval [CI], -2.27 to -0.14; P = 0.026), axial length (ß coefficient = 0.03; 95% CI, 0.01 to 0.05; P = 0.002), and spherical equivalent (ß coefficient = -0.12; 95% CI, -0.22 to -0.02; P = 0.015). No association was detected between age and IOP. The prevalence of cataracts in either eye was 10.83% (95% CI, 6.34-17.89) and was significantly associated with age (odds ratio [OR] = 1.20; 95% CI, 1.06-1.36; P = 0.004). Retinal drusen in either eye was observed in 15.00% (95% CI, 9.60-22.68) of animals, which was also significantly associated with age (OR = 1.14; 95% CI, 1.02-1.27; P = 0.020). Conclusions: Rhesus macaques exhibit age-related ocular associations similar to those observed in human aging, including decreased corneal thickness, increased axial length, myopic shift, and higher prevalence of cataract and retinal drusen.

Early life adversity has sex-dependent effects on survival across the lifespan in rhesus macaques.

Exposure to adversity during early life is linked to lasting detrimental effects on evolutionary fitness across many taxa. However, due to the challenges of collecting longitudinal data, especially in species where one sex disperses, direct evidence from long-lived species remains relatively scarce. Here we test the effects of early life adversity on male and female longevity in a free-ranging population of rhesus macaques (Macaca mulatta) at Cayo Santiago, Puerto Rico. We leveraged six decades of data to quantify the relative importance of ten forms of early life adversity for 6,599 macaques (3,230 male, 3,369 female), with a smaller sample size (N=299) for one form of adversity (maternal social isolation) which required high-resolution behavioral data. We found that individuals who experienced more early life adversity died earlier than those who experienced less adversity. Mortality risk was highest during early life, defined as birth to four years old, suggesting acute survival effects of adversity, but heightened mortality risk was also present in macaques who survived to adulthood. Females and males were affected differently by some forms of adversity, and these differences might be driven by varying energetic demands, female philopatry, and male dispersal. By leveraging data on thousands of macaques collected over decades, our results show that the fitness consequences of early life adversity are not uniform across individuals but vary as a function of the type of adversity, timing, and social context, and thus contribute to our limited but growing understanding of the evolution of early life sensitivities in long-lived species.

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.

A live dengue virus vaccine carrying a chimeric envelope glycoprotein elicits dual DENV2-DENV4 serotype-specific immunity.

The four dengue virus serotypes co-circulate globally and cause significant human disease. Dengue vaccine development is challenging because some virus-specific antibodies are protective, while others are implicated in enhanced viral replication and more severe disease. Current dengue tetravalent vaccines contain four live attenuated serotypes formulated to theoretically induce balanced protective immunity. Among the number of vaccine candidates in clinical trials, only Dengvaxia is licensed for use in DENV seropositive individuals. To simplify live-virus vaccine design, we identify co-evolutionary constraints inherent in flavivirus virion assembly and design chimeric viruses to replace domain II (EDII) of the DENV2 envelope (E) glycoprotein with EDII from DENV4. The chimeric DENV2/4EDII virus replicates efficiently in vitro and in vivo. In male macaques, a single inoculation of DENV2/4EDII induces type-specific neutralizing antibodies to both DENV2 and DENV4, thereby providing a strategy to simplify DENV vaccine design by utilizing a single bivalent E glycoprotein immunogen for two DENV serotypes.

Best practices for genotype imputation from low-coverage sequencing data in natural populations.

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 r2 ≥ 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.

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.

Infection order outweighs the role of CD4+ T cells in tertiary flavivirus exposure.

The link between CD4+ T and B cells during immune responses to DENV and ZIKV and their roles in cross-protection during heterologous infection is an active area of research. Here we used CD4+ lymphocyte depletions to dissect the impact of cellular immunity on humoral responses during a tertiary flavivirus infection in macaques. We show that CD4+ depletion in DENV/ZIKV-primed animals followed by DENV resulted in dysregulated adaptive immune responses. We show a delay in DENV-specific IgM/IgG antibody titers and binding and neutralization in the DENV/ZIKV-primed CD4-depleted animals but not in ZIKV/DENV-primed CD4-depleted animals. This study confirms the critical role of CD4+ cells in priming an early effective humoral response during sequential flavivirus infections. Our work here suggests that the order of flavivirus exposure affects the outcome of a tertiary infection. Our findings have implications for understanding the complex flavivirus immune responses and for the development of effective flavivirus vaccines.

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.

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.

Vitamin D Content in Milk of the Rhesus Monkey.

OBJECTIVE: Milk is the first and continued source of ingested Vitamin D. Extensive studies have been carried out in humans measuring Vitamin D in lactating mothers but to date few values have been obtained for milk of non-human primates and none for rhesus monkeys. Consequently. we have determined Vitamin D and antirachitic activity (ARA) in milk samples obtained from 21 rhesus monkeys. METHODS: Lactating dams were sampled by hand-stripping. 25(OH)D2, Vitamin D2, 25(OH)D3, Vitamin D3 and ARA were assessed in foremilk using LC-MS/MS techniques. RESULTS: 25(OH)D2 and Vitamin D2 were below detectible limits (<0.5 ng/g), 25(OH) D3 =4.2 ± 1.8 ng/ml, Vitamin D3 = 6.1 ± 3.1 ng/ml and ARA = 1080 ± 480 IU/L. CONCLUSION: This is the first report of content of Vitamin D and ARA activity in foremilk of the rhesus monkey and can serve as a reference for future studies.