Skip to main content

Co-occurrence of thyroid and breast cancer is associated with an increased oncogenic SNP burden



Epidemiological evidence suggests that synchronous or metachronous presentation of breast and thyroid cancers exceeds that predicted by chance alone. The following potential explanations have been hypothesized: common environmental or hormonal factors, oncogenic effect of the treatment for the first cancer, closer follow-up of cancer survivors, shared underlying genetic risk factors. While some cases were found to be related to monogenic disorders with autosomal inheritance, the genetic background of most cases of co-occurring breast and thyroid cancer is thought to be polygenic.


In this retrospective case-control study we compared the genetic profile of patients with a history of breast cancer (n = 15) to patients with co-occurring breast and thyroid cancer (n = 19) using next generation sequencing of 112 hereditary cancer risk genes. Identified variants were categorized based on their known association with breast cancer and oncogenesis in general.


No difference between patients with breast and double cancers was observed in clinical and pathological characteristics or the number of neutral SNPs. The unweighted and weighted number of SNPs with an established or potential association with breast cancer was significantly lower in the group with breast cancer only (mean difference − 0.58, BCa 95% CI [− 1.09, − 0.06], p = 0.029, and mean difference − 0.36, BCa 95% CI [− 0.70, − 0.02], p = 0.039, respectively). The difference was also significant when we compared the number of SNPs with potential or known association with any malignancy (mean difference − 1.19, BCa 95% CI [− 2.27, − 0.11], p = 0.032 for unweighted, and mean difference − 0.73, BCa 95% CI [− 1.32, − 0.14], p = 0.017 for weighted scores).


Our findings are compatible with the hypothesis of genetic predisposition in the co-occurrence of breast and thyroid cancer. Further exploration of the underlying genetic mechanisms may help in the identification of patients with an elevated risk for a second cancer at the diagnosis of the first cancer.

Peer Review reports


An increased co-occurrence of breast cancer (BC) and differentiated thyroid cancer (TC) compared to chance alone has been found repeatedly in recent decades [1, 2]. Numerous studies have shown a bidirectional association between these cancers: a more pronounced risk of thyroid cancer in breast cancer survivors and a less pronounced increase in the incidence of breast cancer after thyroid cancer [2,3,4].

A number of hypotheses have been suggested to explain these findings [3, 4]. The one, most frequently cited, describes a common (yet unknown) hormonal etiologic factor. This hypothesis is supported by the observation of a strong female predominance of both tumor types. Estrogen has been implicated in the pathogenesis of both breast and, to a lesser degree, thyroid cancer [5]. Furthermore, thyroid dysfunction has also been tentatively linked to carcinogenesis [6,7,8]. It is also possible that the hormonal milieu of pregnancy (hCG, increased TRH and prolactin) would increase the thyroid and the estrogen-like effects of thyroid hormones, hence the risk of both cancer types [9]. Endocrine disrupting chemicals affecting estrogen receptors have also been theoretically implicated [10].

Another set of hypotheses emphasize carcinogenic effects of previous cancer treatment. Radio- and hormone therapy for breast cancer and I131 treatment for thyroid tumors have both been implicated in carcinogenesis [11,12,13]. It is also possible that the increased surveillance of cancer survivors could lead to earlier recognition or even overdiagnosis of a second primary malignancy.

A shared genetic background of these cancers is also a potential explanation for the observed association. Cowden and Cowden-like syndromes are characterized by the presence of hamartomas and an increased risk of, among others, both thyroid and breast cancer. These conditions are mostly monogenic disorders with an autosomal dominant inheritance. Mutations in the PTEN, SDHB, SDHD, MTHFR and PARP4 genes, and hypermethylation of KLLN have all been implicated [14,15,16]. Much less is known about the genetic background of sporadic cases of metachronous or synchronous breast and thyroid cancer. However, an increased familial risk of other primary malignancies of affected patients points towards underlying germline mutations [3, 17].

In the present paper, we investigated the potential for a shared genetic background of thyroid and breast cancers. The recognition and elucidation of such genetic risk factors could aid in the identification, follow-up, and potential preventive treatment of high-risk individuals with a first primary malignancy.


Patients and setting

