Histone H3.3 lysine 9 and 27 control repressive chromatin at cryptic enhancers and bivalent promoters.

Histone modifications are associated with distinct transcriptional states, but it is unclear whether they instruct gene expression. To investigate this, we mutate histone H3.3 K9 and K27 residues in mouse embryonic stem cells (mESCs). Here, we find that H3.3K9 is essential for controlling specific distal intergenic regions and for proper H3K27me3 deposition at promoters. The H3.3K9A mutation resulted in decreased H3K9me3 at regions encompassing endogenous retroviruses and induced a gain of H3K27ac and nascent transcription. These changes in the chromatin environment unleash cryptic enhancers, resulting in the activation of distinctive transcriptional programs and culminating in protein expression normally restricted to specialized immune cell types. The H3.3K27A mutant disrupts the deposition and spreading of the repressive H3K27me3 mark, particularly impacting bivalent genes with higher basal levels of H3.3 at promoters. Therefore, H3.3K9 and K27 crucially orchestrate repressive chromatin states at cis-regulatory elements and bivalent promoters, respectively, and instruct proper transcription in mESCs.

Basis of gene-specific transcription regulation by the Integrator complex.

The Integrator complex attenuates gene expression via the premature termination of RNA polymerase II (RNAP2) at promoter-proximal pausing sites. It is required for stimulus response, cell differentiation, and neurodevelopment, but how gene-specific and adaptive regulation by Integrator is achieved remains unclear. Here, we identify two sites on human Integrator subunits 13/14 that serve as binding hubs for sequence-specific transcription factors (TFs) and other transcription effector complexes. When Integrator is attached to paused RNAP2, these hubs are positioned upstream of the transcription bubble, consistent with simultaneous TF-promoter tethering. The TFs co-localize with Integrator genome-wide, increase Integrator abundance on target genes, and co-regulate responsive transcriptional programs. For instance, sensory cilia formation induced by glucose starvation depends on Integrator-TF contacts. Our data suggest TF-mediated promoter recruitment of Integrator as a widespread mechanism for targeted transcription regulation.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA.

In placental females, one copy of the two X chromosomes is largely silenced during a narrow developmental time window, in a process mediated by the non-coding RNA Xist1. Here, we demonstrate that Xist can initiate X-chromosome inactivation (XCI) well beyond early embryogenesis. By modifying its endogenous level, we show that Xist has the capacity to actively silence genes that escape XCI both in neuronal progenitor cells (NPCs) and in vivo, in mouse embryos. We also show that Xist plays a direct role in eliminating TAD-like structures associated with clusters of escapee genes on the inactive X chromosome, and that this is dependent on Xist's XCI initiation partner, SPEN2. We further demonstrate that Xist's function in suppressing gene expression of escapees and topological domain formation is reversible for up to seven days post-induction, but that sustained Xist up-regulation leads to progressively irreversible silencing and CpG island DNA methylation of facultative escapees. Thus, the distinctive transcriptional and regulatory topologies of the silenced X chromosome is actively, directly - and reversibly - controlled by Xist RNA throughout life.

APOBEC2 safeguards skeletal muscle cell fate through binding chromatin and regulating transcription of non-muscle genes during myoblast differentiation.

The apolipoprotein B messenger RNA editing enzyme, catalytic polypeptide (APOBEC) family is composed of nucleic acid editors with roles ranging from antibody diversification to RNA editing. APOBEC2, a member of this family with an evolutionarily conserved nucleic acid-binding cytidine deaminase domain, has neither an established substrate nor function. Using a cellular model of muscle differentiation where APOBEC2 is inducibly expressed, we confirmed that APOBEC2 does not have the attributed molecular functions of the APOBEC family, such as RNA editing, DNA demethylation, and DNA mutation. Instead, we found that during muscle differentiation APOBEC2 occupied a specific motif within promoter regions; its removal from those regions resulted in transcriptional changes. Mechanistically, these changes reflect the direct interaction of APOBEC2 with histone deacetylase (HDAC) transcriptional corepressor complexes. We also found that APOBEC2 could bind DNA directly, in a sequence-specific fashion, suggesting that it functions as a recruiter of HDAC to specific genes whose promoters it occupies. These genes are normally suppressed during muscle cell differentiation, and their suppression may contribute to the safeguarding of muscle cell fate. Altogether, our results reveal a unique role for APOBEC2 within the APOBEC family.

Cytotoxic CD8+ Temra cells show loss of chromatin accessibility at genes associated with T cell activation.

