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.

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.

Effects of early daily alcohol exposure on placental function and fetal growth in a rhesus macaque model.

BACKGROUND: Prenatal alcohol exposure is the most common cause of birth defects and intellectual disabilities and can increase the risk of stillbirth and negatively impact fetal growth. OBJECTIVE: To determine the effect of early prenatal alcohol exposure on nonhuman primate placental function and fetal growth. We hypothesized that early chronic prenatal alcohol would alter placental perfusion and oxygen availability that adversely affects fetal growth. STUDY DESIGN: Rhesus macaques self-administered 1.5 g/kg/d of ethanol (n=12) or isocaloric maltose-dextrin (n=12) daily before conception through the first 60 days of gestation (term is approximately 168 days). All animals were serially imaged with Doppler ultrasound to measure fetal biometry, uterine artery volume blood flow, and placental volume blood flow. Following Doppler ultrasound, all animals underwent both blood oxygenation level-dependent magnetic resonance imaging to characterize placental blood oxygenation and dynamic contrast-enhanced magnetic resonance imaging to quantify maternal placental perfusion. Animals were delivered by cesarean delivery for placental collection and fetal necropsy at gestational days 85 (n=8), 110 (n=8), or 135 (n=8). Histologic and RNA-sequencing analyses were performed on collected placental tissue. RESULTS: Placental volume blood flow was decreased at all gestational time points in ethanol-exposed vs control animals, but most significantly at gestational day 110 by Doppler ultrasound (P<.05). A significant decrease in total volumetric blood flow occurred in ethanol-exposed vs control animals on dynamic contrast-enhanced magnetic resonance imaging at both gestation days 110 and 135 (P<.05); moreover, a global reduction in T2∗, high blood deoxyhemoglobin concentration, occurred throughout gestation (P<.05). Similarly, evidence of placental ischemic injury was notable by histologic analysis, which revealed a significant increase in microscopic infarctions in ethanol-exposed, not control, animals, largely present at middle to late gestation. Fetal biometry and weight were decreased in ethanol-exposed vs control animals, but the decrease was not significant. Analysis with RNA sequencing suggested the involvement of the inflammatory and extracellular matrix response pathways. CONCLUSION: Early chronic prenatal alcohol exposure significantly diminished placental perfusion at mid to late gestation and also significantly decreased the oxygen supply to the fetal vasculature throughout pregnancy, these findings were associated with the presence of microscopic placental infarctions in the nonhuman primate. Although placental adaptations may compensate for early environmental perturbations to fetal growth, placental blood flow and oxygenation were reduced, consistent with the evidence of placental ischemic injury.

Disruption of the Francisella noatunensis subsp. orientalis pdpA Gene Results in Virulence Attenuation and Protection in Zebrafish.

Several Francisella spp., including Francisella noatunensis, are regarded as important emerging pathogens of wild and farmed fish. However, very few studies have investigated the virulence factors that allow these bacterial species to be pathogenic in fish. The Francisella pathogenicity island (FPI) is a well-described, gene-dense region encoding major virulence factors for the genus Francisella. pdpA is a member of the pathogenicity-determining protein genes carried by the FPI that are implicated in the ability of the mammalian pathogen Francisella tularensis to escape and replicate in infected host cells. Using a sacB suicide approach, we generated pdpA knockouts to address the role of PdpA as a virulence factor for F. noatunensis. Because polarity can be an issue in gene-dense regions, we generated two different marker-based mutants in opposing polarity (the F. noatunensis subsp. orientalis ΔpdpA1 and ΔpdpA2 strains). Both mutants were attenuated (P < 0.0001) in zebrafish challenges and displayed impaired intracellular replication (P < 0.05) and cytotoxicity (P < 0.05), all of which could be restored to wild-type (WT) levels by complementation for the ΔpdpA1 mutant. Importantly, differences were found for bacterial burden and induction of acute-phase and proinflammatory genes for the F. noatunensis subsp. orientalis ΔpdpA1 and ΔpdpA2 mutants compared to the WT during acute infection. In addition, neither mutant resulted in significant histopathological changes. Finally, immunization with the F. noatunensis subsp. orientalis ΔpdpA1 mutant led to protection (P < 0.012) against an acute 40% lethal dose (LD40) challenge with WT F. noatunensis in the zebrafish model of infection. Taken together, the results from this study further demonstrate physiological similarities within the genus Francisella relative to their phylogenetic relationships and the utility of zebrafish for addressing virulence factors for the genus.