We analyzed the genetic polymorphisms in 112 hereditary cancer risk genes in a case-control study of patients with a history of either breast or both breast and thyroid cancer.

Patients were selected from a pool of n = 274 thyroid cancer patients who had received I131 treatment at our institution between January 2014 and October 2018. Twenty-one of these individuals received treatment for breast cancer before or after the diagnosis of a thyroid cancer. Two patients could not have been reached for consent, resulting in a final sample of 19. Given that the risk of thyroid cancer development after breast cancer is larger than vice versa, the control group consisted of 15 subjects who received treatment for breast cancer at our institution at least 12 years before the present study. Patients with a personal or a family history of thyroid cancer were excluded from this group. These criteria were specified to minimize contamination of the control group with patients with an increased genetic risk for thyroid cancer development. No other inclusion or exclusion criteria were set.

We collected the following clinical characteristics from hospital and outpatient records: age at cancer diagnosis, tumor size, histology, staging and grading, presence or absence of lymph node and distant metastases, details of surgical-, radiation- chemotherapeutic-, hormonal and I131 treatment. Cancer diagnosis for all cases was based on histological data using the WHO classification of breast and endocrine tumors (4th and 3/4th edition respectively). Staging was based on the TNM system (7th edition). Informed consent was obtained from all subjects before any study related procedures were performed. The study was approved by Semmelweis University Regional and Institutional Committee of Science and Research Ethics.

Genetic analysis

We have aimed to cover all the well-established hereditary cancer risk loci. The final list of the 112 investigated genes was compiled based on literature data [18,19,20], and is shown in Table 1.

Table 1 List of cancer risk genes assessed. BC associated genes are underlined

Exome amplicon library was prepared using the Ion AmpliSeq Library Kit Plus combined with the Ion AmpliSeq Exome RDY kit (ThermoFisher, MA, USA). Briefly, 100 ng of genomic DNA was added to dried down, ultra-high multiplexed primer pairs (12 pools) in a 96-well plate and amplified with the following PCR conditions: at 99 °C for 2 min; at 99 °C for 15 s and at 60 °C for 16 min (10 cycles) and holding at 10 °C. Primers were partially digested using a FuPa reagent, and then sequencing adapters and barcodes were ligated to the amplicons. The library was purified using the Agencourt AMPure XP Reagent (Beckmann Coulter, CA, USA). The concentration of the final library was determined by Ion Library TaqMan Quantitation Kit (ThermoFisher, MA, USA) on an ABI 7500 qPCR instrument with absolute quantification method.

Template preparation was performed with Ion 540 OT2 Kit (ThermoFisher, MA, USA) on semi-automated Ion OneTouch 2 instrument using emPCR method. After breaking the emulsion, the non-templated beads were removed from the solution during the semi-automated enrichment process on Ion OneTouch ES (ThermoFisher, MA, USA) machine. Later adding the sequencing primer and polymerase, the fully prepared Ion Sphere Particles (ISPs) were loaded into an Ion 540 chip, and the sequencing runs were performed using the Ion S5 Sequencing kit (ThermoFisher, MA, USA) with 500 flows.

Sequence data from the Ion Torrent run were analyzed using the platform-specific pipeline software Torrent Suite v5.10 for base calling, trim adapter and primer sequences, filtering out poor quality reads, and de-multiplex the reads according to the barcode sequences. Briefly, TMAP algorithm was used to align the reads to the hg19 human reference genome, and then, the variant caller plug-in was executed to search for germline variants in the targeted regions. Integrative genomics viewer (IGV) was used for visualization of the mapped reads. Variants were annotated using the Ion Reporter (ThermoFisher, MA, USA) software. Variant interpretation was restricted to the selected genes with known connection to oncology diseases. Variants were stratified by gene association with breast cancer and by variant association with cancer (no/potential/known) based on available data from the ClinVar database at the time of writing [21]. Synonymous mutations were not excluded from the analysis.

Statistical analysis

Descriptive data are given as means ± standard deviation (SD) for continuous data, and frequencies (percentages) for categorical variables. Between-group differences in clinical parameters were assessed using 2-sample t tests for continuous variables, and χ2 tests or Fisher exact tests for categorical variables.

