A Potential Four-Gene Signature and Nomogram for Predicting the Overall Survival of Papillary Thyroid Cancer

Authors: Ying Ding ⚬, Peng Li ⚬, Wenlong Wang ⚬, Fang Liu ⚬

Affiliations: Thyroid Surgery Department, Xiangya Hospital, Central South University, Changsha 410008, China; National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, 410008 Hunan Province, China; Health Management Center, Xiangya Hospital, Central South University, Changsha 410008, China

License: Copyright © 2022 Ying Ding et al. CC BY 4.0 This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Article links: DOI: 10.1155/2022/8735551  |  PubMed: 36193505  |  PMC: PMC9526076

Relevance: Moderate: mentioned 3+ times in text

Full text: PDF (5.5 MB)

1. Introduction

In the past few decades, the incidence of thyroid cancer has increased rapidly worldwide [ref. 1]. Papillary thyroid carcinoma (PTC) is the most common pathological subtype, which has accounted for more than 85% of cases [ref. 2]. PTCs usually presents an excellent prognosis, and 10-year disease-specific survival rates have been reported to be over 90% via management through the common therapeutic approaches such as thyroidectomy, RAI therapy, and thyroid-stimulating hormone (TSH) suppressive therapy [ref. 3]. However, local recurrence and distant metastasis inevitably occur in up to 20% and 10% in PTC patients [ref. 4]. Moreover, two-thirds of these patients exhibit loss of iodine-131 (131I) uptake initially or gradually, thereby indicating dedifferentiation of the PTC termed as RAI-refractory PTC [ref. 5]. Thus, accurate assessment of the prognosis of PTC is critical to ensure that the high-risk patients receive appropriate treatment and to prevent excessive treatment of the low-risk patients.

The 8^th edition of the AJCC/TNM (tumor node metastasis) manual was released in 2017 [ref. 6, ref. 7]. TNM staging has been identified to be useful in predicting the disease mortality, and it is recommended for all PTC patients [ref. 7, ref. 8]. However, it has been difficult to accurately distinguish the difference between the survival outcomes in PTCs with similar clinicopathological features [ref. 9, ref. 10]. A number of prior studies have proposed that combining BRAF, TERT, and RAS mutations with TNM staging can lead to better prediction of the prognosis of PTC patients [ref. 11–ref. 13]. Our team has previously constructed an integrated nomogram based on clinicopathological factors and the related risk scores, which showed significantly better predictive performance than AJCC stage [ref. 14–ref. 16].

In this study, we aimed to develop a prognosis-predicting model based on nomogram as a prognostic evaluation method. We have obtained data of PTC patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Univariate and multivariate Cox regression analysis and Lasso regression analysis was performed to establish a four-gene signature in the training cohort. Time-dependent receiver-operating characteristic (ROC) and Kaplan-Meier (KM) curve was used to assess the validity of the four-gene signature and the nomogram in the entire patient cohort.

2. Materials and Methods

2.1. Acquisition and Processing of Data

ALL RNA-seq data and the clinical characteristics were extracted from three GEO datasets (GSE33630, GSE3678, and GSE60542) and TCGA-THCA dataset (Table S1). A total of 494 PTC samples were selected from the TCGA-THCA dataset and were then enrolled for the subsequent analysis. The flowchart of the designed study has been illustrated in Figure 1.

2.2. Screening of the Differentially Expressed Genes and Gene Enrichment Analysis

