Genetic diversity of 1,845 rhesus macaques improves genetic variation interpretation and identifies disease models.

Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes.

Mauritian-origin cynomolgus macaques (MCMs) serve as a powerful nonhuman primate model in biomedical research due to their unique genetic homogeneity, which simplifies experimental designs. Despite their extensive use, a comprehensive understanding of crucial immune-regulating gene families, particularly killer Ig-like receptors (KIR) and NK group 2 (NKG2), has been hindered by the lack of detailed genomic reference assemblies. In this study, we employ advanced long-read sequencing techniques to completely assemble eight KIR and seven NKG2 genomic haplotypes, providing an extensive insight into the structural and allelic diversity of these immunoregulatory gene clusters. Leveraging these genomic resources, we prototype a strategy for genotyping KIR and NKG2 using short-read, whole-exome capture data, illustrating the potential for cost-effective multilocus genotyping at colony scale. These results mark a significant enhancement for biomedical research in MCMs and underscore the feasibility of broad-scale genetic investigations.

The CARD8 inflammasome dictates HIV/SIV pathogenesis and disease progression.

While CD4+ T cell depletion is key to disease progression in people living with HIV and SIV-infected macaques, the mechanisms underlying this depletion remain incompletely understood, with most cell death involving uninfected cells. In contrast, SIV infection of "natural" hosts such as sooty mangabeys does not cause CD4+ depletion and AIDS despite high-level viremia. Here, we report that the CARD8 inflammasome is activated immediately after HIV entry by the viral protease encapsulated in incoming virions. Sensing of HIV protease activity by CARD8 leads to rapid pyroptosis of quiescent cells without productive infection, while T cell activation abolishes CARD8 function and increases permissiveness to infection. In humanized mice reconstituted with CARD8-deficient cells, CD4+ depletion is delayed despite high viremia. Finally, we discovered loss-of-function mutations in CARD8 from "natural hosts," which may explain the peculiarly non-pathogenic nature of these infections. Our study suggests that CARD8 drives CD4+ T cell depletion during pathogenic HIV/SIV infections.

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.

Whole Genome Analysis of SNV and Indel Polymorphism in Common Marmosets (Callithrix jacchus).

The common marmoset (Callithrix jacchus) is one of the most widely used nonhuman primate models of human disease. Owing to limitations in sequencing technology, early genome assemblies of this species using short-read sequencing suffered from gaps. In addition, the genetic diversity of the species has not yet been adequately explored. Using long-read genome sequencing and expert annotation, we generated a high-quality genome resource creating a 2.898 Gb marmoset genome in which most of the euchromatin portion is assembled contiguously (contig N50 = 25.23 Mbp, scaffold N50 = 98.2 Mbp). We then performed whole genome sequencing on 84 marmosets sampling the genetic diversity from several marmoset research centers. We identified a total of 19.1 million single nucleotide variants (SNVs), of which 11.9 million can be reliably mapped to orthologous locations in the human genome. We also observed 2.8 million small insertion/deletion variants. This dataset includes an average of 5.4 million SNVs per marmoset individual and a total of 74,088 missense variants in protein-coding genes. Of the 4956 variants orthologous to human ClinVar SNVs (present in the same annotated gene and with the same functional consequence in marmoset and human), 27 have a clinical significance of pathogenic and/or likely pathogenic. This important marmoset genomic resource will help guide genetic analyses of natural variation, the discovery of spontaneous functional variation relevant to human disease models, and the development of genetically engineered marmoset disease models.

Single-Cell RNA-Seq Identifies Pathways and Genes Contributing to the Hyperandrogenemia Associated with Polycystic Ovary Syndrome.

Polycystic ovary syndrome (PCOS) is a common endocrine disorder characterized by hyperandrogenemia of ovarian thecal cell origin, resulting in anovulation/oligo-ovulation and infertility. Our previous studies established that ovarian theca cells isolated and propagated from ovaries of normal ovulatory women and women with PCOS have distinctive molecular and cellular signatures that underlie the increased androgen biosynthesis in PCOS. To evaluate differences between gene expression in single-cells from passaged cultures of theca cells from ovaries of normal ovulatory women and women with PCOS, we performed single-cell RNA sequencing (scRNA-seq). Results from these studies revealed differentially expressed pathways and genes involved in the acquisition of cholesterol, the precursor of steroid hormones, and steroidogenesis. Bulk RNA-seq and microarray studies confirmed the theca cell differential gene expression profiles. The expression profiles appear to be directed largely by increased levels or activity of the transcription factors SREBF1, which regulates genes involved in cholesterol acquisition (LDLR, LIPA, NPC1, CYP11A1, FDX1, and FDXR), and GATA6, which regulates expression of genes encoding steroidogenic enzymes (CYP17A1) in concert with other differentially expressed transcription factors (SP1, NR5A2). This study provides insights into the molecular mechanisms underlying the hyperandrogenemia associated with PCOS and highlights potential targets for molecular diagnosis and therapeutic intervention.

