Within-individual changes reveal increasing social selectivity with age in rhesus macaques.

Accumulating evidence in humans and other mammals suggests older individuals tend to have smaller social networks. Uncovering the cause of these declines can inform how changes in social relationships with age affect health and fitness in later life. While age-based declines in social networks have been thought to be detrimental, physical and physiological limitations associated with age may lead older individuals to adjust their social behavior and be more selective in partner choice. Greater selectivity with age has been shown in humans, but the extent to which this phenomenon occurs across the animal kingdom remains an open question. Using longitudinal data from a population of rhesus macaques on Cayo Santiago, we provide compelling evidence in a nonhuman animal for within-individual increases in social selectivity with age. Our analyses revealed that adult female macaques actively reduced the size of their networks as they aged and focused on partners previously linked to fitness benefits, including kin and partners to whom they were strongly and consistently connected earlier in life. Females spent similar amounts of time socializing as they aged, suggesting that network shrinkage does not result from lack of motivation or ability to engage, nor was this narrowing driven by the deaths of social partners. Furthermore, females remained attractive companions and were not isolated by withdrawal of social partners. Taken together, our results provide rare empirical evidence for social selectivity in nonhumans, suggesting that patterns of increasing selectivity with age may be deeply rooted in primate evolution.

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.

Social history and exposure to pathogen signals modulate social status effects on gene regulation in rhesus macaques.

Social experience is an important predictor of disease susceptibility and survival in humans and other social mammals. Chronic social stress is thought to generate a proinflammatory state characterized by elevated antibacterial defenses and reduced investment in antiviral defense. Here we manipulated long-term social status in female rhesus macaques to show that social subordination alters the gene expression response to ex vivo bacterial and viral challenge. As predicted by current models, bacterial lipopolysaccharide polarizes the immune response such that low status corresponds to higher expression of genes in NF-κB-dependent proinflammatory pathways and lower expression of genes involved in the antiviral response and type I IFN signaling. Counter to predictions, however, low status drives more exaggerated expression of both NF-κB- and IFN-associated genes after cells are exposed to the viral mimic Gardiquimod. Status-driven gene expression patterns are linked not only to social status at the time of sampling, but also to social history (i.e., past social status), especially in unstimulated cells. However, for a subset of genes, we observed interaction effects in which females who fell in rank were more strongly affected by current social status than those who climbed the social hierarchy. Taken together, our results indicate that the effects of social status on immune cell gene expression depend on pathogen exposure, pathogen type, and social history-in support of social experience-mediated biological embedding in adulthood, even in the conventionally memory-less innate immune system.

Social ageing: exploring the drivers of late-life changes in social behaviour in mammals.

Social interactions help group-living organisms cope with socio-environmental challenges and are central to survival and reproductive success. Recent research has shown that social behaviour and relationships can change across the lifespan, a phenomenon referred to as 'social ageing'. Given the importance of social integration for health and well-being, age-dependent changes in social behaviour can modulate how fitness changes with age and may be an important source of unexplained variation in individual patterns of senescence. However, integrating social behaviour into ageing research requires a deeper understanding of the causes and consequences of age-based changes in social behaviour. Here, we provide an overview of the drivers of late-life changes in sociality. We suggest that explanations for social ageing can be categorized into three groups: changes in sociality that (a) occur as a result of senescence; (b) result from adaptations to ameliorate the negative effects of senescence; and/or (c) result from positive effects of age and demographic changes. Quantifying the relative contribution of these processes to late-life changes in sociality will allow us to move towards a more holistic understanding of how and why these patterns emerge and will provide important insights into the potential for social ageing to delay or accelerate other patterns of senescence.

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.

Social status alters chromatin accessibility and the gene regulatory response to glucocorticoid stimulation in rhesus macaques.

Low social status is an important predictor of disease susceptibility and mortality risk in humans and other social mammals. These effects are thought to stem in part from dysregulation of the glucocorticoid (GC)-mediated stress response. However, the molecular mechanisms that connect low social status and GC dysregulation to downstream health outcomes remain elusive. Here, we used an in vitro GC challenge to investigate the consequences of experimentally manipulated social status (i.e., dominance rank) for immune cell gene regulation in female rhesus macaques, using paired control and GC-treated peripheral blood mononuclear cell samples. We show that social status not only influences immune cell gene expression but also chromatin accessibility at hundreds of regions in the genome. Social status effects on gene expression were less pronounced following GC treatment than under control conditions. In contrast, social status effects on chromatin accessibility were stable across conditions, resulting in an attenuated relationship between social status, chromatin accessibility, and gene expression after GC exposure. Regions that were more accessible in high-status animals and regions that become more accessible following GC treatment were enriched for a highly concordant set of transcription factor binding motifs, including motifs for the GC receptor cofactor AP-1. Together, our findings support the hypothesis that social status alters the dynamics of GC-mediated gene regulation and identify chromatin accessibility as a mechanism involved in social stress-driven GC resistance. More broadly, they emphasize the context-dependent nature of social status effects on gene regulation and implicate epigenetic remodeling of chromatin accessibility as a contributing factor.

The Goldilocks effect: female geladas in mid-sized groups have higher fitness.

The cost-benefit ratio of group living is thought to vary with group size: individuals in 'optimally sized' groups should have higher fitness than individuals in groups that are either too large or too small. However, the relationship between group size and individual fitness has been difficult to establish for long-lived species where the number of groups studied is typically quite low. Here, we present evidence for optimal group size that maximizes female fitness in a population of geladas (Theropithecus gelada). Drawing on 14 years of demographic data, we found that females in small groups experienced the highest death rates, while females in mid-sized groups exhibited the highest reproductive performance. This group size effect on female reproductive performance was largely explained by variation in infant mortality (and, in particular, by infanticide from immigrant males) but not by variation in reproductive rates. Taken together, females in mid-sized groups are projected to attain optimal fitness due to conspecific infanticide and, potentially, predation. Our findings provide insight into how and why group size shapes fitness in long-lived species.

Estimating the heritability of cognitive traits across dog breeds reveals highly heritable inhibitory control and communication factors.

Trait heritability is necessary for evolution by both natural and artificial selection, yet we know little about the heritability of cognitive traits. Domestic dogs are a valuable study system for questions regarding the evolution of phenotypic diversity due to their extraordinary intraspecific variation. While previous studies have investigated morphological and behavioral variation across dog breeds, few studies have systematically assessed breed differences in cognition. We integrated data from Dognition.com-a citizen science project on dog cognition-with breed-averaged genetic data from published sources to estimate the among-breed heritability of cognitive traits using mixed models. The resulting dataset included 11 cognitive measures for 1508 adult dogs across 36 breeds. A factor analysis yielded four factors interpreted as reflecting inhibitory control, communication, memory, and physical reasoning. Narrow-sense among-breed heritability estimates-reflecting the proportion of cognitive variance attributable to additive genetic variation-revealed that scores on the inhibitory control and communication factors were highly heritable (inhibitory control: h2 = 0.70; communication: h2 = 0.39), while memory and physical reasoning were less heritable (memory: h2 = 0.17; physical reasoning: h2 = 0.21). Although the heritability of inhibitory control is partially explained by body weight, controlling for breed-average weight still yields a high heritability estimate (h2 = 0.50), while other factors are minimally affected. Our results indicate that cognitive phenotypes in dogs covary with breed relatedness and suggest that cognitive traits have strong potential to undergo selection. The highest heritabilities were observed for inhibitory control and communication, both of which are hypothesized to have been altered by domestication.

Age influences domestic dog cognitive performance independent of average breed lifespan.

Across mammals, increased body size is positively associated with lifespan. However, within species, this relationship is inverted. This is well illustrated in dogs (Canis familiaris), where larger dogs exhibit accelerated life trajectories: growing faster and dying younger than smaller dogs. Similarly, some age-associated traits (e.g., growth rate and physiological pace of aging) exhibit accelerated trajectories in larger breeds. Yet, it is unknown whether cognitive performance also demonstrates an accelerated life course trajectory in larger dogs. Here, we measured cognitive development and aging in a cross-sectional study of over 4000 dogs from 66 breeds using nine memory and decision-making tasks performed by citizen scientists as part of the Dognition project. Specifically, we tested whether cognitive traits follow a compressed (accelerated) trajectory in larger dogs, or the same trajectory for all breeds, which would result in limited cognitive decline in larger breeds. We found that all breeds, regardless of size or lifespan, tended to follow the same quadratic trajectory of cognitive aging-with a period of cognitive development in early life and decline in later life. Taken together, our results suggest that cognitive performance follows similar age-related trajectories across dog breeds, despite remarkable variation in developmental rates and lifespan.

Highly heritable and functionally relevant breed differences in dog behaviour.

Variation across dog breeds presents a unique opportunity to investigate the evolution and biological basis of complex behavioural traits. We integrated behavioural data from more than 14 000 dogs from 101 breeds with breed-averaged genotypic data (n = 5697 dogs) from over 100 000 loci in the dog genome. We found high levels of among-breed heritability for 14 behavioural traits (the proportion of trait variance attributable to genetic similarity among breeds). We next identified 131 single nucleotide polymorphisms associated with breed differences in behaviour, which were found in genes that are highly expressed in the brain and enriched for neurobiological functions and developmental processes, suggesting that they may be functionally associated with behavioural differences. Our results shed light on the heritability and genetic architecture of complex behavioural traits and identify dogs as a powerful model in which to address these questions.

