RNA damage compartmentalization by DHX9 stress granules.

Biomolecules incur damage during stress conditions, and damage partitioning represents a vital survival strategy for cells. Here, we identified a distinct stress granule (SG), marked by dsRNA helicase DHX9, which compartmentalizes ultraviolet (UV)-induced RNA, but not DNA, damage. Our FANCI technology revealed that DHX9 SGs are enriched in damaged intron RNA, in contrast to classical SGs that are composed of mature mRNA. UV exposure causes RNA crosslinking damage, impedes intron splicing and decay, and triggers DHX9 SGs within daughter cells. DHX9 SGs promote cell survival and induce dsRNA-related immune response and translation shutdown, differentiating them from classical SGs that assemble downstream of translation arrest. DHX9 modulates dsRNA abundance in the DHX9 SGs and promotes cell viability. Autophagy receptor p62 is activated and important for DHX9 SG disassembly. Our findings establish non-canonical DHX9 SGs as a dedicated non-membrane-bound cytoplasmic compartment that safeguards daughter cells from parental RNA damage.

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.

Multimodal epigenetic changes and altered NEUROD1 chromatin binding in the mouse hippocampus underlie FOXG1 syndrome.

Forkhead box G1 (FOXG1) has important functions in neuronal differentiation and balances excitatory/inhibitory network activity. Thus far, molecular processes underlying FOXG1 function are largely unexplored. Here, we present a multiomics data set exploring how FOXG1 impacts neuronal maturation at the chromatin level in the mouse hippocampus. At a genome-wide level, FOXG1 i) both represses and activates transcription, ii) binds mainly to enhancer regions, iii) reconfigures the epigenetic landscape through bidirectional alteration of H3K27ac, H3K4me3, and chromatin accessibility, and iv) operates synergistically with NEUROD1. Interestingly, we could not detect a clear hierarchy of FOXG1 and NEUROD1, but instead, provide the evidence that they act in a highly cooperative manner to control neuronal maturation. Genes affected by the chromatin alterations impact synaptogenesis and axonogenesis. Inhibition of histone deacetylases partially rescues transcriptional alterations upon FOXG1 reduction. This integrated multiomics view of changes upon FOXG1 reduction reveals an unprecedented multimodality of FOXG1 functions converging on neuronal maturation. It fuels therapeutic options based on epigenetic drugs to alleviate, at least in part, neuronal dysfunction.

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.

Nuclear Architecture Organized by Rif1 Underpins the Replication-Timing Program.

DNA replication is temporally and spatially organized in all eukaryotes, yet the molecular control and biological function of the replication-timing program are unclear. Rif1 is required for normal genome-wide regulation of replication timing, but its molecular function is poorly understood. Here we show that in mouse embryonic stem cells, Rif1 coats late-replicating domains and, with Lamin B1, identifies most of the late-replicating genome. Rif1 is an essential determinant of replication timing of non-Lamin B1-bound late domains. We further demonstrate that Rif1 defines and restricts the interactions between replication-timing domains during the G1 phase, thereby revealing a function of Rif1 as organizer of nuclear architecture. Rif1 loss affects both number and replication-timing specificity of the interactions between replication-timing domains. In addition, during the S phase, Rif1 ensures that replication of interacting domains is temporally coordinated. In summary, our study identifies Rif1 as the molecular link between nuclear architecture and replication-timing establishment in mammals.

DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome.

Transposable elements are viewed as 'selfish genetic elements', yet they contribute to gene regulation and genome evolution in diverse ways. More than half of the human genome consists of transposable elements. Alu elements belong to the short interspersed nuclear element (SINE) family of repetitive elements, and with over 1 million insertions they make up more than 10% of the human genome. Despite their abundance and the potential evolutionary advantages they confer, Alu elements can be mutagenic to the host as they can act as splice acceptors, inhibit translation of mRNAs and cause genomic instability. Alu elements are the main targets of the RNA-editing enzyme ADAR and the formation of Alu exons is suppressed by the nuclear ribonucleoprotein HNRNPC, but the broad effect of massive secondary structures formed by inverted-repeat Alu elements on RNA processing in the nucleus remains unknown. Here we show that DHX9, an abundant nuclear RNA helicase, binds specifically to inverted-repeat Alu elements that are transcribed as parts of genes. Loss of DHX9 leads to an increase in the number of circular-RNA-producing genes and amount of circular RNAs, translational repression of reporters containing inverted-repeat Alu elements, and transcriptional rewiring (the creation of mostly nonsensical novel connections between exons) of susceptible loci. Biochemical purifications of DHX9 identify the interferon-inducible isoform of ADAR (p150), but not the constitutively expressed ADAR isoform (p110), as an RNA-independent interaction partner. Co-depletion of ADAR and DHX9 augments the double-stranded RNA accumulation defects, leading to increased circular RNA production, revealing a functional link between these two enzymes. Our work uncovers an evolutionarily conserved function of DHX9. We propose that it acts as a nuclear RNA resolvase that neutralizes the immediate threat posed by transposon insertions and allows these elements to evolve as tools for the post-transcriptional regulation of gene expression.

High-Affinity Sites Form an Interaction Network to Facilitate Spreading of the MSL Complex across the X Chromosome in Drosophila.

Dosage compensation mechanisms provide a paradigm to study the contribution of chromosomal conformation toward targeting and spreading of epigenetic regulators over a specific chromosome. By using Hi-C and 4C analyses, we show that high-affinity sites (HAS), landing platforms of the male-specific lethal (MSL) complex, are enriched around topologically associating domain (TAD) boundaries on the X chromosome and harbor more long-range contacts in a sex-independent manner. Ectopically expressed roX1 and roX2 RNAs target HAS on the X chromosome in trans and, via spatial proximity, induce spreading of the MSL complex in cis, leading to increased expression of neighboring autosomal genes. We show that the MSL complex regulates nucleosome positioning at HAS, therefore acting locally rather than influencing the overall chromosomal architecture. We propose that the sex-independent, three-dimensional conformation of the X chromosome poises it for exploitation by the MSL complex, thereby facilitating spreading in males.

pyGenomeTracks: reproducible plots for multivariate genomic datasets.

MOTIVATION: Generating publication ready plots to display multiple genomic tracks can pose a serious challenge. Making desirable and accurate figures requires considerable effort. This is usually done by hand or using a vector graphic software. RESULTS: pyGenomeTracks (PGT) is a modular plotting tool that easily combines multiple tracks. It enables a reproducible and standardized generation of highly customizable and publication ready images. AVAILABILITY AND IMPLEMENTATION: PGT is available through a graphical interface on https://usegalaxy.eu and through the command line. It is provided on conda via the bioconda channel, on pip and it is openly developed on github: https://github.com/deeptools/pyGenomeTracks. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Suv39h-dependent H3K9me3 marks intact retrotransposons and silences LINE elements in mouse embryonic stem cells.

Heterochromatin is required to restrict aberrant expression of retrotransposons, but it remains poorly defined due to the underlying repeat-rich sequences. We dissected Suv39h-dependent histone H3 lysine 9 trimethylation (H3K9me3) by genome-wide ChIP sequencing in mouse embryonic stem cells (ESCs). Refined bioinformatic analyses of repeat subfamilies indicated selective accumulation of Suv39h-dependent H3K9me3 at interspersed repetitive elements that cover ∼5% of the ESC epigenome. The majority of the ∼8,150 intact long interspersed nuclear elements (LINEs) and endogenous retroviruses (ERVs), but only a minor fraction of the >1.8 million degenerate and truncated LINEs/ERVs, are enriched for Suv39h-dependent H3K9me3. Transcriptional repression of intact LINEs and ERVs is differentially regulated by Suv39h and other chromatin modifiers in ESCs but governed by DNA methylation in committed cells. These data provide a function for Suv39h-dependent H3K9me3 chromatin to specifically repress intact LINE elements in the ESC epigenome.

Galaxy HiCExplorer 3: a web server for reproducible Hi-C, capture Hi-C and single-cell Hi-C data analysis, quality control and visualization.

The Galaxy HiCExplorer provides a web service at https://hicexplorer.usegalaxy.eu. It enables the integrative analysis of chromosome conformation by providing tools and computational resources to pre-process, analyse and visualize Hi-C, Capture Hi-C (cHi-C) and single-cell Hi-C (scHi-C) data. Since the last publication, Galaxy HiCExplorer has been expanded considerably with new tools to facilitate the analysis of cHi-C and to provide an in-depth analysis of Hi-C data. Moreover, it supports the analysis of scHi-C data by offering a broad range of tools. With the help of the standard graphical user interface of Galaxy, presented workflows, extensive documentation and tutorials, novices as well as Hi-C experts are supported in their Hi-C data analysis with Galaxy HiCExplorer.

DOT1L promotes progenitor proliferation and primes neuronal layer identity in the developing cerebral cortex.

