scholarly journals Comprehensive analysis of immune cell enrichment in the tumor microenvironment of head and neck squamous cell carcinoma

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Ikko Mito ◽  
Hideyuki Takahashi ◽  
Reika Kawabata-Iwakawa ◽  
Shota Ida ◽  
Hiroe Tada ◽  
...  

AbstractHead and neck squamous carcinoma (HNSCC) is highly infiltrated by immune cells, including tumor-infiltrating lymphocytes and myeloid lineage cells. In the tumor microenvironment, tumor cells orchestrate a highly immunosuppressive microenvironment by secreting immunosuppressive mediators, expressing immune checkpoint ligands, and downregulating human leukocyte antigen expression. In the present study, we aimed to comprehensively profile the immune microenvironment of HNSCC using gene expression data obtained from public database. We calculated enrichment scores of 33 immune cell types based on gene expression data of HNSCC tissues and adjacent non-cancer tissues. Based on these scores, we performed non-supervised clustering and identified three immune signatures—cold, lymphocyte, and myeloid/dendritic cell (DC)—based on the clustering results. We then compared the clinical and biological features of the three signatures. Among HNSCC and non-cancer tissues, human papillomavirus (HPV)-positive HNSCCs exhibited the highest scores in various immune cell types, including CD4+ T cells, CD8+ T cells, B cells, plasma cells, basophils, and their subpopulations. Among the three immune signatures, the proportions of HPV-positive tumors, oropharyngeal cancers, early T tumors, and N factor positive cases were significantly higher in the lymphocyte signature than in other signatures. Among the three signatures, the lymphocyte signature showed the longest overall survival (OS), especially in HPV-positive patients, whereas the myeloid/DC signature demonstrated the shortest OS in these patients. Gene set enrichment analysis revealed the upregulation of several pathways related to inflammatory and proinflammatory responses in the lymphocyte signature. The expression of PRF1, IFNG, GZMB, CXCL9, CXCL10, PDCD1, LAG3, CTLA4, HAVCR2, and TIGIT was the highest in the lymphocyte signature. Meanwhile, the expression of PD-1 ligand genes CD274 and PDCD1LG2 was highest in the myeloid/DC signature. Herein, our findings revealed the transcriptomic landscape of the immune microenvironment that closely reflects the clinical and biological significance of HNSCC, indicating that molecular profiling of the immune microenvironment can be employed to develop novel biomarkers and precision immunotherapies for HNSCC.

2021 ◽  
Author(s):  
Ikko Mito ◽  
Hideyuki Takahashi ◽  
Reika Kawabata-Iwakawa ◽  
Shota Ida ◽  
Hiroe Tada ◽  
...  

Abstract Background: Head and neck squamous carcinoma (HNSCC) is highly infiltrated by immune cells, including tumor-infiltrating lymphocytes and myeloid lineage cells. In the tumor microenvironment, tumor cells orchestrate a highly immunosuppressive microenvironment by secreting immunosuppressive mediators, expressing immune checkpoint ligands, and downregulating human leukocyte antigen expression. In the present study, we aimed to comprehensively profile the immune microenvironment of HNSCC using RNA-sequencing (RNA-seq) data obtained from The Cancer Genome Atlas (TCGA) database.Methods: We calculated enrichment scores of 33 immune cell types based on RNA-seq data of HNSCC tissues and adjacent non-cancer tissues. Based on these scores, we performed non-supervised clustering and identified three immune signatures, i.e., cold, lymphocyte, and myeloid/dendritic cell (DC), using clustering results. We then compared the clinical and biological features of the three signatures.Results: Among HNSCC and non-cancer tissues, human papillomavirus (HPV)-positive HNSCCs exhibited the highest scores in various immune cell types, including CD4+ T cells, CD8+ T cells, B cells, plasma cells, basophils, and their subpopulations. Among the three immune signatures, the proportions of HPV-positive tumors, oropharyngeal cancers, early T tumors, and N factor positive cases were significantly higher in the lymphocyte signature than in other signatures. Among the three signatures, the lymphocyte signature showed the longest overall survival (OS), especially in HPV-positive patients, whereas the myeloid/DC signature demonstrated the shortest OS in these patients. Gene set enrichment analysis revealed the upregulation of several pathways related to inflammatory and proinflammatory responses in the lymphocyte signature. The expression of PRF1, IFNG, GZMB, PDCD1, LAG3, CTLA4, HAVCR2, and TIGIT was the highest in the lymphocyte signature. Meanwhile, the expression of PD-1 ligand genes CD274 and PDCD1LG2 was highest in the myeloid/DC signature. Conclusions: Herein, our findings revealed the transcriptomic landscape of the immune microenvironment that closely reflects the clinical and biological significance of HNSCC, indicating that molecular profiling of the immune microenvironment can be employed to develop novel biomarkers and precision immunotherapies for HNSCC.


