Two distinct auto-regulatory loops operate at the PU.1 locus in B cells and myeloid cells

Chromatin priming elements direct tissue-specific gene activity before hematopoietic specification

Tissue-specific gene regulation during development involves the interplay between transcription factors and epigenetic regulators binding to enhancer and promoter elements. The pattern of active enhancers defines the cellular differentiation state. However, developmental gene activation involves a previous step called chromatin priming which is not fully understood. We recently developed a genome-wide functional assay that allowed us to functionally identify enhancer elements integrated in chromatin regulating five stages spanning the in vitro differentiation of embryonic stem cells to blood. We also measured global chromatin accessibility, histone modifications, and transcription factor binding. The integration of these data identified and characterised cis-regulatory elements which become activated before the onset of gene expression, some of which are primed in a signalling-dependent fashion. Deletion of such a priming element leads to a delay in the up-regulation of its associated gene in development. Our work uncovers the details of a complex network of regulatory interactions with the dynamics of early chromatin opening being at the heart of dynamic tissue-specific gene expression control.

Chromatin structure and 3D architecture define differential functions of PU.1 cis regulatory elements in human blood cell lineages

The precise spatio-temporal expression of the hematopoietic ETS transcription factor PU.1 that determines the hematopoietic cell fates is tightly regulated at the chromatin level. However, it remains elusive as to how chromatin signatures are linked to this dynamic expression pattern of PU.1 across blood cell lineages. Here we performed an unbiased and in-depth analysis of the relationship between human PU.1 expression, the presence of trans-acting factors, and 3D architecture at various cis-regulatory elements (CRE) proximal to the PU.1 locus. We identified multiple novel CREs at the upstream region of the gene following an integrative inspection for conserved DNA elements at the chromatin-accessible regions in primary human blood lineages. We showed that a subset of CREs localize within a 10 kb-wide cluster that exhibits that exhibit molecular features of a myeloid-specific super-enhancer involved in mediating PU.1 autoregulation, including open chromatin, unmethylated DNA, histone enhancer marks, transcription of enhancer RNAs, and occupancy of the PU.1 protein itself. Importantly, we revealed the presence of common 35-kb-wide CTCF-bound insulated neighborhood that contains the CRE cluster, forming the chromatin territory for lineage-specific and CRE-mediated chromatin interactions. These include functional CRE-promoter interactions in myeloid and B cells but not in erythroid and T cells. Our findings also provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in defining certain CREs as enhancers or silencers in chromatin regulation of PU.1 expression. The study lays the groundwork for further examination of PU.1 CREs as well as epigenetic regulation in malignant hematopoiesis.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

PTEN regulates hematopoietic lineage plasticity via PU.1-dependent chromatin accessibility

PTEN loss in fetal liver hematopoietic stem cells (HSCs) leads to alterations in myeloid, T-, and B-lineage potentials and T-lineage acute lymphoblastic leukemia (T-ALL) development. To explore the mechanism underlying PTEN-regulated hematopoietic lineage choices, we carry out integrated assay for transposase-accessible chromatin using sequencing (ATAC-seq), single-cell RNA-seq, and in vitro culture analyses using in vivo-isolated mouse pre-leukemic HSCs and progenitors. We find that PTEN loss alters chromatin accessibility of key lineage transcription factor (TF) binding sites at the prepro-B stage, corresponding to increased myeloid and T-lineage potentials and reduced B-lineage potential. Importantly, we find that PU.1 is an essential TF downstream of PTEN and that altering PU.1 levels can reprogram the chromatin accessibility landscape and myeloid, T-, and B-lineage potentials in Ptennull prepro-B cells. Our study discovers prepro-B as the key developmental stage underlying PTEN-regulated hematopoietic lineage choices and suggests a critical role of PU.1 in modulating the epigenetic state and lineage plasticity of prepro-B progenitors.

Carm1-arginine methylation of the transcription factor C/EBPα regulates transdifferentiation velocity

Here, we describe how the speed of C/EBPα-induced B cell to macrophage transdifferentiation (BMT) can be regulated, using both mouse and human models. The identification of a mutant of C/EBPα (C/EBPαR35A) that greatly accelerates BMT helped to illuminate the mechanism. Thus, incoming C/EBPα binds to PU.1, an obligate partner expressed in B cells, leading to the release of PU.1 from B cell enhancers, chromatin closing and silencing of the B cell program. Released PU.1 redistributes to macrophage enhancers newly occupied by C/EBPα, causing chromatin opening and activation of macrophage genes. All these steps are accelerated by C/EBPαR35A, initiated by its increased affinity for PU.1. Wild-type C/EBPα is methylated by Carm1 at arginine 35 and the enzyme's perturbations modulate BMT velocity as predicted from the observations with the mutant. Increasing the proportion of unmethylated C/EBPα in granulocyte/macrophage progenitors by inhibiting Carm1 biases the cell's differentiation toward macrophages, suggesting that cell fate decision velocity and lineage directionality are closely linked processes.

Unveiling the complexity of transcription factor networks in hematopoietic stem cells: implications for cell therapy and hematological malignancies

The functionality and longevity of hematopoietic tissue is ensured by a tightly controlled balance between self-renewal, quiescence, and differentiation of hematopoietic stem cells (HSCs) into the many different blood lineages. Cell fate determination in HSCs is influenced by signals from extrinsic factors (e.g., cytokines, irradiation, reactive oxygen species, O2 concentration) that are translated and integrated by intrinsic factors such as Transcription Factors (TFs) to establish specific gene regulatory programs. TFs also play a central role in the establishment and/or maintenance of hematological malignancies, highlighting the need to understand their functions in multiple contexts. TFs bind to specific DNA sequences and interact with each other to form transcriptional complexes that directly or indirectly control the expression of multiple genes. Over the past decades, significant research efforts have unraveled molecular programs that control HSC function. This, in turn, led to the identification of more than 50 TF proteins that influence HSC fate. However, much remains to be learned about how these proteins interact to form molecular networks in combination with cofactors (e.g. epigenetics factors) and how they control differentiation, expansion, and maintenance of cellular identity. Understanding these processes is critical for future applications particularly in the field of cell therapy, as this would allow for manipulation of cell fate and induction of expansion, differentiation, or reprogramming of HSCs using specific cocktails of TFs. Here, we review recent findings that have unraveled the complexity of molecular networks controlled by TFs in HSCs and point towards possible applications to obtain functional HSCs ex vivo for therapeutic purposes including hematological malignancies. Furthermore, we discuss the challenges and prospects for the derivation and expansion of functional adult HSCs in the near future.

Cis interactions in the Irf8 locus regulate stage-dependent enhancer activation

Individual elements within a superenhancer can act in a cooperative or temporal manner, but the underlying mechanisms remain obscure. We recently identified an Irf8 superenhancer, within which different elements act at distinct stages of type 1 classical dendritic cell (cDC1) development. The +41-kb Irf8 enhancer is required for pre-cDC1 specification, while the +32-kb Irf8 enhancer acts to support subsequent cDC1 maturation. Here, we found that compound heterozygous Δ32/Δ41 mice, lacking the +32- and +41-kb enhancers on different chromosomes, show normal pre-cDC1 specification but, surprisingly, completely lack mature cDC1 development, suggesting cis dependence of the +32-kb enhancer on the +41-kb enhancer. Transcription of the +32-kb Irf8 enhancer-associated long noncoding RNA (lncRNA) Gm39266 is also dependent on the +41-kb enhancer. However, cDC1 development in mice remained intact when Gm39266 transcripts were eliminated by CRISPR/Cas9-mediated deletion of lncRNA promoters and when transcription across the +32-kb enhancer was blocked by premature polyadenylation. We showed that chromatin accessibility and BATF3 binding at the +32-kb enhancer were dependent on a functional +41-kb enhancer located in cis Thus, the +41-kb Irf8 enhancer controls the subsequent activation of the +32-kb Irf8 enhancer in a manner that is independent of associated lncRNA transcription.

