Search

Exploration of genome organization and function in the HIV infected brain is critical to aid in the development of treatments for HIV-associated neurocognitive disorder (HAND) and HIV cure strategies. Here, we generated a resource comprised of single nuclei transcriptomics, complemented by cell-type-specific Hi-C chromosomal conformation (‘3D genome’) and viral integration site sequencing (IS-seq) in frontal brain tissues from individuals with HIV encephalitis (HIVE), HIV-infected people without encephalitis (HIV+), and HIV uninfected (HIV-) controls. We observed profound 3D genomic reorganization of open/repressive (A/B) compartment structures encompassing 6.4% of the HIVE microglial genome that was associated with transcriptomic reprogramming, including down-regulation of homeostasis and synapse-related functions and robust activation of interferon signaling and cell migratory pathways. HIV RNA was detected in 0.003% of all nuclei in HIVE brain, predominantly in the most activated microglia where it ranked as the second most highly expressed transcript. Microglia from HIV+ brains showed, to a lesser extent, similar transcriptional alterations. IS-seq recovered 1,221 insertion events in glial nuclei that were enriched for chromosomal domains newly mobilized into a permissive chromatin environment in HIVE microglia. Brain and peripheral myeloid cell integration revealed a preference overall for transcription-permissive chromatin, but robust differences in the frequency of recurrent insertions, intergenic integration, and enrichment for pre-integration complex-associated factors at integration sites. Our resource highlights critical differences in the genomic patterns of HIV infection in brain versus blood and points to a dynamic interrelationship between inflammation-associated 3D genome remodeling and successful integration in brain.

Genetic liability associated with schizophrenia most frequently arises from common genetic variation that confers small individual effects and contributes to phenotypes only when considered in aggregate. These risk variants are typically non-coding and act by genetically regulating the expression of one or more gene targets (eGenes), but the mechanisms by which unlinked eGenes interact to contribute to complex genetic risk remains unclear. Here we apply a pooled CRISPR approach to evaluate in parallel ten schizophrenia eGenes in human glutamatergic neurons. Querying the shared neuronal impacts across eGenes reveals shared downstream transcriptomic impacts (“convergence”) concentrated on pathways of brain development and synaptic signaling. The composition and strength of convergent networks is influenced by both the similarity of eGene functional annotation and the strength of eGene co-expression in human postmortem brain tissue. Convergent networks resolve distinct patterns of eGene up-regulation associated with individual-level risk in the post-mortem dorsolateral prefrontal cortex that are broadly enriched for neuropsychiatric disorder risk genes and may represent molecular subtypes of schizophrenia. Convergent gene targets are druggable as novel points of therapeutic intervention. Overall, convergence suggests a model to explain how non-additive interactions arise between risk genes and may explain cross-disorder pleiotropy of genetic risk for psychiatric disorders.

The SET domain containing 2, histone lysine methyltransferase encoded by SETD2 is a dual-function methyltransferase for histones and microtubules and plays an important role for transcriptional regulation, genomic stability, and cytoskeletal functions. Specifically, SETD2 is associated with trimethylation of histone H3 at lysine 36 (H3K36me3) and methylation of α-tubulin at lysine 40. Heterozygous loss of function and missense variants have previously been described with Luscan-Lumish syndrome (LLS), which is characterized by overgrowth, neurodevelopmental features, and absence of overt congenital anomalies. We have identified 15 individuals with de novo variants in codon 1740 of SETD2 whose features differ from those with LLS. Group 1 consists of 12 individuals with heterozygous variant c.5218C>T p.(Arg1740Trp) and Group 2 consists of 3 individuals with heterozygous variant c.5219G>A p.(Arg1740Gln). The phenotype of Group 1 includes microcephaly, profound intellectual disability, congenital anomalies affecting several organ systems, and similar facial features. Individuals in Group 2 had moderate to severe intellectual disability, low normal head circumference, and absence of additional major congenital anomalies. While LLS is likely due to loss of function of SETD2, the clinical features seen in individuals with variants affecting codon 1740 are more severe suggesting an alternative mechanism, such as gain of function, effects on epigenetic regulation, or posttranslational modification of the cytoskeleton. Our report is a prime example of different mutations in the same gene causing diverging phenotypes and the features observed in Group 1 suggest a new clinically recognizable syndrome uniquely associated with the heterozygous variant c.5218C>T p.(Arg1740Trp) in SETD2.

Most risk variants for brain disorders identified by genome-wide association studies reside in the noncoding genome, which makes deciphering biological mechanisms difficult. A commonly used tool, multimarker analysis of genomic annotation (MAGMA), addresses this issue by aggregating single nucleotide polymorphism associations to nearest genes. Here we developed a platform, Hi-C-coupled MAGMA (H-MAGMA), that advances MAGMA by incorporating chromatin interaction profiles from human brain tissue across two developmental epochs and two brain cell types. By analyzing gene regulatory relationships in the disease-relevant tissue, H-MAGMA identified neurobiologically relevant target genes. We applied H-MAGMA to five psychiatric disorders and four neurodegenerative disorders to interrogate biological pathways, developmental windows and cell types implicated for each disorder. Psychiatric-disorder risk genes tended to be expressed during mid-gestation and in excitatory neurons, whereas neurodegenerative-disorder risk genes showed increasing expression over time and more diverse cell-type specificities. H-MAGMA adds to existing analytic frameworks to help identify the neurobiological principles of brain disorders.

