Mammalian Transcription Factor Networks: Recent Advances in Interrogating Biological Complexity

Transcription factor (TF) networks are a key determinant of cell fate decisions in mammalian development and adult tissue homeostasis and are frequently corrupted in disease. However, our inability to experimentally resolve and interrogate the complexity of mammalian TF networks has hampered the progress in this field. Recent technological advances, in particular large-scale genome-wide approaches, single-cell methodologies, live-cell imaging, and genome editing, are emerging as important technologies in TF network biology. Several recent studies even suggest a need to re-evaluate established models of mammalian TF networks. Here, we provide a brief overview of current and emerging methods to define mammalian TF networks. We also discuss how these emerging technologies facilitate new ways to interrogate complex TF networks, consider the current open questions in the field, and comment on potential future directions and biomedical applications.ACW is funded by a Bloodwise Visiting Fellowship. HN is funded by the Japan Science and Technology Agency, the California Institute of Regenerative Medicine and Ludwig Foundation. BG is funded by Bloodwise, Cancer Research UK, the Wellcome Trust, the MRC, NIH-NIDDK and core funding from the Wellcome Trust to the Cambridge Stem Cell Institute

Mbd3/NuRD controls lymphoid cell fate and inhibits tumorigenesis by repressing a B cell transcriptional program

Differentiation of lineage-committed cells from multipotent progenitors requires establishment of accessible chromatin at lineage-specific transcriptional enhancers and promoters, which is mediated by pioneer transcription factors that recruit activating chromatin remodeling complexes. Here we show that the Mbd3/NuRD chromatin remodeling complex opposes this transcriptional pioneering during B cell programming of multipotent lymphoid progenitors by restricting chromatin accessibility at B cell enhancers and promoters. Mbd3/NuRD-deficient lymphoid progenitors therefore prematurely activate a B cell transcriptional program and are biased towards overproduction of pro-B cells at the expense of T cell progenitors. The striking reduction in early thymic T cell progenitors results in compensatory hyperproliferation of immature thymocytes and development of T cell lymphoma. Our results reveal that Mbd3/NuRD can regulate multilineage differentiation by constraining the activation of dormant lineage-specific enhancers and promoters. In this way, Mbd3/NuRD protects the multipotency of lymphoid progenitors, preventing B cell-programming transcription factors from prematurely enacting lineage commitment. Mbd3/NuRD therefore controls the fate of lymphoid progenitors, ensuring appropriate production of lineage-committed progeny and suppressing tumour formation.F.C. was supported by an EMBO long-term fellowship (1305-2015 and Marie Curie Actions LTFCOFUND2013/GA-2013-609409). Work in the Green lab is supported Cancer Research UK (grants refs. C1163/A12765 and C1163/A12526), Bloodwise (grant ref. 13003), and core support grants by the Wellcome Trust to the Cambridge Institute for Medical Research (100140/Z/12/Z) and Wellcome Trust-MRC Cambridge Stem Cell Institute (097922/Z/11/Z)

Single-cell RNA-sequencing reveals a distinct population of proglucagon-expressing cells specific to the mouse upper small intestine

