xCell: Digitally portraying the tissue cellular heterogeneity landscape

Tissues are a complex milieu consisting of numerous cell types. For example, understanding the cellular heterogeneity the tumor microenvironment is an emerging field of research. Numerous methods have been published in recent years for the enumeration of cell subsets from tissue expression profiles. However, the available methods suffer from three major problems: inferring cell subset based on gene sets learned and verified from limited sources; displaying only partial portrayal of the full cellular heterogeneity; and insufficient validation in mixed tissues. To address these issues we developed xCell, a novel gene-signature based method for inferring 64 immune and stroma cell types. We first curated and harmonized 1,822 transcriptomic profiles of pure human cell types from various sources, employed a curve fitting approach for linear comparison of cell types, and introduced a novel spillover compensation technique for separating between closely related cell types. We test the ability of our model learned from pure cell types to infer enrichments of cell types in mixed tissues, using both comprehensive in silico analyses, and by comparison to cytometry immunophenotyping to show that our scores outperform previously published methods. Finally, we explore the cell type enrichments in tumor samples and show that the cellular heterogeneity of the tumor microenvironment uniquely characterizes different cancer types. We provide our method for inferring cell type abundances as a public resource to allow researchers to portray the cellular heterogeneity landscape of tissue expression profiles: http://xCell.ucsf.edu/.

[1]  Derek W Wright,et al.  Gateways to the FANTOM5 promoter level mammalian expression atlas , 2015, Genome Biology.

[2]  Charles H. Yoon,et al.  Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq , 2016, Science.

[3]  R. Weinberg,et al.  Tackling the cancer stem cells — what challenges do they pose? , 2014, Nature Reviews Drug Discovery.

[4]  Z. Modrušan,et al.  Deconvolution of Blood Microarray Data Identifies Cellular Activation Patterns in Systemic Lupus Erythematosus , 2009, PloS one.

[5]  Z. Trajanoski,et al.  Type, Density, and Location of Immune Cells Within Human Colorectal Tumors Predict Clinical Outcome , 2006, Science.

[6]  Douglas Hanahan,et al.  Accessories to the Crime: Functions of Cells Recruited to the Tumor Microenvironment Prospects and Obstacles for Therapeutic Targeting of Function-enabling Stromal Cell Types , 2022 .

[7]  Ludmila V. Danilova,et al.  Tumor immune microenvironment characterization in clear cell renal cell carcinoma identifies prognostic and immunotherapeutically relevant messenger RNA signatures , 2016, Genome Biology.

[8]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[9]  S. Ziegler FOXP3: Not just for regulatory T cells anymore , 2007, European journal of immunology.

[10]  A. Butte,et al.  Systematic pan-cancer analysis of tumour purity , 2015, Nature Communications.

[11]  Ash A. Alizadeh,et al.  Large-Scale and Comprehensive Immune Profiling and Functional Analysis of Normal Human Aging , 2015, PloS one.

[12]  Ash A. Alizadeh,et al.  Robust enumeration of cell subsets from tissue expression profiles , 2015, Nature Methods.

[13]  C. Perou,et al.  Genomic Analysis of Immune Cell Infiltrates Across 11 Tumor Types. , 2016, Journal of the National Cancer Institute.

[14]  N. Hacohen,et al.  Molecular and Genetic Properties of Tumors Associated with Local Immune Cytolytic Activity , 2015, Cell.

[15]  C B Bagwell,et al.  Fluorescence Spectral Overlap Compensation for Any Number of Flow Cytometry Parameters , 1993, Annals of the New York Academy of Sciences.

[16]  N. Friedman,et al.  Densely Interconnected Transcriptional Circuits Control Cell States in Human Hematopoiesis , 2011, Cell.

[17]  Pornpimol Charoentong,et al.  Characterization of the immunophenotypes and antigenomes of colorectal cancers reveals distinct tumor escape mechanisms and novel targets for immunotherapy , 2015, Genome Biology.

[18]  Ash A. Alizadeh,et al.  Abstract PR09: The prognostic landscape of genes and infiltrating immune cells across human cancers , 2015 .

[19]  R. Myers,et al.  Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data , 2005, Nucleic acids research.

[20]  F. Pauler,et al.  Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans , 2016, Genome Biology.

[21]  Adam A. Margolin,et al.  The Cancer Cell Line Encyclopedia enables predictive modeling of anticancer drug sensitivity , 2012, Nature.

[22]  Jun S. Liu,et al.  Comprehensive analyses of tumor immunity: implications for cancer immunotherapy , 2016, Genome Biology.

[23]  Tom C Freeman,et al.  An expression atlas of human primary cells: inference of gene function from coexpression networks , 2013, BMC Genomics.

[24]  Mary Goldman,et al.  Toil enables reproducible, open source, big biomedical data analyses , 2017, Nature Biotechnology.

[25]  S. Turley,et al.  Immunological hallmarks of stromal cells in the tumour microenvironment , 2015, Nature Reviews Immunology.

[26]  Jeffrey A. Wiser,et al.  ImmPort: disseminating data to the public for the future of immunology , 2014, Immunologic Research.

[27]  Andrew H. Beck,et al.  PharmacoGx: an R package for analysis of large pharmacogenomic datasets , 2015, Bioinform..

[28]  Z. Trajanoski,et al.  Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. , 2013, Immunity.

[29]  Pornpimol Charoentong,et al.  Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade , 2016, bioRxiv.

[30]  Andrea Bild,et al.  Alternative preprocessing of RNA-Sequencing data in The Cancer Genome Atlas leads to improved analysis results , 2015, Bioinform..

[31]  A. Butte,et al.  Widespread parainflammation in human cancer , 2016, Genome Biology.

[32]  Ash A. Alizadeh,et al.  High-throughput genomic profiling of tumor-infiltrating leukocytes. , 2016, Current opinion in immunology.

[33]  S. Shen-Orr,et al.  Computational deconvolution: extracting cell type-specific information from heterogeneous samples. , 2013, Current opinion in immunology.

[34]  Ben S. Wittner,et al.  Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1 , 2009, Nature.

[35]  Alexander R. Abbas,et al.  Immune response in silico (IRIS): immune-specific genes identified from a compendium of microarray expression data , 2005, Genes and Immunity.

[36]  ENCODEConsortium,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[37]  Peter S. Linsley,et al.  Copy Number Loss of the Interferon Gene Cluster in Melanomas Is Linked to Reduced T Cell Infiltrate and Poor Patient Prognosis , 2014, PloS one.