Cortical development is controlled by transcriptional programs, which are orchestrated by transcription factors. Yet, stable inheritance of spatio-temporal activity of factors influencing cell fate and localization in different layers is only partly understood. Here we find that deletion of Dot1l in the murine telencephalon leads to cortical layering defects, indicating DOT1L activity and chromatin methylation at H3K79 impact on the cell cycle, and influence transcriptional programs conferring upper layer identity in early progenitors. Specifically, DOT1L prevents premature differentiation by increasing expression of genes that regulate asymmetric cell division (Vangl2, Cenpj). Loss of DOT1L results in reduced numbers of progenitors expressing genes including SoxB1 gene family members. Loss of DOT1L also leads to altered cortical distribution of deep layer neurons that express either TBR1, CTIP2 or SOX5, and less activation of transcriptional programs that are characteristic for upper layer neurons (Satb2, Pou3f3, Cux2, SoxC family members). Data from three different mouse models suggest that DOT1L balances transcriptional programs necessary for proper neuronal composition and distribution in the six cortical layers. Furthermore, because loss of DOT1L in the pre-neurogenic phase of development impairs specifically generation of SATB2-expressing upper layer neurons, our data suggest that DOT1L primes upper layer identity in cortical progenitors.

snakePipes: facilitating flexible, scalable and integrative epigenomic analysis.

