mGAP: the macaque genotype and phenotype resource, a framework for accessing and interpreting macaque variant data, and identifying new models of human disease

Evolutionary constrained genes associated with autism spectrum disorder across 2,054 nonhuman primate genomes

BACKGROUND

Significant progress has been made in elucidating the genetic underpinnings of Autism Spectrum Disorder (ASD). However, there are still significant gaps in our understanding of the link between genomics, neurobiology and clinical phenotype in scientific discovery. New models are therefore needed to address these gaps. Rhesus macaques (Macaca mulatta) have been extensively used for preclinical neurobiological research because of remarkable similarities to humans across biology and behaviour that cannot be captured by other experimental animals.

METHODS

We used the macaque Genotype and Phenotype (mGAP) resource consisting of 2,054 macaque genomes to examine patterns of evolutionary constraint in known human neurodevelopmental genes. Residual variation intolerance scores (RVIS) were calculated for all annotated autosomal genes (N = 18,168) and Gene Set Enrichment Analysis (GSEA) was used to examine patterns of constraint across ASD genes and related neurodevelopmental genes.

RESULTS

We demonstrated that patterns of constraint across autosomal genes are correlated in humans and macaques, and that ASD-associated genes exhibit significant constraint in macaques (p = 9.4 × 10- 27). Among macaques, many key ASD-implicated genes were observed to harbour predicted damaging mutations. A small number of key ASD-implicated genes that are highly intolerant to mutation in humans, however, showed no evidence of similar intolerance in macaques (CACNA1D, MBD5, AUTS2 and NRXN1). Constraint was also observed across genes associated with intellectual disability (p = 1.1 × 10- 46), epilepsy (p = 2.1 × 10- 33) and schizophrenia (p = 4.2 × 10- 45), and for an overlapping neurodevelopmental gene set (p = 4.0 × 10- 10).

LIMITATIONS

The lack of behavioural phenotypes among the macaques whose genotypes were studied means that we are unable to further investigate whether genetic variants have similar phenotypic consequences among nonhuman primates.

CONCLUSION

The presence of pathological mutations in ASD genes among macaques, along with evidence of similar genetic constraints to those in humans, provides a strong rationale for further investigation of genotype-phenotype relationships in macaques. This highlights the importance of developing primate models of ASD to elucidate the neurobiological underpinnings and advance approaches for precision medicine and therapeutic interventions.

Estimating realized relatedness in free-ranging macaques by inferring identity-by-descent segments

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 identical-by-descent DNA segments (IBD) yield the most precise relatedness estimates. Here, we leverage different methods for estimating IBD segments from low-depth whole genome resequencing data to demonstrate the feasibility and value of resolving fine-scaled gradients of relatedness in free-living animals. Using primarily 4 to 6× depth data from a rhesus macaque (Macaca mulatta) population with long-term pedigree data, we show that we can infer the number and length of IBD segments across the genome with high accuracy even at 0.5× sequencing depth. In line with expectations based on simulation, the resulting estimates demonstrate substantial variation in genetic relatedness within kin classes, leading to overlapping distributions between kin classes. By comparing the IBD-based estimates with pedigree and short tandem repeat-based methods, we show that IBD estimates are more reliable and provide more detailed information on kinship. The inferred IBD segments also identify cryptic genetic relatives not represented in the pedigree and reveal elevated recombination rates in females relative to males, which enables the majority of close maternal and paternal kin to be distinguished with genotype data alone. Our findings represent a breakthrough in the ability to study the predictors and consequences of genetic relatedness in natural populations, contributing to our understanding of a fundamental component of population structure in the wild.

Modifications to rhesus macaque TCR constant regions improve TCR cell surface expression

