- 1Medical School of Chinese PLA, Beijing, China
- 2Department of General Surgery, The First Medical Center, Chinese PLA General Hospital., Beijing, China
- 3School of Medicine, Nankai University, Tianjin, China
- 4Department of General Surgery, Beijing Friendship Hospital, Capital Medical University, Beijing, China
Background: Papillary thyroid cancer (PTC), the predominant histologic subtype of thyroid cancer cases, has increased substantially over the past decades. In previous studies, Hashimoto’s thyroiditis (HT) exerts a paradoxical dual role in PTC. However, limited studies have specifically examined the association between HT and the invasion of PTC.
Methods: In this retrospective study, 10329 PTC patients were selected, and the clinicopathological features were retrospectively analyzed. Propensity score matching (PSM) was employed to minimize confounding effects from baseline variables. Univariate analysis and multivariate analysis were performed using binary logistic regression to determine the predictive factor. Odds ratio (OR) and 95% confidence interval (CI) were calculated.
Results: Among 10329 PTC patients, 992 (9.6%) individuals were diagnosed with HT. Compared to the non-HT group, the HT group demonstrated lower rates of extrathyroidal extension (p<0.001), reduced multifocality (p=0.004), decreased bilateral carcinoma involvement (p<0.001), and a greater proportion of pathological N0 stage tumors (p=0.003). Following PSM, a cohort of 970 patients with HT and 2783 non-HT controls without HT were analyzed. HT was significantly associated with lower rates of: central lymph node metastasis (CLNM) in cN0 patients (p<0.001); lateral lymph node metastasis (LLNM) in cN1b patients (p=0.008); and bilateral carcinoma detection in patients with clinically unilateral PTC lesions (p=0.001).
Conclusion: This study found an association between HT and reduced invasiveness of PTC, as evidenced by increased node-negative disease and reduced bilateral pathological involvement.
Introduction
The incidence of thyroid cancer (TC) has increased substantially over the past decades, ranking as the seventh most prevalent cancer globally and the fifth most prevalent malignancy in females (1–3). Papillary thyroid cancer (PTC), the predominant histologic subtype of TC cases, represents the primary driver of this observed epidemiological rise (4). Updated in 2025, the American Thyroid Association (ATA) Management Guidelines for Adult Patients with Thyroid Nodules and Differentiated Thyroid Cancer further mandate precise risk stratification and personalized treatment planning for PTC patients (5).
Hashimoto’s thyroiditis (HT), one of the most prevalent autoimmune thyroid disorders, arises from complex interactions among genetic susceptibility, environmental triggers, and epigenetic modifications (6). It is serologically characterized by elevated anti-thyroglobulin (TgAb) and anti-thyroid peroxidase (TPOAb) antibodies. Histopathologically, HT manifests as lymphoplasmacytic infiltration, germinal center formation within lymphoid follicles, and progressive parenchymal atrophy (7). The relationship between HT and PTC remains contentious. On the one hand, chronic inflammation-driven oncogenic transformation and TSH-mediated follicular epithelial proliferation may promote PTC development (8, 9). Conversely, Anil et al. prospectively demonstrated no increased risk in HT patients with thyroid nodules (10). However, few studies have examined the association between HT and PTC invasion.
In this study, we systematically evaluated the association between clinicopathological characteristics and HT status in PTC patients. Propensity score matching (PSM) was employed to minimize confounding effects from baseline variables. Primary endpoints included lymph node metastasis status and incidence of bilateral carcinomas (initially diagnosed as unilateral disease preoperatively). Our comprehensive analysis provides clinically actionable evidence for refining risk-adapted management strategies in HT-associated PTC.
Materials and methods
Patients’ selection
Consecutive PTC patients who underwent surgical resection for treatment at the First Medical Center of the Chinese People’s Liberation Army (PLA) General Hospital from January 2020 to February 2025 were screened. Related patients meeting the following criteria were selected: (1) an obvious pathologic diagnosis of PTC; (2) primary cases without a history of thyroid tumors; (3) underwent surgical resection. The exclusion criteria were as follows: (1) patients with a history of thyroid tumors; (2) patients with other types of thyroid tumors, including follicular thyroid cancer, medullary thyroid cancer, and anaplastic thyroid cancer; (3) patients with evidence of distant metastatic disease; (4) patients who received ablation of thyroid nodules. This study was conducted in accordance with the Declaration of Helsinki and approved by the Medical Ethics Committee of the First Medical Center of the Chinese PLA General Hospital.
Data collection
The preoperative and postoperative clinicopathological features were retrospectively studied. Baseline variables, such as patient age and gender, were obtained from medical records. The clinical N stage and tumor location were determined from preoperative ultrasound and computed tomography scans. The surgical strategy for PTC, including thyroid lobectomy/total thyroidectomy, central compartment neck dissection, and lateral neck dissection (for cN1b patients), was determined by the surgeon based on preoperative staging. The presence of HT was diagnosed postoperatively by histopathological examination of paraffin-embedded sections, characterized by diffuse lymphoplasmacytic infiltration, and the presence of lymphoid follicles with reactive germinal centers in the background thyroid tissue, irrespective of preoperative antibody status (11). Tumor size, extrathyroidal extension, multifocal carcinoma, bilateral carcinoma, central lymph node metastasis (CLNM) count, central lymph node dissection (CLND) yield, pathological N stage (pN), lateral lymph node metastasis (LLNM) count, and lateral lymph node dissection (LLND) yield were determined through postoperative histopathological examination.
Statistical analysis
IBM SPSS Statistics 26.0 software (SPSS Inc, Chicago, IL, USA) and R software (ver. 4.2.1, R Development Core Team) were adopted for statistical analyses in this study. Numerical data with normal distribution were expressed as the mean ± standard deviation (SD), and numerical data without normal distribution were expressed as median (interquartile range, IQR). Categorical data are presented as absolute numbers and percentages. The MatchIt package was used to conduct the PSM, with a matching ratio of 1:3 between the HT and non-HT groups and a matching capacity of 0.02 (12). The standardized mean differences (SMD) of covariate balance before and after PSM were subsequently calculated. Univariate and multivariate analyses were performed using binary logistic regression to identify the predictive factors for CLNM, LLNM, and bilateral cancer. Odds ratio (OR) and 95% confidence interval (CI) were calculated. P-values were two-tailed, and P < 0.05 was considered statistically significant. Data visualization was performed using GraphPad Prism 8 Software.
Results
Patient selection and characteristics
Initially, a cohort of 10663 PTC patients who underwent thyroidectomy was included. Of these, 151 patients had a history of PTCs, 130 patients had other types of thyroid carcinoma, 2 patients had evidence of distant metastatic disease, and 51 patients had received ablation of thyroid nodules. After rigorous screening, 10329 patients were selected for further analysis, and 992 (9.6%) were diagnosed with HT (Figure 1). Compared with the non-HT group, the HT group exhibited a significantly younger age (p<0.001), greater female predominance (p<0.001), higher rates of clinical N1a and N1b disease (p<0.001), and more frequent bilateral thyroidectomy (p<0.001). Conversely, the HT group demonstrated lower rates of extrathyroidal extension (p<0.001), reduced multifocality (p=0.004), and decreased bilateral carcinoma involvement (p<0.001). Additionally, the HT group had higher numbers of CLND (p<0.001) and a greater proportion of pathological N0 stage tumors (p=0.003) (Table 1; Figure 2A).
Figure 1. The flow chart for patient selection PTC, papillary thyroid carcinoma; HT, Hashimoto’s thyroiditis; PSM, propensity score matching.
Figure 2. Comparison of clinical characteristics in papillary thyroid carcinoma patients with versus without Hashimoto’s thyroiditis before and after propensity score matching. (A) Pre-matching cohort; (B) Post-matching cohort. CLNM, central lymph node metastasis, CLND, central lymph node dissection, LLNM, lateral lymph node metastasis, LLND, lateral lymph node dissection.
Covariates, including age, gender, clinical N stage, clinical tumor location, surgical extent, tumor size, extrathyroidal extension, and multifocal carcinoma, were well-balanced through PSM. Following a 1:3 ratio PSM, the final cohort comprised 970 patients with HT and 2783 non-HT controls. All covariates demonstrated non-significant intergroup differences (P > 0.05; Figure 2B; Supplementary Table 1) with SMD below the 0.1 threshold (Table 2; Figure 3). The propensity score distribution and matching methodology are detailed in Supplementary Figures 2, 3.
Figure 3. The assessment of standardized mean differences (SMD) for covariate balance before and after propensity score matching (PSM).
HT independently associates with lower CLNM in cN0 PTC
After propensity score matching, 3753 patients were included in the final analysis. Among them, 2859 patients (76.2%) were clinically node-negative (cN0), of whom 809 (28.3%) had pathologically confirmed CLNM. Univariate analysis identified male gender (p<0.001; OR = 1.951), right-sided tumors (p=0.024; OR = 1.235), bilateral tumors (p=0.014; OR = 1.332), larger tumor size (p < 0.001), extrathyroidal extension (p<0.001; OR = 1.597), multifocal carcinoma (p=0.011; OR = 1.242), and bilateral carcinoma (p=0.012; OR = 1.264) as potential risk factors for CLNM. In contrast, increasing age (p<0.001; OR = 0.959) and HT (p < 0.001; OR = 0.646) showed potential protective effects. Furthermore, a multivariate analysis was conducted using binary logistic regression to identify the predictive factors for CLNM in cN0 PTC patients (Table 3; Figure 4A). Male gender (p<0.001; OR = 1.698), right-sided tumors (p=0.011; OR = 1.285), bilateral tumors (p=0.032; OR = 1.309), larger tumor size (p<0.001), and extrathyroidal extension (p=0.018; OR = 1.272) as significant risk factors for CLNM, while increasing age (p<0.001; OR = 0.957) and HT (p<0.001; OR = 0.609) were significant associated with lower CLNM rate in cN0 PTC.
Figure 4. Univariate and multivariate logistic regression analyses of risk factors for: (A) CLNM in cN0 PTC patients; (B) CLNM in cN1a PTC patients; (C) LLNM in cN1b PTC patients; (D) bilateral carcinoma identified in clinically unilateral PTC lesions. CLNM, central lymph node metastasis, LLNM, lateral lymph node metastasis, PTC, papillary thyroid cancer.
Among 344 cN1aPTC patients, 291 (84.6%) had pathologically confirmed CLNM. After univariate and multivariate analysis, right-sided tumors (p=0.019; OR = 2.297), larger tumor size (p=0.027), and extrathyroidal extension (p=0.014; OR = 2.509) were identified as significant risk factors for CLNM, while increasing age (p=0.001; OR = 0.952) was a significant protective factor (Table 4; Figure 4B).
HT independently associates with lower LLNM in CN1b PTC
Of the 550 cN1b PTC patients studied, 422 (76.7%) developed pathologically proven LLNM, among whom 54 (12.8%) presented with skip metastasis (Table 5; Figure 4C). Univariate analysis revealed significant associations for LLNM risk: tumor size 1–2 cm (p<0.001, OR = 3.381); tumor size 2–4 cm (p<0.001, OR = 6.679); extrathyroidal extension (p<0.001, OR = 2.611); and CLNM (p<0.001, OR = 8.225). Protective effects were observed for age (p<0.001, OR = 0.952) and HT (p=0.014, OR = 0.593). Multivariate analysis confirmed independent predictors: tumor size 1–2 cm (p=0.047, OR = 2.171); tumor size 2–4 cm (p=0.006, OR = 3.901); CLNM (p<0.001, OR = 5.771); with persistent protection from age (p<0.001, OR = 0.965).HT (p=0.008, OR = 0.514) was significantly associated with lower LLNM rate in cN1b PTC.
HT independently associates with reduced bilateral pathological involvement
In our cohort, 1574 patients had clinically unilateral lesions and underwent bilateral thyroidectomy, involving 1361 cN0 patients and 213 cN1a patients. Among them, 256 patients (16.2%) were pathologically confirmed to have bilateral thyroid carcinoma (Table 6; Figure 4D). Univariate analysis identified cN1a metastasis (p=0.008, OR = 1.613), right-sided tumors (p=0.015, OR = 1.408), tumor size >4cm (p=0.025, OR = 3.866), extrathyroidal extension (p=0.018, OR = 1.443), unilateral multifocal carcinoma (p<0.001, OR = 2.122), and CLNM (p=0.028, OR = 1.360) as potential risk factors for bilateral carcinoma, while HT showed a protective association (p=0.004, OR = 0.614). Multivariate analysis confirmed cN1a (p=0.014, OR = 1.574), right-sided location (p=0.021, OR = 1.389), tumor size >4cm (p=0.031, OR = 3.668), and unilateral multifocal carcinoma (p<0.001, OR = 2.130) as independent risk factors, with HT remaining a significant negative correlate (p=0.001, OR = 0.571).
Table 6. The risk factors for bilateral carcinoma in clinical unilateral lesion PTC patients (N=1574).
Discussion
The imperative for precision risk stratification in PTC management is substantial, which is complicated by concomitant HT. This complexity arises not only from HT-induced diagnostic pitfalls in fine-needle aspiration biopsy (FNAB), but more critically from the paradoxical oncologic implications of HT itself (13, 14). In our cohort, HT coexisted with PTC in 9.6% of cases, demonstrating significant female predominance, consistent with prior epidemiological reports (11, 15–17). Crucially, the presence of HT prompted more aggressive preoperative management, evidenced by significantly higher rates of clinically suspected nodal metastases (cN1a and cN1b stage) and consequent performance of bilateral thyroidectomy. However, postoperative pathological assessment revealed a paradoxical protective association: HT was independently associated with reduced tumor invasiveness, with significantly lower frequencies of extrathyroidal extension, multifocality, and bilateral carcinoma involvement, and a higher proportion of pN0 stage tumors. After PSM to control for key confounders, HT consistently demonstrated an association with reduced PTC invasiveness. This was evidenced by significantly reduced CLNM rates in cN0 patients, lower incidence of LLNM in cN1b cases, and decreased bilateral tumor detection among preoperatively diagnosed unilateral tumors. These robust findings establish a significant association between HT and diminished invasive potential of PTC.
The protective association between HT and attenuated PTC invasiveness observed in our PSM analysis aligns with robust existing evidence. A multicenter cohort of 9210 PTC patients demonstrated HT’s negative correlation with aggressive features, including primary tumors larger than 4 cm, gross extrathyroidal extension, extranodal extension, and distant metastasis, and concurrently correlated with superior 10-year disease-specific and recurrence-free survival (11). This is further substantiated by a meta-analysis of 10648 cases across 38 studies, which found that HT co-occurrence significantly predicted reduced extrathyroidal extension, fewer positive lymph nodes, and improved long-term recurrence-free survival (17). Single-institutional data from 435 surgically treated PTC patients corroborated these findings, showing that PTC-HT cases had smaller tumor dimensions, less lymph node involvement, and earlier-stage disease than HT-negative counterparts (18). Collectively, these findings reinforce HT’s role in contributing to its less aggressive clinicopathological presentation in PTC.
Notwithstanding the predominant protective patterns, certain studies report conflicting associations between Hashimoto’s thyroiditis and papillary thyroid carcinoma behavior. Lee et al.’s aforementioned meta-analysis paradoxically identified HT as a risk factor for multifocal tumor involvement (17), a finding corroborated by Cappellacci et al. in their cohort of 839 PTC patients, in which the HT group demonstrated significantly higher multifocal rates (40.5% vs. 29.6%) (19). Although the latter study confirmed HT’s association with smaller tumor diameters, their multivariate analysis unexpectedly established HT as an independent risk factor for thyroid cancer development. This aligns with Apostolou et al.’s investigation of 3233 patients with multinodular thyroid disease, in which HT independently predicted malignant histopathology (20). Most strikingly, Chen et al.’s cohort demonstrated an 11.8-fold increased thyroid cancer risk in HT patients after full adjustment for demographic and comorbidity confounders (21). Danis et al. similarly reported that HT was associated with elevated PTC risk and consequently recommended more aggressive diagnostic approaches for HT patients with thyroid nodules (22). These contradictory observations suggest HT’s dual role in thyroid carcinogenesis, potentially promoting initial malignant transformation while simultaneously constraining progression through immunomodulatory mechanisms.
Advancements in research methodologies have progressively elucidated the molecular mechanisms underlying HT’s influence on PTC pathogenesis. In vitro experiments by Kim et al. demonstrated that HT attenuates PTC cell invasiveness by upregulated of E-cadherin and TGF-β, key regulators of epithelial-mesenchymal transition (23). Complementarily, Sun et al.’s metabolomic profiling of PTC patients revealed dysregulated aminoacyl-tRNA biosynthesis and serine/threonine metabolism in HT-associated tumors, suggesting metabolic reprogramming as a contributory mechanism (24). Multi-omics analyses further identified HT-specific genetic signatures via The Cancer Genome Atlas, including DMBT1, MET, FAM20A, SERPINA1, CD53, FCER1G, and TYROBP (25–28). Most notably, Ma et al.’s single-cell RNA sequencing identified HT-associated cell clusters that establish a TSH-suppressive microenvironment through interactions among mTE3, nTE0, and nTE2 subpopulations. This immunomodulatory network is functionally mediated by the MIF-(CD74+CXCR4) signaling axis, facilitating stromal-immune communication that may concurrently promote carcinogenesis yet restrain progression (29).
Notably, our study identified preoperative misdiagnosis of contralateral PTC in 16.2% of presumed unilateral cases, with HT emerging as a negative correlate of bilateral involvement, a finding that diverged from the prevailing literature. This contrasts with Lv et al.’s analysis of 1442 patients undergoing total thyroidectomy for unilateral PTC with benign-appearing contralateral nodules, where HT was classified as a risk factor alongside multifocality and capsular invasion (30). Similarly, Wang et al.’s meta-analysis (23% contralateral malignancy misdiagnosis rate) implicated HT as a significant risk predictor (31). Conversely, Wu et al. reported no association between HT and contralateral cancer in their cohort (28.9% misdiagnosis rate) (32). Our large-scale propensity-matched cohort demonstrates HT’s negative correlation against contralateral malignancies, aligning with its established inhibitory roles in protection against CLNM in cN0 patients and LLNM in cN1b patients.
While this study demonstrates an inverse association between HT and markers of PTC aggressiveness, several limitations warrant consideration. Firstly, the retrospective nature inherently introduces selection bias. Secondly, although PSM mitigated major confounders, residual bias persists from unmeasured factors. Besides, the cohort’s restriction to established PTC cases precludes assessment of HT’s impact on de novo carcinogenesis risk, which requires prospective comparison between HT and non-HT populations. Further, our histopathology-based HT diagnosis, while the gold-standard, limits clinical practice where serological criteria (TPOAb/TgAb) or ultrasound features guide preoperative decisions. This methodological discordance may partially explain divergent findings in the literature regarding HT’s dual roles. Finally, while the results may suggest a potentially protective immune microenvironment in HT, the hypothesis remains unproven and requires validation by future molecular and immunopathological studies.
Conclusion
This study demonstrates an inverse association between HT and markers of PTC aggressiveness, evidenced by reduced extrathyroidal extension, multifocality, and bilateral involvement alongside increased node-negative disease. After rigorous PSM, HT demonstrated inverse associations across key clinical scenarios: lowering CLNM rates in clinically node-negative patients, decreasing LLNM in cN1b subgroups, and reducing bilateral tumor detection among preoperatively diagnosed unilateral lesions.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Medical Ethics Committee of the First Medical Center of the Chinese PLA General Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
RL: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing. XS: Conceptualization, Data curation, Methodology, Writing – original draft, Writing – review & editing. LH: Conceptualization, Data curation, Formal analysis, Project administration, Writing – review & editing. ZY: Investigation, Project administration, Writing – review & editing. NL: Investigation, Methodology, Writing – review & editing. PL: Investigation, Project administration, Resources, Writing – review & editing. XZ: Conceptualization, Project administration, Writing – review & editing. WL: Project administration, Supervision, Visualization, Writing – review & editing.
Funding
The author(s) declared that financial support was not received for this work and/or its publication.
Acknowledgments
We thank all the patients whose data were used for the study. Our greatest acknowledgement goes to the authors who made detailed data available for this study and to all our colleagues in this study for their hard work.
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/fimmu.2025.1745452/full#supplementary-material
Abbreviations
TC, thyroid cancer; PTC, papillary thyroid cancer; ATA, American Thyroid Association; HT, Hashimoto’s thyroiditis; TgAb, anti-thyroglobulin antibodies; TPOAb, anti-thyroid peroxidase antibodies; PSM, Propensity score matching; PLA, People’s Liberation Army; CLNM, central lymph node metastasis; CLND, central lymph node dissection; LLNM, lateral lymph node metastasis; LLND, lateral lymph node dissection; SD, standard deviation; IQR, interquartile range; SMD, standardized mean differences; OR, odds ratio; CI, confidence interval; FNAB. fine-needle aspiration biopsy.
References
1. Zhou T, Wang X, Zhang J, Zhou E, Xu C, Shen Y, et al. Global burden of thyroid cancer from 1990 to 2021: a systematic analysis from the Global Burden of Disease Study 2021. J Hematol Oncol. (2024) 17:74. doi: 10.1186/s13045-024-01593-y
2. Li M, Dal Maso L, Pizzato M, and Vaccarella S. Evolving epidemiological patterns of thyroid cancer and estimates of overdiagnosis in 2013–17 in 63 countries worldwide: a population-based study. Lancet Diabetes Endocrinol. (2024) 12:824–36. doi: 10.1016/S2213-8587(24)00223-7
3. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834
4. Filetti S, Durante C, Hartl D, Leboulleux S, Locati LD, Newbold K, et al. Electronic address:Y2xpbmljYWxndWlkZWxpbmVzQGVzbW8ub3Jn. Thyroid cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up†. Ann Oncol. (2019) 30:1856–83. doi: 10.1093/annonc/mdz400
5. Ringel MD, Sosa JA, Baloch Z, Bischoff L, Bloom G, Brent GA, et al. 2025 American thyroid association management guidelines for adult patients with differentiated thyroid cancer. Thyroid. (2025) 35:841–985. doi: 10.1177/10507256251363120
6. Ragusa F, Fallahi P, Elia G, Gonnella D, Paparo SR, Giusti C, et al. Hashimotos’ thyroiditis: Epidemiology, pathogenesis, clinic and therapy. Best Pract Res Clin Endocrinol Metab. (2019) 33:101367. doi: 10.1016/j.beem.2019.101367
7. Ralli M, Angeletti D, Fiore M, D’Aguanno V, Lambiase A, Artico M, et al. Hashimoto’s thyroiditis: An update on pathogenic mechanisms, diagnostic protocols, therapeutic strategies, and potential Malignant transformation. Autoimmun Rev. (2020) 19:102649. doi: 10.1016/j.autrev.2020.102649
8. Dailey M, Lindsay S, and Skahen R. Relation of thyroid neoplasms to Hashimoto disease of the thyroid gland. AMA Arch Surg. (1955) 70:291–7. doi: 10.1001/archsurg.1955.01270080137023
9. Mukasa K, Noh JY, Kunii Y, Matsumoto M, Sato S, Yasuda S, et al. Prevalence of Malignant tumors and adenomatous lesions detected by ultrasonographic screening in patients with autoimmune thyroid diseases. Thyroid. (2011) 21:37–41. doi: 10.1089/thy.2010.0050
10. Anil C, Goksel S, and Gursoy A. Hashimoto’s thyroiditis is not associated with increased risk of thyroid cancer in patients with thyroid nodules: a single-center prospective study. Thyroid. (2010) 20:601–6. doi: 10.1089/thy.2009.0450
11. Xu S, Huang H, Qian J, Liu Y, Huang Y, Wang X, et al. Prevalence of hashimoto thyroiditis in adults with papillary thyroid cancer and its association with cancer recurrence and outcomes. JAMA Netw Open. (2021) 4:e2118526. doi: 10.1001/jamanetworkopen.2021.18526
12. Austin PC. A comparison of 12 algorithms for matching on the propensity score. Stat Med. (2014) 33:1057–69. doi: 10.1002/sim.6004
13. Zhu Y, Song Y, Xu G, Fan Z, and Ren W. Causes of misdiagnoses by thyroid fine-needle aspiration cytology (FNAC): our experience and a systematic review. Diagn Pathol. (2020) 15:1. doi: 10.1186/s13000-019-0924-z
14. Hu F, Yan Z, Ma B, Jiang Y, and Huang H. The impact of concurrent Hashimoto thyroiditis on thyroid nodule cytopathology assessed by ultrasound-guided fine-needle aspiration cytology. Postgrad Med. (2020) 132:506–11. doi: 10.1080/00325481.2020.1739462
15. Ferrari SM, Fallahi P, Elia G, Ragusa F, Ruffilli I, Paparo SR, et al. Thyroid autoimmune disorders and cancer. Semin Cancer Biol. (2020) 64:135–46. doi: 10.1016/j.semcancer.2019.05.019
16. Caturegli P, De Remigis A, Chuang K, Dembele M, Iwama A, and Iwama S. Hashimoto’s thyroiditis: celebrating the centennial through the lens of the Johns Hopkins hospital surgical pathology records. Thyroid. (2013) 23:142–50. doi: 10.1089/thy.2012.0554
17. Lee JH, Kim Y, Choi JW, and Kim YS. The association between papillary thyroid carcinoma and histologically proven Hashimoto’s thyroiditis: a meta-analysis. Eur J Endocrinol. (2013) 168:343–9. doi: 10.1530/EJE-12-0903
18. Battistella E, Pomba L, Costantini A, Scapinello A, and Toniato A. Hashimoto’s thyroiditis and papillary cancer thyroid coexistence exerts a protective effect: a single centre experience. Indian J Surg Oncol. (2022) 13:164–8. doi: 10.1007/s13193-022-01515-9
19. Cappellacci F, Canu GL, Lai ML, Lori E, Biancu M, Boi F, et al. Association between hashimoto thyroiditis and differentiated thyroid cancer: A single-center experience. Front Oncol. (2022) 12:959595. doi: 10.3389/fonc.2022.959595
20. Apostolou K, Zivaljevic V, Tausanovic K, Zoric G, Chelidonis G, Slijepcevic N, et al. Prevalence and risk factors for thyroid cancer in patients with multinodular goitre. BJS Open. (2021) 5:zraa014. doi: 10.1093/bjsopen/zraa014
21. Chen YK, Lin CL, Cheng FT, Sung FC, and Kao CH. Cancer risk in patients with Hashimoto’s thyroiditis: a nationwide cohort study. Br J Cancer. (2013) 109:2496–501. doi: 10.1038/bjc.2013.597
22. Danis N, Comlekci A, Yener S, Durak M, Calan M, Solmaz D, et al. Association between hashimoto’s thyroiditis and papillary thyroid cancer: A single center experience. Acta Endocrinol (Buchar). (2022) 18:74–8. doi: 10.4183/aeb.2022.74
23. Kim SJ, Lee SE, Kim YI, Nam-Goong IS, Jung HW, and Kim ES. Papillary thyroid cancer with Hashimoto’s thyroiditis attenuates the tumour aggressiveness through the up-regulation of E-cadherin and TGF-β expression. Clin Exp Med. (2023) 23:833–40. doi: 10.1007/s10238-022-00857-6
24. Sun D, Zhang Y, Wang D, Zhao X, Han R, Li N, et al. Experimental study on changes in metabolic mechanism of papillary thyroid carcinoma complicated with Hashimoto’s thyroiditis. Heliyon. (2023) 9:e20661. doi: 10.1016/j.heliyon.2023.e20661
25. Gan XX, Li YY, Li SJ, Mo SS, Feng JH, Shen F, et al. Significance of DMBT1 in papillary thyroid carcinoma concurrent with hashimoto’s thyroiditis. Front Oncol. (2021) 11:680873. doi: 10.3389/fonc.2021.680873
26. Li B, Cai Z, Zhang Y, Chen R, Tang S, Kong F, et al. Biomarkers associated with papillary thyroid carcinoma and Hashimoto’s thyroiditis: Bioinformatic analysis and experimental validation. Int Immunopharmacol. (2024) 143:113532. doi: 10.1016/j.intimp.2024.113532
27. Zhang Y, Xie X, Zhou H, Li B, Ding L, Cai Z, et al. Identification of SERPINA1 promoting better prognosis in papillary thyroid carcinoma along with Hashimoto’s thyroiditis through WGCNA analysis. Front Endocrinol (Lausanne). (2023) 14:1131078. doi: 10.3389/fendo.2023.1131078
28. Jiang H, He Y, Lan X, and Xie X. Identification and validation of potential common biomarkers for papillary thyroid carcinoma and Hashimoto’s thyroiditis through bioinformatics analysis and machine learning. Sci Rep. (2024) 14:15578. doi: 10.1038/s41598-024-66162-2
29. Ma H, Li G, Huo D, Su Y, Jin Q, Lu Y, et al. Impact of Hashimoto’s thyroiditis on the tumor microenvironment in papillary thyroid cancer: insights from single-cell analysis. Front Endocrinol (Lausanne). (2024), 15:1339473. doi: 10.3389/fendo.2024.1339473
30. Lv T, Zhu C, and Di Z. Risk factors stratifying Malignancy of nodules in contralateral thyroid lobe in patients with pre-operative ultrasound indicated unilateral papillary thyroid carcinoma: A retrospective analysis from single centre. Clin Endocrinol (Oxf). (2018) 88:279–84. doi: 10.1111/cen.13506
31. Wang W, Kong L, Guo H, and Chen X. Prevalence and predictor for Malignancy of contralateral thyroid nodules in patients with unilateral PTMC: a systematic review and meta-analysis. Endocr Connect. (2021) 10:656–66. doi: 10.1530/EC-21-0164
Keywords: Hashimoto’s thyroiditis, papillary thyroid cancer, predictors, risk factors, thyroid cancer
Citation: Li R, Sun X, Hu L, Yu Z, Liu N, Li P, Zhao X and Liang W (2026) Hashimoto’s thyroiditis is associated with reduced invasiveness in papillary thyroid carcinoma: a propensity score-matched retrospective cohort study. Front. Immunol. 16:1745452. doi: 10.3389/fimmu.2025.1745452
Received: 13 November 2025; Accepted: 31 December 2025; Revised: 23 December 2025;
Published: 21 January 2026.
Edited by:
Wenbo Ma, Southern Medical University, ChinaReviewed by:
Martyna Borowczyk, Poznan University of Medical Sciences, PolandHuizhu Cai, Peking University People’s Hospital, China
Zhouqi Tang, Stanford University, United States
Copyright © 2026 Li, Sun, Hu, Yu, Liu, Li, Zhao and Liang. 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: Xudong Zhao, NjAxNDg5NTU0QHFxLmNvbQ==; Wenquan Liang, bHdxMDYxOEAxNjMuY29t
†These authors have contributed equally to this work
Liai Hu1,2,3†