Functional multi-organelle units control inflammatory lipid metabolism of macrophages.

Eukaryotic cells contain several membrane-separated organelles to compartmentalize distinct metabolic reactions. However, it has remained unclear how these organelle systems are coordinated when cells adapt metabolic pathways to support their development, survival or effector functions. Here we present OrgaPlexing, a multi-spectral organelle imaging approach for the comprehensive mapping of six key metabolic organelles and their interactions. We use this analysis on macrophages, immune cells that undergo rapid metabolic switches upon sensing bacterial and inflammatory stimuli. Our results identify lipid droplets (LDs) as primary inflammatory responder organelle, which forms three- and four-way interactions with other organelles. While clusters with endoplasmic reticulum (ER) and mitochondria (mitochondria-ER-LD unit) help supply fatty acids for LD growth, the additional recruitment of peroxisomes (mitochondria-ER-peroxisome-LD unit) supports fatty acid efflux from LDs. Interference with individual components of these units has direct functional consequences for inflammatory lipid mediator synthesis. Together, we show that macrophages form functional multi-organellar units to support metabolic adaptation and provide an experimental strategy to identify organelle-metabolic signalling hubs.

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.

EBF1 primes B-lymphoid enhancers and limits the myeloid bias in murine multipotent progenitors.

Hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) generate all cells of the blood system. Despite their multipotency, MPPs display poorly understood lineage bias. Here, we examine whether lineage-specifying transcription factors, such as the B-lineage determinant EBF1, regulate lineage preference in early progenitors. We detect low-level EBF1 expression in myeloid-biased MPP3 and lymphoid-biased MPP4 cells, coinciding with expression of the myeloid determinant C/EBPα. Hematopoietic deletion of Ebf1 results in enhanced myelopoiesis and reduced HSC repopulation capacity. Ebf1-deficient MPP3 and MPP4 cells exhibit an augmented myeloid differentiation potential and a transcriptome with an enriched C/EBPα signature. Correspondingly, EBF1 binds the Cebpa enhancer, and the deficiency and overexpression of Ebf1 in MPP3 and MPP4 cells lead to an up- and downregulation of Cebpa expression, respectively. In addition, EBF1 primes the chromatin of B-lymphoid enhancers specifically in MPP3 cells. Thus, our study implicates EBF1 in regulating myeloid/lymphoid fate bias in MPPs by constraining C/EBPα-driven myelopoiesis and priming the B-lymphoid fate.

Hyaluronic acid-GPRC5C signalling promotes dormancy in haematopoietic stem cells.

Bone marrow haematopoietic stem cells (HSCs) are vital for lifelong maintenance of healthy haematopoiesis. In inbred mice housed in gnotobiotic facilities, the top of the haematopoietic hierarchy is occupied by dormant HSCs, which reversibly exit quiescence during stress. Whether HSC dormancy exists in humans remains debatable. Here, using single-cell RNA sequencing, we show a continuous landscape of highly purified human bone marrow HSCs displaying varying degrees of dormancy. We identify the orphan receptor GPRC5C, which enriches for dormant human HSCs. GPRC5C is also essential for HSC function, as demonstrated by genetic loss- and gain-of-function analyses. Through structural modelling and biochemical assays, we show that hyaluronic acid, a bone marrow extracellular matrix component, preserves dormancy through GPRC5C. We identify the hyaluronic acid-GPRC5C signalling axis controlling the state of dormancy in mouse and human HSCs.

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.

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.

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.

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.

MAPCap allows high-resolution detection and differential expression analysis of transcription start sites.

The position, shape and number of transcription start sites (TSS) are critical determinants of gene regulation. Most methods developed to detect TSSs and study promoter usage are, however, of limited use in studies that demand quantification of expression changes between two or more groups. In this study, we combine high-resolution detection of transcription start sites and differential expression analysis using a simplified TSS quantification protocol, MAPCap (Multiplexed Affinity Purification of Capped RNA) along with the software icetea . Applying MAPCap on developing Drosophila melanogaster embryos and larvae, we detected stage and sex-specific promoter and enhancer activity and quantify the effect of mutants of maleless (MLE) helicase at X-chromosomal promoters. We observe that MLE mutation leads to a median 1.9 fold drop in expression of X-chromosome promoters and affects the expression of several TSSs with a sexually dimorphic expression on autosomes. Our results provide quantitative insights into promoter activity during dosage compensation.

CAPRI enables comparison of evolutionarily conserved RNA interacting regions.