T cell immunotherapy success is dependent on effective levels of antigen receptor expressed at the surface of engineered cells. Efforts to optimize surface expression in T cell receptor (TCR)-based therapeutic approaches include optimization of cellular engineering methods and coding sequences, and reducing the likelihood of exogenous TCR α and β chains mispairing with the endogenous TCR chains. Approaches to promote correct human TCR chain pairing include constant region mutations to create an additional disulfide bond between the two chains, full murinization of the constant region of the TCR α and β sequences, and a minimal set of murine mutations to the TCR α and β constant regions. Preclinical animal models are valuable tools to optimize engineering designs and methods, and to evaluate the potential for off-target tissue injury. To further develop rhesus macaque models for TCR based cellular immunotherapy, we tested methods for improving cell surface expression of rhesus macaque TCR in rhesus macaque primary cells by generating five alternative TCRαβ constant region constructs in the context of a SIV Gag-specific TCR: 1. human codon optimized rhesus macaque (RH); 2. RH TCR with an additional disulfide linkage; 3. rhesus macaque constant sequences with minimal murine amino acid substitutions; 4. murinized constant sequences; and 5. murinized constant sequences with a portion of the exposed FG loop in the β constant sequence replaced with rhesus macaque sequence to reduce potential immunogencity. Murinization or mutation of a minimal set of amino acids to the corresponding murine sequence of the constant region resulted in the greatest increase in rhesus macaque TCR surface expression relative to wild type. All novel TCR constructs retained the ability to induce production of cytokines in response to cognate peptide antigen specific stimulation. This work can inform the design of TCRs selected for use in rhesus macaque models of TCR-based cellular immunotherapy.

Karl Landsteiner (1868-1943): A Versatile Blood Scientist

Karl Landsteiner (1868-1943) was an Austrian-American biologist, physician, and immunologist whose groundbreaking discoveries revolutionized the fields of transfusion medicine, immunology, and virology. His most famous work was the identification of the ABO blood group system in 1901, which explained the causes of transfusion reactions and laid the foundation for safe blood transfusions. This discovery earned him the Nobel Prize in Physiology or Medicine in 1930. Landsteiner continued his research, identifying the MN and P blood group system in 1927 and the Rh blood group system in 1940, which addressed the complexities of the Hemolytic Disease of the Fetus and Newborn (HDFN). His work on the poliovirus in 1908, in collaboration with Erwin Popper, established the infectious nature of the disease and laid the groundwork for future vaccine development. Landsteiner's versatility as a scientist is evident in the breadth of his work, spanning hematology, immunology, and virology. His ability to navigate diverse fields and see connections between them allowed him to make pioneering discoveries that have had a lasting impact on medical practice, particularly in blood transfusion, organ transplantation, and immunotherapy. Landsteiner's legacy as the "Father of Transfusion Medicine" is reinforced by the standardization of blood typing procedures, which have saved millions of lives worldwide.

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.

Adaptive functions of structural variants in human brain development

Quantifying the structural variants (SVs) in nonhuman primates could provide a niche to clarify the genetic backgrounds underlying human-specific traits, but such resource is largely lacking. Here, we report an accurate SV map in a population of 562 rhesus macaques, verified by in-house benchmarks of eight macaque genomes with long-read sequencing and another one with genome assembly. This map indicates stronger selective constrains on inversions at regulatory regions, suggesting a strategy for prioritizing them with the most important functions. Accordingly, we identified 75 human-specific inversions and prioritized them. The top-ranked inversions have substantially shaped the human transcriptome, through their dual effects of reconfiguring the ancestral genomic architecture and introducing regional mutation hotspots at the inverted regions. As a proof of concept, we linked APCDD1, located on one of these inversions and down-regulated specifically in humans, to neuronal maturation and cognitive ability. We thus highlight inversions in shaping the human uniqueness in brain development.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BACKGROUND

Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.

METHODS

We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.

RESULTS

Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.

CONCLUSIONS

These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.

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.

A naturally occurring variant of MBD4 causes maternal germline hypermutation in primates

