- Research
- Open access
- Published:
NKD2 is correlated with the occurrence, progression and prognosis of thyroid carcinoma
European Journal of Medical Research volume 27, Article number: 235 (2022)
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.
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 2−△△CT.
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.
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).
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.
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.
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
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.
Paschou SA, Vryonidou A, Goulis DG. Thyroid nodules: Alpha guide to assessment, treatment and follow-up. Maturitas. 2017;96:1–9.
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.
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.
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.
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.
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.
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.
Lebastchi AH, Callender GG. Thyroid cancer. Curr Probl Cancer. 2014;38(2):48–74.
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.
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.
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.
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.
Katoh M. Molecular cloning, gene structure, and expression analyses of NKD1 and NKD2. Int J Oncol. 2001;19(5):963–9.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Roman BR, Morris LG, Davies L. The thyroid cancer epidemic, 2017 perspective. Curr Opin Endocrinol Diabetes Obes. 2017;24(5):332–6.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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
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
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40001-022-00853-2