Identification of constrained sequence elements across 239 primate genomes.

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

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.

Variation in predicted COVID-19 risk among lemurs and lorises.

The novel coronavirus SARS-CoV-2, which in humans leads to the disease COVID-19, has caused global disruption and more than 2 million fatalities since it first emerged in late 2019. As we write, infection rates are at their highest point globally and are rising extremely rapidly in some areas due to more infectious variants. The primary target of SARS-CoV-2 is the cellular receptor angiotensin-converting enzyme-2 (ACE2). Recent sequence analyses of the ACE2 gene predict that many nonhuman primates are also likely to be highly susceptible to infection. However, the anticipated risk is not equal across the Order. Furthermore, some taxonomic groups show high ACE2 amino acid conservation, while others exhibit high variability at this locus. As an example of the latter, analyses of strepsirrhine primate ACE2 sequences to date indicate large variation among lemurs and lorises compared to other primate clades despite low sampling effort. Here, we report ACE2 gene and protein sequences for 71 individual strepsirrhines, spanning 51 species and 19 genera. Our study reinforces previous results while finding additional variability in other strepsirrhine species, and suggests several clades of lemurs have high potential susceptibility to SARS-CoV-2 infection. Troublingly, some species, including the rare and endangered aye-aye (Daubentonia madagascariensis), as well as those in the genera Avahi and Propithecus, may be at high risk. Given that lemurs are endemic to Madagascar and among the primates at highest risk of extinction globally, further understanding of the potential threat of COVID-19 to their health should be a conservation priority. All feasible actions should be taken to limit their exposure to SARS-CoV-2.

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.

Correction to 'Diversity and evolution of the primate skin microbiome'.

[This corrects the article DOI: 10.1098/rspb.2015.2586.].

Diversity and evolution of the primate skin microbiome.

Skin microbes play a role in human body odour, health and disease. Compared with gut microbes, we know little about the changes in the composition of skin microbes in response to evolutionary changes in hosts, or more recent behavioural and cultural changes in humans. No studies have used sequence-based approaches to consider the skin microbe communities of gorillas and chimpanzees, for example. Comparison of the microbial associates of non-human primates with those of humans offers unique insights into both the ancient and modern features of our skin-associated microbes. Here we describe the microbes found on the skin of humans, chimpanzees, gorillas, rhesus macaques and baboons. We focus on the bacterial and archaeal residents in the axilla using high-throughput sequencing of the 16S rRNA gene. We find that human skin microbial communities are unique relative to those of other primates, in terms of both their diversity and their composition. These differences appear to reflect both ancient shifts during millions of years of primate evolution and more recent changes due to modern hygiene.

Analysis of synaptic gene expression in the neocortex of primates reveals evolutionary changes in glutamatergic neurotransmission.

Increased relative brain size characterizes the evolution of primates, suggesting that enhanced cognition plays an important part in the behavioral adaptations of this mammalian order. In addition to changes in brain anatomy, cognition can also be regulated by molecular changes that alter synaptic function, but little is known about modifications of synapses in primate brain evolution. The aim of the current study was to investigate the expression patterns and evolution of 20 synaptic genes from the prefrontal cortex of 12 primate species. The genes investigated included glutamate receptors, scaffolding proteins, synaptic vesicle components, as well as factors involved in synaptic vesicle release and structural components of the nervous system. Our analyses revealed that there have been significant changes during primate brain evolution in the components of the glutamatergic signaling pathway in terms of gene expression, protein expression, and promoter sequence changes. These results could entail functional modifications in the regulation of specific genes related to processes underlying learning and memory.

Comparative analysis of the primate X-inactivation center region and reconstruction of the ancestral primate XIST locus.

Here we provide a detailed comparative analysis across the candidate X-Inactivation Center (XIC) region and the XIST locus in the genomes of six primates and three mammalian outgroup species. Since lemurs and other strepsirrhine primates represent the sister lineage to all other primates, this analysis focuses on lemurs to reconstruct the ancestral primate sequences and to gain insight into the evolution of this region and the genes within it. This comparative evolutionary genomics approach reveals significant expansion in genomic size across the XIC region in higher primates, with minimal size alterations across the XIST locus itself. Reconstructed primate ancestral XIC sequences show that the most dramatic changes during the past 80 million years occurred between the ancestral primate and the lineage leading to Old World monkeys. In contrast, the XIST locus compared between human and the primate ancestor does not indicate any dramatic changes to exons or XIST-specific repeats; rather, evolution of this locus reflects small incremental changes in overall sequence identity and short repeat insertions. While this comparative analysis reinforces that the region around XIST has been subject to significant genomic change, even among primates, our data suggest that evolution of the XIST sequences themselves represents only small lineage-specific changes across the past 80 million years.