As part of an ongoing genome sequencing project at the Oregon National Primate Research Center, we identified a rhesus macaque with a rare homozygous frameshift mutation in the gene methyl-CpG binding domain 4, DNA glycosylase (MBD4). MBD4 is responsible for the repair of C > T deamination mutations at CpG dinucleotides and has been linked to somatic hypermutation and cancer predisposition in humans. We show here that MBD4-associated hypermutation also affects the germline: The six offspring of the MBD4-null dam have a fourfold to sixfold increase in de novo mutation burden. This excess burden was predominantly C > T mutations at CpG dinucleotides consistent with MBD4 loss of function in the dam. There was also a significant excess of C > T at CpA sites, indicating an important, unappreciated role for MBD4 to repair deamination in CpA contexts. The MBD4-null dam developed sustained eosinophilia later in life, but we saw no other signs of neoplastic processes associated with MBD4 loss of function in humans nor any obvious disease in the hypermutated offspring. This work provides the first evidence for a genetic factor causing hypermutation in the maternal germline of a mammal and adds to the very small list of naturally occurring variants known to modulate germline mutation rates in mammals.

Landscape genomics analysis provides insights into future climate change-driven risk in rhesus macaque

Climate change significantly affects the suitability of wildlife habitats. Thus, understanding how animals adapt ecologically and genetically to climate change is important for targeted species protection. Rhesus macaques (Macaca mulatta) are widely distributed and multi-climatically adapted primates. This study explored how rhesus macaques adapt to climate change by integrating ecological and genetic methods and applying species distribution models (SDMs) and a gradient forest (GF) model. The findings suggested that temperature seasonality primarily affects habitat suitability and indicated that climate change will have a dramatic impact on macaque populations in the future. We also applied genotype-environment association (GEA) analyses and selection signature analyses to identify genes associated with climate change and provide possible explanations for the adaptation of rhesus macaques to climatic environments. The population genomics analyses suggested that the Taihang population has the highest genomic vulnerability with inbreeding and low heterozygosity. Combined with the higher ecological vulnerability, additional conservation strategies are required for this population under higher risk of climate change. Our work measured the impact of climate change and enabled the identification of populations that exhibit high vulnerability to severe climate change. Such information is useful for selecting populations of rhesus macaques as subject of long-term monitoring or evolutionary rescue under future climate change.

Genetic variants in melanogenesis proteins TYRP1 and TYR are associated with the golden rhesus macaque phenotype

Nonhuman primates (NHPs) are vital translational research models due to their high genetic, physiological, and anatomical homology with humans. The "golden" rhesus macaque (Macaca mulatta) phenotype is a naturally occurring, inherited trait with a visually distinct pigmentation pattern resulting in light blonde colored fur. Retinal imaging also reveals consistent hypopigmentation and occasional foveal hypoplasia. Here, we describe the use of genome-wide association in 2 distinct NHP populations to identify candidate variants in genes linked to the golden phenotype. Two missense variants were identified in the Tyrosinase-related protein 1 gene (Asp343Gly and Leu415Pro) that segregate with the phenotype. An additional and distinct association was also found with a Tyrosinase variant (His256Gln), indicating the light-colored fur phenotype can result from multiple genetic mechanisms. The implicated genes are related through their contribution to the melanogenesis pathway. Variants in these 2 genes are known to cause pigmentation phenotypes in other species and to be associated with oculocutaneous albinism in humans. The novel associations presented in this study will permit further investigations into the role these proteins and variants play in the melanogenesis pathway and model the effects of genetic hypopigmentation and altered melanogenesis in a naturally occurring nonhuman primate model.

Prenatal delta-9-tetrahydrocannabinol exposure is associated with changes in rhesus macaque DNA methylation enriched for autism genes

BACKGROUND

With the growing availability of cannabis and the popularization of additional routes of cannabis use beyond smoking, including edibles, the prevalence of cannabis use in pregnancy is rapidly increasing. However, the potential effects of prenatal cannabis use on fetal developmental programming remain unknown.

RESULTS

