Background:
Papillary thyroid cancer has been associated with chronic inflammation. A systematic understanding of immune cell infiltration in PTC is essential for subsequent immune research and new diagnostic and therapeutic strategies.
Methods:
Three different algorithms, single-sample gene set enrichment analysis (ssGSEA), immune cell marker and CIBERSORT, were used to evaluate immune cell infiltration levels (abundance and proportion) in 10 data sets (The Cancer Genome Atlas [TCGA], GSE3467, GSE3678, GSE5364, GSE27155, GSE33630, GSE50901, GSE53157, GSE58545, and GSE60542; a total of 799 PTC and 194 normal thyroid samples). Consensus unsupervised clustering divided PTC patients into low-immunity and high-immunity groups. Weighted gene coexpression network analysis (WGCNA) and gene set enrichment analysis (GSEA) were used to analyze the potential mechanisms causing differences in the immune response.
Results:
Compared with normal tissues, PTC tissues had a higher overall immune level and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, dendritic cells (DCs), mast cells (MCs), and M0 macrophages. Compared with early PTC, advanced PTC showed higher immune infiltration and higher abundance levels and proportions of M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages. Compared to the low-immunity group, the high-immunity group exhibited more advanced stages, larger tumor sizes, greater lymph node metastases, higher tall-cell PTCs, lower follicular PTC proportions, more BRAF mutations, and fewer RAS mutations. Epstein-Barr virus (EBV) infection was the most significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for key module genes.
Conclusions:
In human PTC, M2 macrophages, Tregs, monocytes, neutrophils, DCs, MCs, and M0 macrophages appear to play a tumor-promoting role, while M1 macrophages, CD8+ T cells, B cells, NK cells, and T follicular helper (T) cells (including eosinophils, γδ T cells, and Th17 cells with weak supporting evidence) appear to play an antitumor role. During the occurrence and development of PTC, the overall immune level was increased, and the abundance and proportion of tumor-promoting immune cells were significantly increased, indicating that immune escape had been aggravated. Finally, we speculate that EBV may play an important role in changing the immune microenvironment of PTC tumors.
Citing Articles
Thyroid Cancer-The Tumor Immune Microenvironment (TIME) over Time and Space.
Jaume J
Cancers (Basel). 2025; 17(5).
PMID: 40075642
PMC: 11899416.
DOI: 10.3390/cancers17050794.
Circular RNAs in cancer.
Guo Y, Huang Q, Heng Y, Zhou Y, Chen H, Xu C
MedComm (2020). 2025; 6(2):e70079.
PMID: 39901896
PMC: 11788016.
DOI: 10.1002/mco2.70079.
Single-cell RNA-sequencing and spatial transcriptomic analysis reveal a distinct population of APOE cells yielding pathological lymph node metastasis in papillary thyroid cancer.
Xiao G, Xie R, Gu J, Huang Y, Ding M, Shen D
Clin Transl Med. 2025; 15(1):e70172.
PMID: 39810624
PMC: 11733439.
DOI: 10.1002/ctm2.70172.
The expression of CCL17 and potential prognostic value on tumor immunity in thyroid carcinoma based on bioinformatics analysis.
Gu X, Chen B, Zhang S, Zhai X, Hu Y, Ye H
Sci Rep. 2024; 14(1):31580.
PMID: 39738081
PMC: 11686015.
DOI: 10.1038/s41598-024-75750-1.
Prognostic Protein Biomarker Screening for Thyroid Carcinoma Based on Cancer Proteomics Profiles.
Xie P, Yin Q, Wang S, Song D
Biomedicines. 2024; 12(9).
PMID: 39335579
PMC: 11428938.
DOI: 10.3390/biomedicines12092066.
Immune microenvironment in papillary thyroid carcinoma: roles of immune cells and checkpoints in disease progression and therapeutic implications.
Zheng X, Sun R, Wei T
Front Immunol. 2024; 15:1438235.
PMID: 39290709
PMC: 11405226.
DOI: 10.3389/fimmu.2024.1438235.
Exploring markers of immunoresponsiveness in papillary thyroid carcinoma and future treatment strategies.
Mohanty A, Afkhami M, Reyes A, Pharaon R, Yin H, Li H
J Immunother Cancer. 2024; 12(7).
PMID: 39074963
PMC: 11288153.
DOI: 10.1136/jitc-2023-008505.
OIP5-AS1/miR-455-3p/microfibril-associated protein 2 axis exacerbates the progression of thyroid carcinoma.
Huang R, Liu H, Wang C
Transl Cancer Res. 2024; 13(6):3046-3061.
PMID: 38988919
PMC: 11231809.
DOI: 10.21037/tcr-24-630.
Identification and Validation of Immune Implication of R-Spondin 1 and an R-Spondin 1-Related Prognostic Signature in Esophagus Cancer.
Lin Y, Lou X, Li S, Cai W, Che T
Int J Genomics. 2024; 2024:7974277.
PMID: 38962149
PMC: 11222003.
DOI: 10.1155/2024/7974277.
Predominance of CD4+ T cells in metastatic cervical lymph nodes in papillary thyroid carcinoma.
Moma C, Barreto I, Assumpcao L, Zantut-Wittmann D
Endocr Connect. 2024; 13(8).
PMID: 38913547
PMC: 11301539.
DOI: 10.1530/EC-24-0135.
The Potential Role of Virus Infection in the Progression of Thyroid Cancer.
Wu Y, Jiang T, Su Y, Mei L, Sun T, Li Y
World J Oncol. 2024; 15(3):382-393.
PMID: 38751704
PMC: 11092407.
DOI: 10.14740/wjon1830.
Identification of SUMOylation-related biomarkers in papillary thyroid carcinoma.
Li X, Ding Z, Tong Y
Cancer Cell Int. 2024; 24(1):149.
PMID: 38671425
PMC: 11055338.
DOI: 10.1186/s12935-024-03323-3.
The Prognostic Value and Potential Immune Mechanisms of lncRNAs Related to Immunogenic Cell Death in Papillary Thyroid Carcinoma.
Wang Y, Li X, Huang Y, Gang Q, Liu M, Zhang H
J Inflamm Res. 2024; 17:1995-2008.
PMID: 38566983
PMC: 10986630.
DOI: 10.2147/JIR.S456452.
Comprehensive analysis of cellular senescence and immune microenvironment in papillary thyroid carcinoma.
Huang Y, Jiang H, Xu G, Li X, Chen W, Lun Y
Aging (Albany NY). 2024; 16(3):2866-2886.
PMID: 38329430
PMC: 10911381.
DOI: 10.18632/aging.205520.
Bioinformatic identification reveals a m6A-binding protein, IGF2BP2, as a novel tumor-promoting gene signature in thyroid carcinoma.
Xie Y, Xiao J, Ying Y, Liu J, Zhang L, Zeng X
Naunyn Schmiedebergs Arch Pharmacol. 2024; 397(8):5663-5676.
PMID: 38289368
DOI: 10.1007/s00210-024-02961-8.
Autoimmune Thyroid Disease and Differentiated Thyroid Carcinoma: A Review of the Mechanisms That Explain an Intriguing and Exciting Relationship.
Vargas-Uricoechea H
World J Oncol. 2024; 15(1):14-27.
PMID: 38274715
PMC: 10807914.
DOI: 10.14740/wjon1728.
Integrative Analysis of Single-Cell and Bulk RNA Sequencing Reveals Prognostic Characteristics of Macrophage Polarization-Related Genes in Lung Adenocarcinoma.
Mi K, Zeng L, Chen Y, Yang S
Int J Gen Med. 2023; 16:5031-5050.
PMID: 37942473
PMC: 10629586.
DOI: 10.2147/IJGM.S430408.
Cuproptosis-related prognostic signatures predict the prognosis and immunotherapy in HCC patients.
Peng H, Zou Z, Xiang Z, Lu X, Zhang Y, Peng X
Medicine (Baltimore). 2023; 102(34):e34741.
PMID: 37653738
PMC: 10470811.
DOI: 10.1097/MD.0000000000034741.
Mechanism of immune escape mediated by receptor tyrosine kinase KIT in thyroid cancer.
Luo Z, Xu J, Xu D, Xu J, Zhou R, Deng K
Immun Inflamm Dis. 2023; 11(7):e851.
PMID: 37506147
PMC: 10336654.
DOI: 10.1002/iid3.851.
Cuproptosis-related lncRNAs predict prognosis and immune response of thyroid carcinoma.
Shi Y, Sheng P, Guo M, Chen K, Zhou H, Wu M
Front Genet. 2023; 14:1100909.
PMID: 37470034
PMC: 10352785.
DOI: 10.3389/fgene.2023.1100909.