Differentially expressed genes (DEGs) between the normal and PTC samples with absolute log2 fold change (FC) > 1.5 and P < 0.05 were screened from three GEO datasets and the TCGA-THCA dataset by using the “limma” R package. Enrichment analyses of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for DEGs were performed by using the “clusterProfiler” R package and Webgestalt website (http://www.webgestalt.org/).

2.3. Establishment and Validation of the Four-Gene Risk Signature

In total, 344 samples were randomly selected from the 494 PTC samples in the TCGA-THCA dataset as the training cohort. The remaining 150 samples were used as the testing cohort. In the training cohort, 35 OS-related genes were selected based on univariate cox analysis. Thereafter, least absolute shrinkage and selection operator (LASSO) regression was used to screen 12 potential prognostic genes. Finally, only 4 genes (P < 0.05) were included in the risk signature based on the results of multivariate Cox regression. Finally, the risk score for each patient was calculated as follows: ∑_i=1ⁿcoef(i)∗exp(i), where exp(i) represents the expression of genes, and coef(i) is the coefficient of multivariate Cox regression. ROC and KM curves for risk score were drawn with the “survival” and “timeROC” R packages to assess the potential predictive capacity of the risk signature in the training cohort, testing cohort, and the entire cohort.

2.4. Construction and Evaluation of the Nomogram

Univariate and multivariate Cox regression analyses were utilized to screen the various essential clinical characteristics related to OS. A nomogram was established based on the risk score, age, AJCC stage, tumor size, extrathyroidal extension, and history of neoadjuvant treatment with the “rms” R package. The C-index, Akaike information criterion (AIC), and Bayesian Information Criterion (BIC) of the nomogram were calculated, and the ROC, KM, calibration, and decision curves were drawn.

2.5. Immune Analysis

The median nomogram point was used to divide the entire cohort into two distinct groups of good-prognosis and poor-prognosis. Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORTx) (https://cibersortx.stanford.edu/) was then employed to analyze the immune infiltration of 22 types of immune cells in good-prognosis and poor-prognosis groups. The stroma and immune scores were measured by Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) analysis using “estimate” R package.

2.6. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA)

GSEA and GSVA were performed with the GSEA software and the “GSVA” R package. The reference gene sets in GSEA were the c2.cp.kegg.v7.5.1.symbols.gmt, and if the normal P value < 0.05 and FDR (false discovery rate) q value < 0.25, the gene set was considered as significantly enriched.

2.7. Construction of a Competing Endogenous RNA Regulatory Network

The (differentially expressed miRNAs) DEmiRNAs and (differentially expressed lncRNAs) DElncRNAs between the normal and PTC samples in the TCGA-THCA dataset with absolute log2 fold change (FC) > 1.5 and P < 0.05 were screened by using the “limma” R package. The miRcode database was used to match DElncRNAs and DEmiRNAs. ABI3BP, DPT, MRO, and TENM1 target miRNAs were predicted based on the three distinct databases: miRWalk, TargetScan, and miRmap. Subsequently, a ceRNA regulatory network was constructed according to the results obtained above. Cytoscape (version 3.8.2) was used to visualize the competing endogenous RNA (ceRNA) network. OS-related DEmiRNAs and DElncRNAs were screened by the KM curves.

2.8. Data Mining in the GEPIA and HPA Databases

The RNA expression data of ABI3BP, DPT, MRO, and TENM1 genes in normal and PTC samples were extracted from the Gene Expression Profiling Interactive Analysis (GEPIA) database. The immunohistochemistry (IHC) data was obtained from seven patients in the Human Protein Atlas (HPA) database, for whom the basic information is available as shown in Table S2.

2.9. Patients

Ten paired PTC tumors and adjacent normal thyroid tissues were obtained from patients in the Thyroid Surgery Department of Xiangya Hospital from March 2020 to June 2020. An informed consent was obtained from all the participants, and the study was approved by the Ethics Committee of Xiangya Hospital of Central South University (No. 202004192).

2.10. Quantitative RT-PCR

Total RNA was extracted from the normal and PTC tissues by using Trizol Reagent (Accurate Biology, China). cDNA synthesis was performed using Reverse Transcription Kit (Accurate Biology, China). Quantitative RT-PCR was then carried out with the RT-PCR Kit (Accurate Biology, China). The sequence of all the primers used have been described in Table S3.

2.11. Statistical Analysis

Statistical analysis was performed based on R software (version 4.1.0). Chi-square test or Fisher’s exact test were used to analyze the categorical variables. T-test and one-way ANOVA were used to analyze the continuous variables. Univariate and multivariate Cox regression and log-rank test were performed to evaluate OS. Unless otherwise stated, P < 0.05 indicated that the difference was statistically significant.

3. Results

3.1. Identification of DEGs and Gene Enrichment Analysis

The value distribution of the selected samples in the three GEO datasets is uniform which has been shown in Supplementary Figure 1. The various DEGs were obtained in three GEO datasets and TCGA-THCA datasets (Supplementary Figure 2A-D) (453 in GSE33630, 323 in GSE3678, 442 in GSE60542, and 2328 in the TCGA-THCA dataset). The intersection of the four datasets contained 176 DEGs (Supplementary Figure 2E), including 102 upregulated genes and 74 downregulated genes. GO analysis demonstrated that these 176 DEGs were primarily enriched in cell junction assembly, synapse organization, and positive regulation of protein serine in the biological process, collagen-containing extracellular matrix, secretory granule lumen, and cytoplasmic vesicle lumen in the cellular components, serine-type peptidase activity, serine hydrolase activity, and serine-type endopeptidase activity in the molecular functions (Figures 2(a)–2(c)). Additionally, KEGG pathway enrichment analysis indicated that DEGs were significantly enriched in tyrosine metabolism (Figure 2(d)). These results suggested that the genesis and development of PTC might be closely related to these genes or pathways, including pathways that play an important role in development of cancers, such as serine activity, tyrosine metabolism, extracellular matrix, and cellular junction. Therefore, these results can guide our subsequent study related to understanding the molecular mechanisms of PTC.