As expanding genetic and genomic studies continue to implicate a growing list of variants contributing risk to neuropsychiatric disease, an important next step is to understand the functional impact and points of convergence of these risk factors. Here, with a focus on schizophrenia, we survey the most recent findings of the rare and common variants underlying genetic risk for schizophrenia. We discuss the ongoing efforts to validate these variants in post-mortem brain tissue, as well as new approaches to combine CRISPR-based genome engineering with patient-specific human induced pluripotent stem cell (hiPSC)-based models, in order to identify putative causal schizophrenia loci that regulate gene expression and cellular function. We consider the current limitations of hiPSC-based approaches as well as the future advances necessary to improve the fidelity of this human model. With the objective of utilizing patient genotype data to improve diagnosis and predict treatment response, the integration of CRISPR-genome engineering and hiPSC-based models represent an important strategy with which to systematically demonstrate the cell-type-specific effects of schizophrenia-associated variants.

Background Midbrain dopaminergic neurons (MDN) represent 0.0005% of the brain’s neuronal population and mediate cognition, food intake, and metabolism. MDN are also posited to underlay the neurobiological dysfunction of schizophrenia (SCZ), a severe neuropsychiatric disorder that is characterized by psychosis as well as multifactorial medical co-morbidities, including metabolic disease, contributing to markedly increased morbidity and mortality. Paradoxically, however, the genetic risk sequences of psychosis and traits associated with metabolic disease, such as body mass, show very limited overlap. Methods We investigated the genomic interaction of SCZ with medical conditions and traits, including body mass index (BMI), by exploring the MDN’s “spatial genome,” including chromosomal contact landscapes as a critical layer of cell type-specific epigenomic regulation. Low-input Hi-C protocols were applied to 5–10 × 103 dopaminergic and other cell-specific nuclei collected by fluorescence-activated nuclei sorting from the adult human midbrain. Results The Hi-C-reconstructed MDN spatial genome revealed 11 “Euclidean hot spots” of clustered chromatin domains harboring risk sequences for SCZ and elevated BMI. Inter- and intra-chromosomal contacts interconnecting SCZ and BMI risk sequences showed massive enrichment for brain-specific expression quantitative trait loci (eQTL), with gene ontologies, regulatory motifs and proteomic interactions related to adipogenesis and lipid regulation, dopaminergic neurogenesis and neuronal connectivity, and reward- and addiction-related pathways. Conclusions We uncovered shared nuclear topographies of cognitive and metabolic risk variants. More broadly, our PsychENCODE sponsored Hi-C study offers a novel genomic approach for the study of psychiatric and medical co-morbidities constrained by limited overlap of their respective genetic risk architectures on the linear genome.

Social isolation during the juvenile critical window is detrimental to proper functioning of the prefrontal cortex (PFC) and establishment of appropriate adult social behaviors. However, the specific circuits that undergo social experience-dependent maturation to regulate social behavior are poorly understood. We identify a specific activation pattern of parvalbumin-positive interneurons (PVIs) in dorsal-medial PFC (dmPFC) prior to an active bout, or a bout initiated by the focal mouse, but not during a passive bout when mice are explored by a stimulus mouse. Optogenetic and chemogenetic manipulation reveals that brief dmPFC-PVI activation triggers an active social approach to promote sociability. Juvenile social isolation decouples dmPFC-PVI activation from subsequent active social approach by freezing the functional maturation process of dmPFC-PVIs during the juvenile-to-adult transition. Chemogenetic activation of dmPFC-PVI activity in the adult animal mitigates juvenile isolation-induced social deficits. Therefore, social experience-dependent maturation of dmPFC-PVI is linked to long-term impacts on social behavior., Isolation during critical periods of development prevents development of normal social behaviours in mice, and this is thought to involve the prefrontal cortex. Here, the authors identify an activation pattern in parvalbumin-positive interneurons in the dorsal medial prefrontal cortex that when activated promotes sociability behaviours in mice.

The ‘non-linear’ genome, or the spatial proximity of non-contiguous sequences, emerges as an important regulatory layer for genome organization and function, including transcriptional regulation. Here, we review recent genome-scale chromosome conformation mappings (‘Hi-C’) in developing and adult human and mouse brain. Neural differentiation is associated with widespread remodeling of the chromosomal contact map, reflecting dynamic changes in cell-type-specific gene expression programs, with a massive (estimated 20–50%) net loss of chromosomal contacts that is specific for the neuronal lineage. Hi-C datasets provided an unexpected link between locus-specific abnormal expansion of repeat sequences positioned at the boundaries of self-associating topological chromatin domains and associated with monogenic neurodevelopmental and neurodegenerative disease. Furthermore, integrative analyses of cell-type-specific Hi-C and transcriptomes uncovered an expanded genomic risk space interacting with sequences conferring liability for schizophrenia and other cognitive disease. We predict that spatial genome exploration will deliver radically new insights into the brain’s nucleomes in health and disease.

Dopaminergic neurons are critical to movement, mood, addiction, and stress. Current techniques for generating dopaminergic neurons from human induced pluripotent stem cells (hiPSCs) yield heterogenous cell populations with variable purity and inconsistent reproducibility between donors, hiPSC clones, and experiments. Here, we report the rapid (5 weeks) and efficient (~90%) induction of induced dopaminergic neurons (iDANs) through transient overexpression of lineage-promoting transcription factors combined with stringent selection across five donors. We observe maturation-dependent increase in dopamine synthesis and electrophysiological properties consistent with midbrain dopaminergic neuron identity, such as slow-rising after- hyperpolarization potentials, an action potential duration of ~3 ms, tonic sub-threshold oscillatory activity, and spontaneous burst firing at a frequency of ~1.0-1.75 Hz. Transcriptome analysis reveals robust expression of genes involved in fetal midbrain dopaminergic neuron identity. Specifically expressed genes in iDANs, as well as those from isogenic induced GABAergic and glutamatergic neurons, were enriched in loci conferring heritability for cannabis use disorder, schizophrenia, and bipolar disorder; however, each neuronal subtype demonstrated subtype-specific heritability enrichments in biologically relevant pathways, and iDANs alone were uniquely enriched in autism spectrum disorder risk loci. Therefore, iDANs provide a critical tool for modeling midbrain dopaminergic neuron development and dysfunction in psychiatric disease.