RNA-protein complexes play essential regulatory roles at nearly all levels of gene expression. Using in vivo crosslinking and RNA capture, we report a comprehensive RNA-protein interactome in a metazoan at four levels of resolution: single amino acids, domains, proteins and multisubunit complexes. We devise CAPRI, a method to map RNA-binding domains (RBDs) by simultaneous identification of RNA interacting crosslinked peptides and peptides adjacent to such crosslinked sites. CAPRI identifies more than 3000 RNA proximal peptides in Drosophila and human proteins with more than 45% of them forming new interaction interfaces. The comparison of orthologous proteins enables the identification of evolutionary conserved RBDs in globular domains and intrinsically disordered regions (IDRs). By comparing the sequences of IDRs through evolution, we classify them based on the type of motif, accumulation of tandem repeats, conservation of amino acid composition and high sequence divergence.

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.

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.

Differential Methylation of H3K79 Reveals DOT1L Target Genes and Function in the Cerebellum In Vivo.

The disruptor of telomeric silencing 1-like (DOT1L) mediates methylation of histone H3 at position lysine 79 (H3K79). Conditional knockout of Dot1l in mouse cerebellar granule cells (Dot1l-cKOAtoh1) led to a smaller external granular layer with fewer precursors of granule neurons. Dot1l-cKOAtoh1 mice had impaired proliferation and differentiation of granular progenitors, which resulted in a smaller cerebellum. Mutant mice showed mild ataxia in motor behavior tests. In contrast, Purkinje cell-specific conditional knockout mice showed no obvious phenotype. Genome-wide transcription analysis of Dot1l-cKOAtoh1 cerebella using microarrays revealed changes in genes that function in cell cycle, cell migration, axon guidance, and metabolism. To identify direct DOT1L target genes, we used genome-wide profiling of H3K79me2 and transcriptional analysis. Analysis of differentially methylated regions (DR) and differentially expressed genes (DE) revealed in total 12 putative DOT1L target genes in Dot1l-cKOAtoh1 affecting signaling (Tnfaip8l3, B3galt5), transcription (Otx1), cell migration and axon guidance (Sema4a, Sema5a, Robo1), cholesterol and lipid metabolism (Lss, Cyp51), cell cycle (Cdkn1a), calcium-dependent cell-adhesion or exocytosis (Pcdh17, Cadps2), and unknown function (Fam174b). Dysregulated expression of these target genes might be implicated in the ataxia phenotype observed in Dot1l-cKOAtoh1.

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.

Integrative analysis of single-cell expression data reveals distinct regulatory states in bidirectional promoters.

BACKGROUND: Bidirectional promoters (BPs) are prevalent in eukaryotic genomes. However, it is poorly understood how the cell integrates different epigenomic information, such as transcription factor (TF) binding and chromatin marks, to drive gene expression at BPs. Single-cell sequencing technologies are revolutionizing the field of genome biology. Therefore, this study focuses on the integration of single-cell RNA-seq data with bulk ChIP-seq and other epigenetics data, for which single-cell technologies are not yet established, in the context of BPs. RESULTS: We performed integrative analyses of novel human single-cell RNA-seq (scRNA-seq) data with bulk ChIP-seq and other epigenetics data. scRNA-seq data revealed distinct transcription states of BPs that were previously not recognized. We find associations between these transcription states to distinct patterns in structural gene features, DNA accessibility, histone modification, DNA methylation and TF binding profiles. CONCLUSIONS: Our results suggest that a complex interplay of all of these elements is required to achieve BP-specific transcriptional output in this specialized promoter configuration. Further, our study implies that novel statistical methods can be developed to deconvolute masked subpopulations of cells measured with different bulk epigenomic assays using scRNA-seq data.

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.

Community-Driven Data Analysis Training for Biology.

The primary problem with the explosion of biomedical datasets is not the data, not computational resources, and not the required storage space, but the general lack of trained and skilled researchers to manipulate and analyze these data. Eliminating this problem requires development of comprehensive educational resources. Here we present a community-driven framework that enables modern, interactive teaching of data analytics in life sciences and facilitates the development of training materials. The key feature of our system is that it is not a static but a continuously improved collection of tutorials. By coupling tutorials with a web-based analysis framework, biomedical researchers can learn by performing computation themselves through a web browser without the need to install software or search for example datasets. Our ultimate goal is to expand the breadth of training materials to include fundamental statistical and data science topics and to precipitate a complete re-engineering of undergraduate and graduate curricula in life sciences. This project is accessible at https://training.galaxyproject.org.

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.