3.2. Establishment of the Four-Gene Risk Signature

The baseline of the clinical characteristics of the training cohort, testing cohort, and the entire cohort has been presented in Table 1. 35 genes (P < 0.05) related to OS were screened by univariate Cox regression in the training cohort (Table S4). After performing Lasso regression analysis on these 35 genes, 12 genes were obtained (Figures 3(a) and 3(b)). Finally, 12 genes were analyzed by using multivariate Cox regression, and ABI3BP, DPT, MRO, and TENM1 were selected to construct the four-gene risk signature (Figure 3(c)). The risk score = −0.56124 × ABI3BP + (−0.25805 × DPT) + (−0.38946 × MRO) + 0.55421 × TENM1.

3.3. Validation of the Four-Gene Risk Signature

The four-gene risk signature could effectively stratify patients into a high-risk group and a low-risk group with the median risk score. Thus, we could intuitively conclude that the patients who had died were basically the patients with higher ranking (higher risk score) (Figures 4(a)–4(f)). Consistent with this, significant differences of OS between the high-risk and low-risk groups in the training cohort (P = 5.13 × 10⁻³), testing cohort (P = 2.30 × 10⁻²), and the entire cohort (P = 2.56 × 10⁻⁴) were depicted by constructing KM curves (Figures 4(g)–4(i)). Moreover, the AUCs of the risk score corresponding to 1-year, 3-year, and 5-year survival were 1, 0.847, and 0.832 in the training cohort, 1, 0.836, and 0.947 in the testing cohort, and 0.855, 0.829, and 0.841 in the entire cohort, respectively (Figures 4(j)–4(l)). These results indicated that the four-gene risk signature exhibited useful value in accurately predicting the prognosis of PTC patients.

3.4. Building and Validating a Predictive Nomogram

Thereafter, in the entire cohort, univariate and multivariate Cox regression analyses were performed on the risk score and the clinical characteristics (Table 2). The risk score, age, AJCC stage, tumor size, extrathyroidal extension, and history of neoadjuvant treatment were selected to establish a nomogram (Figure 5(a)). The C-index, AIC, and BIC of this nomogram were 0.970, 85.743, and 86.451, respectively (Table 3). Moreover, a calibration curve revealed that the nomogram was excellent at predicting OS (Figure 5(b)). Moreover, the AUCs of the nomogram corresponding to 1-year, 3-year, and 5-year OS were 0.985, 0.962, and 0.973, respectively (Figure 5(c)). They were found to be significantly better than the traditional (Age, Grade, Extrathyroidal extension, Size) AGES score and (Metastases, Age, Completeness of resection, Invasion, Size) MACIS score (Figures 5(d) and 5(e)). Moreover, the decision curves of 3-year OS and 5-year OS of these three models indicated that in terms of predicting OS of PTC patients, the net benefit [ref. 17] of the nomogram model was substantially higher than that of the traditional AGES and MACIS score (Figures 5(f) and 5(g)). Taken together, these results revealed that this nomogram performed well in predicting the prognosis of PTC patients.

3.5. Immune Analysis

The nomogram could effectively stratify patients into a good-prognosis group and a poor-prognosis group based on the median total point. Surprisingly, immune analysis demonstrated that there was no significant difference in the abundance of 22 types of immune cells between the good-prognosis and poor-prognosis groups (Figures 6(a) and 6(b)). Consistent with this, the stroma scores, immune scores, and ESTIMATE scores were observed to be not significantly higher in the poor-prognosis group compared with the good-prognosis group (Figures 6(c)–6(e)), thus suggesting that the proportion of the stromal cells to immune cells in these two groups may be similar, but there was no difference in the tumor purity. These results implied that immune infiltration and immune microenvironment had minimal effect on the OS of PTC patients.

3.6. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA)

