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

SYSTEMATIC REVIEW article

Front. Oncol., 29 January 2026

Sec. Breast Cancer

Volume 16 - 2026 | https://doi.org/10.3389/fonc.2026.1666576

Elevated thymidine kinase 1 expression at baseline predicts poor prognosis in breast cancer patients

Peng LiPeng Li1Yongting Cheng*Yongting Cheng2*Junfeng ZhaoJunfeng Zhao1Yuan FangYuan Fang1Tingting ZhaoTingting Zhao1
  • 1Department of Medical Oncology, The First Affiliated Hospital of Hebei North University, Zhang jia kou, Hebei, China
  • 2Department of Medical Laboratory Science, Hebei North University, Zhang jia kou, Hebei, China

Introduction: Thymidine kinase 1 (TK1), a key enzyme in DNA biosynthesis, has been shown to correlate with breast cancer prognosis and treatment response in dynamic monitoring settings. However, the clinical relevance of baseline TK1 levels remains controversial due to inconsistent evidence across studies. To address this issue, we conducted the first systematic meta-analysis of available studies to investigate the potential association between baseline TK1 levels and prognostic outcomes in breast cancer patients.

Methods: A comprehensive computerized literature search was conducted across major Chinese and English databases to identify studies investigating the association between TK1 expression and breast cancer prognosis. Baseline TK1 expression levels and corresponding patient survival data were systematically extracted for meta-analysis.

Results: The meta-analysis evaluating the association between baseline TK1 expression levels and progression-free survival (PFS) in breast cancer included 2,887 patients from 11 studies. Significant heterogeneity was observed across the included studies (I2 = 87.9%, p = 0.099), which persisted even after subgroup analyses. Therefore, a random-effects model was employed, yielding a pooled hazard ratio (HR) of 1.63 (95% confidence interval [CI]: 1.28–2.10, p = 0.000, Z = 3.88). The meta-analysis evaluating the association between baseline TK1 expression levels and OS in breast cancer included 2,233 patients from six studies. Significant heterogeneity was initially observed (I2 = 72.3%, p = 0.003), which was resolved through subgroup stratification by treatment status (treatment-naive versus recurrent disease). In the treatment-naive subgroup, the HR was 1.30 (95% CI: 1.11–1.52, p = 0.001, Z = 3.26). For the recurrent disease subgroup, the HR was 2.10 (95% CI: 1.74–2.54, p = 0.000, Z = 7.64).

Conclusion: Breast cancer patients presenting with high baseline TK1 expression are associated with significantly worse prognostic outcomes. Collectively, these findings support the clinical potential of TK1 assessment for prognostic risk stratification and treatment guidance, which merits further verification in large-scale, multicenter clinical trials.

1 Introduction

Invasive breast cancer (IBC) is the most commonly diagnosed malignancy and the second leading cause of cancer-related mortality in women worldwide. Although standardized multimodal therapies have significantly improved disease-free survival (DFS) and overall survival (OS) in breast cancer patients, substantial interindividual variability in treatment response persists, driven by tumor biological heterogeneity and complex clinical factors (1). This heterogeneity underscores the critical need for robust prognostic biomarkers to enable precise risk stratification and personalized therapeutic decision-making.

Current clinical practice relies on established prognostic indicators, including hormone receptor status, human epidermal growth factor receptor 2 (HER2) expression, and Ki-67 proliferation index. The integration of comprehensive genomic profiling has further enhanced risk assessment and therapeutic selection, particularly for early-stage disease. However, despite these advancements, many proposed biomarkers exhibit inconsistent clinical utility and have not been widely adopted into routine clinical practice (2, 3).

A major clinical challenge remains the significant proportion of patients who experience disease recurrence or metastasis despite receiving optimal standard therapy, highlighting the urgent need for more accurate prognostic tools. The discovery of novel, reliable biomarkers could substantially improve risk prediction and facilitate the development of tailored treatment strategies.

In humans, thymidine kinase exists as two distinct isoenzymes: cytosolic thymidine kinase 1 (TK1) and mitochondrial thymidine kinase 2 (TK2) (4). As a cell cycle-regulated proliferation marker, TK1 is markedly upregulated during the S-phase to facilitate DNA synthesis. It plays a pivotal role in the pyrimidine salvage pathway—one of the two major routes for DNA precursor biosynthesis (alongside de novo synthesis)—thereby contributing to DNA replication and repair (5). Under normal physiological conditions, serum TK1 activity in healthy individuals is negligible or undetectable. By contrast, it is substantially elevated in malignant tumors, with strong correlations to tumor proliferation, progression, and metastatic potential (6). These findings have positioned TK1 as a promising biomarker for prognostic assessment and therapeutic decision-making in oncology.

