Skip to main content

NKD2 is correlated with the occurrence, progression and prognosis of thyroid carcinoma

Abstract

Background

Thyroid carcinoma (THCA) is the most prevalent type of tumor in endocrine system. NKD2 has been increasingly evidenced to play crucial roles in many cancers, except for THCA. We herein aimed to explore the potential role of NKD2 in THCA.

Methods

Totally 502 THCA patient data were downloaded from TCGA (The Cancer Genome Atlas) database. Overall survival was estimated by Kaplan–Meier method. Gene set enrichment analysis was conducted to obtain significant functional pathways. Wilcoxon rank sum test was used to determine the NKD2 expression differences among various groups. The NKD2 expression was validated in cell lines and tissue microarray.

Results

Significantly higher NKD2 expression was observed in THCA samples compared with adjacent samples, which were successfully verified in cell lines and tissue microarray. Moreover, NKD2 expression gradually elevated along with the increase of TNM Stage, and NKD2 expression was significantly higher in elder THCA patients compared with young patients. NKD2 highly expressed THCA patients had worse prognosis compared with NKD2 low-expressed patients. Furthermore, 53 pathways were significantly activated in the high NKD2 expression patients compared with low NKD2 expression THCA patients.

Conclusions

In summary, high NKD2 expression was probably related to the progression and poor prognosis of THCA. NKD2 is a promising prognostic biomarker and pathogenic target of THCA.

Background

Thyroid carcinoma (THCA) is the most prevalent type of tumor in endocrine system, accounting for approximately 4% of the newly occurred tumors [1, 2]. Over the past decades, growing THCA incidence rate has been reported in China and many other countries [3,4,5], which has caught our attention. Moreover, the region, race, and gender usually result in differential THCA incidence [6,7,8]. Generally, according to histopathology and clinical characteristics, there are four subtypes of THCA, including papillary, follicular, medullary, and undifferentiated THCA [9]. Currently, most THCA patients undergo conventional treatments like surgery, thyroid hormone replacement, or radioiodine therapy [10], but these methods are still less effective for patients with some specific THCA subtypes [11]. Furthermore, regarding early detection and diagnosis of THCA, the limitation of core needle biopsy is lack of widely recognized diagnostic criteria now [12], leading to a limited diagnostic accuracy or delayed diagnosis of THCA. Not only that, most poor prognosis or death of THCA patients result from the recurrence of tumor, metastasis, or radioiodine resistance [10]. Accordingly, it is imperative to find novel sensitively diagnostic and prognostic biomarkers for THCA.

Naked cuticle homolog 2 (NKD2) belongs to the NKD family, and it is located on chromosome 5p15.3 [13], encoding a 451-amino acid protein [14]. NKD2 or NKD2 protein mainly works by acting as a suppressor of WNT signaling pathway [15, 16]. Aberrant expression and epigenetic alterations of NKD2 have been increasingly evidenced to play crucial roles in many cancers. For example, it has been demonstrated that in many solid tumors, aberrant NKD2 expression resulted from promoter hypermethylation [17,18,19]. Besides, NKD2 has been frequently methylated in human breast cancer and gastric cancer [20]. In addition, it has been documented that NKD2 played important roles in the development of several tumors, like glioma, lung cancer, osteosarcoma, and gastric cancer [21]. Decreased NKD2 expression has been reported to be associated with the poor outcome in cytogenetically normal acute myeloid leukemia patients [22]. NKD2 has been indicated to be a potential prognostic biomarker and therapeutic target in hepatocellular carcinoma [15]. NKD2 could suppress the cell growth and tumor metastasis of osteosarcoma [23]. However, as far as we know, NKD2 has not been systematically explored in THCA. Therefore, it would provide more reference information for NKD2 and THCA research if we could preliminarily explore the role of NKD2 in THCA.

In this work, through comprehensive bioinformatics analyses of the THCA-related data obtained from The Cancer Genome Atlas (TCGA) database and further experimental validation, we aimed to study the possible roles of NKD2 in the occurrence, progression, or prognosis of THCA patients, in order to provide more optional choices for THCA clinical treatment strategies.

Materials and methods

Patient data

