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.We thank J. Downing (St. Jude Children's Research Hospital, Memphis, TN, USA) for the Runx1-ires-GFP mouse. Research in the authors' laboratory is supported by the Medical Research Council, Biotechnology and Biological Sciences Research Council, Leukaemia and Lymphoma Research, the Leukemia and Lymphoma Society, Microsoft Research and core support grants by the Wellcome Trust to the Cambridge Institute for Medical Research and Wellcome Trust - MRC Cambridge Stem Cell Institute. V.M. is supported by a Medical Research Council Studentship and Centenary Award and S.W. by a Microsoft Research PhD Scholarship.This is the accepted manuscript for a paper published in Nature Biotechnology 33, 269–276 (2015) doi:10.1038/nbt.315

GFI1 proteins regulate stem cell formation in the AGM

In vertebrates, the first haematopoietic stem cells (HSCs) with multi-lineage and long-term repopulating potential arise in the AGM (aorta-gonad-mesonephros) region. These HSCs are generated from a rare and transient subset of endothelial cells, called haemogenic endothelium (HE), through an endothelial-to-haematopoietic transition (EHT). Here, we establish the absolute requirement of the transcriptional repressors GFI1 and GFI1B (growth factor independence 1 and 1B) in this unique trans-differentiation process. We first demonstrate that Gfi1 expression specifically defines the rare population of HE that generates emerging HSCs. We further establish that in the absence of GFI1 proteins, HSCs and haematopoietic progenitor cells are not produced in the AGM, revealing the critical requirement for GFI1 proteins in intra-embryonic EHT. Finally, we demonstrate that GFI1 proteins recruit the chromatin-modifying protein LSD1, a member of the CoREST repressive complex, to epigenetically silence the endothelial program in HE and allow the emergence of blood cells.We thank the staff at the Advanced Imaging, animal facility, Molecular Biology Core facilities and Flow Cytometry of CRUK Manchester Institute for technical support and Michael Lie-A-Ling and Elli Marinopoulou for initiating the DamID-PIP bioinformatics project. We thank members of the Stem Cell Biology group, the Stem Cell Haematopoiesis groups and Martin Gering for valuable advice and critical reading of the manuscript. Work in our laboratory is supported by the Leukaemia and Lymphoma Research Foundation (LLR), Cancer Research UK (CRUK) and the Biotechnology and Biological Sciences Research Council (BBSRC). SC is the recipient of an MRC senior fellowship (MR/J009202/1).This is the author accepted manuscript. The final version is available from NPG via http://dx.doi.org/10.1038/ncb327

Epigenetic reprogramming of breast cancer cells with oocyte extracts

<p>Abstract</p> <p>Background</p> <p>Breast cancer is a disease characterised by both genetic and epigenetic alterations. Epigenetic silencing of tumour suppressor genes is an early event in breast carcinogenesis and reversion of gene silencing by epigenetic reprogramming can provide clues to the mechanisms responsible for tumour initiation and progression. In this study we apply the reprogramming capacity of oocytes to cancer cells in order to study breast oncogenesis.</p> <p>Results</p> <p>We show that breast cancer cells can be directly reprogrammed by amphibian oocyte extracts. The reprogramming effect, after six hours of treatment, in the absence of DNA replication, includes DNA demethylation and removal of repressive histone marks at the promoters of tumour suppressor genes; also, expression of the silenced genes is re-activated in response to treatment. This activity is specific to oocytes as it is not elicited by extracts from ovulated eggs, and is present at very limited levels in extracts from mouse embryonic stem cells. Epigenetic reprogramming in oocyte extracts results in reduction of cancer cell growth under anchorage independent conditions and a reduction in tumour growth in mouse xenografts.</p> <p>Conclusions</p> <p>This study presents a new method to investigate tumour reversion by epigenetic reprogramming. After testing extracts from different sources, we found that axolotl oocyte extracts possess superior reprogramming ability, which reverses epigenetic silencing of tumour suppressor genes and tumorigenicity of breast cancer cells in a mouse xenograft model. Therefore this system can be extremely valuable for dissecting the mechanisms involved in tumour suppressor gene silencing and identifying molecular activities capable of arresting tumour growth. These applications can ultimately shed light on the contribution of epigenetic alterations in breast cancer and advance the development of epigenetic therapies.</p

Conserved Regulation of p53 Network Dosage by MicroRNA–125b Occurs through Evolving miRNA–Target Gene Pairs

MicroRNAs regulate networks of genes to orchestrate cellular functions. MiR-125b, the vertebrate homologue of the Caenorhabditis elegans microRNA lin-4, has been implicated in the regulation of neural and hematopoietic stem cell homeostasis, analogous to how lin-4 regulates stem cells in C. elegans. Depending on the cell context, miR-125b has been proposed to regulate both apoptosis and proliferation. Because the p53 network is a central regulator of both apoptosis and proliferation, the dual roles of miR-125b raise the question of what genes in the p53 network might be regulated by miR-125b. By using a gain- and loss-of-function screen for miR-125b targets in humans, mice, and zebrafish and by validating these targets with the luciferase assay and a novel miRNA pull-down assay, we demonstrate that miR-125b directly represses 20 novel targets in the p53 network. These targets include both apoptosis regulators like Bak1, Igfbp3, Itch, Puma, Prkra, Tp53inp1, Tp53, Zac1, and also cell-cycle regulators like cyclin C, Cdc25c, Cdkn2c, Edn1, Ppp1ca, Sel1l, in the p53 network. We found that, although each miRNA–target pair was seldom conserved, miR-125b regulation of the p53 pathway is conserved at the network level. Our results lead us to propose that miR-125b buffers and fine-tunes p53 network activity by regulating the dose of both proliferative and apoptotic regulators, with implications for tissue stem cell homeostasis and oncogenesis

Toxicity Testing in the 21st Century: Defining New Risk Assessment Approaches Based on Perturbation of Intracellular Toxicity Pathways

The approaches to quantitatively assessing the health risks of chemical exposure have not changed appreciably in the past 50 to 80 years, the focus remaining on high-dose studies that measure adverse outcomes in homogeneous animal populations. This expensive, low-throughput approach relies on conservative extrapolations to relate animal studies to much lower-dose human exposures and is of questionable relevance to predicting risks to humans at their typical low exposures. It makes little use of a mechanistic understanding of the mode of action by which chemicals perturb biological processes in human cells and tissues. An alternative vision, proposed by the U.S. National Research Council (NRC) report Toxicity Testing in the 21st Century: A Vision and a Strategy, called for moving away from traditional high-dose animal studies to an approach based on perturbation of cellular responses using well-designed in vitro assays. Central to this vision are (a) “toxicity pathways” (the innate cellular pathways that may be perturbed by chemicals) and (b) the determination of chemical concentration ranges where those perturbations are likely to be excessive, thereby leading to adverse health effects if present for a prolonged duration in an intact organism. In this paper we briefly review the original NRC report and responses to that report over the past 3 years, and discuss how the change in testing might be achieved in the U.S. and in the European Union (EU). EU initiatives in developing alternatives to animal testing of cosmetic ingredients have run very much in parallel with the NRC report. Moving from current practice to the NRC vision would require using prototype toxicity pathways to develop case studies showing the new vision in action. In this vein, we also discuss how the proposed strategy for toxicity testing might be applied to the toxicity pathways associated with DNA damage and repair

Characterization of transcriptional networks in blood stem and progenitor cells using high-throughput single-cell gene expression analysis

Cellular decision-making is mediated by a complex interplay of external stimuli with the intracellular environment, in particular transcription factor regulatory networks. Here we have determined the expression of a network of 18 key haematopoietic transcription factors in 597 single primary blood stem and progenitor cells isolated from mouse bone marrow. We demonstrate that different stem/progenitor populations are characterized by distinctive transcription factor expression states, and through comprehensive bioinformatic analysis reveal positively and negatively correlated transcription factor pairings, including previously unrecognized relationships between Gata2, Gfi1 and Gfi1b. Validation using transcriptional and transgenic assays confirmed direct regulatory interactions consistent with a regulatory triad in immature blood stem cells, where Gata2 may function to modulate cross-inhibition between Gfi1 and Gfi1b. Single-cell expression profiling therefore identifies network states and allows reconstruction of network hierarchies involved in controlling stem cell fate choices, and provides a blueprint for studying both normal development and human disease

Virtual genome walking across the 32 Gb Ambystoma mexicanum genome; assembling gene models and intronic sequence

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n. The software pipeline is available from https://github.com/LooseLab/iterassemble

Computational Modeling of the Hematopoietic Erythroid-Myeloid Switch Reveals Insights into Cooperativity, Priming, and Irreversibility

Hematopoietic stem cell lineage choices are decided by genetic networks that are turned ON/OFF in a switch-like manner. However, prior to lineage commitment, genes are primed at low expression levels. Understanding the underlying molecular circuitry in terms of how it governs both a primed state and, at the other extreme, a committed state is of relevance not only to hematopoiesis but also to developmental systems in general. We develop a computational model for the hematopoietic erythroid-myeloid lineage decision, which is determined by a genetic switch involving the genes PU.1 and GATA-1. Dynamical models based upon known interactions between these master genes, such as mutual antagonism and autoregulation, fail to make the system bistable, a desired feature for robust lineage determination. We therefore suggest a new mechanism involving a cofactor that is regulated as well as recruited by one of the master genes to bind to the antagonistic partner that is necessary for bistability and hence switch-like behavior. An interesting fallout from this architecture is that suppression of the cofactor through external means can lead to a loss of cooperativity, and hence to a primed state for PU.1 and GATA-1. The PU. 1-GATA-1 switch also interacts with another mutually antagonistic pair, C/EBP alpha-FOG-1. The latter pair inherits the state of its upstream master genes and further reinforces the decision due to several feedback loops, thereby leading to irreversible commitment. The genetic switch, which handles the erythroid-myeloid lineage decision, is an example of a network that implements both a primed and a committed state by regulating cooperativity through recruitment of cofactors. Perturbing the feedback between the master regulators and downstream targets suggests potential reprogramming strategies. The approach points to a framework for lineage commitment studies in general and could aid the search for lineage-determining genes