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.

TiO2 nanorod arrays@PDA/Ag with biomimetic polydopamine as binary mediators for duplex SERS detection of illegal food dyes

Based on TiO2 nanorod arrays@PDA/Ag (TNRs@PDA/Ag), a better surface-enhanced Raman scattering (SERS) sensor with effective enrichment and enhancement was investigated for duplex SERS detection of illicit food dyes. Biomimetic PDA functions as binary mediators by utilizing the structural characteristics of polydopamine (PDA), which include the conjugated structure and abundant hydrophilic groups. One PDA functioned as an electron transfer mediator to enhance the efficiency of electron transfer, and the other as an enrichment mediator to effectively enrich rhodamine B (RhB) and crystal violet (CV) through hydrogen bonding, π-π stacking, and electrostatic interactions. Individual and duplex detection of illicit food dyes (RhB and CV) was performed using TNRs@PDA/Ag to estimate SERS applications. Their linear equations and limits of detection of 1 nM for RhB and 5 nM for CV were derived. Individual and duplex food colour detection was successfully accomplished even in genuine chili meal with good results. The bifunctional TNRs@PDA/Ag-based highly sensitive and duplex SERS dye detection will have enormous potential for food safety monitoring.

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.

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.

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.

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.

Loci on chromosome 12q13.2 encompassing ERBB3, PA2G4 and RAB5B are associated with polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is characterized by hyperandrogenemia of ovarian theca cell origin. We report significant association of androgen production with 15 single nucleotide variants (SNVs) identified by exome sequencing of theca cells from women with PCOS and normal ovulatory women. Ten SNVs are located within a 150 kbp region on 12q13.2 which encompasses loci identified in PCOS genome-wide association studies (GWAS) and contains PCOS candidate genes ERBB3 and RAB5B. The region also contains PA2G4 which encodes a transcriptional corepressor of androgen receptor and androgen receptor-regulated genes. PA2G4 has not previously been recognized as related to PCOS in published GWAS studies. Two of the SNVs are predicted to have functional consequences (ERBB3 missense SNV, PA2G4 promoter SNV). PA2G4 interacts with the ERBB3 cytoplasmic domain containing the missense variant, suggesting a potential signaling pathway disruption that could lead to the PCOS ovarian phenotype. Single cell RNA sequencing of theca cells showed significantly less expression of PA2G4 after forskolin treatment in PCOS cells compared to normal cells (padj = 3.82E-30) and in cells heterozygous for the PA2G4 promoter SNV compared to those without the SNV (padj = 2.16E-11). This is consistent with a functional effect of the PA2G4 promoter SNV. No individual SNV was significantly associated with PCOS in an independent family cohort, but a haplotype with minor alleles of three SNVs was found preferentially in women with PCOS. These findings suggest a functional role for 12q13.2 variants in PCOS and implicate variants in ERBB3 and PA2G4 in the pathophysiology of PCOS.

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 Collapse

Key Words Collapse

MESH Headings

• Animals
• Colonic Neoplasms/genetics
• Colorectal Neoplasms/genetics
• Colorectal Neoplasms/pathology
• Colorectal Neoplasms, Hereditary Nonpolyposis/genetics
• Colorectal Neoplasms, Hereditary Nonpolyposis/pathology
• DNA Methylation/genetics
• DNA Mismatch Repair/genetics
• Genomics
• Humans
• Macaca mulatta/genetics
• Mice
• Microsatellite Instability
• MutL Protein Homolog 1/genetics

Collapse

Grants

• P30 CA016672 NCI NIH HHS
• R24 OD011173 NIH HHS
• Feinberg Family Foundation
• MDACC Institutional Research Grant Program
• Cattleman for Cancer Research
• National Institutes of Health
• National Cancer Institute

Collapse

Affiliation(s)

Nejla Ozirmak Lermi Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America School of Health Professions, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Stanton B. Gray Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Charles M. Bowen Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Laura Reyes-Uribe Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Beth K. Dray Charles River Laboratories, Ashland, Ohio, United States of America
Nan Deng Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
R. Alan Harris Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Muthuswamy Raveendran Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Fernando Benavides Department of Epigenetics and Molecular Carcinogenesis, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Carolyn L. Hodo Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Melissa W. Taggart Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Karen Colbert Maresso Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Krishna M. Sinha Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Jeffrey Rogers Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Eduardo Vilar Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America Clinical Cancer Genetics Program, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America