A genome-wide relay of signalling-responsive enhancers drives hematopoietic specification

Developmental control of gene expression critically depends on distal cis-regulatory elements including enhancers which interact with promoters to activate gene expression. To date no global experiments have been conducted that identify their cell type and cell stage-specific activity within one developmental pathway and in a chromatin context. Here, we describe a high-throughput method that identifies thousands of differentially active cis-elements able to stimulate a minimal promoter at five stages of hematopoietic progenitor development from embryonic stem (ES) cells, which can be adapted to any ES cell derived cell type. We show that blood cell-specific gene expression is controlled by the concerted action of thousands of differentiation stage-specific sets of cis-elements which respond to cytokine signals terminating at signalling responsive transcription factors. Our work provides an important resource for studies of hematopoietic specification and highlights the mechanisms of how and where extrinsic signals program a cell type-specific chromatin landscape driving hematopoietic differentiation.

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.

Combined Approach to Leukemic Differentiation Using Transcription Factor PU.1-Enhancing Agents

The transcription factor PU.1 (Purine-rich DNA binding, SPI1) is a key regulator of hematopoiesis, whose level is influenced by transcription through its enhancers and its post-transcriptional degradation via microRNA-155 (miR-155). The degree of transcriptional regulation of the PU.1 gene is influenced by repression via DNA methylation, as well as other epigenetic factors, such as those related to progenitor maturation status, which is modulated by the transcription factor Myeloblastosis oncogene (MYB). In this work, we show that combinatorial treatment of acute myeloid leukemia (AML) cells with DNA methylation inhibitors (5-Azacytidine), MYB inhibitors (Celastrol), and anti-miR-155 (AM155) ideally leads to overproduction of PU.1. We also show that PU.1 reactivation can be compensated by miR-155 and that only a combined approach leads to sustained PU.1 derepression, even at the protein level. The triple effect on increasing PU.1 levels in myeloblasts stimulates the myeloid transcriptional program while inhibiting cell survival and proliferation, leading to partial leukemic differentiation.

Spi1 -14 Kb upstream regulatory element (URE) is not required for maintenance of PU.1 expression in macrophages

Background: Previous work suggested an upstream regulatory element (URE) of Spi1 was required to maintain constant expression of the PU.1 transcription factor in bone marrow and foetal liver cells. PU.1, encoded by Spi1, is essential for development and maintenance of myeloid and B-lymphocyte populations in mice. Deletion of this (-14 Kb) URE potentially reduces expression of PU.1 and therefore provides a way to investigate its role in myeloid populations in development and disease. This study aimed to examine the impact of removal of the -14 Kb Spi1 URE in Cx3cr1+ cells on the myeloid lineage formation and maintenance. Methods: B6;129-Spi1tm1.2Dgt/J mice, whose -14 Kb Spi1 URE mice is flanked by LoxP sites (‘floxed’), were bred to a strain with constitutively active Cre expressed under the Cx3cr1 promoter (B6J.B6N(Cg)-Cx3cr1tm1.1(cre)Jung/J) to delete the Spi1 URE in myeloid cells. The floxed mice were also bred to mice with a tamoxifen-inducible Cre expressed under the Cx3cr1 promoter (B6.129P2(C)-Cx3cr1tm2.1(cre/ERT2)Jung/J) to be used as URE intact controls and to permit temporally-controlled deletion of the URE if required. PU.1 protein expression was measured in the peritoneal macrophages and microglia by flow cytometry. Additionally, a Cre-encoding lentiviral vector was used to assess the impact on PU.1 expression in bone-marrow derived macrophages from these mice in vitro. Results: Expression of the PU.1 transcription factor was not significantly altered in the peritoneal macrophages or microglia in mice lacking the -14 Kb Spi1 URE. Moreover, initial experiments utilising Cre encoding lentivirus did not reduce PU.1 protein in bone-marrow derived macrophages differentiated from the -14 Kb Spi1 URE floxed mice. Conclusions: These observations suggest that the -14 Kb URE does not play a major role in PU.1 protein expression in either mature peritoneal macrophages or microglia.

Blood stem cell PU.1 upregulation is a consequence of differentiation without fast autoregulation

Transcription factors (TFs) regulate cell fates, and their expression must be tightly regulated. Autoregulation is assumed to regulate many TFs' own expression to control cell fates. Here, we manipulate and quantify the (auto)regulation of PU.1, a TF controlling hematopoietic stem and progenitor cells (HSPCs), and correlate it to their future fates. We generate transgenic mice allowing both inducible activation of PU.1 and noninvasive quantification of endogenous PU.1 protein expression. The quantified HSPC PU.1 dynamics show that PU.1 up-regulation occurs as a consequence of hematopoietic differentiation independently of direct fast autoregulation. In contrast, inflammatory signaling induces fast PU.1 up-regulation, which does not require PU.1 expression or its binding to its own autoregulatory enhancer. However, the increased PU.1 levels induced by inflammatory signaling cannot be sustained via autoregulation after removal of the signaling stimulus. We conclude that PU.1 overexpression induces HSC differentiation before PU.1 up-regulation, only later generating cell types with intrinsically higher PU.1.

Data-driven modeling predicts gene regulatory network dynamics during the differentiation of multipotential hematopoietic progenitors

Cellular differentiation during hematopoiesis is guided by gene regulatory networks (GRNs) comprising transcription factors (TFs) and the effectors of cytokine signaling. Based largely on analyses conducted at steady state, these GRNs are thought to be organized as a hierarchy of bistable switches, with antagonism between Gata1 and PU.1 driving red- and white-blood cell differentiation. Here, we utilize transient gene expression patterns to infer the genetic architecture—the type and strength of regulatory interconnections—and dynamics of a twelve-gene GRN including key TFs and cytokine receptors. We trained gene circuits, dynamical models that learn genetic architecture, on high temporal-resolution gene-expression data from the differentiation of an inducible cell line into erythrocytes and neutrophils. The model is able to predict the consequences of gene knockout, knockdown, and overexpression experiments and the inferred interconnections are largely consistent with prior empirical evidence. The inferred genetic architecture is densely interconnected rather than hierarchical, featuring extensive cross-antagonism between genes from alternative lineages and positive feedback from cytokine receptors. The analysis of the dynamics of gene regulation in the model reveals that PU.1 is one of the last genes to be upregulated in neutrophil conditions and that the upregulation of PU.1 and other neutrophil genes is driven by Cebpa and Gfi1 instead. This model inference is confirmed in an independent single-cell RNA-Seq dataset from mouse bone marrow in which Cebpa and Gfi1 expression precedes the neutrophil-specific upregulation of PU.1 during differentiation. These results demonstrate that full PU.1 upregulation during neutrophil development involves regulatory influences extrinsic to the Gata1-PU.1 bistable switch. Furthermore, although there is extensive cross-antagonism between erythroid and neutrophil genes, it does not have a hierarchical structure. More generally, we show that the combination of high-resolution time series data and data-driven dynamical modeling can uncover the dynamics and causality of developmental events that might otherwise be obscured.

