Single-cell nanobiopsy enables multigenerational longitudinal transcriptomics of cancer cells.

Single-cell RNA sequencing has revolutionized our understanding of cellular heterogeneity, but routine methods require cell lysis and fail to probe the dynamic trajectories responsible for cellular state transitions, which can only be inferred. Here, we present a nanobiopsy platform that enables the injection of exogenous molecules and multigenerational longitudinal cytoplasmic sampling from a single cell and its progeny. The technique is based on scanning ion conductance microscopy (SICM) and, as a proof of concept, was applied to longitudinally profile the transcriptome of single glioblastoma (GBM) brain tumor cells in vitro over 72 hours. The GBM cells were biopsied before and after exposure to chemotherapy and radiotherapy, and our results suggest that treatment either induces or selects for more transcriptionally stable cells. We envision the nanobiopsy will contribute to transforming standard single-cell transcriptomics from a static analysis into a dynamic assay.

Perivascular niche cells sense thrombocytopenia and activate hematopoietic stem cells in an IL-1 dependent manner.

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr+ perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion.

BACKGROUND: Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored. RESULTS: We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells. CONCLUSIONS: Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Computational prediction and experimental validation of Salmonella Typhimurium SopE-mediated fine-tuning of autophagy in intestinal epithelial cells.

Macroautophagy is a ubiquitous homeostasis and health-promoting recycling process of eukaryotic cells, targeting misfolded proteins, damaged organelles and intracellular infectious agents. Some intracellular pathogens such as Salmonella enterica serovar Typhimurium hijack this process during pathogenesis. Here we investigate potential protein-protein interactions between host transcription factors and secreted effector proteins of Salmonella and their effect on host gene transcription. A systems-level analysis identified Salmonella effector proteins that had the potential to affect core autophagy gene regulation. The effect of a SPI-1 effector protein, SopE, that was predicted to interact with regulatory proteins of the autophagy process, was investigated to validate our approach. We then confirmed experimentally that SopE can directly bind to SP1, a host transcription factor, which modulates the expression of the autophagy gene MAP1LC3B. We also revealed that SopE might have a double role in the modulation of autophagy: Following initial increase of MAP1LC3B transcription triggered by Salmonella infection, subsequent decrease in MAP1LC3B transcription at 6h post-infection was SopE-dependent. SopE also played a role in modulation of the autophagy flux machinery, in particular MAP1LC3B and p62 autophagy proteins, depending on the level of autophagy already taking place. Upon typical infection of epithelial cells, the autophagic flux is increased. However, when autophagy was chemically induced prior to infection, SopE dampened the autophagic flux. The same was also observed when most of the intracellular Salmonella cells were not associated with the SCV (strain lacking sifA) regardless of the autophagy induction status before infection. We demonstrated how regulatory network analysis can be used to better characterise the impact of pathogenic effector proteins, in this case, Salmonella. This study complements previous work in which we had demonstrated that specific pathogen effectors can affect the autophagy process through direct interaction with autophagy proteins. Here we show that effector proteins can also influence the upstream regulation of the process. Such interdisciplinary studies can increase our understanding of the infection process and point out targets important in intestinal epithelial cell defense.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella.

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Free fatty-acid transport via CD36 drives ß-oxidation-mediated hematopoietic stem cell response to infection.

Acute infection is known to induce rapid expansion of hematopoietic stem cells (HSCs), but the mechanisms supporting this expansion remain incomplete. Using mouse models, we show that inducible CD36 is required for free fatty acid uptake by HSCs during acute infection, allowing the metabolic transition from glycolysis towards ß-oxidation. Mechanistically, high CD36 levels promote FFA uptake, which enables CPT1A to transport fatty acyl chains from the cytosol into the mitochondria. Without CD36-mediated FFA uptake, the HSCs are unable to enter the cell cycle, subsequently enhancing mortality in response to bacterial infection. These findings enhance our understanding of HSC metabolism in the bone marrow microenvironment, which supports the expansion of HSCs during pathogenic challenge.

Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia.

Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We, therefore, investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

Common clonal origin of conventional T cells and induced regulatory T cells in breast cancer patients.

Regulatory CD4+ T cells (Treg) prevent tumor clearance by conventional T cells (Tconv) comprising a major obstacle of cancer immune-surveillance. Hitherto, the mechanisms of Treg repertoire formation in human cancers remain largely unclear. Here, we analyze Treg clonal origin in breast cancer patients using T-Cell Receptor and single-cell transcriptome sequencing. While Treg in peripheral blood and breast tumors are clonally distinct, Tconv clones, including tumor-antigen reactive effectors (Teff), are detected in both compartments. Tumor-infiltrating CD4+ cells accumulate into distinct transcriptome clusters, including early activated Tconv, uncommitted Teff, Th1 Teff, suppressive Treg and pro-tumorigenic Treg. Trajectory analysis suggests early activated Tconv differentiation either into Th1 Teff or into suppressive and pro-tumorigenic Treg. Importantly, Tconv, activated Tconv and Treg share highly-expanded clones contributing up to 65% of intratumoral Treg. Here we show that Treg in human breast cancer may considerably stem from antigen-experienced Tconv converting into secondary induced Treg through intratumoral activation.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment.

BACKGROUND: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs. OBJECTIVES: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis. RESULTS: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state. CONCLUSIONS: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

ROS-mediated PI3K activation drives mitochondrial transfer from stromal cells to hematopoietic stem cells in response to infection.

Hematopoietic stem cells (HSCs) undergo rapid expansion in response to stress stimuli. Here we investigate the bioenergetic processes which facilitate the HSC expansion in response to infection. We find that infection by Gram-negative bacteria drives an increase in mitochondrial mass in mammalian HSCs, which results in a metabolic transition from glycolysis toward oxidative phosphorylation. The initial increase in mitochondrial mass occurs as a result of mitochondrial transfer from the bone marrow stromal cells (BMSCs) to HSCs through a reactive oxygen species (ROS)-dependent mechanism. Mechanistically, ROS-induced oxidative stress regulates the opening of connexin channels in a system mediated by phosphoinositide 3-kinase (PI3K) activation, which allows the mitochondria to transfer from BMSCs into HSCs. Moreover, mitochondria transfer from BMSCs into HSCs, in the response to bacterial infection, occurs before the HSCs activate their own transcriptional program for mitochondrial biogenesis. Our discovery demonstrates that mitochondrial transfer from the bone marrow microenvironment to HSCs is an early physiologic event in the mammalian response to acute bacterial infection and results in bioenergetic changes which underpin emergency granulopoiesis.

Acute myeloid leukemia induces protumoral p16INK4a-driven senescence in the bone marrow microenvironment.

Acute myeloid leukemia (AML) is an age-related disease that is highly dependent on the bone marrow (BM) microenvironment. With increasing age, tissues accumulate senescent cells, characterized by an irreversible arrest of cell proliferation and the secretion of a set of proinflammatory cytokines, chemokines, and growth factors, collectively known as the senescence-associated secretory phenotype (SASP). Here, we report that AML blasts induce a senescent phenotype in the stromal cells within the BM microenvironment and that the BM stromal cell senescence is driven by p16INK4a expression. The p16INK4a-expressing senescent stromal cells then feed back to promote AML blast survival and proliferation via the SASP. Importantly, selective elimination of p16INK4a+ senescent BM stromal cells in vivo improved the survival of mice with leukemia. Next, we find that the leukemia-driven senescent tumor microenvironment is caused by AML-induced NOX2-derived superoxide. Finally, using the p16-3MR mouse model, we show that by targeting NOX2 we reduced BM stromal cell senescence and consequently reduced AML proliferation. Together, these data identify leukemia-generated NOX2-derived superoxide as a driver of protumoral p16INK4a-dependent senescence in BM stromal cells. Our findings reveal the importance of a senescent microenvironment for the pathophysiology of leukemia. These data now open the door to investigate drugs that specifically target the "benign" senescent cells that surround and support AML.

Defining Cell Identity with Single-Cell Omics.

Cells are a fundamental unit of life, and the ability to study the phenotypes and behaviors of individual cells is crucial to understanding the workings of complex biological systems. Cell phenotypes (epigenomic, transcriptomic, proteomic, and metabolomic) exhibit dramatic heterogeneity between and within the different cell types and states underlying cellular functional diversity. Cell genotypes can also display heterogeneity throughout an organism, in the form of somatic genetic variation-most notably in the emergence and evolution of tumors. Recent technical advances in single-cell isolation and the development of omics approaches sensitive enough to reveal these aspects of cell identity have enabled a revolution in the study of multicellular systems. In this review, we discuss the technologies available to resolve the genomes, epigenomes, transcriptomes, proteomes, and metabolomes of single cells from a wide variety of living systems.

Proliferation Drives Aging-Related Functional Decline in a Subpopulation of the Hematopoietic Stem Cell Compartment.

Aging of the hematopoietic stem cell (HSC) compartment is characterized by lineage bias and reduced stem cell function, the molecular basis of which is largely unknown. Using single-cell transcriptomics, we identified a distinct subpopulation of old HSCs carrying a p53 signature indicative of stem cell decline alongside pro-proliferative JAK/STAT signaling. To investigate the relationship between JAK/STAT and p53 signaling, we challenged HSCs with a constitutively active form of JAK2 (V617F) and observed an expansion of the p53-positive subpopulation in old mice. Our results reveal cellular heterogeneity in the onset of HSC aging and implicate a role for JAK2V617F-driven proliferation in the p53-mediated functional decline of old HSCs.

Genome-wide in vivo screen identifies novel host regulators of metastatic colonization.

Metastasis is the leading cause of death for cancer patients. This multi-stage process requires tumour cells to survive in the circulation, extravasate at distant sites, then proliferate; it involves contributions from both the tumour cell and tumour microenvironment ('host', which includes stromal cells and the immune system). Studies suggest the early steps of the metastatic process are relatively efficient, with the post-extravasation regulation of tumour growth ('colonization') being critical in determining metastatic outcome. Here we show the results of screening 810 mutant mouse lines using an in vivo assay to identify microenvironmental regulators of metastatic colonization. We identify 23 genes that, when disrupted in mouse, modify the ability of tumour cells to establish metastatic foci, with 19 of these genes not previously demonstrated to play a role in host control of metastasis. The largest reduction in pulmonary metastasis was observed in sphingosine-1-phosphate (S1P) transporter spinster homologue 2 (Spns2)-deficient mice. We demonstrate a novel outcome of S1P-mediated regulation of lymphocyte trafficking, whereby deletion of Spns2, either globally or in a lymphatic endothelial-specific manner, creates a circulating lymphopenia and a higher percentage of effector T cells and natural killer (NK) cells present in the lung. This allows for potent tumour cell killing, and an overall decreased metastatic burden.

Resolving early mesoderm diversification through single-cell expression profiling.

In mammals, specification of the three major germ layers occurs during gastrulation, when cells ingressing through the primitive streak differentiate into the precursor cells of major organ systems. However, the molecular mechanisms underlying this process remain unclear, as numbers of gastrulating cells are very limited. In the mouse embryo at embryonic day 6.5, cells located at the junction between the extra-embryonic region and the epiblast on the posterior side of the embryo undergo an epithelial-to-mesenchymal transition and ingress through the primitive streak. Subsequently, cells migrate, either surrounding the prospective ectoderm contributing to the embryo proper, or into the extra-embryonic region to form the yolk sac, umbilical cord and placenta. Fate mapping has shown that mature tissues such as blood and heart originate from specific regions of the pre-gastrula epiblast, but the plasticity of cells within the embryo and the function of key cell-type-specific transcription factors remain unclear. Here we analyse 1,205 cells from the epiblast and nascent Flk1(+) mesoderm of gastrulating mouse embryos using single-cell RNA sequencing, representing the first transcriptome-wide in vivo view of early mesoderm formation during mammalian gastrulation. Additionally, using knockout mice, we study the function of Tal1, a key haematopoietic transcription factor, and demonstrate, contrary to previous studies performed using retrospective assays, that Tal1 knockout does not immediately bias precursor cells towards a cardiac fate.

Single-cell RNA sequencing reveals molecular and functional platelet bias of aged haematopoietic stem cells.

Aged haematopoietic stem cells (HSCs) generate more myeloid cells and fewer lymphoid cells compared with young HSCs, contributing to decreased adaptive immunity in aged individuals. However, it is not known how intrinsic changes to HSCs and shifts in the balance between biased HSC subsets each contribute to the altered lineage output. Here, by analysing HSC transcriptomes and HSC function at the single-cell level, we identify increased molecular platelet priming and functional platelet bias as the predominant age-dependent change to HSCs, including a significant increase in a previously unrecognized class of HSCs that exclusively produce platelets. Depletion of HSC platelet programming through loss of the FOG-1 transcription factor is accompanied by increased lymphoid output. Therefore, increased platelet bias may contribute to the age-associated decrease in lymphopoiesis.

Single-Cell RNA-Sequencing Reveals a Continuous Spectrum of Differentiation in Hematopoietic Cells.

The transcriptional programs that govern hematopoiesis have been investigated primarily by population-level analysis of hematopoietic stem and progenitor cells, which cannot reveal the continuous nature of the differentiation process. Here we applied single-cell RNA-sequencing to a population of hematopoietic cells in zebrafish as they undergo thrombocyte lineage commitment. By reconstructing their developmental chronology computationally, we were able to place each cell along a continuum from stem cell to mature cell, refining the traditional lineage tree. The progression of cells along this continuum is characterized by a highly coordinated transcriptional program, displaying simultaneous suppression of genes involved in cell proliferation and ribosomal biogenesis as the expression of lineage specific genes increases. Within this program, there is substantial heterogeneity in the expression of the key lineage regulators. Overall, the total number of genes expressed, as well as the total mRNA content of the cell, decreases as the cells undergo lineage commitment.

Dynamic spatio-temporal contribution of single ß5t+ cortical epithelial precursors to the thymus medulla.

Intrathymic T-cell development is critically dependent on cortical and medullary thymic epithelial cells (TECs). Both epithelial subsets originate during early thymus organogenesis from progenitor cells that express the thymoproteasome subunit ß5t, a typical feature of cortical TECs. Using in vivo lineage fate mapping, we demonstrate in mice that ß5t(+) TEC progenitors give rise to the medullary TEC compartment early in life but significantly limit their contribution once the medulla has completely formed. Lineage-tracing studies at single cell resolution demonstrate for young mice that the postnatal medulla is expanded from individual ß5t(+) cortical progenitors located at the cortico-medullary junction. These results therefore not only define a developmental window during which the expansion of medulla is efficiently enabled by progenitors resident in the thymic cortex, but also reveal the spatio-temporal dynamics that control the growth of the thymic medulla.

Integrin Based Isolation Enables Purification of Murine Lineage Committed Cardiomyocytes.

In contrast to mature cardiomyocytes which have limited regenerative capacity, pluripotent stem cells represent a promising source for the generation of new cardiomyocytes. The tendency of pluripotent stem cells to form teratomas and the heterogeneity from various differentiation stages and cardiomyocyte cell sub-types, however, are major obstacles to overcome before this type of therapy could be applied in a clinical setting. Thus, the identification of extracellular markers for specific cardiomyocyte progenitors and mature subpopulations is of particular importance. The delineation of cardiomyocyte surface marker patterns not only serves as a means to derive homogeneous cell populations by FACS, but is also an essential tool to understand cardiac development. By using single-cell expression profiling in early mouse embryonic hearts, we found that a combination of integrin alpha-1, alpha-5, alpha-6 and N-cadherin enables isolation of lineage committed murine cardiomyocytes. Additionally, we were able to separate trabecular cardiomyocytes from solid ventricular myocardium and atrial murine cells. These cells exhibit expected subtype specific phenotype confirmed by electrophysiological analysis. We show that integrin expression can be used for the isolation of living, functional and lineage-specific murine cardiomyocytes.

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

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.