Collapse

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.

Multiplexed Simian Immunodeficiency Virus-Specific Paired RNA-Guided Cas9 Nickases Inactivate Proviral DNA

Human and simian immunodeficiency virus (HIV and SIV) infections establish lifelong reservoirs of cells harboring an integrated proviral genome. Genome editing CRISPR-associated Cas9 nucleases, combined with SIV-specific guiding RNA (gRNA) molecules, inactivate integrated provirus DNA in vitro and in animal models. We generated RNA-guided Cas9 nucleases (RGNu) and nickases (RGNi) targeting conserved SIV regions with no homology in the human or rhesus macaque genome. Assays in cells cotransfected with SIV provirus and plasmids coding for RGNus identified SIV long terminal repeat (LTR), trans-activation response (TAR) element, and ribosome slip site (RSS) regions as the most effective at virus suppression; RGNi targeting these regions inhibited virus production significantly. Multiplex plasmids that coexpressed these three RGNu (Nu3), or six (three pairs) RGNi (Ni6), were more efficient at virus suppression than any combination of individual RGNu and RGNi plasmids. Both Nu3 and Ni6 plasmids were tested in lymphoid cells chronically infected with SIVmac239, and whole-genome sequencing was used to determine on- and off-target mutations. Treatment with these all-in-one plasmids resulted in similar levels of mutations of viral sequences from the cellular genome; Nu3 induced indels at the 3 SIV-specific sites, whereas for Ni6 indels were present at the LTR and TAR sites. Levels of off-target effects detected by two different algorithms were indistinguishable from background mutations. In summary, we demonstrate that Cas9 nickase in association with gRNA pairs can specifically eliminate parts of the integrated provirus DNA; also, we show that careful design of an all-in-one plasmid coding for 3 gRNAs and Cas9 nuclease inhibits SIV production with undetectable off-target mutations, making these tools a desirable prospect for moving into animal studies. IMPORTANCE Our approach to HIV cure, utilizing the translatable SIV/rhesus macaque model system, aims at provirus inactivation and its removal with the least possible off-target side effects. We developed single molecules that delivered either three truncated SIV-specific gRNAs along with Cas9 nuclease or three pairs of SIV-specific gRNAs (six individual gRNAs) along with Cas9 nickase to enhance efficacy of on-target mutagenesis. Whole-genome sequencing demonstrated effective SIV sequence mutation and inactivation and the absence of demonstrable off-target mutations. These results open the possibility to employ Cas9 variants that introduce single-strand DNA breaks to eliminate integrated proviral DNA.

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.

Origins and Long-Term Patterns of Copy-Number Variation in Rhesus Macaques

Mutations play a key role in the development of disease in an individual and the evolution of traits within species. Recent work in humans and other primates has clarified the origins and patterns of single-nucleotide variants, showing that most arise in the father’s germline during spermatogenesis. It remains unknown whether larger mutations, such as deletions and duplications of hundreds or thousands of nucleotides, follow similar patterns. Such mutations lead to copy-number variation (CNV) within and between species, and can have profound effects by deleting or duplicating genes. Here, we analyze patterns of CNV mutations in 32 rhesus macaque individuals from 14 parent–offspring trios. We find the rate of CNV mutations per generation is low (less than one per genome) and we observe no correlation between parental age and the number of CNVs that are passed on to offspring. We also examine segregating CNVs within the rhesus macaque sample and compare them to a similar data set from humans, finding that both species have far more segregating deletions than duplications. We contrast this with long-term patterns of gene copy-number evolution between 17 mammals, where the proportion of deletions that become fixed along the macaque lineage is much smaller than the proportion of segregating deletions. These results suggest purifying selection acting on deletions, such that the majority of them are removed from the population over time. Rhesus macaques are an important biomedical model organism, so these results will aid in our understanding of this species and the disease models it supports.

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.

Comparative genomic analysis of sifakas (Propithecus) reveals selection for folivory and high heterozygosity despite endangered status

Sifakas (genus Propithecus) are critically endangered, large-bodied diurnal lemurs that eat leaf-based diets and show corresponding anatomical and microbial adaptations to folivory. We report on the genome assembly of Coquerel's sifaka (P. coquereli) and the resequenced genomes of Verreaux's (P. verreauxi), the golden-crowned (P. tattersalli), and the diademed (P. diadema) sifakas. We find high heterozygosity in all sifakas compared with other primates and endangered mammals. Demographic reconstructions nevertheless suggest declines in effective population size beginning before human arrival on Madagascar. Comparative genomic analyses indicate pervasive accelerated evolution in the ancestral sifaka lineage affecting genes in several complementary pathways relevant to folivory, including nutrient absorption and xenobiotic and fatty acid metabolism. Sifakas show convergent evolution at the level of the pathway, gene family, gene, and amino acid substitution with other folivores. Although sifakas have relatively generalized diets, the physiological challenges of habitual folivory likely led to strong selection.

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.

Pedigree reconstruction and distant pairwise relatedness estimation from genome sequence data: A demonstration in a population of rhesus macaques (Macaca mulatta)

A primary challenge in the analysis of free‐ranging animal populations is the accurate estimation of relatedness among individuals. Many aspects of population analysis rely on knowledge of relatedness patterns, including socioecology, demography, heritability and gene mapping analyses, wildlife conservation and the management of breeding colonies. Methods for determining relatedness using genome‐wide data have improved our ability to determine kinship and reconstruct pedigrees in humans. However, methods for reconstructing complex pedigree structures and estimating distant relatedness (beyond third‐degree) have not been widely applied to other species. We sequenced the genomes of 150 male rhesus macaques from the Tulane National Primate Research Center colony to estimate pairwise relatedness, reconstruct closely related pedigrees, estimate more distant relationships and augment colony records. Methods for determining relatedness developed for human genetic data were applied and evaluated in the analysis of nonhuman primates, including identity‐by‐descent‐based methods for pedigree reconstruction and shared segment‐based inference of more distant relatedness. We compared the genotype‐based pedigrees and estimated relationships to available colony pedigree records and found high concordance (95.5% agreement) between expected and identified relationships for close relatives. In addition, we detected distant relationships not captured in colony records, including some as distant as twelfth‐degree. Furthermore, while deep sequence coverage is preferable, we show that this approach can also provide valuable information when only low‐coverage (5×) sequence data is available. Our findings demonstrate the value of these methods for determination of relatedness in various animal populations, with diverse applications to conservation biology, evolutionary and ecological research and biomedical studies.

Sequence diversity analyses of an improved rhesus macaque genome enhance its biomedical utility

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

Role of toll-like receptor 7 (TLR7) in voluntary alcohol consumption

Overactivation of neuroimmune signaling has been linked to excessive ethanol consumption. Toll-like receptors (TLRs) are a major component of innate immune signaling and initiate anti- and pro-inflammatory responses via intracellular signal transduction cascades. TLR7 is upregulated in post-mortem brain tissue from humans with alcohol use disorder (AUD) and animals with prior exposure to ethanol. Despite this evidence, the role of TLR7 in the regulation of voluntary ethanol consumption has not been studied. We test the hypothesis that TLR7 activation regulates voluntary ethanol drinking behavior by administering a TLR7 agonist (R848) during an intermittent access drinking procedure in mice. Acute activation of TLR7 reduced ethanol intake, preference, and total fluid intake due, at least in part, to an acute sickness response. However, chronic pre-treatment with R848 resulted in tolerance to the adverse effects of the drug and a subsequent increase in ethanol consumption. To determine the molecular machinery that mediates these behavioral changes, we evaluated gene expression after acute and chronic TLR7 activation. We found that acute TLR7 activation produces brain region specific changes in expression of immune pathway genes, whereas chronic TLR7 activation causes downregulation of TLRs and blunted cytokine induction, suggesting molecular tolerance. Our results demonstrate a novel role for TLR7 signaling in regulating voluntary ethanol consumption. Taken together, our findings suggest TLR7 may be a viable target for development of therapies to treat AUD.