When comparing the number of different variants across groups, the assumptions of normality and homogeneity of variances were assessed using the Shapiro–Wilk test and Levene’s test respectively. Given that in some cases these assumptions were not met, Welch’s t-test was used to compare the number of variants with known, potential or no association with breast cancer and with any cancer. We also used the bootstrap procedure to get robust confidence intervals and p-values. To increase statistical power, we lumped together variants with potential and known associations using both weighted and unweighted scores. To generate the weighted scores, we gave double weight for variants with known association compared to those variants with potential association. To correct for multiple testing we applied the Benjamini-Hochberg procedure targeting a false discovery rate < 15%. A corrected two-tailed p-value less than 0.05 was considered statistically significant. Effect sizes were computed as adjusted standardized mean differences (Hedge’s g). Analyses were conducted using IBM SPSS Statistics for Windows, Version 23.0 (IBM Corp. Released 2015 Armonk, New York, U.S.A.).


Table 2 displays patients’ clinical characteristics. BC only patients had a slightly younger age at the time of breast cancer diagnosis compared to the synchronous/metachronous cancer group (47.7 vs 54.4, p = 0.079). All other clinical parameters evaluated were similar in the BC and BC-TC groups. It is worth noting that no difference was found in the frequency of the different BC treatment modalities (irradiation, hormone or chemotherapy).

Table 2 Patient characteristics at BC diagnosis and at follow-up

The results of the genetic analyses for the BC-TC patients and the control group are presented in Tables 3 and 4 respectively. Details of the statistical analysis are reported in Table 5.

Table 3 List of any cancer (and breast cancer) related genes and unique variants identified by individual patients in the BC-TC group. The strength of known association of each variant with oncogenesis is marked as known (**), potential (*) or no (no marker)
Table 4 List of any cancer (and breast cancer) related genes and unique variants identified by individual patients in the control group. The strength of known association of each variant with oncogenesis is marked as known (**), potential (*) or no (no marker)
Table 5 Comparison of the two patient groups by the number of genetic variants

We found a significantly lower number of SNPs with potential or known association with breast cancer in the BC only group (mean difference − 0.58, BCa 95% CI [− 1.09, − 0.06], t (27) = − 2.30, p = 0.029, g = 0.74 for unweighted and mean difference − 0.36, BCa 95% CI [− 0.70, − 0.02], t(27.26) = − 2.17, p = 0.039, g = 0.70 for weighted scores).

The individual number of SNPs with no known genetic association with cancer was nominally higher in the BC only group (mean difference 0.46, BCa 95% CI [− 0.27, 1.20], t(25.56) = 1.30, p = 0.206, g = 0.39), while the number of potentially associated SNP or those with a known association was nominally lower in the BC only group (mean difference − 0.92, BCa 95% CI [− 2.01, 0.17], t(30.98) = − 1.72, p = 0.096, g = 0.59 and mean difference − 0.27, BCa 95% CI [− 0.63, 0.09], t(31.37) = − 1.55, p = 0.131, g = 0.51, respectively).

When we lumped together SNPs with potential or known association with carcinogenesis to increase statistical power, we found significantly higher values in the synchronous/metachronous cancer group compared to the BC only group using both unweighted (mean difference − 1.19, BCa 95% CI [− 2.27, − 0.11], t(31.88) = − 2.24, p = 0.032, g = 0.76) and weighted (mean difference − 0.73, BCa 95% CI [− 1.32, − 0.14], t(31.80) = − 2.51, p = 0.017, g = 0.83) scores.


In our study, we found a significantly higher burden of established and potential germline cancer risk variants among subjects with a history of metachronous thyroid and breast cancer compared to patients with a history of breast cancer only. The trend was similar and effect sizes were still considerable when we assessed the number of SNPs with known and potential cancer risk separately. However, in these cases, probably due to lack of power, statistical significance was not reached. Our results are compatible with the hypothesis that the sporadic co-occurrence of thyroid and breast cancer is of multigenetic origin and probably related to the burden of the carcinogenic SNPs rather than an individual gene alteration.

Comparison of clinical features of breast cancer, such as onset, stage, histology and treatment also yielded no significant differences. There was a non-significant tendency for later age of breast cancer onset in the study group which is contrary to previous findings [22].