Effect of chronic ethanol consumption in rhesus macaques on the nucleus accumbens core transcriptome.

The nucleus accumbens core (NAcc) has been repeatedly demonstrated to be a key component of the circuitry associated with excessive ethanol consumption. Previous studies have illustrated that in a nonhuman primate (NHP) model of chronic ethanol consumption, there is significant epigenetic remodeling of the NAcc. In the current study, RNA-Seq was used to examine genome-wide gene expression in eight each of control, low/binge (LD*), and high/very high (HD*) rhesus macaque drinkers. Using an FDR < 0.05, zero genes were significantly differentially expressed (DE) between LD* and controls, six genes between HD* and LD*, and 734 genes between HD* and controls. Focusing on HD* versus control DE genes, the upregulated genes (N = 366) were enriched in genes with annotations associated with signal recognition particle (SRP)-dependent co-translational protein targeting to membrane (FDR < 3 × 10-59 ), structural constituent of ribosome (FDR < 3 × 10-47 ), and ribosomal subunit (FDR < 5 × 10-48 ). Downregulated genes (N = 363) were enriched in annotations associated with behavior (FDR < 2 × 10-4 ), membrane organization (FDR < 1 × 10-4 ), inorganic cation transmembrane transporter activity (FDR < 2 × 10-3 ), synapse part (FDR < 4 × 10-10 ), glutamatergic synapse (FDR < 1 × 10-6 ), and GABAergic synapse (FDR < 6 × 10-4 ). Ingenuity Pathway Analysis (IPA) revealed that EIF2 signaling and mTOR pathways were significantly upregulated in HD* animals (FDR < 3 × 10-33 and <2 × 10-16 , respectively). Overall, the data supported our working hypothesis; excessive consumption would be associated with transcriptional differences in GABA/glutamate-related genes.

Effects of repeated alcohol abstinence on within-subject prefrontal cortical gene expression in rhesus macaques.

Male rhesus monkeys (n = 24) had a biopsy of prefrontal cortical area 46 prior to chronic ethanol self-administration (n = 17) or caloric control (n = 7). Fourteen months of daily self-administration (water vs. 4% alcohol, 22 h access/day termed "open-access") was followed by two cycles of prolonged abstinence (5 weeks) each followed by 3 months of open-access alcohol and a final abstinence followed by necropsy. At necropsy, a biopsy of Area 46, contralateral to the original biopsy, was obtained. Gene expression data (RNA-Seq) were collected comparing biopsy/necropsy samples. Monkeys were categorized by drinking status during the final post-abstinent drinking phase as light (LD), binge (BD), heavy (HD) and very heavy (VHD drinkers). Comparing pre-ethanol to post-abstinent biopsies, four animals that converted from HD to VHD status had significant ontology enrichments in downregulated genes (necropsy minus biopsy n = 286) that included immune response (FDR < 9 × 10-7) and plasma membrane changes (FDR < 1 × 10-7). Genes in the immune response category included IL16 and 18, CCR1, B2M, TLR3, 6 and 7, SP2 and CX3CR1. Upregulated genes (N = 388) were particularly enriched in genes associated with the negative regulation of MAP kinase activity (FDR < 3 × 10-5), including DUSP 1, 4, 5, 6 and 18, SPRY 2, 3, and 4, SPRED2, BMP4 and RGS2. Overall, these data illustrate the power of the NHP model and the within-subject design of genomic changes due to alcohol and suggest new targets for treating severe escalated drinking following repeated alcohol abstinence attempts.

Transcriptional profiling of mucus production and modification in rhesus macaque endocervical cells under hormonal regulation.