We have downloaded the THCA patient data from The Cancer Genome Atlas (TCGA) (https://tcga-data.nci.nih.gov/tcga/) database. Totally 502 THCA patients’ mRNA profiles and clinical data were downloaded, which included 510 THCA tissues and 58 paired adjacent tissues. Among which, one patient’s data were excluded as his/her survival information was incomplete. The the other 501 patients’ survival information was recorded completely and their data were used for further analyses. Taking the median of NKD2 expression level as the cut-off value, all THCA samples were divided into two groups, high-NKD2 expression group and low-NKD2 expression group. The clinical information of the 501 THCA patients is summarized in Table 1.

Table1 Clinicopathological characteristics of THCA patients from TCGA database

Cell lines and culture

Totally three cell lines were used in our study, including normal human thyroid follicular epithelial cell line NTHY-ORI3-1 (purchased from BeNa Culture Collection, Beijing, China), THCA cell line SW579 (purchased from Cell Resource Center, Beijing, China), and human anaplastic thyroid carcinoma 8305C (purchased from National Collection of Authenticated Cell Cultures, Shanghai, China). Then, NTHY-ORI3-1, SW579, and 8305C were cultured in RPMI-1640 medium (GIBCO, 31870082), L-15 medium (GIBCO, 41300039), and Eagle’s minimal essential medium (EMEM) complete medium (Invitrogen, 11090081), respectively. All cell lines were cultured in 90% medium and 10% FBS at 37 °C with an atmosphere of 5% CO2.

Gene set enrichment analysis (GSEA)

GSEA analysis was conducted in GSEA software (version: #4.0), on the basis of c2.cp.kegg.v7.0.symbols in Molecular Signatures Database (MSigDB). The P value was adjusted according to Benjamini and Hochberg (BH) method, and the pathways with adjust P value < 0.05 were taken as significantly enriched KEGG pathways.

qRT-PCR

The total RNA was extracted with TRIzol (Tiangen, DP431, Beijing, China), which was measured with a spectrophotometer (Thermo, New York, USA). The reagents used in qPCR experiments: reverse transcription kit (Tiangen Biochemical, KR118, Beijing, China) and SYBR detection reagent (Tiangen, FP205, Beijing, China). The following thermocycling conditions were applied for qPCR: denaturation at 95˚C for 15 min; 40 cycles of 95 °C for 10 s, 60 °C for 20 s, 72 °C for 20 s. The internal reference was β-actin. The primer sequences are listed in Table 2. There were three repeats per sample. The relative mRNA expression was calculated based on formula 2CT.

Table 2 Primer sequences for RT-PCR

Western blot analysis

The methods we used in western blot analysis were consistent with the previous common methods [24]. The reagents we used were shown as below: BCA protein concentration detection kit (Solarbio, cat#PC0020, Beijing, China), microplate reader (SpectraMax, M5, USA), primary antibody NaKed 2 polycIonal AnTibody (Abcam, ab153717, 1:300, UK), internal reference GAPDH (Bioss, cat#bs-2188R, 1:1000, Beijing, China), and secondary antibody Goat anti-rabbit IgG-HRP (Bioss, cat#bs-0295G-HRP, 1:4000, Beijing, China). Software Gelpro32 was used for gray value analyses.

Tissue microarray and immunohistochemical analyses

The THCA tissue microarray (HThyP120CS02, OUTDO Biotech, Shanghai, China) contained 4 normal thyroid tissue samples, and 116 cases of papillary thyroid carcinoma and the paired paracancerous tissues. The clinical information of the tissue microarray is shown in Additional file 1: Table S1 (all detailed information could be available at https://www.superchip.com.cn/biology/tissue.html). The sections were scanned using Pannoramic 250FLASH (3DHISTECH, Hungary) and captured with Aipathwell (Servicebio, China), then the integral optic density (IOD) was calculated. The primary antibody NaKed 2 polycIonal AnTibody (ab33654A36, 1:30, Abcam, UK), the secondary antibody Goat Anti-rabbit IgG/HRP (SE134, Solarbio, 1:150, Beijing, China), and histochemistry kit DAB (DAKO, K5007, USA) were used in our experiment.

Statistical analyses

The NKD2 expression differences between THCA samples, adjacent samples, and samples with various clinicopathological characteristics were determined by Wilcoxon rank sum test. The multivariate Cox regression analysis was done to determine the independent prognostic indicators for THCA. All statistical analyses were conducted using R software v3.5.2. P value < 0.05 was considered as statistically significant difference.

Results

Higher expression of NKD2 in THCA patients

Firstly, the expression of NKD2 in THCA samples and adjacent samples were compared to explore the potential correlation between NKD2 expression and THCA onset. The results showed that compared with in adjacent samples, the expression of NKD2 in paired THCA samples (N = 58) was significantly higher (P = 0.00015) (Fig. 1A). Moreover, significantly higher NKD2 expression was also observed in all THCA samples (N = 510) compared with adjacent samples (N = 58) (P = 3.4e–08) (Fig. 1B). Our results indicated that high expression of NKD2 was related to the occurrence of THCA.

Fig. 1
figure 1

The NKD2 expression levels in THCA samples and adjacent samples. A The NKD2 expression levels in paired THCA and adjacent samples. B NKD2 expression levels in all THCA and adjacent samples

Correlation between NKD2 expression and clinicopathologic features

Subsequently, the correlation between NKD2 expression and various corresponding THCA clinicopathological characteristics was analyzed based on Wilcoxon rank sum test. We found that the expression of NKD2 gradually increased along with the increase of TNM Stage, besides, there were significant differences between Stage I vs. Stage III, Stage II vs. Stage III, Stage II vs. Stage IV, and Stage I vs. Stage IV (P value < 0.05, Fig. 2A). Moreover, elevated NKD2 expression was significantly correlated with Age of THCA patients. NKD2 expression was significantly higher in elder THCA patients compared with young patients (P value < 0.05, Fig. 2C). However, there was no significant correlation between NKD2 expression and sex (P value > 0.05, Fig. 2B).

Fig. 2
figure 2

The correlation between NKD2 expression and clinicopathological characteristics, prognosis of THCA patients. A NKD2 expression levels in different TNM stages of THCA samples. B NKD2 expression levels in male and female THCA patients. C NKD2 expression levels in different ages. D Kaplan–Meier survival curves of high and low NKD2 expression THCA patients. P value was calculated based on log-rank test. X-axis: time (days); Y-axis: survival probability. E The results of multivariate Cox regression analysis

THCA patients with high NKD2 expression had poor prognosis

To further study the prognostic value of NKD2 in THCA patients, the survival analysis was conducted on high and low NKD2 expression THCA patients, based on the data in TCGA database. Comparing with patients with low NKD2 expression, patients with high NKD2 expression had a poorer overall survival (p = 0.036, HR = 0.32, 95% CI 0.12–0.85) (Fig. 2D). Additionally, a multivariate Cox regression analysis was conducted including NKD2 expression and sex, which suggested that NKD2 was still significantly associated with the prognosis of THCA patients (Fig. 2E).

NKD2-related results of GSEA

Furthermore, signaling pathways significantly enriched in NKD2 high-expressed THCA patients were identified via GSEA enrichment analysis, compared with NKD2 low-expressed patients. The KEGG pathways with adjusted P value < 0.05 were screened as significantly enriched pathways. There were totally 53 pathways significantly enriched in the NKD2 high-expressed group, and the top 9 most significant pathways included T cell receptor signaling pathway, B cell receptor signaling pathway, thyroid cancer, viral myocarditis, antigen processing and presentation, cytosolic DNA sensing pathway, pathways in cancer, acute myeloid leukemia, and JAK/STAT signaling pathway (Fig. 3A–I). The results suggested that thyroid cancer pathway and several other cancer-related pathways were significantly activated in NKD2 highly expressed THCA patients. The detailed results of GSEA are shown in Additional file 2: Table S2.

Fig. 3
figure 3

Significantly activated KEGG pathways in NKD2 high-expressed THCA patients, based on GSEA. A T cell receptor signaling pathway, (B) B cell receptor signaling pathway, (C) thyroid cancer, (D) viral myocarditis, (E) antigen processing and presentation, (F) cytosolic DNA sensing pathway, (G) pathways in cancer, (H) acute myeloid leukemia, (I) jak/stat signaling pathway. The top 9 significantly enriched pathways (P < 0.05)

High NKD2 expression validated in THCA cell lines and tissue microarray

Subsequently, we have further validated our results in cell lines and tissue microarray. Firstly, in three cell lines, the relative NKD2 mRNA expression was significantly higher in THCA cell lines compared with in normal cells (Fig. 4A), and the tendency of protein expression was consistent (Fig. 4B–C, Additional file 3). Moreover, the results of IHC analyses suggested that NKD2 expression was still significantly higher in THCA tissue samples compared with adjacent samples (Fig. 5A–B). Our findings have been successfully validated in cell lines and clinical samples.

Fig. 4
figure 4

The NKD2 expression in cell lines. A NKD2 mRNA was upregulated in THCA cell lines (* P < 0.05, vs. SW579; ## P < 0.001, vs. 8305C). B, C The NKD2 protein expression levels in THCA cell lines (*** P < 0.001, vs. SW579; ### P < 0.001, vs. 8305C)

Fig. 5
figure 5

The expression of A NKD2 in B THCA tissue microarray.

Discussion

Over the past decades, increasing incidence of THCA has been widely reported, especially in females [25]. Despite significant advances have been achieved in THCA therapy, effective treatment strategies for some specific types of THCA are still lacking. In the present study, we have explored the potential role of NKD2 in the progression and prognosis of THCA. High NKD2 expression was probably related to the progression and occurrence of THCA. Moreover, NKD2 highly expressed THCA patients had worse prognosis compared with patients with low NKD2 expression.

The role of NKD2 in many tumors has been reported [26, 27], except for THCA. Consequently, we aimed to better understand the role of NKD2 in THCA, on the basis of data in TCGA database. Firstly, we found that significantly higher NKD2 expression was observed in THCA samples (including paired THCA vs. adjacent samples and all THCA vs. adjacent samples), indicating the potential pathogenic role of NKD2 in THCA. Our results were also successfully validated in cell lines and clinical samples, and were compatible with some previous reports directly or indirectly. NKD2 has been evidenced to serve as a suppressor of WNT signaling pathway. Additionally, the role of NKD2 involving development, progression, or prognosis in several tumors has been demonstrated, like glioma, lung cancer, osteosarcoma, and gastric cancer [19, 21]. Despite we have explored the role of NKD2 in THCA for the first time, only a small part of the NKD2 function was found in our study. NKD2 was previously reported to be regulated by promoter region methylation [18, 20] and it was the target of some lncRNAs or miRNAs in the regulation of several cancers [26, 28], which should be further studied in our future research. Furthermore, elevated NKD2 expression was associated with the increase of TNM Stage, which suggested indirectly that NKD2 might contribute to the progression of THCA, as a pathogenic factor. On this point, NKD2-related underlying mechanisms in THCA deserved further investigation, and its detailed function is probably conducive to studying potential therapy strategies targeting NKD2. Moreover, there was no significant correlation between NKD2 expression and gender, although it has been widely reported that a higher incidence of THCA was observed in females than in males [29]. While NKD2 expression seems to be not correlated with gender, which might result from limited samples in this research. Additionally, THCA patients with high NKD2 expression had a worse prognosis compared with patients with low NKD2 expression. The prognostic value of NKD2 has been demonstrated in some other tumors, for example, hepatocellular carcinoma [15], ovarian cancer [30], acute myeloid leukemia [22]. More recently, Wu et al. have documented the prognostic value of several WNT signaling pathway molecules in colorectal cancer (CRC), including NKD2, and the inhibition of NKD2 could inhibit the tumor cell growth of CRC [31]. Higher NKD2 expression was significantly correlated with poorer survival of CRC [31], thereby their results partly supported our findings in THCA. Collectively, NKD2 is a poor prognostic factor for THCA, which is firstly reported in our research.

Furthermore, the GSEA analysis was conducted to better understand the potential function of NKD2 in THCA. Totally, 53 pathways were significantly activated in the NKD2 high-expressed group compared to the NKD2 low-expressed group. Notably, 9 various cancer-related pathways, including Thyroid cancer pathway, were significantly enriched, which implied that NKD2 was actually involved in cancer. Moreover, some of the enriched pathways have been recently evidenced to involve in the anticancer activities of THCA, for instance, JAK/STAT3 signaling pathway [32, 33]. Not only that, several known pathways that were involved in oncogenic mechanisms in thyroid cancer initiation or progression were also found to be significantly activated, including VEGF signaling pathway [34], Toll-like receptor signaling pathway [35], MAPK signaling pathway [36]. Nonetheless, the precise role of NKD2 in THCA still needs to be clarified.

Although the pathogenic role and prognostic value of NKD2 in THCA have been demonstrated in our study for the first time, there were still some limitations in our work. The THCA-related data in TCGA database are limited, and more public data will be helpful to strengthen our results. Secondly, we have preliminarily explored the role of NKD2 via in silico analysis and wet-lab validation, but deepening underlying mechanisms regarding NKD2 are not well investigated herein.

Conclusions

In conclusion, our study explored the important role of NKD2 in THCA patients for the first time. On the basis of integrated bioinformatics analyses of THCA-related data in TCGA database, we also conducted further validation experiments in cell lines and clinical tissues to further confirm that high NKD2 expression was correlated with the progression and prognosis of THCA. Collectively, our findings indicate that NKD2 is a promising prognostic biomarker and pathogenic target of THCA.

Availability of data and materials

The datasets generated and analyzed for this study can be found in the TCGA database (https://tcga-data.nci.nih.gov/tcga/).

Abbreviations

THCA:

Thyroid carcinoma

NKD2:

Naked cuticle homolog 2

TCGA:

The Cancer Genome Atlas

GSEA:

Gene set enrichment analysis

References

  1. Salazar-Vega J, Ortiz-Prado E, Solis-Pazmino P, Gomez-Barreno L, Simbana-Rivera K, Henriquez-Trujillo AR, et al. Thyroid cancer in ecuador, a 16 years population-based analysis (2001–2016). BMC Cancer. 2019;19(1):294.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Paschou SA, Vryonidou A, Goulis DG. Thyroid nodules: Alpha guide to assessment, treatment and follow-up. Maturitas. 2017;96:1–9.

    Article  PubMed  Google Scholar 

  3. Kilfoy BA, Zheng T, Holford TR, Han X, Ward MH, Sjodin A, et al. International patterns and trends in thyroid cancer incidence, 1973–2002. Cancer Causes Control. 2009;20(5):525–31.

    Article  PubMed  Google Scholar 

  4. Enewold L, Zhu K, Ron E, Marrogi AJ, Stojadinovic A, Peoples GE, et al. Rising thyroid cancer incidence in the United States by demographic and tumor characteristics, 1980–2005. Cancer Epidemiol Biomarkers Prev. 2009;18(3):784–91.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Huk DJ, Ashtekar A, Magner A, La Perle K, Kirschner LS. Deletion of Rap1b, but not Rap1a or Epac1, reduces protein kinase a-mediated thyroid cancer. Thyroid. 2018;28(9):1153–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ward EM, Sherman RL, Henley SJ, Jemal A, Siegel DA, Feuer EJ, et al. Annual report to the nation on the status of cancer, featuring cancer in men and women age 20–49 years. J Natl Cancer Inst. 2019;111(12):1279–97.

    Article  PubMed  PubMed Central  Google Scholar 

  7. La Vecchia C, Malvezzi M, Bosetti C, Garavello W, Bertuccio P, Levi F, et al. Thyroid cancer mortality and incidence: a global overview. Int J Cancer. 2015;136(9):2187–95.

    Article  PubMed  Google Scholar 

  8. Colonna M, Uhry Z, Guizard AV, Delafosse P, Schvartz C, Belot A, et al. Recent trends in incidence, geographical distribution, and survival of papillary thyroid cancer in France. Cancer Epidemiol. 2015;39(4):511–8.

    Article  CAS  PubMed  Google Scholar 

  9. Lebastchi AH, Callender GG. Thyroid cancer. Curr Probl Cancer. 2014;38(2):48–74.

    Article  PubMed  Google Scholar 

  10. Heiden KB, Williamson AJ, Doscas ME, Ye J, Wang Y, Liu D, et al. The sonic hedgehog signaling pathway maintains the cancer stem cell self-renewal of anaplastic thyroid cancer by inducing snail expression. J Clin Endocrinol Metab. 2014;99(11):E2178–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Conzo G, Avenia N, Bellastella G, Candela G, de Bellis A, Esposito K, et al. The role of surgery in the current management of differentiated thyroid cancer. Endocrine. 2014;47(2):380–8.

    Article  CAS  PubMed  Google Scholar 

  12. Yoon JH, Kim EK, Kwak JY, Moon HJ. Effectiveness and limitations of core needle biopsy in the diagnosis of thyroid nodules: review of current literature. J Pathol Transl Med. 2015;49(3):230–5.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lu Y, Yu Y, Zhu Z, Xu H, Ji J, Bu L, et al. Identification of a new target region by loss of heterozygosity at 5p1533 in sporadic gastric carcinomas: genotype and phenotype related. Cancer Lett. 2005;224(2):329–37.

    Article  CAS  PubMed  Google Scholar 

  14. Katoh M. Molecular cloning, gene structure, and expression analyses of NKD1 and NKD2. Int J Oncol. 2001;19(5):963–9.

    CAS  PubMed  Google Scholar 

  15. Wang D, Zhang S, Chen Y, Hu B, Lu C. Low expression of NKD2 is associated with enhanced cell proliferation and poor prognosis in human hepatocellular carcinoma. Hum Pathol. 2018;72:80–90.

    Article  CAS  PubMed  Google Scholar 

  16. Hu T, Krezel AM, Li C, Coffey RJ. Structural studies of human Naked2: a biologically active intrinsically unstructured protein. Biochem Biophys Res Commun. 2006;350(4):911–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Gotze S, Wolter M, Reifenberger G, Muller O, Sievers S. Frequent promoter hypermethylation of Wnt pathway inhibitor genes in malignant astrocytic gliomas. Int J Cancer. 2010;126(11):2584–93.

    PubMed  Google Scholar 

  18. Dong Y, Cao B, Zhang M, Han W, Herman JG, Fuks F, et al. Epigenetic silencing of NKD2, a major component of Wnt signaling, promotes breast cancer growth. Oncotarget. 2015;6(26):22126–38.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Jia Y, Cao B, Yang Y, Linghu E, Zhan Q, Lu Y, et al. Silencing NKD2 by promoter region hypermethylation promotes gastric cancer invasion and metastasis by up-regulating SOX18 in human gastric cancer. Oncotarget. 2015;6(32):33470–85.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Cao B, Yang W, Jin Y, Zhang M, He T, Zhan Q, et al. Silencing NKD2 by promoter region hypermethylation promotes esophageal cancer progression by activating wnt signaling. J Thorac Oncol. 2016;11(11):1912–26.

    Article  PubMed  Google Scholar 

  21. Tabatabai R, Linhares Y, Bolos D, Mita M, Mita A. Targeting the wnt pathway in cancer: a review of novel therapeutics. Target Oncol. 2017;12(5):623–41.

    Article  PubMed  Google Scholar 

  22. Li XX, Zhou JD, Zhang TJ, Yang L, Wen XM, Ma JC, et al. Epigenetic dysregulation of NKD2 is a valuable predictor assessing treatment outcome in acute myeloid leukemia. J Cancer. 2017;8(3):460–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zhao S, Kurenbekova L, Gao Y, Roos A, Creighton CJ, Rao P, et al. NKD2, a negative regulator of Wnt signaling, suppresses tumor growth and metastasis in osteosarcoma. Oncogene. 2015;34(39):5069–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Wang YX, Liu L, Zeng QX, Fan TY, Jiang JD, Deng HB, et al. Synthesis and Identification of novel berberine derivatives as potent inhibitors against TNF-alpha-induced NF-kappaB activation. Molecules. 2017;22(8):1257.

    Article  PubMed Central  Google Scholar 

  25. Roman BR, Morris LG, Davies L. The thyroid cancer epidemic, 2017 perspective. Curr Opin Endocrinol Diabetes Obes. 2017;24(5):332–6.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Tang H, Dou Y, Meng Y, Lu Q, Liang L, Luo Y. LINC00538 promotes the progression of colon cancer through inhibiting NKD2 expression. J BUON. 2020;25(6):2657–64.

    PubMed  Google Scholar 

  27. Wang Y, Dong T, Wang P, Li S, Wu G, Zhou J, et al. LINC00922 regulates epithelial-mesenchymal transition, invasive and migratory capacities in breast cancer through promoting NKD2 methylation. Cell Signal. 2021;77: 109808.

    Article  CAS  PubMed  Google Scholar 

  28. Li Z, Li Y, Wang N, Yang L, Zhao W, Zeng X. miR-130b targets NKD2 and regulates the Wnt signaling to promote proliferation and inhibit apoptosis in osteosarcoma cells. Biochem Biophys Res Commun. 2016;471(4):479–85.

    Article  CAS  PubMed  Google Scholar 

  29. Zhang D, Tang J, Kong D, Cui Q, Wang K, Gong Y, et al. Impact of gender and age on the prognosis of differentiated thyroid carcinoma: a retrospective analysis based on SEER. Horm Cancer. 2018;9(5):361–70.

    Article  CAS  PubMed  Google Scholar 

  30. Wei W, Zheng L, Gao Y, He M, Yang F. Expression and prognostic significance of NKD2 in ovarian cancer. Jpn J Clin Oncol. 2021;51(3):459–68.

    Article  PubMed  Google Scholar 

  31. Wu YF, Wang CY, Tang WC, Lee YC, Ta HDK, Lin LC, et al. Expression profile and prognostic value of wnt signaling pathway molecules in colorectal cancer. Biomedicines. 2021;9(10):1331.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Khan AQ, Ahmed EI, Elareer N, Fathima H, Prabhu KS, Siveen KS, et al. Curcumin-mediated apoptotic cell death in papillary thyroid cancer and cancer stem-like cells through targeting of the JAK/STAT3 signaling pathway. Int J Mol Sci. 2020;21(2):438.

    Article  CAS  PubMed Central  Google Scholar 

  33. Bian P, Hu W, Liu C, Li L. Resveratrol potentiates the anti-tumor effects of rapamycin in papillary thyroid cancer: PI3K/AKT/mTOR pathway involved. Arch Biochem Biophys. 2020;689: 108461.

    Article  CAS  PubMed  Google Scholar 

  34. Du Y, Zhu J, Chu BF, Yang YP, Zhang SL. MiR-548c-3p suppressed the progression of papillary thyroid carcinoma via inhibition of the HIF1alpha-mediated VEGF signaling pathway. Eur Rev Med Pharmacol Sci. 2019;23(15):6570–8.

    CAS  PubMed  Google Scholar 

  35. Peyret V, Nazar M, Martin M, Quintar AA, Fernandez EA, Geysels RC, et al. Functional toll-like receptor 4 overexpression in papillary thyroid cancer by MAPK/ERK-induced ETS1 transcriptional activity. Mol Cancer Res. 2018;16(5):833–45.

    Article  CAS  PubMed  Google Scholar 

  36. Lu ZW, Wen D, Wei WJ, Han LT, Xiang J, Wang YL, et al. Silencing of PPM1D inhibits cell proliferation and invasion through the p38 MAPK and p53 signaling pathway in papillary thyroid carcinoma. Oncol Rep. 2020;43(3):783–94.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Contributions

YG contributed to the study conception and design. Data collection and analysis were performed by all authors. Validation was conducted by YG and WW. The first draft of the manuscript was written by YG and DG and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Rende Guo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Clinical features of thyroid cancer patients.

Additional file 2:

Table S2. GSEA result.

Additional file 3.

The original results of western blot.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gao, Y., Wang, Y. & Guo, R. NKD2 is correlated with the occurrence, progression and prognosis of thyroid carcinoma. Eur J Med Res 27, 235 (2022). https://doi.org/10.1186/s40001-022-00853-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-022-00853-2

Keywords