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.

Mediterranean diet protects against a neuroinflammatory cortical transcriptome: Associations with brain volumetrics, peripheral inflammation, social isolation, and anxiety in nonhuman primates (Macaca fascicularis).

Mediterranean diets may be neuroprotective and prevent cognitive decline relative to Western diets; however, the underlying biology is poorly understood. We assessed the effects of Western versus Mediterranean-like diets on RNAseq-generated transcriptional profiles in lateral temporal cortex and their relationships with longitudinal changes in neuroanatomy, circulating monocyte gene expression, and observations of social isolation and anxiety in 38 socially-housed, middle-aged female cynomolgus macaques (Macaca fascicularis). Diet resulted in differential expression of seven transcripts (FDR < 0.05). Cyclin dependent kinase 14 (CDK14), a proinflammatory regulator, was lower in the Mediterranean group. The remaining six transcripts [i.e., "lunatic fringe" (LFNG), mannose receptor C type 2 (MRC2), solute carrier family 3 member 2 (SLCA32), butyrophilin subfamily 2 member A1 (BTN2A1), katanin regulatory subunit B1 (KATNB1), and transmembrane protein 268 (TMEM268)] were higher in cortex of the Mediterranean group and generally associated with anti-inflammatory/neuroprotective pathways. KATNB1 encodes a subcomponent of katanin, important in maintaining microtubule homeostasis. BTN2A1 is involved in immunomodulation of γδ T-cells which have anti-neuroinflammatory and neuroprotective effects. CDK14, LFNG, MRC2, and SLCA32 are associated with inflammatory pathways. The latter four differentially expressed cortex transcripts were associated with peripheral monocyte transcript levels, neuroanatomical changes determined by MRI, and with social isolation and anxiety. These results provide important insights into the potential mechanistic processes linking diet, peripheral and central inflammation, and behavior. Collectively, our results provide evidence that, relative to Western diets, Mediterranean diets confer protection against peripheral and central inflammation which is reflected in preserved brain structure and socioemotional behavior. Ultimately, such protective effects may confer resilience to the development of neuropathology and associated disease.

Vascularization underlies differences in sexually selected skin coloration in a wild primate.

Male reproductive competition can select for condition-dependent, conspicuous traits that signal some aspect of fighting ability and facilitate assessment of potential rivals. However, the underlying mechanisms that link the signal to a male's current condition are difficult to investigate in wild populations, often requiring invasive experimental manipulation. Here, we use digital photographs and chest skin samples to investigate the mechanisms of a visual signal used in male competition in a wild primate, the red chest patch in geladas (Theropithecus gelada). We analysed photographs collected during natural (n = 144) and anaesthetized conditions (n = 38) to understand variability in male and female chest redness, and we used chest skin biopsies (n = 38) to explore sex differences in gene expression. Male and female geladas showed similar average redness, but males exhibited a wider within-individual range in redness under natural conditions. These sex differences were also reflected at the molecular level, with 10.5% of genes exhibiting significant sex differences in expression. Subadult males exhibited intermediate gene expression patterns between adult males and females, pointing to mechanisms underlying the development of the red chest patch. We found that genes more highly expressed in males were associated with blood vessel development and maintenance but not with androgen or oestrogen activity. Together, our results suggest male gelada redness variability is driven by increased blood vessel branching in the chest skin, providing a potential link between male chest redness and current condition as increased blood circulation to exposed skin could lead to heat loss in the cold, high-altitude environment of geladas.

A novel statovirus identified in fecal samples from wild geladas in the Ethiopian highlands.

We present a novel statovirus in geladas (Theropithecus gelada), graminivorous primates endemic to the Ethiopian highlands. Using a high-throughput sequencing approach, we identified contiguous sequences in feces from two adult female geladas in the Simien Mountains National Park, Ethiopia, that share similarities to statoviruses. Our phylogenetic analysis of the whole genome, as well as the RNA-dependent RNA polymerase (RdRp) and capsid protein (CP) amino acid sequences, revealed that the gelada statoviruses cluster with those from other primates (laboratory populations of Macaca nemestrina and Macaca mulatta). As the first report of statovirus in wild primates, this finding contributes to our understanding of the phylogenetic and geographic distribution of statoviruses and their hosts.

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

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

Trichromacy increases fruit intake rates of wild capuchins (Cebus capucinus imitator).

Intraspecific color vision variation is prevalent among nearly all diurnal monkeys in the neotropics and is seemingly a textbook case of balancing selection acting to maintain genetic polymorphism. Clear foraging advantages to monkeys with trichromatic vision over those with dichromatic "red-green colorblind" vision have been observed in captive studies; however, evidence of trichromatic advantage during close-range foraging has been surprisingly scarce in field studies, perhaps as a result of small sample sizes and strong impacts of environmental or individual variation on foraging performance. To robustly test the effects of color vision type on foraging efficiency in the wild, we conducted an extensive study of dichromatic and trichromatic white-faced capuchin monkeys (Cebus capucinus imitator), controlling for plant-level and monkey-level variables that may affect fruit intake rates. Over the course of 14 months, we collected behavioral data from 72 monkeys in Sector Santa Rosa, Costa Rica. We analyzed 19,043 fruit feeding events within 1,602 foraging bouts across 27 plant species. We find that plant species, color conspicuity category, and monkey age class significantly impact intake rates, while sex does not. When plant species and age are controlled for, we observe that trichromats have higher intake rates than dichromats for plant species with conspicuously colored fruits. This study provides clear evidence of trichromatic advantage in close-range fruit feeding in wild monkeys. Taken together with previous reports of dichromatic advantage for finding cryptic foods, our results illuminate an important aspect of balancing selection maintaining primate opsin polymorphism.

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.

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.

Mediterranean Diet Protects Against a Neuroinflammatory Cortical Transcriptome: Associations with Brain Volumetrics, Peripheral Inflammation, Social Isolation and Anxiety.

INTRODUCTION: Mediterranean diets may be neuroprotective and prevent cognitive decline relative to Western diets, however the underlying biology is poorly understood. METHODS: We assessed the effects of Western vs. Mediterranean-like diets on RNAseq generated transcriptional profiles in temporal cortex and their relationships with changes in MRI neuroimaging phenotypes, circulating monocyte gene expression, and observations of social isolation and anxiety in 38 socially-housed, middle-aged female cynomolgus macaques. RESULTS: Diet resulted in differential expression of seven transcripts (FDR<0.05). Cyclin dependent kinase 14 ( CDK14 ), a proinflammatory regulator, was lower in the Mediterranean group. The remaining six transcripts [i.e., "lunatic fringe" ( LFNG ), mannose receptor C type 2 ( MRC2 ), solute carrier family 3 member 2 ( SLCA32 ), butyrophilin subfamily 2 member A1 ( BTN2A1 ), katanin regulatory subunit B1 ( KATNB1 ), and transmembrane protein 268 ( TMEM268 )] were higher in cortex of the Mediterranean group and generally associated with anti-inflammatory/neuroprotective pathways. KATNB1 encodes a subcomponent of katanin, important in maintaining microtubule homeostasis. BTN2A1 is involved in immunomodulation of γδ T-cells which have anti-neuroinflammatory and neuroprotective effects. CDK14 , LFNG , MRC2, and SLCA32 are associated with inflammatory pathways. The latter four differentially expressed cortex transcripts were associated with monocyte transcript levels, changes in AD-relevant brain volumes determined by MRI over the course of the study, and social isolation and anxiety. CDK14 was positively correlated with monocyte inflammatory transcripts, changes in total brain, gray matter, cortical gray matter volumes, and time alone and anxious behavior, and negatively correlated with changes in total white matter and cerebrospinal fluid (CSF) volumes. In contrast, LFNG , MRC2 , and SLCA32 were negatively correlated with monocyte inflammatory transcripts and changes in total gray matter volume, and positively correlated with CSF volume changes, and SLCA32 was negatively correlated with time alone. DISCUSSION: Collectively, our results suggest that relative to Western diets, Mediterranean diets confer protection against peripheral and central inflammation which is reflected in preserved brain structure and behavior.

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.

Human microglia show unique transcriptional changes in Alzheimer's disease.

Microglia, the innate immune cells of the brain, influence Alzheimer's disease (AD) progression and are potential therapeutic targets. However, microglia exhibit diverse functions, the regulation of which is not fully understood, complicating therapeutics development. To better define the transcriptomic phenotypes and gene regulatory networks associated with AD, we enriched for microglia nuclei from 12 AD and 10 control human dorsolateral prefrontal cortices (7 males and 15 females, all aged >60 years) before single-nucleus RNA sequencing. Here we describe both established and previously unrecognized microglial molecular phenotypes, the inferred gene networks driving observed transcriptomic change, and apply trajectory analysis to reveal the putative relationships between microglial phenotypes. We identify microglial phenotypes more prevalent in AD cases compared with controls. Further, we describe the heterogeneity in microglia subclusters expressing homeostatic markers. Our study demonstrates that deep profiling of microglia in human AD brain can provide insight into microglial transcriptional changes associated with AD.

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

Genome-wide coancestry reveals details of ancient and recent male-driven reticulation in baboons.

Baboons (genus Papio ) are a morphologically and behaviorally diverse clade of catarrhine monkeys that have experienced hybridization between phenotypically and genetically distinct phylogenetic species. We used high coverage whole genome sequences from 225 wild baboons representing 19 geographic localities to investigate population genomics and inter-species gene flow. Our analyses provide an expanded picture of evolutionary reticulation among species and reveal novel patterns of population structure within and among species, including differential admixture among conspecific populations. We describe the first example of a baboon population with a genetic composition that is derived from three distinct lineages. The results reveal processes, both ancient and recent, that produced the observed mismatch between phylogenetic relationships based on matrilineal, patrilineal, and biparental inheritance. We also identified several candidate genes that may contribute to species-specific phenotypes. One-Sentence Summary: Genomic data for 225 baboons reveal novel sites of inter-species gene flow and local effects due to differences in admixture.

Genome-wide coancestry reveals details of ancient and recent male-driven reticulation in baboons.

Baboons (genus Papio) are a morphologically and behaviorally diverse clade of catarrhine monkeys that have experienced hybridization between phenotypically and genetically distinct phylogenetic species. We used high-coverage whole-genome sequences from 225 wild baboons representing 19 geographic localities to investigate population genomics and interspecies gene flow. Our analyses provide an expanded picture of evolutionary reticulation among species and reveal patterns of population structure within and among species, including differential admixture among conspecific populations. We describe the first example of a baboon population with a genetic composition that is derived from three distinct lineages. The results reveal processes, both ancient and recent, that produced the observed mismatch between phylogenetic relationships based on matrilineal, patrilineal, and biparental inheritance. We also identified several candidate genes that may contribute to species-specific phenotypes.

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.

Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas.

Primates have adapted to numerous environments and lifestyles but very few species are native to high elevations. Here we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of haematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared with baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult haemoglobin but found that gelada haemoglobin does not exhibit markedly altered oxygenation properties compared with lowland primates. We also found that geladas at high altitude do not exhibit elevated blood haemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythaemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia research.

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.

Pleistocene diversification of living squirrel monkeys (Saimiri spp.) inferred from complete mitochondrial genome sequences.

In order to enhance our understanding of the evolutionary history of squirrel monkeys (Saimiri spp.), we newly sequenced and analyzed data from seven complete mitochondrial genomes representing six squirrel monkey taxa. While previous studies have lent insights into the taxonomy and phylogeny of the genus, phylogenetic relationships and divergence date estimates among major squirrel monkey clades remain unclear. Using maximum likelihood and Bayesian procedures, we inferred a highly resolved phylogenetic tree with strong support for a sister relationship between Saimiri boliviensis and all other Saimiri, for monophyly of Saimiri oerstedii and Saimiri sciureus sciureus, and for Saimiri sciureus macrodon as the sister lineage to the S. oerstedii/S. s. sciureus clade. We inferred that crown lineages for extant squirrel monkeys diverged around 1.5 million years ago (MYA) in the Pleistocene Epoch, with other major clades diverging between 0.9 and 1.1 MYA. Our results suggest a relatively recent timeline of squirrel monkey evolution and challenge previous conceptions about the diversification of the genus and its expansion into Central America.