Reduced meiotic recombination in rhesus macaques and the origin of the human recombination landscape

Characterizing meiotic recombination rates across the genomes of nonhuman primates is important for understanding the genetics of primate populations, performing genetic analyses of phenotypic variation and reconstructing the evolution of human recombination. Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primates in biomedical research. We constructed a high-resolution genetic map of the rhesus genome based on whole genome sequence data from Indian-origin rhesus macaques. The genetic markers used were approximately 18 million SNPs, with marker density 6.93 per kb across the autosomes. We report that the genome-wide recombination rate in rhesus macaques is significantly lower than rates observed in apes or humans, while the distribution of recombination across the macaque genome is more uniform. These observations provide new comparative information regarding the evolution of recombination in primates.

Paternal age in rhesus macaques is positively associated with germline mutation accumulation but not with measures of offspring sociability

Mutation is the ultimate source of all genetic novelty and the cause of heritable genetic disorders. Mutational burden has been linked to complex disease, including neurodevelopmental disorders such as schizophrenia and autism. The rate of mutation is a fundamental genomic parameter and direct estimates of this parameter have been enabled by accurate comparisons of whole-genome sequences between parents and offspring. Studies in humans have revealed that the paternal age at conception explains most of the variation in mutation rate: Each additional year of paternal age in humans leads to approximately 1.5 additional inherited mutations. Here, we present an estimate of the de novo mutation rate in the rhesus macaque (Macaca mulatta) using whole-genome sequence data from 32 individuals in four large pedigrees. We estimated an average mutation rate of 0.58 × 10⁻⁸ per base pair per generation (at an average parental age of 7.5 yr), much lower than found in direct estimates from great apes. As in humans, older macaque fathers transmit more mutations to their offspring, increasing the per generation mutation rate by 4.27 × 10⁻¹⁰ per base pair per year. We found that the rate of mutation accumulation after puberty is similar between macaques and humans, but that a smaller number of mutations accumulate before puberty in macaques. We additionally investigated the role of paternal age on offspring sociability, a proxy for normal neurodevelopment, by studying 203 male macaques in large social groups.

Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta)

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

Unusual sequence characteristics of human chromosome 19 are conserved across 11 nonhuman primates

BACKGROUND

Human chromosome 19 has many unique characteristics including gene density more than double the genome-wide average and 20 large tandemly clustered gene families. It also has the highest GC content of any chromosome, especially outside gene clusters. The high GC content and concomitant high content of hypermutable CpG sites raises the possibility chromosome 19 exhibits higher levels of nucleotide diversity both within and between species, and may possess greater variation in DNA methylation that regulates gene expression.

RESULTS

We examined GC and CpG content of chromosome 19 orthologs across representatives of the primate order. In all 12 primate species with suitable genome assemblies, chromosome 19 orthologs have the highest GC content of any chromosome. CpG dinucleotides and CpG islands are also more prevalent in chromosome 19 orthologs than other chromosomes. GC and CpG content are generally higher outside the gene clusters. Intra-species variation based on SNPs in human common dbSNP, rhesus, crab eating macaque, baboon and marmoset datasets is most prevalent on chromosome 19 and its orthologs. Inter-species comparisons based on phyloP conservation show accelerated nucleotide evolution for chromosome 19 promoter flanking and enhancer regions. These same regulatory regions show the highest CpG density of any chromosome suggesting they possess considerable methylome regulatory potential.

CONCLUSIONS

The pattern of high GC and CpG content in chromosome 19 orthologs, particularly outside gene clusters, is present from human to mouse lemur representing 74 million years of primate evolution. Much CpG variation exists both within and between primate species with a portion of this variation occurring in regulatory regions.

MHC genotyping from rhesus macaque exome sequences

