Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data.
Ontology highlight
ABSTRACT:
SUBMITTER: Racle J
PROVIDER: S-EPMC5718706 | biostudies-literature | 2017 Nov
REPOSITORIES: biostudies-literature
Publications
Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data.
Racle Julien J de Jonge Kaat K Baumgaertner Petra P Speiser Daniel E DE Gfeller David D
eLife 20171113
Immune cells infiltrating tumors can have important impact on tumor progression and response to therapy. We present an efficient algorithm to simultaneously estimate the fraction of cancer and immune cell types from bulk tumor gene expression data. Our method integrates novel gene expression profiles from each major non-malignant cell type found in tumors, renormalization based on cell-type-specific mRNA content, and the ability to consider uncharacterized and possibly highly variable cell types ...[more]
Similar Datasets
Project description:Cancers are composed of various cell types. We present an efficient algorithm to simultaneously Estimate the Proportion of Immune and Cancer cells (EPIC) from bulk tumor gene expression data. Our method integrates novel gene expression profiles from each major non-malignant cell type found in tumors, renormalization based on cell-type specific mRNA content, and the ability to consider uncharacterized and possibly highly variable cell types. Feasibility is demonstrated by validation with flow cytometry, immunohistochemistry and single-cell RNA-seq analyses of human melanoma and colorectal tumor specimens. Altogether, our work not only improves accuracy but also broadens the scope of absolute cell fraction predictions from tumor gene expression data, and provides a unique novel experimental benchmark for immunogenomics analyses in cancer research.
2017-11-15 | GSE93722 | GEO
Project description:BackgroundDifferential expression analysis is a standard approach in molecular biology. For example, genes whose expression levels differ between diseased and non-diseased samples are considered to be associated with that disease. On the other hand, differential variability analysis focuses on the differences of the variances of gene expression between sample groups. Although differential variability is also known to capture biological information, its interpretation remains unclear and controversial. Recent single-cell analyses have revealed that differences between sample groups can affect gene expression in a cellular subset-specific manner or by altering the proportion of a particular cellular subset. The aim of this study is to clarify the interpretation of mean and variance of bulk gene expression data.MethodWe developed a mathematical model in which the bulk gene expression value is proportional to the mean value of the single-cell gene expression profile. Based on this model, we performed theoretical, simulated and real single-cell RNA-seq data analyses.Result and conclusionWe identified how differences in single-cell gene expression profiles affect the differences in the mean and the variance of bulk gene expression. It is shown that differential expression analysis of bulk expression data can overlook significant changes in gene expression at the single-cell level. Further, differential variability analysis capture the complex feature affected by different gene expression shifts for each subset, changes in the proportions of cellular subsets, and variation in single-cell distribution parameters among samples.
| S-EPMC9474327 | biostudies-literature
Project description:BackgroundTranscript enumeration methods such as SAGE, MPSS, and sequencing-by-synthesis EST "digital northern", are important high-throughput techniques for digital gene expression measurement. As other counting or voting processes, these measurements constitute compositional data exhibiting properties particular to the simplex space where the summation of the components is constrained. These properties are not present on regular Euclidean spaces, on which hybridization-based microarray data is often modeled. Therefore, pattern recognition methods commonly used for microarray data analysis may be non-informative for the data generated by transcript enumeration techniques since they ignore certain fundamental properties of this space.ResultsHere we present a software tool, Simcluster, designed to perform clustering analysis for data on the simplex space. We present Simcluster as a stand-alone command-line C package and as a user-friendly on-line tool. Both versions are available at: http://xerad.systemsbiology.net/simcluster.ConclusionSimcluster is designed in accordance with a well-established mathematical framework for compositional data analysis, which provides principled procedures for dealing with the simplex space, and is thus applicable in a number of contexts, including enumeration-based gene expression data.
| S-EPMC2147035 | biostudies-literature
Project description:Computational deconvolution is a time and cost-efficient approach to obtain cell type-specific information from bulk gene expression of heterogeneous tissues like blood. Deconvolution can aim to either estimate cell type proportions or abundances in samples, or estimate how strongly each present cell type expresses different genes, or both tasks simultaneously. Among the two separate goals, the estimation of cell type proportions/abundances is widely studied, but less attention has been paid on defining the cell type-specific expression profiles. Here, we address this gap by introducing a novel method Rodeo and empirically evaluating it and the other available tools from multiple perspectives utilizing diverse datasets.
| S-EPMC7803005 | biostudies-literature
Project description:BackgroundNew technologies are focusing on characterizing cell types to better understand their heterogeneity. With large volumes of cellular data being generated, innovative methods are needed to structure the resulting data analyses. Here, we describe an 'Ontologically BAsed Molecular Signature' (OBAMS) method that identifies novel cellular biomarkers and infers biological functions as characteristics of particular cell types. This method finds molecular signatures for immune cell types based on mapping biological samples to the Cell Ontology (CL) and navigating the space of all possible pairwise comparisons between cell types to find genes whose expression is core to a particular cell type's identity.ResultsWe illustrate this ontological approach by evaluating expression data available from the Immunological Genome project (IGP) to identify unique biomarkers of mature B cell subtypes. We find that using OBAMS, candidate biomarkers can be identified at every strata of cellular identity from broad classifications to very granular. Furthermore, we show that Gene Ontology can be used to cluster cell types by shared biological processes in order to find candidate genes responsible for somatic hypermutation in germinal center B cells. Moreover, through in silico experiments based on this approach, we have identified genes sets that represent genes overexpressed in germinal center B cells and identify genes uniquely expressed in these B cells compared to other B cell types.ConclusionsThis work demonstrates the utility of incorporating structured ontological knowledge into biological data analysis - providing a new method for defining novel biomarkers and providing an opportunity for new biological insights.
| S-EPMC3844401 | biostudies-literature
Project description:Single-cell RNA sequencing (scRNA-seq) facilitates the study of cell type heterogeneity and the construction of cell atlas. However, due to its limitations, many genes may be detected to have zero expressions, i.e. dropout events, leading to bias in downstream analyses and hindering the identification and characterization of cell types and cell functions. Although many imputation methods have been developed, their performances are generally lower than expected across different kinds and dimensions of data and application scenarios. Therefore, developing an accurate and robust single-cell gene expression data imputation method is still essential. Considering to maintain the original cell-cell and gene-gene correlations and leverage bulk RNA sequencing (bulk RNA-seq) data information, we propose scINRB, a single-cell gene expression imputation method with network regularization and bulk RNA-seq data. scINRB adopts network-regularized non-negative matrix factorization to ensure that the imputed data maintains the cell-cell and gene-gene similarities and also approaches the gene average expression calculated from bulk RNA-seq data. To evaluate the performance, we test scINRB on simulated and experimental datasets and compare it with other commonly used imputation methods. The results show that scINRB recovers gene expression accurately even in the case of high dropout rates and dimensions, preserves cell-cell and gene-gene similarities and improves various downstream analyses including visualization, clustering and trajectory inference.
| S-EPMC11006796 | biostudies-literature
Project description:BACKGROUND:Clinical benefit from checkpoint inhibitors has been associated in a tumor-agnostic manner with two main tumor traits. The first is tumor antigenicity, which is typically measured by tumor mutation burden, microsatellite instability (MSI), or Mismatch Repair Deficiency using gene sequence platforms and/or immunohistochemistry. The second is the presence of a pre-existing adaptive immune response, typically measured by immunohistochemistry (e.g. single analyte PD-L1 expression) and/or gene expression signatures (e.g. tumor "inflamed" phenotype). These two traits have been shown to provide independent predictive information. Here we investigated the potential of using gene expression to predict tumor MSI, thus enabling the measurement of both tumor antigenicity and the level of tumor inflammation in a single assay, possibly reducing sample requirement, turn-around time, and overall cost. METHODS:Using The Cancer Genome Atlas RNA-seq datasets with the greatest MSI-H incidence, i.e. those from colon (n?=?208), stomach (n?=?269), and endometrial (n?=?241) cancers, we trained an algorithm to predict tumor MSI from under-expression of the mismatch repair genes MLH1, PMS2, MSH2, and MSH6 and from 10 additional genes with strong pan-cancer associations with tumor hypermutation. The algorithms were validated on the NanoString nCounter™ platform in independent cohorts of colorectal (n?=?52), endometrial (n?=?11), and neuroendocrine (n?=?4) tumors pre-characterized using the MMR immunohistochemistry assay. RESULTS:In the validation cohorts, the algorithm showed high prediction accuracy of tumor MSI status, with sensitivity of at least 88% attained at thresholds chosen to achieve 100% specificity. Furthermore, MSI status was compared to the Tumor Inflammation Signature (TIS), an analytically validated diagnostic assay which measures a suppressed adaptive immune response in the tumor and enriches for response to immune checkpoint blockade. TIS score was largely independent of MSI status, suggesting that measuring both parameters may identify more patients that would respond to immune checkpoint blockade than either assay alone. CONCLUSIONS:Development of a gene expression signature of MSI status raises the possibility of a combined diagnostic assay on a single platform which measures both tumor antigenicity and presence of a suppressed adaptive immune response. Such an assay would have significant advantages over multi-platform assays for both ease of use and turnaround time and could lead to a diagnostic test with improved clinical performance.
| S-EPMC6341623 | biostudies-literature
Project description:Immune infiltration of the tumor microenvironment has been associated with improved survival for some patients with solid tumors. The precise makeup and prognostic relevance of immune infiltrates across a broad spectrum of tumors remain unclear.Using mRNA sequencing data from The Cancer Genome Atlas (TCGA) from 11 tumor types representing 3485 tumors, we evaluated lymphocyte and macrophage gene expression by tissue type and by genomic subtypes defined within and across tumor tissue of origin (Cox proportional hazards, Pearson correlation). We investigated clonal diversity of B-cell infiltrates through calculating B-cell receptor (BCR) repertoire sequence diversity. All statistical tests were two-sided.High expression of T-cell and B-cell signatures predicted improved overall survival across many tumor types including breast, lung, and melanoma (breast CD8_T_Cells hazard ratio [HR] = 0.36, 95% confidence interval [CI] = 0.16 to 0.81, P = .01; lung adenocarcinoma B_Cell_60gene HR?=?0.71, 95% CI?=?0.58 to 0.87, P = 7.80E-04; melanoma LCK HR?=?0.86, 95% CI?=?0.79 to 0.94, P = 6.75E-04). Macrophage signatures predicted worse survival in GBM, as did B-cell signatures in renal tumors (Glioblastoma Multiforme [GBM]: macrophages HR?=?1.62, 95% CI?=?1.17 to 2.26, P = .004; renal: B_Cell_60gene HR?=?1.17, 95% CI?=?1.04 to 1.32, P = .009). BCR diversity was associated with survival beyond gene segment expression in melanoma (HR?=?2.67, 95% CI?=?1.32 to 5.40, P = .02) and renal cell carcinoma (HR?=?0.36, 95% CI?=?0.15 to 0.87, P = .006).These data support existing studies suggesting that in diverse tissue types, heterogeneous immune infiltrates are present and typically portend an improved prognosis. In some tumor types, BCR diversity was also associated with survival. Quantitative genomic signatures of immune cells warrant further testing as prognostic markers and potential biomarkers of response to cancer immunotherapy.
| S-EPMC5241901 | biostudies-other
Project description:ArrayExpress (https://www.ebi.ac.uk/arrayexpress) is an archive of functional genomics data from a variety of technologies assaying functional modalities of a genome, such as gene expression or promoter occupancy. The number of experiments based on sequencing technologies, in particular RNA-seq experiments, has been increasing over the last few years and submissions of sequencing data have overtaken microarray experiments in the last 12 months. Additionally, there is a significant increase in experiments investigating single cells, rather than bulk samples, known as single-cell RNA-seq. To accommodate these trends, we have substantially changed our submission tool Annotare which, along with raw and processed data, collects all metadata necessary to interpret these experiments. Selected datasets are re-processed and loaded into our sister resource, the value-added Expression Atlas (and its component Single Cell Expression Atlas), which not only enables users to interpret the data easily but also serves as a test for data quality. With an increasing number of studies that combine different assay modalities (multi-omics experiments), a new more general archival resource the BioStudies Database has been developed, which will eventually supersede ArrayExpress. Data submissions will continue unchanged; all existing ArrayExpress data will be incorporated into BioStudies and the existing accession numbers and application programming interfaces will be maintained.
| S-EPMC6323929 | biostudies-literature
Project description:BackgroundTo facilitate the investigation of the pathogenic roles played by various immune cells in complex tissues such as tumors, a few computational methods for deconvoluting bulk gene expression profiles to predict cell composition have been created. However, available methods were usually developed along with a set of reference gene expression profiles consisting of imbalanced replicates across different cell types. Therefore, the objective of this study was to create a new deconvolution method equipped with a new set of reference gene expression profiles that incorporate more microarray replicates of the immune cells that have been frequently implicated in the poor prognosis of cancers, such as T helper cells, regulatory T cells and macrophage M1/M2 cells.MethodsOur deconvolution method was developed by choosing ε-support vector regression (ε-SVR) as the core algorithm assigned with a loss function subject to the L1-norm penalty. To construct the reference gene expression signature matrix for regression, a subset of differentially expressed genes were chosen from 148 microarray-based gene expression profiles for 9 types of immune cells by using ANOVA and minimizing condition number. Agreement analyses including mean absolute percentage errors and Bland-Altman plots were carried out to compare the performances of our method and CIBERSORT.ResultsIn silico cell mixtures, simulated bulk tissues, and real human samples with known immune-cell fractions were used as the test datasets for benchmarking. Our method outperformed CIBERSORT in the benchmarks using in silico breast tissue-immune cell mixtures in the proportions of 30:70 and 50:50, and in the benchmark using 164 human PBMC samples. Our results suggest that the performance of our method was at least comparable to that of a state-of-the-art tool, CIBERSORT.ConclusionsWe developed a new cell composition deconvolution method and the implementation was entirely based on the publicly available R and Python packages. In addition, we compiled a new set of reference gene expression profiles, which might allow for a more robust prediction of the immune cell fractions from the expression profiles of cell mixtures. The source code of our method could be downloaded from https://github.com/holiday01/deconvolution-to-estimate-immune-cell-subsets.
| S-EPMC6923925 | biostudies-literature