Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

MSL2 ensures biallelic gene expression in mammals.

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling.

Genome rearrangements that occur during evolution impose major challenges on regulatory mechanisms that rely on three-dimensional genome architecture. Here, we developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. We observe extensive genome shuffling between these species with one synteny breakpoint after approximately every six genes. A/B compartments, a set of large gene-dense topologically associating domains (TADs), and spatial contacts between high-affinity sites (HAS) located on the X chromosome are maintained over 40 million years, indicating architectural conservation at various hierarchies. Evolutionary conserved genes cluster in the vicinity of HAS, while HAS locations appear evolutionarily flexible, thus uncoupling functional requirement of dosage compensation from individual positions on the linear X chromosome. Therefore, 3D architecture is preserved even in scenarios of thousands of rearrangements highlighting its relevance for essential processes such as dosage compensation of the X chromosome.

Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

Mechanism suppressing H3K9 trimethylation in pluripotent stem cells and its demise by polyQ-expanded huntingtin mutations.

Pluripotent stem cells are invaluable resources to study development and disease, holding a great promise for regenerative medicine. Here we use human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) from patients with Huntington's disease (HD-iPSCs) to shed light into the normal function of huntingtin (HTT) and its demise in disease. We find that HTT binds ATF7IP, a regulator of the histone H3 methyltransferase SETDB1. HTT inhibits the interaction of the ATF7IP-SETDB1 complex with other heterochromatin regulators and transcriptional repressors, maintaining low levels of H3K9 trimethylation (H3K9me3) in hESCs. Loss of HTT promotes global increased H3K9me3 levels and enrichment of H3K9me3 marks at distinct genes, including transcriptional regulators of neuronal differentiation. Although these genes are normally expressed at low amounts in hESCs, HTT knockdown (KD) reduces their induction during neural differentiation. Notably, mutant expanded polyglutamine repeats in HTT diminish its interaction with ATF7IP-SETDB1 complex and trigger H3K9me3 in HD-iPSCs. Conversely, KD of ATF7IP in HD-iPSCs reduces H3K9me3 alterations and ameliorates gene expression changes in their neural counterparts. Taken together, our results indicate ATF7IP as a potential target to correct aberrant H3K9me3 levels induced by mutant HTT.

Standardizing chromatin research: a simple and universal method for ChIP-seq.

Chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) is a key technique in chromatin research. Although heavily applied, existing ChIP-seq protocols are often highly fine-tuned workflows, optimized for specific experimental requirements. Especially the initial steps of ChIP-seq, particularly chromatin shearing, are deemed to be exceedingly cell-type-specific, thus impeding any protocol standardization efforts. Here we demonstrate that harmonization of ChIP-seq workflows across cell types and conditions is possible when obtaining chromatin from properly isolated nuclei. We established an ultrasound-based nuclei extraction method (NEXSON: Nuclei EXtraction by SONication) that is highly effective across various organisms, cell types and cell numbers. The described method has the potential to replace complex cell-type-specific, but largely ineffective, nuclei isolation protocols. By including NEXSON in ChIP-seq workflows, we completely eliminate the need for extensive optimization and sample-dependent adjustments. Apart from this significant simplification, our approach also provides the basis for a fully standardized ChIP-seq and yields highly reproducible transcription factor and histone modifications maps for a wide range of different cell types. Even small cell numbers (∼10,000 cells per ChIP) can be easily processed without application of modified chromatin or library preparation protocols.

Hydrogen peroxide-mediated induction of SOD1 gene transcription is independent from Nrf2 in a cellular model of neurodegeneration.

BACKGROUND: It is still unclear whether oxidative stress (OS) is a disease consequence or is directly involved in the etiology of neurodegenerative disorders (NDs) onset and/or progression; however, many of these conditions are associated with increased levels of oxidation markers and damaged cell components. Previously we demonstrated the accumulation of reactive oxygen species (ROS) and increased SOD1 gene expression in H2O2-treated SH-SY5Y cells, recapitulating pathological features of Amyotrophic Lateral Sclerosis (ALS). Since we observed a post-transcriptional regulation of SOD1 gene in this cellular model, we investigated the transcriptional regulation of SOD1 mRNA under oxidative stress (OS). RESULTS: In response to H2O2 treatment, PolII increased its association to SOD1 promoter. Electrophoretic mobility shift assays and mass spectrometry analyses on SOD1 promoter highlighted the formation of a transcriptional complex bound to the ARE sequences. Western Blotting experiments showed that in our in vitro model, H2O2 exposure increases Nrf2 expression in the nuclear fraction while immunoprecipitation confirmed its phosphorylation and release from Keap1 inhibition. However, H2O2 treatment did not modify Nrf2 binding on SOD1 promoter, which seems to be regulated by different transcription factors (TFs). CONCLUSIONS: Although our data suggest that SOD1 is transcriptionally regulated in response to OS, Nrf2 does not appear to associate with SOD1 promoter in this cellular model of neurodegeneration. Our results open new perspectives in the comprehension of two key antioxidant pathways involved in neurodegenerative disorders.

