RNA editing underlies genetic risk of common inflammatory diseases.

A major challenge in human genetics is to identify the molecular mechanisms of trait-associated and disease-associated variants. To achieve this, quantitative trait locus (QTL) mapping of genetic variants with intermediate molecular phenotypes such as gene expression and splicing have been widely adopted1,2. However, despite successes, the molecular basis for a considerable fraction of trait-associated and disease-associated variants remains unclear3,4. Here we show that ADAR-mediated adenosine-to-inosine RNA editing, a post-transcriptional event vital for suppressing cellular double-stranded RNA (dsRNA)-mediated innate immune interferon responses5-11, is an important potential mechanism underlying genetic variants associated with common inflammatory diseases. We identified and characterized 30,319 cis-RNA editing QTLs (edQTLs) across 49 human tissues. These edQTLs were significantly enriched in genome-wide association study signals for autoimmune and immune-mediated diseases. Colocalization analysis of edQTLs with disease risk loci further pinpointed key, putatively immunogenic dsRNAs formed by expected inverted repeat Alu elements as well as unexpected, highly over-represented cis-natural antisense transcripts. Furthermore, inflammatory disease risk variants, in aggregate, were associated with reduced editing of nearby dsRNAs and induced interferon responses in inflammatory diseases. This unique directional effect agrees with the established mechanism that lack of RNA editing by ADAR1 leads to the specific activation of the dsRNA sensor MDA5 and subsequent interferon responses and inflammation7-9. Our findings implicate cellular dsRNA editing and sensing as a previously underappreciated mechanism of common inflammatory diseases.

Learning cis-regulatory principles of ADAR-based RNA editing from CRISPR-mediated mutagenesis.

Adenosine-to-inosine (A-to-I) RNA editing catalyzed by ADAR enzymes occurs in double-stranded RNAs. Despite a compelling need towards predictive understanding of natural and engineered editing events, how the RNA sequence and structure determine the editing efficiency and specificity (i.e., cis-regulation) is poorly understood. We apply a CRISPR/Cas9-mediated saturation mutagenesis approach to generate libraries of mutations near three natural editing substrates at their endogenous genomic loci. We use machine learning to integrate diverse RNA sequence and structure features to model editing levels measured by deep sequencing. We confirm known features and identify new features important for RNA editing. Training and testing XGBoost algorithm within the same substrate yield models that explain 68 to 86 percent of substrate-specific variation in editing levels. However, the models do not generalize across substrates, suggesting complex and context-dependent regulation patterns. Our integrative approach can be applied to larger scale experiments towards deciphering the RNA editing code.

p53 is a central regulator driving neurodegeneration caused by C9orf72 poly(PR).

The most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is a GGGGCC repeat expansion in the C9orf72 gene. We developed a platform to interrogate the chromatin accessibility landscape and transcriptional program within neurons during degeneration. We provide evidence that neurons expressing the dipeptide repeat protein poly(proline-arginine), translated from the C9orf72 repeat expansion, activate a highly specific transcriptional program, exemplified by a single transcription factor, p53. Ablating p53 in mice completely rescued neurons from degeneration and markedly increased survival in a C9orf72 mouse model. p53 reduction also rescued axonal degeneration caused by poly(glycine-arginine), increased survival of C9orf72 ALS/FTD-patient-induced pluripotent stem cell (iPSC)-derived motor neurons, and mitigated neurodegeneration in a C9orf72 fly model. We show that p53 activates a downstream transcriptional program, including Puma, which drives neurodegeneration. These data demonstrate a neurodegenerative mechanism dynamically regulated through transcription-factor-binding events and provide a framework to apply chromatin accessibility and transcription program profiles to neurodegeneration.

Widespread RNA editing dysregulation in brains from autistic individuals.

Transcriptomic analyses of postmortem brains have begun to elucidate molecular abnormalities in autism spectrum disorder (ASD). However, a crucial pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here we profiled global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of postmortem brains of people with ASD. We observed a global bias for hypoediting in ASD brains, which was shared across brain regions and involved many synaptic genes. We show that the Fragile X proteins FMRP and FXR1P interact with RNA-editing enzymes (ADAR proteins) and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA-editing alterations in ASD and Fragile X syndrome, establishing this as a molecular link between these related diseases. Our findings, which are corroborated across multiple data sets, including dup15q (genomic duplication of 15q11.2-13.1) cases associated with intellectual disability, highlight RNA-editing dysregulation in ASD and reveal new mechanisms underlying this disorder.