Since the immune infiltration analysis did not yield significant results, GSEA and GSVA were performed to further explore the potential differences in the molecular mechanisms between the good-prognosis and poor-prognosis groups. As shown in Figures 7(a)–7(h), the “ribosome,” “intestinal immune network for IgA production,” “systemic lupus erythematosus,” “asthma,” “nod like receptor signaling,” “glycosaminoglycan degradation,” “viral myocarditis,” and “cell adhesion molecules cams” pathways were found to be significantly enriched in the poor-prognosis group. Moreover, GSVA demonstrated that “peroxisome,” “myc targets V2,” “uv response up,” and “cholesterol homeostasis” were the hallmark pathways modulated in the poor-prognosis group (t value > 2) (Figure 7(i)). This implied that these pathways or the genes regulating them could play an important role in promoting the progression of PTC.

3.7. Construction of a ceRNA Regulatory Network

95 DEmiRNAs (Table S5) and 839 DElncRNAs (Table S6) between the normal and PTC samples in the TCGA-THCA dataset were screened (Figures 8(a) and 8(b)). A total of 54 DElncRNAs and 10 DEmiRNAs were paired into 151 DElncRNA-DEmiRNA interactions, whereas 10 DEmiRNAs and 4 DEmRNAs were matched to form 15 pairs of DEmiRNA-DEmRNA interactions (Table S7). Consequently, the lncRNA-miRNA-mRNA ceRNA regulatory network, which contained 68 distinct nodes and 166 edges was constructed (Figure 8(c)). Among these, DElncRNAs and DEmiRNAs, KM curves showed that one lncRNA (MIR181A2HG) and one miRNA (hsa−mir−375) could serve as potential protective biomarkers in PTC patients (Figures 8(d) and 8(e)).

3.8. Validation of the Expression of the Four Genes in PTC Tissues

Among the four genes in four-gene risk signature, the levels of ABI3BP, DPT, and MRO were downregulated, whereas that of TENM1 was upregulated in PTC samples (Figures 9(a)–9(d)). Thereafter, we used quantitative RT-PCR to verify the expression of these four genes in 10 pairs of PTC tissues (Figures 9(e)–9(h)). Furthermore, immunohistochemical staining data of TENM1 in the normal and tumor tissues were obtained by searching the HPA database, and higher protein expression of TENM1 was observed in analyzed tumor tissues (Figures 9(i)–9(k)).

4. Discussion

It has been established that even when the patients receive the standardized treatment, about 5%-23% of PTC patients display a poor prognosis [ref. 18]. Therefore, prediction of the prognosis of PTC patients can not only promote the active implementation of individualized treatment but also aid to avoid the various negative effects associated with excessive medical treatment. A number of prognostic markers obtained from gene expression profiles can accurately predict the prognosis of a single patient at the molecular level and can be complementary to the traditional clinical staging prediction system such as TNM staging [ref. 19].

In the present study, a novel four-gene risk signature comprising ABI3BP, DPT, MRO, and TENM1 to predict the OS of PTC was identified. The efficacy of this signature was validated in the study cohort. Among these four genes, the levels of ABI3BP, DPT, and MRO were downregulated whereas that of TENM1 was upregulated in PTC. ABI3BP is an ArgBP/E3B1/Avi2/NESH family protein, which can participate in the negative regulation of the cell movement and metastasis through its influence on membrane folding and layer formation [ref. 20, ref. 21]. It has been proven to be an src-homologous 3(SH3) adapter molecule and can exhibit a tumor-suppressive effect in thyroid cancer [ref. 22, ref. 23]. ABI3BP, which is reexpressed in the thyroid cells, has been reported to trigger cellular senescence through affecting the p21 pathway, resulting in a reduction in transformation activity, cell growth, viability, migration, invasion, and tumor growth in nude mice [ref. 22]. Moreover, the loss of ABI3BP expression may be functionally involved in the pathogenesis of several types of cancer such as gallbladder cancer [ref. 24] and esophageal cancer [ref. 25]. DPT is a tyrosine-rich noncollagenous extracellular matrix component, and the depletion of DPT has been associated with hyperproliferation of scars, skin fibrosis, systemic sclerosis as well as some cancers [ref. 26]. DPT has been reported to regulate cell proliferation and invasiveness of a variety of tumors like endometrial cancer [ref. 27], prostate cancer [ref. 28], hepatocellular carcinoma [ref. 29], and oral cancer [ref. 30]. Moreover, low DPT expression in PTC has been related to higher T classification. DPT can regulate CDK4, CDK6, and p21 through modulating MEK-ERK-MYC signaling to inhibit PTC proliferation [ref. 31]. The expression of TENM1 as a cell signal sensor in neurons has been positively correlated with the growth and invasion of PTC. TENM1 has been identified as the direct functional target of miR-486 in PTC cells. The restoration of miR-486 can significantly inhibit the growth of PTC in vivo [ref. 32, ref. 33]. Similar to the previous two genes, TENM1 can also play an important role in different types of cancer especially in the brain tumors such as prolactin pituitary tumor [ref. 34] and glioblastoma [ref. 35]. MRO belongs to a novel gene family, named “maestro heat-like repeat family (MROH)” [ref. 36]. Expression of MRO in lean-type polycystic ovarian syndrome has been found to be increased [ref. 37]. Nevertheless, MRO has been poorly studied in the tumors, and the role of MRO in thyroid cancer has not been previously reported. We speculate that MRO, as a sex-determining gene affecting the prognosis of PTC, might be related to the sex differences between PTC patients, but this observation needs further analysis.

