YY1-mediated enhancer-promoter communication in the immunoglobulin µ locus is regulated by MSL/MOF recruitment.

The rearrangement and expression of the immunoglobulin µ heavy chain (Igh) gene require communication of the intragenic Eµ and 3' regulatory region (RR) enhancers with the variable (VH) gene promoter. Eµ binding of the transcription factor YY1 has been implicated in enhancer-promoter communication, but the YY1 protein network remains obscure. By analyzing the comprehensive proteome of the 1-kb Eµ wild-type enhancer and that of Eµ lacking the YY1 binding site, we identified the male-specific lethal (MSL)/MOF complex as a component of the YY1 protein network. We found that MSL2 recruitment depends on YY1 and that gene knockout of Msl2 in primary pre-B cells reduces µ gene expression and chromatin looping of Eµ to the 3' RR enhancer and VH promoter. Moreover, Mof heterozygosity in mice impaired µ expression and early B cell differentiation. Together, these data suggest that the MSL/MOF complex regulates Igh gene expression by augmenting YY1-mediated enhancer-promoter communication.

Antagonistic interactions safeguard mitotic propagation of genetic and epigenetic information in zebrafish.

The stability of cellular phenotypes in developing organisms depends on error-free transmission of epigenetic and genetic information during mitosis. Methylation of cytosine residues in genomic DNA is a key epigenetic mark that modulates gene expression and prevents genome instability. Here, we report on a genetic test of the relationship between DNA replication and methylation in the context of the developing vertebrate organism instead of cell lines. Our analysis is based on the identification of hypomorphic alleles of dnmt1, encoding the DNA maintenance methylase Dnmt1, and pole1, encoding the catalytic subunit of leading-strand DNA polymerase epsilon holoenzyme (Pole). Homozygous dnmt1 mutants exhibit genome-wide DNA hypomethylation, whereas the pole1 mutation is associated with increased DNA methylation levels. In dnmt1/pole1 double-mutant zebrafish larvae, DNA methylation levels are restored to near normal values, associated with partial rescue of mutant-associated transcriptional changes and phenotypes. Hence, a balancing antagonism between DNA replication and maintenance methylation buffers against replicative errors contributing to the robustness of vertebrate development.

Single-cell deconvolution reveals high lineage- and location-dependent heterogeneity in mesenchymal multivisceral stage 4 colorectal cancer.

Metastasized colorectal cancer (CRC) is associated with a poor prognosis and rapid disease progression. Besides hepatic metastasis, peritoneal carcinomatosis is the major cause of death in Union for International Cancer Control (UICC) stage IV CRC patients. Insights into differential site-specific reconstitution of tumor cells and the corresponding tumor microenvironment are still missing. Here, we analyzed the transcriptome of single cells derived from murine multivisceral CRC and delineated the intermetastatic cellular heterogeneity regarding tumor epithelium, stroma, and immune cells. Interestingly, we found an intercellular site-specific network of cancer-associated fibroblasts and tumor epithelium during peritoneal metastasis as well as an autologous feed-forward loop in cancer stem cells. We furthermore deciphered a metastatic dysfunctional adaptive immunity by a loss of B cell-dependent antigen presentation and consecutive effector T cell exhaustion. Furthermore, we demonstrated major similarities of this murine metastatic CRC model with human disease and - based on the results of our analysis - provided an auspicious site-specific immunomodulatory treatment approach for stage IV CRC by intraperitoneal checkpoint inhibition.

Transcriptional reprogramming by mutated IRF4 in lymphoma.

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

Tnpo3 enables EBF1 function in conditions of antagonistic Notch signaling.

Transcription factor EBF1 (early B cell factor 1) acts as a key regulator of B cell specification. The transcriptional network in which EBF1 operates has been extensively studied; however, the regulation of EBF1 function remains poorly defined. By mass spectrometric analysis of proteins associated with endogenous EBF1 in pro-B cells, we identified the nuclear import receptor Transportin-3 (Tnpo3) and found that it interacts with the immunoglobulin-like fold domain of EBF1. We delineated glutamic acid 271 of EBF1 as a critical residue for the association with Tnpo3. EBF1E271A showed normal nuclear localization; however, it had an impaired B cell programming ability in conditions of Notch signaling, as determined by retroviral transduction of Ebf1 -/- progenitors. By RNA-seq analysis of EBF1E271A-expressing progenitors, we found an up-regulation of T lineage determinants and down-regulation of early B genes, although similar chromatin binding of EBF1E271A and EBF1wt was detected in pro-B cells expressing activated Notch1. B lineage-specific inactivation of Tnpo3 in mice resulted in a block of early B cell differentiation, accompanied by a down-regulation of B lineage genes and up-regulation of T and NK lineage genes. Taken together, our observations suggest that Tnpo3 ensures B cell programming by EBF1 in nonpermissive conditions.

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.

Senescent Tumor Cells in the Peritoneal Carcinomatosis Drive Immunosenescence in the Tumor Microenvironment.

More than half of all patients with colorectal cancer (CRC) develop distant metastasis and, depending on the local stage of the primary tumor, up to 48% of patients present peritoneal carcinomatosis (PC). PC is often considered as a widespread metastatic disease, which is almost resistant to current systemic therapies like chemotherapeutic and immunotherapeutic regimens. Here we could show that tumor cells of PC besides being senescent also exhibit stem cell features. To investigate these surprising findings in more detail, we established a murine model based on tumor organoids that resembles the clinical setting. In this murine orthotopic transplantation model for peritoneal carcinomatosis, we could show that the metastatic site in the peritoneum is responsible for senescence and stemness induction in tumor cells and that induction of senescence is not due to oncogene activation or therapy. In both mouse and human PC, senescence is associated with a senescence-associated secretory phenotype (SASP) influencing the tumor microenvironment (TME) of PC. SASP factors are able to induce a senescence phenotype in neighbouring cells. Here we could show that SASP leads to enhanced immunosenescence in the TME of PC. Our results provide a new immunoescape mechanism in PC explaining the resistance of PC to known chemo- and immunotherapeutic approaches. Therefore, senolytic approaches may represent a novel roadmap to target this terminal stage of CRC.

Epigenetic regulator genes direct lineage switching in MLL/AF4 leukemia.

The fusion gene MLL/AF4 defines a high-risk subtype of pro-B acute lymphoblastic leukemia. Relapse can be associated with a lineage switch from acute lymphoblastic to acute myeloid leukemia, resulting in poor clinical outcomes caused by resistance to chemotherapies and immunotherapies. In this study, the myeloid relapses shared oncogene fusion breakpoints with their matched lymphoid presentations and originated from various differentiation stages from immature progenitors through to committed B-cell precursors. Lineage switching is linked to substantial changes in chromatin accessibility and rewiring of transcriptional programs, including alternative splicing. These findings indicate that the execution and maintenance of lymphoid lineage differentiation is impaired. The relapsed myeloid phenotype is recurrently associated with the altered expression, splicing, or mutation of chromatin modifiers, including CHD4 coding for the ATPase/helicase of the nucleosome remodelling and deacetylation complex. Perturbation of CHD4 alone or in combination with other mutated epigenetic modifiers induces myeloid gene expression in MLL/AF4+ cell models, indicating that lineage switching in MLL/AF4 leukemia is driven and maintained by disrupted epigenetic regulation.

Sequences to Differences in Gene Expression: Analysis of RNA-Seq Data.

RNA-Seq is now a routinely employed assay to measure gene expression. As the technique matured over the last decade, so have dedicated analytic tools. In this chapter, we first describe the mainstream as well as the most up-to-date protocols and their implications on downstream analysis. We then detail the steps entailing RNA-Seq analysis in three main stages: (i) preprocessing and data preparation, (ii) upstream processing, and (iii) high-level analyses. We review the most recent and relevant tools as one workflow following a stepwise order. The chapter further encompasses in-depth features of these tools. Details of the required code are made available throughout the chapter, as well as of the underlying statistics. We illustrate these steps with analysis of publicly available RNA-Seq data.