While the utility of serial TK1 monitoring for treatment response and prognosis in breast cancer has been established, significant discrepancies persist regarding the diagnostic, prognostic, and predictive value of baseline TK1 levels (79). The correlation between baseline TK1 levels and distinct prognostic outcomes in breast cancer patients, as well as its potential utility in guiding diagnosis and treatment, remains controversial due to inconsistent conclusions across studies. To address this evidence gap, we conducted the first meta-analysis to systematically assess the relationship between baseline TK1 levels and prognosis in breast cancer.

2 Methods

2.1 Literature search

A comprehensive literature search was performed in both Chinese and international databases up to 31 May 2025. Chinese databases included China National Knowledge Infrastructure (CNKI), Weipu Information (VIP), Wanfang, and Chinese Biomedical Literature Database (CBMDisc), while English databases comprised PubMed (US National Library of Medicine), EMBASE, and the Cochrane Library. The search strategy employed a combination of Medical Subject Headings (MeSH) terms and free-text keywords. The following search terms were used: “thymidine kinase”, “thymidine kinase 1”, “TK”, or “TK1”, in combination with “breast cancer”. Both Chinese and English search terms were adapted to the respective databases to ensure comprehensive coverage.

2.2 Inclusion criteria

1. Original research articles published in peer-reviewed journals, with full text available in English or Chinese, focusing on primary breast cancer patients.

2. Studies investigating the association between TK1 expression levels and breast cancer prognosis, including OS, DFS, progression-free survival (PFS), or recurrence-free survival (RFS).

3. Clear definition of TK1 expression status, with patients stratified into high- and low-TK1 expression groups based on standardized thresholds (e.g., immunohistochemistry cutoff values, mRNA levels, or enzyme activity).

4. TK1 measurements obtained at initial diagnosis or after a treatment-free interval of ≥6 months, ensuring that TK1 levels were assessed prior to any therapeutic intervention (e.g., surgery, chemotherapy, or radiotherapy).

5. Reported survival outcomes (e.g., OS, DFS/PFS, RFS) with adjusted hazard ratios (HRs) and 95% confidence intervals (CIs) derived from multivariate Cox regression analyses (preferred) or univariate analyses if multivariate data were unavailable.

6. Sufficient statistical data to directly extract or calculate HRs and 95% CIs (e.g., Kaplan-Meier survival curves with log-rank P-values, event numbers, or regression coefficients).

2.3 Exclusion criteria

1. Non-English/non-Chinese articles; cross-sectional studies, experimental studies, review articles, conference abstracts, case reports, and letters.

2. Studies focusing on nonprimary breast cancer.

3. Studies lacking survival outcomes (e.g., OS, DFS) or with insufficient data to calculate/extract HRs with 95% CIs.

4. Studies unable to compare baseline TK1 levels (e.g., missing two-group HR data for high/low expression) or with insufficient data for meta-analysis pooling.

5. Studies with overlapping patient cohorts or duplicate data sources.

2.4 Information extraction

Two independent researchers systematically extracted data from eligible studies, with any discrepancies resolved through discussion or by consulting a third reviewer. Key extracted information included study characteristics (first author, publication year, country, design), patient demographics (sample size, age, stage), TK1 detection methods (specimen type, assay technique, cutoff values), and survival outcomes (HRs with 95% CIs for OS/DFS/PFS). For studies with multiple reports, the most comprehensive dataset was selected. Study quality was assessed using the Newcastle-Ottawa Scale (NOS), evaluating cohort selection, comparability, and outcome assessment. All extracted data were cross-verified to ensure accuracy prior to analysis.

2.5 Statistical analysis

Statistical analyses were performed using Stata 12.0 (StataCorp College Station, TX, United States), with HRs and 95% CIs for OS and PFS serving as the primary effect measures. Heterogeneity was assessed using I2 statistics and Cochran’s Q test (p < 0.10), with fixed-effects models applied when I2 < 50% and random-effects models applied when I2 ≥ 50%. For studies lacking reported HRs, estimates were derived from available survival data (10). Sensitivity analyses included sequential study exclusion and subgroup analyses. Publication bias was evaluated through funnel plots, Begg’s test, Egger’s test, and the trim-and-fill method. All tests were two-tailed, with statistical significance set at p < 0.05.

3 Results

3.1 Literature search

The systematic literature search initially identified 790 potential studies from Chinese (CNKI, VIP, Wanfang, CBMDisc) and English (PubMed) databases. After removing 732 irrelevant records through title/abstract screening, 58 full-text articles were assessed for eligibility. Based on our predefined criteria, 47 studies were excluded, yielding 11 qualified studies for final meta-analysis (Figure 1).