The 8th edition of the AJCC/TNM staging system (TNM-8th) was released in 2017. A number of studies have proved that it was more appropriate for the prediction of the survival and recurrence than TNM-7th [ref. 7, ref. 38]. However, studies also have shown that approximately 30%-40% [ref. 39, ref. 40] of patients were downstaged upon reclassification. In addition, TNM-8th oversimplifies the survival deterioration that primarily occurs with increasing age at the time of diagnosis and underestimates the prognosis for the younger patients, particularly those aged 45–55 years [ref. 9]. What is more, the molecular profile, which is important for the precision medicine, has not been contained in TNM-8th. Consequently, as knowledge of cancer biology evolves, innovative diagnostic tools and treatment modalities need to be developed and improved. Novel nomograms for predicting the survival of thyroid cancer patients have been previously formulated in several studies. For instance, Pathak and colleagues developed nomograms to predict the likelihood of the relapse and death from thyroid cancer in an individual patient (C-indices were 0.92 and 0.76, respectively) based on patients’ characteristics but without identification of any specific gene signature [ref. 41]. A nomogram that first included a gene signature to predict 1-year, 3-year, and 5-year DFS of DTC patients was established (C‐index = 0.801) by Pan Ruchong [ref. 42]. However, five-gene signature appears to cause significant expenditure in the healthcare compared to the four-gene signature. In addition, to the best of our knowledge, our study is the first study to confirm that the developed nomogram performed better than the traditional AGES score and MACIS score for prediction of OS. Besides, our study demonstrated that BRAF mutations had minimal effect on OS of PTC. Consistent with this, unlike previous studies suggesting that BRAF mutations can imply poor prognosis [ref. 43], some recent studies have also indicated that BRAF mutations cannot be used as an independent prognostic and predictive factor in PTC [ref. 44]. Moreover, Wang et al. found that compared to male PTC patients, BRAF mutations cannot be considered as robust independent risk factor for female patients [ref. 45]. In general, the nomogram constructed in this study exhibited better predictive efficacy for the OS of PTC patients and clinical applicability.

To understand the potential mechanisms affecting the prognosis of patients with PTC, immune analysis, GSEA, and GSVA were performed between good-prognosis and poor-prognosis groups. Surprisingly, in contrast to some other prior studies [ref. 46], immune infiltration and immune microenvironment displayed little effect on the OS of PTC patients. These results can be explained, in part, by the fact that the four key genes we identified in this study are not directly involved in the process of tumor immunity. Therefore, no significant biological difference in immune infiltration was observed when PTC samples were divided into good and poor prognosis groups. The “cell adhesion molecules cams” [ref. 47], “Myc” pathways [ref. 48], and glycosaminoglycans (GAGs) [ref. 49], which can play an essential role in the tumor pathogenesis and distant metastasis were found to be enriched in poor-prognosis groups. The effects of intestinal immunity caused by intestinal flora on the various cancers are currently being elucidated, including gastrointestinal tumors, liver cancer [ref. 50], and breast cancer [ref. 51]. In addition, intestinal immunity can also affect the thyroid function [ref. 52]. However, the effect of intestinal immunity on PTC remains to be further studied. Moreover, a number of studies have reported that systemic lupus erythematosus (SLE) and asthma were associated with an increased risk of overall cancers including non-Hodgkin’s lymphoma, Hodgkin’s lymphoma, leukemia, multiple myeloma, and thyroid cancer. Therefore, the increased expression of genes associated with lupus and asthma in the poor-prognosis group is understandable. What is more, the key molecules were presented in the ceRNA network including MIR181A2HG and hsa-mir-375, which need to be studied in detail in PTC.

