Comprehensive Analysis of the Functions and Prognostic Value of RNA-Binding Proteins in Thyroid Cancer

RNA binding proteins (RBPs) have been proved to play pivotal roles in a variety types of tumors. However, there is no convincible evidence disclosing the functions of RBPs in thyroid cancer (THCA) thoroughly and systematically. Integrated analysis of the functional and prognostic effect of RBPs help better understanding tumorigenesis and development in thyroid and may provide a novel therapeutic method for THCA. In this study, we obtained a list of human RBPs from Gerstberger database, which covered 1,542 genes encoding RBPs. Gene expression data of THCA was downloaded from The Cancer Genome Atlas (TCGA, n = 567), from which we extracted 1,491 RBPs' gene expression data. We analyzed differentially expressed RBPs using R package "limma". Based on differentially expressed RBPs, we constructed protein-protein interaction network and the GO and KEGG pathway enrichment analyses were carried out. We found six RBPs (AZGP1, IGF2BP2, MEX3A, NUDT16, NUP153, USB1) independently associated with prognosis of patients with thyroid cancer according to univariate and multivariate Cox proportional hazards regression models. The survival analysis and risk score analysis achieved good performances from this six-gene prognostic model. Nomogram was constructed to guide clinical decision in practice. Finally, biological experiments disclosed that NUP153 and USB1 can significantly impact cancer cell proliferation and migration. In conclusion, our research provided a new insight of thyroid tumorigenesis and development based on analyses of RBPs. More importantly, the six-gene model may play an important role in clinical practice in the future.