Evangelia Diamanti

Esrrb Is a Pivotal Target of the Gsk3/Tcf3 Axis Regulating Embryonic Stem Cell Self-Renewal

Summary Inhibition of glycogen synthase kinase-3 (Gsk3) supports mouse embryonic stem cells (ESCs) by modulating Tcf3, but the critical targets downstream of Tcf3 are unclear. We analyzed the intersection between genome localization and transcriptome data sets to identify genes repressed by Tcf3. Among these, manipulations of Esrrb gave distinctive phenotypes in functional assays. Knockdown and knockout eliminated response to Gsk3 inhibition, causing extinction of pluripotency markers and loss of colony forming capability. Conversely, forced expression phenocopied Gsk3 inhibition or Tcf3 deletion by suppressing differentiation and sustaining self-renewal. Thus the nuclear receptor Esrrb is necessary and sufficient to mediate self-renewal downstream of Gsk3 inhibition. Leukaemia inhibitory factor (LIF) regulates ESCs through Stat3, independently of Gsk3 inhibition. Consistent with parallel operation, ESCs in LIF accommodated Esrrb deletion and remained pluripotent. These findings highlight a key role for Esrrb in regulating the naive pluripotent state and illustrate compensation among the core pluripotency factors.

Decoding the regulatory network of early blood development from single-cell gene expression measurements

Reconstruction of the molecular pathways controlling organ development has been hampered by a lack of methods to resolve embryonic progenitor cells. Here we describe a strategy to address this problem that combines gene expression profiling of large numbers of single cells with data analysis based on diffusion maps for dimensionality reduction and network synthesis from state transition graphs. Applying the approach to hematopoietic development in the mouse embryo, we map the progression of mesoderm toward blood using single-cell gene expression analysis of 3,934 cells with blood-forming potential captured at four time points between E7.0 and E8.5. Transitions between individual cellular states are then used as input to develop a single-cell network synthesis toolkit to generate a computationally executable transcriptional regulatory network model of blood development. Several model predictions concerning the roles of Sox and Hox factors are validated experimentally. Our results demonstrate that single-cell analysis of a developing organ coupled with computational approaches can reveal the transcriptional programs that underpin organogenesis.

Combined Single-Cell Functional and Gene Expression Analysis Resolves Heterogeneity within Stem Cell Populations

Summary Heterogeneity within the self-renewal durability of adult hematopoietic stem cells (HSCs) challenges our understanding of the molecular framework underlying HSC function. Gene expression studies have been hampered by the presence of multiple HSC subtypes and contaminating non-HSCs in bulk HSC populations. To gain deeper insight into the gene expression program of murine HSCs, we combined single-cell functional assays with flow cytometric index sorting and single-cell gene expression assays. Through bioinformatic integration of these datasets, we designed an unbiased sorting strategy that separates non-HSCs away from HSCs, and single-cell transplantation experiments using the enriched population were combined with RNA-seq data to identify key molecules that associate with long-term durable self-renewal, producing a single-cell molecular dataset that is linked to functional stem cell activity. Finally, we demonstrated the broader applicability of this approach for linking key molecules with defined cellular functions in another stem cell system.

A tripartite transcription factor network regulates primordial germ cell specification in mice

Transitions in cell states are controlled by combinatorial actions of transcription factors. BLIMP1, the key regulator of primordial germ cell (PGC) specification, apparently acts together with PRDM14 and AP2γ. To investigate their individual and combinatorial functions, we first sought an in vitro system for transcriptional readouts and chromatin immunoprecipitation sequencing analysis. We then integrated this data with information from single-cell transcriptome analysis of normal and mutant PGCs. Here we show that BLIMP1 binds directly to repress somatic and cell proliferation genes. It also directly induces AP2γ, which together with PRDM14 initiates the PGC-specific fate. We determined the occupancy of critical genes by AP2γ—which, when computed altogether with those of BLIMP1 and PRDM14 (both individually and cooperatively), reveals a tripartite mutually interdependent transcriptional network for PGCs. We also demonstrate that, in principle, BLIMP1, AP2γ and PRDM14 are sufficient for PGC specification, and the unprecedented resetting of the epigenome towards a basal state.

A Single Cell Resolution Map of Mouse Haematopoietic Stem and Progenitor Cell Differentiation

Maintenance of the blood system requires balanced cell fate decisions by hematopoietic stem and progenitor cells (HSPCs). Because cell fate choices are executed at the individual cell level, new single-cell profiling technologies offer exciting possibilities for mapping the dynamic molecular changes underlying HSPC differentiation. Here, we have used single-cell RNA sequencing to profile more than 1600 single HSPCs, and deep sequencing has enabled detection of an average of 6558 protein-coding genes per cell. Index sorting, in combination with broad sorting gates, allowed us to retrospectively assign cells to 12 commonly sorted HSPC phenotypes while also capturing intermediate cells typically excluded by conventional gating. We further show that independently generated single-cell data sets can be projected onto the single-cell resolution expression map to directly compare data from multiple groups and to build and refine new hypotheses. Reconstruction of differentiation trajectories reveals dynamic expression changes associated with early lymphoid, erythroid, and granulocyte-macrophage differentiation. The latter two trajectories were characterized by common upregulation of cell cycle and oxidative phosphorylation transcriptional programs. By using external spike-in controls, we estimate absolute messenger RNA (mRNA) levels per cell, showing for the first time that despite a general reduction in total mRNA, a subset of genes shows higher expression levels in immature stem cells consistent with active maintenance of the stem-cell state. Finally, we report the development of an intuitive Web interface as a new community resource to permit visualization of gene expression in HSPCs at single-cell resolution for any gene of choice.

CODEX: a next-generation sequencing experiment database for the haematopoietic and embryonic stem cell communities

CODEX (http://codex.stemcells.cam.ac.uk/) is a user-friendly database for the direct access and interrogation of publicly available next-generation sequencing (NGS) data, specifically aimed at experimental biologists. In an era of multi-centre genomic dataset generation, CODEX provides a single database where these samples are collected, uniformly processed and vetted. The main drive of CODEX is to provide the wider scientific community with instant access to high-quality NGS data, which, irrespective of the publishing laboratory, is directly comparable. CODEX allows users to immediately visualize or download processed datasets, or compare user-generated data against the databases cumulative knowledge-base. CODEX contains four types of NGS experiments: transcription factor chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq), histone modification ChIP-Seq, DNase-Seq and RNA-Seq. These are largely encompassed within two specialized repositories, HAEMCODE and ESCODE, which are focused on haematopoiesis and embryonic stem cell samples, respectively. To date, CODEX contains over 1000 samples, including 221 unique TFs and 93 unique cell types. CODEX therefore provides one of the most complete resources of publicly available NGS data for the direct interrogation of transcriptional programmes that regulate cellular identity and fate in the context of mammalian development, homeostasis and disease.

Brachyury and SMAD signalling collaboratively orchestrate distinct mesoderm and endoderm gene regulatory networks in differentiating human embryonic stem cells

The transcription factor brachyury (T, BRA) is one of the first markers of gastrulation and lineage specification in vertebrates. Despite its wide use and importance in stem cell and developmental biology, its functional genomic targets in human cells are largely unknown. Here, we use differentiating human embryonic stem cells to study the role of BRA in activin A-induced endoderm and BMP4-induced mesoderm progenitors. We show that BRA has distinct genome-wide binding landscapes in these two cell populations, and that BRA interacts and collaborates with SMAD1 or SMAD2/3 signalling to regulate the expression of its target genes in a cell-specific manner. Importantly, by manipulating the levels of BRA in cells exposed to different signalling environments, we demonstrate that BRA is essential for mesoderm but not for endoderm formation. Together, our data illuminate the function of BRA in the context of human embryonic development and show that the regulatory role of BRA is context dependent. Our study reinforces the importance of analysing the functions of a transcription factor in different cellular and signalling environments. Summary: In differentiating hESCs, BRACHYURY has cell type-specific targets and functions: it acts with SMAD1 to promote mesoderm specification, while in the endoderm it interacts with SMAD2/3.

Key regulators control distinct transcriptional programmes in blood progenitor and mast cells

Despite major advances in the generation of genome‐wide binding maps, the mechanisms by which transcription factors (TFs) regulate cell type identity have remained largely obscure. Through comparative analysis of 10 key haematopoietic TFs in both mast cells and blood progenitors, we demonstrate that the largely cell type‐specific binding profiles are not opportunistic, but instead contribute to cell type‐specific transcriptional control, because (i) mathematical modelling of differential binding of shared TFs can explain differential gene expression, (ii) consensus binding sites are important for cell type‐specific binding and (iii) knock‐down of blood stem cell regulators in mast cells reveals mast cell‐specific genes as direct targets. Finally, we show that the known mast cell regulators Mitf and c‐fos likely contribute to the global reorganisation of TF binding profiles. Taken together therefore, our study elucidates how key regulatory TFs contribute to transcriptional programmes in several distinct mammalian cell types.

The human GFI136N variant induces epigenetic changes at the Hoxa9 locus and accelerates K-RAS driven myeloproliferative disorder in mice

The coding single nucleotide polymorphism GFI136N in the human gene growth factor independence 1 (GFI1) is present in 3%-7% of whites and increases the risk for acute myeloid leukemia (AML) by 60%. We show here that GFI136N, in contrast to GFI136S, lacks the ability to bind to the Gfi1 target gene that encodes the leukemia-associated transcription factor Hoxa9 and fails to initiate histone modifications that regulate HoxA9 expression. Consistent with this, AML patients heterozygous for the GFI136N variant show increased HOXA9 expression compared with normal controls. Using ChipSeq, we demonstrate that GFI136N specific epigenetic changes are also present in other genes involved in the development of AML. Moreover, granulomonocytic progenitors, a bone marrow subset from which AML can arise in humans and mice, show a proliferative expansion in the presence of the GFI136N variant. In addition, granulomonocytic progenitors carrying the GFI136N variant allele have altered gene expression patterns and differ in their ability to grow after transplantation. Finally, GFI136N can accelerate a K-RAS driven fatal myeloproliferative disease in mice. Our data suggest that the presence of a GFI136N variant allele induces a preleukemic state in myeloid precursors by deregulating the expression of Hoxa9 and other AML-related genes.

Gene set control analysis predicts hematopoietic control mechanisms from genome-wide transcription factor binding data

Transcription factors are key regulators of both normal and malignant hematopoiesis. Chromatin immunoprecipitation (ChIP) coupled with high-throughput sequencing (ChIP-Seq) has become the method of choice to interrogate the genome-wide effect of transcription factors. We have collected and integrated 142 publicly available ChIP-Seq datasets for both normal and leukemic murine blood cell types. In addition, we introduce the new bioinformatic tool Gene Set Control Analysis (GSCA). GSCA predicts likely upstream regulators for lists of genes based on statistical significance of binding event enrichment within the gene loci of a user-supplied gene set. We show that GSCA analysis of lineage-restricted gene sets reveals expected and previously unrecognized candidate upstream regulators. Moreover, application of GSCA to leukemic gene sets allowed us to predict the reactivation of blood stem cell control mechanisms as a likely contributor to LMO2 driven leukemia. It also allowed us to clarify the recent debate on the role of Myc in leukemia stem cell transcriptional programs. As a result, GSCA provides a valuable new addition to analyzing gene sets of interest, complementary to Gene Ontology and Gene Set Enrichment analyses. To facilitate access to the wider research community, we have implemented GSCA as a freely accessible web tool (http://bioinformatics.cscr.cam.ac.uk/GSCA/GSCA.html).