Figure 1
Flowchart illustrating the literature selection process for a meta-analysis. Initially, 86 articles were retrieved from Chinese databases and 704 from English databases. After exclusions based on citations and abstracts (732), 58 articles were evaluated in full text. Of these, 38 were excluded, leaving 20 for re-evaluation. Following re-evaluation, eight were excluded, with reasons including overlapping data sources and insufficient data. Ultimately, 11 articles were included in the meta-analysis.

Figure 1. Literature identification and screening process.

3.2 Characteristics of included studies and literature quality evaluation

A total of 11 studies were included, with sample sizes ranging from 31 to 1,310. Based on the general disease characteristics of enrolled patients, five studies included primary treatment, while the remaining six included patients with recurrent disease. All 11 studies evaluated the correlation between TK1 levels and breast cancer prognosis. Although all reported PFS or DFS data comparisons, six also provided OS outcome comparisons. Geographically, four studies were conducted in Italy, four in Sweden, and the remaining three in France, the USA, and Israel. Regarding detection methods, nine studies used enzyme-linked immunosorbent assay (ELISA), one used immunohistochemistry, and one employed the bicinchoninic acid (BCA) method. Notably, six studies conducted multivariate analyses of factors associated with breast cancer prognosis. All eligible studies were identified from English-language databases and demonstrated high quality, as reflected in NOS scores of 8 to 9 (Table 1).

Table 1
www.frontiersin.org

Table 1. Characteristics of included studies.

3.3 Meta-analysis of baseline TK1 expression and prognosis in breast cancer patients

A meta-analysis was first performed to explore the correlation between baseline TK1 expression levels and PFS in breast cancer patients. For the PFS-focused studies (Table 2), heterogeneity testing revealed significant heterogeneity among the included studies (I2 = 87.9%, p = 0.099). A meta-regression analysis was conducted to assess the impact of treatment-naive/recurrent treatment status (as part of enrollment characteristics) on heterogeneity. The results showed that the treatment-naive/recurrent variable explained 34.03% of the heterogeneity, which was not statistically significant (Adj R-squared = 34.03%, p = 0.175). Considering that, in addition to treatment-naive/recurrent treatment status, baseline characteristics such as patient pathology and treatment modality across studies might also contribute to substantial heterogeneity—and that subgroup analyses could not fully eliminate this heterogeneity—a random-effects model was selected for the meta-analysis. The results indicated that baseline high TK1 expression predicted poorer PFS (HR: 1.63, 95% CI: 1.28–2.10, p < 0.001, Z = 3.88) (Figure 2).

Table 2
www.frontiersin.org

Table 2. Association between baseline TK1 expression and PFS in breast cancer.

Figure 2
Forest plot showing hazard ratios (HR) with 95% confidence intervals for eleven studies from 2001 to 2024, listed alongside percentage weights. The HRs are depicted by squares with lines representing confidence intervals. An overall analysis is shown at the bottom with a diamond, indicating a pooled HR of 1.63. A vertical line at HR equals 1 serves as a reference point. Weights are from random effects analysis.

Figure 2. Meta-analysis of PFS in 11 studies (forest plot).

A meta-analysis was performed to investigate the association between baseline TK1 expression and OS in breast cancer patients. For OS-focused studies (Table 3), heterogeneity testing revealed significant interstudy heterogeneity (I2 = 72.3%, p = 0.003). A meta-regression analysis was conducted to assess the impact of treatment-naive/recurrence status (as part of enrollment characteristics) on heterogeneity. The results showed that the treatment-naive/recurrence variable explained 100% of the heterogeneity, which was statistically significant (Adj R-squared = 100%, p = 0.019). Based on baseline treatment-naive/recurrence status, relevant studies were divided into two subgroups, and fixed-effects models were applied for separate meta-analyses. No significant heterogeneity was observed in either the treatment-naive subgroup (I2 = 21.4%, p = 0.28) or the recurrence subgroup (I2 = 0.0%, p = 0.645). The meta-analysis results indicated that baseline high TK1 expression predicted poorer OS in both subgroups (treatment-naive: HR = 1.30, 95% CI: 1.11–1.52, p = 0.001, Z = 3.26; recurrence: HR = 2.10, 95% CI: 1.74–2.54, p < 0.001, Z = 7.64) (Figure 3).

Table 3
www.frontiersin.org

Table 3. Association between baseline TK1 expression and OS in breast cancer.

Figure 3
Forest plot depicting hazard ratios (HRs) with 95% confidence intervals (CIs) from multiple studies. Group 1 includes studies by Broet (2001), Fanelli (2021), and Matikas (2021) with a subtotal HR of 1.30. Group 2 includes studies by Bjøhle (2013), Larsson (2020), and Bergqvist (2023) with a subtotal HR of 2.10. Overall HR is 1.58. Weight percentages are shown for each study, and heterogeneity statistics are provided.