The metachronous and synchronous occurrence of breast and thyroid cancer has been first described nearly four decades ago [1]. The relationship between the two tumor types has since been found to be bidirectional. In recent meta-analyses, the odds ratio for thyroid cancer following breast cancer treatment was 1.55 while the odds ratio of developing BC after TC was reported between 1.18–1.32 [3, 4]. Studies exploring the cause behind this association emphasize either the role of common underlying hormonal and environmental factors, or the role of the first malignancy in the development of the second. The potential role of cancer treatment and the increased surveillance of cancer survivors fall in this latter category. With breast cancer being the most common type of malignancy in women [23] and with the increasing incidence of thyroid cancer worldwide [24, 25], the importance of this topic is clear.

One potential hypothesis explaining this bidirectional relationship implicates common hormonal factors and is rooted in the overwhelming predominance of both cancer types in reproductive-aged women [4]. Estrogen, progesterone and androgen receptors have been shown to be overexpressed not only in breast cancer but also in thyroid neoplasms [26]. In vitro and animal studies also point towards the potential oncogenic effect of estrogens in thyroid cells [26, 27]. An increased TSH secretion in response to estrogens has also been suggested as a potential pathomechanism for thyroid cancer development [3], as serum TSH levels have been shown to correlate with thyroid cancer risk [28]. TSH and thyroid hormones also have been implicated in breast carcinogenesis by in vitro animal and observational studies. T3 seems to stimulate proliferation in breast cancer cells in vitro at least in part through interactions with the estrogen signalling system [29, 30]. Tumor suppressor pathways downstream of the TRβ nuclear thyroid receptor and oncogenic pathways downstream of the membrane receptor αVβ3 have been suggested to mediate the role of thyroid hormones in carcinogenesis [31,32,33]. While some observational studies in humans found a decreased breast cancer incidence in hypothyroid and an increased incidence in hyperthyroid patients [7, 8], other studies have failed to replicate these findings [34, 35].

Endocrine disrupting chemicals (EDC) have been linked to numerous types of hormone dependent malignancies including both thyroid and breast cancer [36, 37]. In addition to EDCs, obesity, a shared risk factor for both tumor types has also been postulated to increase the risk of synchronous/metachronous cancer development [38, 39]. Selective estrogen receptor modulators (SERMs) are widely used in breast cancer treatment. While there is some evidence that these drugs due to their partial estrogen effect increase TSH secretion and thyroid cell proliferation, their potential role in thyroid cancer development has not yet been studied [13]. The similar frequency of hormone therapy in the BC only, and the BC-TC groups in our study point toward a limited role of sex hormones in the development of co-occurring breast and thyroid cancer.

Other theories suggest the role of prior cancer treatment in the development of the second primary malignancy. Despite several improvements minimizing radiation scatter to surrounding tissues, there is a definitely increased risk for certain malignancies following external beam radiation for breast cancer [40]. While some reports are conflicting [41], most available data do not substantiate such a connection between thyroid cancer and previous adjuvant breast irradiation [42,43,44,45]. Similarly, radioactive iodine treatment given for thyroid cancer does not seem to play a role in subsequent breast cancer development [46,47,48]. Our data with similar frequency of radiation therapy for breast cancer in the investigated groups also argue against the role of external or internal radiation in the development of synchronous/metachronous breast and thyroid cancer.

The increased co-occurrence of these tumor types could also be related to surveillance bias. The higher compliance with, and higher rate of screening efforts among cancer survivors could also affect the time of diagnosis and could lead the substantial overdiagnosis of clinically irrelevant second primary malignancy [49]. This seems to be especially true for differentiated thyroid cancer [50]. However, as there are no screening programs for breast cancer among thyroid cancer survivors or vice versa, surveillance bias is unlikely to account for the preferential association of these tumor types.

