Schizophrenia Risk Mapping and Functional Engineering of the 3D Genome in Three Neuronal Subtypes.

Common variants associated with schizophrenia are concentrated in non-coding regulatory sequences, but their precise target genes are context-dependent and impacted by cell-type-specific three-dimensional spatial chromatin organization. Here, we map long-range chromosomal conformations in isogenic human dopaminergic, GABAergic, and glutamatergic neurons to track developmentally programmed shifts in the regulatory activity of schizophrenia risk loci. Massive repressive compartmentalization, concomitant with the emergence of hundreds of neuron-specific multi-valent chromatin architectural stripes, occurs during neuronal differentiation, with genes interconnected to genetic risk loci through these long-range chromatin structures differing in their biological roles from genes more proximal to sequences conferring heritable risk. Chemically induced CRISPR-guided chromosomal loop-engineering for the proximal risk gene SNAP91 and distal risk gene BHLHE22 profoundly alters synaptic development and functional activity. Our findings highlight the large-scale cell-type-specific reorganization of chromosomal conformations at schizophrenia risk loci during neurodevelopment and establish a causal link between risk-associated gene-regulatory loop structures and neuronal function.

R-loop landscapes in the developing human brain are linked to neural differentiation and cell-type specific transcription.

Here, we construct genome-scale maps for R-loops, three-stranded nucleic acid structures comprised of a DNA/RNA hybrid and a displaced single strand of DNA, in the proliferative and differentiated zones of the human prenatal brain. We show that R-loops are abundant in the progenitor-rich germinal matrix, with preferential formation at promoters slated for upregulated expression at later stages of differentiation, including numerous neurodevelopmental risk genes. RNase H1-mediated contraction of the genomic R-loop space in neural progenitors shifted differentiation toward the neuronal lineage and was associated with transcriptomic alterations and defective functional and structural neuronal connectivity in vivo and in vitro. Therefore, R-loops are important for fine-tuning differentiation-sensitive gene expression programs of neural progenitor cells.

Induction of dopaminergic neurons for neuronal subtype-specific modeling of psychiatric disease risk.

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.

Identifying synergistic high-order 3D chromatin conformations from genome-scale nanopore concatemer sequencing.

High-order three-dimensional (3D) interactions between more than two genomic loci are common in human chromatin, but their role in gene regulation is unclear. Previous high-order 3D chromatin assays either measure distant interactions across the genome or proximal interactions at selected targets. To address this gap, we developed Pore-C, which combines chromatin conformation capture with nanopore sequencing of concatemers to profile proximal high-order chromatin contacts at the genome scale. We also developed the statistical method Chromunity to identify sets of genomic loci with frequencies of high-order contacts significantly higher than background ('synergies'). Applying these methods to human cell lines, we found that synergies were enriched in enhancers and promoters in active chromatin and in highly transcribed and lineage-defining genes. In prostate cancer cells, these included binding sites of androgen-driven transcription factors and the promoters of androgen-regulated genes. Concatemers of high-order contacts in highly expressed genes were demethylated relative to pairwise contacts at the same loci. Synergies in breast cancer cells were associated with tyfonas, a class of complex DNA amplicons. These results rigorously link genome-wide high-order 3D interactions to lineage-defining transcriptional programs and establish Pore-C and Chromunity as scalable approaches to assess high-order genome structure.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains.

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

The SEQC2 epigenomics quality control (EpiQC) study.

BACKGROUND: Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group. RESULTS: Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms. CONCLUSIONS: The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Toll-like Receptor (TLR)-induced Rasgef1b expression in macrophages is regulated by NF-κB through its proximal promoter.

Ras Guanine Exchange Factor (RasGEF) domain family member 1b is encoded by a Toll-like receptor (TLR)-inducible gene expressed in macrophages, but transcriptional mechanisms that govern its expression are still unknown. Here, we have functionally characterized the 5' flanking Rasgef1b sequence and analyzed its transcriptional activation. We have identified that the inflammation-responsive promoter is contained within a short sequence (-183 to +119) surrounding the transcriptional start site. The promoter sequence is evolutionarily conserved and harbors a cluster of five NF-κB binding sites. Luciferase reporter gene assay showed that the promoter is responsive to TLR activation and RelA or cRel, but not RelB, transcription factors. Besides, site-directed mutagenesis showed that the κB binding sites are required for maximal promoter activation induced by LPS. Analysis by Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) revealed that the promoter is located in an accessible chromatin region. More important, Chromatin Immunoprecipitation sequencing (ChIP-seq) showed that RelA is recruited to the promoter region upon LPS stimulation of bone marrow-derived macrophages. Finally, studies with Rela-deficient macrophages or pharmacological inhibition by Bay11-7082 showed that NF-κB is required for optimal Rasgef1b expression induced by TLR agonists. Our data provide evidence of the regulatory mechanism mediated by NF-κB that facilitates Rasgef1b expression after TLR activation in macrophages.

