Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 11 November 2025

Sec. Thyroid Endocrinology

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1689674

This article is part of the Research TopicMolecular Characterization of Thyroid Lesions in the Era of “Next Generation” Techniques: Volume IIIView all 6 articles

A Nomogram for Predicting Lung Metastasis in Papillary Thyroid Cancer Patients Aged Less Than 55 Years

Huiyun Yang&#x;Huiyun Yang1†Qiuqin Qian,&#x;Qiuqin Qian1,2†Yanling LiYanling Li1Wenjie PanWenjie Pan1Jing PengJing Peng1Zhenyu ZouZhenyu Zou1Haiqing ZhuHaiqing Zhu1Feng Shi*Feng Shi1*
  • 1Department of Nuclear Medicine, The Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Hunan Cancer Hospital, Changsha, China
  • 2Department of Radiology, Second Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, Hunan, China

Highlights: A nomogram predicting lung metastasis in patients with papillary thyroid cancer was developed and validated. Our study focused on patients under 55 years old, revealing the genetic factors linked to lung metastasis in this younger group. Young patients with the BRAF V600E mutation are less likely to develop lung metastasis.

Objective: To develop and validate a nomogram for predicting the risk of lung metastasis in patients under 55 years old with papillary thyroid cancer (PTC).

Methods: A total of 243 patients under 55 years old with PTC were retrospectively collected from January 2017 to June 2020 and randomly divided into a training group (n = 170) and a validation group (n = 73) in a 7:3 ratio. Genetic testing data and clinical information were compiled, and univariate and multivariate binary logistic regression analyses were performed. Based on the results, a nomogram predicting the risk of lung metastasis was constructed using the training cohort. The nomogram’s performance was assessed using calibration curves, decision curve analysis (DCA), the concordance index (C-index), and the area under the receiver operating characteristic (ROC) curve (AUC) in both the training and validation groups.

Results: T stage, unilateral thyroid involvement, TERT mutation, and BRAF mutation were identified as independent prognostic factors and were used to construct the nomogram. The C-index of the model was 0.89 in the training group and 0.88 in the validation group. The AUC, DCA, and calibration curves demonstrated favorable predictive accuracy. Using a cut-off value of 0.229, the nomogram achieved a sensitivity of 0.829 and a specificity of 0.793.

Conclusion: A nomogram with strong predictive performance has been successfully developed and validated, which may assist clinicians in estimating the risk of lung metastasis in young patients with PTC.

Introduction

Papillary thyroid cancer (PTC) is the most common endocrine malignancy of the head and neck. Differentiated thyroid cancer accounts for over 90% of thyroid cancer cases. The standard treatment typically involves radioactive iodine-131 therapy or active surveillance following surgery. Although most patients have a favorable prognosis, 10-15% still develop distant metastases, and the 10-year overall survival rate in this subgroup is only about 40% (1, 2). Therefore, identifying factors associated with distant metastasis is critical for the early recognition of high-risk patients and the timely initiation of individualized treatment strategies.

Previous studies have investigated risk factors for distant metastasis in differentiated thyroid cancer. Khan U et al. (3) reported that tumor size ≥2 cm, nodal metastasis, and histological subtype are significantly associated with distant spread. Other researchers have identified tumor size, age, surgical intervention, N stage, T stage, and pathological type as independent predictors of distant metastasis in female patients with differentiated thyroid cancer (4). In recent years, the role of genetic mutations in metastatic risk has gained increasing attention. Gene fusions involving RET, ALK, or NTRK1, among others, have been shown to correlate with metastatic potential in thyroid cancer (5). However, few studies have integrated clinical features with genetic testing data to predict the risk of distant metastasis in these patients.

The lungs are the most common site of distant metastasis in PTC, accounting for approximately 45% of all metastatic cases (6). Pulmonary metastasis is also a leading cause of mortality in patients with PTC (7). Hence, early identification of individuals at risk for lung metastasis is essential to improve outcomes through targeted and timely intervention. While metastasis in thyroid cancer has been associated with age (4), lung metastases can still occur in younger patients under 55 years of age. Therefore, this study aims to develop a nomogram that combines clinical characteristics with genetic testing results to predict the likelihood of lung metastasis in patients under 55 years old with PTC.

Patients and methods

This study was approved by the Ethics Committee of the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University/Hunan Cancer Hospital (Approval No. KYJJ-2020-195).

Study population

A total of 243 patients under 55 years old with PTC were selected from the Affiliated Cancer Hospital of Xiangya School of Medicine between January 2017 and June 2020 and followed up until July 2023. The inclusion criteria were as follows: (I) postoperative histopathological confirmation of PTC; (II) age <55 years; (III) availability of pathological tissue and complete clinical data; (IV) no evidence of metastasis at the time of initial diagnosis; and (V) development of pulmonary metastases only during the follow-up period. For patients diagnosed solely with “lung metastasis,” at least one of the following criteria was required to confirm the diagnosis (1): pulmonary pathology confirmed independently by two pathologists; (2) radioactive iodine uptake in the lungs; (3) positive findings on positron emission tomography-computed tomography (PET-CT); or (4) chest computed tomography (CT) findings reviewed and confirmed by two radiologists. The exclusion criteria were: (I) the presence of other primary tumors in addition to PTC; (II) metastasis to organs other than the lungs; and (III) other pathological subtypes of thyroid cancer. “unilateral thyroid involvement” is determined based on postoperative histopathological reports. It is defined as: the primary tumor is confined to a single thyroid lobe, with no evidence of tumor involvement in the contralateral lobe on pathological examination.

Genetic testing

All patients underwent whole-genome sequencing using next-generation sequencing (NGS) with the Geneseeq Prime™ 425-gene panel. The specific genes analyzed are listed in Supplementary Table 1. The sequencing procedure was as follows: DNA was extracted from pathological tissue samples. The extracted DNA was then end-repaired, and poly(A) tails were added. Illumina sequencing adapters were ligated to both ends of the DNA fragments, followed by library construction through iterative optimization of conditions. Custom-designed DNA probes targeting oncogenes were synthesized to enrich the library for relevant genes.

Enriched libraries of tumor tissue oncogene exons and negative controls were prepared, pooled based on sequencing throughput, and subjected to high-throughput sequencing using the PE75 kit (Illumina, California, USA). Bioinformatics analysis was subsequently performed to obtain comprehensive mutation data, including single nucleotide variants (SNVs), gene fusions, amplifications, deletions, and insertions in tumor tissues.

Construction of the nomogram

The data were organized, and univariate binary logistic regression analysis was performed to identify variables with P < 0.05. These significant variables were then included in a multivariate binary logistic regression analysis with stepwise selection to identify independent predictors (P < 0.05) for constructing the nomogram. Statistical analyses were conducted using SPSS software and R language.

The Hosmer-Lemeshow test was used to evaluate the goodness of fit between the observed outcomes and the nomogram’s predictions. A calibration curve was generated to compare the predicted probability of lung metastasis with the actual observed outcomes. Decision curve analysis (DCA) was employed to assess the clinical utility and net benefit of the nomogram. The apparent performance of the nomogram is evaluated using the concordance index (C-index), and its optimistic bias is corrected through 1000 bootstrap resamples to obtain the bias-corrected C-index. The C-index and the area under the receiver operating characteristic (ROC) curve (AUC) were calculated to evaluate the nomogram’s predictive accuracy. Using SPSS software, the probability of lung metastasis was calculated for each patient in the training group. The optimal cut-off value, also known as Youden index, was determined using MedCalc software (8).

Statistical analysis

Statistical analyses were performed using SPSS software and R version 4.3.2 (http://www.R-project.org/), with P < 0.05 considered statistically significant. MedCalc software (Ostend, Belgium) was used to determine the optimal cut-off value.

Results

Patients’ characteristics

From January 2017 to June 2020, a total of 243 patients under 55 years old with PTC were enrolled in the study and randomly assigned to a training group (n = 170) and a validation group (n = 73) in a 7:3 ratio. The clinical characteristics of the patients are summarized in Table 1. The median age was 36 years in the training group and 38 years in the validation group. Lung metastases were observed in 35 patients in the training group and 15 patients in the validation group. A total of 151 patients had unilateral thyroid lobe involvement. The mean maximum diameter of primary thyroid lesions was 2.3 cm in the training group and 1.8 cm in the validation group. Extraglandular invasion was observed in 51.2% of patients in the training group and 46.6% in the validation group.

Table 1
www.frontiersin.org

Table 1. The characteristics of patients in the training and validation cohorts.

All patients underwent whole-genome sequencing using NGS. BRAF mutations were detected in 147 of 243 patients (60.4%), representing the most frequently observed mutation, followed by RET mutations (35/243) and TERT mutations (16/243). Detailed results are presented in Table 1.

Identifying predictors

Univariate binary logistic regression analysis identified age, T stage, N stage, unilateral thyroid involvement, maximum tumor diameter, extra-glandular invasion, and RET, TERT, and BRAF gene mutations as factors significantly associated with lung metastasis in thyroid cancer. In the multivariate analysis, T stage (OR = 5.564, 95% CI: 1.769-17.502, P = 0.003), unilateral thyroid involvement (OR = 8.706, 95% CI: 3.002-25.248, P < 0.001), TERT mutation (OR = 7.585, 95% CI: 1.245-46.199, P = 0.028), and BRAF mutation (OR = 0.070, 95% CI: 0.022-0.225, P < 0.001) were confirmed as independent predictors of lung metastasis and were used to construct the nomogram in R (Table 2). The Hosmer-Lemeshow test yielded a P-value of 0.986, indicating a good model fit.

Table 2
www.frontiersin.org

Table 2. Univariate and multivariate binary logistic regression analyses in the training cohort.

Construction and evaluation of nomogram

Based on the multivariate regression analysis results, each factor was assigned a specific score. The detailed regression coefficients, odds ratios, and the score for each predictor in the final nomogram are provided in Supplementary Table 3. According to Supplementary Table 3, for young patients suspected of lung metastasis, we can assign scores based on BRAF mutation, TERT promoter mutation, unilateral thyroid involvement, and T stage. These scores are then summed to obtain the “total points” in our nomogram. A vertical line is drawn from the total points to the corresponding risk of lung metastasis, thereby determining the patient’s probability of developing lung metastasis (Figure 1).

Figure 1
Nomogram for predicting lung metastasis risk using various factors. Points range from zero to one hundred. Factors include T stage (T1/T2, T3/T4), unilateral thyroid involvement (yes, no), TERT mutation (no, yes), and BRAF(V600E) mutation (yes, no). Total points range from zero to four hundred, correlating with metastasis risk from zero to 0.9.

Figure 1. Nomogram for predicting lung metastases in patients with PTC under 55 years old.

The nomogram demonstrated strong predictive performance, with a C-index of 0.89 in the training cohort and 0.88 in the validation cohort. Calibration curves showed good agreement between predicted and observed outcomes in both the training (Figure 2A) and validation (Figure 2B) groups. This indicates the model’s strong predictive accuracy. Such reliability enables clinicians to trust the model outputs, accurately communicate risks, and formulate appropriate treatment decisions—avoiding both overestimation that leads to overtreatment and underestimation that results in delayed intervention.

Figure 2
Two calibration plots (A and B) compare predicted probability to actual probability. Both plots show three lines: ideal (dashed), apparent (red), and bias-corrected (blue). Plot A demonstrates close alignment of all lines, while in Plot B, the apparent and bias-corrected lines deviate more significantly from the ideal line.

Figure 2. Calibration plots of the nomogram in the training (A) and validation (B) cohorts.

DCA indicated that the nomogram has valuable clinical utility (Figures 3A, B). Specifically, decision analysis indicates that within a broad threshold probability range of approximately 10% to 90%, this predictive map delivers superior net benefits compared to both “all interventions” and “no interventions” strategies, which suggests that when clinicians consider intervention thresholds (such as enhanced follow-up or treatment) within this range, the nomogram can guide decision-making. The area under the ROC curve (AUC) was 0.89 for the training group and 0.88 for the validation group (Figures 4A, B). To assess the added value of the integrated model, we compared its performance against each independent predictor in the training cohort. As shown in Supplementary Figure 1, the nomogram demonstrated superior discriminatory accuracy (Supplementary Figure 1A) and provided a greater net benefit across a wide range of threshold probabilities (Supplementary Figure 1B) compared to any single variable (T stage, unilateral thyroid involvement, TERT mutation, or BRAF mutation).The analysis confirms that the integrated model outperforms any single predictor. The optimal cut-off value, defined by the Youden index, was 0.229, corresponding to a sensitivity of 82.9% and specificity of 79.3% (Figure 5).

Figure 3
Two decision curves labeled A and B show net benefit versus high-risk threshold. Both graphs have three lines representing “Model” in red, “All” in light gray, and “None” in dark gray. The red lines show decreasing net benefit as the threshold increases, with fluctuations closer to zero. The gray lines serve as baselines. Legends are present for clarity.

Figure 3. DCA of the nomogram in the training (A) and validation (B) cohorts.

Figure 4
Two ROC curves labeled A and B. Curve A has an area under the curve of 0.89 with a 95% confidence interval of 0.833 to 0.933. Curve B has an area of 0.88 with a 95% confidence interval of 0.783 to 0.945. Both curves show increasing sensitivity with decreasing specificity.

Figure 4. AUC of the nomogram in the training (A) and validation (B) cohorts.

Figure 5
ROC curve illustrating diagnostic accuracy with the y-axis labeled “Sensitivity” and the x-axis labeled “100-Specificity.” A blue line indicates performance, showing a particular point with sensitivity of 82.9 percent, specificity of 79.3 percent, and a criterion of greater than 0.2286.

Figure 5. The cut-off value determined by the ROC curve in the training cohort.

Discussion

This study developed and validated a nomogram to predict the risk of lung metastasis in patients under 55 years old with PTC by integrating clinical data and genomic sequencing information. T stage, unilateral thyroid involvement, TERT promoter mutation, and BRAF V600E mutation were identified as independent predictors of lung metastasis and incorporated into the nomogram. Patients were stratified into high-risk and low-risk groups using a risk score threshold of 0.229, with the high-risk group showing a significantly higher probability of developing lung metastases.

Several previous studies have developed nomograms to predict the prognosis of thyroid cancer (3, 4, 9), but our study differs from them in key aspects. Wen W et al. (4) constructed a nomogram including age at diagnosis, pathological type, N1 stage, T3–4 stage, thyroidectomy, and tumor size to predict distant metastasis in thyroid cancer. Peng G et al. (10) developed a nomogram based on the number of lymph node metastases, lymph node ratio, and preoperative thyroglobulin to predict survival prognosis in PTC patients. However, these models lack genetic background information. NGS has been widely recommended for molecular classification and predicting recurrence and metastasis in PTC (11). In our study, all patients underwent NGS, and our model combined clinical and genetic data to predict distant metastasis, which could reflect more genetic background information, differing from their studies.

Additionally, previous studies did not specifically focus on age stratification (4, 10). Our study targeted patients younger than 55 years old. According to the 8th edition of the AJCC staging system for thyroid cancer, patients younger than 55 years are staged differently from those 55 or older. Patients under 55 with distant metastases are classified as stage II, while those 55 or older with distant metastases are classified as stage IVB, highlighting age 55 as a critical prognostic threshold (1214). Although younger patients generally have better prognoses, a subset still develops distant metastases and experiences poor outcomes (15, 16). Therefore, identifying risk factors for distant metastases in young patients is essential for enabling timely individualized treatment.

Third, our study specifically focuses on single-organ (lung)-specific metastasis of thyroid cancer, which is rarely addressed in previous research. The lungs are the most common and often the earliest site of metastasis in thyroid cancer patients (1719), particularly in younger patients (15, 20). Kuang HF et al. (21) previously developed a nomogram to predict lung metastasis in young thyroid cancer patients using the SEER database; however, his target population was primarily under 19 years old, and the studied pathologies included both papillary and follicular carcinomas. Importantly, our study excluded patients with metastases to other sites, differentiating it from prior work. Given that lung metastases are typically the most frequent and initial metastatic site in thyroid cancer, this focus enhances clinical relevance.

Using our nomogram, clinicians may identify young patients at high risk for lung metastasis early and adjust clinical management accordingly. For instance, high-risk patients should undergo more frequent and stringent follow-ups, especially those presenting with multiple pulmonary nodules on imaging but without clinical confirmation of metastasis. Postoperative radioactive iodine (I-131) therapy doses typically vary based on individual patient conditions (22); however, for patients identified as high-risk, an increased therapeutic dose may be considered. Furthermore, comprehensive genetic testing is recommended to identify actionable mutations for targeted therapies.

T stage, unilateral thyroid involvement, and TERT and BRAF V600E mutations were all independent predictors of lung metastasis. Advanced T stage (T3/T4; OR: 5.564, 95% CI: 1.769–17.502, P = 0.003) correlated with metastasis consistent with its role in local invasion (21, 23). Similarly, TERT promoter mutations, which promote uncontrolled growth by enhancing telomerase activity (24, 25), demonstrated strong predictive value (OR: 7.585, 95% CI: 1.245–46.199, P = 0.028), further confirming their established association with distant metastasis (2629).

Although extensively implicated in thyroid cancer tumorigenesis and progression (3033), the role of the BRAF V600E mutation in distant metastasis remains complex. Contrasting with its association with aggressive local features (3437), our study identified the BRAF V600E as an independent factor inversely associated with lung metastasis risk (OR: 0.070, 95% CI: 0.022-0.225, P < 0.001). This inverse association is an intriguing phenomenon that finds support in other studies. Some research suggests that BRAF mutation may primarily influence local rather than distant metastasis (27, 3840). For instance, Miguel Melo M et al. (27) reported that in PTC, compared to primary tumors, TERT promoter mutations increased while BRAF mutations decreased in metastatic lesions, indicating that BRAF may play a limited role during distant metastasis, whereas TERT mutations could be more critical. Similarly, Cheng X et al. (41) found a significantly higher frequency of BRAF mutations in PTC patients without distant metastases compared to those with metastases (84.8% vs. 51.5%, P < 0.0001), aligning with our observations.

The molecular mechanisms underlying this inverse association between BRAF mutation and lung metastasis remain unclear and warrant further investigation. One hypothesis is based on the concept of mutual exclusivity among oncogenic drivers, which has been reported in other studies (42, 43). We speculate that in this cohort of young patients with papillary thyroid carcinoma, the genetic drivers that result in distant metastasis may be other genes that are mutually exclusive with the BRAF V600E mutation. According to this hypothesis, when the BRAF V600E mutation is present, key genetic alterations required to complete the lung metastasis cascade may be absent or suppressed, potentially explaining the observed negative correlation. This mechanistic hypothesis remains speculative and requires direct validation through future studies.

However, the coexistence of BRAF and TERT mutations creates a powerful synergistic effect that profoundly enhances metastatic potential (44, 45). In our training cohort, the lung metastasis rate was 5.56% in patients harboring only the BRAF mutation, but this rate rose dramatically to 28.57% when both BRAF and TERT mutations were present. A similar pattern was observed in the validation cohort (Supplementary Table 2). This synergy can be explained by their mechanistic interaction. the BRAF-driven MAPK pathway can upregulate TERT promoter activity (46), while TERT cooperates with mutant BRAF to drive tumor dedifferentiation and progression through mechanisms such as the modulation of ribosomal biogenesis (47). Therefore, BRAF and TERT mutations may synergistically enhance metastasis.

Furthermore, our study characterized genetic mutations in a binary manner (presence/absence), without incorporating quantitative metrics such as variant allele frequency (VAF) or mutational burden. Emerging evidence suggests that VAF may reflect the clonal dominance of a mutation within the tumor, which could be associated with more aggressive behavior and poorer outcomes in thyroid cancer (48). Future studies incorporating VAF and other quantitative genetic features could provide a more nuanced understanding of tumor biology and further refine risk stratification models for lung metastasis in young patients with PTC.

We also conducted a comprehensive analysis of other genetic markers, including RAS family mutations (NRAS, HRAS, KRAS), TP53, and PIK3CA. As shown in Table 2, univariate analyses conducted within the training cohort confirmed that none of these mutations demonstrated a significant association with prognosis (RAS: OR = 1.488, P = 0.573; TP53: OR = 0.632, P = 0.676; PIK3CA: OR = 0.765, P = 0.809) and were therefore excluded from the final model. Future studies with larger sample sizes are needed to clearly evaluate the contribution of rare genetic events.

Not all regions or patients have access to genetic testing. Therefore, the clinical utility of this nomogram should be considered. We compared the predictive performance of the model using only two clinical variables (T stage and unilateral thyroid involvement) with that of the model utilizing all variables. As shown in Supplementary Figure 2, the model incorporating only clinical variables demonstrated reasonable predictive value, though its performance was markedly inferior to that of the full model integrating genetic data. Thus, to achieve optimal predictive accuracy, inclusion of BRAF and TERT mutation status is strongly recommended.

Although the 425-gene panel provides more genetic information compared to specific targeted gene testing, the 425-gene panel may prove difficult to implement widely in clinical practice due to its high cost. As our multivariate analysis revealed, only a subset of genes—specifically TERT and BRAF—along with clinical factors (T stage and unilateral thyroid involvement) were independently predictive of lung metastasis. This suggests that by detecting few specific genes(such as BRAF and TERT promoter mutation) and combining them with key clinical variables, we can predict lung metastasis, thereby reducing testing costs and enhancing feasibility.

There are some limitations to our study. First, as a retrospective analysis, it may be subject to selection bias during patient enrollment. Second, the model has only been validated internally. Although it demonstrates satisfactory performance within our cohort, external validation using independent, multi-institutional data is necessary to confirm its generalizability before clinical application, particularly regarding the cut-off value. We are actively collaborating with multiple leading cancer centers in China to establish an independent, multicenter cohort of young PTC patients for external validation, which will optimize the model and enhance its potential for translational application. Third, although other candidate genes such as RAS have been implicated in thyroid cancer metastasis (49), too few patients (only 3 with lung metastasis) presented with these mutations in our cohort to allow meaningful analysis. Larger samples are needed to evaluate other target genes and improve model accuracy. Fourth, we lacked data on patients’ thyroglobulin(Tg) levels. Tg serves as a core indicator for postoperative monitoring of differentiated thyroid cancer, demonstrating high sensitivity in detecting distant metastasis (50). Consequently, the absence of this data may diminish the practical value of this model. However, based on baseline genetic and pathological assessments, our model can identify young patients at high risk for lung metastasis during the early diagnostic stages—even before Tg levels change. Simultaneously, for specific patients with difficult-to-interpret Tg results (potentially affected by Tg antibodies), this model still provides supplementary predictive value to some extent. Therefore, we will integrate Tg data in the future to construct a more comprehensive predictive model and enhance its forecasting capabilities.

Conclusion

We developed and validated a nomogram that integrates T stage, unilateral thyroid involvement, BRAF V600E, and TERT promoter mutations to individually predict the probability of lung metastasis in PTC patients under 55 years of age. This tool may assist clinicians in early risk stratification and personalized management.

Data availability statement

The original data are available from the corresponding author.

Ethics statement

The study was approved by the Ethics Committee of The Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University/Hunan Cancer Hospital (approval number: KYJJ-2020-195). Written informed consent from the patients was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author contributions

HY: Funding acquisition, Writing – original draft, Writing – review & editing. QQ: Data curation, Methodology, Writing – review & editing. YL: Data curation, Formal analysis, Writing – review & editing. WP: Conceptualization, Formal analysis, Investigation, Writing – review & editing. JP: Software, Validation, Writing – review & editing. ZZ: Project administration, Visualization, Writing – review & editing. HZ: Formal analysis, Methodology, Writing – review & editing. FS: Funding acquisition, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. The study was supported by the Science and Technology Innovation Program of Hunan Province (2023SK4043),the Research Program Project of Hunan Provincial Health and Wellness Commission (202103060919) and the National Natural Science Foundation of China (82404040). High-Level talent Support Program of Hunan Cancer Hospital (20250805-1001). High-Level talent Support Program of Hunan Cancer Hospital(20250731-1040).

Acknowledgments

The authors sincerely thank all the patients for their consent to publish this information.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1689674/full#supplementary-material

References

1. Cabanillas ME, McFadden DG, and Durante C. Thyroid cancer. Lancet. (2016) 388:2783–95. doi: 10.1016/S0140-6736(16)30172-6

PubMed Abstract | Crossref Full Text | Google Scholar

2. O'Neill CJ, Oucharek J, Learoyd D, and Sidhu SB. Standard and emerging therapies for metastatic differentiated thyroid cancer. Oncologist. (2010) 15:146–56. doi: 10.1634/theoncologist.2009-0190

PubMed Abstract | Crossref Full Text | Google Scholar

3. Khan U, Al Afif A, Aldaihani A, MacKay C, Rigby MH, Rajaraman M, et al. Patient and tumor factors contributing to distant metastasis in well-differentiated thyroid cancer: a retrospective cohort study. J Otolaryngol Head Neck Surg. (2020) 49:78. doi: 10.1186/s40463-020-00469-8

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wang W, Shen C, and Yang Z. Nomogram individually predicts the risk for distant metastasis and prognosis value in female differentiated thyroid cancer patients: A SEER-based study. Front Oncol. (2022) 12:800639. doi: 10.3389/fonc.2022.800639

PubMed Abstract | Crossref Full Text | Google Scholar

5. Lan X, Bao H, Ge X, Cao J, Fan X, Zhang Q, et al. Genomic landscape of metastatic papillary thyroid carcinoma and novel biomarkers for predicting distant metastasis. Cancer Sci. (2020) 111:2163–73. doi: 10.1111/cas.14389

PubMed Abstract | Crossref Full Text | Google Scholar

6. Toraih EA, Hussein MH, Zerfaoui M, Attia AS, Marzouk Ellythy A, Mostafa A, et al. Site-specific metastasis and survival in papillary thyroid cancer: the importance of brain and multi-organ disease. Cancers (Basel). (2021) 13:1625. doi: 10.3390/cancers13071625

PubMed Abstract | Crossref Full Text | Google Scholar

7. Sohn SY, Kim HI, Kim YN, Kim TH, Kim SW, and Chung JH. Prognostic indicators of outcomes in patients with lung metastases from differentiated thyroid carcinoma during long-term follow-up. Clin Endocrinol (Oxf). (2018) 88:318–26. doi: 10.1111/cen.13489

PubMed Abstract | Crossref Full Text | Google Scholar

8. Thedsawad A, Wanachiwanawin W, Taka O, and Hantaweepant C. Cut-off values for diagnosis of G6PD deficiency by flow cytometry in Thai population. Ann Hematol. (2022) 101:2149–57. doi: 10.1007/s00277-022-04923-7

PubMed Abstract | Crossref Full Text | Google Scholar

9. Li Y, Gao X, Guo T, and Liu J. Development and validation of a nomogram for risk of pulmonary metastasis in non-papillary thyroid carcinoma: A SEER-based study. Med (Baltimore). (2023) 102:e34581. doi: 10.1097/MD.0000000000034581

PubMed Abstract | Crossref Full Text | Google Scholar

10. Guo P, Wang X, Xia L, Shawureding N, and Hu Z. Analysis of factors associated with the prognosis of papillary thyroid cancer and the construction of a survival model. Cancer Med. (2023) 12:7868–76. doi: 10.1002/cam4.5555

PubMed Abstract | Crossref Full Text | Google Scholar

11. Mosele MF, Westphalen CB, Stenzinger A, Barlesi F, Bayle A, Bièche I, et al. Recommendations for the use of next-generation sequencing (NGS) for patients with advanced cancer in 2024: a report from the ESMO Precision Medicine Working Group. Ann Oncol. (2024) 35:588–606. doi: 10.1016/j.annonc.2024.04.005

PubMed Abstract | Crossref Full Text | Google Scholar

12. Lamartina L, Grani G, Arvat E, Nervo A, Zatelli MC, Rossi R, et al. 8th edition of the AJCC/TNM staging system of thyroid cancer: what to expect (ITCO2). Endocr Relat Can. (2018) 25:L7–L11. doi: 10.1530/ERC-17-0453

PubMed Abstract | Crossref Full Text | Google Scholar

13. Mazurat A, Torroni A, Hendrickson-Rebizant J, Benning H, Nason RW, and Pathak KA. The age factor in survival of a population cohort of well-differentiated thyroid cancer. Endocr Con. (2013) 2:154–60. doi: 10.1530/EC-13-0056

PubMed Abstract | Crossref Full Text | Google Scholar

14. Alzahrani AS, Albalawi L, Mazi S, Mukhtar N, Aljamei H, Moria Y, et al. How does the AJCC/TNM staging system eighth edition perform in thyroid cancer at A major middle eastern medical center? Endocr Pract. (2021) 27:607–13. doi: 10.1016/j.eprac.2020.11.018

PubMed Abstract | Crossref Full Text | Google Scholar

15. Sugino K, Nagahama M, Kitagawa W, Ohkuwa K, Uruno T, Matsuzu K, et al. Distant metastasis in pediatric and adolescent differentiated thyroid cancer: clinical outcomes and risk factor analyses. J Clin Endocrinol Metab. (2020) 105:dgaa545. doi: 10.1210/clinem/dgaa545

PubMed Abstract | Crossref Full Text | Google Scholar

16. Parvathareddy SK, Siraj AK, Annaiyappanaidu P, Siraj N, Al-Rasheed M, Al-Haqawi W, et al. Predictive risk factors for distant metastasis in pediatric differentiated thyroid cancer from Saudi Arabia. Front Endocrinol (Lausanne). (2023) 14:1228049. doi: 10.3389/fendo.2023.1228049

PubMed Abstract | Crossref Full Text | Google Scholar

17. Kim H, Kim HI, Kim SW, Jung J, Jeon MJ, Kim WG, et al. Prognosis of differentiated thyroid carcinoma with initial distant metastasis: A multicenter study in korea. Endocrinol Metab (Seoul). (2018) 33:287–95. doi: 10.3803/EnM.2018.33.2.287

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hirsch D, Levy S, Tsvetov G, Gorshtein A, Slutzky-Shraga I, Akirov A, et al. Long-term outcomes and prognostic factors in patients with differentiated thyroid cancer and distant metastases. Endocr Pract. (2017) 23:1193–200. doi: 10.4158/EP171924

Crossref Full Text | Google Scholar

19. Kim M, Kim WG, Park S, Kwon H, Jeon MJ, Lee JJ, et al. Initial size of metastatic lesions is best prognostic factor in patients with metastatic differentiated thyroid carcinoma confined to the lung. Thyroid. (2017) 27:49–58. doi: 10.1089/thy.2016.0347

PubMed Abstract | Crossref Full Text | Google Scholar

20. Chesover AD, Vali R, Hemmati SH, and Wasserman JD. Lung metastasis in children with differentiated thyroid cancer: factors associated with diagnosis and outcomes of therapy. Thyroid. (2021) 31:50–60. doi: 10.1089/thy.2020.0002

PubMed Abstract | Crossref Full Text | Google Scholar

21. Kuang HF and Lu WL. Predictive factors for lung metastasis in pediatric differentiated thyroid cancer: a clinical prediction study. J Pediatr Endocrinol Metab. (2024) 37:250–9. doi: 10.1515/jpem-2023-0425

PubMed Abstract | Crossref Full Text | Google Scholar

22. Nguyen NC, Anigati EM, Desai NB, and Öz OK. Radioactive iodine therapy in differentiated thyroid cancer: an update on dose recommendations and risk of secondary primary Malignancies. Semin Nucl Med. (2024) 54:488–96. doi: 10.1053/j.semnuclmed.2024.05.002

PubMed Abstract | Crossref Full Text | Google Scholar

23. Omi Y, Yanagida J, Yoshida Y, Horiuchi K, and Okamoto T. Risk factors for mediastinal lymph node metastasis and lung metastasis in papillary thyroid carcinoma patients: who benefits from preoperative computed tomography? Endocr J. (2024) 71:357–62. doi: 10.1507/endocrj.EJ23-0413

PubMed Abstract | Crossref Full Text | Google Scholar

24. Liu X, Bishop J, Shan Y, Pai S, Liu D, Murugan AK, et al. Highly prevalent TERT promoter mutations in aggressive thyroid cancers. Endocr Relat Can. (2013) 20:603–10. doi: 10.1530/ERC-13-0210

PubMed Abstract | Crossref Full Text | Google Scholar

25. Hafezi F and Perez Bercoff D. The solo play of TERT promoter mutations. Cells. (2020) 9:749. doi: 10.3390/cells9030749

PubMed Abstract | Crossref Full Text | Google Scholar

26. Yang L, Wang M, Li N, Yan LD, Zhou W, Yu ZQ, et al. TERT mutations in non-small cell lung cancer: clinicopathologic features and prognostic implications. Clin Med Insights Oncol. (2023) 17:11795549221140781. doi: 10.1177/11795549221140781

PubMed Abstract | Crossref Full Text | Google Scholar

27. Melo M, Gaspar da Rocha A, Batista R, Vinagre J, Martins MJ, Costa G, et al. TERT, BRAF, and NRAS in primary thyroid cancer and metastatic disease. J Clin Endocrinol Metab. (2017) 102(6):1898–907. doi: 10.1210/jc.2016-2785

PubMed Abstract | Crossref Full Text | Google Scholar

28. Liu C, Liu Z, Chen T, Zeng W, Guo Y, and Huang T. TERT promoter mutation and its association with clinicopathological features and prognosis of papillary thyroid cancer: A meta-analysis. Sci Rep. (2016) 6:36990. doi: 10.1038/srep36990

PubMed Abstract | Crossref Full Text | Google Scholar

29. Parvathareddy SK, Siraj AK, Iqbal K, Qadri Z, Ahmed SO, Al-Rasheed M, et al. TERT promoter mutations are an independent predictor of distant metastasis in middle eastern papillary thyroid microcarcinoma. Front Endocrinol (Lausanne). (2022) 13:808298. doi: 10.3389/fendo.2022.808298

PubMed Abstract | Crossref Full Text | Google Scholar

30. Rossi M, Buratto M, Tagliati F, Rossi R, Lupo S, Trasforini G, et al. Relevance of BRAF(V600E) mutation testing versus RAS point mutations and RET/PTC rearrangements evaluation in the diagnosis of thyroid cancer. Thyroid. (2015) 25:221–8. doi: 10.1089/thy.2014.0338

PubMed Abstract | Crossref Full Text | Google Scholar

31. Pappa T, Ahmadi S, Marqusee E, Johnson HL, Nehs MA, Cho NL, et al. Oncogenic mutations in PI3K/AKT/mTOR pathway effectors associate with worse prognosis in BRAFV600E -driven papillary thyroid cancer patients. Clin Cancer Res. (2021) 27:4256–64. doi: 10.1158/1078-0432.CCR-21-0874

PubMed Abstract | Crossref Full Text | Google Scholar

32. Xing M. BRAF mutation and thyroid cancer recurrence. J Clin Oncol. (2015) 33:2482–3. doi: 10.1200/JCO.2015.61.4016

PubMed Abstract | Crossref Full Text | Google Scholar

33. Subbiah V, Baik C, and Kirkwood JM. Clinical development of BRAF plus MEK inhibitor combinations. Trends Can. (2020) 6:797–810. doi: 10.1016/j.trecan.2020.05.009

PubMed Abstract | Crossref Full Text | Google Scholar

34. Yan C, Huang M, Li X, Wang T, and Ling R. Relationship between BRAF V600E and clinical features in papillary thyroid carcinoma. Endocr Con. (2019) 8:988–96. doi: 10.1530/EC-19-0246

PubMed Abstract | Crossref Full Text | Google Scholar

35. Geng J, Wang H, Liu Y, Tai J, Jin Y, Zhang J, et al. Correlation between BRAF V600E mutation and clinicopathological features in pediatric papillary thyroid carcinoma. Sci China Life Sci. (2017) 60(7):729–38. doi: 10.1007/s11427-017-9083-8

PubMed Abstract | Crossref Full Text | Google Scholar

36. Nucera C, Lawler J, and Parangi S. BRAF(V600E) and microenvironment in thyroid cancer: a functional link to drive cancer progression. Cancer Res. (2011) 71:2417–22. doi: 10.1158/0008-5472.CAN-10-3844

PubMed Abstract | Crossref Full Text | Google Scholar

37. Tufano RP, Teixeira GV, Bishop J, Carson KA, and Xing M. BRAF mutation in papillary thyroid cancer and its value in tailoring initial treatment: a systematic review and meta-analysis. Med (Baltimore). (2012) 91:274–86. doi: 10.1097/MD.0b013e31826a9c71

PubMed Abstract | Crossref Full Text | Google Scholar

38. Sancisi V, Nicoli D, Ragazzi M, Piana S, and Ciarrocchi A. BRAFV600E mutation does not mean distant metastasis in thyroid papillary carcinomas. J Clin Endocrinol Metab. (2012) 97:E1745–9. doi: 10.1210/jc.2012-1526

PubMed Abstract | Crossref Full Text | Google Scholar

39. Chakravarty D, Santos E, Ryder M, Knauf JA, Liao XH, West BL, et al. Small-molecule MAPK inhibitors restore radioiodine incorporation in mouse thyroid cancers with conditional BRAF activation. J Clin Invest. (2011) 121:4700–11. doi: 10.1172/JCI46382

PubMed Abstract | Crossref Full Text | Google Scholar

40. Ito Y, Yoshida H, Maruo R, Morita S, Takano T, Hirokawa M, et al. BRAF mutation in papillary thyroid carcinoma in a Japanese population: its lack of correlation with high-risk clinicopathological features and disease-free survival of patients. Endocr J. (2009) 56:89–97. doi: 10.1507/endocrj.k08e-208

PubMed Abstract | Crossref Full Text | Google Scholar

41. Cheng X, Zhou Y, Xu S, Yu H, Wu J, Bao J, et al. Risk-stratified distant metastatic thyroid cancer with clinicopathological factors and BRAF/TERT promoter mutations. Exp Clin Endocrinol Diabet. (2023) 131:577–82. doi: 10.1055/a-2177-1051

PubMed Abstract | Crossref Full Text | Google Scholar

42. Di Cristofaro J, Marcy M, Vasko V, Sebag F, Fakhry N, Wynford-Thomas D, et al. Molecular genetic study comparing follicular variant versus classic papillary thyroid carcinomas: association of N-ras mutation in codon 61 with follicular variant. Hum Pathol. (2006) 37:824–30. doi: 10.1016/j.humpath.2006.01.030

PubMed Abstract | Crossref Full Text | Google Scholar

43. Bazhenova L, Lokker A, Snider J, Castellanos E, Fisher V, Fellous M, et al. TRK fusion cancer: patient characteristics and survival analysis in the real-world setting. Tar Oncol. (2021) 16:389–99. doi: 10.1007/s11523-021-00815-4

PubMed Abstract | Crossref Full Text | Google Scholar

44. Xing M, Liu R, Liu X, Murugan AK, Zhu G, Zeiger MA, et al. BRAF V600E and TERT promoter mutations cooperatively identify the most aggressive papillary thyroid cancer with highest recurrence. J Clin Oncol. (2014) 32:2718–26. doi: 10.1200/JCO.2014.55.5094

PubMed Abstract | Crossref Full Text | Google Scholar

45. Ngeow J and Eng C. TERT and BRAF in thyroid cancer: teaming up for trouble. J Clin Oncol. (2014) 32:2683–4. doi: 10.1200/JCO.2014.56.5614

PubMed Abstract | Crossref Full Text | Google Scholar

46. Liu X, Qu S, Liu R, Sheng C, Shi X, Zhu G, et al. TERT promoter mutations and their association with BRAF V600E mutation and aggressive clinicopathological characteristics of thyroid cancer. J Clin Endocrinol Metab. (2014) 99:E1130–6. doi: 10.1210/jc.2013-4048

PubMed Abstract | Crossref Full Text | Google Scholar

47. Yu P, Qu N, Zhu R, Hu J, Han P, Wu J, et al. TERT accelerates BRAF mutant-induced thyroid cancer dedifferentiation and progression by regulating ribosome biogenesis. Sci Adv. (2023) 9:eadg7125. doi: 10.1126/sciadv.adg7125

PubMed Abstract | Crossref Full Text | Google Scholar

48. Fu G, Chazen RS, Monteiro E, Vescan A, Freeman JL, Witterick IJ, et al. Facilitation of definitive cancer diagnosis with quantitative molecular assays of BRAF V600E and TERT promoter variants in patients with thyroid nodules. JAMA Netw Open. (2023) 6:e2323500. doi: 10.1001/jamanetworkopen.2023.23500

PubMed Abstract | Crossref Full Text | Google Scholar

49. Riccio I, Laforteza A, Landau MB, Hussein MH, Linhuber J, Staav J, et al. Decoding RAS mutations in thyroid cancer: A meta-analysis unveils specific links to distant metastasis and increased mortality. Am J Otolaryngol. (2025) 46:104570. doi: 10.1016/j.amjoto.2024.104570

PubMed Abstract | Crossref Full Text | Google Scholar

50. Guo P, Wang X, Xia L, Shawureding N, and Hu Z. Analysis of factors associated with the prognosis of papillary thyroid cancer and the construction of a survival model. Cancer Med. (2023) 12:7868–76. doi: 10.1002/cam4.5555

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: thyroid cancer, lung metastasis, nomogram, next-generation sequencing, prognosis

Citation: Yang H, Qian Q, Li Y, Pan W, Peng J, Zou Z, Zhu H and Shi F (2025) A Nomogram for Predicting Lung Metastasis in Papillary Thyroid Cancer Patients Aged Less Than 55 Years. Front. Endocrinol. 16:1689674. doi: 10.3389/fendo.2025.1689674

Received: 20 August 2025; Accepted: 24 October 2025;
Published: 11 November 2025.

Edited by:

Dario de Biase, University of Bologna, Italy

Reviewed by:

Mohd Amir, Aligarh Muslim University, India
Sijia Zhang, Heidelberg University, Germany

Copyright © 2025 Yang, Qian, Li, Pan, Peng, Zou, Zhu and Shi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Feng Shi, c2hpZmVuZ0BobmNhLm9yZy5jbg==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.