Indian rhesus macaque major histocompatibility complex (MHC) variation can influence the outcomes of transplantation and infectious disease studies. Frequently, rhesus macaques are MHC genotyped to identify variants that could account for unexpected results. Since the MHC is only one region in the genome where variation could impact experimental outcomes, strategies for simultaneously profiling variation in the macaque MHC and the remainder of the protein coding genome would be useful. Here we determine MHC class I and class II genotypes using target-capture probes enriched for MHC sequences, a method we term macaque exome sequence (MES) genotyping. For a cohort of 27 Indian rhesus macaques, we describe two methods for obtaining MHC genotypes from MES data and demonstrate that the MHC class I and class II genotyping results obtained with these methods are 98.1% and 98.7% concordant, respectively, with expected MHC genotypes. In contrast, conventional MHC genotyping results obtained by deep sequencing of short multiplex PCR amplicons were only 92.6% concordant with expectations for this cohort.

Population-Based Estimation of Dental Caries and Periodontal Disease Rates of Gravid and Recently Postpartum Women in Lilongwe, Malawi

Objectives  The objective of this study was to determine the rate of dental caries and periodontal disease among gravid and recently postpartum women at five delivery centers within and surrounding Lilongwe, Malawi. Study Design  We partnered with obstetric specialists, community health workers, and dentists to perform dental history interviews and dental examinations during the study period from December 2012 to May 2014. Dental examinations were performed according to World Health Organization standards to assess periodontal and oral health status. Results  Among the 387 gravid and recently postpartum women, the rate of dental caries was 69.3% and the rate of composite dental disease (caries and periodontal disease) was 76.7%. The majority (69.5%) of women examined had a decayed-missing-filled (DMF) index greater than or equal to one; the average DMF Index was 2.48. The majority of women had never seen a dentist (62.8%). However, most did perform oral hygiene, two or more times per day (90.2%); most women reported brushing with toothpaste (88.1%). Conclusion  When assessing this population for dental caries and periodontal disease, the rate of dental disease was high. Therefore, this may be an ideal setting to test for impactful interventions aimed at reducing caries and periodontal disease.

Associations between muscle mass, physical activity and dietary behaviour in adolescents

BACKGROUND AND OBJECTIVES

Investigations in older individuals suggest that adequate nutrition and physical activity are particularly important to skeletal muscle health; however, data in adolescents are scant and equivocal. The objective was to determine the associations among diet, physical activity and skeletal muscle mass in adolescents.

METHODS

We assessed diet with four to seven 24-h recalls and physical activity by accelerometry in 640 adolescents. Using total body measures of fat-free soft tissue mass and fat mass assessed by dual-energy X-ray absorptiometry, the skeletal muscle mass index (SMMI) was derived by adjusting fat-free soft tissue mass for fat mass in addition to height.

RESULTS

Skeletal muscle mass index (SMMI) was negatively associated with consumption of sugar-sweetened beverages (standardized beta coefficient [β] = -0.10, P = 0.001) and saturated fats (β = -0.28, P < 0.001). SMMI was positively associated with physical activity (moderate + vigorous) (β = 0.20, P < 0.001). In further analysis, we observed a significant interaction between physical activity and sugar-sweetened beverage intake on SMMI (P = 0.002).

CONCLUSION

Our study in adolescence suggests that physical activity and consumption of both sugar-sweetened beverages and saturated fats are associated with skeletal muscle mass. More importantly, our findings suggest that sugar-sweetened beverage intake may attenuate the beneficial effects of physical activity on skeletal muscle mass.

A nonhuman primate model of inherited retinal disease

Inherited retinal degenerations are a common cause of untreatable blindness worldwide, with retinitis pigmentosa and cone dystrophy affecting approximately 1 in 3500 and 1 in 10,000 individuals, respectively. A major limitation to the development of effective therapies is the lack of availability of animal models that fully replicate the human condition. Particularly for cone disorders, rodent, canine, and feline models with no true macula have substantive limitations. By contrast, the cone-rich macula of a nonhuman primate (NHP) closely mirrors that of the human retina. Consequently, well-defined NHP models of heritable retinal diseases, particularly cone disorders that are predictive of human conditions, are necessary to more efficiently advance new therapies for patients. We have identified 4 related NHPs at the California National Primate Research Center with visual impairment and findings from clinical ophthalmic examination, advanced retinal imaging, and electrophysiology consistent with achromatopsia. Genetic sequencing confirmed a homozygous R565Q missense mutation in the catalytic domain of PDE6C, a cone-specific phototransduction enzyme associated with achromatopsia in humans. Biochemical studies demonstrate that the mutant mRNA is translated into a stable protein that displays normal cellular localization but is unable to hydrolyze cyclic GMP (cGMP). This NHP model of a cone disorder will not only serve as a therapeutic testing ground for achromatopsia gene replacement, but also for optimization of gene editing in the macula and of cone cell replacement in general.