We designed this study to determine whether the use of edible cannabis during pregnancy is deleterious to the fetal and placental epigenome. Pregnant rhesus macaques consumed a daily edible containing either delta-9-tetrahydrocannabinol (THC) (2.5 mg/7 kg/day) or placebo. DNA methylation was measured in 5 tissues collected at cesarean delivery (placenta, lung, cerebellum, prefrontal cortex, and right ventricle of the heart) using the Illumina MethylationEPIC platform and filtering for probes previously validated in rhesus macaque. In utero exposure to THC was associated with differential methylation at 581 CpGs, with 573 (98%) identified in placenta. Loci differentially methylated with THC were enriched for candidate autism spectrum disorder (ASD) genes from the Simons Foundation Autism Research Initiative (SFARI) database in all tissues. The placenta demonstrated greatest SFARI gene enrichment, including genes differentially methylated in placentas from a prospective ASD study.

CONCLUSIONS

Overall, our findings reveal that prenatal THC exposure alters placental and fetal DNA methylation at genes involved in neurobehavioral development that may influence longer-term offspring outcomes. The data from this study add to the limited existing literature to help guide patient counseling and public health polices focused on prenatal cannabis use in the future.

Nonhuman primate genetic models for the study of rare diseases

Pre-clinical research and development relies heavily upon translationally valid models of disease. A major difficulty in understanding the biology of, and developing treatments for, rare disease is the lack of animal models. It is important that these models not only recapitulate the presentation of the disease in humans, but also that they share functionally equivalent underlying genetic causes. Nonhuman primates share physiological, anatomical, and behavioral similarities with humans resulting from close evolutionary relationships and high genetic homology. As the post-genomic era develops and next generation sequencing allows for the resequencing and screening of large populations of research animals, naturally occurring genetic variation in nonhuman primates with clinically relevant phenotypes is regularly emerging. Here we review nonhuman primate models of multiple rare genetic diseases with a focus on the similarities and differences in manifestation and etiologies across species. We discuss how these models are being developed and how they can offer new tools and opportunities for researchers interested in exploring novel therapeutics for these and other genetic diseases. Modeling human genetic diseases in translationally relevant nonhuman primates presents new prospects for development of therapeutics and a better understanding of rare diseases. The post-genomic era offers the opportunity for the discovery and further development of more models like those discussed here.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Cyclopia in a newborn rhesus macaque born to a dam infected with SIV and receiving antiretroviral therapy during pregnancy

Cyclopia, a rare genetic anomaly and birth defect, was recently observed in our nonhuman primate study. A newborn rhesus macaque, delivered via cesarean section, exhibited facial abnormalities, including a single eye in the middle of the forehead. This macaque was born to a dam who had been inoculated with SIV in the first trimester and received antiretroviral therapy (ART) in the early third trimester of pregnancy. Prenatal ultrasound detected fetal defects, including the fusion of the thalami and absence of third ventricle during the third trimester of fetal development. Remarkably, the newborn macaque was diagnosed with severe alobar holoprosencephaly, characterized by a single eye located on the facial midline and proboscises positioned above and below the eye. This condition was accompanied by the absence of a nose, mouth, mandible, maxilla, nasal and oral cavities, tongue, as well as the esophagus. Subsequent genetic screening identified a significant down-regulation of craniofacial development-associated genes, although genetic mutations in the sonic hedgehog gene (SHH) were not present. As the fetal defects were identified prior to the initiation of antiretroviral therapy, it is possible that other environmental factors may have contributed to the development of cyclopia in this rhesus case. However, the etiology of this congenital HPE case remains essentially unknown.

Spontaneous, naturally occurring cancers in non-human primates as a translational model for cancer immunotherapy