Genome-wide cell-free DNA mutational integration enables ultra-sensitive cancer monitoring.

In many areas of oncology, we lack sensitive tools to track low-burden disease. Although cell-free DNA (cfDNA) shows promise in detecting cancer mutations, we found that the combination of low tumor fraction (TF) and limited number of DNA fragments restricts low-disease-burden monitoring through the prevailing deep targeted sequencing paradigm. We reasoned that breadth may supplant depth of sequencing to overcome the barrier of cfDNA abundance. Whole-genome sequencing (WGS) of cfDNA allowed ultra-sensitive detection, capitalizing on the cumulative signal of thousands of somatic mutations observed in solid malignancies, with TF detection sensitivity as low as 10-5. The WGS approach enabled dynamic tumor burden tracking and postoperative residual disease detection, associated with adverse outcome. Thus, we present an orthogonal framework for cfDNA cancer monitoring via genome-wide mutational integration, enabling ultra-sensitive detection, overcoming the limitation of cfDNA abundance and empowering treatment optimization in low-disease-burden oncology care.

A chromosomal connectome for psychiatric and metabolic risk variants in adult dopaminergic neurons.

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.

Dynamic Lamin B1-Gene Association During Oligodendrocyte Progenitor Differentiation.

Differentiation of oligodendrocytes (OL) from progenitor cells (OPC) is the result of a unique program of gene expression, which is further regulated by the formation of topological domains of association with the nuclear lamina. In this study, we show that cultured OPC were characterized by progressively declining levels of endogenous Lamin B1 (LMNB1) during differentiation into OL. We then identify the genes dynamically associated to the nuclear lamina component LMNB1 during this transition, using a well established technique called DamID, which is based on the ability of a bacterially-derived deoxyadenosine methylase (Dam), to modify genomic regions in close proximity. We expressed a fusion protein containing Dam and LMNB1 in OPC (OPCLMNB1-Dam) and either kept them proliferating or differentiated them into OL (OLLMNB1-Dam) and identified genes that were dynamically associated to LMNB1 with differentiation. Importantly, we identified Lss, the gene encoding for lanosterol synthase, a key enzyme in cholesterol synthesis, as associated to the nuclear lamina in OLLMNB1-Dam. This finding could at least in part explain the lipid dysregulation previously reported for mouse models of ADLD characterized by persistent LMNB1 expression in oligodendrocytes.

Spatial genome exploration in the context of cognitive and neurological disease.

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 monogenic neurodevelopmental and neurodegenerative disease. Furthermore, integrative cell-type-specific Hi-C and transcriptomic analysis uncovered an expanded genomic risk space for sequences conferring liability for schizophrenia and other cognitive disease. We predict that spatial genome exploration will deliver radically new insights into the brain nucleome in health and disease.

The histone chaperone CAF-1 cooperates with the DNA methyltransferases to maintain Cd4 silencing in cytotoxic T cells.

The transcriptional repression of alternative lineage genes is critical for cell fate commitment. Mechanisms by which locus-specific gene silencing is initiated and heritably maintained during cell division are not clearly understood. To study the maintenance of silent gene states, we investigated how the Cd4 gene is stably repressed in CD8+ T cells. Through CRISPR and shRNA screening, we identified the histone chaperone CAF-1 as a critical component for Cd4 repression. We found that the large subunit of CAF-1, Chaf1a, requires the N-terminal KER domain to associate with the histone deacetylases HDAC1/2 and the histone demethylase LSD1, enzymes that also participate in Cd4 silencing. When CAF-1 was lacking, Cd4 derepression was markedly enhanced in the absence of the de novo DNA methyltransferase Dnmt3a but not the maintenance DNA methyltransferase Dnmt1. In contrast to Dnmt1, Dnmt3a deficiency did not significantly alter levels of DNA methylation at the Cd4 locus. Instead, Dnmt3a deficiency sensitized CD8+ T cells to Cd4 derepression mediated by compromised functions of histone-modifying factors, including the enzymes associated with CAF-1. Thus, we propose that the heritable silencing of the Cd4 gene in CD8+ T cells exploits cooperative functions among the DNA methyltransferases, CAF-1, and histone-modifying enzymes.

Neuron-specific signatures in the chromosomal connectome associated with schizophrenia risk.

To explore the developmental reorganization of the three-dimensional genome of the brain in the context of neuropsychiatric disease, we monitored chromosomal conformations in differentiating neural progenitor cells. Neuronal and glial differentiation was associated with widespread developmental remodeling of the chromosomal contact map and included interactions anchored in common variant sequences that confer heritable risk for schizophrenia. We describe cell type-specific chromosomal connectomes composed of schizophrenia risk variants and their distal targets, which altogether show enrichment for genes that regulate neuronal connectivity and chromatin remodeling, and evidence for coordinated transcriptional regulation and proteomic interaction of the participating genes. Developmentally regulated chromosomal conformation changes at schizophrenia-relevant sequences disproportionally occurred in neurons, highlighting the existence of cell type-specific disease risk vulnerabilities in spatial genome organization.

Whole-genome bisulfite sequencing with improved accuracy and cost.

DNA methylation patterns in the genome both reflect and help to mediate transcriptional regulatory processes. The digital nature of DNA methylation, present or absent on each allele, makes this assay capable of quantifying events in subpopulations of cells, whereas genome-wide chromatin studies lack the same quantitative capacity. Testing DNA methylation throughout the genome is possible using whole-genome bisulfite sequencing (WGBS), but the high costs associated with the assay have made it impractical for studies involving more than limited numbers of samples. We have optimized a new transposase-based library preparation assay for the Illumina HiSeq X platform suitable for limited amounts of DNA and providing a major cost reduction for WGBS. By incorporating methylated cytosines during fragment end repair, we reveal an end-repair artifact affecting 1%-2% of reads that we can remove analytically. We show that the use of a high (G + C) content spike-in performs better than PhiX in terms of bisulfite sequencing quality. As expected, the loci with transposase-accessible chromatin are DNA hypomethylated and enriched in flanking regions by post-translational modifications of histones usually associated with positive effects on gene expression. Using these transposase-accessible loci to represent the cis-regulatory loci in the genome, we compared the representation of these loci between WGBS and other genome-wide DNA methylation assays, showing WGBS to outperform substantially all of the alternatives. We conclude that it is now technologically and financially feasible to perform WGBS in larger numbers of samples with greater accuracy than previously possible.

ATRX is a regulator of therapy induced senescence in human cells.

Senescence is a state of stable cell cycle exit with important implications for development and disease. Here, we demonstrate that the chromatin remodeling enzyme ATRX is required for therapy-induced senescence. ATRX accumulates in nuclear foci and is required for therapy-induced senescence in multiple types of transformed cells exposed to either DNA damaging agents or CDK4 inhibitors. Mobilization into foci depends on the ability of ATRX to interact with H3K9me3 histone and HP1. Foci form soon after cells exit the cell cycle, before other hallmarks of senescence appear. Eliminating ATRX in senescent cells destabilizes the senescence-associated heterochromatic foci. Additionally, ATRX binds to and suppresses expression from the HRAS locus; repression of HRAS is sufficient to promote the transition of quiescent cells into senescence and preventing repression blocks progression into senescence. Thus ATRX is a critical regulator of therapy-induced senescence and acts in multiple ways to drive cells into this state.Therapy induced senescence (TIS) is a growth suppressive program activated by cytostatic agents in some cancer cells. Here the authors show that the chromatin remodeling enzyme ATRX is a regulator of TIS and drives cells into this state via multiple mechanisms.

An NF-κB Transcription-Factor-Dependent Lineage-Specific Transcriptional Program Promotes Regulatory T Cell Identity and Function.

Both conventional T (Tconv) cells and regulatory T (Treg) cells are activated through ligation of the T cell receptor (TCR) complex, leading to the induction of the transcription factor NF-κB. In Tconv cells, NF-κB regulates expression of genes essential for T cell activation, proliferation, and function. However the role of NF-κB in Treg function remains unclear. We conditionally deleted canonical NF-κB members p65 and c-Rel in developing and mature Treg cells and found they have unique but partially redundant roles. c-Rel was critical for thymic Treg development while p65 was essential for mature Treg identity and maintenance of immune tolerance. Transcriptome and NF-κB p65 binding analyses demonstrated a lineage specific, NF-κB-dependent transcriptional program, enabled by enhanced chromatin accessibility. These dual roles of canonical NF-κB in Tconv and Treg cells highlight the functional plasticity of the NF-κB signaling pathway and underscores the need for more selective strategies to therapeutically target NF-κB.

