seq-ImmuCC: Cell-Centric View of Tissue Transcriptome Measuring Cellular Compositions of Immune Microenvironment From Mouse RNA-Seq Data

The RNA sequencing approach has been broadly used to provide gene-, pathway-, and network-centric analyses for various cell and tissue samples. However, thus far, rich cellular information carried in tissue samples has not been thoroughly characterized from RNA-Seq data. Therefore, it would expand our horizons to better understand the biological processes of the body by incorporating a cell-centric view of tissue transcriptome. Here, a computational model named seq-ImmuCC was developed to infer the relative proportions of 10 major immune cells in mouse tissues from RNA-Seq data. The performance of seq-ImmuCC was evaluated among multiple computational algorithms, transcriptional platforms, and simulated and experimental datasets. The test results showed its stable performance and superb consistency with experimental observations under different conditions. With seq-ImmuCC, we generated the comprehensive landscape of immune cell compositions in 27 normal mouse tissues and extracted the distinct signatures of immune cell proportion among various tissue types. Furthermore, we quantitatively characterized and compared 18 different types of mouse tumor tissues of distinct cell origins with their immune cell compositions, which provided a comprehensive and informative measurement for the immune microenvironment inside tumor tissues. The online server of seq-ImmuCC are freely available at http://wap-lab.org:3200/immune/.

[1]  Ting Gong,et al.  DeconRNASeq: a statistical framework for deconvolution of heterogeneous tissue samples based on mRNA-Seq data , 2013, Bioinform..

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

[3]  R. Xavier,et al.  Differential Effects of Environmental and Genetic Factors on T and B Cell Immune Traits , 2016, Cell reports.

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

[5]  H. Stunnenberg,et al.  Transcriptional Landscape of Human Tissue Lymphocytes Unveils Uniqueness of Tumor-Infiltrating T Regulatory Cells , 2016, Immunity.

[6]  J. Szustakowski,et al.  Optimal Deconvolution of Transcriptional Profiling Data Using Quadratic Programming with Application to Complex Clinical Blood Samples , 2011, PloS one.

[7]  I. Amit,et al.  Digital cell quantification identifies global immune cell dynamics during influenza infection , 2014, Molecular systems biology.

[8]  R. Lahesmaa,et al.  From a gene-centric to whole-proteome view of differentiation of T helper cell subsets , 2013, Briefings in functional genomics.

[9]  M. Vetizou,et al.  Resistance Mechanisms to Immune-Checkpoint Blockade in Cancer: Tumor-Intrinsic and -Extrinsic Factors. , 2016, Immunity.

[10]  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.

[11]  Sean R. Davis,et al.  SRAdb: query and use public next-generation sequencing data from within R , 2013, BMC Bioinformatics.

[12]  Yi-jun Wang,et al.  Tumor-Infiltrating FoxP3+ Tregs and CD8+ T Cells Affect the Prognosis of Hepatocellular Carcinoma Patients , 2012, Digestion.

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

[14]  T. Vyse,et al.  Profiling RNA-Seq at multiple resolutions markedly increases the number of causal eQTLs in autoimmune disease , 2017, bioRxiv.

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

[16]  Kun Huang,et al.  MMAD: microarray microdissection with analysis of differences is a computational tool for deconvoluting cell type-specific contributions from tissue samples , 2014, Bioinform..

[17]  Trey Ideker,et al.  Genotype to phenotype via network analysis. , 2013, Current opinion in genetics & development.

[18]  Young-Mo Kim,et al.  Integrated Omics Analysis of Pathogenic Host Responses during Pandemic H1N1 Influenza Virus Infection: The Crucial Role of Lipid Metabolism. , 2016, Cell host & microbe.

[19]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

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

[21]  T. Jiang,et al.  Inference of immune cell composition on the expression profiles of mouse tissue , 2017, Scientific Reports.

[22]  D. Speiser,et al.  Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data , 2017, bioRxiv.

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

[24]  T. Chan,et al.  Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab , 2017, Cell.

[25]  M. Girvan,et al.  A pathway-centric view of spatial proximity in the 3D nucleome across cell lines , 2015, Scientific Reports.

[26]  Mei Yu,et al.  PERT: A Method for Expression Deconvolution of Human Blood Samples from Varied Microenvironmental and Developmental Conditions , 2012, PLoS Comput. Biol..