Deconstructing sociality: the types of social connections that predict longevity in a group-living primate.

Many species use social interactions to cope with challenges in their environment and a growing number of studies show that individuals which are well-connected to their group have higher fitness than socially isolated individuals. However, there are many ways to be 'well-connected' and it is unclear which aspects of sociality drive fitness benefits. Being well-connected can be conceptualized in four main ways: individuals can be socially integrated by engaging in a high rate of social behaviour or having many partners; they can have strong and stable connections to favoured partners; they can indirectly connect to the broader group structure; or directly engage in a high rate of beneficial behaviours, such as grooming. In this study, we use survival models and long-term data in adult female rhesus macaques (Macaca mulatta) to compare the fitness outcomes of multiple measures of social connectedness. Females that maintained strong connections to favoured partners had the highest relative survival probability, as did females well-integrated owing to forming many weak connections. We found no survival benefits to being structurally well-connected or engaging in high rates of grooming. Being well-connected to favoured partners could provide fitness benefits by, for example, increasing the efficacy of coordinated or mutualistic behaviours.

Social affiliation predicts mitochondrial DNA copy number in female rhesus macaques.

In many social mammals, social adversity predicts compromised health and reduced fitness. These effects are thought to be driven in part by chronic social stress, but their molecular underpinnings are not well understood. Recent work suggests that chronic stress can affect mitochondrial copy number, heteroplasmy rates and function. Here, we tested the first two possibilities for the first time in non-human primates. We manipulated dominance rank in captive female rhesus macaques ( n = 45), where low rank induces chronic social stress, and measured mitochondrial DNA (mtDNA) copy number and heteroplasmy in five peripheral blood mononuclear cell types from each study subject. We found no effect of dominance rank on either mtDNA copy number or heteroplasmy rates. However, grooming rate, a measure of affiliative social behaviour predicted by high social status, was positively associated with mtDNA copy number in B cells, cytotoxic T cells and monocytes. Our results suggest that social interactions can influence mtDNA regulation in immune cells. Further, they indicate the importance of considering both affiliative and competitive interactions in investigating this relationship.

Developmental responses to early-life adversity: Evolutionary and mechanistic perspectives.

Adverse ecological and social conditions during early life are known to influence development, with rippling effects that may explain variation in adult health and fitness. The adaptive function of such developmental plasticity, however, remains relatively untested in long-lived animals, resulting in much debate over which evolutionary models are most applicable. Furthermore, despite the promise of clinical interventions that might alleviate the health consequences of early-life adversity, research on the proximate mechanisms governing phenotypic responses to adversity have been largely limited to studies on glucocorticoids. Here, we synthesize the current state of research on developmental plasticity, discussing both ultimate and proximate mechanisms. First, we evaluate the utility of adaptive models proposed to explain developmental responses to early-life adversity, particularly for long-lived mammals such as humans. In doing so, we highlight how parent-offspring conflict complicates our understanding of whether mothers or offspring benefit from these responses. Second, we discuss the role of glucocorticoids and a second physiological system-the gut microbiome-that has emerged as an additional, clinically relevant mechanism by which early-life adversity can influence development. Finally, we suggest ways in which nonhuman primates can serve as models to study the effects of early-life adversity, both from evolutionary and clinical perspectives.

Composition and stability of the vervet monkey milk microbiome.

The human milk microbiome is vertically transmitted to offspring during the postnatal period and has emerged as a critical driver of infant immune and metabolic development. Despite this importance in humans, the milk microbiome of nonhuman primates remains largely unexplored. This dearth of comparative work precludes our ability to understand how species-specific differences in the milk microbiome may differentially drive maternal effects and limits how translational models can be used to understand the role of vertically transmitted milk microbes in human development. Here, we present the first culture-independent data on the milk microbiome of a nonhuman primate. We collected milk and matched fecal microbiome samples at early and late lactation from a cohort of captive lactating vervet monkeys (N = 15). We found that, similar to humans, the vervet monkey milk microbiome comprises a shared community of taxa that are universally present across individuals. However, unlike in humans, this shared community is dominated by the genera Lactobacillus, Bacteroides, and Prevotella. We also found that, in contrast to previous culture-dependent studies in humans, the vervet milk microbiome exhibits greater alpha-diversity than the gut microbiome across lactation. Finally, we did not find support for the translocation of microbes from the gut to the mammary gland within females (i.e., "entero-mammary pathway"). Taken together, our results show that the vervet monkey milk microbiome is taxonomically diverse, distinct from the gut microbiome, and largely stable. These findings demonstrate that the milk microbiome is a unique substrate that may selectively favor the establishment and persistence of particular microbes across lactation and highlights the need for future experimental studies on the origin of microbes in milk.