Genetic comparisons yield insight into the evolution of enamel thickness during human evolution.

Enamel thickness varies substantially among extant hominoids and is a key trait with significance for interpreting dietary adaptation, life history trajectory, and phylogenetic relationships. There is a strong link in humans between enamel formation and mutations in the exons of the four genes that code for the enamel matrix proteins and the associated protease. The evolution of thick enamel in humans may have included changes in the regulation of these genes during tooth development. The cis-regulatory region in the 5' flank (upstream non-coding region) of MMP20, which codes for enamelysin, the predominant protease active during enamel secretion, has previously been shown to be under strong positive selection in the lineages leading to both humans and chimpanzees. Here we examine evidence for positive selection in the 5' flank and 3' flank of AMELX, AMBN, ENAM, and MMP20. We contrast the human sequence changes with other hominoids (chimpanzees, gorillas, orangutans, gibbons) and rhesus macaques (outgroup), a sample comprising a range of enamel thickness. We find no evidence for positive selection in the protein-coding regions of any of these genes. In contrast, we find strong evidence for positive selection in the 5' flank region of MMP20 and ENAM along the lineage leading to humans, and in both the 5' flank and 3' flank regions of MMP20 along the lineage leading to chimpanzees. We also identify putative transcription factor binding sites overlapping some of the species-specific nucleotide sites and we refine which sections of the up- and downstream putative regulatory regions are most likely to harbor important changes. These non-coding changes and their potential for differential regulation by transcription factors known to regulate tooth development may offer insight into the mechanisms that allow for rapid evolutionary changes in enamel thickness across closely-related species, and contribute to our understanding of the enamel phenotype in hominoids.

A molecular phylogeny of living primates.

Comparative genomic analyses of primates offer considerable potential to define and understand the processes that mold, shape, and transform the human genome. However, primate taxonomy is both complex and controversial, with marginal unifying consensus of the evolutionary hierarchy of extant primate species. Here we provide new genomic sequence (~8 Mb) from 186 primates representing 61 (~90%) of the described genera, and we include outgroup species from Dermoptera, Scandentia, and Lagomorpha. The resultant phylogeny is exceptionally robust and illuminates events in primate evolution from ancient to recent, clarifying numerous taxonomic controversies and providing new data on human evolution. Ongoing speciation, reticulate evolution, ancient relic lineages, unequal rates of evolution, and disparate distributions of insertions/deletions among the reconstructed primate lineages are uncovered. Our resolution of the primate phylogeny provides an essential evolutionary framework with far-reaching applications including: human selection and adaptation, global emergence of zoonotic diseases, mammalian comparative genomics, primate taxonomy, and conservation of endangered species.

Taking identity-by-descent analysis into the wild: Estimating realized relatedness in free-ranging macaques.

Biological relatedness is a key consideration in studies of behavior, population structure, and trait evolution. Except for parent-offspring dyads, pedigrees capture relatedness imperfectly. The number and length of DNA segments that are identical-by-descent (IBD) yield the most precise estimates of relatedness. Here, we leverage novel methods for estimating locus-specific IBD from low coverage whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4-6× coverage data from a rhesus macaque (Macaca mulatta) population with available long-term pedigree data, we show that we can call the number and length of IBD segments across the genome with high accuracy even at 0.5× coverage. The resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. They identify cryptic genetic relatives that are not represented in the pedigree and reveal elevated recombination rates in females relative to males, which allows us to discriminate maternal and paternal kin using genotype data alone. Our findings represent a breakthrough in the ability to understand the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Yearly variation coupled with social interactions shape the skin microbiome in free-ranging rhesus macaques.