2021 ◽  
Author(s):  
Ikko Mito ◽  
Hideyuki Takahashi ◽  
Reika Kawabata-Iwakawa ◽  
Shota Ida ◽  
Hiroe Tada ◽  
...  

Abstract Background: Head and neck squamous carcinoma (HNSCC) is highly infiltrated by immune cells, including tumor-infiltrating lymphocytes and myeloid lineage cells. In the tumor microenvironment, tumor cells orchestrate a highly immunosuppressive microenvironment by secreting immunosuppressive mediators, expressing immune checkpoint ligands, and downregulating human leukocyte antigen expression. In the present study, we aimed to comprehensively profile the immune microenvironment of HNSCC using RNA-sequencing (RNA-seq) data obtained from The Cancer Genome Atlas (TCGA) database.Methods: We calculated enrichment scores of 33 immune cell types based on RNA-seq data of HNSCC tissues and adjacent non-cancer tissues. Based on these scores, we performed non-supervised clustering and identified three immune signatures, i.e., cold, lymphocyte, and myeloid/dendritic cell (DC), using clustering results. We then compared the clinical and biological features of the three signatures.Results: Among HNSCC and non-cancer tissues, human papillomavirus (HPV)-positive HNSCCs exhibited the highest scores in various immune cell types, including CD4+ T cells, CD8+ T cells, B cells, plasma cells, basophils, and their subpopulations. Among the three immune signatures, the proportions of HPV-positive tumors, oropharyngeal cancers, early T tumors, and N factor positive cases were significantly higher in the lymphocyte signature than in other signatures. Among the three signatures, the lymphocyte signature showed the longest overall survival (OS), especially in HPV-positive patients, whereas the myeloid/DC signature demonstrated the shortest OS in these patients. Gene set enrichment analysis revealed the upregulation of several pathways related to inflammatory and proinflammatory responses in the lymphocyte signature. The expression of PRF1, IFNG, GZMB, PDCD1, LAG3, CTLA4, HAVCR2, and TIGIT was the highest in the lymphocyte signature. Meanwhile, the expression of PD-1 ligand genes CD274 and PDCD1LG2 was highest in the myeloid/DC signature. Conclusions: Herein, our findings revealed the transcriptomic landscape of the immune microenvironment that closely reflects the clinical and biological significance of HNSCC, indicating that molecular profiling of the immune microenvironment can be employed to develop novel biomarkers and precision immunotherapies for HNSCC.


eLife ◽  
2017 ◽  
Vol 6 ◽  
Author(s):  
Julien Racle ◽  
Kaat de Jonge ◽  
Petra Baumgaertner ◽  
Daniel E Speiser ◽  
David Gfeller

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. 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 (http://epic.gfellerlab.org).


2019 ◽  
Author(s):  
Samuel A Danziger ◽  
David L Gibbs ◽  
Ilya Shmulevich ◽  
Mark McConnell ◽  
Matthew WB Trotter ◽  
...  

AbstractImmune cell infiltration of tumors can be an important component for determining patient outcomes, e.g. by inferring immune cell presence by deconvolving gene expression data drawn from a heterogenous mix of cell types. One particularly powerful family of deconvolution techniques uses signature matrices of genes that uniquely identify each cell type as determined from cell type purified gene expression data. Many methods of this type have been recently published, often including new signature matrices appropriate for a single purpose, such as investigating a specific type of tumor. The package ADAPTS helps users make the most of this expanding knowledge base by introducing a framework for cell type deconvolution. ADAPTS implements modular tools for customizing signature matrices for new tissue types by adding custom cell types or building new matrices de novo, including from single cell RNAseq data. It includes a common interface to several popular deconvolution algorithms that use a signature matrix to estimate the proportion of cell types present in heterogenous samples. ADAPTS also implements a novel method for clustering cell types into groups that are hard to distinguish by deconvolution and then re-splitting those clusters using hierarchical deconvolution. We demonstrate that the techniques implemented in ADAPTS improve the ability to reconstruct the cell types present in a single cell RNAseq data set in a blind predictive analysis. ADAPTS is currently available for use in R on CRAN and GitHub.


2017 ◽  
Author(s):  
Julien Racle ◽  
Kaat de Jonge ◽  
Petra Baumgaertner ◽  
Daniel E. Speiser ◽  
David Gfeller

AbstractImmune 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. 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.


2021 ◽  
Author(s):  
Francesca Petralia ◽  
Anna P Calinawan ◽  
Song Feng ◽  
Sara JC Gosline ◽  
Pietro Pugliese ◽  
...  

Characterizing the tumor microenvironment is crucial in order to improve responsiveness to immunotherapy and develop new therapeutic strategies. The fraction of different cell-types in the tumor microenvironment can be estimated based on transcriptomic profiling of bulk tumor data via deconvolution algorithms. One class of such algorithms, known as reference-based, rely on a reference signature containing gene expression data for various cell-types. The limitation of these methods is that such a signature is derived from the gene expression of pure cell-types, which might not be consistent with the transcriptomic profiling in solid tumors. On the other hand, reference-free methods usually require only a set of cell-specific markers to perform deconvolution; however, once the different components have been estimated from the data, their labeling can be problematic. To overcome these limitations, we propose BayesDeBulk - a new reference-free Bayesian method for bulk deconvolution based on gene expression data. Given a list of markers expressed in each cell-type (cell-specific markers), a repulsive prior is placed on the mean of gene expression in different cell-types to ensure that cell-specific markers are upregulated in a particular component. Contrary to existing reference-free methods, the labeling of different components is decided a priori through a repulsive prior. Furthermore, the advantage over reference-based algorithms is that the cell fractions as well as the gene expression of different cells are estimated from the data, simultaneously. Given its flexibility, BayesDeBulk can be utilized to perform bulk deconvolution beyond transcriptomic data, based on other data types such as proteomic profiles or the integration of both transcriptomic and proteomic profiles.


2022 ◽  
Author(s):  
Yongsheng Zhang ◽  
Yunlong Wang ◽  
Jichuang Wang ◽  
Kaixiang Zhang

Abstract Bladder cancer (BLCA) is among the most frequent types of cancer. Patients with BLCA have a significant recurrence rate and a poor post-surgery survival rate. Recent research has found a link between tumor immune cell infiltration (ICI) and the prognosis of BLCA patients. However, the ICI picture of BLCA remains unclear. Common gene expression data was obtained by combining the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) expression databases. Two computational algorithms were proposed to unravel the ICI landscape of BLCA patients. The R package "limma" was applied to find differentially expressed genes (DEGs). Principal-component analysis (PCA) was used to calculate the ICI score. A total of 569 common gene expression data were retrieved from TCGA and GEO cohorts. CD8+ T cells were found to have a substantial positive connection with activated memory CD4+ T cells and immune score. On the contrary, CD8+ T cells were found to have a substantial negative connection with Macrophages M0. Thirty-eight DEGs were selected. Two ICI patterns were defined by unsupervised clustering method. Patients of BLCA were separated into two groups. The high ICI score group exhibits better outcome than the low one (p < 0.001). Finally, the group with a high tumor mutation burden (TMB) as well as a high ICI score had the best outcome. (p <0.001). Combining TMB and ICI score resulted in a more accurate survival prediction, suggesting that ICI score could be used as a prognostic marker for BLCA patients.


2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Gaia Griguolo ◽  
Maria Vittoria Dieci ◽  
Laia Paré ◽  
Federica Miglietta ◽  
Daniele Giulio Generali ◽  
...  

AbstractLittle is known regarding the interaction between immune microenvironment and tumor biology in hormone receptor (HR)+/HER2− breast cancer (BC). We here assess pretreatment gene-expression data from 66 HR+/HER2− early BCs from the LETLOB trial and show that non-luminal tumors (HER2-enriched, Basal-like) present higher tumor-infiltrating lymphocyte levels than luminal tumors. Moreover, significant differences in immune infiltrate composition, assessed by CIBERSORT, were observed: non-luminal tumors showed a more proinflammatory antitumor immune infiltrate composition than luminal ones.


Sign in / Sign up

Export Citation Format

Share Document