The comparative genomics and complex population history of Papio baboons

Recent studies suggest that closely related species can accumulate substantial genetic and phenotypic differences despite ongoing gene flow, thus challenging traditional ideas regarding the genetics of speciation. Baboons (genus Papio) are Old World monkeys consisting of six readily distinguishable species. Baboon species hybridize in the wild, and prior data imply a complex history of differentiation and introgression. We produced a reference genome assembly for the olive baboon (Papio anubis) and whole-genome sequence data for all six extant species. We document multiple episodes of admixture and introgression during the radiation of Papio baboons, thus demonstrating their value as a model of complex evolutionary divergence, hybridization, and reticulation. These results help inform our understanding of similar cases, including modern humans, Neanderthals, Denisovans, and other ancient hominins.

Maresin 1 attenuates neuroinflammation in a mouse model of perioperative neurocognitive disorders

Background

Resolution of inflammation is an active and dynamic process after surgery. Maresin 1 (MaR1) is one of a growing number of specialised pro-resolving lipids biosynthesised by macrophages that regulates acute inflammation. We investigated the effects of MaR1 on postoperative neuroinflammation, macrophage activity, and cognitive function in mice.

Methods

Adult male C57BL/6 (n=111) and Ccr2^RFP/+Cx3cr1^GFP/+ (n=54) mice were treated with MaR1 before undergoing anaesthesia and orthopaedic surgery. Systemic inflammatory changes, bone healing, neuroinflammation, and cognition were assessed at different time points. MaR1 protective effects were also evaluated using bone marrow derived macrophage cultures.

Results

MaR1 exerted potent systemic anti-inflammatory effects without impairing fracture healing. Prophylaxis with MaR1 prevented surgery-induced glial activation and opening of the blood–brain barrier. In Ccr2^RFP/+Cx3cr1^GFP/+ mice, fewer infiltrating macrophages were detected in the hippocampus after surgery with MaR1 prophylaxis, which resulted in improved memory function. MaR1 treatment also reduced expression of pro-inflammatory cell surface markers and cytokines by in vitro cultured macrophages. MaR1 was detectable in the cerebrospinal fluid of older adults before and after surgery.

Conclusions

MaR1 exerts distinct anti-inflammatory and pro-resolving effects through regulation of macrophage infiltration, NF-κB signalling, and cytokine release after surgery. Future studies on the use of pro-resolving lipid mediators may inform novel approaches to treat neuroinflammation and postoperative neurocognitive disorders.

Neuroimmune signaling in alcohol use disorder

Alcohol use disorder (AUD) is a widespread disease with limited treatment options. Targeting the neuroimmune system is a new avenue for developing or repurposing effective pharmacotherapies. Alcohol modulates innate immune signaling in different cell types in the brain by altering gene expression and the molecular pathways that regulate neuroinflammation. Chronic alcohol abuse may cause an imbalance in neuroimmune function, resulting in prolonged perturbations in brain function. Likewise, manipulating the neuroimmune system may change alcohol-related behaviors. Psychiatric disorders that are comorbid with AUD, such as post-traumatic stress disorder, major depressive disorder, and other substance use disorders, may also have underlying neuroimmune mechanisms; current evidence suggests that convergent immune pathways may be involved in AUD and in these comorbid disorders. In this review, we provide an overview of major neuroimmune cell-types and pathways involved in mediating alcohol behaviors, discuss potential mechanisms of alcohol-induced neuroimmune activation, and present recent clinical evidence for candidate immune-related drugs to treat AUD.