Figure 3. Meta-analysis of OS in six studies (forest plot).

3.4 Sensitivity analysis

In this study, sensitivity analyses were performed using leave-one-out exclusion to evaluate the robustness of the meta-analysis results. For PFS-related studies, sequential exclusion of each included study showed that the direction of pooled effect sizes in the remaining studies remained consistent, and the 95% CIs of the pooled effect sizes did not include the null value of 1.0. These results indicated that the absence of any individual study did not significantly alter the overall conclusions, suggesting that the meta-analysis findings were not substantially influenced by single studies (Figure 4).

Figure 4
Forest plot depicting meta-analysis estimates with specific studies omitted. Each row represents a study with circles indicating estimates and lines representing confidence intervals. Studies listed include Broet, Nisman, BJohle, Bonechi, McCartney (2019 and 2020), Larsson, Fanelli, Matikas, Bergqvist, and Zhu Y. The x-axis spans from 1.22 to 2.22 with a center line at 1.63, showing variations in confidence intervals across studies.

Figure 4. Sensitivity analysis of PFS data.

In the sensitivity analysis of the OS treatment-naive subgroup, exclusion of the 2001 Broet study—which enrolled a patient cohort with a notably distinct treatment background compared to those in subsequent investigations—resulted in the 95% CIs of the HR encompassing 1.0. Furthermore, whereas TK1 was measured in tumor tissue in this study, the majority of subsequent investigations utilized serum samples. This biological discrepancy in sample sources may introduce systematic bias when comparing the prognostic utility of TK1 across these studies. Combined with the sensitivity analysis results and the observations mentioned above, this indicates that the Broet study critically influenced the statistical significance of the overall effect, reflecting the poor stability of the pooled effect size (Figure 5).

Figure 5
Forest plot showing meta-analysis estimates when specific studies are omitted. Studies listed are Broet (2001), Fanelli (2021), and Matikas (2021). Circles represent estimates, with horizontal lines indicating confidence intervals. X-axis ranges from 0.67 to 2.02.

Figure 5. Sensitivity analysis of OS in treatment-naive subgroups.

In the sensitivity analysis of OS recurrence subgroup studies, sequential exclusion of each included study showed that the direction of pooled effect sizes in the remaining studies remained consistent, and none of the 95% CIs for the pooled effect sizes included the null value of 1.0. These results indicated that the absence of any individual study did not significantly alter the overall conclusions, suggesting that the meta-analysis findings were not substantially influenced by single studies (Figure 6).

Figure 6
Plot illustrating the meta-analysis estimates with certain studies omitted, showing estimates and confidence intervals for BJohle (2013), Larsson (2020), and Bergqvist (2023). Estimates are marked as circles on a horizontal scale from 1.40 to 2.77, with confidence intervals depicted as horizontal lines.

Figure 6. Sensitivity analysis of OS in recurrence subgroups.

3.5 Publication bias

For the meta-analysis of baseline TK1 expression and PFS involving 11 studies, publication bias was assessed using Begg’s and Egger’s tests. The results showed no significant publication bias (p > 0.05). Additionally, the trim-and-fill method detected no studies requiring trimming or filling. Collectively, these findings indicate that the included studies were not substantially affected by publication bias (Figure 7).

Figure 7
Funnel plot showing points representing data, with pseudo ninety-five percent confidence limits. The x-axis is the standard error of theta, filled, and the y-axis is theta, filled. Points are dispersed within the funnel-shaped confidence limits.

Figure 7. Funnel plot assessing publication bias in PFS studies (posttrim-and-fill adjustment).

Publication bias was assessed for the six studies included in the meta-analysis of baseline TK1 expression and OS. Results from Begg’s and Egger’s tests indicated no significant publication bias (p > 0.05). The trim-and-fill method also detected no studies requiring adjustment. Given the limited number of studies included in the OS analysis, the statistical tests used to detect publication bias may be underpowered. Therefore, the potential for publication bias cannot be entirely ruled out (Figure 8).

Figure 8
Funnel plot showing the standard error of theta (x-axis) against theta filled (y-axis). The plot includes a funnel shape with pseudo 95% confidence limits and several data points distributed within and outside the funnel region.

Figure 8. Funnel plot assessing publication bias in OS studies (posttrim-and-fill adjustment).

4 Discussion

This meta-analysis evaluated the prognostic significance of baseline TK1 expression in breast cancer, with strict inclusion criteria requiring TK1 assessment prior to treatment initiation. Pooled analysis of 11 high-quality studies (NOS scores 8–9) encompassing diverse clinical subtypes (early/advanced stage, various receptor statuses, and operable/inoperable cases) revealed significantly worse progression-free survival in patients with elevated baseline TK1 (HR = 1.63, 95% CI: 1.28–2.10, p < 0.001; random-effects model, I2 = 87.9%).

The significant heterogeneity observed may originate from several potential sources. First, clinical heterogeneity across patient populations likely contributed substantially. The included studies involved a broad spectrum of breast cancer patients, ranging from those with early-stage operable disease (9, 11, 18) to those with advanced metastatic disease (1317, 19), as well as patients receiving neoadjuvant therapy (9, 20). Notable differences exist between early- and advanced-stage patients regarding tumor burden and the tumor microenvironment. For example, TK1 may more directly reflect high proliferative activity and disease burden in advanced-stage patients, whereas in early-stage patients, it could be associated with micrometastatic potential. Although a subgroup analysis based on treatment-naive versus recurrent disease was performed, this stratification did not fully account for the observed heterogeneity(Adj R-squared = 34.03%, p = 0.175), suggesting that additional clinicopathological factors—such as hormone receptor status, HER-2 expression, prior lines of therapy, and metastatic sites—may also influence the results.

Second, methodological variations in detection assays and specimen types represent another key source of heterogeneity. As summarized in Table 1, most studies (nine of 11) measured serum TK1 levels using ELISA. However, considerable differences existed among studies in the commercial kits used and in the definitions of positive thresholds (e.g., median, upper quartile, or prespecified cutoff values). Moreover, one study employed immunohistochemistry (IHC) to assess TK1 protein expression in tumor tissue (18), whereas an earlier study utilized the BCA method (11). Importantly, IHC captures in situ TK1 expression within tumor tissue, whereas ELISA quantifies circulating TK1 in peripheral blood. These two approaches reflect biologically distinct dimensions of TK1, differing in quantitative units, analytical sensitivity, and clinical interpretation—differences that inevitably introduce substantial heterogeneity across studies. Due to the limited overall number of available studies, performing univariate subgroup analyses for heterogeneity factors beyond treatment−naive/recurrent status would further reduce the number of studies in each subgroup, thereby precluding a meaningful examination or reduction of heterogeneity through subgroup comparisons. For these reasons, a random−effects model was ultimately employed to conduct the meta−analysis of PFS.

Analysis of six studies reporting overall survival demonstrated consistent prognostic value across subgroups (treatment-naive: HR = 1.30; recurrence: HR = 2.10; pooled HR = 1.58, 95% CI: 1.40–1.78, p = 0.003), with heterogeneity resolved through stratification by disease status. Based on the observed disparity in HR between the two subgroups, the prognostic significance of TK1 appears to be disease stage-dependent. These findings confirm that baseline TK1 expression is consistently associated with adverse outcomes in breast cancer patients, establishing it as a reliable prognostic indicator.

As a key enzyme in the DNA salvage pathway, TK1 demonstrates peak expression during the S-phase of the cell cycle, serving as a precise indicator of tumor proliferative activity. Its biological properties have been well characterized. Current research on TK1 spans diverse solid and hematological malignancies (2124), with the most clinically significant positive findings concentrated in breast cancer (25). Accumulating evidence has established TK1 as a key mediator of breast cancer progression. Recent mechanistic studies have further defined its specific function as a downstream effector of the cyclin-dependent kinase 4/6 (CDK4/6) pathway. In patients with hormone receptor-positive/human epidermal growth factor receptor 2-negative (HR+/HER2−) advanced breast cancer, an early elevation in TK1 activity following CDK4/6 inhibitor administration is consistently correlated with decreased endocrine therapy sensitivity or the development of acquired drug resistance (26, 27). Our findings substantiate the clinical utility of baseline TK1 levels in prognostic assessment for breast cancer. Importantly, baseline TK1 maintained stable prognostic performance across diverse clinical subtypes, suggesting its potential to overcome the limitations of conventional pathological classification and provide more broadly applicable clinical guidance.

Nevertheless, this study has several limitations that should be acknowledged. The primary limitation stems from substantial heterogeneity in the progression-free survival analysis (I2 = 87.9%), which persisted despite employing disease status stratification and other statistical adjustments. Interpretation of PFS outcomes from this meta-analysis requires caution, as the high heterogeneity observed suggests that the prognostic value of TK1 may be influenced by variations in patient characteristics and detection methodologies. To consolidate the clinical utility of TK1, subsequent investigations are recommended to establish standardized inclusion criteria and uniform TK1 detection protocols, coupled with rigorous prospective validation studies, thereby yielding robust evidence to advance the development and optimization of TK1-driven therapeutic strategies. Second, our analysis may be susceptible to publication bias, in which studies demonstrating a significant association between TK1 expression and poor prognosis are more likely to be published, potentially leading to an overestimation of its true prognostic value. Although an exhaustive literature search was conducted, this limitation cannot be entirely excluded.

Building on our findings, we propose three critical directions for future investigation: (1) establishing standardized TK1 detection protocols to enhance interlaboratory reproducibility and clinical applicability; (2) conducting mechanistic studies to elucidate TK1’s prognostic role in breast cancer, with particular emphasis on its crosstalk with core signaling pathways to inform targeted therapy development; and (3) performing biomarker-driven clinical trials to evaluate TK1-guided therapeutic strategies, especially in patients exhibiting elevated TK1 expression.

An optimal prognostic biomarker should fulfill a tripartite role: outcome prediction, treatment response monitoring, and guidance for therapeutic decision-making. Our meta-analysis demonstrates a significant association between elevated baseline TK1 expression and adverse prognostic outcomes in breast cancer. Although high baseline TK1 expression indicates poor prognosis, whether adjusting treatment strategies in advance or implementing more active interventions can improve clinical benefits for these patients under the current treatment system requires further support from high-quality clinical research data.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

PL: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing – original draft. YC: Data curation, Formal Analysis, Investigation, Project administration, Resources, Supervision, Validation, Writing – review & editing. JZ: Funding acquisition, Resources, Writing – review & editing. YF: Formal Analysis, Validation, Writing – original draft. TZ: Formal Analysis, Visualization, Writing – original draft.

Funding

The author(s) declared that financial support was received for this work and/or its publication. The Key Project of Zhangjiakou Science and Technology Program (Grant No. 2021026D) and the Basic Scientific Research Operating Fund for Provincial Universities (Grant No. JYT2023016).

Acknowledgments

The authors gratefully thank Yonggang Huang and Junfeng Zhao for their support.

Conflict of interest

The author(s) declared that this work 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) declared that generative AI was not 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/fonc.2026.1666576/full#supplementary-material

Supplementary Table 1 | Detailed quality assessment of included cohort studies using the Newcastle-Ottawa Scale (NOS). Scoring Guide: ★: One star awarded for fulfilling the criterion. Selection (S): S1:​ Representativeness of the Exposed Cohort. S2:​ Selection of the Non-Exposed Cohort. S3:​ Ascertainment of Exposure (e.g., TK1 measurement method). S4:​ Demonstration that outcome of interest was not present at start. Comparability (C): C1:​ Comparability of cohorts on the basis of the design or analysis (controlling for confounders). The most important factor(s) controlled for are noted (e.g., age, tumor stage). Outcome (O): O1:​ Assessment of Outcome (e.g., progression, death). O2:​ Was follow-up long enough for outcomes to occur? O3:​ Adequacy of follow-up of cohorts.

References

1. Tarantino P, Viale G, Press MF, Hu X, Penault-Llorca F, Bardia A, et al. ESMO expert consensus statements (ECS) on the definition, diagnosis, and management of HER2-low breast cancer. Ann Oncol. (2023) 34:645–59. doi: 10.1016/j.annonc.2023.05.008

PubMed Abstract | Crossref Full Text | Google Scholar

2. Wolff AC, Somerfield MR, Dowsett M, Hammond M, Hayes DF, McShane LM, et al. Human epidermal growth factor receptor 2 testing in breast cancer: ASCO-college of American pathologists guideline update. J Clin Oncol. (2023) 41:3867–72. doi: 10.1200/JCO.22.02864

PubMed Abstract | Crossref Full Text | Google Scholar

3. Clusan L, Ferriere F, Flouriot G, and Pakdel F. A basic review on estrogen receptor signaling pathways in breast cancer. Int J Mol Sci. (2023) 24:6834. doi: 10.3390/ijms24076834

PubMed Abstract | Crossref Full Text | Google Scholar

4. Chen YL, Eriksson S, and Chang ZF. Regulation and functional contribution of thymidine kinase 1 in repair of DNA damage. J Biol Chem. (2010) 285:27327–35. doi: 10.1074/jbc.M110.137042

PubMed Abstract | Crossref Full Text | Google Scholar

5. Ke PY and Chang ZF. Mitotic degradation of human thymidine kinase 1 is dependent on the anaphase-promoting complex/cyclosome-CDH1-mediated pathway. Mol Cell Biol. (2004) 24:514–26. doi: 10.1128/MCB.24.2.514-526.2004

PubMed Abstract | Crossref Full Text | Google Scholar