Objective: Endocervical mucus production is a key regulator of fertility throughout the menstrual cycle. With cycle-dependent variability in mucus quality and quantity, cervical mucus can either facilitate or block sperm ascension into the upper female reproductive tract. This study seeks to identify genes involved in the hormonal regulation of mucus production, modification, and regulation through profiling the transcriptome of endocervical cells from the non-human primate, the Rhesus Macaque (Macaca mulatta). Design: Experimental. Setting: Translational science laboratory. Intervention: We treated differentiated primary endocervical cultures with estradiol (E2) and progesterone (P4) to mimic peri-ovulatory and luteal-phase hormonal changes. Using RNA-sequencing, we identified differential expression of gene pathways and mucus producing and modifying genes in cells treated with E2 compared to hormone-free conditions and E2 compared to E2-primed cells treated with P4. Main Outcome Measures: We pursued differential gene expression analysis on RNA-sequenced cells. Sequence validation was done using qPCR. Results: Our study identified 158 genes that show significant differential expression in E2-only conditions compared to hormone-free control, and 250 genes that show significant differential expression in P4-treated conditions compared to E2-only conditions. From this list, we found hormone-induced changes in transcriptional profiles for genes across several classes of mucus production, including ion channels and enzymes involved in post-translational mucin modification that have not previously been described as hormonally regulated. Conclusion: Our study is the first to use an in vitro culture system to create an epithelial-cell specific transcriptome of the endocervix. As a result, our study identifies new genes and pathways that are altered by sex-steroids in cervical mucus production.

Mayaro virus pathogenesis and immunity in rhesus macaques.

Mayaro virus (MAYV) is a mosquito-transmitted alphavirus that causes debilitating and persistent arthritogenic disease. While MAYV was previously reported to infect non-human primates (NHP), characterization of MAYV pathogenesis is currently lacking. Therefore, in this study we characterized MAYV infection and immunity in rhesus macaques. To inform the selection of a viral strain for NHP experiments, we evaluated five MAYV strains in C57BL/6 mice and showed that MAYV strain BeAr505411 induced robust tissue dissemination and disease. Three male rhesus macaques were subcutaneously challenged with 105 plaque-forming units of this strain into the arms. Peak plasma viremia occurred at 2 days post-infection (dpi). NHPs were taken to necropsy at 10 dpi to assess viral dissemination, which included the muscles and joints, lymphoid tissues, major organs, male reproductive tissues, as well as peripheral and central nervous system tissues. Histological examination demonstrated that MAYV infection was associated with appendicular joint and muscle inflammation as well as presence of perivascular inflammation in a wide variety of tissues. One animal developed a maculopapular rash and two NHP had viral RNA detected in upper torso skin samples, which was associated with the presence of perivascular and perifollicular lymphocytic aggregation. Analysis of longitudinal peripheral blood samples indicated a robust innate and adaptive immune activation, including the presence of anti-MAYV neutralizing antibodies with activity against related Una virus and chikungunya virus. Inflammatory cytokines and monocyte activation also peaked coincident with viremia, which was well supported by our transcriptomic analysis highlighting enrichment of interferon signaling and other antiviral processes at 2 days post MAYV infection. The rhesus macaque model of MAYV infection recapitulates many of the aspects of human infection and is poised to facilitate the evaluation of novel therapies and vaccines targeting this re-emerging virus.

Transcriptomic profiles conducive to immune-mediated tumor rejection in human breast cancer skin metastases treated with Imiquimod.

Imiquimod is a topical toll-like-receptor-7 agonist currently used for treating basal cell carcinoma. Recently, imiquimod has demonstrated tumor regression in melanoma and breast cancer skin metastases. However, the molecular perturbations induced by imiquimod in breast cancer metastases have not been previously characterized. Here, we describe transcriptomic profiles associated with responsiveness to imiquimod in breast cancer skin metastases. Baseline and post-treatment tumor samples from patients treated with imiquimod in a clinical trial were profiled using Nanostring technology. Through an integrative analytic pipeline, we showed that tumors from patients who achieved a durable clinical response displayed a permissive microenvironment, substantiated by the upregulation of transcripts encoding for molecules involved in leukocyte adhesion and migration, cytotoxic functions, and antigen presentation. In responding patients, Imiquimod triggered a strong T-helper-1 (Th-1)/cytotoxic immune response, characterized by the coordinated upregulation of Th-1 chemokines, migration of Th-1 and cytotoxic T cells into the tumor, and activation of immune-effector functions, ultimately mediating tumor destruction. In conclusion, we have shown that topical imiquimod can induce a robust immune response in breast cancer metastases, and this response is more likely to occur in tumors with a pre-activated microenvironment. In this setting, imiquimod could be utilized in combination with other targeted immunotherapies to increase therapeutic efficacy.