Here, we mapped cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic ‘B2NeuN+’ megabase-scaling subcompartment. Comparative chromosomal conformation mapping in Mus spretus and Mus musculus revealed neuron-specific reconfigurations tracking recent ERV2 retrotransposon expansions in the murine germline, with significantly higher B2 megadomain contact frequencies at sites with ongoing ERV2 insertions in Mus musculus. Ablation of the retrotransposon silencer Kmt1e/Setdb1 triggered B2 megadomain disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly of ERV2/Intracisternal-A-Particles (IAPs) infiltrating dendrites and spines. We conclude that neuronal megadomain architectures include evolutionarily adaptive heterochromatic organization which, upon perturbation, unleashes ERV proviruses with strong tropism within mature neurons.

Non-random spatial organization of the chromosomal material inside the nuclei of brain cells emerges as an important regulatory layer of genome organization and function in health and disease. Here, we discuss how integrative approaches assessing chromatin in context of the 3D genome is providing new insights into normal and diseased neurodevelopment. Studies in primate (incl. human) and rodent brain have confirmed that chromosomal organization in neurons and glia undergoes highly dynamic changes during pre- and early postnatal development, with potential for plasticity across a much wider age window. For example, neuronal 3D genomes from juvenile and adult cerebral cortex and hippocampus undergo chromosomal conformation changes at hundreds of loci in the context of learning and environmental enrichment, viral infection, and neuroinflammation. Furthermore, locus-specific structural DNA variations, such as micro-deletions, duplications, repeat expansions, and retroelement insertions carry the potential to disrupt the broader epigenomic and transcriptional landscape far beyond the boundaries of the site-specific variation, highlighting the critical importance of long-range intra- and inter-chromosomal contacts for neuronal and glial function.

SUMMARYPluripotent stem cell-derived cardiomyocytes (PSC-CMs) provide an unprecedented opportunity to study human heart development and disease. A major caveat however is that they remain functionally and structurally immature in culture, limiting their potential for disease modeling and regenerative approaches. Here, we address the question of how different metabolic pathways can be modulated in order to induce efficient hPSC-CM maturation. We show that PPAR signaling acts in an isoform-specific manner to balance glycolysis and fatty acid oxidation (FAO). PPARD activation or inhibition results in efficient respective up- or down-regulation of the gene regulatory networks underlying FAO in hPSC-CMs. PPARD induction further increases mitochondrial and peroxisome content, enhances mitochondrial cristae formation and augments FAO flux. Lastly PPARD activation results in enhanced myofibril organization and improved contractility. Transient lactate exposure, commonly used in hPSC-CM purification protocols, induces an independent program of cardiac maturation, but when combined with PPARD activation equally results in a metabolic switch to FAO. In summary, we identify multiple axes of metabolic modifications of hPSC-CMs and a role for PPARD signaling in inducing the metabolic switch to FAO in hPSC-CMs. Our findings provide new and easily implemented opportunities to generate mature hPSC-CMs for disease modeling and regenerative therapy.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Overproduction or suppression of certain microRNAs (miRNAs) in Alzheimer’s disease (AD) brains promote alterations in tau proteostasis and neurodegeneration. However, the mechanisms governing how specific miRNAs are dysregulated in AD brains are still under investigation. Epitranscriptomic regulation adds a layer of post-transcriptional control to brain function during development and adulthood. NOP2/Sun RNA methyltransferase 2 (NSun2) is one of the few known brain-enriched methyltransferases able to modify mammalian non-coding RNAs and loss of function autosomal-recessive mutations in NSUN2 have been associated with neurological abnormalities in humans. Here, we provide evidence that NSun2 is expressed in adult human neurons in the hippocampal formation and prefrontal cortex. When we evaluated NSun2 protein expression in post-mortem brain tissue from AD patients we find is dysregulated which was also found in mice and human cellular AD models. To probe these observed alterations were unique to AD we further evaluated brain tissue from other tauopathies, observing NSun2 protein levels were similar between cases and controls. In a well-established Drosophila melanogaster model of tau-induced toxicity we investigated the pathological role of NSun2 observing that reduction of NSun2 protein levels exacerbated tau toxicity, while overexpression of NSun2 partially abrogated the toxic effects. We further show using human induced pluripotent stem cell (iPSC) derived neuronal cultures that NSun2 deficiency results in tau hyperphosphorylation and we found in primary hippocampal neuronal cultures NSun2 levels decrease in response to amyloid-beta oligomers (AβO). Furthermore, in mice, we observed that NSun2 deficiency promotes aberrant levels of m6A methylated miR-125b and tau hyperphosphorylation. Altogether, our study supports that neuronal NSun2 deficiency in AD promotes neurodegeneration by altering tau phosphorylation and tau toxicity through an epitranscriptomic regulatory mechanism and highlights a novel avenue for therapeutic targeting.

Opioid dependence is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. We interrogated the effects of opioid dependence on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicts case-control status with high accuracy. We focus on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or “plexus”: ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid dependence and captures genetic and environmental factors perpetuating the opioid epidemic.