6. He Q, Zhang P, Zou L, Li H, Wang X, Zhou S, et al. Concentration of thymidine kinase 1 in serum (S-TK1) is a more sensitive proliferation marker in human solid tumors than its activity. Oncol Rep. (2005) 14:1013–9. doi: 10.3892/or.14.4.1013

PubMed Abstract | Crossref Full Text | Google Scholar

7. Paoletti C, Barlow WE, Cobain EF, Bergqvist M, Mehta RS, Gralow JR, et al. Evaluating serum thymidine kinase 1 in patients with hormone receptor-positive metastatic breast cancer receiving first-line endocrine therapy in the SWOG S0226 trial. Clin Cancer Res. (2021) 27:6115–23. doi: 10.1158/1078-0432.CCR-21-1562

PubMed Abstract | Crossref Full Text | Google Scholar

8. Malorni L, Bianchini G, Caputo R, Zambelli A, Puglisi F, Bianchi GV, et al. Serum thymidine kinase activity in patients with HR-positive/HER2-negative advanced breast cancer treated with ribociclib plus letrozole: Results from the prospective BioItaLEE trial. Eur J Cancer. (2023) 186:1–11. doi: 10.1016/j.ejca.2023.03.001

PubMed Abstract | Crossref Full Text | Google Scholar

9. Matikas A, Wang K, Lagoudaki E, Acs B, Zerdes I, Hartman J, et al. Prognostic role of serum thymidine kinase 1 kinetics during neoadjuvant chemotherapy for early breast cancer. ESMO Open. (2021) 6:100076. doi: 10.1016/j.esmoop.2021.100076

PubMed Abstract | Crossref Full Text | Google Scholar

10. Altman DG and Bland JM. How to obtain the confidence interval from a P value. Bmj-Brit Med J. (2011) 343:d2090. doi: 10.1136/bmj.d2090

PubMed Abstract | Crossref Full Text | Google Scholar

11. Broet P, Romain S, Daver A, Ricolleau G, Quillien V, Rallet A, et al. Thymidine kinase as a proliferative marker: clinical relevance in 1,692 primary breast cancer patients. J Clin Oncol. (2001) 19:2778–87. doi: 10.1200/JCO.2001.19.11.2778

PubMed Abstract | Crossref Full Text | Google Scholar

12. Nisman B, Allweis T, Kaduri L, Maly B, Gronowitz S, Hamburger T, et al. Serum thymidine kinase 1 activity in breast cancer. Cancer biomark. (2010) 7:65–72. doi: 10.3233/CBM-2010-0148

PubMed Abstract | Crossref Full Text | Google Scholar

13. Bjohle J, Bergqvist J, Gronowitz JS, Johansson H, Carlsson L, Einbeigi Z, et al. Serum thymidine kinase activity compared with CA 15–3 in locally advanced and metastatic breast cancer within a randomized trial. Breast Cancer Res TR. (2013) 139:751–8. doi: 10.1007/s10549-013-2579-x

PubMed Abstract | Crossref Full Text | Google Scholar

14. Bonechi M, Galardi F, Biagioni C, De Luca F, Bergqvist M, Neumuller M, et al. Plasma thymidine kinase-1 activity predicts outcome in patients with hormone receptor positive and HER2 negative metastatic breast cancer treated with endocrine therapy. Oncotarget. (2018) 9:16389–99. doi: 10.18632/oncotarget.24700

PubMed Abstract | Crossref Full Text | Google Scholar

15. McCartney A, Biagioni C, Schiavon G, Bergqvist M, Mattsson K, Migliaccio I, et al. Prognostic role of serum thymidine kinase 1 activity in patients with hormone receptor-positive metastatic breast cancer: Analysis of the randomised phase III Evaluation of Faslodex versus Exemestane Clinical Trial (EFECT). Eur J Cancer. (2019) 114:55–66. doi: 10.1016/j.ejca.2019.04.002

PubMed Abstract | Crossref Full Text | Google Scholar

16. Larsson AM, Bendahl PO, Aaltonen K, Jansson S, Forsare C, Bergqvist M, et al. Serial evaluation of serum thymidine kinase activity is prognostic in women with newly diagnosed metastatic breast cancer. Sci Rep-UK. (2020) 10:4484. doi: 10.1038/s41598-020-61416-1

PubMed Abstract | Crossref Full Text | Google Scholar

17. McCartney A, Bonechi M, De Luca F, Biagioni C, Curigliano G, Moretti E, et al. Plasma thymidine kinase activity as a biomarker in patients with luminal metastatic breast cancer treated with palbociclib within the TREnd trial. Clin Cancer Res. (2020) 26:2131–9. doi: 10.1158/1078-0432.CCR-19-3271

PubMed Abstract | Crossref Full Text | Google Scholar

