The stepwise commitment from hematopoietic stem cells in the bone marrow to T lymphocyte-restricted progenitors in the thymus represents a paradigm for understanding the requirement for distinct extrinsic cues during different stages of lineage restriction from multipotent to lineage-restricted progenitors. However, the commitment stage at which progenitors migrate from the bone marrow to the thymus remains unclear. Here we provide functional and molecular evidence at the single-cell level that the earliest progenitors in the neonatal thymus had combined granulocyte-monocyte, T lymphocyte and B lymphocyte lineage potential but not megakaryocyte-erythroid lineage potential. These potentials were identical to those of candidate thymus-seeding progenitors in the bone marrow, which were closely related at the molecular level. Our findings establish the distinct lineage-restriction stage at which the T cell lineage-commitment process transits from the bone marrow to the remote thymus.

Identifying novel players of the pluripotency gene regulatory network centered on Oct4, Sox2, and Nanog as well as delineating the interactions within the complex network is key to understanding self-renewal and early cell fate commitment of embryonic stem cells (ESC). While overexpression of the transcriptional regulator Cited2 sustains ESC pluripotency, its role in ESC functions remains unclear. Here, we show that Cited2 is important for proliferation, survival, and self-renewal of mouse ESC. We position Cited2 within the pluripotency gene regulatory network by defining Nanog, Tbx3, and Klf4 as its direct targets. We also demonstrate that the defects caused by Cited2 depletion are, at least in part, rescued by Nanog constitutive expression. Finally, we demonstrate that Cited2 is required for and enhances reprogramming of mouse embryonic fibroblasts to induced pluripotent stem cells.

The ability to quantify differentiation potential of single cells is a task of critical importance. Here we demonstrate, using over 7,000 single-cell RNA-Seq profiles, that differentiation potency of a single cell can be approximated by computing the signalling promiscuity, or entropy, of a cell's transcriptome in the context of an interaction network, without the need for feature selection. We show that signalling entropy provides a more accurate and robust potency estimate than other entropy-based measures, driven in part by a subtle positive correlation between the transcriptome and connectome. Signalling entropy identifies known cell subpopulations of varying potency and drug resistant cancer stem-cell phenotypes, including those derived from circulating tumour cells. It further reveals that expression heterogeneity within single-cell populations is regulated. In summary, signalling entropy allows in silico estimation of the differentiation potency and plasticity of single cells and bulk samples, providing a means to identify normal and cancer stem-cell phenotypes.

The relationships between normal and leukemic stem/progenitor cells are unclear. We show that in ∼80% of primary human CD34+ acute myeloid leukemia (AML), two expanded populations with hemopoietic progenitor immunophenotype coexist in most patients. Both populations have leukemic stem cell (LSC) activity and are hierarchically ordered; one LSC population gives rise to the other. Global gene expression profiling shows the LSC populations are molecularly distinct and resemble normal progenitors but not stem cells. The more mature LSC population most closely mirrors normal granulocyte-macrophage progenitors (GMP) and the immature LSC population a previously uncharacterized progenitor functionally similar to lymphoid-primed multipotential progenitors (LMPPs). This suggests that in most cases primary CD34+ AML is a progenitor disease where LSCs acquire abnormal self-renewal potential.

Molecular mechanisms employed by individual multipotent cells at the point of lineage commitment remain largely uncharacterized. Current paradigms span from instructive to noise-driven mechanisms. Of considerable interest is also whether commitment involves a limited set of genes or the entire transcriptional program, and to what extent gene expression configures multiple trajectories into commitment. Importantly, the transient nature of the commitment transition confounds the experimental capture of committing cells. We develop a computational framework that simulates stochastic commitment events, and affords mechanistic exploration of the fate transition. We use a combined modeling approach guided by gene expression classifier methods that infers a time-series of stochastic commitment events from experimental growth characteristics and gene expression profiling of individual hematopoietic cells captured immediately before and after commitment. We define putative regulators of commitment and probabilistic rules of transition through machine learning methods, and employ clustering and correlation analyses to interrogate gene regulatory interactions in multipotent cells. Against this background, we develop a Monte Carlo time-series stochastic model of transcription where the parameters governing promoter status, mRNA production and mRNA decay in multipotent cells are fitted to experimental static gene expression distributions. Monte Carlo time is converted to physical time using cell culture kinetic data. Probability of commitment in time is a function of gene expression as defined by a logistic regression model obtained from experimental single-cell expression data. Our approach should be applicable to similar differentiating systems where single cell data is available. Within our system, we identify robust model solutions for the multipotent population within physiologically reasonable values and explore model predictions with regard to molecular scenarios of entry into commitment. The model suggests distinct dependencies of different commitment-associated genes on mRNA dynamics and promoter activity, which globally influence the probability of lineage commitment.

Kinase signaling fuels growth of B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Yet its role in leukemia initiation is unclear and has not been shown in primary human hematopoietic cells. We previously described activating mutations in interleukin-7 receptor alpha (IL7RA) in poor-prognosis "ph-like" BCP-ALL. Here we show that expression of activated mutant IL7RA in human CD34+ hematopoietic stem and progenitor cells induces a preleukemic state in transplanted immunodeficient NOD/LtSz-scid IL2Rγnull mice, characterized by persistence of self-renewing Pro-B cells with non-productive V(D)J gene rearrangements. Preleukemic CD34+CD10highCD19+ cells evolve into BCP-ALL with spontaneously acquired Cyclin Dependent Kinase Inhibitor 2 A (CDKN2A) deletions, as commonly observed in primary human BCP-ALL. CRISPR mediated gene silencing of CDKN2A in primary human CD34+ cells transduced with activated IL7RA results in robust development of BCP-ALLs in-vivo. Thus, we demonstrate that constitutive activation of IL7RA can initiate preleukemia in primary human hematopoietic progenitors and cooperates with CDKN2A silencing in progression into BCP-ALL.

Checkpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

ETV6-RUNX1 is almost exclusively associated with childhood B-cell acute lymphoblastic leukemia (B-ALL), but the consequences of ETV6-RUNX1 expression on cell lineage decisions during B-cell leukemogenesis are completely unknown. Clinically silent ETV6-RUNX1 preleukemic clones are frequently found in neonatal cord blood, but few carriers develop B-ALL as a result of secondary genetic alterations. The understanding of the mechanisms underlying the first transforming steps could greatly advance the development of non-toxic prophylactic interventions. Using genetic lineage tracing, we examined the capacity of ETV6-RUNX1 to instruct a malignant phenotype in the hematopoietic lineage by cell-specific Cre-mediated activation of ETV6-RUNX1 from the endogenous Etv6 gene locus. Here we show that, while ETV6-RUNX1 has the propensity to trigger both T- and B-lymphoid malignancies, it is the second hit that determines tumor cell identity. To instigate leukemia, both oncogenic hits must place early in the development of hematopoietic/precursor cells, not in already committed B-cells. Depending on the nature of the second hit, the resulting B-ALLs presented distinct entities that were clearly separable based on their gene expression profiles. Our findings give a novel mechanistic insight into the early steps of ETV6-RUNX1+ B-ALL development and might have major implications for the potential development of ETV6-RUNX1+ B-ALL prevention strategies.

Prospective isolation is critical for understanding the cellular and molecular aspects of stem cell heterogeneity. Here, we identify the cell surface antigen CD9 as a positive marker that provides a simple alternative for hematopoietic stem cell isolation at high purity. Crucially, CD9 affords the capture of all hematopoietic stem cells in murine bone marrow in the absence of contaminating populations that lack authentic stem cell function. Using CD9 as a tool to subdivide hematopoietic stem-cell-containing populations, we provide evidence for heterogeneity at the cellular, functional, and molecular levels.

Transcription factories are nuclear domains where gene transcription takes place although the molecular basis for their formation and maintenance are unknown. In this study, we explored how the properties of chromatin as a polymer may contribute to the structure of transcription factories. We found that transcriptional active chromatin contains modifications like histone H4 acetylated at Lysine 16 (H4K16ac). Single fibre analysis showed that this modification spans the entire body of the gene. Furthermore, H4K16ac genes cluster in regions up to 500 Kb alternating active and inactive chromatin. The introduction of H4K16ac in chromatin induces stiffness in the chromatin fibre. The result of this change in flexibility is that chromatin could behave like a multi-block copolymer with repetitions of stiff-flexible (active-inactive chromatin) components. Copolymers with such structure self-organize through spontaneous phase separation into microdomains. Consistent with such model H4K16ac chromatin form foci that associates with nascent transcripts. We propose that transcription factories are the result of the spontaneous concentration of H4K16ac chromatin that are in proximity, mainly in cis.

Gene loci are found in nuclear subcompartments that are related to their expression status. For instance, silent genes are often localized to heterochromatin and the nuclear periphery, whereas active genes tend to be found in the nuclear center. Evidence also suggests that chromosomes may be specifically positioned within the nucleus; however, the nature of this organization and how it is achieved are not yet fully understood. To examine whether gene regulation is related to a discernible pattern of genomic organization, we analyzed the linear arrangement of co-regulated genes along chromosomes and determined the organization of chromosomes during the differentiation of a hematopoietic progenitor to erythroid and neutrophil cell types. Our analysis reveals that there is a significant tendency for co-regulated genes to be proximal, which is related to the association of homologous chromosomes and the spatial juxtaposition of lineage-specific gene domains. We suggest that proximity in the form of chromosomal gene distribution and homolog association may be the basis for organizing the genome for coordinate gene regulation during cellular differentiation.

AML1/RUNX1 is the most frequently mutated gene in leukaemia and is central to the normal biology of hematopoietic stem and progenitor cells. However, the role of different AML1 isoforms within these primitive compartments is unclear. Here we investigate whether altering relative expression of AML1 isoforms impacts the balance between cell self-renewal and differentiation in vitro and in vivo.

Umbilical cord blood (UCB) has had considerable impact in pediatric stem cell transplantation, but its wider use is limited in part by unit size. Long-term ex vivo culture offers one approach to increase engraftment capacity by seeking to expand stem and progenitor cells. Here, we show brief incubation (8 h) of UCB CD34+ cells with the matricellular regulator Nov (CCN3) increases the frequency of serially transplantable hematopoietic stem cells (HSCs) 6-fold. This rapid response suggests recruitment rather than expansion of stem cells; accordingly, in single-cell assays, Nov increases the clonogenicity of phenotypic HSCs without increasing their number through cell division. Recruitment is associated with both metabolic and transcriptional changes, and tracing of cell divisions demonstrates that the increased clonogenic activity resides within the undivided fraction of cells. Harnessing latent stem cell potential through recruitment-based approaches will inform understanding of stem cell state transitions with implications for translation to the clinic.

Retinoids regulate a wide spectrum of cellular functions from the embryo throughout adulthood, including cell differentiation, metabolic regulation, and inflammation. These traits make retinoids very attractive molecules for medical purposes. In light of some of the physicochemical limitations of retinoids, the development of drug delivery systems offers several advantages for clinical translation of retinoid-based therapies, including improved solubilization, prolonged circulation, reduced toxicity, sustained release, and improved efficacy. In this Review, we discuss advances in preclinical and clinical tests regarding retinoid formulations, specifically the ones based in natural retinoids, evaluated in the context of regenerative medicine, brain, cancer, skin, and immune diseases. Advantages and limitations of retinoid formulations, as well as prospects to push the field forward, will be presented.

Tumor heterogeneity of human colorectal cancer (CRC)-initiating cells (CRCICs) in cancer tissues often represents aggressive features of cancer progression. For high-resolution examination of CRCICs, we performed single-cell whole-exome sequencing (scWES) and bulk cell targeted exome sequencing (TES) of CRCICs to investigate stemness-specific somatic alterations or clonal evolution.

Hematopoietic stem cells (HSCs) emerge during development from the vascular wall of the main embryonic arteries. The onset of circulation triggers several processes that provide critical external factors for HSC generation. Nevertheless, it is not fully understood how and when the onset of circulation affects HSC emergence. Here we show that in Ncx1-/- mouse embryos devoid of circulation the HSC lineage develops until the phenotypic pro-HSC stage. However, these cells reside in an abnormal microenvironment, fail to activate the hematopoietic program downstream of Runx1, and are functionally impaired. Single-cell transcriptomics shows that during the endothelial-to-hematopoietic transition, Ncx1-/- cells fail to undergo a glycolysis to oxidative phosphorylation metabolic switch present in wild-type cells. Interestingly, experimental activation of glycolysis results in decreased intraembryonic hematopoiesis. Our results suggest that the onset of circulation triggers metabolic changes that allow HSC generation to proceed.

Defective silencing of retrotransposable elements has been linked to inflammageing, cancer and autoimmune diseases. However, the underlying mechanisms are only partially understood. Here we implicate the histone H3.3 chaperone Daxx, a retrotransposable element repressor inactivated in myeloid leukaemia and other neoplasms, in protection from inflammatory disease. Loss of Daxx alters the chromatin landscape, H3.3 distribution and histone marks of haematopoietic progenitors, leading to engagement of a Pu.1-dependent transcriptional programme for myelopoiesis at the expense of B-cell differentiation. This causes neutrophilia and inflammation, predisposing mice to develop an autoinflammatory skin disease. While these molecular and phenotypic perturbations are in part reverted in animals lacking both Pu.1 and Daxx, haematopoietic progenitors in these mice show unique chromatin and transcriptome alterations, suggesting an interaction between these two pathways. Overall, our findings implicate retrotransposable element silencing in haematopoiesis and suggest a cross-talk between the H3.3 loading machinery and the pioneer transcription factor Pu.1.

Innate lymphoid cells are central to the regulation of immunity at mucosal barrier sites, with group 2 innate lymphoid cells (ILC2s) being particularly important in type 2 immunity. In this study, we demonstrate that microRNA(miR)-142 plays a critical, cell-intrinsic role in the homeostasis and function of ILC2s. Mice deficient for miR-142 expression demonstrate an ILC2 progenitor-biased development in the bone marrow, and along with peripheral ILC2s at mucosal sites, these cells display a greatly altered phenotype based on surface marker expression. ILC2 proliferative and effector functions are severely dysfunctional following Nippostrongylus brasiliensis infection, revealing a critical role for miR-142 isoforms in ILC2-mediated immune responses. Mechanistically, Socs1 and Gfi1 expression are regulated by miR-142 isoforms in ILC2s, impacting ILC2 phenotypes as well as the proliferative and effector capacity of these cells. The identification of these novel pathways opens potential new avenues to modulate ILC2-dependent immune functions.

Altered cell metabolism is inherently connected with pathological conditions including cancer and viral infections. Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiological agent of Kaposi's sarcoma (KS). KS tumour cells display features of lymphatic endothelial differentiation and in their vast majority are latently infected with KSHV, while a small number are lytically infected, producing virions. Latently infected cells express only a subset of viral genes, mainly located within the latency-associated region, among them 12 microRNAs. Notably, the metabolic properties of KSHV-infected cells closely resemble the metabolic hallmarks of cancer cells. However, how and why KSHV alters host cell metabolism remains poorly understood. Here, we investigated the effect of KSHV infection on the metabolic profile of primary dermal microvascular lymphatic endothelial cells (LEC) and the functional relevance of this effect. We found that the KSHV microRNAs within the oncogenic cluster collaborate to decrease mitochondria biogenesis and to induce aerobic glycolysis in infected cells. KSHV microRNAs expression decreases oxygen consumption, increase lactate secretion and glucose uptake, stabilize HIF1α and decreases mitochondria copy number. Importantly this metabolic shift is important for latency maintenance and provides a growth advantage. Mechanistically we show that KSHV alters host cell energy metabolism through microRNA-mediated down regulation of EGLN2 and HSPA9. Our data suggest that the KSHV microRNAs induce a metabolic transformation by concurrent regulation of two independent pathways; transcriptional reprograming via HIF1 activation and reduction of mitochondria biogenesis through down regulation of the mitochondrial import machinery. These findings implicate viral microRNAs in the regulation of the cellular metabolism and highlight new potential avenues to inhibit viral latency.

Jarid1b/KDM5b is a histone demethylase that regulates self-renewal and differentiation in stem cells and cancer; however, its function in hematopoiesis is unclear. Here, we find that Jarid1b is highly expressed in primitive hematopoietic compartments and is overexpressed in acute myeloid leukemias. Constitutive genetic deletion of Jarid1b did not impact steady-state hematopoiesis. In contrast, acute deletion of Jarid1b from bone marrow increased peripheral blood T cells and, following secondary transplantation, resulted in loss of bone marrow reconstitution. Our results reveal that deletion of Jarid1b compromises hematopoietic stem cell (HSC) self-renewal capacity and suggest that Jarid1b is a positive regulator of HSC potential.