Multilayer omics analysis reveals a non-classical retinoic acid signaling axis that regulates hematopoietic stem cell identity.

Hematopoietic stem cells (HSCs) rely on complex regulatory networks to preserve stemness. Due to the scarcity of HSCs, technical challenges have limited our insights into the interplay between metabolites, transcription, and the epigenome. In this study, we generated low-input metabolomics, transcriptomics, chromatin accessibility, and chromatin immunoprecipitation data, revealing distinct metabolic hubs that are enriched in HSCs and their downstream multipotent progenitors. Mechanistically, we uncover a non-classical retinoic acid (RA) signaling axis that regulates HSC function. We show that HSCs rely on Cyp26b1, an enzyme conventionally considered to limit RA effects in the cell. In contrast to the traditional view, we demonstrate that Cyp26b1 is indispensable for production of the active metabolite 4-oxo-RA. Further, RA receptor beta (Rarb) is required for complete transmission of 4-oxo-RA-mediated signaling to maintain stem cells. Our findings emphasize that a single metabolite controls stem cell fate by instructing epigenetic and transcriptional attributes.

ZFP451-mediated SUMOylation of SATB2 drives embryonic stem cell differentiation.

The establishment of cell fates involves alterations of transcription factor repertoires and repurposing of transcription factors by post-translational modifications. In embryonic stem cells (ESCs), the chromatin organizers SATB2 and SATB1 balance pluripotency and differentiation by activating and repressing pluripotency genes, respectively. Here, we show that conditional Satb2 gene inactivation weakens ESC pluripotency, and we identify SUMO2 modification of SATB2 by the E3 ligase ZFP451 as a potential driver of ESC differentiation. Mutations of two SUMO-acceptor lysines of Satb2 (Satb2K →R ) or knockout of Zfp451 impair the ability of ESCs to silence pluripotency genes and activate differentiation-associated genes in response to retinoic acid (RA) treatment. Notably, the forced expression of a SUMO2-SATB2 fusion protein in either Satb2K →R or Zfp451-/- ESCs rescues, in part, their impaired differentiation potential and enhances the down-regulation of Nanog The differentiation defect of Satb2K →R ESCs correlates with altered higher-order chromatin interactions relative to Satb2wt ESCs. Upon RA treatment of Satb2wt ESCs, SATB2 interacts with ZFP451 and the LSD1/CoREST complex and gains binding at differentiation genes, which is not observed in RA-treated Satb2K →R cells. Thus, SATB2 SUMOylation may contribute to the rewiring of transcriptional networks and the chromatin interactome of ESCs in the transition of pluripotency to differentiation.

Repetitive Elements Trigger RIG-I-like Receptor Signaling that Regulates the Emergence of Hematopoietic Stem and Progenitor Cells.

Inflammatory signaling is required for hematopoietic stem and progenitor cell (HSPC) development. Here, we studied the involvement of RIG-I-like receptors (RLRs) in HSPC formation. Rig-I or Mda5 deficiency impaired, while Lgp2 deficiency enhanced, HSPC emergence in zebrafish embryos. Rig-I or Mda5 deficiency reduced HSPC numbers by inhibiting inflammatory signals that were in turn enhanced in Lgp2 deficient embryos. Simultaneous reduction of Lgp2 and either Rig-I or Mda5 rescued inflammatory signals and HSPC numbers. Modulating the expression of the signaling mediator Traf6 in RLR deficient embryos restored HSPC numbers. Repetitive element transcripts could be detected in hemogenic endothelial cells and HSPCs, suggesting a role as RLR ligands. Indeed, ectopic expression of repetitive elements enhanced HSPC formation in wild-type, but not in Rig-I or Mda5 deficient embryos. Manipulation of RLR expression in mouse fetal liver HSPCs indicated functional conservation among species. Thus, repetitive elements transcribed during development drive RLR-mediated inflammatory signals that regulate HSPC formation.

Interactions between lineage-associated transcription factors govern haematopoietic progenitor states.