As humans age, their memory T cell compartment expands due to the lifelong exposure to antigens. This expansion is characterized by terminally differentiated CD8+ T cells (Temra), which possess NK cell-like phenotype and are associated with chronic inflammatory conditions. Temra cells are predominantly driven by the sporadic reactivation of cytomegalovirus (CMV), yet their epigenomic patterns and cellular heterogeneity remain understudied. To address this gap, we correlated their gene expression profiles with chromatin openness and conducted single-cell transcriptome analysis, comparing them to other CD8+ subsets and CMV-responses. We confirmed that Temra cells exhibit high expression of genes associated with cytotoxicity and lower expression of costimulatory and chemokine genes. The data revealed that CMV-responsive CD8+ T cells (Tcmv) were predominantly derived from a mixed population of Temra and memory cells (Tcm/em) and shared their transcriptomic profiles. Using ATAC-seq analysis, we identified 1449 differentially accessible chromatin regions between CD8+ Temra and Tcm/em cells, of which only 127 sites gained chromatin accessibility in Temra cells. We further identified 51 gene loci, including costimulatory CD27, CD28, and ICOS genes, whose chromatin accessibility correlated with their gene expression. The differential chromatin regions Tcm/em cells were enriched in motifs that bind multiple transcriptional activators, such as Jun/Fos, NFkappaB, and STAT, whereas the open regions in Temra cells mainly contained binding sites of T-box transcription factors. Our single-cell analysis of CD8+CCR7loCD45RAhi sorted Temra population showed several subsets of Temra and NKT-like cells and CMC1+ Temra populations in older individuals that were shifted towards decreased cytotoxicity. Among CD8+CCR7loCD45RAhi sorted cells, we found a decreased proportion of IL7R+ Tcm/em-like and MAIT cells in individuals with high levels of CMV antibodies (CMVhi). These results shed new light on the molecular and cellular heterogeneity of CD8+ Temra cells and their relationship to aging and CMV infection.

The remission status of AML patients after allo-HCT is associated with a distinct single-cell bone marrow T-cell signature.

ABSTRACT: Acute myeloid leukemia (AML) is a hematologic malignancy for which allogeneic hematopoietic cell transplantation (allo-HCT) often remains the only curative therapeutic approach. However, incapability of T cells to recognize and eliminate residual leukemia stem cells might lead to an insufficient graft-versus-leukemia (GVL) effect and relapse. Here, we performed single-cell RNA-sequencing (scRNA-seq) on bone marrow (BM) T lymphocytes and CD34+ cells of 6 patients with AML 100 days after allo-HCT to identify T-cell signatures associated with either imminent relapse (REL) or durable complete remission (CR). We observed a higher frequency of cytotoxic CD8+ effector and gamma delta (γδ) T cells in CR vs REL samples. Pseudotime and gene regulatory network analyses revealed that CR CD8+ T cells were more advanced in maturation and had a stronger cytotoxicity signature, whereas REL samples were characterized by inflammatory tumor necrosis factor/NF-κB signaling and an immunosuppressive milieu. We identified ADGRG1/GPR56 as a surface marker enriched in CR CD8+ T cells and confirmed in a CD33-directed chimeric antigen receptor T cell/AML coculture model that GPR56 becomes upregulated on T cells upon antigen encounter and elimination of AML cells. We show that GPR56 continuously increases at the protein level on CD8+ T cells after allo-HCT and confirm faster interferon gamma (IFN-γ) secretion upon re-exposure to matched, but not unmatched, recipient AML cells in the GPR56+ vs GPR56- CD8+ T-cell fraction. Together, our data provide a single-cell reference map of BM-derived T cells after allo-HCT and propose GPR56 expression dynamics as a surrogate for antigen encounter after allo-HCT.

GRaNIE and GRaNPA: inference and evaluation of enhancer-mediated gene regulatory networks.

Enhancers play a vital role in gene regulation and are critical in mediating the impact of noncoding genetic variants associated with complex traits. Enhancer activity is a cell-type-specific process regulated by transcription factors (TFs), epigenetic mechanisms and genetic variants. Despite the strong mechanistic link between TFs and enhancers, we currently lack a framework for jointly analysing them in cell-type-specific gene regulatory networks (GRN). Equally important, we lack an unbiased way of assessing the biological significance of inferred GRNs since no complete ground truth exists. To address these gaps, we present GRaNIE (Gene Regulatory Network Inference including Enhancers) and GRaNPA (Gene Regulatory Network Performance Analysis). GRaNIE (https://git.embl.de/grp-zaugg/GRaNIE) builds enhancer-mediated GRNs based on covariation of chromatin accessibility and RNA-seq across samples (e.g. individuals), while GRaNPA (https://git.embl.de/grp-zaugg/GRaNPA) assesses the performance of GRNs for predicting cell-type-specific differential expression. We demonstrate their power by investigating gene regulatory mechanisms underlying the response of macrophages to infection, cancer and common genetic traits including autoimmune diseases. Finally, our methods identify the TF PURA as a putative regulator of pro-inflammatory macrophage polarisation.

MYT1L haploinsufficiency in human neurons and mice causes autism-associated phenotypes that can be reversed by genetic and pharmacologic intervention.

MYT1L is an autism spectrum disorder (ASD)-associated transcription factor that is expressed in virtually all neurons throughout life. How MYT1L mutations cause neurological phenotypes and whether they can be targeted remains enigmatic. Here, we examine the effects of MYT1L deficiency in human neurons and mice. Mutant mice exhibit neurodevelopmental delays with thinner cortices, behavioural phenotypes, and gene expression changes that resemble those of ASD patients. MYT1L target genes, including WNT and NOTCH, are activated upon MYT1L depletion and their chemical inhibition can rescue delayed neurogenesis in vitro. MYT1L deficiency also causes upregulation of the main cardiac sodium channel, SCN5A, and neuronal hyperactivity, which could be restored by shRNA-mediated knockdown of SCN5A or MYT1L overexpression in postmitotic neurons. Acute application of the sodium channel blocker, lamotrigine, also rescued electrophysiological defects in vitro and behaviour phenotypes in vivo. Hence, MYT1L mutation causes both developmental and postmitotic neurological defects. However, acute intervention can normalise resulting electrophysiological and behavioural phenotypes in adulthood.

Convolutional networks for supervised mining of molecular patterns within cellular context.

Cryo-electron tomograms capture a wealth of structural information on the molecular constituents of cells and tissues. We present DeePiCt (deep picker in context), an open-source deep-learning framework for supervised segmentation and macromolecular complex localization in cryo-electron tomography. To train and benchmark DeePiCt on experimental data, we comprehensively annotated 20 tomograms of Schizosaccharomyces pombe for ribosomes, fatty acid synthases, membranes, nuclear pore complexes, organelles, and cytosol. By comparing DeePiCt to state-of-the-art approaches on this dataset, we show its unique ability to identify low-abundance and low-density complexes. We use DeePiCt to study compositionally distinct subpopulations of cellular ribosomes, with emphasis on their contextual association with mitochondria and the endoplasmic reticulum. Finally, applying pre-trained networks to a HeLa cell tomogram demonstrates that DeePiCt achieves high-quality predictions in unseen datasets from different biological species in a matter of minutes. The comprehensively annotated experimental data and pre-trained networks are provided for immediate use by the community.

Transcriptional profiling unveils molecular subgroups of adaptive and maladaptive right ventricular remodeling in pulmonary hypertension.

Right ventricular (RV) function is critical to prognosis in all forms of pulmonary hypertension. Here we perform molecular phenotyping of RV remodeling by transcriptome analysis of RV tissue obtained from 40 individuals, and two animal models of RV dysfunction of both sexes. Our unsupervised clustering analysis identified 'early' and 'late' subgroups within compensated and decompensated states, characterized by the expression of distinct signaling pathways, while fatty acid metabolism and estrogen response appeared to underlie sex-specific differences in RV adaptation. The circulating levels of several extracellular matrix proteins deregulated in decompensated RV subgroups were assessed in two independent cohorts of individuals with pulmonary arterial hypertension, revealing that NID1, C1QTNF1 and CRTAC1 predicted the development of a maladaptive RV state, as defined by magnetic resonance imaging parameters, and were associated with worse clinical outcomes. Our study provides a resource for subphenotyping RV states, identifying state-specific biomarkers, and potential therapeutic targets for RV dysfunction.

Current challenges in understanding the role of enhancers in disease.

Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.

Identifying a novel role for the master regulator Tal1 in the Endothelial to Hematopoietic Transition.

Progress in the generation of Hematopoietic Stem and Progenitor Cells (HSPCs) in vitro and ex vivo has been built on the knowledge of developmental hematopoiesis, underscoring the importance of understanding this process. HSPCs emerge within the embryonic vasculature through an Endothelial-to-Hematopoietic Transition (EHT). The transcriptional regulator Tal1 exerts essential functions in the earliest stages of blood development, but is considered dispensable for the EHT. Nevertheless, Tal1 is expressed with its binding partner Lmo2 and it homologous Lyl1 in endothelial and transitioning cells at the time of EHT. Here, we investigated the function of these genes using a mouse embryonic-stem cell (mESC)-based differentiation system to model hematopoietic development. We showed for the first time that the expression of TAL1 in endothelial cells is crucial to ensure the efficiency of the EHT process and a sustained hematopoietic output. Our findings uncover an important function of Tal1 during the EHT, thus filling the current gap in the knowledge of the role of this master gene throughout the whole process of hematopoietic development.

Drug-microenvironment perturbations reveal resistance mechanisms and prognostic subgroups in CLL.

The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited and systematic analyses are lacking. Using chronic lymphocytic leukaemia (CLL) as a model disease, we investigated the influence of 17 microenvironmental stimuli on 12 drugs in 192 genetically characterised patient samples. Based on microenvironmental response, we identified four subgroups with distinct clinical outcomes beyond known prognostic markers. Response to multiple microenvironmental stimuli was amplified in trisomy 12 samples. Trisomy 12 was associated with a distinct epigenetic signature. Bromodomain inhibition reversed this epigenetic profile and could be used to target microenvironmental signalling in trisomy 12 CLL. We quantified the impact of microenvironmental stimuli on drug response and their dependence on genetic alterations, identifying interleukin 4 (IL4) and Toll-like receptor (TLR) stimulation as the strongest actuators of drug resistance. IL4 and TLR signalling activity was increased in CLL-infiltrated lymph nodes compared with healthy samples. High IL4 activity correlated with faster disease progression. The publicly available dataset can facilitate the investigation of cell-extrinsic mechanisms of drug resistance and disease progression.

Operation of a TCA cycle subnetwork in the mammalian nucleus.

Nucleic acid and histone modifications critically depend on the tricarboxylic acid (TCA) cycle for substrates and cofactors. Although a few TCA cycle enzymes have been reported in the nucleus, the corresponding pathways are considered to operate in mitochondria. Here, we show that a part of the TCA cycle is operational also in the nucleus. Using 13C-tracer analysis, we identified activity of glutamine-to-fumarate, citrate-to-succinate, and glutamine-to-aspartate routes in the nuclei of HeLa cells. Proximity labeling mass spectrometry revealed a spatial vicinity of the involved enzymes with core nuclear proteins. We further show nuclear localization of aconitase 2 and 2-oxoglutarate dehydrogenase in mouse embryonic stem cells. Nuclear localization of the latter enzyme, which produces succinyl-CoA, changed from pluripotency to a differentiated state with accompanying changes in the nuclear protein succinylation. Together, our results demonstrate operation of an extended metabolic pathway in the nucleus, warranting a revision of the canonical view on metabolic compartmentalization.

Comparative chromatin accessibility upon BDNF stimulation delineates neuronal regulatory elements.

Neuronal stimulation induced by the brain-derived neurotrophic factor (BDNF) triggers gene expression, which is crucial for neuronal survival, differentiation, synaptic plasticity, memory formation, and neurocognitive health. However, its role in chromatin regulation is unclear. Here, using temporal profiling of chromatin accessibility and transcription in mouse primary cortical neurons upon either BDNF stimulation or depolarization (KCl), we identify features that define BDNF-specific chromatin-to-gene expression programs. Enhancer activation is an early event in the regulatory control of BDNF-treated neurons, where the bZIP motif-binding Fos protein pioneered chromatin opening and cooperated with co-regulatory transcription factors (Homeobox, EGRs, and CTCF) to induce transcription. Deleting cis-regulatory sequences affect BDNF-mediated Arc expression, a regulator of synaptic plasticity. BDNF-induced accessible regions are linked to preferential exon usage by neurodevelopmental disorder-related genes and the heritability of neuronal complex traits, which were validated in human iPSC-derived neurons. Thus, we provide a comprehensive view of BDNF-mediated genome regulatory features using comparative genomic approaches to dissect mammalian neuronal stimulation.

Glucose Metabolism and Aging of Hematopoietic Stem and Progenitor Cells.

Comprehensive proteomics studies of human hematopoietic stem and progenitor cells (HSPC) have revealed that aging of the HSPC compartment is characterized by elevated glycolysis. This is in addition to deregulations found in murine transcriptomics studies, such as an increased differentiation bias towards the myeloid lineage, alterations in DNA repair, and a decrease in lymphoid development. The increase in glycolytic enzyme activity is caused by the expansion of a more glycolytic HSPC subset. We therefore developed a method to isolate HSPC into three distinct categories according to their glucose uptake (GU) levels, namely the GUhigh, GUinter and GUlow subsets. Single-cell transcriptomics studies showed that the GUhigh subset is highly enriched for HSPC with a differentiation bias towards myeloid lineages. Gene set enrichment analysis (GSEA) demonstrated that the gene sets for cell cycle arrest, senescence-associated secretory phenotype, and the anti-apoptosis and P53 pathways are significantly upregulated in the GUhigh population. With this series of studies, we have produced a comprehensive proteomics and single-cell transcriptomics atlas of molecular changes in human HSPC upon aging. Although many of the molecular deregulations are similar to those found in mice, there are significant differences. The most unique finding is the association of elevated central carbon metabolism with senescence. Due to the lack of specific markers, the isolation and collection of senescent cells have yet to be developed, especially for human HSPC. The GUhigh subset from the human HSPC compartment possesses all the transcriptome characteristics of senescence. This property may be exploited to accurately enrich, visualize, and trace senescence development in human bone marrow.

CDK7/12/13 inhibition targets an oscillating leukemia stem cell network and synergizes with venetoclax in acute myeloid leukemia.

The heterogeneous response of acute myeloid leukemia (AML) to current anti-leukemic therapies is only partially explained by mutational heterogeneity. We previously identified GPR56 as a surface marker associated with poor outcome across genetic groups, which characterizes two leukemia stem cell (LSC)-enriched compartments with different self-renewal capacities. How these compartments self-renew remained unclear. Here, we show that GPR56+ LSC compartments are promoted in a complex network involving epithelial-to-mesenchymal transition (EMT) regulators besides Rho, Wnt, and Hedgehog (Hh) signaling. Unexpectedly, Wnt pathway inhibition increased the more immature, slowly cycling GPR56+ CD34+ fraction and Hh/EMT gene expression, while Wnt activation caused opposite effects. Our data suggest that the crucial role of GPR56 lies in its ability to co-activate these opposing signals, thus ensuring the constant supply of both LSC subsets. We show that CDK7 inhibitors suppress both LSC-enriched subsets in vivo and synergize with the Bcl-2 inhibitor venetoclax. Our data establish reciprocal transition between LSC compartments as a novel concept underlying the poor outcome in GPR56high AML and propose combined CDK7 and Bcl-2 inhibition as LSC-directed therapy in this disease.

Genome-wide quantification of transcription factor binding at single-DNA-molecule resolution using methyl-transferase footprinting.

Precise control of gene expression requires the coordinated action of multiple factors at cis-regulatory elements. We recently developed single-molecule footprinting to simultaneously resolve the occupancy of multiple proteins including transcription factors, RNA polymerase II and nucleosomes on single DNA molecules genome-wide. The technique combines the use of cytosine methyltransferases to footprint the genome with bisulfite sequencing to resolve transcription factor binding patterns at cis-regulatory elements. DNA footprinting is performed by incubating permeabilized nuclei with recombinant methyltransferases. Upon DNA extraction, whole-genome or targeted bisulfite libraries are prepared and loaded on Illumina sequencers. The protocol can be completed in 4-5 d in any laboratory with access to high-throughput sequencing. Analysis can be performed in 2 d using a dedicated R package and requires access to a high-performance computing system. Our method can be used to analyze how transcription factors cooperate and antagonize to regulate transcription.

Enhancers in disease: molecular basis and emerging treatment strategies.

Enhancers are genomic sequences that play a key role in regulating tissue-specific gene expression levels. An increasing number of diseases are linked to impaired enhancer function through chromosomal rearrangement, genetic variation within enhancers, or epigenetic modulation. Here, we review how these enhancer disruptions have recently been implicated in congenital disorders, cancers, and common complex diseases and address the implications for diagnosis and treatment. Although further fundamental research into enhancer function, target genes, and context is required, enhancer-targeting drugs and gene editing approaches show great therapeutic promise for a range of diseases.

Transcription factors: Bridge between cell signaling and gene regulation.

Transcription factors (TFs) are key regulators of intrinsic cellular processes, such as differentiation and development, and of the cellular response to external perturbation through signaling pathways. In this review we focus on the role of TFs as a link between signaling pathways and gene regulation. Cell signaling tends to result in the modulation of a set of TFs that then lead to changes in the cell's transcriptional program. We highlight the molecular layers at which TF activity can be measured and the associated technical and conceptual challenges. These layers include post-translational modifications (PTMs) of the TF, regulation of TF binding to DNA through chromatin accessibility and epigenetics, and expression of target genes. We highlight that a large number of TFs are understudied in both signaling and gene regulation studies, and that our knowledge about known TF targets has a strong literature bias. We argue that TFs serve as a perfect bridge between the fields of gene regulation and signaling, and that separating these fields hinders our understanding of cell functions. Multi-omics approaches that measure multiple dimensions of TF activity are ideally suited to study the interplay of cell signaling and gene regulation using TFs as the anchor to link the two fields.