While skin microbes are known to mediate human health and disease, there has been minimal research on the interactions between skin microbiota, social behavior, and year-to-year effects in non-human primates-important animal models for translational biomedical research. To examine these relationships, we analyzed skin microbes from 78 rhesus macaques living on Cayo Santiago Island, Puerto Rico. We considered age, sex, and social group membership, and characterized social behavior by assessing dominance rank and patterns of grooming as compared to nonsocial behaviors. To measure the effects of a shifting environment, we sampled skin microbiota (based on sequence analysis of the 16S rRNA V4 region) and assessed weather across sampling periods between 2013 and 2015. We hypothesized that, first, monkeys with similar social behavior and/or in the same social group would possess similar skin microbial composition due, in part, to physical contact, and, second, microbial diversity would differ across sampling periods. We found significant phylum-level differences between social groups in the core microbiome as well as an association between total grooming rates and alpha diversity in the complete microbiome, but no association between microbial diversity and measures of rank or other nonsocial behaviors. We also identified alpha and beta diversity differences in microbiota and differential taxa abundance across two sampling periods. Our findings indicate that social dynamics interact with yearly environmental changes to shape the skin microbiota in rhesus macaques, with potential implications for understanding the factors affecting the microbiome in humans, which share many biological and social characteristics with these animals. IMPORTANCE Primate studies are valuable for translational and evolutionary insights into the human microbiome. The majority of primate microbiome studies focus on the gut, so less is known about the factors impacting the microbes on skin and how their links affect health and behavior. Here, we probe the impact of social interactions and the yearly environmental changes on food-provisioned, free-ranging monkeys living on a small island. We expected animals that lived together and groomed each other would have more similar microbes on their skin, but surprisingly found that the external environment was a stronger influence on skin microbiome composition. These findings have implications for our understanding of the human skin microbiome, including potential manipulations to improve health and treat disease.

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.

A global catalog of whole-genome diversity from 233 primate species.

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

The landscape of tolerated genetic variation in humans and primates.

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole-genome sequencing data for 809 individuals from 233 primate species and identified 4.3 million common protein-altering variants with orthologs in humans. We show that these variants can be inferred to have nondeleterious effects in humans based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

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.

Variation in predicted COVID-19 risk among lemurs and lorises.

The novel coronavirus SARS-CoV-2, which in humans leads to the disease COVID-19, has caused global disruption and more than 1.5 million fatalities since it first emerged in late 2019. As we write, infection rates are currently at their highest point globally and are rising extremely rapidly in some areas due to more infectious variants. The primary viral target is the cellular receptor angiotensin-converting enzyme-2 (ACE2). Recent sequence analyses of the ACE2 gene predicts that many nonhuman primates are also likely to be highly susceptible to infection. However, the anticipated risk is not equal across the Order. Furthermore, some taxonomic groups show high ACE2 amino acid conservation, while others exhibit high variability at this locus. As an example of the latter, analyses of strepsirrhine primate ACE2 sequences to date indicate large variation among lemurs and lorises compared to other primate clades despite low sampling effort. Here, we report ACE2 gene and protein sequences for 71 individual strepsirrhines, spanning 51 species and 19 genera. Our study reinforces previous results and finds additional variability in other strepsirrhine species, and suggests several clades of lemurs have high potential susceptibility to SARS-CoV-2 infection. Troublingly, some species, including the rare and Endangered aye-aye (Daubentonia madagascariensis), as well as those in the genera Avahi and Propithecus, may be at high risk. Given that lemurs are endemic to Madagascar and among the primates at highest risk of extinction globally, further understanding of the potential threat of COVID-19 to their health should be a conservation priority. All feasible actions should be taken to limit their exposure to SARS-CoV-2.

Primate comparative genomics: lemur biology and evolution.

Comparative genome sequencing projects are providing insight into aspects of genome biology that raise new questions and challenge existing paradigms. Placement in the phylogenetic tree can often be a major determinant of which organism to choose for study. Lemurs hold a key position at the base of the primate evolutionary tree and will be highly informative for the genomics community by offering comparisons of primate-specific characteristics and processes. Combining research in chromosome evolution, genome evolution and behavior with lemur comparative genomic sequencing will offer insights into many levels of primate evolution. We discuss the current state of lemur cytogenetic and phylogenetic analyses, and suggest how focusing more genomic efforts on lemurs will be beneficial to understanding human and primate evolution, as well as disease, and will contribute to conservation efforts.

The structure and evolution of centromeric transition regions within the human genome.

An understanding of how centromeric transition regions are organized is a critical aspect of chromosome structure and function; however, the sequence context of these regions has been difficult to resolve on the basis of the draft genome sequence. We present a detailed analysis of the structure and assembly of all human pericentromeric regions (5 megabases). Most chromosome arms (35 out of 43) show a gradient of dwindling transcriptional diversity accompanied by an increasing number of interchromosomal duplications in proximity to the centromere. At least 30% of the centromeric transition region structure originates from euchromatic gene-containing segments of DNA that were duplicatively transposed towards pericentromeric regions at a rate of six-seven events per million years during primate evolution. This process has led to the formation of a minimum of 28 new transcripts by exon exaptation and exon shuffling, many of which are primarily expressed in the testis. The distribution of these duplicated segments is nonrandom among pericentromeric regions, suggesting that some regions have served as preferential acceptors of euchromatic DNA.