Spontaneous Spongiform Brainstem Degeneration in a Young Mouse Lemur (Microcebus murinus) with Conspicuous Behavioral, Motor, Growth, and Ocular Pathologies

Here we report a case of severe growth retardation and neurologic abnormalities in a female gray mouse lemur (Microcebus murinus), a small NHP species for which the genomic sequence recently became available. The female lemur we present here died on postnatal day 125. This lemur had impaired development of motor skills and showed severe ataxia and tremors. In addition, hearing seemed normal whereas ophthalmic examination revealed incipient bilateral cataracts, abnormal pigmentation in the lens of the left eye, and a missing optokinetic nystagmus, which indicated impaired vision. Most prominently, the lemur showed severe growth retardation. Necropsy revealed maldevelopment of the left reproductive organs and unilateral dilation of the right lateral ventricle, which was confirmed on brain MRI. Brain histology further revealed large, bilateral areas of vacuolation within the brainstem, but immunohistochemistry indicated no sign of pathologic prion protein deposition. Full genomic sequencing of the lemur revealed a probably pathologic mutation in LARGE2 of the LARGE gene family, which has been associated with congenital muscular dystrophies. However, potentially functional mutations in other genes were also present. The observed behavioral and motor signs in the presented animal might have been linked to spongiform degeneration and resulting brainstem dysfunction and progressive muscle weakness. The macroscopic developmental abnormalities and ophthalmic findings might be genetic in origin and linked to the mutation in LARGE2.

Allele-specific epigenome maps reveal sequence-dependent stochastic switching at regulatory loci

To assess the impact of genetic variation in regulatory loci on human health, we constructed a high-resolution map of allelic imbalances in DNA methylation, histone marks, and gene transcription in 71 epigenomes from 36 distinct cell and tissue types from 13 donors. Deep whole-genome bisulfite sequencing of 49 methylomes revealed sequence-dependent CpG methylation imbalances at thousands of heterozygous regulatory loci. Such loci are enriched for stochastic switching, which is defined as random transitions between fully methylated and unmethylated states of DNA. The methylation imbalances at thousands of loci are explainable by different relative frequencies of the methylated and unmethylated states for the two alleles. Further analyses provided a unifying model that links sequence-dependent allelic imbalances of the epigenome, stochastic switching at gene regulatory loci, and disease-associated genetic variation.

Reproductive Longevity Predicts Mutation Rates in Primates

Mutation rates vary between species across several orders of magnitude, with larger organisms having the highest per-generation mutation rates. Hypotheses for this pattern typically invoke physiological or population-genetic constraints imposed on the molecular machinery preventing mutations [1]. However, continuing germline cell division in multicellular eukaryotes means that organisms with longer generation times and of larger size will leave more mutations to their offspring simply as a byproduct of their increased lifespan [2, 3]. Here, we deeply sequence the genomes of 30 owl monkeys (Aotus nancymaae) from six multi-generation pedigrees to demonstrate that paternal age is the major factor determining the number of de novo mutations in this species. We find that owl monkeys have an average mutation rate of 0.81 × 10^-8 per site per generation, roughly 32% lower than the estimate in humans. Based on a simple model of reproductive longevity that does not require any changes to the mutational machinery, we show that this is the expected mutation rate in owl monkeys. We further demonstrate that our model predicts species-specific mutation rates in other primates, including study-specific mutation rates in humans based on the average paternal age. Our results suggest that variation in life history traits alone can explain variation in the per-generation mutation rate among primates, and perhaps among a wide range of multicellular organisms.

Sooty mangabey genome sequence provides insight into AIDS resistance in a natural SIV host

In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia¹. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3–4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

Hybrid de novo genome assembly and centromere characterization of the gray mouse lemur (Microcebus murinus)

BACKGROUND

The de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly.

METHODS

We used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome.

RESULTS

We improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies.

CONCLUSIONS

We describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes.

Improved full-length killer cell immunoglobulin-like receptor transcript discovery in Mauritian cynomolgus macaques