Quantifying uncertainty due to fission-fusion dynamics as a component of social complexity.

Groups of animals (including humans) may show flexible grouping patterns, in which temporary aggregations or subgroups come together and split, changing composition over short temporal scales, (i.e. fission and fusion). A high degree of fission-fusion dynamics may constrain the regulation of social relationships, introducing uncertainty in interactions between group members. Here we use Shannon's entropy to quantify the predictability of subgroup composition for three species known to differ in the way their subgroups come together and split over time: spider monkeys (Ateles geoffroyi), chimpanzees (Pan troglodytes) and geladas (Theropithecus gelada). We formulate a random expectation of entropy that considers subgroup size variation and sample size, against which the observed entropy in subgroup composition can be compared. Using the theory of set partitioning, we also develop a method to estimate the number of subgroups that the group is likely to be divided into, based on the composition and size of single focal subgroups. Our results indicate that Shannon's entropy and the estimated number of subgroups present at a given time provide quantitative metrics of uncertainty in the social environment (within which social relationships must be regulated) for groups with different degrees of fission-fusion dynamics. These metrics also represent an indirect quantification of the cognitive challenges posed by socially dynamic environments. Overall, our novel methodological approach provides new insight for understanding the evolution of social complexity and the mechanisms to cope with the uncertainty that results from fission-fusion dynamics.

Weak effects of common genetic variation in oxytocin and vasopressin receptor genes on rhesus macaque social behavior.

The neuropeptides oxytocin (OT) and arginine vasopressin (AVP) influence pair bonding, attachment, and sociality, as well as anxiety and stress responses in humans and other mammals. The effects of these peptides are mediated by genetic variability in their associated receptors, OXTR and the AVPR gene family. However, the role of these genes in regulating social behaviors in non-human primates is not well understood. To address this question, we examined whether genetic variation in the OT receptor gene OXTR and the AVP receptor genes AVPR1A and AVPR1B influence naturally-occurring social behavior in free-ranging rhesus macaques-gregarious primates that share many features of their biology and social behavior with humans. We assessed rates of social behavior across 3,250 hr of observational behavioral data from 201 free-ranging rhesus macaques on Cayo Santiago island in Puerto Rico, and used genetic sequence data to identify 25 OXTR, AVPR1A, and AVPR1B single-nucleotide variants (SNVs) in the population. We used an animal model to estimate the effects of 12 SNVs (n = 3 OXTR; n = 5 AVPR1A; n = 4 AVPR1B) on rates of grooming, approaches, passive contact, contact aggression, and non-contact aggression, given and received. Though we found evidence for modest heritability of these behaviors, estimates of effect sizes of the selected SNVs were close to zero, indicating that common OXTR and AVPR variation contributed little to social behavior in these animals. Our results are consistent with recent findings in human genetics that the effects of individual common genetic variants on complex phenotypes are generally small.

Pervasive Effects of Aging on Gene Expression in Wild Wolves.

Gene expression levels change as an individual ages and responds to environmental conditions. With the exception of humans, such patterns have principally been studied under controlled conditions, overlooking the array of developmental and environmental influences that organisms encounter under conditions in which natural selection operates. We used high-throughput RNA sequencing (RNA-Seq) of whole blood to assess the relative impacts of social status, age, disease, and sex on gene expression levels in a natural population of gray wolves (Canis lupus). Our findings suggest that age is broadly associated with gene expression levels, whereas other examined factors have minimal effects on gene expression patterns. Further, our results reveal evolutionarily conserved signatures of senescence, such as immunosenescence and metabolic aging, between wolves and humans despite major differences in life history and environment. The effects of aging on gene expression levels in wolves exhibit conservation with humans, but the more rapid expression differences observed in aging wolves is evolutionarily appropriate given the species' high level of extrinsic mortality due to intraspecific aggression. Some expression changes that occur with age can facilitate physical age-related changes that may enhance fitness in older wolves. However, the expression of these ancestral patterns of aging in descendant modern dogs living in highly modified domestic environments may be maladaptive and cause disease. This work provides evolutionary insight into aging patterns observed in domestic dogs and demonstrates the applicability of studying natural populations to investigate the mechanisms of aging.