VirtualCytometry: a webserver for evaluating immune cell differentiation using single-cell RNA sequencing data

MOTIVATION The immune system has diverse types of cells that are differentiated or activated via various signaling pathways and transcriptional regulation upon challenging conditions. Immunophenotyping by flow and mass cytometry are the major approaches for identifying key signaling molecules and transcription factors directing the transition between the functional states of immune cells. However, few proteins can be evaluated by flow cytometry in a single experiment, preventing researchers from obtaining a comprehensive picture of the molecular programs involved in immune cell differentiation. Recent advances in single-cell RNA sequencing (scRNA-seq) have enabled unbiased genome-wide quantification of gene expression in individual cells on a large scale, providing a new and versatile analytical pipeline for studying immune cell differentiation. RESULTS We present VirtualCytometry, a web-based computational pipeline for evaluating immune cell differentiation by exploiting cell-to-cell variation in gene expression with scRNA-seq data. Differentiating cells often show a continuous spectrum of cellular states rather than distinct populations. VirtualCytometry enables the identification of cellular subsets for different functional states of differentiation based on the expression of marker genes. Case studies have highlighted the usefulness of this subset analysis strategy for discovering signaling molecules and transcription factors for human T cell exhaustion, a state of T cell dysfunction, in tumor and mouse dendritic cells activated by pathogens. With more than 226 scRNA-seq datasets precompiled from public repositories covering diverse mouse and human immune cell types in normal and disease tissues, VirtualCytometry is a useful resource for the molecular dissection of immune cell differentiation. AVAILABILITY AND IMPLEMENTATION www.grnpedia.org/cytometry.

[1]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[2]  Matheus C. Bürger,et al.  Defining CD8+ T cells that provide the proliferative burst after PD-1 therapy , 2016, Nature.

[3]  Zhongming Zhao,et al.  scRNASeqDB: A Database for RNA-Seq Based Gene Expression Profiles in Human Single Cells , 2017, Genes.

[4]  M. Daha,et al.  Immune modulation of human dendritic cells by complement , 2007, European journal of immunology.

[5]  C. Rice,et al.  Interferon-stimulated genes: a complex web of host defenses. , 2014, Annual review of immunology.

[6]  Toshihisa Takagi,et al.  DNA Data Bank of Japan , 2016, Nucleic Acids Res..

[7]  I. Amit,et al.  Massively Parallel Single-Cell RNA-Seq for Marker-Free Decomposition of Tissues into Cell Types , 2014, Science.

[8]  R. Satija,et al.  Single-cell RNA sequencing to explore immune cell heterogeneity , 2017, Nature Reviews Immunology.

[9]  Fabian J Theis,et al.  Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells , 2015, Nature Biotechnology.

[10]  J. Marioni,et al.  Using single‐cell genomics to understand developmental processes and cell fate decisions , 2018, Molecular systems biology.

[11]  Kerstin B. Meyer,et al.  Single-cell reconstruction of the early maternal–fetal interface in humans , 2018, Nature.

[12]  Fabian J Theis,et al.  The Human Cell Atlas , 2017, bioRxiv.

[13]  John D. Storey,et al.  Statistical significance of variables driving systematic variation in high-dimensional data , 2013, Bioinform..

[14]  T. Schumacher,et al.  T Cell Dysfunction in Cancer. , 2018, Cancer cell.

[15]  Tal Shay,et al.  JingleBells: A Repository of Immune-Related Single-Cell RNA–Sequencing Datasets , 2017, The Journal of Immunology.

[16]  A. Kamphorst,et al.  CD8 T Cell Exhaustion in Chronic Infection and Cancer: Opportunities for Interventions. , 2018, Annual review of medicine.

[17]  Sean R. Davis,et al.  NCBI GEO: archive for functional genomics data sets—update , 2012, Nucleic Acids Res..

[18]  A. Regev,et al.  Scaling single-cell genomics from phenomenology to mechanism , 2017, Nature.

[19]  Rugang Zhang,et al.  SATB1 Expression Governs Epigenetic Repression of PD‐1 in Tumor‐Reactive T Cells , 2017, Immunity.

[20]  I. Mellman,et al.  Oncology meets immunology: the cancer-immunity cycle. , 2013, Immunity.

[21]  E. Clark,et al.  The role of CD40 and CD154/CD40L in dendritic cells. , 2009, Seminars in immunology.

[22]  Boxi Kang,et al.  Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing , 2017, Cell.

[23]  T. Hughes,et al.  The Human Transcription Factors , 2018, Cell.

[24]  Aysun Adan,et al.  Flow cytometry: basic principles and applications , 2017, Critical reviews in biotechnology.

[25]  Sachi Kato,et al.  SCPortalen: human and mouse single-cell centric database , 2017, Nucleic Acids Res..

[26]  G. Sen,et al.  Interferon-Induced Ifit Proteins: Their Role in Viral Pathogenesis , 2014, Journal of Virology.

[27]  Ariel S. Schwartz,et al.  An Atlas of Combinatorial Transcriptional Regulation in Mouse and Man , 2010, Cell.

[28]  G. Freeman,et al.  Selective expansion of a subset of exhausted CD8 T cells by αPD-L1 blockade , 2008, Proceedings of the National Academy of Sciences.

[29]  G. Nolan,et al.  Mass Cytometry: Single Cells, Many Features , 2016, Cell.

[30]  Thomas M. Keane,et al.  The European Nucleotide Archive in 2018 , 2018, Nucleic Acids Res..

[31]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[32]  Lai Guan Ng,et al.  Dimensionality reduction for visualizing single-cell data using UMAP , 2018, Nature Biotechnology.