Exome Sequencing Implicates DGKZ , ESRRA , and GXYLT1 for Modulating Granuloma Formation in Crohn Disease.

Non-caseating granulomas may indicate a more aggressive phenotype of Crohn disease (CD). Genetic associations of granulomatous CD (GCD) may help elucidate disease pathogenesis. Whole-exome sequencing was performed on peripheral blood-derived DNA from 17 pediatric patients with GCD and 19 with non-GCD (NGCD), and from an independent validation cohort of 44 GCD and 19 NGCD cases. PLINK (a tool set for whole-genome association and population-based linkage analyses) analysis was used to identify single nucleotide polymorphisms (SNPs) differentiating between groups, and subgroup allele frequencies were also compared to a public genomic database (gnomAD). The Combined Annotation Dependent Depletion scoring tool was used to predict deleteriousness of SNPs. Human leukocyte antigen (HLA) haplotype findings were compared to a control group (n = 8496). PLINK-based analysis between GCD and NGCD groups did not find consistently significant hits. gnomAD control comparisons, however, showed consistent subgroup associations with DGKZ , ESRRA , and GXYLT1 , genes that have been implicated in mammalian granulomatous inflammation. Our findings may guide future research and precision medicine.

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.

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.

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.

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 human. We show that these variants can be inferred to have non-deleterious effects in human 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. One Sentence Summary: Deep learning classifier trained on 4.3 million common primate missense variants predicts variant pathogenicity in humans.

SVhound: detection of regions that harbor yet undetected structural variation.

BACKGROUND: Recent population studies are ever growing in number of samples to investigate the diversity of a population or species. These studies reveal new polymorphism that lead to important insights into the mechanisms of evolution, but are also important for the interpretation of these variations. Nevertheless, while the full catalog of variations across entire species remains unknown, we can predict which regions harbor additional not yet detected variations and investigate their properties, thereby enhancing the analysis for potentially missed variants. RESULTS: To achieve this we developed SVhound ( https://github.com/lfpaulin/SVhound ), which based on a population level SVs dataset can predict regions that harbor unseen SV alleles. We tested SVhound using subsets of the 1000 genomes project data and showed that its correlation (average correlation of 2800 tests r = 0.7136) is high to the full data set. Next, we utilized SVhound to investigate potentially missed or understudied regions across 1KGP and CCDG. Lastly we also apply SVhound on a small and novel SV call set for rhesus macaque (Macaca mulatta) and discuss the impact and choice of parameters for SVhound. CONCLUSIONS: SVhound is a unique method to identify potential regions that harbor hidden diversity in model and non model organisms and can also be potentially used to ensure high quality of SV call sets.

Examining the Effects of Hibernation on Germline Mutation Rates in Grizzly Bears.

A male mutation bias is observed across vertebrates, and, where data are available, this bias is accompanied by increased per-generation mutation rates with parental age. While continuing mitotic cell division in the male germline post puberty has been proposed as the major cellular mechanism underlying both patterns, little direct evidence for this role has been found. Understanding the evolution of the per-generation mutation rate among species requires that we identify the molecular mechanisms that change between species. Here, we study the per-generation mutation rate in an extended pedigree of the brown (grizzly) bear, Ursus arctos horribilis. Brown bears hibernate for one-third of the year, a period during which spermatogenesis slows or stops altogether. The reduction of spermatogenesis is predicted to lessen the male mutation bias and to lower the per-generation mutation rate in this species. However, using whole-genome sequencing, we find that both male bias and per-generation mutation rates are highly similar to that expected for a non-hibernating species. We also carry out a phylogenetic comparison of substitution rates along the lineage leading to brown bear and panda (a non-hibernating species) and find no slowing of the substitution rate in the hibernator. Our results contribute to accumulating evidence that suggests that male germline cell division is not the major determinant of mutation rates and mutation biases. The results also provide a quantitative basis for improved estimates of the timing of carnivore evolution.

De novo Mutations in Domestic Cat are Consistent with an Effect of Reproductive Longevity on Both the Rate and Spectrum of Mutations.

The mutation rate is a fundamental evolutionary parameter with direct and appreciable effects on the health and function of individuals. Here, we examine this important parameter in the domestic cat, a beloved companion animal as well as a valuable biomedical model. We estimate a mutation rate of 0.86 × 10-8 per bp per generation for the domestic cat (at an average parental age of 3.8 years). We find evidence for a significant paternal age effect, with more mutations transmitted by older sires. Our analyses suggest that the cat and the human have accrued similar numbers of mutations in the germline before reaching sexual maturity. The per-generation mutation rate in the cat is 28% lower than what has been observed in humans, but is consistent with the shorter generation time in the cat. Using a model of reproductive longevity, which takes into account differences in the reproductive age and time to sexual maturity, we are able to explain much of the difference in per-generation rates between species. We further apply our reproductive longevity model in a novel analysis of mutation spectra and find that the spectrum for the cat resembles the human mutation spectrum at a younger age of reproduction. Together, these results implicate changes in life-history as a driver of mutation rate evolution between species. As the first direct observation of the paternal age effect outside of rodents and primates, our results also suggest a phenomenon that may be universal among mammals.

Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus.

BACKGROUND: The Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity. FINDINGS: Here, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gb, similar to the 2.50-Gb length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity, with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein-coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein-coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where ∼17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0, in which the number of unresolved bases is reduced to 3.00%. CONCLUSIONS: Access to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models.

Comparative molecular genomic analyses of a spontaneous rhesus macaque model of mismatch repair-deficient colorectal cancer.

Colorectal cancer (CRC) remains the third most common cancer in the US with 15% of cases displaying Microsatellite Instability (MSI) secondary to Lynch Syndrome (LS) or somatic hypermethylation of the MLH1 promoter. A cohort of rhesus macaques from our institution developed spontaneous mismatch repair deficient (MMRd) CRC with a notable fraction harboring a pathogenic germline mutation in MLH1 (c.1029C

Heritability of social behavioral phenotypes and preliminary associations with autism spectrum disorder risk genes in rhesus macaques: A whole exome sequencing study.

Nonhuman primates and especially rhesus macaques (Macaca mulatta) have been indispensable animal models for studies of various aspects of neurobiology, developmental psychology, and other aspects of neuroscience. While remarkable progress has been made in our understanding of influences on atypical human social behavior, such as that observed in autism spectrum disorders (ASD), many significant questions remain. Improved understanding of the relationships among variation in specific genes and variation in expressed social behavior in a nonhuman primate would benefit efforts to investigate risk factors, developmental mechanisms, and potential therapies for behavioral disorders including ASD. To study genetic influences on key aspects of social behavior and interactions-individual competence and/or motivation for specific aspects of social behavior-we quantified individual variation in social interactions among juvenile rhesus macaques using both a standard macaque ethogram and a macaque-relevant modification of the human Social Responsiveness Scale. Our analyses demonstrate that various aspects of juvenile social behavior exhibit significant genetic heritability, with estimated quantitative genetic effects similar to that described for ASD in human children. We also performed exome sequencing and analyzed variants in 143 genes previously suggested to influence risk for human ASD. We find preliminary evidence for genetic association between specific variants and both individual behaviors and multi-behavioral factor scores. To our knowledge, this is the first demonstration that spontaneous social behaviors performed by free-ranging juvenile rhesus macaques display significant genetic heritability and then to use exome sequencing data to examine potential macaque genetic associations in genes associated with human ASD.

SULT1A1 and SULT1A2 Associated Extensive Prolapse-Type Inflammatory Polyposis in Crohn's Colitis.

Mucosal prolapse syndrome most commonly involves the rectum and presents as solitary rectal ulcer syndrome and proctitis cystica profunda. Symptoms and endoscopic appearances are nonspecific. Histologically, mucosal prolapse is characterized by fibromuscular obliteration of the lamina propria, and displacement of crypts into submucosa and muscularis mucosae. Mucosal prolapse presenting as polyposis is rare and has only been reported involving the rectosigmoid colon. In this report, we describe a case of mucosal prolapse syndrome presenting as diffuse polyposis and colitis cystica profunda involving the hepatic, splenic flexures and descending colon in a teenage boy suffering from refractory fibrostenosing Crohn's disease. This patient was found to have possibly deleterious homozygous single nucleotide polymorphisms in both SULT1A1 and SULT1A2 genes within a unique polygenic variation of altered cell adhesion.