Transcriptional characterization of iPSC-derived microglia as a model for therapeutic development in neurodegeneration.

Microglia are the resident immune cells in the brain that play a key role in driving neuroinflammation, a hallmark of neurodegenerative disorders. Inducible microglia-like cells have been developed as an in vitro platform for molecular and therapeutic hypothesis generation and testing. However, there has been no systematic assessment of similarity of these cells to primary human microglia along with their responsiveness to external cues expected of primary cells in the brain. In this study, we performed transcriptional characterization of commercially available human inducible pluripotent stem cell (iPSC)-derived microglia-like (iMGL) cells by bulk and single cell RNA sequencing to assess their similarity with primary human microglia. To evaluate their stimulation responsiveness, iMGL cells were treated with Liver X Receptor (LXR) pathway agonists and their transcriptional responses characterized by bulk and single cell RNA sequencing. Bulk transcriptome analyses demonstrate that iMGL cells have a similar overall expression profile to freshly isolated human primary microglia and express many key microglial transcription factors and functional and disease-associated genes. Notably, at the single-cell level, iMGL cells exhibit distinct transcriptional subpopulations, representing both homeostatic and activated states present in normal and diseased primary microglia. Treatment of iMGL cells with LXR pathway agonists induces robust transcriptional changes in lipid metabolism and cell cycle at the bulk level. At the single cell level, we observe heterogeneity in responses between cell subpopulations in homeostatic and activated states and deconvolute bulk expression changes into their corresponding single cell states. In summary, our results demonstrate that iMGL cells exhibit a complex transcriptional profile and responsiveness, reminiscent of in vivo microglia, and thus represent a promising model system for therapeutic development in neurodegeneration.

Broad transcriptomic dysregulation occurs across the cerebral cortex in ASD.

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations1-3. In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex1,2,4-6. However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions.

Developmentally regulated alternate 3' end cleavage of nascent transcripts controls dynamic changes in protein expression in an adult stem cell lineage.

Alternative polyadenylation (APA) generates transcript isoforms that differ in the position of the 3' cleavage site, resulting in the production of mRNA isoforms with different length 3' UTRs. Although widespread, the role of APA in the biology of cells, tissues, and organisms has been controversial. We identified >500 Drosophila genes that express mRNA isoforms with a long 3' UTR in proliferating spermatogonia but a short 3' UTR in differentiating spermatocytes due to APA. We show that the stage-specific choice of the 3' end cleavage site can be regulated by the arrangement of a canonical polyadenylation signal (PAS) near the distal cleavage site but a variant or no recognizable PAS near the proximal cleavage site. The emergence of transcripts with shorter 3' UTRs in differentiating cells correlated with changes in expression of the encoded proteins, either from off in spermatogonia to on in spermatocytes or vice versa. Polysome gradient fractionation revealed >250 genes where the long 3' UTR versus short 3' UTR mRNA isoforms migrated differently, consistent with dramatic stage-specific changes in translation state. Thus, the developmentally regulated choice of an alternative site at which to make the 3' end cut that terminates nascent transcripts can profoundly affect the suite of proteins expressed as cells advance through sequential steps in a differentiation lineage.

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.

Coexpression network architecture reveals the brain-wide and multiregional basis of disease susceptibility.

Gene networks have yielded numerous neurobiological insights, yet an integrated view across brain regions is lacking. We leverage RNA sequencing in 864 samples representing 12 brain regions to robustly identify 12 brain-wide, 50 cross-regional and 114 region-specific coexpression modules. Nearly 40% of genes fall into brain-wide modules, while 25% comprise region-specific modules reflecting regional biology, such as oxytocin signaling in the hypothalamus, or addiction pathways in the nucleus accumbens. Schizophrenia and autism genetic risk are enriched in brain-wide and multiregional modules, indicative of broad impact; these modules implicate neuronal proliferation and activity-dependent processes, including endocytosis and splicing, in disease pathophysiology. We find that cell-type-specific long noncoding RNA and gene isoforms contribute substantially to regional synaptic diversity and that constrained, mutation-intolerant genes are primarily enriched in neurons. We leverage these data using an omnigenic-inspired network framework to characterize how coexpression and gene regulatory networks reflect neuropsychiatric disease risk, supporting polygenic models.

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.

Gene co-expression network analysis in human spinal cord highlights mechanisms underlying amyotrophic lateral sclerosis susceptibility.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease defined by motor neuron (MN) loss. Multiple genetic risk factors have been identified, implicating RNA and protein metabolism and intracellular transport, among other biological mechanisms. To achieve a systems-level understanding of the mechanisms governing ALS pathophysiology, we built gene co-expression networks using RNA-sequencing data from control human spinal cord samples, identifying 13 gene co-expression modules, each of which represents a distinct biological process or cell type. Analysis of four RNA-seq datasets from a range of ALS disease-associated contexts reveal dysregulation in numerous modules related to ribosomal function, wound response, and leukocyte activation, implicating astrocytes, oligodendrocytes, endothelia, and microglia in ALS pathophysiology. To identify potentially causal processes, we partitioned heritability across the genome, finding that ALS common genetic risk is enriched within two specific modules, SC.M4, representing genes related to RNA processing and gene regulation, and SC.M2, representing genes related to intracellular transport and autophagy and enriched in oligodendrocyte markers. Top hub genes of this latter module include ALS-implicated risk genes such as KPNA3, TMED2, and NCOA4, the latter of which regulates ferritin autophagy, implicating this process in ALS pathophysiology. These unbiased, genome-wide analyses confirm the utility of a systems approach to understanding the causes and drivers of ALS.

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.

Alterations in Retrotransposition, Synaptic Connectivity, and Myelination Implicated by Transcriptomic Changes Following Maternal Immune Activation in Nonhuman Primates.

BACKGROUND: Maternal immune activation (MIA) is a proposed risk factor for multiple neuropsychiatric disorders, including schizophrenia. However, the molecular mechanisms through which MIA imparts risk remain poorly understood. A recently developed nonhuman primate model of exposure to the viral mimic poly:ICLC during pregnancy shows abnormal social and repetitive behaviors and elevated striatal dopamine, a molecular hallmark of human psychosis, providing an unprecedented opportunity for studying underlying molecular correlates. METHODS: We performed RNA sequencing across psychiatrically relevant brain regions (prefrontal cortex, anterior cingulate, hippocampus) and primary visual cortex for comparison from 3.5- to 4-year-old male MIA-exposed and control offspring-an age comparable to mid adolescence in humans. RESULTS: We identify 266 unique genes differentially expressed in at least one brain region, with the greatest number observed in hippocampus. Co-expression networks identified region-specific alterations in synaptic signaling and oligodendrocytes. Although we observed temporal and regional differences, transcriptomic changes were shared across first- and second-trimester exposures, including for the top differentially expressed genes-PIWIL2 and MGARP. In addition to PIWIL2, several other regulators of retrotransposition and endogenous transposable elements were dysregulated following MIA, potentially connecting MIA to retrotransposition. CONCLUSIONS: Together, these results begin to elucidate the brain-level molecular processes through which MIA may impart risk for psychiatric disease.

Integrative genomics identifies a convergent molecular subtype that links epigenomic with transcriptomic differences in autism.

Autism spectrum disorder (ASD) is a phenotypically and genetically heterogeneous neurodevelopmental disorder. Despite this heterogeneity, previous studies have shown patterns of molecular convergence in post-mortem brain tissue from autistic subjects. Here, we integrate genome-wide measures of mRNA expression, miRNA expression, DNA methylation, and histone acetylation from ASD and control brains to identify a convergent molecular subtype of ASD with shared dysregulation across both the epigenome and transcriptome. Focusing on this convergent subtype, we substantially expand the repertoire of differentially expressed genes in ASD and identify a component of upregulated immune processes that are associated with hypomethylation. We utilize eQTL and chromosome conformation datasets to link differentially acetylated regions with their cognate genes and identify an enrichment of ASD genetic risk variants in hyperacetylated noncoding regulatory regions linked to neuronal genes. These findings help elucidate how diverse genetic risk factors converge onto specific molecular processes in ASD.

Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms.

Tissue-specific regulatory regions harbor substantial genetic risk for disease. Because brain development is a critical epoch for neuropsychiatric disease susceptibility, we characterized the genetic control of the transcriptome in 201 mid-gestational human brains, identifying 7,962 expression quantitative trait loci (eQTL) and 4,635 spliceQTL (sQTL), including several thousand prenatal-specific regulatory regions. We show that significant genetic liability for neuropsychiatric disease lies within prenatal eQTL and sQTL. Integration of eQTL and sQTL with genome-wide association studies (GWAS) via transcriptome-wide association identified dozens of novel candidate risk genes, highlighting shared and stage-specific mechanisms in schizophrenia (SCZ). Gene network analysis revealed that SCZ and autism spectrum disorder (ASD) affect distinct developmental gene co-expression modules. Yet, in each disorder, common and rare genetic variation converges within modules, which in ASD implicates superficial cortical neurons. More broadly, these data, available as a web browser and our analyses, demonstrate the genetic mechanisms by which developmental events have a widespread influence on adult anatomical and behavioral phenotypes.

Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic autism in post-mortem human brain tissue.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

Shared Molecular Neuropathology Across Major Psychiatric Disorders Parallels Polygenic Overlap.

(Gandal et al., "Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap" Science 09 Feb 2018:Vol. 359, Issue 6376, pp. 693-697 (DOI: 10.1126/science.aad6469). Reprinted with permission from AAAS).

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.

Author Correction: Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism.

Change history: In this Letter, the labels for splicing events A3SS and A5SS were swapped in column D of Supplementary Table 3a and b. This has been corrected online.

Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap.

The predisposition to neuropsychiatric disease involves a complex, polygenic, and pleiotropic genetic architecture. However, little is known about how genetic variants impart brain dysfunction or pathology. We used transcriptomic profiling as a quantitative readout of molecular brain-based phenotypes across five major psychiatric disorders-autism, schizophrenia, bipolar disorder, depression, and alcoholism-compared with matched controls. We identified patterns of shared and distinct gene-expression perturbations across these conditions. The degree of sharing of transcriptional dysregulation is related to polygenic (single-nucleotide polymorphism-based) overlap across disorders, suggesting a substantial causal genetic component. This comprehensive systems-level view of the neurobiological architecture of major neuropsychiatric illness demonstrates pathways of molecular convergence and specificity.