The current research aimed to development and validation in signature immune genes for lymphatic metastasis in papillary thyroid cancer (PTC).
Weighted correlation network analysis (WGCNA) was performed to identify genes closely correlated with lymphatic metastasis in PTC from TCGA database. Information on immune-related genes (IRGs) was obtained from the ImmPort database. Crossover genes were used with the R package clusterProfiler for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment. Key genes in the protein–protein interaction network of cross-targets were obtained using Cytoscape. Lasso and Random Forest (RF) models were utilized to identify pivotal genes. We constructed a nomogram based on the hub genes. The correlation between hub genes and immune cell infiltration was explored. We collected and assessed clinical samples
In total, 122 IRGs were correlated with lymphatic metastases from PTC. There are 10 key IRGs in the protein–protein interaction network. Then, three hub genes including PTGS2, MET, and ICAM1 were established using the LASSO and RF models. The expression of these hub genes was upregulated in samples collected from patients with lymphatic metastases. The average area under the curve of the model reached 0.83 after a 10-fold and 200-time cross-validation, which had a good prediction ability. Immuno-infiltration analysis showed that the three hub genes were significantly positively correlated with resting dendritic cells and were negatively correlated with activated natural cells, monocytes, and eosinophils. Immunohistochemistry results revealed that lymph node metastasis samples had a higher expression of the three hub genes than non-metastasis samples.