The methyltransferase SETDB1 regulates a large neuron-specific topological chromatin domain.

We report locus-specific disintegration of megabase-scale chromosomal conformations in brain after neuronal ablation of Setdb1 (also known as Kmt1e; encodes a histone H3 lysine 9 methyltransferase), including a large topologically associated 1.2-Mb domain conserved in humans and mice that encompasses >70 genes at the clustered protocadherin locus (hereafter referred to as cPcdh). The cPcdh topologically associated domain (TADcPcdh) in neurons from mutant mice showed abnormal accumulation of the transcriptional regulator and three-dimensional (3D) genome organizer CTCF at cryptic binding sites, in conjunction with DNA cytosine hypomethylation, histone hyperacetylation and upregulated expression. Genes encoding stochastically expressed protocadherins were transcribed by increased numbers of cortical neurons, indicating relaxation of single-cell constraint. SETDB1-dependent loop formations bypassed 0.2-1 Mb of linear genome and radiated from the TADcPcdh fringes toward cis-regulatory sequences within the cPcdh locus, counterbalanced shorter-range facilitative promoter-enhancer contacts and carried loop-bound polymorphisms that were associated with genetic risk for schizophrenia. We show that the SETDB1 repressor complex, which involves multiple KRAB zinc finger proteins, shields neuronal genomes from excess CTCF binding and is critically required for structural maintenance of TADcPcdh.

Graded Foxo1 activity in Treg cells differentiates tumour immunity from spontaneous autoimmunity.

Regulatory T (Treg) cells expressing the transcription factor Foxp3 have a pivotal role in maintaining immunological self-tolerance; yet, excessive Treg cell activities suppress anti-tumour immune responses. Compared to the resting Treg (rTreg) cell phenotype in secondary lymphoid organs, Treg cells in non-lymphoid tissues exhibit an activated Treg (aTreg) cell phenotype. However, the function of aTreg cells and whether their generation can be manipulated are largely unexplored. Here we show that the transcription factor Foxo1, previously demonstrated to promote Treg cell suppression of lymphoproliferative diseases, has an unexpected function in inhibiting aTreg-cell-mediated immune tolerance in mice. We find that aTreg cells turned over at a slower rate than rTreg cells, but were not locally maintained in tissues. aTreg cell differentiation was associated with repression of Foxo1-dependent gene transcription, concomitant with reduced Foxo1 expression, cytoplasmic localization and enhanced phosphorylation at the Akt sites. Treg-cell-specific expression of an Akt-insensitive Foxo1 mutant prevented downregulation of lymphoid organ homing molecules, and impeded Treg cell homing to non-lymphoid organs, causing CD8(+) T-cell-mediated autoimmune diseases. Compared to Treg cells from healthy tissues, tumour-infiltrating Treg cells downregulated Foxo1 target genes more substantially. Expression of the Foxo1 mutant at a lower dose was sufficient to deplete tumour-associated Treg cells, activate effector CD8(+) T cells, and inhibit tumour growth without inflicting autoimmunity. Thus, Foxo1 inactivation is essential for the migration of aTreg cells that have a crucial function in suppressing CD8(+) T-cell responses; and the Foxo signalling pathway in Treg cells can be titrated to break tumour immune tolerance preferentially.

Genome Wide Mapping of Foxo1 Binding-sites in Murine T Lymphocytes.

The Forkhead box O (Foxo) family of transcription factors has a critical role in controlling the development, differentiation, and function of T cells. However, the direct target genes of Foxo transcription factors in T cells have not been well characterized. In this study, we focused on mapping the genome wide Foxo1-binding sites in naïve CD4+ T cells, CD8+ T cells, and Foxp3+ regulatory T (Treg) cells. By using chromatin immunoprecipitation coupled with deep sequencing (ChIP-Seq), we identified Foxo1 binding sites that were shared among or specific to the three T cell populations. Here we describe the experiments, quality controls, as well as the deep sequencing data. Part of the data analysis has been published by Ouyang W et al. in Nature 2012[1] and Kim MV et al. in Immunity 2013[2], and the associated data set were uploaded to NCBI Gene Expression Omnibus.