BackgroundRecent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes and readily available compounds that reduce COVID-19 host susceptibility is a critical next step.MethodsWe integrate COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX) and perturbargen signatures to identify candidate genes and compounds that reverse the predicted gene expression dysregulation associated with COVID-19 susceptibility. The top candidate gene is validated by testing both its GReX and observed blood transcriptome association with COVID-19 severity, as well as by in vitro perturbation to quantify effects on viral load and molecular pathway dysregulation. We validate the in silico drug repositioning analysis by examining whether the top candidate compounds decrease COVID-19 incidence based on epidemiological evidence.ResultsWe identify IL10RB as the top key regulator of COVID-19 host susceptibility. Predicted GReX up-regulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes. In vitro IL10RB overexpression is associated with increased viral load and activation of immune-related molecular pathways. Azathioprine and retinol are prioritized as candidate compounds to reduce the likelihood of testing positive for COVID-19.ConclusionsWe establish an integrative data-driven approach for gene target prioritization. We identify and validate IL10RB as a suitable molecular target for modulation of COVID-19 host susceptibility. Finally, we provide evidence for a few readily available medications that would warrant further investigation as drug repositioning candidates.

To explore modular organization of chromosomes in schizophrenia (SCZ) and bipolar disorder (BD), we applied ‘population-scale’ correlational structuring of 739 histone H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from the prefrontal cortex (PFC) of 568 cases and controls. Neuronal histone acetylomes and methylomes assembled as thousands of cis-regulatory domains (CRDs), revealing fine-grained, kilo-to megabase scale chromatin organization at higher resolution but firmly integrated into Hi-C chromosomal conformations. Large clusters of domains that were hyperacetylated in disease shared spatial positioning within the nucleus, predominantly regulating PFC projection neuron function and excitatory neurotransmission. Hypoacetylated domains were linked to inhibitory interneuron- and myelination-relevant genes. Chromosomal modular architecture is affected in SCZ and BD, with hyperacetylated domains showing unexpectedly strong convergences defined by cell type, nuclear topography, genetic risk, and active chromatin state across a wide developmental window.

Human prefrontal cortex (PFC) is associated with broad individual variabilities in functions linked to personality, social behaviors, and cognitive functions. The phenotype variabilities associated with brain functions can be caused by genetic or epigenetic factors. The interactions between these factors in human subjects is, as of yet, poorly understood. The heterogeneity of cerebral tissue, consisting of neuronal and nonneuronal cells, complicates the comparative analysis of gene activities in brain specimens. To approach the underlying neurogenomic determinants, we performed a deep analysis of open chromatin–associated histone methylation in PFC neurons sorted from multiple human individuals in conjunction with whole-genome and transcriptome sequencing. Integrative analyses produced novel unannotated neuronal genes and revealed individual-specific chromatin “blueprints” of neurons that, in part, relate to genetic background. Surprisingly, we observed gender-dependent epigenetic signals, implying that gender may contribute to the chromatin variabilities in neurons. Finally, we found epigenetic, allele-specific activation of the testis-specific gene nucleoporin 210 like (NUP210L) in brain in some individuals, which we link to a genetic variant occurring in

Schizophrenia is a debilitating psychiatric disorder with a complex genetic architecture and limited understanding of its neuropathology, reflected by the lack of diagnostic measures and effective pharmacological treatments. Geneticists have recently identified more than 145 risk loci comprising hundreds of common variants of small effect sizes, most of which lie in noncoding genomic regions. This review will discuss how the epigenetic toolbox can be applied to contextualize genetic findings in schizophrenia. Progress in next-generation sequencing, along with increasing methodological complexity, has led to the compilation of genome-wide maps of DNA methylation, histone modifications, RNA expression, and more. Integration of chromatin conformation datasets is one of the latest efforts in deciphering schizophrenia risk, allowing the identification of genes in contact with regulatory variants across 100s of kilobases. Large-scale multiomics studies will facilitate the prioritization of putative causal risk variants and gene networks that contribute to schizophrenia etiology, informing clinical diagnostics and treatment downstream. .La esquizofrenia es un trastorno psiquiátrico invalidante con una arquitectura genética compleja y una comprensión limitada de su neuropatología, lo que se traduce en la ausencia de mediciones diagnósticas y tratamientos farmacológicos eficaces. Recientemente los genetistas han identificado más de 145 loci de riesgo que comprenden cientos de variantes comunes de pequeño efecto, la mayoría de los cuales se encuentran en regiones genómicas no codificantes. Este artículo revisará cómo se pueden aplicar las herramientas epigenéticas para contextualizar los hallazgos genéticos en la esquizofrenia. El progreso en la secuenciación de segunda generación, junto con la creciente complejidad metodológica, ha llevado a la recopilación de mapas de metilación de ADN, modificaciones de histonas, expresión de ARN y otros más. Uno de los últimos esfuerzos para descifrar el riesgo de presentar esquizofrenia es la integración de los conjuntos de datos de conformación de cromatina, lo que permite la identificación de genes en contacto con variantes reguladoras en cientos de kilobases. Los estudios multiómicos a gran escala facilitarán la priorización de las posibles variantes de riesgo causal y las redes de genes que contribuyen a la etiología de la esquizofrenia, lo que incidirá en el diagnóstico clínico y en el tratamiento.La schizophrénie est une maladie psychiatrique invalidante dont l’architecture génétique est complexe : sa neuropathologie est mal comprise comme le montre le manque de mesures diagnostiques et de traitements médicamenteux efficaces. Les généticiens ont récemment identifié plus de 145 loci de risque comprenant des centaines de variants courants de faible effet, la plupart situés dans des régions non codantes du génome. Cet article discute de la façon d’appliquer les outils épigénétiques pour replacer les résultats génétiques dans le contexte de la schizophrénie. Les avancées du séquençage de nouvelle génération associées à une méthodologie de plus en plus complexe, ont permis de regrouper les cartes pangénomiques de la méthylation de l’ADN, des modifications de l’histone, de l’expression de l’ARN et plus encore. L’un des derniers progrès en date dans le déchiffrage du risque schizophrénique porte sur l’intégration de séries de données sur l’architecture de la chromatine, permettant l’identification de gènes en contact avec les variants de régulation parmi des centaines de kilobases. Des études multiomiques à grande échelle aideront à prioriser les variants du risque causal putatif et les réseaux de gènes qui contribuent à l’étiologie de la schizophrénie, pour renseigner en aval les diagnostics cliniques et les traitements.

The complexity of human astrocytes remains poorly defined in primary human tissue, requiring better tools for their isolation and molecular characterization. Fluorescence-activated nuclei sorting (FANS) can be used to successfully isolate and study human neuronal nuclei (NeuN+) populations from frozen archival tissue, thereby avoiding problems associated with handling fresh tissue. However, efforts to similarly isolate astroglia from the non-neuronal (NeuN−) element are lacking. A recently developed and validated immunotagging strategy uses three transcription factor antibodies to simultaneously isolate enriched neuronal (NeuN+), astrocyte (paired box protein 6 (PAX6)+NeuN−), and oligodendrocyte progenitor (OLIG2+NeuN−) nuclei populations from non-diseased, fresh (unfixed) snap-frozen postmortem human temporal neocortex tissue. This technique was shown to be useful for the characterization of cell type-specific transcriptome alterations in primary pathological epilepsy neocortex. Transcriptomic analyses confirmed that PAX6+NeuN− sorted populations are robustly enriched for pan-astrocyte markers and capture astrocytes in both resting and reactive conditions. This paper describes the FANS methodology for the isolation of astrocyte-enriched nuclei populations from fresh-frozen human cortex, including tissue dissociation into single-nucleus (sn) suspension; immunotagging of nuclei with anti-NeuN and anti-PAX6 fluorescently conjugated antibodies; FANS gating strategies and quality control metrics for optimizing sensitivity and specificity during sorting and for confirming astrocyte enrichment; and recommended procurement for downstream transcriptome and chromatin accessibility sequencing at bulk or sn resolution. This protocol is applicable for non-necrotic, fresh-frozen, human cortical specimens with various pathologies and recommended postmortem tissue collection within 24 h.

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration—particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs., This report discusses gaps in knowledge and possibilities for the next phase of functional discovery for addiction.

Cigarette smoking and alcohol use are among the most prevalent substances used worldwide and account for a substantial proportion of preventable morbidity and mortality, underscoring the public health significance of understanding their etiology. Genome-wide association studies (GWAS) have successfully identified genetic variants associated with cigarette smoking and alcohol use traits. However, the vast majority of risk variants reside in non-coding regions of the genome, and their target genes and neurobiological mechanisms are unknown. Chromosomal conformation mappings can address this knowledge gap by charting the interaction profiles of risk-associated regulatory variants with target genes. To investigate the functional impact of common variants associated with cigarette smoking and alcohol use traits, we applied Hi-C coupled MAGMA (H-MAGMA) built upon cortical and midbrain dopaminergic neuronal Hi-C datasets to GWAS summary statistics of nicotine dependence, cigarettes per day, problematic alcohol use, and drinks per week. The identified risk genes mapped to key pathways associated with cigarette smoking and alcohol use traits, including drug metabolic processes and neuronal apoptosis. Risk genes were highly expressed in cortical glutamatergic, midbrain dopaminergic, GABAergic, and serotonergic neurons, suggesting them as relevant cell types in understanding the mechanisms by which genetic risk factors influence cigarette smoking and alcohol use. Lastly, we identified pleiotropic genes between cigarette smoking and alcohol use traits under the assumption that they may reveal substance-agnostic, shared neurobiological mechanisms of addiction. The number of pleiotropic genes was ∼26-fold higher in dopaminergic neurons than in cortical neurons, emphasizing the critical role of ascending dopaminergic pathways in mediating general addiction phenotypes. Collectively, brain region- and neuronal subtype-specific 3D genome architecture refines neurobiological hypotheses for smoking, alcohol, and general addiction phenotypes by linking genetic risk factors to their target genes.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant interindividual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly afflict healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants influence vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single-nucleotide polymorphism (rs4702), common in the population and located in the 3′ UTR of the protease FURIN, influences alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can have an impact on viral infection and thus contribute to clinical heterogeneity in COVID-19. Ongoing genetic studies will help to identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs., Graphical abstract, Brennand, Akbarian, and colleagues combine human induced pluripotent stem cell-based models and CRISPR engineering to explore the host genetics of SARS-CoV-2 in alveolar cells, intestinal cells, and neurons. A SNP(rs4702) located in the 3′ UTR of the protease FURIN influences alveolar and neuron infection by SARS-CoV-2 in vitro.

SummaryIn early development, the environment triggers mnemonic epigenomic programs resulting in memory and learning experiences to confer cognitive phenotypes into adulthood. To uncover how environmental stimulation impacts the epigenome and genome organization, we used the paradigm of environmental enrichment (EE) in young mice constantly receiving novel stimulation. We profiled epigenome and chromatin architecture in whole cortex and sorted neurons by deep-sequencing techniques. Specifically, we studied chromatin accessibility, gene and protein regulation, and 3D genome conformation, combined with predicted enhancer and chromatin interactions. We identified increased chromatin accessibility, transcription factor binding including CTCF-mediated insulation, differential occupancy of H3K36me3 and H3K79me2, and changes in transcriptional programs required for neuronal development. EE stimuli led to local genome re-organization by inducing increased contacts between chromosomes 7 and 17 (inter-chromosomal). Our findings support the notion that EE-induced learning and memory processes are directly associated with the epigenome and genome organization.Highlights-Environmental enrichment (EE) alters chromatin conformation, CTCF binding, and spatially 3D genome changes, thereby regulating cognitive function during the first steps of life after birth.-Transcription-associated gene body marks H3K79me2 and H3K36me3 are differently influenced by EE in cortical brain cells and binding is exacerbated upon stimulation in an age-dependent manner.-EE-induced changes of 3D genome organization increase inter-chromosomal interactions of genes associated with synaptic transmission and AMPA receptor genes on chromosomes 7 and 17.

In coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, the relationship between brain tropism, neuroinflammation and host immune response has not been well characterized. We analyzed 68,557 single-nucleus transcriptomes from three brain regions (dorsolateral prefrontal cortex, medulla oblongata and choroid plexus) and identified an increased proportion of stromal cells and monocytes in the choroid plexus of COVID-19 patients. Differential gene expression, pseudo-temporal trajectory and gene regulatory network analyses revealed microglial transcriptome perturbations, mediating a range of biological processes, including cellular activation, mobility and phagocytosis. Quantification of viral spike S1 protein and SARS-CoV-2 transcripts did not support the notion of brain tropism. Overall, our findings suggest extensive neuroinflammation in patients with acute COVID-19.One Sentence SummarySingle-nucleus transcriptome analysis suggests extensive neuroinflammation in human brain tissue of patients with acute coronavirus disease 2019.

Additional file 1:. Figure S1-S13, Table S1,S3,S4.

Additional file 5. BSMN_Membership_List.

As the complex genetics underlying risk for psychiatric disorders are uncovered, the translation of genetic findings to inform our understanding of disease pathophysiology and treatment is an urgent and important challenge. Human induced pluripotent stem cell (hiPSC)-based models provide a source of live human donor-specific neurons and glia, while CRISPR-based methods make it possible to precisely engineer and manipulate the genome and epigenome. Together, these emerging technologies make possible the functional exploration of the fundamental roles of risk variants and genes linked to psychiatric disease. Here, we provide a brief overview of the genetic risk architecture associated with psychiatric disorders, highlight relevant technical advancements in hiPSC models and CRISPR engineering that make possible new avenues to empirically evaluate the functional impact of risk variants in disease-relevant cell types. We next assess novel applications combining hiPSC-based approaches and genetic screening to uncover the effect of risk variants at unparalleled scale and culminate with a critical appraisal of the developments necessary to translate new insights into clinically actionable information and novel therapeutic approaches.

Additional file 6. Review history.

Pluripotent stem-cell-derived cardiomyocytes (PSC-CMs) provide an unprecedented opportunity to study human heart development and disease, but they are functionally and structurally immature. Here, we induce efficient human PSC-CM (hPSC-CM) maturation through metabolic-pathway modulations. Specifically, we find that peroxisome-proliferator-associated receptor (PPAR) signaling regulates glycolysis and fatty acid oxidation (FAO) in an isoform-specific manner. While PPARalpha (PPARa) is the most active isoform in hPSC-CMs, PPARdelta (PPARd) activation efficiently upregulates the gene regulatory networks underlying FAO, increases mitochondrial and peroxisome content, enhances mitochondrial cristae formation, and augments FAO flux. PPARd activation further increases binucleation, enhances myofibril organization, and improves contractility. Transient lactate exposure, which is frequently used for hPSC-CM purification, induces an independent cardiac maturation program but, when combined with PPARd activation, still enhances oxidative metabolism. In summary, we investigate multiple metabolic modifications in hPSC-CMs and identify a role for PPARd signaling in inducing the metabolic switch from glycolysis to FAO in hPSC-CMs.

Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. Here, the Brain Somatic Mosaicism Network conducted a coordinated, multi-institutional study to: (i) examine the ability of existing methods to detect simulated somatic single nucleotide variants (SNVs) in DNA mixing experiments; (ii) generate multiple replicates of whole genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual; (iii) devise strategies to discover somatic SNVs; and (iv) apply various approaches to validate somatic SNVs. These efforts led to the identification of 43 bona fide somatic SNVs that ranged in variant allele fractions from ~0.005 to ~0.28. Guided by these results, we devised best practices for calling mosaic SNVs from 250X whole genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrated that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. Thus, this study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Epigenetic decoration of the genome is a highly regulated cell-specific and life-long process including DNA methylation, many different types of covalent and posttranslational histone modifications, and complex nonrandom chromosomal conformations including “loopings” and self-folding “domains.” These epigenomic determinants remain “plastic” throughout all periods of development and aging, with ongoing dynamic regulation even in neurons and other highly differentiated brain cells. There can be little doubt that DNA methylation landscapes show substantial reorganization during the course of normal development and aging of the human cerebral cortex and participate in alterations in gene expression with disease. With the underlying DNA sequence left intact, it is possible to imagine that simultaneous epigenomic targeting of enhancer and promoter sequences within multiple risk haplotypes could offer a promising approach to effectively alter cognition and behavior.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant inter-individual variability. In addition to showing more disease in males, the elderly, and individuals with underlying co-morbidities, SARS-CoV-2 can seemingly render healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants also impact vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR-engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single nucleotide polymorphism (rs4702), common in the population at large, and located in the 3’UTR of the protease FURIN, impacts alveolar and neuron infection by SARS-CoV-2in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can impact viral infection, and thus contribute to clinical heterogeneity in SARS-CoV-2. Ongoing genetic studies will help to better identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs that might treat disease.