SUMMARY: Due to the rapidly increasing scale and diversity of epigenomic data, modular and scalable analysis workflows are of wide interest. Here we present snakePipes, a workflow package for processing and downstream analysis of data from common epigenomic assays: ChIP-seq, RNA-seq, Bisulfite-seq, ATAC-seq, Hi-C and single-cell RNA-seq. snakePipes enables users to assemble variants of each workflow and to easily install and upgrade the underlying tools, via its simple command-line wrappers and yaml files. AVAILABILITY AND IMPLEMENTATION: snakePipes can be installed via conda: `conda install -c mpi-ie -c bioconda -c conda-forge snakePipes'. Source code (https://github.com/maxplanck-ie/snakepipes) and documentation (https://snakepipes.readthedocs.io/en/latest/) are available online. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Galaxy HiCExplorer: a web server for reproducible Hi-C data analysis, quality control and visualization.

Galaxy HiCExplorer is a web server that facilitates the study of the 3D conformation of chromatin by allowing Hi-C data processing, analysis and visualization. With the Galaxy HiCExplorer web server, users with little bioinformatic background can perform every step of the analysis in one workflow: mapping of the raw sequence data, creation of Hi-C contact matrices, quality assessment, correction of contact matrices and identification of topological associated domains (TADs) and A/B compartments. Users can create publication ready plots of the contact matrix, A/B compartments, and TADs on a selected genomic locus, along with additional information like gene tracks or ChIP-seq signals. Galaxy HiCExplorer is freely usable at: https://hicexplorer.usegalaxy.eu and is available as a Docker container: https://github.com/deeptools/docker-galaxy-hicexplorer.

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.

deepTools2: a next generation web server for deep-sequencing data analysis.

We present an update to our Galaxy-based web server for processing and visualizing deeply sequenced data. Its core tool set, deepTools, allows users to perform complete bioinformatic workflows ranging from quality controls and normalizations of aligned reads to integrative analyses, including clustering and visualization approaches. Since we first described our deepTools Galaxy server in 2014, we have implemented new solutions for many requests from the community and our users. Here, we introduce significant enhancements and new tools to further improve data visualization and interpretation. deepTools continue to be open to all users and freely available as a web service at deeptools.ie-freiburg.mpg.de The new deepTools2 suite can be easily deployed within any Galaxy framework via the toolshed repository, and we also provide source code for command line usage under Linux and Mac OS X. A public and documented API for access to deepTools functionality is also available.

Putative Prostate Cancer Risk SNP in an Androgen Receptor-Binding Site of the Melanophilin Gene Illustrates Enrichment of Risk SNPs in Androgen Receptor Target Sites.

Genome-wide association studies have identified genomic loci, whose single-nucleotide polymorphisms (SNPs) predispose to prostate cancer (PCa). However, the mechanisms of most of these variants are largely unknown. We integrated chromatin-immunoprecipitation-coupled sequencing and microarray expression profiling in TMPRSS2-ERG gene rearrangement positive DUCaP cells with the GWAS PCa risk SNPs catalog to identify disease susceptibility SNPs localized within functional androgen receptor-binding sites (ARBSs). Among the 48 GWAS index risk SNPs and 3,917 linked SNPs, 80 were found located in ARBSs. Of these, rs11891426:T>G in an intron of the melanophilin gene (MLPH) was within a novel putative auxiliary AR-binding motif, which is enriched in the neighborhood of canonical androgen-responsive elements. T→G exchange attenuated the transcriptional activity of the ARBS in an AR reporter gene assay. The expression of MLPH in primary prostate tumors was significantly lower in those with the G compared with the T allele and correlated significantly with AR protein. Higher melanophilin level in prostate tissue of patients with a favorable PCa risk profile points out a tumor-suppressive effect. These results unravel a hidden link between AR and a functional putative PCa risk SNP, whose allele alteration affects androgen regulation of its host gene MLPH.

deepTools: a flexible platform for exploring deep-sequencing data.

We present a Galaxy based web server for processing and visualizing deeply sequenced data. The web server's core functionality consists of a suite of newly developed tools, called deepTools, that enable users with little bioinformatic background to explore the results of their sequencing experiments in a standardized setting. Users can upload pre-processed files with continuous data in standard formats and generate heatmaps and summary plots in a straight-forward, yet highly customizable manner. In addition, we offer several tools for the analysis of files containing aligned reads and enable efficient and reproducible generation of normalized coverage files. As a modular and open-source platform, deepTools can easily be expanded and customized to future demands and developments. The deepTools webserver is freely available at http://deeptools.ie-freiburg.mpg.de and is accompanied by extensive documentation and tutorials aimed at conveying the principles of deep-sequencing data analysis. The web server can be used without registration. deepTools can be installed locally either stand-alone or as part of Galaxy.

Nontransformed, GM-CSF-dependent macrophage lines are a unique model to study tissue macrophage functions.

Macrophages are diverse cell types in the first line of antimicrobial defense. Only a limited number of primary mouse models exist to study their function. Bone marrow-derived, macrophage-CSF-induced cells with a limited life span are the most common source. We report here a simple method yielding self-renewing, nontransformed, GM-CSF/signal transducer and activator of transcription 5-dependent macrophages (Max Planck Institute cells) from mouse fetal liver, which reflect the innate immune characteristics of alveolar macrophages. Max Planck Institute cells are exquisitely sensitive to selected microbial agents, including bacterial LPS, lipopeptide, Mycobacterium tuberculosis, cord factor, and adenovirus and mount highly proinflammatory but no anti-inflammatory IL-10 responses. They show a unique pattern of innate responses not yet observed in other mononuclear phagocytes. This includes differential LPS sensing and an unprecedented regulation of IL-1α production upon LPS exposure, which likely plays a key role in lung inflammation in vivo. In conclusion, Max Planck Institute cells offer an useful tool to study macrophage biology and for biomedical science.

The NSL complex regulates housekeeping genes in Drosophila.

MOF is the major histone H4 lysine 16-specific (H4K16) acetyltransferase in mammals and Drosophila. In flies, it is involved in the regulation of X-chromosomal and autosomal genes as part of the MSL and the NSL complexes, respectively. While the function of the MSL complex as a dosage compensation regulator is fairly well understood, the role of the NSL complex in gene regulation is still poorly characterized. Here we report a comprehensive ChIP-seq analysis of four NSL complex members (NSL1, NSL3, MBD-R2, and MCRS2) throughout the Drosophila melanogaster genome. Strikingly, the majority (85.5%) of NSL-bound genes are constitutively expressed across different cell types. We find that an increased abundance of the histone modifications H4K16ac, H3K4me2, H3K4me3, and H3K9ac in gene promoter regions is characteristic of NSL-targeted genes. Furthermore, we show that these genes have a well-defined nucleosome free region and broad transcription initiation patterns. Finally, by performing ChIP-seq analyses of RNA polymerase II (Pol II) in NSL1- and NSL3-depleted cells, we demonstrate that both NSL proteins are required for efficient recruitment of Pol II to NSL target gene promoters. The observed Pol II reduction coincides with compromised binding of TBP and TFIIB to target promoters, indicating that the NSL complex is required for optimal recruitment of the pre-initiation complex on target genes. Moreover, genes that undergo the most dramatic loss of Pol II upon NSL knockdowns tend to be enriched in DNA Replication-related Element (DRE). Taken together, our findings show that the MOF-containing NSL complex acts as a major regulator of housekeeping genes in flies by modulating initiation of Pol II transcription.

ChroMoS: an integrated web tool for SNP classification, prioritization and functional interpretation.

UNLABELLED: Genome-wide association studies and re-sequencing projects are revealing an increasing number of disease-associated SNPs, a large fraction of which are non-coding. Although they could have relevance for disease susceptibility and progression, the lack of information about regulatory regions impedes the assessment of their functionality. Here we present a web server, ChroMoS (Chromatin Modified SNPs), which combines genetic and epigenetic data with the goal of facilitating SNPs' classification, prioritization and prediction of their functional consequences. ChroMoS uses a large database of SNPs and chromatin states, but allows a user to provide his/her own genetic information. Based on the SNP classification and interactive prioritization, a user can compute the functional impact of multiple SNPs using two prediction tools, one for differential analysis of transcription factor binding (sTRAP) and another for SNPs with potential impact on binding of miRNAs (MicroSNiPer). AVAILABILITY: Web server, ChroMoS, is freely available at http://epicenter.immunbio.mpg.de/services/chromos.