Pre-reproductive stress and fluoxetine treatment in rats affect offspring A-to-I RNA editing, gene expression and social behavior.

Adenosine to inosine RNA editing is an epigenetic process that entails site-specific modifications in double-stranded RNA molecules, catalyzed by adenosine deaminases acting on RNA (ADARs). Using the multiplex microfluidic PCR and deep sequencing technique, we recently showed that exposing adolescent female rats to chronic unpredictable stress before reproduction affects editing in the prefrontal cortex and amygdala of their newborn offspring, particularly at the serotonin receptor 5-HT2c (encoded by Htr2c). Here, we used the same technique to determine whether post-stress, pre-reproductive maternal treatment with fluoxetine (5 mg/kg, 7 days) reverses the effects of stress on editing. We also examined the mRNA expression of ADAR enzymes in these regions, and asked whether social behavior in adult offspring would be altered by maternal exposure to stress and/or fluoxetine. Maternal treatment with fluoxetine altered Htr2c editing in offspring amygdala at birth, enhanced the expression of Htr2c mRNA and RNA editing enzymes in the prefrontal cortex, and reversed the effects of pre-reproductive stress on Htr2c editing in this region. Furthermore, maternal fluoxetine treatment enhanced differences in editing of glutamate receptors between offspring of control and stress-exposed rats, and led to enhanced social preference in adult offspring. Our findings indicate that pre-gestational fluoxetine treatment affects patterns of RNA editing and editing enzyme expression in neonatal offspring brain in a region-specific manner, in interaction with pre-reproductive stress. Overall, these findings imply that fluoxetine treatment affects serotonergic signaling in offspring brain even when treatment is discontinued before gestation, and its effects may depend upon prior exposure to stress.

A-to-I RNA editing in the rat brain is age-dependent, region-specific and sensitive to environmental stress across generations.

BACKGROUND: Adenosine-to-inosine (A-to-I) RNA editing is an epigenetic modification catalyzed by adenosine deaminases acting on RNA (ADARs), and is especially prevalent in the brain. We used the highly accurate microfluidics-based multiplex PCR sequencing (mmPCR-seq) technique to assess the effects of development and environmental stress on A-to-I editing at 146 pre-selected, conserved sites in the rat prefrontal cortex and amygdala. Furthermore, we asked whether changes in editing can be observed in offspring of stress-exposed rats. In parallel, we assessed changes in ADARs expression levels. RESULTS: In agreement with previous studies, we found editing to be generally higher in adult compared to neonatal rat brain. At birth, editing was generally lower in prefrontal cortex than in amygdala. Stress affected editing at the serotonin receptor 2c (Htr2c), and editing at this site was significantly altered in offspring of rats exposed to prereproductive stress across two generations. Stress-induced changes in Htr2c editing measured with mmPCR-seq were comparable to changes measured with Sanger and Illumina sequencing. Developmental and stress-induced changes in Adar and Adarb1 mRNA expression were observed but did not correlate with editing changes. CONCLUSIONS: Our findings indicate that mmPCR-seq can accurately detect A-to-I RNA editing in rat brain samples, and confirm previous accounts of a developmental increase in RNA editing rates. Our findings also point to stress in adolescence as an environmental factor that alters RNA editing patterns several generations forward, joining a growing body of literature describing the transgenerational effects of stress.

Dynamic landscape and regulation of RNA editing in mammals.

Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules. Although many editing sites have recently been discovered, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

Molecular definition of a metastatic lung cancer state reveals a targetable CD109-Janus kinase-Stat axis.

Lung cancer is the leading cause of cancer deaths worldwide, with the majority of mortality resulting from metastatic spread. However, the molecular mechanism by which cancer cells acquire the ability to disseminate from primary tumors, seed distant organs, and grow into tissue-destructive metastases remains incompletely understood. We combined tumor barcoding in a mouse model of human lung adenocarcinoma with unbiased genomic approaches to identify a transcriptional program that confers metastatic ability and predicts patient survival. Small-scale in vivo screening identified several genes, including Cd109, that encode novel pro-metastatic factors. We uncovered signaling mediated by Janus kinases (Jaks) and the transcription factor Stat3 as a critical, pharmacologically targetable effector of CD109-driven lung cancer metastasis. In summary, by coupling the systematic genomic analysis of purified cancer cells in distinct malignant states from mouse models with extensive human validation, we uncovered several key regulators of metastatic ability, including an actionable pro-metastatic CD109-Jak-Stat3 axis.

Adenosine-to-inosine RNA editing by ADAR1 is essential for normal murine erythropoiesis.

Adenosine deaminases that act on RNA (ADARs) convert adenosine residues to inosine in double-stranded RNA. In vivo, ADAR1 is essential for the maintenance of hematopoietic stem/progenitors. Whether other hematopoietic cell types also require ADAR1 has not been assessed. Using erythroid- and myeloid-restricted deletion of Adar1, we demonstrate that ADAR1 is dispensable for myelopoiesis but is essential for normal erythropoiesis. Adar1-deficient erythroid cells display a profound activation of innate immune signaling and high levels of cell death. No changes in microRNA levels were found in ADAR1-deficient erythroid cells. Using an editing-deficient allele, we demonstrate that RNA editing is the essential function of ADAR1 during erythropoiesis. Mapping of adenosine-to-inosine editing in purified erythroid cells identified clusters of hyperedited adenosines located in long 3'-untranslated regions of erythroid-specific transcripts and these are ADAR1-specific editing events. ADAR1-mediated RNA editing is essential for normal erythropoiesis.

Editing of Cellular Self-RNAs by Adenosine Deaminase ADAR1 Suppresses Innate Immune Stress Responses.

Adenosine deaminases acting on double-stranded RNA (ADARs) catalyze the deamination of adenosine (A) to produce inosine (I) in double-stranded (ds) RNA structures, a process known as A-to-I RNA editing. dsRNA is an important trigger of innate immune responses, including interferon (IFN) production and action. We examined the role of A-to-I RNA editing by two ADARs, ADAR1 and ADAR2, in the sensing of self-RNA in the absence of pathogen infection, leading to activation of IFN-induced, RNA-mediated responses in mouse embryo fibroblasts. IFN treatment of Adar1(-/-) cells lacking both the p110 constitutive and p150 IFN-inducible ADAR1 proteins induced formation of stress granules, whereas neither wild-type (WT) nor Adar2(-/-) cells displayed a comparable stress granule response following IFN treatment. Phosphorylation of protein synthesis initiation factor eIF2α at serine 51 was increased in IFN-treated Adar1(-/-) cells but not in either WT or Adar2(-/-) cells following IFN treatment. Analysis by deep sequencing of mouse exonic loci containing A-to-I-editing sites revealed that the majority of editing in mouse embryo fibroblasts was carried out by ADAR1. IFN treatment increased editing in both WT and Adar2(-/-) cells but not in either Adar1(-/-) or Adar1(-/-) (p150) cells or Stat1(-/-) or Stat2(-/-) cells. Hyper-edited sites found in predicted duplex structures showed strand bias of editing for some RNAs. These results implicate ADAR1 p150 as the major A-to-I editor in mouse embryo fibroblasts, acting as a feedback suppressor of innate immune responses otherwise triggered by self-RNAs possessing regions of double-stranded character.

Genetic mapping uncovers cis-regulatory landscape of RNA editing.

Adenosine-to-inosine (A-to-I) RNA editing, catalysed by ADAR enzymes conserved in metazoans, plays an important role in neurological functions. Although the fine-tuning mechanism provided by A-to-I RNA editing is important, the underlying rules governing ADAR substrate recognition are not well understood. We apply a quantitative trait loci (QTL) mapping approach to identify genetic variants associated with variability in RNA editing. With very accurate measurement of RNA editing levels at 789 sites in 131 Drosophila melanogaster strains, here we identify 545 editing QTLs (edQTLs) associated with differences in RNA editing. We demonstrate that many edQTLs can act through changes in the local secondary structure for edited dsRNAs. Furthermore, we find that edQTLs located outside of the edited dsRNA duplex are enriched in secondary structure, suggesting that distal dsRNA structure beyond the editing site duplex affects RNA editing efficiency. Our work will facilitate the understanding of the cis-regulatory code of RNA editing.

RNA editing by ADAR1 prevents MDA5 sensing of endogenous dsRNA as nonself.

