Social status is one of the strongest predictors of human disease risk and mortality, and it also influences Darwinian fitness in social mammals more generally. To understand the biological basis of these effects, we combined genomics with a social status manipulation in female rhesus macaques to investigate how status alters immune function. We demonstrate causal but largely plastic social status effects on immune cell proportions, cell type–specific gene expression levels, and the gene expression response to immune challenge. Further, we identify specific transcription factor signaling pathways that explain these differences, including low-status–associated polarization of the Toll-like receptor 4 signaling pathway toward a proinflammatory response. Our findings provide insight into the direct biological effects of social inequality on immune function, thus improving our understanding of social gradients in health.

Variation in the social environment is a fundamental component of many vertebrate societies. In humans and other primates, adverse social environments often translate into lasting physiological costs. The biological mechanisms associated with these effects are therefore of great interest, both for understanding the evolutionary impacts of social behavior and in the context of human health. However, large gaps remain in our understanding of the mechanisms that mediate these effects at the molecular level. Here we addressed these questions by leveraging the power of an experimental system that consisted of 10 social groups of female macaques, in which each individual's social status (i.e., dominance rank) could be experimentally controlled. Using this paradigm, we show that dominance rank results in a widespread, yet plastic, imprint on gene regulation, such that peripheral blood mononuclear cell gene expression data alone predict social status with 80% accuracy. We investigated the mechanistic basis of these effects using cell type-specific gene expression profiling and glucocorticoid resistance assays, which together contributed to rank effects on gene expression levels for 694 (70%) of the 987 rankrelated genes. We also explored the possible contribution of DNA methylation levels to these effects, and identified global associations between dominance rank and methylation profiles that suggest epigenetic flexibility in response to status-related behavioral cues. Together, these results illuminate the importance of the molecular response to social conditions, particularly in the immune system, and demonstrate a key role for gene regulation in linking the social environment to individual physiology.inflammation | social gradient | differential gene expression | sociogenomics

BackgroundThe rhesus macaque (Macaca mulatta) is a key species for advancing biomedical research. Like all draft mammalian genomes, the draft rhesus assembly (rheMac2) has gaps, sequencing errors and misassemblies that have prevented automated annotation pipelines from functioning correctly. Another rhesus macaque assembly, CR_1.0, is also available but is substantially more fragmented than rheMac2 with smaller contigs and scaffolds. Annotations for these two assemblies are limited in completeness and accuracy. High quality assembly and annotation files are required for a wide range of studies including expression, genetic and evolutionary analyses.ResultsWe report a new de novo assembly of the rhesus macaque genome (MacaM) that incorporates both the original Sanger sequences used to assemble rheMac2 and new Illumina sequences from the same animal. MacaM has a weighted average (N50) contig size of 64 kilobases, more than twice the size of the rheMac2 assembly and almost five times the size of the CR_1.0 assembly. The MacaM chromosome assembly incorporates information from previously unutilized mapping data and preliminary annotation of scaffolds. Independent assessment of the assemblies using Ion Torrent read alignments indicates that MacaM is more complete and accurate than rheMac2 and CR_1.0. We assembled messenger RNA sequences from several rhesus tissues into transcripts which allowed us to identify a total of 11,712 complete proteins representing 9,524 distinct genes. Using a combination of our assembled rhesus macaque transcripts and human transcripts, we annotated 18,757 transcripts and 16,050 genes with complete coding sequences in the MacaM assembly. Further, we demonstrate that the new annotations provide greatly improved accuracy as compared to the current annotations of rheMac2. Finally, we show that the MacaM genome provides an accurate resource for alignment of reads produced by RNA sequence expression studies.ConclusionsThe MacaM assembly and annotation files provide a substantially more complete and accurate representation of the rhesus macaque genome than rheMac2 or CR_1.0 and will serve as an important resource for investigators conducting next-generation sequencing studies with nonhuman primates.ReviewersThis article was reviewed by Dr. Lutz Walter, Dr. Soojin Yi and Dr. Kateryna Makova.

The contribution of pre-mRNA processing mechanisms to the regulation of immune responses remains poorly studied despite emerging examples of their role as regulators of immune defenses. We sought to investigate the role of mRNA processing in the cellular responses of human macrophages to live bacterial infections. Here, we used mRNA sequencing to quantify gene expression and isoform abundances in primary macrophages from 60 individuals, before and after infection with Listeria monocytogenes and Salmonella typhimurium. In response to both bacteria we identified thousands of genes that significantly change isoform usage in response to infection, characterized by an overall increase in isoform diversity after infection. In response to both bacteria, we found global shifts towards (i) the inclusion of cassette exons and (ii) shorter 3’ UTRs, with near-universal shifts towards usage of more upstream polyadenylation sites. Using complementary data collected in non-human primates, we show that these features are evolutionarily conserved among primates. Following infection, we identify candidate RNA processing factors whose expression is associated with individual-specific variation in isoform abundance. Finally, by profiling microRNA levels, we show that 3’ UTRs with reduced abundance after infection are significantly enriched for target sites for particular miRNAs. These results suggest that the pervasive usage of shorter 3’ UTRs is a mechanism for particular genes to evade repression by immune-activated miRNAs. Collectively, our results suggest that dynamic changes in RNA processing may play key roles in the regulation of innate immune responses.