The complexity of cancer immunotherapy (CIT) demands reliable preclinical models to successfully translate study findings to the clinics. Non-human primates (NHPs; here referring to rhesus and cynomolgus macaques) share broad similarities with humans including physiology, genetic homology, and importantly also immune cell populations, immune regulatory mechanisms, and protein targets for CIT. Furthermore, NHP naturally develop cancers such as colorectal and breast cancer with an incidence, pathology, and age pattern comparable to humans. Thus, these tumor-bearing monkeys (TBMs) have the potential to bridge the experimental gap between early preclinical cancer models and patients with human cancer.This review presents our current knowledge of NHP immunology, the incidence and features of naturally-occurring cancers in NHP, and recent TBM trials investigating CIT to provide a scientific rationale for this unique model for human cancer.

Genome-Wide DNA Methylation Profile Indicates Potential Epigenetic Regulation of Aging in the Rhesus Macaque Thymus

We analyzed whole-genome bisulfite sequencing (WGBS) and RNA sequencing data of two young (1 year old) and two adult (9 years old) rhesus macaques (Macaca mulatta) to characterize the genomic DNA methylation profile of the thymus and explore the molecular mechanism of age-related changes in the thymus. Combining the two-omics data, we identified correlations between DNA methylation and gene expression and found that DNA methylation played an essential role in the functional changes of the aging thymus, especially in immunity and coagulation. The hypomethylation levels of C3 and C5AR2 and the hypermethylation level of C7 may lead to the high expressions of these genes in adult rhesus macaque thymuses, thus activating the classical complement pathway and the alternative pathway and enhancing their innate immune function. Adult thymuses had an enhanced coagulation pathway, which may have resulted from the hypomethylation and upregulated expressions of seven coagulation-promoting factor genes (F13A1, CLEC4D, CLEC4E, FCN3, PDGFRA, FGF2 and FGF7) and the hypomethylation and low expression of CPB2 to inhibit the degradation of blood clots. Furthermore, the functional decline in differentiation, activation and maturation of T cells in adult thymuses was also closely related to the changes in methylation levels and gene expression levels of T cell development genes (CD3G, GAD2, ADAMDEC1 and LCK) and the thymogenic hormone gene TMPO. A comparison of the age-related methylated genes among four mammal species revealed that most of the epigenetic clocks were species-specific. Furthermore, based on the genomic landscape of allele-specific DNA methylation, we identified several age-related clustered sequence-dependent allele-specific DNA methylated (cS-ASM) genes. Overall, these DNA methylation patterns may also help to assist with understanding the mechanisms of the aging thymus with the epigenome.

Oxytocin and vasotocin receptor variation and the evolution of human prosociality

Modern human lifestyle strongly depends on complex social traits like empathy, tolerance and cooperation. These diverse facets of social cognition have been associated with variation in the oxytocin receptor (OTR) and its sister genes, the vasotocin/vasopressin receptors (VTR1A/AVPR1A and AVPR1B/VTR1B). Here, we compared the available genomic sequences of these receptors between modern humans, archaic humans, and 12 non-human primate species, and identified sites that show heterozygous variation in modern humans and archaic humans distinct from variation in other primates, and for which we could find association studies with clinical implications. On these sites, we performed a range of analyses (variant clustering, pathogenicity prediction, regulation, linkage disequilibrium frequency), and reviewed the literature on selection data in different modern-human populations. We found five sites with modern human specific variation, where the modern human allele is the major allele in the global population (OTR: rs1042778, rs237885, rs6770632; VTR1A: rs10877969; VTR1B: rs33985287). Among them, variation in the OTR-rs6770632 site was predicted to be the most functional. Two alleles (OTR: rs59190448 and rs237888) present only in modern humans and archaic humans were putatively under positive selection in modern humans, with rs237888 predicted to be a highly functional site. Three sites showed convergent evolution between modern humans and bonobos (OTR: rs2228485 and rs237897; VTR1A: rs1042615), with OTR-rs2228485 ranking highly in terms of functionality and reported to be under balancing selection in modern humans (Schaschl, 2015) [1]. Our findings have implications for understanding hominid prosociality, as well as the similarities between modern human and bonobo social behavior.

Genomic resources for rhesus macaques (Macaca mulatta)