Prognostic role of elevated mir-24-3p in breast cancer and its association with the metastatic process.

MicroRNAs have been shown to play important roles in breast cancer progression and can serve as biomarkers. To assess the prognostic role of a panel of miRNAs in breast cancer, we collected plasma prospectively at the time of initial diagnosis from 1,780 patients with stage I-III breast cancer prior to definitive treatment. We identified plasma from 115 patients who subsequently developed distant metastases and 115 patients without metastatic disease. Both groups were matched by: age at blood collection, year of blood collection, breast cancer subtype, and stage. The median follow up was 3.4 years (range, 1-9 years). We extracted RNA from plasma and analyzed the expression of 800 miRNAs using Nanostring technology. We then assessed the expression of miRNAs in primary and metastatic breast cancer samples from The Cancer Genome Atlas (TCGA). We found that, miR-24-3p was upregulated in patients with metastases, both in plasma and in breast cancer tissues. Patients whose primary tumors expressed high levels of miR-24-3p had a significantly lower survival rate compared to patients with low mir-24-3p levels in the TCGA cohort (n=1,024). RNA-Seq data of the samples with the highest miR-24-3p expression versus those with the lowest miR-24-3p in the TCGA cohort identified a specific gene expression signature for those tumors with high miR-24-3p. Possible target genes for miR-24-3p were predicted based on gene expression and binding site, and their effects on cancer pathways were evaluated. Cancer, breast cancer and proteoglycans were the top three pathways affected by miR-24-3p overexpression.

Recurrent homozygous deletion of DROSHA and microduplication of PDE4DIP in pineoblastoma.

Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. Here, we analyzed pediatric and adult pineoblastoma samples (n = 23) using a combination of genome-wide DNA methylation profiling and whole-exome sequencing or whole-genome sequencing. Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower-grade pineal tumors and normal pineal gland. Recurrent variants were found in genes involved in PKA- and NF-κB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expresion of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain.

Apolipoprotein L1 risk variants associate with prevalent atherosclerotic disease in African American systemic lupus erythematosus patients.

OBJECTIVE: Atherosclerosis is exaggerated in African American (AA) systemic lupus erythematosus (SLE) patients, with doubled cardiovascular disease (CVD) risk compared to White patients. The extent to which common Apolipoprotein L1 (APOL1) risk alleles (RA) contribute to this trend is unknown. This retrospective cohort study assessed prevalent atherosclerotic disease across APOL1 genotypes in AA SLE patients. METHODS: One hundred thirteen AA SLE subjects were APOL1-genotyped and stratified as having: zero risk alleles, one risk allele, or two risk alleles. Chart review assessed CVD manifestations including abdominal aortic aneurysm, angina, carotid artery disease, coronary artery disease, myocardial infarction, peripheral vascular disease, stroke, and vascular calcifications. Associations between the genotypes and a composite endpoint defined as one or more CVD manifestations were calculated using logistic regression. Symptomatic atherosclerotic disease, excluding incidental vascular calcifications, was also assessed. RESULTS: The 0-risk-allele, 1-risk-allele and 2-risk-allele groups, respectively, comprised 34%, 53%, and 13% of the cohort. Respectively, 13.2%, 41.7%, and 60.0% of the 0-risk allele, 1-risk-allele, and 2-risk-allele groups met the composite endpoint of atherosclerotic CVD (p = 0.001). Adjusting for risk factors-including smoking, ESRD, BMI >25 and hypertension-we observed an association between carrying one or more RA and atherosclerotic CVD (OR = 7.1; p = 0.002). For symptomatic disease, the OR was 3.5 (p = 0.02). In a time-to-event analysis, the proportion of subjects free from the composite primary endpoint, symptomatic atherosclerotic CVD, was higher in the 0-risk-allele group compared to the 1-risk-allele and 2-risk-allele groups (χ2 = 6.5; p = 0.04). CONCLUSIONS: Taken together, the APOL1 RAs associate with prevalent atherosclerotic CVD in this cohort of AA SLE patients, perhaps reflecting a potentiating effect of SLE on APOL1-related cardiovascular phenotypes.