Objectives: To identify sub-populations of intestinal preproglucagon-expressing (PPG) cells producing Glucagon-like Peptide-1, and their associated expression profiles of sensory receptors, thereby enabling the discovery of therapeutic strategies that target these cell populations for the treatment of diabetes and obesity. Methods: We performed single cell RNA sequencing of PPG-cells purified by flow cytometry from the upper small intestine of 3 GLU-Venus mice. Cells from 2 mice were sequenced at low depth, and from the third mouse at high depth. High quality sequencing data from 234 PPG-cells were used to identify clusters by tSNE analysis. qPCR was performed to compare the longitudinal and crypt/villus locations of cluster-specific genes. Immunofluorescence and mass spectrometry were used to confirm protein expression. Results: PPG-cells formed 3 major clusters: a group with typical characteristics of classical L-cells, including high expression of Gcg and Pyy (comprising 51% of all PPG-cells); a cell type overlapping with Gip-expressing K-cells (14%); and a unique cluster expressing Tph1 and Pzp that was predominantly located in proximal small intestine villi and co-produced 5-HT (35%). Expression of G-protein coupled receptors differed between clusters, suggesting the cell types are differentially regulated, and would be differentially targetable. Conclusions: Our findings support the emerging concept that many enteroendocrine cell populations are highly overlapping, with individual cells producing a range of peptides previously assigned to distinct cell types. Different receptor expression profiles across the clusters highlight potential drug targets to increase gut hormone secretion for the treatment of diabetes and obesity.Research in the FR/FMG lab is funded by a Wellcome joint investigator award (106262/Z/14/Z and 106263/Z/14/Z) and a joint MRC programme within the Metabolic Diseases Unit (MRC_MC_UU_12012/3). Single cell collection and analysis was supported through MRC Clinical Research Infrastructure funds for the Cambridge Single Cell Analysis Clinical Core Facility. Work in the Göttgens laboratory is supported by grants from the Wellcome Trust, Bloodwise, Cancer Research UK, NIDDK and core support grants by the Wellcome Trust to the Wellcome Trust-MRC Cambridge Stem Cell Institute

Neurogenin3 phosphorylation controls reprogramming efficiency of pancreatic ductal organoids into endocrine cells

β-cell replacement has been proposed as an effective treatment for some forms of diabetes, and in vitro methods for β-cell generation are being extensively explored. A potential source of β-cells comes from fate conversion of exocrine pancreatic cells into the endocrine lineage, by overexpression of three regulators of pancreatic endocrine formation and β-cell identity, Ngn3, Pdx1 and MafA. Pancreatic ductal organoid cultures have recently been developed that can be expanded indefinitely, while maintaining the potential to differentiate into the endocrine lineage. Here, using mouse pancreatic ductal organoids, we see that co-expression of Ngn3, Pdx1 and MafA are required and sufficient to generate cells that express insulin and resemble β-cells transcriptome-wide. Efficiency of β-like cell generation can be significantly enhanced by preventing phosphorylation of Ngn3 protein and further augmented by conditions promoting differentiation. Taken together, our new findings underline the potential of ductal organoid cultures as a source material for generation of β-like cells and demonstrate that post-translational regulation of reprogramming factors can be exploited to enhance β-cell generation

Defining Lineage Potential and Fate Behavior of Precursors during Pancreas Development

Pancreas development involves a coordinated process in which an early phase of cell segregation is followed by a longer phase of lineage restriction, expansion, and tissue remodeling. By combining clonal tracing and whole-mount reconstruction with proliferation kinetics and single-cell transcriptional profiling, we define the functional basis of pancreas morphogenesis. We show that the large-scale organization of mouse pancreas can be traced to the activity of self-renewing precursors positioned at the termini of growing ducts, which act collectively to drive serial rounds of stochastic ductal bifurcation balanced by termination. During this phase of branching morphogenesis, multipotent precursors become progressively fate-restricted, giving rise to self-renewing acinar-committed precursors that are conveyed with growing ducts, as well as ductal progenitors that expand the trailing ducts and give rise to delaminating endocrine cells. These findings define quantitatively how the functional behavior and lineage progression of precursor pools determine the large-scale patterning of pancreatic sub-compartments

Mapping and Functional Characterisation of a CTCF-Dependent Insulator Element at the 3′ Border of the Murine Scl Transcriptional Domain

The Scl gene encodes a transcription factor essential for haematopoietic development. Scl transcription is regulated by a panel of cis-elements spread over 55 kb with the most distal 3′ element being located downstream of the neighbouring gene Map17, which is co-regulated with Scl in haematopoietic cells. The Scl/Map17 domain is flanked upstream by the ubiquitously expressed Sil gene and downstream by a cluster of Cyp genes active in liver, but the mechanisms responsible for delineating the domain boundaries remain unclear. Here we report identification of a DNaseI hypersensitive site at the 3′ end of the Scl/Map17 domain and 45 kb downstream of the Scl transcription start site. This element is located at the boundary of active and inactive chromatin, does not function as a classical tissue-specific enhancer, binds CTCF and is both necessary and sufficient for insulator function in haematopoietic cells in vitro. Moreover, in a transgenic reporter assay, tissue-specific expression of the Scl promoter in brain was increased by incorporation of 350 bp flanking fragments from the +45 element. Our data suggests that the +45 region functions as a boundary element that separates the Scl/Map17 and Cyp transcriptional domains, and raise the possibility that this element may be useful for improving tissue-specific expression of transgenic constructs

The stem/progenitor landscape is reshaped in a mouse model of essential thrombocythemia and causes excess megakaryocyte production

Frameshift mutations in CALR (calreticulin) are associated with essential thrombocythemia (ET), but the stages at and mechanisms by which mutant CALR drives transformation remain incompletely defined. Here, we use single-cell approaches to examine the hematopoietic stem/progenitor cell landscape in a mouse model of mutant CALR-driven ET. We identify a trajectory linking hematopoietic stem cells (HSCs) with megakaryocytes and prospectively identify a previously unknown intermediate population that is overrepresented in the disease state. We also show that mutant CALR drives transformation primarily from the earliest stem cell compartment, with some contribution from megakaryocyte progenitors. Last, relative to wild-type HSCs, mutant CALR HSCs show increases in JAK-STAT signaling, the unfolded protein response, cell cycle, and a previously undescribed up-regulation of cholesterol biosynthesis. Overall, we have identified a novel megakaryocyte-biased cell population that is increased in a mouse model of ET and described transcriptomic changes linking CALR mutations to increased HSC proliferation and megakaryopoiesis

Cell-Specific Computational Modeling of the PIM Pathway in Acute Myeloid Leukemia

Personalized therapy is a major goal of modern oncology, as patient responses vary greatly even within a histologically defined cancer subtype. This is especially true in acute myeloid leukemia (AML), which exhibits striking heterogeneity in molecular segmentation. When calibrated to cell-specific data, executable network models can reveal subtle differences in signaling that help explain differences in drug response. Furthermore, they can suggest drug combinations to increase efficacy and combat acquired resistance. Here, we experimentally tested dynamic proteomic changes and phenotypic responses in diverse AML cell lines treated with pan-PIM kinase inhibitor and fms-related tyrosine kinase 3 ($\textit{FLT3}$) inhibitor as single agents and in combination. We constructed cell-specific executable models of the signaling axis, connecting genetic aberrations in $\textit{FLT3}$, tyrosine kinase 2 ($\textit{TYK2}$), platelet-derived growth factor receptor alpha ($\textit{PDGFRA}$), and fibroblast growth factor receptor 1 ($\textit{FGFR1}$) to cell proliferation and apoptosis via the PIM and PI3K kinases. The models capture key differences in signaling that later enabled them to accurately predict the unique proteomic changes and phenotypic responses of each cell line. Furthermore, using cell-specific models, we tailored combination therapies to individual cell lines and successfully validated their efficacy experimentally. Specifically, we showed that cells mildly responsive to PIM inhibition exhibited increased sensitivity in combination with PIK3CA inhibition. We also used the model to infer the origin of PIM resistance engineered through prolonged drug treatment of MOLM16 cell lines and successfully validated experimentally our prediction that this resistance can be overcome with AKT1/2 inhibition. Cancer Res; 77(4); 827–38.We would also like to thank Bloodwise for supporting BG, and the Israeli ministry of science, technology and space and Edmond J. Safra Center for Bioinformatics at Tel-Aviv University for supporting DS