Rhesus macaques (Macaca mulatta) are among the most extensively studied of nonhuman primates. This species has been the subject of many investigations concerning basic primate biology and behavior, including studies of social organization, developmental psychology, physiology, endocrinology, and neurodevelopment. Rhesus macaques are also critically important as a nonhuman primate model of human health and disease, including use in studies of infectious diseases, metabolic diseases, aging, and drug or alcohol abuse. Current research addressing fundamental biological and/or applied biomedical questions benefits from various genetic and genomic analyses. As a result, the genome of rhesus macaques has been the subject of more study than most nonhuman primates. This paper briefly discusses a number of information resources that can provide interested researchers with access to genetic and genomic data describing the content of the rhesus macaque genome, available information regarding genetic variation within the species, results from studies of gene expression, and other aspects of genomic analysis. Specific online databases are discussed, including the US National Center for Biotechnology Information, the University of California Santa Cruz genome browser, Ensembl genome browser, the Macaque Genotype and Phenotype database (mGAP), Rhesusbase, and others.

Improving rigor and reproducibility in nonhuman primate research

Nonhuman primates (NHPs) are a critical component of translational/preclinical biomedical research due to the strong similarities between NHP and human physiology and disease pathology. In some cases, NHPs represent the most appropriate, or even the only, animal model for complex metabolic, neurological, and infectious diseases. The increased demand for and limited availability of these valuable research subjects requires that rigor and reproducibility be a prime consideration to ensure the maximal utility of this scarce resource. Here, we discuss a number of approaches that collectively can contribute to enhanced rigor and reproducibility in NHP research.

A novel non-human primate model of Pelizaeus-Merzbacher disease

Pelizaeus-Merzbacher disease (PMD) is a severe hypomyelinating disorder of the central nervous system (CNS) linked to mutations in the proteolipid protein-1 (PLP1) gene. Although there are multiple animal models of PMD, few of them fully mimic the human disease. Here, we report three spontaneous cases of male neonatal rhesus macaques with the clinical symptoms of hypomyelinating disease, including intention tremors, progressively worsening motor dysfunction, and nystagmus. These animals demonstrated a paucity of CNS myelination accompanied by reactive astrogliosis, and a lack of PLP1 expression throughout white matter. Genetic analysis revealed that these animals were related to one another and that their parents carried a rare, hemizygous missense variant in exon 5 of the PLP1 gene. These animals therefore represent the first reported non-human primate model of PMD, providing a novel and valuable opportunity for preclinical studies that aim to promote myelination in pediatric hypomyelinating diseases.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases

PURPOSE

Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful.

METHODS

We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols.

RESULTS

We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases.

CONCLUSION

The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

MACSNVdb: a high-quality SNV database for interspecies genetic divergence investigation among macaques

Macaques are the most widely used non-human primates in biomedical research. The genetic divergence between these animal models is responsible for their phenotypic differences in response to certain diseases. However, the macaque single nucleotide polymorphism resources mainly focused on rhesus macaque (Macaca mulatta), which hinders the broad research and biomedical application of other macaques. In order to overcome these limitations, we constructed a database named MACSNVdb that focuses on the interspecies genetic diversity among macaque genomes. MACSNVdb is a web-enabled database comprising ~74.51 million high-quality non-redundant single nucleotide variants (SNVs) identified among 20 macaque individuals from six species groups (muttla, fascicularis, sinica, arctoides, silenus, sylvanus). In addition to individual SNVs, MACSNVdb also allows users to browse and retrieve groups of user-defined SNVs. In particular, users can retrieve non-synonymous SNVs that may have deleterious effects on protein structure or function within macaque orthologs of human disease and drug-target genes. Besides position, alleles and flanking sequences, MACSNVdb integrated additional genomic information including SNV annotations and gene functional annotations. MACSNVdb will facilitate biomedical researchers to discover molecular mechanisms of diverse responses to diseases as well as primatologist to perform population genetic studies. We will continue updating MACSNVdb with newly available sequencing data and annotation to keep the resource up to date.