Adenosine-to-inosine (A-to-I) editing is a highly prevalent posttranscriptional modification of RNA, mediated by ADAR (adenosine deaminase acting on RNA) enzymes. In addition to RNA editing, additional functions have been proposed for ADAR1. To determine the specific role of RNA editing by ADAR1, we generated mice with an editing-deficient knock-in mutation (Adar1(E861A), where E861A denotes Glu(861)→Ala(861)). Adar1(E861A/E861A) embryos died at ~E13.5 (embryonic day 13.5), with activated interferon and double-stranded RNA (dsRNA)-sensing pathways. Genome-wide analysis of the in vivo substrates of ADAR1 identified clustered hyperediting within long dsRNA stem loops within 3' untranslated regions of endogenous transcripts. Finally, embryonic death and phenotypes of Adar1(E861A/E861A) were rescued by concurrent deletion of the cytosolic sensor of dsRNA, MDA5. A-to-I editing of endogenous dsRNA is the essential function of ADAR1, preventing the activation of the cytosolic dsRNA response by endogenous transcripts.

RADAR: a rigorously annotated database of A-to-I RNA editing.

We present RADAR--a rigorously annotated database of A-to-I RNA editing (available at http://RNAedit.com). The identification of A-to-I RNA editing sites has been dramatically accelerated in the past few years by high-throughput RNA sequencing studies. RADAR includes a comprehensive collection of A-to-I RNA editing sites identified in humans (Homo sapiens), mice (Mus musculus) and flies (Drosophila melanogaster), together with extensive manually curated annotations for each editing site. RADAR also includes an expandable listing of tissue-specific editing levels for each editing site, which will facilitate the assignment of biological functions to specific editing sites.

Whole-exome resequencing distinguishes cystic kidney diseases from phenocopies in renal ciliopathies.

Rare single-gene disorders cause chronic disease. However, half of the 6000 recessive single gene causes of disease are still unknown. Because recessive disease genes can illuminate, at least in part, disease pathomechanism, their identification offers direct opportunities for improved clinical management and potentially treatment. Rare diseases comprise the majority of chronic kidney disease (CKD) in children but are notoriously difficult to diagnose. Whole-exome resequencing facilitates identification of recessive disease genes. However, its utility is impeded by the large number of genetic variants detected. We here overcome this limitation by combining homozygosity mapping with whole-exome resequencing in 10 sib pairs with a nephronophthisis-related ciliopathy, which represents the most frequent genetic cause of CKD in the first three decades of life. In 7 of 10 sibships with a histologic or ultrasonographic diagnosis of nephronophthisis-related ciliopathy, we detect the causative gene. In six sibships, we identify mutations of known nephronophthisis-related ciliopathy genes, while in two additional sibships we found mutations in the known CKD-causing genes SLC4A1 and AGXT as phenocopies of nephronophthisis-related ciliopathy. Thus, whole-exome resequencing establishes an efficient, noninvasive approach towards early detection and causation-based diagnosis of rare kidney diseases. This approach can be extended to other rare recessive disorders, thereby providing accurate diagnosis and facilitating the study of disease mechanisms.

Identifying RNA editing sites using RNA sequencing data alone.

We show that RNA editing sites can be called with high confidence using RNA sequencing data from multiple samples across either individuals or species, without the need for matched genomic DNA sequence. We identified many previously unidentified editing sites in both humans and Drosophila; our results nearly double the known number of human protein recoding events. We also found that human genes harboring conserved editing sites within Alu repeats are enriched for neuronal functions.

Exome capture reveals ZNF423 and CEP164 mutations, linking renal ciliopathies to DNA damage response signaling.

Nephronophthisis-related ciliopathies (NPHP-RC) are degenerative recessive diseases that affect kidney, retina, and brain. Genetic defects in NPHP gene products that localize to cilia and centrosomes defined them as "ciliopathies." However, disease mechanisms remain poorly understood. Here, we identify by whole-exome resequencing, mutations of MRE11, ZNF423, and CEP164 as causing NPHP-RC. All three genes function within the DNA damage response (DDR) pathway. We demonstrate that, upon induced DNA damage, the NPHP-RC proteins ZNF423, CEP164, and NPHP10 colocalize to nuclear foci positive for TIP60, known to activate ATM at sites of DNA damage. We show that knockdown of CEP164 or ZNF423 causes sensitivity to DNA damaging agents and that cep164 knockdown in zebrafish results in dysregulated DDR and an NPHP-RC phenotype. Our findings link degenerative diseases of the kidney and retina, disorders of increasing prevalence, to mechanisms of DDR.

FAN1 mutations cause karyomegalic interstitial nephritis, linking chronic kidney failure to defective DNA damage repair.

Chronic kidney disease (CKD) represents a major health burden. Its central feature of renal fibrosis is not well understood. By exome sequencing, we identified mutations in FAN1 as a cause of karyomegalic interstitial nephritis (KIN), a disorder that serves as a model for renal fibrosis. Renal histology in KIN is indistinguishable from that of nephronophthisis, except for the presence of karyomegaly. The FAN1 protein has nuclease activity and acts in DNA interstrand cross-link (ICL) repair within the Fanconi anemia DNA damage response (DDR) pathway. We show that cells from individuals with FAN1 mutations have sensitivity to the ICL-inducing agent mitomycin C but do not exhibit chromosome breakage or cell cycle arrest after diepoxybutane treatment, unlike cells from individuals with Fanconi anemia. We complemented ICL sensitivity with wild-type FAN1 but not with cDNA having mutations found in individuals with KIN. Depletion of fan1 in zebrafish caused increased DDR, apoptosis and kidney cysts. Our findings implicate susceptibility to environmental genotoxins and inadequate DNA repair as novel mechanisms contributing to renal fibrosis and CKD.

Accurate identification of human Alu and non-Alu RNA editing sites.

We developed a computational framework to robustly identify RNA editing sites using transcriptome and genome deep-sequencing data from the same individual. As compared with previous methods, our approach identified a large number of Alu and non-Alu RNA editing sites with high specificity. We also found that editing of non-Alu sites appears to be dependent on nearby edited Alu sites, possibly through the locally formed double-stranded RNA structure.

Genotype-phenotype correlation in 440 patients with NPHP-related ciliopathies.

Nephronophthisis (NPHP), an autosomal recessive cystic kidney disease, is the most frequent genetic cause for end-stage renal failure in the first three decades of life. Mutations in 13 genes (NPHP1-NPHP11, AHI1, and CC2D2A) cause NPHP with ubiquitous expression of the corresponding proteins consistent with the multiorgan involvement of NPHP-related diseases. The genotype-phenotype correlation in these ciliopathies can be explained by gene locus heterogeneity, allelism, and the impact of modifier genes. In some NPHP-related ciliopathies, the nature of the recessive mutations determines disease severity. In order to define the genotype-phenotype correlation more clearly, we evaluated a worldwide cohort of 440 patients from 365 families with NPHP-related ciliopathies, in whom both disease-causing alleles were identified. The phenotypes were ranked in the order of severity from degenerative to degenerative/dysplastic to dysplastic. A genotype of two null alleles caused a range of phenotypes, with an increasing order of severity of NPHP1, NPHP3, NPHP4, NPHP5, NPHP2, NPHP10, NPHP6, to AHI1. Only NPHP6 showed allelic influences on the phenotypes; the presence of two null mutations caused dysplastic phenotypes, whereas at least one missense allele rescued it to a milder degenerative phenotype. We also found nine novel mutations in the NPHP genes. Thus, our studies have important implications for genetic counseling and planning of renal replacement therapy.

A transition zone complex regulates mammalian ciliogenesis and ciliary membrane composition.

Mutations affecting ciliary components cause ciliopathies. As described here, we investigated Tectonic1 (Tctn1), a regulator of mouse Hedgehog signaling, and found that it is essential for ciliogenesis in some, but not all, tissues. Cell types that do not require Tctn1 for ciliogenesis require it to localize select membrane-associated proteins to the cilium, including Arl13b, AC3, Smoothened and Pkd2. Tctn1 forms a complex with multiple ciliopathy proteins associated with Meckel and Joubert syndromes, including Mks1, Tmem216, Tmem67, Cep290, B9d1, Tctn2 and Cc2d2a. Components of this complex co-localize at the transition zone, a region between the basal body and ciliary axoneme. Like Tctn1, loss of Tctn2, Tmem67 or Cc2d2a causes tissue-specific defects in ciliogenesis and ciliary membrane composition. Consistent with a shared function for complex components, we identified a mutation in TCTN1 that causes Joubert syndrome. Thus, a transition zone complex of Meckel and Joubert syndrome proteins regulates ciliary assembly and trafficking, suggesting that transition zone dysfunction is the cause of these ciliopathies.