Cowden syndrome (CS) and Cowden-like syndromes (CLS) are characterised by hamartomas and an extremely increased risk for several types of malignancies including breast cancer, thyroid cancer, endometrial cancer, colorectal cancer, and melanoma with standardized incidence ratios for breast and thyroid cancer in the range of 6 to 9 [14]. CS and most forms of CLS have an autosomal dominant inheritance. However, these conditions, referred to as PTEN hamartoma tumor syndromes, make up only a small fraction of synchronous/metachronous thyroid and breast cancer cases. This is consistent with our finding that no patient in our samples had a mutation in the PTEN gene. Thus, the existence of other shared genetic risk factors is highly probable. To the best of our knowledge, genetic analyses exploring this notion have not been conducted and genetic factors underlying sporadic metachronous cancer cases are yet to be elucidated.

The main strengths of this study were the state-of-the-art genetic analysis using NGS technology, and the fact that both cases and controls came from the same source population. The main limitation of our research was the relatively low number of patients limiting statistical power even in analyses of genetic scores. Furthermore, our study was unable to identify or investigate any single genetic risk locus behind synchronous/metachronous cancer development. Two types of bias could also have affected our findings. First, participants were recruited years after diagnosis potentially leading to survivor bias. Second, the exclusion of controls with less than 12 years of follow-up lead to an age difference between our study groups that is contrary to previous findings [22], and probably reflects selection bias. However, both of these selections were deemed to be necessary to minimize contamination of the control group with thyroid cancer.


In conclusion, we reported an increased burden of carcinogenic SNPs in people with both thyroid and breast cancer compared to individuals with breast cancer only, based on whole exome sequencing of 112 known hereditary cancer risk genes. We found no differences in clinicopathologic parameters between the groups suggesting that both groups have similar presentation at the time of diagnosis of breast cancer. While our study was not powered to identify specific risk loci for metachronous cancer development, our findings further support the multigenetic etiology of co-occurring breast and thyroid cancer. Our findings do not directly contradict any of the other, previously detailed theories explaining the association between these tumor types. Nevertheless, our results do underline the need for further genetic research in this field, which are lacking at the moment.

Availability of data and materials

All data generated or analysed during this study are included in this published article.



Breast cancer

BCa 95% CI:

Bias-corrected and accelerated 95% bootstrap confidence interval


Cowden-like syndromes


Cowden syndrome


Endocrine disrupting chemical


Human chorionic gonadotropin


Integrative genomics viewer


Ion sphere particles


Standard deviation


Selective estrogen receptor modulator


Thyroid cancer