Database URL: http://big.cdu.edu.cn/macsnvdb/

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

Naturally Occurring and Experimentally Induced Rhesus Macaque Models for Polycystic Ovary Syndrome: Translational Gateways to Clinical Application

Indian rhesus macaque nonhuman primate models for polycystic ovary syndrome (PCOS) implicate both female hyperandrogenism and developmental molecular origins as core components of PCOS etiopathogenesis. Establishing and exploiting macaque models for translational impact into the clinic, however, has required multi-year, integrated basic-clinical science collaborations. Paradigm shifting insight has accrued from such concerted investment, leading to novel mechanistic understanding of PCOS, including hyperandrogenic fetal and peripubertal origins, epigenetic programming, altered neural function, defective oocytes and embryos, adipogenic constraint enhancing progression to insulin resistance, pancreatic decompensation and type 2 diabetes, together with placental compromise, all contributing to transgenerational transmission of traits likely to manifest in adult PCOS phenotypes. Our recent demonstration of PCOS-related traits in naturally hyperandrogenic (High T) female macaques additionally creates opportunities to employ whole genome sequencing to enable exploration of gene variants within human PCOS candidate genes contributing to PCOS-related traits in macaque models. This review will therefore consider Indian macaque model contributions to various aspects of PCOS-related pathophysiology, as well as the benefits of using macaque models with compellingly close homologies to the human genome, phenotype, development and aging.

Bardet-Biedl Syndrome in rhesus macaques: A nonhuman primate model of retinitis pigmentosa

The development of therapies for retinal disorders is hampered by a lack of appropriate animal models. Higher nonhuman primates are the only animals with retinal structure similar to humans, including the presence of a macula and fovea. However, few nonhuman primate models of genetic retinal disease are known. We identified a lineage of rhesus macaques with a frameshift mutation in exon 3 of the BBS7 gene c.160delG (p.Ala54fs) that is predicted to produce a non-functional protein. In humans, mutations in this and other BBS genes cause Bardet-Biedl syndrome, a ciliopathy and a syndromic form of retinitis pigmentosa generally occurring in conjunction with kidney dysfunction, polydactyly, obesity, and/or hypogonadism. Three full- or half-sibling monkeys homozygous for the BBS7 c.160delG variant, at ages 3.5, 4 and 6 years old, displayed a combination of severe photoreceptor degeneration and progressive kidney disease. In vivo retinal imaging revealed features of severe macular degeneration, including absence of photoreceptor layers, degeneration of the retinal pigment epithelium, and retinal vasculature atrophy. Electroretinography in the 3.5-year-old case demonstrated loss of scotopic and photopic a-waves and markedly reduced and delayed b-waves. Histological assessments in the 4- and 6-year-old cases confirmed profound loss of photoreceptors and inner retinal neurons across the posterior retina, with dramatic thinning and disorganization of all cell layers, abundant microglia, absent or displaced RPE cells, and significant gliosis in the subretinal space. Retinal structure, including presence of photoreceptors, was preserved only in the far periphery. Ultrasound imaging of the kidneys revealed deranged architecture, and renal histopathology identified distorted contours with depressed, fibrotic foci and firmly adhered renal capsules; renal failure occurred in the 6-year-old case. Magnetic resonance imaging obtained in one case revealed abnormally low total brain volume and unilateral ventricular enlargement. The one male had abnormally small testes at 4 years of age, but polydactyly and obesity were not observed. Thus, monkeys homozygous for the BBS7 c.160delG variant closely mirrored several key features of the human BBS syndrome. This finding represents the first identification of a naturally-occurring nonhuman primate model of BBS, and more broadly the first such model of retinitis pigmentosa and a ciliopathy with an associated genetic mutation. This important new preclinical model will provide the basis for better understanding of disease progression and for the testing of new therapeutic options, including gene and cell-based therapies, not only for BBS but also for multiple forms of photoreceptor degeneration.