Recent advances in molecular profiling provide descriptive datasets of complex differentiation landscapes including the haematopoietic system, but the molecular mechanisms defining progenitor states and lineage choice remain ill-defined. Here, we employed a cellular model of murine multipotent haematopoietic progenitors (Hoxb8-FL) to knock out 39 transcription factors (TFs) followed by RNA-Seq analysis, to functionally define a regulatory network of 16,992 regulator/target gene links. Focussed analysis of the subnetworks regulated by the B-lymphoid TF Ebf1 and T-lymphoid TF Gata3 revealed a surprising role in common activation of an early myeloid programme. Moreover, Gata3-mediated repression of Pax5 emerges as a mechanism to prevent precocious B-lymphoid differentiation, while Hox-mediated activation of Meis1 suppresses myeloid differentiation. To aid interpretation of large transcriptomics datasets, we also report a new method that visualises likely transitions that a progenitor will undergo following regulatory network perturbations. Taken together, this study reveals how molecular network wiring helps to establish a multipotent progenitor state, with experimental approaches and analysis tools applicable to dissecting a broad range of both normal and perturbed cellular differentiation landscapes.

ERBB2b mRNA isoform encodes a nuclear variant of the ERBB2 oncogene in breast cancer.

The presence of nuclear ERBB2 receptor-type tyrosine kinase is one of the causes of the resistance to membrane ERBB2-targeted therapy in breast cancers. It has been previously reported that this nuclear location arises through at least two different mechanisms: proteolytic shedding of the extracellular domain of the full-length receptor and translation of the messenger RNA (mRNA)-encoding ERBB2 from internal initiation codons. Here, we report a new mechanism and function where a significant portion of nuclear ERBB2 results from the translation of the variant ERBB2 mRNA under the transcriptional control of a distal promoter that is actively used in breast cancer cells. We show that both membrane ERBB2a and nuclear ERBB2b isoforms are prevalently expressed in breast cancer cell lines and carcinoma samples. The ERBB2b isoform, which is translated from mRNA variant 2, can directly translocate into the nucleus due to the lack of the signal peptide which is required for an intermediate membrane location. Small interfering RNA-mediated gene silencing showed that ERBB2b can repress ERBB2a expression, encoded by variant 1, whereas ERBB2a activates ERBB2b. Nuclear ERBB2 binding to its own promoter was revealed by chromatin immunoprecipitation assay. Altogether, our results provide new insights into the origin and function of nuclear ERBB2 where it can participate at the same time in a positive or a negative feedback autoregulatory loop, dependent on which of its promoters this bona fide transcription factor is acting. They also provide a new understanding for the resistance to therapies targeting the membrane-anchored ERBB2 in breast cancer.

EBF1-deficient bone marrow stroma elicits persistent changes in HSC potential.

Crosstalk between mesenchymal stromal cells (MSCs) and hematopoietic stem cells (HSCs) is essential for hematopoietic homeostasis and lineage output. Here, we investigate how transcriptional changes in bone marrow (BM) MSCs result in long-lasting effects on HSCs. Single-cell analysis of Cxcl12-abundant reticular (CAR) cells and PDGFRα+Sca1+ (PαS) cells revealed an extensive cellular heterogeneity but uniform expression of the transcription factor gene Ebf1. Conditional deletion of Ebf1 in these MSCs altered their cellular composition, chromatin structure and gene expression profiles, including the reduced expression of adhesion-related genes. Functionally, the stromal-specific Ebf1 inactivation results in impaired adhesion of HSCs, leading to reduced quiescence and diminished myeloid output. Most notably, HSCs residing in the Ebf1-deficient niche underwent changes in their cellular composition and chromatin structure that persist in serial transplantations. Thus, genetic alterations in the BM niche lead to long-term functional changes of HSCs.

2HybridTools, a handy software to facilitate clone identification and mutation mapping from yeast two-hybrid screening.

Yeast Two-Hybrid (Y2H) and reverse Two-Hybrid (RY2H) are powerful protein-protein interaction screening methods that rely on the interaction of bait and prey proteins fused to DNA binding (DB) and activation domains (AD), respectively. Y2H allows identification of protein interaction partners using screening libraries, while RY2H is used to determine residues critical to a given protein-protein interaction by exploiting site-directed mutagenesis. Currently, both these techniques still rely on sequencing of positive clones using conventional Sanger sequencing. For Y2H, a screen can yield several positives; the identification of such clones is further complicated by the fact that sequencing products usually contain vector sequence. For RY2H, obtaining a complete sequence is required to identify the full range of residues involved in protein-protein interactions. However, with Sanger sequencing limited to 500-800 nucleotides, sequencing is usually carried from both ends for clones greater than this length. Analysis of such RY2H data thus requires assembly of sequencing products combined with trimming of vector sequences and of low-quality bases at the beginning and ends of sequencing products. Further, RY2H analysis requires collation of mutations that abrogate a DB/AD interaction. Here, we present 2HybridTools, a Java program with a user-friendly interface that allows addressing all these issues inherent to both Y2H and RY2H. Specifically, for Y2H, 2HybridTools enables automated identification of positive clones, while for RY2H, 2HybridTools provides detailed mutation reports as a basis for further investigation of given protein-protein interactions.

Dependence on Myb expression is attenuated in myeloid leukaemia with N-terminal CEBPA mutations.

Mutations at the N- or C-terminus of C/EBPα are frequent in acute myeloid leukaemia (AML) with normal karyotype. Here, we investigate the role of the transcription factor Myb in AMLs driven by different combinations of CEBPA mutations. Using knockdown of Myb in murine cell lines modelling the spectrum of CEBPA mutations, we show that the effect of reduced Myb depends on the mutational status of the two Cebpa alleles. Importantly, Myb knockdown fails to override the block in myeloid differentiation in cells with biallelic N-terminal C/EBPα mutations, demonstrating for the first time that the dependency on Myb is much lower in AML with this mutational profile. By comparing gene expression following Myb knockdown and chromatin immunoprecipitation sequencing data for the binding of C/EBPα isoforms, we provide evidence for a functional cooperation between C/EBPα and Myb in the maintenance of AML. This co-dependency breaks down when both alleles of CEBPA harbour N-terminal mutations, as a subset of C/EBPα-regulated genes only bind the short p30 C/EBPα isoform and, unlike other C/EBPα-regulated genes, do so without a requirement for Myb.

Subtype-specific regulatory network rewiring in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a heterogeneous disease caused by a variety of alterations in transcription factors, epigenetic regulators and signaling molecules. To determine how different mutant regulators establish AML subtype-specific transcriptional networks, we performed a comprehensive global analysis of cis-regulatory element activity and interaction, transcription factor occupancy and gene expression patterns in purified leukemic blast cells. Here, we focused on specific subgroups of subjects carrying mutations in genes encoding transcription factors (RUNX1, CEBPα), signaling molecules (FTL3-ITD, RAS) and the nuclear protein NPM1). Integrated analysis of these data demonstrates that each mutant regulator establishes a specific transcriptional and signaling network unrelated to that seen in normal cells, sustaining the expression of unique sets of genes required for AML growth and maintenance.

Global long terminal repeat activation participates in establishing the unique gene expression programme of classical Hodgkin lymphoma.

Long terminal repeat (LTR) elements are wide-spread in the human genome and have the potential to act as promoters and enhancers. Their expression is therefore under tight epigenetic control. We previously reported in classical Hodgkin Lymphoma (cHL) that a member of the THE1B class of LTR elements acted as a promoter for the proto-oncogene and growth factor receptor gene CSF1R and that expression of this gene is required for cHL tumour survival. However, to which extent and how such elements participate in globally shaping the unique cHL gene expression programme is unknown. To address this question we mapped the genome-wide activation of THE1-LTRs in cHL cells using a targeted next generation sequencing approach (RACE-Seq). Integration of these data with global gene expression data from cHL and control B cell lines showed a unique pattern of LTR activation impacting on gene expression, including genes associated with the cHL phenotype. We also show that global LTR activation is induced by strong inflammatory stimuli. Together these results demonstrate that LTR activation provides an additional layer of gene deregulation in classical Hodgkin lymphoma and highlight the potential impact of genome-wide LTR activation in other inflammatory diseases.