XPD mutations in trichothiodystrophy hamper collagen VI expression and reveal a role of TFIIH in transcription derepression.

Mutations in the XPD subunit of the transcription/DNA repair factor (TFIIH) give rise to trichothiodystrophy (TTD), a rare hereditary multisystem disorder with skin abnormalities. Here, we show that TTD primary dermal fibroblasts contain low amounts of collagen type VI alpha1 subunit (COL6A1), a fundamental component of soft connective tissues. We demonstrate that COL6A1 expression is downregulated by the sterol regulatory element-binding protein-1 (SREBP-1) whose removal from the promoter is a key step in COL6A1 transcription upregulation in response to cell confluence. We provide evidence for TFIIH being involved in transcription derepression, thus highlighting a new function of TFIIH in gene expression regulation. The lack of COL6A1 upregulation in TTD is caused by the inability of the mutated TFIIH complexes to remove SREBP-1 from COL6A1 promoter and to sustain the subsequent high rate of COL6A1 transcription. This defect might account for the pathologic features that TTD shares with hereditary disorders because of mutations in COL6A genes.

Identification of regulatory links between transcription and RNA processing with long-read sequencing.

We present a detailed protocol for sequencing full-length mRNA isoforms using the Oxford nanopore long-read sequencing technology. We describe steps for poly(A) RNA isolation, library preparation, and cDNA size selection. We then detail procedures for sequencing and processing and a computational framework to identify exon couplings and assign mRNA 5' ends and 3' ends to each other. Our approach enables the identification of links between transcription initiation and co-transcriptional RNA processing events. For complete details on the use and execution of this protocol, please refer to Alfonso-Gonzalez et al.1.

Isolation and Processing of Murine White Adipocytes for Transcriptome and Epigenome Analyses.

Obesity is a complex disease influenced by genetics, epigenetics, the environment, and their interactions. Mature adipocytes represent the major cell type in white adipose tissue. Understanding how adipocytes function and respond to (epi)genetic and environmental signals is essential for identifying the cause(s) of obesity. RNA and chromatin have previously been isolated from adipocytes using enzymatic digestion. In addition, protocols have been developed for nuclear isolation, where purification is achieved by fluorescence-activated cell sorting (FACS) of adipocyte-specific transgenic reporters. One of the greatest challenges to achieving high yield and quality during such protocols is the substantial amount of lipid contained in adipose tissue. The present protocol describes an optimized procedure for isolating mature adipocytes that leverages heptane to separate lipids from the targets of interest (RNA/chromatin). The resulting RNA has high integrity and generates high-quality RNA-seq results. Likewise, the procedure improves nuclei yield rate and generates reproducible ChIP-seq results across samples. Therefore, the current study provides a reliable and universal murine adipocyte isolation protocol suitable for whole-genome transcriptome and epigenome studies.

The glucose transporter GLUT3 controls T helper 17 cell responses through glycolytic-epigenetic reprogramming.

Metabolic reprogramming is a hallmark of activated T cells. The switch from oxidative phosphorylation to aerobic glycolysis provides energy and intermediary metabolites for the biosynthesis of macromolecules to support clonal expansion and effector function. Here, we show that glycolytic reprogramming additionally controls inflammatory gene expression via epigenetic remodeling. We found that the glucose transporter GLUT3 is essential for the effector functions of Th17 cells in models of autoimmune colitis and encephalomyelitis. At the molecular level, we show that GLUT3-dependent glucose uptake controls a metabolic-transcriptional circuit that regulates the pathogenicity of Th17 cells. Metabolomic, epigenetic, and transcriptomic analyses linked GLUT3 to mitochondrial glucose oxidation and ACLY-dependent acetyl-CoA generation as a rate-limiting step in the epigenetic regulation of inflammatory gene expression. Our findings are also important from a translational perspective because inhibiting GLUT3-dependent acetyl-CoA generation is a promising metabolic checkpoint to mitigate Th17-cell-mediated inflammatory diseases.