Thyrotropin-releasing hormone


  1. McTiernan A, Weiss NS, Daling JR. Incidence of thyroid cancer in women in relation to known or suspected risk factors for breast cancer. Cancer Res. 1987;47(1):292–5.

    CAS  PubMed  Google Scholar 

  2. An JH, Hwangbo Y, Ahn HY, Keam B, Lee KE, Han W, et al. A possible association between thyroid Cancer and breast Cancer. Thyroid. 2015;25(12):1330–8.

    Article  CAS  PubMed  Google Scholar 

  3. Nielsen SM, White MG, Hong S, Aschebrook-Kilfoy B, Kaplan EL, Angelos P, et al. The breast-thyroid Cancer link: a systematic review and meta-analysis. Cancer Epidemiol Biomark Prev. 2016;25(2):231–8.

    Article  CAS  Google Scholar 

  4. Joseph KR, Edirimanne S, Eslick GD. The association between breast cancer and thyroid cancer: a meta-analysis. Breast Cancer Res Treat. 2015;152(1):173–81.

    Article  PubMed  Google Scholar 

  5. Moleti M, Sturniolo G, Di Mauro M, Russo M, Vermiglio F. Female Reproductive Factors and Differentiated Thyroid Cancer. Front Endocrinol (Lausanne). 2017;8:111.

    Article  Google Scholar 

  6. Hardefeldt PJ, Eslick GD, Edirimanne S. Benign thyroid disease is associated with breast cancer: a meta-analysis. Breast Cancer Res Treat. 2012;133(3):1169–77.

    Article  CAS  PubMed  Google Scholar 

  7. Cristofanilli M, Yamamura Y, Kau SW, Bevers T, Strom S, Patangan M, et al. Thyroid hormone and breast carcinoma. Primary hypothyroidism is associated with a reduced incidence of primary breast carcinoma. Cancer. 2005;103(6):1122–8.

    Article  CAS  PubMed  Google Scholar 

  8. Shi XZ, Jin X, Xu P, Shen HM. Relationship between breast cancer and levels of serum thyroid hormones and antibodies: a meta-analysis. Asian Pac J Cancer Prev. 2014;15(16):6643–7.

    Article  PubMed  Google Scholar 

  9. Kurata A. Differentiated thyroid cancer: why does it affect predominantly women during the reproductive period and have higher incidence of mutual association with breast cancer? Med Hypotheses. 2019;122:5–7.

    Article  PubMed  Google Scholar 

  10. Hoffman K, Lorenzo A, Butt CM, Hammel SC, Henderson BB, Roman SA, et al. Exposure to flame retardant chemicals and occurrence and severity of papillary thyroid cancer: a case-control study. Environ Int. 2017;107:235–42.

    Article  CAS  PubMed  Google Scholar 

  11. Tunio MA, Al Asiri M, Bayoumi Y, Stanciu LG, Al Johani N, Al Saeed EF. Is thyroid gland an organ at risk in breast cancer patients treated with locoregional radiotherapy? Results of a pilot study. J Cancer Res Ther. 2015;11(4):684–9.

    Article  CAS  PubMed  Google Scholar 

  12. Zhang Y, Liang J, Li H, Cong H, Lin Y. Risk of second primary breast cancer after radioactive iodine treatment in thyroid cancer: a systematic review and meta-analysis. Nucl Med Commun. 2016;37(2):110–5.

    Article  CAS  PubMed  Google Scholar 

  13. Zidan J, Rubenstein W. Effect of adjuvant tamoxifen therapy on thyroid function in postmenopausal women with breast cancer. Oncology. 1999;56(1):43–5.

    Article  CAS  PubMed  Google Scholar 

  14. Ngeow J, Stanuch K, Mester JL, Barnholtz-Sloan JS, Eng C. Second malignant neoplasms in patients with Cowden syndrome with underlying germline PTEN mutations. J Clin Oncol. 2014;32(17):1818–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Ni Y, He X, Chen J, Moline J, Mester J, Orloff MS, et al. Germline SDHx variants modify breast and thyroid cancer risks in Cowden and Cowden-like syndrome via FAD/NAD-dependant destabilization of p53. Hum Mol Genet. 2012;21(2):300–10.

    Article  CAS  PubMed  Google Scholar 

  16. Bolf EL, Sprague BL, Carr FE. A linkage between thyroid and breast Cancer: a common etiology? Cancer Epidemiol Biomark Prev. 2019;28(4):643–9.

    Article  Google Scholar 

  17. Goldgar DE, Easton DF, Cannon-Albright LA, Skolnick MH. Systematic population-based assessment of cancer risk in first-degree relatives of cancer probands. J Natl Cancer Inst. 1994;86(21):1600–8.

    Article  CAS  PubMed  Google Scholar 

  18. Bailey MH, Tokheim C, Porta-Pardo E, Sengupta S, Bertrand D, Weerasinghe A, et al. Comprehensive characterization of Cancer driver genes and mutations. Cell. 2018;174(4):1034–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Foulkes WD. Inherited susceptibility to common cancers. N Engl J Med. 2008;359(20):2143–53.

    Article  CAS  PubMed  Google Scholar 

  20. Shah PD, Nathanson KL. Application of panel-based tests for inherited risk of Cancer. Annu Rev Genomics Hum Genet. 2017;18(1):201–27.

    Article  CAS  PubMed  Google Scholar 

  21. Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, et al. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016;44(D1):D862–8.

    Article  CAS  PubMed  Google Scholar 

  22. Kuo JH, Chabot JA, Lee JA. Breast cancer in thyroid cancer survivors: An analysis of the surveillance, epidemiology, and end Results-9 database. Surgery. 2016;159(1):23–9.

    Article  PubMed  Google Scholar 

  23. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  25. Davies L, Welch HG. Current thyroid cancer trends in the United States. JAMA Otolaryngol Head Neck Surg. 2014;140(4):317–22.

    Article  PubMed  Google Scholar 

  26. Yane K, Kitahori Y, Konishi N, Okaichi K, Ohnishi T, Miyahara H, et al. Expression of the estrogen receptor in human thyroid neoplasms. Cancer Lett. 1994;84(1):59–66.

    Article  CAS  PubMed  Google Scholar 

  27. Rajoria S, Suriano R, Shanmugam A, Wilson YL, Schantz SP, Geliebter J, et al. Metastatic phenotype is regulated by estrogen in thyroid cells. Thyroid. 2010;20(1):33–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Haymart MR, Repplinger DJ, Leverson GE, Elson DF, Sippel RS, Jaume JC, et al. Higher serum thyroid stimulating hormone level in thyroid nodule patients is associated with greater risks of differentiated thyroid cancer and advanced tumor stage. J Clin Endocrinol Metab. 2008;93(3):809–14.

    Article  CAS  PubMed  Google Scholar 

  29. Hall LC, Salazar EP, Kane SR, Liu N. Effects of thyroid hormones on human breast cancer cell proliferation. J Steroid Biochem Mol Biol. 2008;109(1–2):57–66.

    Article  CAS  PubMed  Google Scholar 

  30. Silva JM, Dominguez G, Gonzalez-Sancho JM, Garcia JM, Silva J, Garcia-Andrade C, et al. Expression of thyroid hormone receptor/erbA genes is altered in human breast cancer. Oncogene. 2002;21(27):4307–16.

    Article  CAS  PubMed  Google Scholar 

  31. Kim WG, Cheng SY. Thyroid hormone receptors and cancer. Biochim Biophys Acta. 2013;1830(7):3928–36.

    Article  CAS  PubMed  Google Scholar 

  32. Martinez-Iglesias O, Garcia-Silva S, Tenbaum SP, Regadera J, Larcher F, Paramio JM, et al. Thyroid hormone receptor beta1 acts as a potent suppressor of tumor invasiveness and metastasis. Cancer Res. 2009;69(2):501–9.

    Article  CAS  PubMed  Google Scholar 

  33. Carr FE, Tai PW, Barnum MS, Gillis NE, Evans KG, Taber TH, et al. Thyroid hormone receptor-beta (TRbeta) mediates runt-related transcription factor 2 (Runx2) expression in thyroid Cancer cells: a novel signaling pathway in thyroid Cancer. Endocrinology. 2016;157(8):3278–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Muller I, Kilburn LS, Taylor PN, Barrett-Lee PJ, Bliss JM, Ellis P, et al. TPOAb and thyroid function are not associated with breast Cancer outcome: evidence from a large-scale study using data from the Taxotere as adjuvant chemotherapy trial (TACT, CRUK01/001). Eur Thyroid J. 2017;6(4):197–207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Fang Y, Yao L, Sun J, Yang R, Chen Y, Tian J, et al. Does thyroid dysfunction increase the risk of breast cancer? A systematic review and meta-analysis. J Endocrinol Investig. 2017;40(10):1035–47.

    Article  CAS  Google Scholar 

  36. Gore AC, Chappell VA, Fenton SE, Flaws JA, Nadal A, Prins GS, et al. EDC-2: the Endocrine Society's second scientific statement on endocrine-disrupting chemicals. Endocr Rev. 2015;36(6):E1–E150.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Soto AM, Vandenberg LN, Maffini MV, Sonnenschein C. Does breast cancer start in the womb? Basic Clin Pharmacol Toxicol. 2008;102(2):125–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Harvie M, Howell A, Evans DG. Can diet and lifestyle prevent breast cancer: what is the evidence? Am Soc Clin Oncol Educ Book. 2015;(35):e66–73.

  39. Ma J, Huang M, Wang L, Ye W, Tong Y, Wang H. Obesity and risk of thyroid cancer: evidence from a meta-analysis of 21 observational studies. Med Sci Monit. 2015;21:283–91.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Matesich SM, Shapiro CL. Second cancers after breast cancer treatment. Semin Oncol. 2003;30(6):740–8.

    Article  CAS  PubMed  Google Scholar 

  41. Grantzau T, Overgaard J. Risk of second non-breast cancer among patients treated with and without postoperative radiotherapy for primary breast cancer: a systematic review and meta-analysis of population-based studies including 522,739 patients. Radiother Oncol. 2016;121(3):402–13.

    Article  PubMed  Google Scholar 

  42. Li CI, Rossing MA, Voigt LF, Daling JR. Multiple primary breast and thyroid cancers: role of age at diagnosis and cancer treatments (United States). Cancer Causes Control. 2000;11(9):805–11.

    Article  CAS  PubMed  Google Scholar 

  43. Huang J, Walker R, Groome PG, Shelley W, Mackillop WJ. Risk of thyroid carcinoma in a female population after radiotherapy for breast carcinoma. Cancer. 2001;92(6):1411–8.<1411::AID-CNCR1464>3.0.CO;2-9.

    Article  CAS  PubMed  Google Scholar 

  44. Lal G, Groff M, Howe JR, Weigel RJ, Sugg SL, Lynch CF. Risk of subsequent primary thyroid cancer after another malignancy: latency trends in a population-based study. Ann Surg Oncol. 2012;19(6):1887–96.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Kim SS, Kim SJ, Bae YT, Lee JY, Kim BH, Kim YK, et al. Factors associated with the development of new onset diffuse thyroid F18-fluorodeoxyglucose uptake after treatment of breast cancer in patients without a history of thyroid disease or thyroid dysfunction. Thyroid. 2012;22(1):53–8.

    Article  CAS  PubMed  Google Scholar 

  46. Brown AP, Chen J, Hitchcock YJ, Szabo A, Shrieve DC, Tward JD. The risk of second primary malignancies up to three decades after the treatment of differentiated thyroid cancer. J Clin Endocrinol Metab. 2008;93(2):504–15.

    Article  CAS  PubMed  Google Scholar 

  47. Sawka AM, Thabane L, Parlea L, Ibrahim-Zada I, Tsang RW, Brierley JD, et al. Second primary malignancy risk after radioactive iodine treatment for thyroid cancer: a systematic review and meta-analysis. Thyroid. 2009;19(5):451–7.

    Article  CAS  PubMed  Google Scholar 

  48. Verkooijen RB, Smit JW, Romijn JA, Stokkel MP. The incidence of second primary tumors in thyroid cancer patients is increased, but not related to treatment of thyroid cancer. Eur J Endocrinol. 2006;155(6):801–6.

    Article  CAS  PubMed  Google Scholar 

  49. Ito Y, Miyauchi A, Oda H. Low-risk papillary microcarcinoma of the thyroid: a review of active surveillance trials. Eur J Surg Oncol. 2018;44(3):307–15.

    Article  CAS  PubMed  Google Scholar 

  50. Ahn HS, Kim HJ, Welch HG. Korea's thyroid-cancer "epidemic"--screening and overdiagnosis. N Engl J Med. 2014;371(19):1765–7.

    Article  PubMed  Google Scholar 

