Reveal the potential molecular mechanism of circRNA regulating immune-related mRNA through sponge miRNA in the occurrence and immune regulation of papillary thyroid cancer

Abstract

Background Papillary thyroid cancer (PTC) is the most common endocrine malignant tumour. The purpose of this study was to explore the potential molecular mechanism of circRNA regulating immune-related mRNA through sponge miRNA in the occurrence and immune regulation of PTC. Methods All data were downloaded from public databases, such as GEO, Immport and TCGA. Differentially expressed (DE) mRNAs (DEmRNAs), DEmiRNAs and DEcircRNAs were identified using metaMA and limma packages. Subsequently, immune-related DEmRNAs were screened, and circRNA-miRNA-mRNA (ceRNA) regulatory network was constructed. In addition, functional annotation, protein-protein interaction (PPI) network construction, immune cell infiltration analysis and Pearson correlation analysis were performed. Finally, qRT-PCR validation and cell experiments were also performed. Results In total, 2962 DEmRNAs, 78 DEmiRNAs and 51 DEcircRNAs were obtained. Subsequently, 195 immune-related DEmRNAs were obtained based on Immport database. Cytokine-cytokine receptor interaction was the only signalling pathway obtained in KEGG analysis. Then, 8 hub immune-related DEmRNAs were identified based on PPI network and CytoHubba plug-in. Subsequently, ceRNA sub-network containing hub immune-related DEmRNAs was extracted from ceRNA regulatory network. In ceRNA sub-network, hsa_circ_0082182-hsa-miR-18b-5p-FGF1/PDGFC, hsa_circ_0016404-hsa-miR-1275-FGF1/CTSB/IL13RA1, hsa_circ_0070100-hsa-miR-27a-3p/hsa-miR-27b-3p-TGFBR3, hsa_circ_0060055/hsa_circ_0038718-hsa-miR-150-3p-CXCL14, hsa_circ_0030427/hsa_circ_0002917-hsa-miR-22-3p-BMP7 and hsa_circ_0030427/hsa_circ_0002917-hsa-miR-125a-5p-LIFR axes were identified. Moreover, FGF1, PDGFC, CTSB, IL13RA1, TGFBR3, CXCL14, BMP7, LIFR, hsa-miR-125a-5p, hsa-miR-1275, hsa-miR-150-3p, hsa-miR-18b-5p and hsa-miR-27b-3p were also found to have good diagnostic accuracy and may be potential novel diagnostic markers for PTC. XCell analysis showed that the levels of immune cell infiltration (including Tregs, HSC, DC and Monocytes) were significantly different between the PTC and the control groups. Knockdown of the expression of hsa_circ_0082182 significantly inhibits the activity, proliferation, migration and invasion of TPC-1 cells. Conclusion Several circRNA-miRNA-mRNA axes identified in this study may be related to the occurrence, progression and survival of PTC. This lays a theoretical foundation for further understanding the molecular mechanism of PTC, and also contributes to clinical management and research.