Schizophrenia is a chronic and severe neuropsychiatric condition manifested by cognitive, emotional, affective, perceptual, and behavioral abnormalities. Despite decades of research, the biological substrates driving the signs and symptoms of the disorder remain elusive, thus hampering progress in the development of treatments aimed at disease etiologies. The recent emergence of human induced pluripotent stem cell (hiPSC)-based models has provided the field with a highly innovative approach to generate, study, and manipulate living neural tissue derived from patients, making possible the exploration of fundamental roles of genes and early-life stressors in disease-relevant cell types. Here, we begin with a brief overview of the clinical, epidemiological, and genetic aspects of the condition, with a focus on schizophrenia as a neurodevelopmental disorder. We then highlight relevant technical advancements in hiPSC models and assess novel findings attained using hiPSC-based approaches and their implications for disease biology and treatment innovation. We close with a critical appraisal of the developments necessary for both further expanding knowledge of schizophrenia and the translation of new insights into therapeutic innovations.

Cellular heterogeneity in the human brain obscures the identification of robust cellular regulatory networks. Here we integrated genome-wide chromosome conformation in purified neurons and glia with transcriptomic and enhancer profiles to build the gene regulatory landscape of two major cell classes in the human brain. Within glutamatergic and GABAergic neurons, we were able to link enhancers to their cognate genes via neuronal chromatin interaction profiles. These cell-type-specific regulatory landscapes were then leveraged to gain insight into the cellular etiology of several brain disorders. We found that Alzheimer’s disease (AD)-associated epigenetic dysregulation was linked to neurons and oligodendrocytes, whereas genetic risk factors for AD highlighted microglia as a central cell type, suggesting that different cell types may confer risk to the disease via different genetic mechanisms. Moreover, neuronal subtype-specific annotation of genetic risk factors for schizophrenia and bipolar disorder identified shared (parvalbumin-expressing interneurons) and distinct cellular etiology (upper layer neurons for bipolar and deeper layer projection neurons for schizophrenia) between these two closely related psychiatric illnesses. Collectively, these findings shed new light on cell-type-specific gene regulatory networks in brain disorders.

Cellular heterogeneity in the human brain obscures the identification of robust cellular regulatory networks, which is necessary to understand the function of non-coding elements and the impact of non-coding genetic variation. Here we integrate genome-wide chromosome conformation data from purified neurons and glia with transcriptomic and enhancer profiles, to characterize the gene regulatory landscape of two major cell classes in the human brain. We then leverage cell-type-specific regulatory landscapes to gain insight into the cellular etiology of several brain disorders. We find that Alzheimer’s disease (AD)-associated epigenetic dysregulation is linked to neurons and oligodendrocytes, whereas genetic risk factors for AD highlighted microglia, suggesting that different cell types may contribute to disease risk, via different mechanisms. Moreover, integration of glutamatergic and GABAergic regulatory maps with genetic risk factors for schizophrenia (SCZ) and bipolar disorder (BD) identifies shared (parvalbumin-expressing interneurons) and distinct cellular etiologies (upper layer neurons for BD, and deeper layer projection neurons for SCZ). Collectively, these findings shed new light on cell-type-specific gene regulatory networks in brain disorders., The cellular heterogeneity in brain obscures the identification of robust cellular regulatory networks. Here the authors integrate genome-wide chromosome conformation data from sorted neurons and glia, with transcriptomic and enhancer profiles, to characterize cell-type-specific gene regulatory landscapes in the human brain, and provide insights into cell-type-specific gene regulatory networks in brain disorders.

Juvenile social isolation reduces sociability in adulthood, but the underlying neural circuit mechanisms are poorly understood. We found that, in male mice, 2 weeks of social isolation immediately following weaning leads to a failure to activate medial prefrontal cortex neurons projecting to the posterior paraventricular thalamus (mPFC→pPVT) during social exposure in adulthood. Chemogenetic or optogenetic suppression of mPFC→pPVT activity in adulthood was sufficient to induce sociability deficits without affecting anxiety-related behaviors or preference toward rewarding food. Juvenile isolation led to both reduced excitability of mPFC→pPVT neurons and increased inhibitory input drive from low-threshold-spiking somatostatin interneurons in adulthood, suggesting a circuit mechanism underlying sociability deficits. Chemogenetic or optogenetic stimulation of mPFC→pPVT neurons in adulthood could rescue the sociability deficits caused by juvenile isolation. Our study identifies a pair of specific medial prefrontal cortex excitatory and inhibitory neuron populations required for sociability that are profoundly affected by juvenile social experience. Yamamuro et al. show that juvenile social isolation disrupts prefrontal neurons projecting to the paraventricular thalamus and associated prefrontal somatostatin interneurons, and thereby impairs sociability in adulthood.

BACKGROUND: Schizophrenia (SZ) is associated with increased all-cause mortality, smoking, and age-associated proteins, yet multiple previous studies found no association between SZ and biological age using Horvath’s epigenetic clock, a well-established aging biomarker based on DNA methylation (DNAm). However, numerous epigenetic clocks have been developed that may capture distinct aspects of aging. This study tested the hypothesis that altered aging in SZ manifests in these other clocks. METHODS: We performed a comprehensive analysis of 14 epigenetic clocks, categorized according to what they were trained to predict: chronological age, mortality, mitotic divisions, or telomere length. To understand the etiology of biological age differences, we also examined DNAm predictors of smoking, alcohol, BMI, serum proteins, and cell proportions. We independently analyzed three publicly available multi-ethnic DNAm datasets from whole blood, a total of 567 SZ cases and 594 non-psychiatric controls. RESULTS: All data sets showed accelerations in SZ for the three mortality clocks up to 5 years, driven by smoking and elevated levels of 6 age-associated proteins. The two mitotic clocks were decelerated in SZ related to anti-tumor NK and CD8T cells, which may help explain conflicting reports about low cancer rates in epidemiological studies of SZ. One cohort with available medication data showed clozapine is associated with male-specific decelerations up to 7 years in multiple chronological age clocks. CONCLUSIONS: Our study demonstrates the utility of studying the various epigenetic clocks in tandem and highlights potential mechanisms by which mental illness influences long-term outcomes including cancer and early mortality.