Download references


Not applicable.


The research was financed by the FIKP and the Thematic Excellence Program (Tématerületi Kiválósági Program, 2020–4.1.1.-TKP2020) of the Ministry for Innovation and Technology in Hungary, within the framework of the Molecular Biology thematic program of Semmelweis University. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



BB1, AK, PL, and IT conceived of the presented idea. All authors were involved in planning and supervising the work. BB1, BB2, AK, PL, MD, ZsV supervised sample and data collection. Genetic data collection, analysis and interpretation were carried out by KÁ, BB2, BT, and JK. Statistical analysis was done by ÁT, BSz, BDK, and RÁ. Manuscript was drafted by BB1, AK, PL, ZsP, and RA, and was substantively revised by MD and ZsV. All authors discussed the results and contributed to the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bence Bakos.

Ethics declarations

Ethics approval and consent to participate

This research was approved and licenced by the Semmelweis University Regional and Institutional Committee of Science and Research Ethics, and was conducted in accordance with the World Medical Association’s Declaration of Helsinki. Written informed consent was obtained from all participants before entering the study and prior to any study related procedures. This consent extended to the study team accessing their relevant medical history and genetic data. All participants were 18 years or older at the time of the study.

Consent for publication

Not applicable.

Competing interests

The authors state that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bakos, B., Kiss, A., Árvai, K. et al. Co-occurrence of thyroid and breast cancer is associated with an increased oncogenic SNP burden. BMC Cancer 21, 706 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Thyroid cancer
  • Breast cancer
  • Metachronous cancer
  • Oncogenesis