Abstract
OBJECTIVE: Since there are few studies exploring genes associated with lymphatic metastasis of thyroid carcinoma (THCA), this study was conducted to explore genes associated with lymphatic metastasis of THCA and to investigate the relationship with immune infiltration. METHODS: Differentially expressed genes associated with THCA lymphatic metastasis were analyzed based on The Cancer Genome Atlas Program (TCGA) database; a protein-protein interaction(PPI)network was constructed to screen for pivotal genes. Based on the identified hub genes, their expression in THCA with and without lymphatic metastasis were determined. Functional enrichment analysis was performed. The correlation between the identified genes and immune cell infiltration was explored. LASSO logistic regression analysis was performed to determine the risk score of the most relevant gene constructs and multifactor COX regression analysis based on genes in the risk score formula. RESULTS: A total of 115 genes were differentially expressed in THCA with and without lymphatic metastasis, including 28 upregulated genes and 87 downregulated genes. The PPI network identified seven hub genes (EVA1A, TIMP1, SERPINA1, FAM20A, FN1, TNC, MXRA8); the expression of all seven genes was upregulated in the group with lymphatic metastasis; Immuno-infiltration analysis showed that all seven genes were significantly positively correlated with macrophage M1 and NK cells and negatively correlated with T-cell CD4+ and myeloid dendritic cells. LASSO logistic regression analysis identified the five most relevant genes (EVA1A, SERPINA1, FN1, TNC, MXRA8), and multi-factor COX regression analysis showed EVA1A, SERPINA1 and FN1 as independent prognostic factors. CONCLUSION: Seven genes were associated with lymphatic metastasis of THCA and with tumor immune cell infiltration.