18. Fanelli GN, Scarpitta R, Cinacchi P, Fuochi B, Szumera-Cieckiewicz A, De Ieso K, et al. Immunohistochemistry for thymidine kinase-1 (TK1): A potential tool for the prognostic stratification of breast cancer patients. J Clin Med. (2021) 10:5416. doi: 10.3390/jcm10225416

PubMed Abstract | Crossref Full Text | Google Scholar

19. Bergqvist M, Nordmark A, Williams A, Paoletti C, Barlow W, Cobain EF, et al. Thymidine kinase activity levels in serum can identify HR+ metastatic breast cancer patients with a low risk of early progression (SWOG S0226). Biomarkers. (2023) 28:313–22. doi: 10.1080/1354750X.2023.2168063

PubMed Abstract | Crossref Full Text | Google Scholar

20. Zhu Y, Zerdes I, Matikas A, Cruz IR, Bergqvist M, Elinder E, et al. The role of serum thymidine kinase 1 activity in neoadjuvant-treated HER2-positive breast cancer: biomarker analysis from the Swedish phase II randomized PREDIX HER2 trial. Breast Cancer Res TR. (2024) 204:299–308. doi: 10.1007/s10549-023-07200-x

PubMed Abstract | Crossref Full Text | Google Scholar

21. Gasparri F, Wang N, Skog S, Galvani A, and Eriksson S. Thymidine kinase 1 expression defines an activated G1 state of the cell cycle as revealed with site-specific antibodies and ArrayScan assays. Eur J Cell Biol. (2009) 88:779–85. doi: 10.1016/j.ejcb.2009.06.005

PubMed Abstract | Crossref Full Text | Google Scholar

22. Nisman B, Appelbaum L, Yutkin V, Nechushtan H, Hubert A, Uziely B, et al. Serum thymidine kinase 1 activity following nephrectomy for renal cell carcinoma and radiofrequency ablation of metastases to lung and liver. Anticancer Res. (2016) 36:1791–7. doi: 10.21873/anticanres.10602

PubMed Abstract | Crossref Full Text | Google Scholar

23. Chen G, He C, Li L, Lin A, Zheng X, He E, et al. Nuclear TK1 expression is an independent prognostic factor for survival in pre-malignant and Malignant lesions of the cervix. BMC Cancer. (2013) 13:249. doi: 10.1186/1471-2407-13-249

PubMed Abstract | Crossref Full Text | Google Scholar

24. Tao J, Wang Z, Shi R, Lin L, Li M, Meng Y, et al. ERK-USP9X-coupled regulation of thymidine kinase 1 promotes both its enzyme activity-dependent and its enzyme activity-independent functions for tumor growth. Nat Struct Mol Biol. (2025) 32:853–63. doi: 10.1038/s41594-024-01473-6

PubMed Abstract | Crossref Full Text | Google Scholar

25. McCartney A and Malorni L. Thymidine kinase-1 as a biomarker in breast cancer: estimating prognosis and early recognition of treatment resistance. biomark Med. (2020) 14:495–8. doi: 10.2217/bmm-2020-0072

PubMed Abstract | Crossref Full Text | Google Scholar

26. Stender JD, Frasor J, Komm B, Chang KC, Kraus WL, and Katzenellenbogen BS. Estrogen-regulated gene networks in human breast cancer cells: involvement of E2F1 in the regulation of cell proliferation. Mol Endocrinol. (2007) 21:2112–23. doi: 10.1210/me.2006-0474

PubMed Abstract | Crossref Full Text | Google Scholar

27. Jia X, Wang K, Wang X, Luo Y, Liu H, Zhao Z, et al. Options after progression on first-line CDK4/6 inhibitors in advanced breast cancer patients. Ther Adv Med Oncol. (2025) 17:22781359. doi: 10.1177/17588359251336623

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: baseline TK1 expression, breast cancer, meta-analysis, OS, PFS, thymidine kinase 1

Citation: Li P, Cheng Y, Zhao J, Fang Y and Zhao T (2026) Elevated thymidine kinase 1 expression at baseline predicts poor prognosis in breast cancer patients. Front. Oncol. 16:1666576. doi: 10.3389/fonc.2026.1666576

Received: 15 July 2025; Accepted: 12 January 2026; Revised: 17 December 2025;
Published: 29 January 2026.

Edited by:

Haiyan Li, The Sixth Affiliated Hospital of Sun Yat-sen University, China

Reviewed by:

Henrik Rönnberg, Swedish University of Agricultural Sciences, Sweden
Vijay Paul Samuel, RAK Medical and Health Sciences University, United Arab Emirates

Copyright © 2026 Li, Cheng, Zhao, Fang and Zhao. 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: Yongting Cheng, Y3l0aGJudUAxNjMuY29t

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.