Killer cell immunoglobulin-like receptors (KIRs) modulate disease progression of pathogens including HIV, malaria, and hepatitis C. Cynomolgus and rhesus macaques are widely used as nonhuman primate models to study human pathogens, and so, considerable effort has been put into characterizing their KIR genetics. However, previous studies have relied on cDNA cloning and Sanger sequencing that lack the throughput of current sequencing platforms. In this study, we present a high throughput, full-length allele discovery method utilizing Pacific Biosciences circular consensus sequencing (CCS). We also describe a new approach to Macaque Exome Sequencing (MES) and the development of the Rhexome1.0, an adapted target capture reagent that includes macaque-specific capture probe sets. By using sequence reads generated by whole genome sequencing (WGS) and MES to inform primer design, we were able to increase the sensitivity of KIR allele discovery. We demonstrate this increased sensitivity by defining nine novel alleles within a cohort of Mauritian cynomolgus macaques (MCM), a geographically isolated population with restricted KIR genetics that was thought to be completely characterized. Finally, we describe an approach to genotyping KIRs directly from sequence reads generated using WGS/MES reads. The findings presented here expand our understanding of KIR genetics in MCM by associating new genes with all eight KIR haplotypes and demonstrating the existence of at least one KIR3DS gene associated with every haplotype.

Standardization of a Videofluoroscopic Swallow Study Protocol to Investigate Dysphagia in Dogs

Background

Videofluoroscopic swallow study (VFSS) is the gold standard for diagnosis of dysphagia in veterinary medicine but lacks standardized protocols that emulate physiologic feeding practices. Age impacts swallow function in humans but has not been evaluated by VFSS in dogs.

Hypothesis/Objectives

To develop a protocol with custom kennels designed to allow free‐feeding of 3 optimized formulations of contrast media and diets that address limitations of current VFSS protocols. We hypothesized that dogs evaluated by a free‐feeding VFSS protocol would show differences in objective swallow metrics based on age.

Animals

Healthy juvenile, adult, and geriatric dogs (n = 24).

Methods

Prospective, experimental study. Custom kennels were developed to maintain natural feeding behaviors during VFSS. Three food consistencies (thin liquid, pureed food, and dry kibble) were formulated with either iohexol or barium to maximize palatability and voluntary prehension. Dogs were evaluated by 16 swallow metrics and compared across age groups.

Results

Development of a standardized VFSS protocol resulted in successful collection of swallow data in healthy dogs. No significant differences in swallow metrics were observed among age groups. Substantial variability was observed in healthy dogs when evaluated under these physiologic conditions. Features typically attributed to pathologic states, such as gastric reflux, were seen in healthy dogs.

Conclusions and Clinical Importance

Development of a VFSS protocol that reflects natural feeding practices may allow emulation of physiology resulting in clinical signs of dysphagia. Age did not result in significant changes in swallow metrics, but additional studies are needed, particularly in light of substantial normal variation.

Primary Human Placental Trophoblasts are Permissive for Zika Virus (ZIKV) Replication

Zika virus (ZIKV) is an emerging mosquito-borne (Aedes genus) arbovirus of the Flaviviridae family. Although ZIKV has been predominately associated with a mild or asymptomatic dengue-like disease, its appearance in the Americas has been accompanied by a multi-fold increase in reported incidence of fetal microcephaly and brain malformations. The source and mode of vertical transmission from mother to fetus is presumptively transplacental, although a causal link explaining the interval delay between maternal symptoms and observed fetal malformations following infection has been missing. In this study, we show that primary human placental trophoblasts from non-exposed donors (n = 20) can be infected by primary passage ZIKV-FLR isolate, and uniquely allowed for ZIKV viral RNA replication when compared to dengue virus (DENV). Consistent with their being permissive for ZIKV infection, primary trophoblasts expressed multiple putative ZIKV cell entry receptors, and cellular function and differentiation were preserved. These findings suggest that ZIKV-FLR strain can replicate in human placental trophoblasts without host cell destruction, thereby serving as a likely permissive reservoir and portal of fetal transmission with risk of latent microcephaly and malformations.