Although this study revealed numerous important findings, but there are several limitations associated with it. First, selection bias and confounding bias were inevitable due to the retrospective design of this study. Second, the clinical characteristics were mainly derived from TCGA database, and thus, caution should be exercised when expanding our results to patients of other ethnicities. Besides, in future studies, the nomogram should be validated in different external datasets. Finally, additional in vitro and in vivo functional experiments need to be performed to further elucidate the detailed molecular mechanisms affecting the prognosis of PTC.

5. Conclusion

To conclude, our study established a four-gene risk signature and developed a novel prognostic nomogram in combination with prognosis-related clinical characteristics to predict the OS of PTC. The four DEGs were found closely related to the prognosis of PTC and thus can act as potential therapeutic targets. These results might be beneficial for individualized treatment and medical decision-making during the management of PTC patients.

References

1. C. D. Seib, J. A. Sosa. Evolving understanding of the epidemiology of thyroid cancer. Endocrinology and Metabolism Clinics of North America, 2019. [DOI | PubMed]
2. M. E. Cabanillas, D. G. McFadden, C. Durante. Thyroid cancer. Lancet, 2016. [DOI | PubMed]
3. W. Wang, N. Bai, Q. Ouyang, B. Sun, C. Shen, X. Li. Prediction of level V metastases in papillary thyroid microcarcinoma: a single center analysis. Gland Surgery, 2020. [DOI | PubMed]
4. W. Wang, C. Meng, Q. Ouyang, J. Xie, X. Li. Magnesemia: an independent risk factor of hypocalcemia after thyroidectomy. Cancer Management and Research, 2019. [DOI | PubMed]
5. J. M. Oh, B. C. Ahn. Molecular mechanisms of radioactive iodine refractoriness in differentiated thyroid cancer: impaired sodium iodide symporter (NIS) expression owing to altered signaling pathway activity and intracellular localization of NIS. Theranostics, 2021. [DOI | PubMed]
6. E. F. S. van Velsen, W. E. Visser, M. T. Stegenga. Finding the optimal age cutoff for the UICC/AJCC TNM staging system in patients with papillary or follicular thyroid cancer. Thyroid, 2021. [DOI | PubMed]
7. S. H. Nam, M. R. Bae, J. L. Roh. A comparison of the 7th and 8th editions of the AJCC staging system in terms of predicting recurrence and survival in patients with papillary thyroid carcinoma. Oral Oncology, 2018. [DOI | PubMed]
8. A. S. Alzahrani, L. Albalawi, S. Mazi. How does the AJCC/TNM staging system eighth edition perform in thyroid cancer at a major middle eastern medical center?. Endocrine Practice, 2021. [DOI | PubMed]
9. Y. N. Kim, M. Kim, H. S. Ahn. Refining the tumor-node-metastasis staging system for individualized treatment of differentiated thyroid carcinoma. Oral Oncology, 2019. [DOI | PubMed]
10. O. A. Manzardo, M. Cellini, R. Indirli. TNM 8th edition in thyroid cancer staging: is there an improvement in predicting recurrence?. Endocrine-Related Cancer, 2020. [DOI | PubMed]
11. R. Liu, M. Xing. TERT promoter mutations in thyroid cancer. Endocrine-Related Cancer, 2016. [DOI | PubMed]
12. V. Marotta, M. Bifulco, M. Vitale. Significance of RAS mutations in thyroid benign nodules and non-medullary thyroid cancer. Cancers (Basel), 2021. [DOI]
13. J. A. Silver, M. Bogatchenko, M. Pusztaszeri. BRAF V600E mutation is associated with aggressive features in papillary thyroid carcinomas ≤1.5 cm. Journal of Otolaryngology – Head & Neck Surgery, 2021. [DOI]
14. W. Wang, C. Shen, Y. Zhao, B. Sun, N. Bai, X. Li. Identification and validation of potential novel biomarkers to predict distant metastasis in differentiated thyroid cancer. Annals of Translational Medicine, 2021. [DOI | PubMed]
15. W. Wang, Z. Yang, Q. Ouyang. A nomogram to predict skip metastasis in papillary thyroid cancer. World Journal of Surgical Oncology, 2020. [DOI | PubMed]
16. W. Wang, C. Shen, Y. Zhao. The role of m6A RNA methylation-related lncRNAs in the prognosis and tumor immune microenvironment of papillary thyroid carcinoma. Frontiers in Cell and Development Biology, 2021. [DOI]
17. B. Van Calster, L. Wynants, J. F. M. Verbeek. Reporting and interpreting decision curve analysis: a guide for investigators. European Urology, 2018. [DOI | PubMed]
18. M. Haq, C. Harmer. Differentiated thyroid carcinoma with distant metastases at presentation: prognostic factors and outcome. Clinical Endocrinology, 2005. [DOI | PubMed]
19. W. Wang, N. Bai, X. Li. Comprehensive analysis of the prognosis and drug sensitivity of differentiation-related lncRNAs in papillary thyroid cancer. Cancers (Basel), 2022. [DOI]
20. Y. Ichigotani, K. Fujii, M. Hamaguchi, S. Matsuda. In search of a function for the E3B1/Abi2/Argbp1/NESH family (review). International Journal of Molecular Medicine, 2002. [DOI | PubMed]
21. K. Miyazaki, S. Matsuda, Y. Ichigotani. Isolation and characterization of a novel human gene (NESH) which encodes a putative signaling molecule similar to e3B1 protein1. Biochimica et Biophysica Acta, 2000. [DOI | PubMed]
22. F. R. Latini, J. P. Hemerly, G. Oler, G. J. Riggins, J. M. Cerutti. Re-expression of ABI3-binding protein suppresses thyroid tumor growth by promoting senescence and inhibiting invasion. Endocrine-Related Cancer, 2008. [DOI | PubMed]
23. F. R. Latini, J. P. Hemerly, B. C. Freitas, G. Oler, G. J. Riggins, J. M. Cerutti. ABI3 ectopic expression reduces in vitro and in vivo cell growth properties while inducing senescence. BMC Cancer, 2011. [DOI | PubMed]
24. N. Lin, Z. Yao, M. Xu. Long noncoding RNA MALAT1 potentiates growth and inhibits senescence by antagonizing ABI3BP in gallbladder cancer cells. Journal of Experimental & Clinical Cancer Research, 2019. [DOI | PubMed]
25. H. Cai, Y. Li, D. Qin. The depletion of ABI3BP by microRNA-183 promotes the development of esophageal carcinoma. Mediators of Inflammation, 2020. [DOI]
26. G. P. Sidgwick, A. Bayat. Extracellular matrix molecules implicated in hypertrophic and keloid scarring. Journal of the European Academy of Dermatology and Venereology : JEADV, 2012. [DOI | PubMed]
27. H. Huang, Z. Hao, L. Long. Dermatopontin as a potential pathogenic factor in endometrial cancer. Oncology Letters, 2021. [DOI | PubMed]
28. T. Takeuchi, M. Suzuki, J. Kumagai, T. Kamijo, M. Sakai, T. Kitamura. Extracellular matrix dermatopontin modulates prostate cell growth in vivo. The Journal of Endocrinology, 2006. [DOI | PubMed]
29. S. Liu, J. Qiu, G. He. Dermatopontin inhibits WNT signaling pathway via CXXC finger protein 4 in hepatocellular carcinoma. Journal of Cancer, 2020. [DOI | PubMed]
30. M. Yamatoji, A. Kasamatsu, Y. Kouzu. Dermatopontin: a potential predictor for metastasis of human oral cancer. International Journal of Cancer, 2012. [DOI | PubMed]
31. Y. Guo, H. Li, H. Guan. Dermatopontin inhibits papillary thyroid cancer cell proliferation through MYC repression. Molecular and Cellular Endocrinology, 2019. [DOI | PubMed]
32. Y. H. Sun, Z. F. Liu, B. B. Yang, B. Yu. MicroRNA-486 inhibits cell proliferation, invasion and migration via down-regulating the TENM1 expressions and affecting ERK and Akt signaling pathways and epithelial-to-mesenchymal transition in papillary thyroid carcinoma. European Review for Medical and Pharmacological Sciences, 2019. [DOI | PubMed]
33. S. P. Cheng, M. J. Chen, M. N. Chien, C. H. Lin, J. J. Lee, C. L. Liu. Overexpression of teneurin transmembrane protein 1 is a potential marker of disease progression in papillary thyroid carcinoma. Clinical and Experimental Medicine, 2017. [DOI | PubMed]
34. W. Zhang, Z. Zang, Y. Song, H. Yang, Q. Yin. Co-expression network analysis of differentially expressed genes associated with metastasis in prolactin pituitary tumors. Molecular Medicine Reports, 2014. [DOI | PubMed]
35. O. Kadioglu, M. E. M. Saeed, N. Mahmoud. Identification of novel drug resistance mechanisms by genomic and transcriptomic profiling of glioblastoma cells with mutation-activated EGFR. Life Sciences, 2021. [DOI | PubMed]
36. S. Kenigsberg, P. D. Lima, L. Maghen. The elusive MAESTRO gene: its human reproductive tissue-specific expression pattern. PLoS One, 2017. [DOI | PubMed]
37. A. Mitchell, H. Y. Chang, L. Daugherty. The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Research, 2015. [DOI | PubMed]
38. B. G. Cavalheiro, L. L. de Matos, A. K. N. Leite, M. A. V. Kulcsar, C. R. Cernea, L. P. Kowalski. Survival in differentiated thyroid carcinoma: comparison between the 7th and 8th editions of the AJCC/UICC TNM staging system and the ATA initial risk stratification system. Head & Neck, 2021. [DOI | PubMed]
39. R. M. Tuttle, B. Haugen, N. D. Perrier. Updated American Joint Committee on cancer/tumor-node-metastasis staging system for differentiated and anaplastic thyroid cancer (eighth edition): what changed and why?. Thyroid, 2017. [DOI | PubMed]
40. L. Lamartina, G. Grani, E. Arvat. 8th edition of the AJCC/TNM staging system of thyroid cancer: what to expect (ITCO#2). Endocrine-Related Cancer, 2018. [DOI | PubMed]
41. K. A. Pathak, P. Lambert, R. W. Nason, T. Klonisch. Comparing a thyroid prognostic nomogram to the existing staging systems for prediction risk of death from thyroid cancers. European Journal of Surgical Oncology, 2016. [DOI | PubMed]
42. M. Wu, H. Yuan, X. Li, Q. Liao, Z. Liu. Identification of a five-gene signature and establishment of a prognostic nomogram to predict progression-free interval of papillary thyroid carcinoma. Frontiers in Endocrinology, 2019. [DOI | PubMed]
43. M. Xing, B. R. Haugen, M. Schlumberger. Progress in molecular-based management of differentiated thyroid cancer. Lancet, 2013. [DOI | PubMed]
44. A. Czarniecka, M. Oczko-Wojciechowska, M. Barczyński. BRAF V600E mutation in prognostication of papillary thyroid cancer (PTC) recurrence. Gland Surgery, 2016. [DOI | PubMed]
45. F. Wang, S. Zhao, X. Shen. BRAF V600E confers male sex disease-specific mortality risk in patients with papillary thyroid cancer. Journal of Clinical Oncology, 2018. [DOI | PubMed]
46. J. Shi, P. Wu, L. Sheng, W. Sun, H. Zhang. Ferroptosis-related gene signature predicts the prognosis of papillary thyroid carcinoma. Cancer Cell International, 2021. [DOI | PubMed]
47. M. Xin, X. W. Dong, X. L. Guo. Role of the interaction between galectin-3 and cell adhesion molecules in cancer metastasis. Biomedicine & Pharmacotherapy, 2015. [DOI | PubMed]
48. H. Chen, H. Liu, G. Qing. Targeting oncogenic Myc as a strategy for cancer treatment. Signal Transduction and Targeted Therapy, 2018. [DOI | PubMed]
49. S. Morla. Glycosaminoglycans and glycosaminoglycan mimetics in cancer and inflammation. International Journal of Molecular Sciences, 2019. [DOI | PubMed]
50. C. Bi, G. Xiao, C. Liu. Molecular immune mechanism of intestinal microbiota and their metabolites in the occurrence and development of liver cancer. Frontiers in Cell and Development Biology, 2021. [DOI | PubMed]
51. S. Terrisse, L. Derosa, V. Iebba. Intestinal microbiota influences clinical outcome and side effects of early breast cancer treatment. Cell Death and Differentiation, 2021. [DOI | PubMed]
52. J. Knezevic, C. Starchl, A. Tmava Berisha, K. Amrein. Thyroid-gut-axis: how does the microbiota influence thyroid function?. Nutrients, 2020. [DOI | PubMed]