DOT1L-mediated murine neuronal differentiation associates with H3K79me2 accumulation and preserves SOX2-enhancer accessibility.

During neuronal differentiation, the transcriptional profile and the epigenetic context of neural committed cells is subject to significant rearrangements, but a systematic quantification of global histone modification changes is still missing. Here, we show that H3K79me2 increases and H3K27ac decreases globally during in-vitro neuronal differentiation of murine embryonic stem cells. DOT1L mediates all three degrees of methylation of H3K79 and its enzymatic activity is critical to modulate cellular differentiation and reprogramming. In this context, we find that inhibition of DOT1L in neural progenitor cells biases the transcriptional state towards neuronal differentiation, resulting in transcriptional upregulation of genes marked with H3K27me3 on the promoter region. We further show that DOT1L inhibition affects accessibility of SOX2-bound enhancers and impairs SOX2 binding in neural progenitors. Our work provides evidence that DOT1L activity gates differentiation of progenitors by allowing SOX2-dependent transcription of stemness programs.

CRUP: a comprehensive framework to predict condition-specific regulatory units.

We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

A comprehensive analysis of 195 DNA methylomes reveals shared and cell-specific features of partially methylated domains.

BACKGROUND: Partially methylated domains are extended regions in the genome exhibiting a reduced average DNA methylation level. They cover gene-poor and transcriptionally inactive regions and tend to be heterochromatic. We present a comprehensive comparative analysis of partially methylated domains in human and mouse cells, to identify structural and functional features associated with them. RESULTS: Partially methylated domains are present in up to 75% of the genome in human and mouse cells irrespective of their tissue or cell origin. Each cell type has a distinct set of partially methylated domains, and genes expressed in such domains show a strong cell type effect. The methylation level varies between cell types with a more pronounced effect in differentiating and replicating cells. The lowest level of methylation is observed in highly proliferating and immortal cancer cell lines. A decrease of DNA methylation within partially methylated domains tends to be linked to an increase in heterochromatic histone marks and a decrease of gene expression. Characteristic combinations of heterochromatic signatures in partially methylated domains are linked to domains of early and middle S-phase and late S-G2 phases of DNA replication. CONCLUSIONS: Partially methylated domains are prominent signatures of long-range epigenomic organization. Integrative analysis identifies them as important general, lineage- and cell type-specific topological features. Changes in partially methylated domains are hallmarks of cell differentiation, with decreased methylation levels and increased heterochromatic marks being linked to enhanced cell proliferation. In combination with broad histone marks, partially methylated domains demarcate distinct domains of late DNA replication.

RELACS nuclei barcoding enables high-throughput ChIP-seq.

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) is an invaluable tool for mapping chromatin-associated proteins. Current barcoding strategies aim to improve assay throughput and scalability but intense sample handling and lack of standardization over cell types, cell numbers and epitopes hinder wide-spread use in the field. Here, we present a barcoding method to enable high-throughput ChIP-seq using common molecular biology techniques. The method, called RELACS (restriction enzyme-based labeling of chromatin in situ) relies on standardized nuclei extraction from any source and employs chromatin cutting and barcoding within intact nuclei. Barcoded nuclei are pooled and processed within the same ChIP reaction, for maximal comparability and workload reduction. The innovative barcoding concept is particularly user-friendly and suitable for implementation to standardized large-scale clinical studies and scarce samples. Aiming to maximize universality and scalability, RELACS can generate ChIP-seq libraries for transcription factors and histone modifications from hundreds of samples within three days.

High-resolution TADs reveal DNA sequences underlying genome organization in flies.

Despite an abundance of new studies about topologically associating domains (TADs), the role of genetic information in TAD formation is still not fully understood. Here we use our software, HiCExplorer (hicexplorer.readthedocs.io) to annotate >2800 high-resolution (570 bp) TAD boundaries in Drosophila melanogaster. We identify eight DNA motifs enriched at boundaries, including a motif bound by the M1BP protein, and two new boundary motifs. In contrast to mammals, the CTCF motif is only enriched on a small fraction of boundaries flanking inactive chromatin while most active boundaries contain the motifs bound by the M1BP or Beaf-32 proteins. We demonstrate that boundaries can be accurately predicted using only the motif sequences at open chromatin sites. We propose that DNA sequence guides the genome architecture by allocation of boundary proteins in the genome. Finally, we present an interactive online database to access and explore the spatial organization of fly, mouse and human genomes, available at http://chorogenome.ie-freiburg.mpg.de .