The supply of blood cells is replenished by the maturation of hematopoietic progenitor cells into different cell types. Which cell type a progenitor cell develops into is determined by a complex network of genes whose protein products directly or indirectly regulate each others’ expression and that of downstream genes characteristic of the cell type. We inferred the nature and causality of the regulatory connections in a 12-gene network known to affect the decision between erythrocyte and neutrophil cell fates using a predictive machine-learning approach. Our analysis showed that the overall architecture of the network is densely interconnected and not hierarchical. Furthermore, the model inferred that PU.1, considered a master regulator of all white-blood cell lineages, is upregulated during neutrophil development by two other proteins, Cebpa and Gfi1. We validated this prediction by showing that Cebpa and Gfi1 expression precedes that of PU.1 in single-cell gene expression data from mouse bone marrow. These results revise the architecture of the gene network and the causality of regulatory events guiding hematopoiesis. The results also show that combining machine learning approaches with time course data can help resolve causality during development.

Genomic and epigenomic adaptation in SP-R210 (Myo18A) isoform-deficient macrophages

Macrophages play an important role in maintaining tissue homeostasis, from regulating the inflammatory response to pathogens to resolving inflammation and aiding tissue repair. The surfactant protein A (SP-A) receptor SP-R210 (MYO18A) has been shown to affect basal and inflammatory macrophage states. Specifically, disruption of the longer splice isoform SP-R210_L/MYO18Aα renders macrophages hyper-inflammatory, although the mechanism by which this occurs is not well understood. We asked whether disruption of the L isoform led to the hyper-inflammatory state via alteration of global genomic responses. RNA sequencing analysis of L isoform-deficient macrophages (SP-R210_L(DN)) revealed basal and influenza-induced upregulation of genes associated with inflammatory pathways, such as TLR, RIG-I, NOD, and cytoplasmic DNA signaling, whereas knockout of both SP-R210 isoforms (L and S) only resulted in increased RIG-I and NOD signaling. Chromatin immunoprecipitation sequencing (ChIP-seq) analysis showed increased genome-wide deposition of the pioneer transcription factor PU.1 in SP-R210_L(DN) cells, with increased representation around genes relevant to inflammatory pathways. Additional ChIP-seq analysis of histone H3 methylation marks showed decreases in both repressive H3K9me3 and H3K27me3 marks with a commensurate increase in transcriptionally active (H3K4me3) histone marks in the L isoform deficient macrophages. Influenza A virus (IAV) infection, known to stimulate a wide array of anti-viral responses, caused a differential redistribution of PU.1 binding between proximal promoter and distal sites and decoupling from Toll-like receptor regulated gene promoters in SP-R210_L(DN) cells. These finding suggest that the inflammatory differences seen in SP-R210_L-deficient macrophages are a result of transcriptional differences that are mediated by epigenetic changes brought about by differential expression of the SP-R210 isoforms. This provides an avenue to explore how the signaling pathways downstream of the receptor and the ligands can modulate the macrophage inflammatory response.

Modest changes in Spi1 dosage reveal the potential for altered microglial function as seen in Alzheimer's disease

Genetic association studies have identified multiple variants at the SPI1 locus that modify risk and age of onset for Alzheimer's Disease (AD). Reports linking risk variants to gene expression suggest that variants denoting higher SPI1 expression are likely to have an earlier AD onset, and several other AD risk genes contain PU.1 binding sites in the promoter region. Overall, this suggests the level of SPI1 may alter microglial phenotype potentially impacting AD. This study determined how the microglial transcriptome was altered following modest changes to Spi1 expression in primary mouse microglia. RNA-sequencing was performed on microglia with reduced or increased Spi1/PU.1 expression to provide an unbiased approach to determine transcriptomic changes affected by Spi1. In summary, a reduction in microglial Spi1 resulted in the dysregulation of transcripts encoding proteins involved in DNA replication pathways while an increased Spi1 results in an upregulation of genes associated with immune response pathways. Additionally, a subset of 194 Spi1 dose-sensitive genes was identified and pathway analysis suggests that several innate immune and interferon response pathways are impacted by the concentration of Spi1. Together these results suggest Spi1 levels can alter the microglial transcriptome and suggests interferon pathways may be altered in individuals with AD related Spi1 risk SNPs.

Expression of RUNX1-ETO Rapidly Alters the Chromatin Landscape and Growth of Early Human Myeloid Precursor Cells

Acute myeloid leukemia (AML) is a hematopoietic malignancy caused by recurrent mutations in genes encoding transcriptional, chromatin, and/or signaling regulators. The t(8;21) translocation generates the aberrant transcription factor RUNX1-ETO (RUNX1-RUNX1T1), which by itself is insufficient to cause disease. t(8;21) AML patients show extensive chromatin reprogramming and have acquired additional mutations. Therefore, the genomic and developmental effects directly and solely attributable to RUNX1-ETO expression are unclear. To address this, we employ a human embryonic stem cell differentiation system capable of forming definitive myeloid progenitor cells to express RUNX1-ETO in an inducible fashion. Induction of RUNX1-ETO causes extensive chromatin reprogramming by interfering with RUNX1 binding, blocks differentiation, and arrests cellular growth, whereby growth arrest is reversible following RUNX1-ETO removal. Single-cell gene expression analyses show that RUNX1-ETO induction alters the differentiation of early myeloid progenitors, but not of other progenitor types, indicating that oncoprotein-mediated transcriptional reprogramming is highly target cell specific.

•

RUNX1-ETO reversibly arrests the growth of human ESC-derived early myeloid cells

•

RUNX1-ETO disrupts global RUNX1 binding and deregulates RUNX1 target genes

•

RUNX1-ETO blocks myeloid differentiation by rapidly downregulating SPI1 and CEBPA

•

The impact of RUNX1-ETO induction is cell type specific

Logic and lineage impacts on functional transcription factor deployment for T-cell fate commitment

Transcription factors are the major agents that read the regulatory sequence information in the genome to initiate changes in expression of specific genes, both in development and in physiological activation responses. Their actions depend on site-specific DNA binding and are largely guided by their individual DNA target sequence specificities. However, their action is far more conditional in a real developmental context than would be expected for simple reading of local genomic DNA sequence, which is common to all cells in the organism. They are constrained by slow-changing chromatin states and by interactions with other transcription factors, which affect their occupancy patterns of potential sites across the genome. These mechanisms lead to emergent discontinuities in function even for transcription factors with minimally changing expression. This is well revealed by diverse lineages of blood cells developing throughout life from hematopoietic stem cells, which use overlapping combinations of transcription factors to drive strongly divergent gene regulation programs. Here, using development of T lymphocytes from hematopoietic multipotent progenitor cells as a focus, recent evidence is reviewed on how binding specificity and dynamics, transcription factor cooperativity, and chromatin state changes impact the effective regulatory functions of key transcription factors including PU.1, Runx1, Notch-RBPJ, and Bcl11b.

How haematopoiesis research became a fertile ground for regulatory network biology as pioneered by Eric Davidson

PURPOSE OF REVIEW

This historical perspective reviews how work of Eric H. Davidson was a catalyst and exemplar for explaining haematopoietic cell fate determination through gene regulation.

RECENT FINDINGS

Researchers studying blood and immune cells pioneered many of the early mechanistic investigations of mammalian gene regulatory processes. These efforts included the characterization of complex gene regulatory sequences exemplified by the globin and T-cell/B-cell receptor gene loci, as well as the identification of many key regulatory transcription factors through the fine mapping of chromosome translocation breakpoints in leukaemia patients. As the repertoire of known regulators expanded, assembly into gene regulatory network models became increasingly important, not only to account for the truism that regulatory genes do not function in isolation but also to devise new ways of extracting biologically meaningful insights from even more complex information. Here we explore how Eric H. Davidson's pioneering studies of gene regulatory network control in nonvertebrate model organisms have had an important and lasting impact on research into blood and immune cell development.

SUMMARY

The intellectual framework developed by Davidson continues to contribute to haematopoietic research, and his insistence on demonstrating logic and causality still challenges the frontier of research today.

Multi-scale Dynamical Modeling of T Cell Development from an Early Thymic Progenitor State to Lineage Commitment

Intrathymic development of committed progenitor (pro)-T cells from multipotent hematopoietic precursors offers an opportunity to dissect the molecular circuitry establishing cell identity in response to environmental signals. This transition encompasses programmed shutoff of stem/progenitor genes, upregulation of T cell specification genes, proliferation, and ultimately commitment. To explain these features in light of reported cis-acting chromatin effects and experimental kinetic data, we develop a three-level dynamic model of commitment based upon regulation of the commitment-linked gene Bcl11b. The levels are (1) a core gene regulatory network (GRN) architecture from transcription factor (TF) perturbation data, (2) a stochastically controlled chromatin-state gate, and (3) a single-cell proliferation model validated by experimental clonal growth and commitment kinetic assays. Using RNA fluorescence in situ hybridization (FISH) measurements of genes encoding key TFs and measured bulk population dynamics, this single-cell model predicts state-switching kinetics validated by measured clonal proliferation and commitment times. The resulting multi-scale model provides a mechanistic framework for dissecting commitment dynamics.

Regulatory sharing between estrogen receptor α bound enhancers

The human genome encodes an order of magnitude more gene expression enhancers than promoters, suggesting that most genes are regulated by the combined action of multiple enhancers. We have previously shown that neighboring estrogen-responsive enhancers exhibit complex synergistic contributions to the production of an estrogenic transcriptional response. Here we sought to determine the molecular underpinnings of this enhancer cooperativity. We generated genetic deletions of four estrogen receptor α (ER) bound enhancers that regulate two genes and found that enhancers containing full estrogen response element (ERE) motifs control ER binding at neighboring sites, while enhancers with pre-existing histone acetylation/accessibility confer a permissible chromatin environment to the neighboring enhancers. Genome engineering revealed that two enhancers with half EREs could not compensate for the lack of a full ERE site within the cluster. In contrast, two enhancers with full EREs produced a transcriptional response greater than the wild-type locus. By swapping genomic sequences, we found that the genomic location of a full ERE strongly influences enhancer activity. Our results lead to a model in which a full ERE is required for ER recruitment, but the presence of a pre-existing permissible chromatin environment can also be needed for estrogen-driven gene regulation to occur.

Interconversion between Tumorigenic and Differentiated States in Acute Myeloid Leukemia

Tumors are composed of phenotypically heterogeneous cancer cells that often resemble various differentiation states of their lineage of origin. Within this hierarchy, it is thought that an immature subpopulation of tumor-propagating cancer stem cells (CSCs) differentiates into non-tumorigenic progeny, providing a rationale for therapeutic strategies that specifically eradicate CSCs or induce their differentiation. The clinical success of these approaches depends on CSC differentiation being unidirectional rather than reversible, yet this question remains unresolved even in prototypically hierarchical malignancies, such as acute myeloid leukemia (AML). Here, we show in murine and human models of AML that, upon perturbation of endogenous expression of the lineage-determining transcription factor PU.1 or withdrawal of established differentiation therapies, some mature leukemia cells can de-differentiate and reacquire clonogenic and leukemogenic properties. Our results reveal plasticity of CSC maturation in AML, highlighting the need to therapeutically eradicate cancer cells across a range of differentiation states.

A network-based computational framework to predict and differentiate functions for gene isoforms using exon-level expression data

MOTIVATION

Alternative splicing makes significant contributions to functional diversity of transcripts and proteins. Many alternatively spliced gene isoforms have been shown to perform specific biological functions under different contexts. In addition to gene-level expression, the advances of high-throughput sequencing offer a chance to estimate isoform-specific exon expression with a high resolution, which is informative for studying splice variants with network analysis.

RESULTS

In this study, we propose a novel network-based analysis framework to predict isoform-specific functions from exon-level RNA-Seq data. In particular, based on exon-level expression data, we firstly propose a unified framework, referred to as Iso-Net, to integrate two new mathematical methods (named MINet and RVNet) that infer co-expression networks at different data scenarios. We demonstrate the superior prediction accuracy of Iso-Net over the existing methods for most simulation data, especially in two extreme cases: sample size is very small and exon numbers of two isoforms are quite different. Furthermore, by defining relevant quantitative measures (e.g., Jaccard correlation coefficient) and combining differential co-expression network analysis and GO functional enrichment analysis, a co-expression network analysis framework is developed to predict functions of isoforms and further, to discover their distinct functions within the same gene. We apply Iso-Net to study gene isoforms for several important transcription factors in human myeloid differentiation with the exon-level RNA-Seq data from three different cell lines.

AVAILABILITY AND IMPLEMENTATION

Iso-Net is open source and freely available from https://github.com/Dingjie-Wang/Iso-Net.

Labelled regulatory elements are pervasive features of the macrophage genome and are dynamically utilized by classical and alternative polarization signals

The concept of tissue-specific gene expression posits that lineage-determining transcription factors (LDTFs) determine the open chromatin profile of a cell via collaborative binding, providing molecular beacons to signal-dependent transcription factors (SDTFs). However, the guiding principles of LDTF binding, chromatin accessibility and enhancer activity have not yet been systematically evaluated. We sought to study these features of the macrophage genome by the combination of experimental (ChIP-seq, ATAC-seq and GRO-seq) and computational approaches. We show that Random Forest and Support Vector Regression machine learning methods can accurately predict chromatin accessibility using the binding patterns of the LDTF PU.1 and four other key TFs of macrophages (IRF8, JUNB, CEBPA and RUNX1). Any of these TFs alone were not sufficient to predict open chromatin, indicating that TF binding is widespread at closed or weakly opened chromatin regions. Analysis of the PU.1 cistrome revealed that two-thirds of PU.1 binding occurs at low accessible chromatin. We termed these sites labelled regulatory elements (LREs), which may represent a dormant state of a future enhancer and contribute to macrophage cellular plasticity. Collectively, our work demonstrates the existence of LREs occupied by various key TFs, regulating specific gene expression programs triggered by divergent macrophage polarizing stimuli.

Dynamic control of the T-cell specification gene regulatory network

Specification of multipotent blood precursor cells in postnatal mice to become committed T-cell precursors involves a gene regulatory network of several interacting but functionally distinct modules. Many links of this network have been defined by perturbation tests and by functional genomics. However, using the network model to predict real-life kinetics of the commitment process is still difficult, partly due to the tenacity of repressive chromatin states, and to the ability of transcription factors to affect each other's binding site choices through competitive recruitment to alternative sites ("coregulator theft"). To predict kinetics, future models will need to incorporate mechanistic information about chromatin state change dynamics and more sophisticated understanding of the proteomics and cooperative DNA site choices of transcription factor complexes.

Models of dendritic cell development correlate ontogeny with function

Rapid advances have been made to uncover the mechanisms that regulate dendritic cell (DC) development, and in turn, how models of development can be employed to define dendritic cell function. Models of DC development have been used to define the unique functions of DC subsets during immune responses to distinct pathogens. More recently, models of DC function have expanded to include their homeostatic and inflammatory physiology, modes of communication with various innate and adaptive immune lineages, and specialized functions across different lymphoid organs. New models of DC development call for revisions of previously accepted paradigms with respect to the ontogeny of plasmacytoid DC (pDC) and classical DC (cDC) subsets. By far, development of the cDC1 subset is best understood, and models have now been developed that can separate deficiencies in development from deficiencies in function. Such models are lacking for pDCs and cDC2s, limiting the depth of our understanding of their unique and essential roles during immune responses. If novel immunotherapies aim to harness the functions of human DCs, understanding of DC development will be essential to develop models DC function. Here we review emerging models of DC development and function.

Runx1 promotes murine erythroid progenitor proliferation and inhibits differentiation by preventing Pu.1 downregulation

Pu.1 is an ETS family transcription factor (TF) that plays critical roles in erythroid progenitors by promoting proliferation and blocking terminal differentiation. However, the mechanisms controlling expression and down-regulation of Pu.1 during early erythropoiesis have not been defined. In this study, we identify the actions of Runx1 and Pu.1 itself at the Pu.1 gene Upstream Regulatory Element (URE) as major regulators of Pu.1 expression in Burst-Forming Unit erythrocytes (BFUe). During early erythropoiesis, Runx1 and Pu.1 levels decline, and chromatin accessibility at the URE is lost. Ectopic expression of Runx1 or Pu.1, both of which bind the URE, prevents Pu.1 down-regulation and blocks terminal erythroid differentiation, resulting in extensive ex vivo proliferation and immortalization of erythroid progenitors. Ectopic expression of Runx1 in BFUe lacking a URE fails to block terminal erythroid differentiation. Thus, Runx1, acting at the URE, and Pu.1 itself directly regulate Pu.1 levels in erythroid cells, and loss of both factors is critical for Pu.1 down-regulation during terminal differentiation. The molecular mechanism of URE inactivation in erythroid cells through loss of TF binding represents a distinct pattern of Pu.1 regulation from those described in other hematopoietic cell types such as T cells which down-regulate Pu.1 through active repression. The importance of down-regulation of Runx1 and Pu.1 in erythropoiesis is further supported by genome-wide analyses showing that their DNA-binding motifs are highly overrepresented in regions that lose chromatin accessibility during early erythroid development.

An enriched network motif family regulates multistep cell fate transitions with restricted reversibility

Multistep cell fate transitions with stepwise changes of transcriptional profiles are common to many developmental, regenerative and pathological processes. The multiple intermediate cell lineage states can serve as differentiation checkpoints or branching points for channeling cells to more than one lineages. However, mechanisms underlying these transitions remain elusive. Here, we explored gene regulatory circuits that can generate multiple intermediate cellular states with stepwise modulations of transcription factors. With unbiased searching in the network topology space, we found a motif family containing a large set of networks can give rise to four attractors with the stepwise regulations of transcription factors, which limit the reversibility of three consecutive steps of the lineage transition. We found that there is an enrichment of these motifs in a transcriptional network controlling the early T cell development, and a mathematical model based on this network recapitulates multistep transitions in the early T cell lineage commitment. By calculating the energy landscape and minimum action paths for the T cell model, we quantified the stochastic dynamics of the critical factors in response to the differentiation signal with fluctuations. These results are in good agreement with experimental observations and they suggest the stable characteristics of the intermediate states in the T cell differentiation. These dynamical features may help to direct the cells to correct lineages during development. Our findings provide general design principles for multistep cell linage transitions and new insights into the early T cell development. The network motifs containing a large family of topologies can be useful for analyzing diverse biological systems with multistep transitions.

The functions of cells are dynamically controlled in many biological processes including development, regeneration and disease progression. Cell fate transition, or the switch of cellular functions, often involves multiple steps. The intermediate stages of the transition provide the biological systems with the opportunities to regulate the transitions in a precise manner. These transitions are controlled by key regulatory genes of which the expression shows stepwise patterns, but how the interactions of these genes can determine the multistep processes was unclear. Here, we present a comprehensive analysis on the design principles of gene circuits that govern multistep cell fate transition. We found a large network family with common structural features that can generate systems with the ability to control three consecutive steps of the transition. We found that this type of networks is enriched in a gene circuit controlling the development of T lymphocyte, a crucial type of immune cells. We performed mathematical modeling using this gene circuit and we recapitulated the stepwise and irreversible loss of stem cell properties of the developing T lymphocytes. Our findings can be useful to analyze a wide range of gene regulatory networks controlling multistep cell fate transitions.

STAT5-Driven Enhancers Tightly Control Temporal Expression of Mammary-Specific Genes

The de novo formation of milk-secreting mammary epithelium during pregnancy is regulated by prolactin through activation of the transcription factor STAT5, which stimulates the expression of several hundred mammary-specific genes. In addition to its key role in activating gene expression in mammary tissue, STAT5, which is ubiquitously expressed in most cell types, implements T cell-specific programs controlled by interleukins. However, the mechanisms by which STAT5 controls cell-specific genetic programs activated by distinct cytokines remain relatively unknown. Integration of data from genome-wide surveys of chromatin markers and transcription factor binding at regulatory elements may shed light on the mechanisms that drive cell-specific programs. Here, we have illustrated how STAT5 controls cell-specific gene expression through its concentration and an auto-regulatory enhancer supporting its high levels in mammary tissue. The unique genomic features of STAT5-driven enhancers or super-enhancers that regulate mammary-specific genes and their dynamic remodeling in response to pregnancy hormone levels are described. We have further provided biological evidence supporting the in vivo function of a STAT5-driven super-enhancer with the aid of CRISPR/Cas9 genome editing. Finally, we discuss how the functions of mammary-specific super-enhancers are confined by the zinc finger protein, CTCF, to allow exclusive activation of mammary-specific genes without affecting common neighboring genes. This review comprehensively summarizes the molecular pathways underlying differential control of cell-specific gene sets by STAT5 and provides novel insights into STAT5-dependent mammary physiology.

Mechanisms of Action of Hematopoietic Transcription Factor PU.1 in Initiation of T-Cell Development

PU.1 is an ETS-family transcription factor that plays a broad range of roles in hematopoiesis. A direct regulator of myeloid, dendritic-cell, and B cell functional programs, and a well-known antagonist of terminal erythroid cell differentiation, it is also expressed in the earliest stages of T-cell development of each cohort of intrathymic pro-T cells. Its expression in this context appears to give T-cell precursors initial, transient access to myeloid and dendritic cell developmental competence and therefore to represent a source of antagonism or delay of T-cell lineage commitment. However, it has remained uncertain until recently why T-cell development is also intensely dependent upon PU.1. Here, we review recent work that sheds light on the molecular biology of PU.1 action across the genome in pro-T cells and identifies the genes that depend on PU.1 for their correct regulation. This work indicates modes of chromatin engagement, pioneering, and cofactor recruitment (“coregulator theft”) by PU.1 as well as gene network interactions that not only affect specific target genes but also have system-wide regulatory consequences, amplifying the impact of PU.1 beyond its own direct binding targets. The genes directly regulated by PU.1 also suggest a far-reaching transformation of cell biology and signaling potential between the early stages of T-cell development when PU.1 is expressed and when it is silenced. These cell-biological functions can be important to distinguish fetal from adult T-cell development and have the potential to illuminate aspects of thymic function that have so far remained the most mysterious.

Temporal autoregulation during human PU.1 locus SubTAD formation

Epigenetic control of gene expression occurs within discrete spatial chromosomal units called topologically associating domains (TADs), but the exact spatial requirements of most genes are unknown; this is of particular interest for genes involved in cancer. We therefore applied high-resolution chromosomal conformation capture sequencing to map the three-dimensional (3D) organization of the human locus encoding the key myeloid transcription factor PU.1 in healthy monocytes and acute myeloid leukemia (AML) cells. We identified a dynamic ∼75-kb unit (SubTAD) as the genomic region in which spatial interactions between PU.1 gene regulatory elements occur during myeloid differentiation and are interrupted in AML. Within this SubTAD, proper initiation of the spatial chromosomal interactions requires PU.1 autoregulation and recruitment of the chromatin-adaptor protein LDB1 (LIM domain-binding protein 1). However, once these spatial interactions have occurred, LDB1 stabilizes them independently of PU.1 autoregulation. Thus, our data support that PU.1 autoregulates its expression in a "hit-and-run" manner by initiating stable chromosomal loops that result in a transcriptionally active chromatin architecture.

Illuminating stem cell transcription factor dynamics: long-term single-cell imaging of fluorescent protein fusions

Most single-cell approaches to date are based on destructive snapshot measurements which do not permit to correlate a current molecular state with future fate. However, to understand how cell fate choices are established by transcription factor networks (TFNs) regulating cell fates, TFN dynamics must be continuously monitored in single cells. Here we review how quantitative time-lapse imaging can contribute to understanding TFN dependent cell fate regulation at the single-cell level. We outline potentials of the technology and highlight challenges for interpreting the dynamics of fluorescent protein reporters that may interfere with endogenous TF function. We provide an outlook on how continuous observation of TF dynamics and single-cell fates may be complemented by perturbation studies and be linked to multidimensional molecular profiling in the future.

LMO2 is required for TAL1 DNA binding activity and initiation of definitive haematopoiesis at the haemangioblast stage

LMO2 is a bridging factor within a DNA binding complex and is required for definitive haematopoiesis to occur. The developmental stage of the block in haematopoietic specification is not known. We show that Lmo2^−/− mouse embryonic stem cells differentiated to Flk-1⁺ haemangioblasts, but less efficiently to haemogenic endothelium, which only produced primitive haematopoietic progenitors. Genome-wide approaches indicated that LMO2 is required at the haemangioblast stage to position the TAL1/LMO2/LDB1 complex to regulatory elements that are important for the establishment of the haematopoietic developmental program. In the absence of LMO2, the target site recognition of TAL1 is impaired. The lack of LMO2 resulted in altered gene expression levels already at the haemangioblast stage, with transcription factor genes accounting for ∼15% of affected genes. Comparison of Lmo2^−/− with Tal1^−/− Flk-1⁺ cells further showed that TAL1 was required to initiate or sustain Lmo2 expression.

Multiplex Enhancer Interference Reveals Collaborative Control of Gene Regulation by Estrogen Receptor α-Bound Enhancers

Multiple regulatory regions have the potential to regulate a single gene, yet how these elements combine to affect gene expression remains unclear. To uncover the combinatorial relationships between enhancers, we developed Enhancer-interference (Enhancer-i), a CRISPR interference-based approach that uses 2 different repressive domains, KRAB and SID, to prevent enhancer activation simultaneously at multiple regulatory regions. We applied Enhancer-i to promoter-distal estrogen receptor α binding sites (ERBS), which cluster around estradiol-responsive genes and therefore may collaborate to regulate gene expression. Targeting individual sites revealed predominant ERBS that are completely required for the transcriptional response, indicating a lack of redundancy. Simultaneous interference of different ERBS combinations identified supportive ERBS that contribute only when predominant sites are active. Using mathematical modeling, we find strong evidence for collaboration between predominant and supportive ERBS. Overall, our findings expose a complex functional hierarchy of enhancers, where multiple loci bound by the same transcription factor combine to fine-tune the expression of target genes.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

The formation of multicellular organisms from the egg to the adult stage is largely under genetic control. The activation of specific genes is governed by regulatory DNA sequences present nearby on the chromosome. Most of these sequences promote activation and are called enhancers. In this paper, we study two enhancers governing the expression of a gene involved in the formation of the posterior brain in vertebrates. One of these enhancers is involved in a positive feedback loop: it is itself activated by the protein product of the gene that it regulates. The other enhancer was thought to be only involved in the initial accumulation of the protein, necessary for the subsequent activation of the feedback loop. Here we show that the second enhancer directly cooperates with the autoregulatory enhancer to increase its accessibility and its activity. Our work uncovers a novel, long-range mode of cooperation between enhancers that restricts the domain of action of autoregulatory enhancers within embryos and might be essential to avoid their inappropriate activation.

Dynamic Gene Regulatory Networks of Human Myeloid Differentiation

The reconstruction of gene regulatory networks underlying cell differentiation from high-throughput gene expression and chromatin data remains a challenge. Here, we derive dynamic gene regulatory networks for human myeloid differentiation using a 5-day time series of RNA-seq and ATAC-seq data. We profile HL-60 promyelocytes differentiating into macrophages, neutrophils, monocytes, and monocyte-derived macrophages. We find a rapid response in the expression of key transcription factors and lineage markers that only regulate a subset of their targets at a given time, which is followed by chromatin accessibility changes that occur later along with further gene expression changes. We observe differences between promyelocyte- and monocyte-derived macrophages at both the transcriptional and chromatin landscape level, despite using the same differentiation stimulus, which suggest that the path taken by cells in the differentiation landscape defines their end cell state. More generally, our approach of combining neighboring time points and replicates to achieve greater sequencing depth can efficiently infer footprint-based regulatory networks from long series data.

Applying attractor dynamics to infer gene regulatory interactions involved in cellular differentiation

The dynamics of gene regulatory networks (GRNs) guide cellular differentiation. Determining the ways regulatory genes control expression of their targets is essential to understand and control cellular differentiation. The way a regulatory gene controls its target can be expressed as a gene regulatory function. Manual derivation of these regulatory functions is slow, error-prone and difficult to update as new information arises. Automating this process is a significant challenge and the subject of intensive effort. This work presents a novel approach to discovering biologically plausible gene regulatory interactions that control cellular differentiation. This method integrates known cell type expression data, genetic interactions, and knowledge of the effects of gene knockouts to determine likely GRN regulatory functions. We employ a genetic algorithm to search for candidate GRNs that use a set of transcription factors that control differentiation within a lineage. Nested canalyzing functions are used to constrain the search space to biologically plausible networks. The method identifies an ensemble of GRNs whose dynamics reproduce the gene expression pattern for each cell type within a particular lineage. The method's effectiveness was tested by inferring consensus GRNs for myeloid and pancreatic cell differentiation and comparing the predicted gene regulatory interactions to manually derived interactions. We identified many regulatory interactions reported in the literature and also found differences from published reports. These discrepancies suggest areas for biological studies of myeloid and pancreatic differentiation. We also performed a study that used defined synthetic networks to evaluate the accuracy of the automated search method and found that the search algorithm was able to discover the regulatory interactions in these defined networks with high accuracy. We suggest that the GRN functions derived from the methods described here can be used to fill gaps in knowledge about regulatory interactions and to offer hypotheses for experimental testing of GRNs that control differentiation and other biological processes.

Chromatin priming of genes in development: Concepts, mechanisms and consequences

During ontogeny, cells progress through multiple alternate differentiation states by activating distinct gene regulatory networks. In this review, we highlight the important role of chromatin priming in facilitating gene activation during lineage specification and in maintaining an epigenetic memory of previous gene activation. We show that chromatin priming is part of a hugely diverse repertoire of regulatory mechanisms that genes use to ensure that they are expressed at the correct time, in the correct cell type, and at the correct level, but also that they react to signals. We also emphasize how increasing our knowledge of these principles could inform our understanding of developmental failure and disease.

Hematopoiesis and T-cell specification as a model developmental system

The pathway to generate T cells from hematopoietic stem cells guides progenitors through a succession of fate choices while balancing differentiation progression against proliferation, stage to stage. Many elements of the regulatory system that controls this process are known, but the requirement for multiple, functionally distinct transcription factors needs clarification in terms of gene network architecture. Here, we compare the features of the T-cell specification system with the rule sets underlying two other influential types of gene network models: first, the combinatorial, hierarchical regulatory systems that generate the orderly, synchronized increases in complexity in most invertebrate embryos; second, the dueling 'master regulator' systems that are commonly used to explain bistability in microbial systems and in many fate choices in terminal differentiation. The T-cell specification process shares certain features with each of these prevalent models but differs from both of them in central respects. The T-cell system is highly combinatorial but also highly dose-sensitive in its use of crucial regulatory factors. The roles of these factors are not always T-lineage-specific, but they balance and modulate each other's activities long before any mutually exclusive silencing occurs. T-cell specification may provide a new hybrid model for gene networks in vertebrate developmental systems.

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

Distinct Genetic Networks Orchestrate the Emergence of Specific Waves of Fetal and Adult B-1 and B-2 Development

B cell development is often depicted as a linear process initiating in the fetus and continuing postnatally. Using a PU.1 hypomorphic mouse model, we found that B-1 and B-2 lymphopoiesis occurred in distinct fetal and adult waves differentially dependent on the Sfpi1 14 kB upstream regulatory element. The initial wave of fetal B-1 development was absent in PU.1 hypomorphic mice, while subsequent fetal and adult waves emerged. In contrast, B-2 lymphopoiesis occurred in distinct fetal and adult waves. Whole-transcriptome profiling of fetal and adult B cell progenitors supported the existence of three waves of B-1 and two waves of B-2 development and revealed that the network of transcription factors governing B lineage specification and commitment was highly divergent between B-1 and B-2 progenitors. These findings support the view that the B-1 and B-2 lineages are distinct and provide a genetic basis for layering of immune system development.

Irreversibility of T-Cell Specification: Insights from Computational Modelling of a Minimal Network Architecture

BACKGROUND/OBJECTIVES

A cascade of gene activations under the control of Notch signalling is required during T-cell specification, when T-cell precursors gradually lose the potential to undertake other fates and become fully committed to the T-cell lineage. We elucidate how the gene/protein dynamics for a core transcriptional module governs this important process by computational means.

METHODS

We first assembled existing knowledge about transcription factors known to be important for T-cell specification to form a minimal core module consisting of TCF-1, GATA-3, BCL11B, and PU.1 aiming at dynamical modeling. Model architecture was based on published experimental measurements of the effects on each factor when each of the others is perturbed. While several studies provided gene expression measurements at different stages of T-cell development, pure time series are not available, thus precluding a straightforward study of the dynamical interactions among these genes. We therefore translate stage dependent data into time series. A feed-forward motif with multiple positive feed-backs can account for the observed delay between BCL11B versus TCF-1 and GATA-3 activation by Notch signalling. With a novel computational approach, all 32 possible interactions among Notch signalling, TCF-1, and GATA-3 are explored by translating combinatorial logic expressions into differential equations for BCL11B production rate.

RESULTS

Our analysis reveals that only 3 of 32 possible configurations, where GATA-3 works as a dimer, are able to explain not only the time delay, but very importantly, also give rise to irreversibility. The winning models explain the data within the 95% confidence region and are consistent with regard to decay rates.

CONCLUSIONS

This first generation model for early T-cell specification has relatively few players. Yet it explains the gradual transition into a committed state with no return. Encoding logics in a rate equation setting allows determination of binding properties beyond what is possible in a Boolean network.

CD19 CAR immune pressure induces B-precursor acute lymphoblastic leukaemia lineage switch exposing inherent leukaemic plasticity

Adoptive immunotherapy using chimeric antigen receptor (CAR) expressing T cells targeting the CD19 B lineage receptor has demonstrated marked success in relapsed pre-B-cell acute lymphoblastic leukaemia (ALL). Persisting CAR-T cells generate sustained pressure against CD19 that may drive unique mechanisms of resistance. Pre-B ALL originates from a committed pre-B cell or an earlier progenitor, with potential to reprogram into other hematopoietic lineages. Here we report changes in lineage markers including myeloid conversion in patients following CD19 CAR therapy. Using murine ALL models we study the long-term effects of CD19 CAR-T cells and demonstrate partial or complete lineage switch as a consistent mechanism of CAR resistance depending on the underlying genetic oncogenic driver. Deletion of Pax5 or Ebf1 recapitulates lineage reprogramming occurring during CD19 CAR pressure. Our findings establish lineage switch as a mechanism of CAR resistance exposing inherent plasticity in genetic subtypes of pre-B-cell ALL.

CAR-T targeting CD19 have been successfully used in a variety of B-cell malignancies but patients may eventually relapse. Here, the authors show that CD19 CAR-T resistance in pre-B cell ALL can be due to the induction of a myeloid lineage switch through an epigenetic alterations in master regulators of B cell development.

GATA-1 Inhibits PU.1 Gene via DNA and Histone H3K9 Methylation of Its Distal Enhancer in Erythroleukemia

GATA-1 and PU.1 are two important hematopoietic transcription factors that mutually inhibit each other in progenitor cells to guide entrance into the erythroid or myeloid lineage, respectively. PU.1 controls its own expression during myelopoiesis by binding to the distal URE enhancer, whose deletion leads to acute myeloid leukemia (AML). We herein present evidence that GATA-1 binds to the PU.1 gene and inhibits its expression in human AML-erythroleukemias (EL). Furthermore, GATA-1 together with DNA methyl Transferase I (DNMT1) mediate repression of the PU.1 gene through the URE. Repression of the PU.1 gene involves both DNA methylation at the URE and its histone H3 lysine-K9 methylation and deacetylation as well as the H3K27 methylation at additional DNA elements and the promoter. The GATA-1-mediated inhibition of PU.1 gene transcription in human AML-EL mediated through the URE represents important mechanism that contributes to PU.1 downregulation and leukemogenesis that is sensitive to DNA demethylation therapy.

Dynamic Gene Regulatory Networks Drive Hematopoietic Specification and Differentiation

Metazoan development involves the successive activation and silencing of specific gene expression programs and is driven by tissue-specific transcription factors programming the chromatin landscape. To understand how this process executes an entire developmental pathway, we generated global gene expression, chromatin accessibility, histone modification, and transcription factor binding data from purified embryonic stem cell-derived cells representing six sequential stages of hematopoietic specification and differentiation. Our data reveal the nature of regulatory elements driving differential gene expression and inform how transcription factor binding impacts on promoter activity. We present a dynamic core regulatory network model for hematopoietic specification and demonstrate its utility for the design of reprogramming experiments. Functional studies motivated by our genome-wide data uncovered a stage-specific role for TEAD/YAP factors in mammalian hematopoietic specification. Our study presents a powerful resource for studying hematopoiesis and demonstrates how such data advance our understanding of mammalian development.

The analysis of novel distal Cebpa enhancers and silencers using a transcriptional model reveals the complex regulatory logic of hematopoietic lineage specification

C/EBPα plays an instructive role in the macrophage-neutrophil cell-fate decision and its expression is necessary for neutrophil development. How Cebpa itself is regulated in the myeloid lineage is not known. We decoded the cis-regulatory logic of Cebpa, and two other myeloid transcription factors, Egr1 and Egr2, using a combined experimental-computational approach. With a reporter design capable of detecting both distal enhancers and silencers, we analyzed 46 putative cis-regulatory modules (CRMs) in cells representing myeloid progenitors, and derived early macrophages or neutrophils. In addition to novel enhancers, this analysis revealed a surprisingly large number of silencers. We determined the regulatory roles of 15 potential transcriptional regulators by testing 32,768 alternative sequence-based transcriptional models against CRM activity data. This comprehensive analysis allowed us to infer the cis-regulatory logic for most of the CRMs. Silencer-mediated repression of Cebpa was found to be effected mainly by TFs expressed in non-myeloid lineages, highlighting a previously unappreciated contribution of long-distance silencing to hematopoietic lineage resolution. The repression of Cebpa by multiple factors expressed in alternative lineages suggests that hematopoietic genes are organized into densely interconnected repressive networks instead of hierarchies of mutually repressive pairs of pivotal TFs. More generally, our results demonstrate that de novo cis-regulatory dissection is feasible on a large scale with the aid of transcriptional modeling. Current address: Department of Biology, University of North Dakota, 10 Cornell Street, Stop 9019, Grand Forks, ND 58202-9019, USA.

Wellington-bootstrap: differential DNase-seq footprinting identifies cell-type determining transcription factors

Background

The analysis of differential gene expression is a fundamental tool to relate gene regulation with specific biological processes. Differential binding of transcription factors (TFs) can drive differential gene expression. While DNase-seq data can provide global snapshots of TF binding, tools for detecting differential binding from pairs of DNase-seq data sets are lacking.

Results

In order to link expression changes with changes in TF binding we introduce the concept of differential footprinting alongside a computational tool. We demonstrate that differential footprinting is associated with differential gene expression and can be used to define cell types by their specific TF occupancy patterns.

Conclusions

Our new tool, Wellington-bootstrap, will enable the detection of differential TF binding facilitating the study of gene regulatory systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2081-4) contains supplementary material, which is available to authorized users.

Epigenetic program and transcription factor circuitry of dendritic cell development

Dendritic cells (DC) are professional antigen presenting cells that develop from hematopoietic stem cells through successive steps of lineage commitment and differentiation. Multipotent progenitors (MPP) are committed to DC restricted common DC progenitors (CDP), which differentiate into specific DC subsets, classical DC (cDC) and plasmacytoid DC (pDC). To determine epigenetic states and regulatory circuitries during DC differentiation, we measured consecutive changes of genome-wide gene expression, histone modification and transcription factor occupancy during the sequel MPP-CDP-cDC/pDC. Specific histone marks in CDP reveal a DC-primed epigenetic signature, which is maintained and reinforced during DC differentiation. Epigenetic marks and transcription factor PU.1 occupancy increasingly coincide upon DC differentiation. By integrating PU.1 occupancy and gene expression we devised a transcription factor regulatory circuitry for DC commitment and subset specification. The circuitry provides the transcription factor hierarchy that drives the sequel MPP-CDP-cDC/pDC, including Irf4, Irf8, Tcf4, Spib and Stat factors. The circuitry also includes feedback loops inferred for individual or multiple factors, which stabilize distinct stages of DC development and DC subsets. In summary, here we describe the basic regulatory circuitry of transcription factors that drives DC development.

An autoregulatory enhancer controls mammary-specific STAT5 functions

Signal Transducers and Activators of Transcription (STATs) are principal transcription factors downstream of cytokine receptors. Although STAT5A is expressed in most tissues it remains to be understood why its premier, non-redundant functions are restricted to prolactin-induced mammary gland development and function. We report that the ubiquitously expressed Stat5a/b locus is subject to additional lineage-specific transcriptional control in mammary epithelium. Genome-wide surveys of epigenetic status and transcription factor occupancy uncovered a putative mammary-specific enhancer within the intergenic sequences separating the two Stat5 genes. This region exhibited several hallmarks of genomic enhancers, including DNaseI hypersensitivity, H3K27 acetylation and binding by GR, NFIB, ELF5 and MED1. Mammary-specific STAT5 binding was obtained at two canonical STAT5 binding motifs. CRISPR/Cas9-mediated genome editing was used to delete these sites in mice and determine their biological function. Mutant animals exhibited an 80% reduction of Stat5 levels in mammary epithelium and a concomitant reduction of STAT5-dependent gene expression. Transcriptome analysis identified a class of mammary-restricted genes that was particularly dependent on high STAT5 levels as a result of the intergenic enhancer. Taken together, the mammary-specific enhancer enables a positive feedback circuit that contributes to the remarkable abundance of STAT5 and, in turn, to the efficacy of STAT5-dependent mammary physiology.

Isocitrate treatment of acute anemia of inflammation in a mouse model

Acute and severe anemia of inflammation (AI) is a common complication of various clinical syndromes, including fulminant infections, critical illness with multiorgan failure, and exacerbations of autoimmune diseases. Building on recent data showing beneficial results with isocitrate treatment for chronic low-grade AI in a rat model, we used a mouse model of acute and severe AI induced by intraperitoneal heat-killed Brucella abortus to determine if isocitrate would be effective in this more stringent application. Inflamed mice treated with isocitrate developed an early but transient improvement in hemoglobin compared to solvent-treated controls, with a robust improvement on day 7, and only a trend towards improvement by day 14. Reticulocyte counts were increased in treated mice transiently, with no significant difference by day 21. Serum erythropoietin (EPO) levels were similar in treated versus control mice, indicating that isocitrate increased sensitivity to EPO. Serum and tissue iron levels showed no significant differences between the treated and control mice, ruling out improved iron availability as the cause of the increased response to endogenous EPO. Compared to the milder rat model, much higher doses of isocitrate were required for a relatively modest benefit.