Causal genes and variants within genome-wide association study (GWAS) loci can be identified by integrating GWAS statistics with expression quantitative trait loci (eQTL) and determining which variants underlie both GWAS and eQTL signals. Most analyses, however, consider only the marginal eQTL signal, rather than dissect this signal into multiple conditionally independent signals for each gene. Here we show that analyzing conditional eQTL signatures, which could be important under specific cellular or temporal contexts, leads to improved fine mapping of GWAS associations. Using genotypes and gene expression levels from post-mortem human brain samples (n = 467) reported by the CommonMind Consortium (CMC), we find that conditional eQTL are widespread; 63% of genes with primary eQTL also have conditional eQTL. In addition, genomic features associated with conditional eQTL are consistent with context-specific (e.g., tissue-, cell type-, or developmental time point-specific) regulation of gene expression. Integrating the 2014 Psychiatric Genomics Consortium schizophrenia (SCZ) GWAS and CMC primary and conditional eQTL data reveals 40 loci with strong evidence for co-localization (posterior probability > 0.8), including six loci with co-localization of conditional eQTL. Our co-localization analyses support previously reported genes, identify novel genes associated with schizophrenia risk, and provide specific hypotheses for their functional follow-up.

Large-scale consortia including the Psychiatric Genomics Consortium, the Common Minds Consortium, BrainSeq and PsychENCODE, and many other studies taken together provide increasingly detailed insights into the genetic and epigenetic risk architectures of schizophrenia and offer vast amounts of molecular information, but with largely unexplored therapeutic potential. Here we discuss how epigenomic studies in human brain could guide animal work to test the impact of disease-associated alterations in chromatin structure and function on cognition and behavior. For example, transcription factors such as MYOCYTE-SPECIFIC ENHANCER FACTOR (MEF2C), or multiple regulators of the open chromatin mark, methyl-histone H3-lysine 4, are associated with the genetic risk architectures of common psychiatric disease and alterations in chromatin structure and function in diseased brain tissue. Importantly, these molecules also affect cognition and behavior in genetically engineered mice, including virus-mediated expression changes in prefrontal cortex and other key nodes in the circuitry underlying psychosis. Therefore, preclinical and small laboratory animal work could target genomic sequences affected by chromatin alterations in schizophrenia. To this end, in vivo editing of enhancer and other regulatory non-coding DNA by RNA-guided nucleases including CRISPR-Cas, and designer transcription factors, could be expected to deliver pipelines for novel therapeutic approaches aimed at improving cognitive dysfunction and other core symptoms of schizophrenia.

CRISPR/Cas9 technology has transformed our abilities to manipulate the genome and epigenome, as applications have expanded from efficient genomic editing to include the targeted localization of effectors to specific genomic loci. By facilitating the manipulation of DNA- and histone-modifying enzyme activities, activation or repression of gene expression, and targeting of transcriptional regulators to defined loci, it is now possible to directly probe the role of gene-regulatory and epigenetic pathways in order to examine their roles in basic biology and disease processes. Here we discuss these emerging CRISPR-based methodologies, with specific consideration of neurobiological applications using human induced pluripotent stem cell (hiPSC)-based models.

PHF8 is a histone demethylase with specificity for repressive modifications. While mutations of PHF8 have been associated with cognitive defects and cleft lip/palate, its role in mammalian development and physiology remains unexplored. Here, we have generated a Phf8 knockout allele in mice to examine the consequences of Phf8 loss for development and behaviour. Phf8 deficient mice neither display obvious developmental defects nor signs of cognitive impairment. However, we report a striking resiliency to stress-induced anxiety- and depression-like behaviour on loss of Phf8. We further observe misregulation of serotonin signalling within the prefrontal cortex of Phf8 deficient mice and identify the serotonin receptors Htr1a and Htr2a as direct targets of PHF8. Our results clarify the functional role of Phf8 in mammalian development and behaviour and establish a direct link between Phf8 expression and serotonin signalling, identifying this histone demethylase as a potential target for the treatment of anxiety and depression., Mutation of the human gene PHF8, encoding a histone demethylase, is linked to cognitive defects but its role in development is unclear. Here, the authors show that Phf8 deletion in mice causes no overt developmental defects but confers resilience to depression, likely through increased serotonin signalling.

Both heritability and environment contribute to risk for schizophrenia. However, the molecular mechanisms of interactions between genetic and non-genetic factors remain unclear. Epigenetic regulation of neuronal genome may be a presumable mechanism in pathogenesis of schizophrenia. Here, we performed analysis of open chromatin landscape of gene promoters in prefrontal cortical (PFC) neurons from schizophrenic patients. We cataloged cell-type-based epigenetic signals of transcriptional start sites (TSS) marked by histone H3-K4 trimethylation (H3K4me3) across the genome in PFC from multiple schizophrenia subjects and age-matched control individuals. One of the top-ranked chromatin alterations was found in the major histocompatibility (MHC) locus on chromosome 6 highlighting the overlap between genetic and epigenetic risk factors in schizophrenia. The chromosome conformation capture (3C) analysis in human brain cells revealed the architecture of multipoint chromatin interactions between the schizophrenia-associated genetic and epigenetic polymorphic sites and distantly located HLA-DRB5 and BTNL2 genes. In addition, schizophrenia-specific chromatin modifications in neurons were particularly prominent for non-coding RNA genes, including an uncharacterized LINC01115 gene and recently identified BNRNA_052780. Notably, protein-coding genes with altered epigenetic state in schizophrenia are enriched for oxidative stress and cell motility pathways. Our results imply the rare individual epigenetic alterations in brain neurons are involved in the pathogenesis of schizophrenia.