The Polycomb-Dependent Epigenome Controls ß Cell Dysfunction, Dedifferentiation, and Diabetes.

To date, it remains largely unclear to what extent chromatin machinery contributes to the susceptibility and progression of complex diseases. Here, we combine deep epigenome mapping with single-cell transcriptomics to mine for evidence of chromatin dysregulation in type 2 diabetes. We find two chromatin-state signatures that track ß cell dysfunction in mice and humans: ectopic activation of bivalent Polycomb-silenced domains and loss of expression at an epigenomically unique class of lineage-defining genes. ß cell-specific Polycomb (Eed/PRC2) loss of function in mice triggers diabetes-mimicking transcriptional signatures and highly penetrant, hyperglycemia-independent dedifferentiation, indicating that PRC2 dysregulation contributes to disease. The work provides novel resources for exploring ß cell transcriptional regulation and identifies PRC2 as necessary for long-term maintenance of ß cell identity. Importantly, the data suggest a two-hit (chromatin and hyperglycemia) model for loss of ß cell identity in diabetes.

The histone code reader Spin1 controls skeletal muscle development.

While several studies correlated increased expression of the histone code reader Spin1 with tumor formation or growth, little is known about physiological functions of the protein. We generated Spin1M5 mice with ablation of Spin1 in myoblast precursors using the Myf5-Cre deleter strain. Most Spin1M5 mice die shortly after birth displaying severe sarcomere disorganization and necrosis. Surviving Spin1M5 mice are growth-retarded and exhibit the most prominent defects in soleus, tibialis anterior, and diaphragm muscle. Transcriptome analyses of limb muscle at embryonic day (E) 15.5, E16.5, and at three weeks of age provided evidence for aberrant fetal myogenesis and identified deregulated skeletal muscle (SkM) functional networks. Determination of genome-wide chromatin occupancy in primary myoblast revealed direct Spin1 target genes and suggested that deregulated basic helix-loop-helix transcription factor networks account for developmental defects in Spin1M5 fetuses. Furthermore, correlating histological and transcriptome analyses, we show that aberrant expression of titin-associated proteins, abnormal glycogen metabolism, and neuromuscular junction defects contribute to SkM pathology in Spin1M5 mice. Together, we describe the first example of a histone code reader controlling SkM development in mice, which hints at Spin1 as a potential player in human SkM disease.

Inactivation of Lsd1 triggers senescence in trophoblast stem cells by induction of Sirt4.

Coordination of energy metabolism is essential for homeostasis of stem cells, whereas an imbalance in energy homeostasis causes disease and accelerated aging. Here we show that deletion or enzymatic inactivation of lysine-specific demethylase 1 (Lsd1) triggers senescence in trophoblast stem cells (TSCs). Genome-wide transcriptional profiling of TSCs following Lsd1 inhibition shows gene set enrichment of aging and metabolic pathways. Consistently, global metabolomic and phenotypic analyses disclose an unbalanced redox status, decreased glutamine anaplerosis and mitochondrial function. Loss of homeostasis is caused by increased expression of sirtuin 4 (Sirt4), a Lsd1-repressed direct target gene. Accordingly, Sirt4 overexpression in wild-type TSCs recapitulates the senescence phenotype initiated by Lsd1 deletion or inhibition. Inversely, absence of Lsd1 enzymatic activity concomitant with knockdown of Sirt4 reestablishes normal glutamine anaplerosis, redox balance and mitochondrial function. In conclusion, by repression of Sirt4, Lsd1 directs the epigenetic control of TSC immortality via maintenance of metabolic flexibility.

Lsd1 Ablation Triggers Metabolic Reprogramming of Brown Adipose Tissue.

Previous work indicated that lysine-specific demethylase 1 (Lsd1) can positively regulate the oxidative and thermogenic capacities of white and beige adipocytes. Here we investigate the role of Lsd1 in brown adipose tissue (BAT) and find that BAT-selective Lsd1 ablation induces a shift from oxidative to glycolytic metabolism. This shift is associated with downregulation of BAT-specific and upregulation of white adipose tissue (WAT)-selective gene expression. This results in the accumulation of di- and triacylglycerides and culminates in a profound whitening of BAT in aged Lsd1-deficient mice. Further studies show that Lsd1 maintains BAT properties via a dual role. It activates BAT-selective gene expression in concert with the transcription factor Nrf1 and represses WAT-selective genes through recruitment of the CoREST complex. In conclusion, our data uncover Lsd1 as a key regulator of gene expression and metabolic function in BAT.