eTable 1. Standardized Mortality Ratios (SMRs) and 95% Confidence Intervals (CIs) for Solid Cancer Mortality by Treatment Combination, Based on Ratio of Observed Numbers of Deaths Among Patients With Hyperthyroidism Versus Expected Numbers of Deaths in the General Population, by Time Between Study Entry and Death
eTable 2. Standardized Mortality Ratios (SMRs) by Treatment Combination, Based on Ratio of Observed Numbers of Deaths Among Patients With Hyperthyroidism Versus Expected Numbers of Deaths in the General Population, Excluding Patients With Prior Cancers
Customize your JAMA Network experience by selecting one or more topics from the list below.
Identify all potential conflicts of interest that might be relevant to your comment.
Conflicts of interest comprise financial interests, activities, and relationships within the past 3 years including but not limited to employment, affiliation, grants or funding, consultancies, honoraria or payment, speaker's bureaus, stock ownership or options, expert testimony, royalties, donation of medical equipment, or patents planned, pending, or issued.
Err on the side of full disclosure.
If you have no conflicts of interest, check "No potential conflicts of interest" in the box below. The information will be posted with your response.
Not all submitted comments are published. Please see our commenting policy for details.
Kitahara CM, Preston DL, Sosa JA, Berrington de Gonzalez A. Association of Radioactive Iodine, Antithyroid Drug, and Surgical Treatments With Solid Cancer Mortality in Patients With Hyperthyroidism. JAMA Netw Open. 2020;3(7):e209660. doi:10.1001/jamanetworkopen.2020.9660
这项队列研究涉及 31,363 名在基线时没有癌症的患者。研究在治疗组（放射性碘、药物和/或手术）与实体癌死亡风险之间未发现关联。在接受放射性碘治疗的患者中，与实体癌死亡率的关联性随着总给药活动的增加而增加。
The long-term health effects of radioactive iodine (RAI) and antithyroid drug (ATD) treatments compared with surgery for hyperthyroidism remain uncertain.
To compare solid cancer mortality rates associated with RAI and ATD treatments vs surgical management for hyperthyroidism.
Design, Setting, and Participants
This multicenter cohort study assessed patients treated for hyperthyroidism from January 1, 1946, to December 31, 1964, with follow-up through December 31, 2014. Data analysis was performed from August 1, 2019, to April 23, 2020.
Management with RAI, ATDs, surgical intervention, or combinations of these treatments.
Main Outcomes and Measures
Comparisons of solid cancer mortality rates in each treatment group with expected rates from the general population were assessed using standardized mortality ratios (SMRs), and internal comparisons were assessed using hazard ratios (HRs) adjusted for age, sex, and underlying diagnosis (Graves disease or toxic nodular goiter).
Of 31 363 patients (24 894 [79.4%] female; mean [SD] age, 46.9 [14.8] years) included in the study, 28 523 (90.9%) had Graves disease. The median follow-up time was 26.0 years (interquartile range, 12.3-41.9 years). Important differences in patient characteristics existed across treatment groups at study entry. Notably, the drug-only group (3.6% of the cohort) included a higher proportion of patients with prior cancers (7.3% vs 1.9%-4.0%), contributing to an elevated SMR for solid cancer mortality. After excluding prior cancers, solid cancer SMRs were not elevated in any of the treatment groups (SMR for surgery only, 0.82 [95% CI, 0.66-1.00]; SMR for drugs only, 0.90 [95% CI, 0.74-1.09]; SMR for drugs and surgery, 0.88 [95% CI, 0.84-0.94]; SMR for RAI only, 0.90 [95% CI, 0.84-0.96]; SMR for surgery and RAI, 0.66 [95% CI, 0.52-0.85]; SMR for drugs and RAI, 0.94 [95% CI, 0.89-1.00]; and SMR for drugs, surgery, and RAI, 0.85 [95% CI, 0.75-0.96]), and no significant HRs for solid cancer death were observed across treatment groups. Among RAI-treated patients, HRs for solid cancer mortality increased significantly across levels of total administered activity (1.08 per 370 MBq; 95% CI, 1.03-1.13 per 370 MBq); this association was stronger among patients treated with only RAI (HR, 1.19 per 370 MBq; 95% CI, 1.09-1.30 per 370 MBq).
Conclusions and Relevance
After controlling for known sources of confounding, the study found no significant differences in the risk of solid cancer mortality by treatment group. However, among RAI-treated patients, a modest positive association was observed between total administered activity and solid cancer mortality, providing further evidence in support of a dose-dependent association between RAI and solid cancer mortality.
Overt hyperthyroidism is diagnosed in approximately 0.5% of the US population.1 The most common diagnosis is Graves disease, with a lifetime prevalence of approximately 3% for women and less than 1% for men. Radioactive iodine (RAI), antithyroid drugs (ATDs), and thyroidectomy alone or in combination are the 3 main modalities used in the management of hyperthyroidism; however, treatment practices vary substantially over time and by geographic region. For toxic nodular goiter (TNG), RAI and thyroidectomy are the primary treatment options, with long-term ATD treatment offered less frequently.1 For uncomplicated Graves disease, RAI has been the treatment of choice in the US, whereas ATDs have been the more commonly preferred option in most other regions of the world.2-4 However, in the US and Europe, RAI treatment of Graves disease has been decreasing in recent decades in favor of ATDs.2 Evidence of long-term risks of leukemia and other secondary malignant tumors after postsurgical adjuvant treatment of differentiated thyroid cancer with RAI may have contributed to this trend.5-9
Studies evaluating the long-term adverse effects of RAI and other treatment options for hyperthyroidism may help to better inform clinical treatment practices and guidelines. Because of the lack of large randomized clinical trials, information on long-term associations of treatment with cancer and other outcomes has, thus far, relied on results from observational studies. However, results from the few large observational studies10-16 on this topic have been inconsistent, and summarizing and comparing the findings across studies are challenging because of differences in sample size, duration of follow-up, information on treatment types and doses, statistical methods, and ability to control for confounding and other sources of bias. Confounding is especially of concern for studies that rely on an external comparison group (eg, general population rates) because patients can differ in many important ways from the general public in terms of genetic, environmental, and behavioral factors and access to care. In addition, the underlying disease itself may be directly associated with risk of the outcome under study independently of any treatment effects.17,18 Another source of bias is confounding by indication (or treatment allocation bias) in which the reason or choice of one treatment vs another is associated with the baseline risk of the outcome under study (eg, patients’ age, sex, general health status, and disease origin and severity).19,20 Thus, although some of these studies have found increased incidence or mortality from cancers of the stomach,11-13 kidney,12 breast,12 and thyroid10,14 among RAI-treated patients, the direct effects of treatment are difficult to discern in the presence of uncontrolled confounding.
With use of data from the Cooperative Thyrotoxicosis Therapy Follow-up Study (CTTFUS), a large cohort study of patients treated for hyperthyroidism between the 1940s and 1960s in the US and UK, a positive association between organ-absorbed dose and mortality from solid cancers (including breast cancer) was previously reported.21 The focus on RAI-treated patients was a strength because it helped to minimize confounding caused by imbalances in risk factors across treatment groups and provided a direct assessment of the risk per unit dose.22 Nonetheless, comparisons of risks across treatment groups are of interest for practitioners involved in the management of hyperthyroidism who must select the treatment modality and the extent of treatment.
Thus, for the current study, we used data from the CTTFUS cohort to evaluate the association between treatment received (RAI, ATDs, surgery, and various combinations of these) and risk of solid cancer death. We also investigated these risks across levels of RAI-administered activity (as a more intuitive measure and proxy for organ-absorbed doses). Results presented in the current article informed the methodologic approach used in a recent analysis21 of organ-absorbed dose. With an additional 24 years of follow-up, the current study also provides an update of the standardized mortality ratios (SMRs) for solid cancer published by Ron et al10 more than 2 decades ago.
This cohort study included data from the original CTTFUS cohort, which included 35 593 patients with hyperthyroidism treated with RAI, ATDS, surgery, or some combination of these at 25 US and 1 British medical centers between January 1, 1946, and December 31, 1964.10 Comprehensive clinical data were abstracted from medical records, and early follow-up information was obtained from treating physicians, clinics, or patients when no medical sources were available. In 1984, the National Cancer Institute and 4 regional study centers reassembled patient and treatment data from computer lists, microfiche, microfilm cassettes, and handwritten material; all but 4010 patients (n = 31 583) were eligible for subsequent tracing and mortality follow-up through 1990.10 Follow-up of US patients was recently extended through December 31, 2014.21 Person-time for all patients was defined as the date of the first study examination through the date of death, date of loss to follow-up or last known vital status, or end of study, whichever occurred first. This study was approved by the institutional review board of the National Cancer Institute of the National Institutes of Health. Because mortality follow-up was based on linkages with available databases and involved no direct contact with study participants, the requirement for informed consent was waived. All data were deidentified. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline.
For the current analyses, we additionally excluded patients with no follow-up time (n = 206), those with missing entry or exit dates (n = 8), and those whose exit dates occurred on or before the study entry date (n = 6) (Figure). Of the 31 363 remaining patients, death records were available for 25 445 (81.1%), whereas 2829 (9.0%) were censored at their last known vital status date and the others were found or presumed to be alive by the end of follow-up on December 31, 2014 (3089 [9.8%]). Data analysis was performed from August 1, 2019, to April 23, 2020.
For external comparisons, SMRs and 95% CIs were used to compare the observed number of solid cancer (including breast and nonbreast solid cancer) deaths in the cohort with the expected number of deaths based on age-, sex-, and calendar year–specific mortality rates in the general populations of the US and UK. The SMRs were calculated separately by treatment group. Influence of possible sources of bias were examined through exclusion of specific patient groups (eg, those with comorbidities) or person-time (to account for potential exposure-disease latency periods).
For internal cohort comparisons, Cox proportional hazards regression models were used to estimate multivariable-adjusted hazard ratios (HRs) and 95% CIs for the association between treatment groups (or categories of total administered activity) and risk of solid cancer death. All models excluded patients with prior cancer at study entry (n = 1013) or between entry and date of last RAI treatment (n = 36), used attained age as the time scale (incorporating a 5-year latency period), and were adjusted for sex and underlying diagnosis (Graves disease, TNG, and other or unknown). Statistical analyses were conducted with Stata, version 16.0 (StataCorp) and Epicure, version 2.00.02 (Risk Sciences International). A 2-sided P < .05 was considered to be statistically significant.
A total of 31 363 patients (24 894 [79.4%] female; mean [SD] age, 46.9 [14.8] years) were included in the study. Baseline characteristics of CTTFUS patients (1946-1964) are described in Table 1. A total of 28 523 patients (90.9%) received a diagnosis of Graves disease. Use of ATDs (alone or in combination with other treatments) was recorded for 22 357 patients (71.3%), RAI for 19 589 patients (62.5%), and surgery for 13 676 patients (43.6%). Whereas 7474 (23.8%) were treated with only RAI, 1138 (3.6%) were treated using only drugs and 800 (2.6%) underwent surgery. Data on ATD types, dosages, and treatment dates were unavailable or largely incomplete. For 2718 patients (89.0%) treated with both surgery and RAI, the date of surgery preceded the date of first RAI treatment, suggesting that these were mostly patients undergoing subtotal thyroidectomy who subsequently experienced relapse. Among RAI-treated patients, 12 979 (66.3%) received 1 round of treatment, whereas 6610 (33.7%) received multiple rounds. Median administered activity was 270.8 MBq (interquartile range [IQR], 187.8-422.5 MBq) for Graves disease, 488.1 MBq (IQR, 300.1-817.7 MBq) for TNG, 373.4 MBq (IQR, 233.3-569.9 MBq) for other or unknown, and 296.2 MBq (IQR, 188.7-446.2 MBq) for RAI-treated patients overall. Treatment groups differed by age, underlying diagnosis, and baseline prevalence of comorbidities. A total of 4990 patients (50.8%) in the drugs and surgery group were younger than 40 years compared with 10 487 (33.4%) overall. The surgery-only group included a high proportion of patients with TNG (356 [44.5%] vs 2491 [7.9%] overall), which accounted for the higher mean (SD) age at entry (50.8 [13.9] vs 46.9 [14.8] years). The drug-only group had the highest proportion with comorbidities, including cancer (83 [7.3%] vs 1015 [3.2%] overall).
Median follow-up time was 26.0 years (IQR, 12.3-41.9 years). Duration of follow-up and proportion of early censoring differed across treatment groups (Table 2), with the shortest follow-up occurring in the drug-only group (median, 20.2 years; IQR, 5.4-41.3 years) and the longest in the drugs and surgery group (median, 33.6 years; IQR, 18.4-50.5 years). The drug-only group had the highest proportion of early censoring (135 [11.9%] vs 27 [3.4%] to 916 [10.1%]). The drug-only and drugs and surgery groups also had the highest proportion of cancer-related deaths (160 [18.3%] to 1413 [18.6%] vs 78 [13.0%] to 1253 [16.7%]). However, 59 cancer deaths (36.9%) in the drug-only group occurred in the first 5 years of follow-up compared with 74 (5.2%) in the drugs and surgery group and 10 (12.8%) to 193 (15.4%) in the other treatment groups. Overall, 59 (5.2%) in the drug-only group died of cancer in the first 5 years of follow-up vs 74 (0.8%) to 193 (2.1%) in the other groups.
The SMRs by treatment group are given in Table 3 and eTable 1 in the Supplement. The drug-only group had elevated mortality rates for solid cancer (SMR, 1.31; 95% CI, 1.11-1.53), breast cancer (SMR, 1.85; 95% CI, 1.31-2.62), and nonbreast solid cancer (SMR, 1.22; 95% CI, 1.01-1.46). Elevated breast cancer mortality was observed in the drugs and RAI group (SMR, 1.16; 95% CI, 1.00-1.35). Excluding patients with prior cancers attenuated these risks. Additional exclusions (patients with TNG in the first 5 years of follow-up) had little additional influence on these results. The SMRs were not elevated in the other treatment groups. The elevated solid cancer mortality rate in the drug-only group was restricted to the first 5 years of follow-up and was no longer apparent after excluding prior cancers (eTable 1 in the Supplement). After excluding prior cancers, no consistent patterns by treatment type were observed for the other solid cancers (eTable 2 in the Supplement) apart from elevated thyroid cancer mortality in RAI-treated patients, which was no longer apparent after excluding the first 5 years of follow-up and patients with TNG.
After excluding prior cancers and adjusting for patient age, sex, and underlying diagnosis, we found no significant differences in the HR of solid cancer death across treatment groups (Table 4). Adjustment for physician-reported clinical severity of disease had almost no effect on these results (HR for drugs only, 0.98; 95% CI, 0.79-1.21; HR for surgery only, 0.85, 95% CI, 0.68-1.07; and HR for drugs and surgery, 0.98, 95% CI, 0.91-1.05). Among RAI-treated patients, risk of solid cancer death increased with total administered activity (HR, 1.08 per 370 MBq; 95% CI, 1.03-1.13 per 370 MBq; P = .001 for trend). Similar exposure-response patterns were observed for mortality from breast cancer (HR, 1.10 per 370 MBq, 95% CI, 0.99-1.22 per 370 MBq; P = .08 for trend) and nonbreast solid cancer (HR, 1.08 per 370 MBq; 95% CI, 1.03-1.13 per 370 MBq; P = .003 for trend). Results were stronger among patients receiving only RAI (no additional ATDs or surgery): HRs were 1.19 per 370 MBq (95% CI, 1.09-1.30 per 370 MBq) for solid cancer, 1.28 per 370 MBq (95% CI, 1.03-1.59 per 370 MBq) for breast cancer, and 1.18 (95% CI, 1.07-1.29 per 370 MBq) for nonbreast solid cancer mortality. None of the findings above changed notably after restricting the findings to patients with Graves disease.
Using data from a large cohort of patients treated for hyperthyroidism between the mid-1940s and mid-1960s, we evaluated the association between treatment received (RAI, ATDs, surgery, or combinations) and long-term risk of solid cancer death. After accounting for known sources of bias, we observed no significant differences in risk across treatment groups based on external and internal cohort comparisons, and notably, we found no evidence of an increased risk of solid cancer death associated with use of ATDs. Among RAI-treated patients, however, total administered activity was positively associated with risks of all solid cancers, female breast cancer, and nonbreast solid cancer death, with increases of 8% to 10% per 370 MBq. The direction and magnitude of these associations were consistent with previous findings based on estimated organ-absorbed doses.21
Caution is needed when comparing risks across treatment groups using observational study data, particularly when relying on general population rates as the comparison group. SMRs may be less than, equal to, or greater than unity (1.0) for reasons other than the treatment (eg, genetic factors, the disease, environmental and behavioral factors, and access to care) and should be interpreted in the context of other SMRs. Internal comparisons (eg, HR estimates) provide greater reassurance against bias by restricting comparisons to individuals with the same disease and allowing for covariate adjustment. Nonetheless, these approaches yielded similar results after accounting for major known sources of confounding.
Although some have argued that the lack of an unexposed comparison group was an important limitation of the previous analysis of organ and tissue absorbed doses and cancer mortality,23-25 the results of the current analysis suggest that the focus on RAI-treated patients was actually a strength because their inclusion could have biased our estimates of the radiation-associated risks.26 For instance, the surgery-only group included a higher proportion of TNG than Graves disease diagnoses (44.5% vs 7.9% in the full cohort), and the drug-only group included a high proportion of patients with baseline comorbidities, including prior cancer diagnoses, raising the possibility that these groups may have had different baseline risks of cancer death for reasons unrelated to the treatments. The elevated solid cancer mortality rate in the drug-only group, observed previously by Ron et al,10 has been suggested to indicate direct association of ATDs with cancer risk23-25; however, confounding by indication seems to be the more likely explanation. As was demonstrated in the earlier study,10 this association was no longer apparent after excluding patients with prior cancer diagnoses, and no association was observed after the first 5 years of follow-up. In addition, similar elevations in solid cancer mortality were not observed in the drugs and surgery, drugs and RAI, or drugs, surgery, and RAI treatment groups. Although numerous changes have occurred in the management of hyperthyroidism since the 1940s to 1960s, it seems plausible that patients presenting with multiple comorbidities or short life expectancies were more likely to have been managed using only ATDs rather than RAI and/or surgery. Although excluding prior cancers helped minimize this bias, it also negatively biased all the SMRs to some extent (as can be observed in Table 3) because the same exclusion could not be performed for the expected (general population) rates.
More generally, the CTTFUS cohort is not well suited to the evaluation of long-term risks associated with ATDs. The ATD formulations have changed over time, as methimazole has replaced propylthioracil as the preferred option for most patients with Graves disease because of evidence of adverse effects of propylthioracil on the liver.27,28 Propylthioracil was likely the most commonly prescribed ATD at the start of cohort follow-up. As described previously,10 no information was abstracted or retained from the medical records of these patients that would allow us to distinguish between the specific ATD types or to reconstruct the prescribed doses, and we had incomplete (and sometimes inconsistent) information on drug treatment starting and stopping dates. Considering the temporal changes in ATD treatment, more contemporary cohort studies are needed to evaluate the long-term adverse effects of methimazole specifically.
Ionizing radiation, however, is an established carcinogen.29 The breast is 1 of the most radiosensitive cancer sites, with epidemiologic studies29-34 finding consistent positive associations between absorbed dose to the breast, even in the low-dose range, and subsequent risk of breast cancer. In contrast to ATDs, data on RAI treatment dates and administered activity were nearly complete for this cohort. A previous study21 estimated that risks of solid cancer increased by 6%, breast cancer by 12%, and nonbreast solid cancer death by 5% for every 100-mGy absorbed dose to the stomach or breast. Similar estimates of 8% to 10% increased risk per 370 MBq were obtained in the current study using administered activity instead of organ doses (considering that 370 MBq corresponds to a median stomach dose of approximately 150 mGy [IQR, 140-170 mGy] and median breast dose of approximately 140 mGy [IQR, 120-160 mGy]).21 However, risk estimates based on organ doses should be more reliable because organ dose estimates accounted for clinical parameters, such as thyroid mass and percentage of uptake, in addition to administered activity.35 Previous estimates of risk per unit organ dose were consistent with those from other observational studies of radiation-exposed populations.30-34 In the current study, the lack of discernible differences in solid cancer risk across treatment groups may be partly explained by the low and narrow range of administered activities (median, 270.8 MBq [IQR, 187.8-422.5 MBq] for Graves disease and 488.1 MBq [IQR, 300.1-817.7 MBq] for TNG). Applying our estimates to currently recommended levels of RAI treatment (typically 370-555 MBq from a single application),1 the recent previous analysis21 estimated that 20 to 30 excess solid cancers may occur per 1000 patients treated as a direct result of the radiation exposure, although most excess deaths would likely occur more than 20 years after initial treatment. Thus, for many patients, the benefits of RAI may outweigh these risks.
This study has limitations. Those not already mentioned include lack of information on some potential confounding factors, such as cigarette smoking, obesity, and reproductive factors21,26; however, we found no association between treatment type (or radiation-absorbed dose to the lung21) and lung cancer mortality, suggesting that a major smoking-related bias did not occur. We lacked laboratory measures of thyroid function to assess severity of the underlying disease. The reliance on cancer mortality, as opposed to incidence, follow-up was another limitation because we could not distinguish between factors associated with cancer development vs survival. In addition, mortality follow-up is not ideal for capturing cancers with high survival rates (eg, thyroid or breast) or for studying cancer subtypes. Findings related to RAI treatment may not be generalizable to patients with thyroid cancer, who typically receive higher administered activities (>3700 MBq) in the context of postsurgical remnant ablation or therapy for presumed or known advanced or metastatic disease, potentially yielding higher or lower doses to individual organs or tissues, depending on the size and location of the residual thyroid tissue.36
This study found positive associations between total administered activity of RAI and risks of death from total solid cancer, female breast cancer, and nonbreast solid cancers. These findings were consistent with earlier results from this cohort based on organ-absorbed doses.21 The low, narrow range of total administered activity and the potential for residual confounding by risk factors not captured in this study may have hindered our ability to detect differences in risk between RAI-treated and non–RAI-treated patients. No evidence supporting an association between ATD treatment and risk of solid cancer death was found in this cohort.
Accepted for Publication: April 26, 2020.
Published: July 23, 2020. doi:10.1001/jamanetworkopen.2020.9660
Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2020 Kitahara CM et al. JAMA Network Open.
Corresponding Author: Cari M. Kitahara, PhD, MHS, Radiation Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, 9609 Medical Center Dr, Room 7E-456, Bethesda, MD 20892 (email@example.com).
Author Contributions: Dr Kitahara had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Kitahara, Preston, Berrington de Gonzalez.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Kitahara, Preston.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Kitahara, Preston, Berrington de Gonzalez.
Obtained funding: Kitahara.
Conflict of Interest Disclosures: Dr Sosa reported being a member of the data monitoring committee of the Medullary Thyroid Cancer Consortium Registry supported by GlaxoSmithKline, Novo Nordisk, Astra Zeneca, and Eli Lilly and receiving institutional research funding from Exelixis and Eli Lilly. No other disclosures were reported.
Funding/Support: This research was funded by the Intramural Research Program of the National Cancer Institute, National Institutes of Health (Drs Kitahara, Preston, and Berrington de Gonzalez).
Role of the Funder/Sponsor: The funding source had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Create a personal account or sign in to: