Into the multiverse: Advances in single-cell multiomic profiling

Single-cell transcriptomic approaches have revolutionised the study of complex biological systems, with the routine measurement of gene expression in thousands of cells enabling construction of whole-organism cell atlases. However, the transcriptome is just one layer amongst many that coordinate to define cell type and state and, ultimately, function. In parallel with the widespread uptake of single-cell RNA-seq (scRNA-seq), there has been a rapid emergence of methods that enable multiomic profiling of individual cells, enabling parallel measurement of intercellular heterogeneity in the genome, epigenome, transcriptome, and proteomes. Linking measurements from each of these layers has the potential to reveal regulatory and functional mechanisms underlying cell behaviour in healthy development and disease

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.The study was supported by Cancer Research UK grant number C45041/A14953 to A.C., C.L. and L.F and a core support grant from the Wellcome Trust and MRC to the Wellcome Trust–Medical Research Council Cambridge Stem Cell Institute. S.T would like to acknowledge the Lister Research Prize from the Lister Institute. The authors declare no competing financial interestsThis is the final version of the article. It first appeared from Cell Press via http://dx.doi.org/10.1016/j.celrep.2015.12.08

Single-cell gene and isoform expression analysis reveals signatures of ageing in haematopoietic stem and progenitor cells

Single-cell approaches have revealed that the haematopoietic hierarchy is a continuum of differentiation, from stem cell to committed progenitor, marked by changes in gene expression. However, many of these approaches neglect isoform-level information and thus do not capture the extent of alternative splicing within the system. Here, we present an integrated short- and long-read single-cell RNA-seq analysis of haematopoietic stem and progenitor cells. We demonstrate that over half of genes detected in standard short-read single-cell analyses are expressed as multiple, often functionally distinct, isoforms, including many transcription factors and key cytokine receptors. We observe global and HSC-specific changes in gene expression with ageing but limited impact of ageing on isoform usage. Integrating single-cell and cell-type-specific isoform landscape in haematopoiesis thus provides a new reference for comprehensive molecular profiling of heterogeneous tissues, as well as novel insights into transcriptional complexity, cell-type-specific splicing events and consequences of ageing

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

Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos

The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation

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.Work in the author’s laboratory is supported by grants from the Leukaemia and Lymphoma Research, the Medical Research Council, Cancer Research UK, Biotechnology and Biological Sciences Research Council, Leukemia Lymphoma Society, and the National Institute for Health Research Cambridge Biomedical Research Centre and core support grants by the Wellcome Trust to the Cambridge Institute for Medical Research and Wellcome Trust-MRC Cambridge Stem Cell Institute. D.G.K. is the recipient of a Canadian Institutes of Health Research Postdoctoral Fellowship. F.B. and F.J.T. are funded by the European Research Council (starting grant “LatentCauses”). For funding for the open access charge, the core support grant was provided by the Wellcome Trust-MRC Cambridge Stem Cell Institute. We acknowledge the support of the University of Cambridge, Cancer Research UK Institute (core grant C14303/A17197), and Hutchison Whampoa Limited.This is the final published version. It first appeared at http://www.cell.com/cell-stem-cell/abstract/S1934-5909%2815%2900162-9

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

Non-neuronal expression of SARS-CoV-2 entry genes in the olfactory system suggests mechanisms underlying COVID-19-associated anosmia

Abstract:Altered olfactory function is a common symptom of COVID-19, but its etiology is unknown. A key question is whether SARS-CoV-2 (CoV-2) – the causal agent in COVID-19 – affects olfaction directly, by infecting olfactory sensory neurons or their targets in the olfactory bulb, or indirectly, through perturbation of supporting cells. Here we identify cell types in the olfactory epithelium and olfactory bulb that express SARS-CoV-2 cell entry molecules. Bulk sequencing demonstrated that mouse, non-human primate and human olfactory mucosa expresses two key genes involved in CoV-2 entry, ACE2 and TMPRSS2. However, single cell sequencing revealed that ACE2 is expressed in support cells, stem cells, and perivascular cells, rather than in neurons. Immunostaining confirmed these results and revealed pervasive expression of ACE2 protein in dorsally-located olfactory epithelial sustentacular cells and olfactory bulb pericytes in the mouse. These findings suggest that CoV-2 infection of non-neuronal cell types leads to anosmia and related disturbances in odor perception in COVID-19 patients

Eliciting health state utilities for Dupuytren's contracture using a discrete choice experiment

Background and purpose An internet-based discrete choice experiment (DCE) was conducted to elicit preferences for a wide range of Dupuytren’s contracture (DC)-related health states. An algorithm was subsequently developed to convert these preferences into health state utilities that can be used to assess DC’s impact on quality of life and the value of its treatments. Methods Health state preferences for varying levels of DC hand severity were elicited via an internet survey from a sample of the UK adult population. Severity levels were deined using a combination of contractures (0, 45, or 90 degrees) in 8 proximal interphalangeal and metacarpophalangeal joints of the index, middle, ring, and little ingers. Right-handed, left-handed, and ambidextrous respondents indicated which hand was preferable in each of the 10 randomly-selected hand-pairings comparing different DC severity levels. For consistency across comparisons, anatomically precise digital hand drawings were used. To anchor preferences onto the traditional 0–1 utility scale used in health economic evaluations, unaffected hands were assigned a utility of 1.0 whereas the utility for a maximally affected hand (i.e., all 8 joints set at 90 degrees of contracture) was derived by asking respondents to indicate what combination of attributes and levels of the EQ-5D-5L proile most accurately relects the impact of living with such hand. Conditional logistic models were used to estimate indirect utilities, then rescaled to the anchor points on the EQ-5D-5L. Results Estimated utilities based on the responses of 1,745 qualiied respondents were 0.49, 0.57, and 0.63 for completely affected dominant hands, non-dominant hands, or ambidextrous hands, respectively. Utility for a dominant hand with 90-degree contracture in t h e metacarpophalangeal joints of the ring and little ingers was estimated to be 0.89. Separately, reducing the contracture of metacarpophalangeal joint for a little inger from 50 to 12 degrees would improve utility by 0.02. Interpretation DC is associated with substantial utility decrements. The algorithms presented herein provide a robust and lexible framework to assess utility for varying degrees of DC severity

Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch, and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here, we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify three clusters of neurons located predominantly in laminae I–III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 and ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons