Identification of Transcriptional Pattern Related to Immune Cell Infiltration With Gene Co-Expression Network in Papillary Thyroid Cancer

Abstract

BackgroundA growing body of evidence suggests that immune cell infiltration in cancer is closely related to clinical outcomes. However, there is still a lack of research on papillary thyroid cancer (PTC).MethodsBased on single-sample gene set enrichment analysis (SSGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) tool, the infiltration level of immune cell and key modules and genes associated with the level of immune cell infiltration were identified using PTC gene expression data from The Cancer Genome Atlas (TCGA) database. In addition, the co-expression network and protein-protein interactions network analysis were used to identify the hub genes. Moreover, the immunological and clinical characteristics of these hub genes were verified in TCGA and GSE35570 datasets and quantitative real-time polymerase chain reaction (PCR). Finally, receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic value of hub genes.ResultsActivated B cell, activated dendritic cell, CD56bright natural killer cell, CD56dim natural killer cell, Eosinophil, Gamma delta T cell, Immature dendritic cell, Macrophage, Mast cell, Monocyte, Natural killer cell, Neutrophil and Type 17 T helper cell were significantly changed between PTC and adjacent normal groups. WGCNA results showed that the black model had the highest correlation with the infiltration level of activated dendritic cells. We found 14 hub genes whose expression correlated to the infiltration level of activated dendritic cells in both TCGA and GSE35570 datasets. After validation in the TCGA dataset, the expression level of only 5 genes (C1QA, HCK, HLA-DRA, ITGB2 and TYROBP) in 14 hub genes were differentially expressed between PTC and adjacent normal groups. Meanwhile, the expression levels of these 5 hub genes were successfully validated in GSE35570 dataset. Quantitative real-time PCR results showed the expression of these 4 hub genes (except C1QA) was consistent with the results in TCGA and GSE35570 dataset. Finally, these 4 hub genes had diagnostic value to distinguish PTC and adjacent normal controls.ConclusionsHCK, HLA-DRA, ITGB2 and TYROBP may be key diagnostic biomarkers and immunotherapy targets in PTC.

Bioinformatics analysis of the prognostic and clinical value of senescence-related gene signature in papillary thyroid cancer

“…[18] The infiltrating score of 16 immune cells and the activity of 13 immune-related pathways were performed by the “single-sample gene set enrichment analysis” function in “gsva” R package. [19]…”

Construction of ferroptosis-related prediction model for pathogenesis, diagnosis and treatment of ruptured abdominal aortic aneurysm

“…[16] The "single-sample gene set enrichment analysis (ssGSEA)" was performed to analyzed infiltrating score of 16 immune cells and activity of 13 immune-related functions between stable group and rupture group. [17] 2.9. Drug analysis DGIdb (https://dgidb.genome.wustl.edu/) is an open-source project that offers user-friendly browsing, searching, and filtering of information on drug-gene interactions and the druggable genome.…”

Prognosis and therapy in thyroid cancer by gene signatures related to natural killer cells

“…Additionally, previous studies have shown a strong positive association between improved clinical outcomes across various types of tumors and higher levels of tumor‐infiltrating NK cells 12,13 . Although the molecular properties of NK cells in malignant tumors and infectious disorders have been studied, the association between NK cells and the prognosis of THCA patients remains unclear 14,15 …”

“…12,13 Although the molecular properties of NK cells in malignant tumors and infectious disorders have been studied, the association between NK cells and the prognosis of THCA patients remains unclear. 14,15 With the development of bioinformatics, biomarkers have been defined in various ways. On the basis of RNA-seq expression profiles, there is significant evidence that analyzing the characteristics of infiltrating immune cells could be valuable in predicting the prognosis and response to immunotherapy in patients with malignant tumors.…”