Abstract

Context

The waiting time for radioactive iodine therapy (WRAIT) after total thyroidectomy (TT) in patients with papillary thyroid cancer (PTC) and lymph node metastases (N1) has not been sufficiently investigated for risk of adverse outcomes.

Objective

This work aimed to estimate the effect of WRAIT on the outcomes of disease persistence and recurrence among patients with N1 PTC and investigate factors predictive of delayed radioactive iodine therapy (RAIT).

Methods

This retrospective cohort study was conducted in a university hospital. A total of 909 patients with N1 PTC were referred for RAIT between 2014 and 2018. WRAIT is the duration between TT and initial RAIT. The optimal WRAIT threshold determined using recursive partitioning analysis was used to define early and delayed RAIT. The primary end point was tumor persistence/recurrence. We compared the outcomes of patients with early and delayed RAIT using inverse probability weighting based on the propensity score.

Results

The WRAIT threshold that optimally differentiated worse long-term remission/excellent response outcomes was greater than 88 days (51% of our cohort; n = 464). WRAIT exceeding 88 days was associated with an augmented risk of disease persistence/recurrence (odds ratio, 2.47; 95% CI, 1.60-3.82) after adjustment. Predictors of delayed RAIT included residence in lower-income areas, reoperation before the initial RAIT, TT at a nonuniversity-affiliated hospital, multifocality, extrathyroidal extension, N1b disease, and pre–RAIT-stimulated thyroglobulin level less than 1 ng/mL.

Conclusion

Delayed RAIT beyond 88 days after TT in patients with N1 PTC independently increased the risk of disease persistence/recurrence. Evaluation of the predictive determinants of prolonged WRAIT may help target at-risk patients and facilitate interventions.

Thyroid cancer is the most commonly diagnosed endocrine malignancy, and its incidence is rapidly increasing worldwide, with more than 586 000 new cases in 2020 (1). Among them, 80% to 85% are papillary thyroid cancers (PTCs) (2, 3). While regional cervical lymph node (LN) metastasis occurs in approximately half of the patients with PTC and represents an adverse pathologic finding associated with a detrimental effect on oncologic outcomes (4-7), in select individuals with LN involvement, the addition of radioactive iodine therapy (RAIT) to total thyroidectomy (TT) has been shown to result in high rates of disease-free and overall survival (8-11).

In real-world practice, however, patients may experience longer waiting time for radioactive iodine therapy (WRAIT) following surgery because of various patient-, disease-, or system-related factors, including comorbidities, postoperative complications, nonadherence, increased pretreatment radiologic/pathological workup, low-efficiency health care process, or limited medical facilities for RAIT (12-14). Presently, the issue of delays for patients scheduled for RAIT has been further exacerbated by the ongoing COVID-2019 pandemic. However, the established recommendations offer little guidance regarding the appropriate waiting time for postoperative administration of RAI, regardless of low-, intermediate-, or high-risk status (15-17).

Previous retrospective studies in low-, intermediate-, or high-risk PTC patients have characterized the prognostic significance of WRAIT by using an arbitrary time-to-RAI cutoff, such as 1, 2, 3, 6, or 9 months, and different primary outcome measures (14, 18-23); thus, conclusions drawn from these studies have varied, and further investigation is warranted. An optimal WRAIT threshold in patients with node-positive PTC has not been determined. Identifying a threshold number of days precisely defining a prolonged WRAIT would help guide cancer treatment centers on the appropriate allocation and increased use of resources. Moreover, data regarding determinants of the extended delay to RAIT for PTC are scarce, and research is needed to identify susceptible patients and intervene prospectively.

In the present study, we sought to establish a time-to-RAI threshold based on clinical outcomes in patients with PTC and pathological N1 disease in a university medical center. We also analyzed the association of delayed RAIT with the risk of disease persistence or recurrence. In addition, we investigated factors predictive of extended WRAIT in this population.

Materials and Methods

Study Population

This study was approved by the institutional review board of the Zhujiang Hospital, Southern Medical University (approval No. 2021-KY-059-01); the institutional review board waived the requirement for written informed consent. We used the institution's electronic medical records system to retrieve the records of all consecutive patients with thyroid cancer admitted to our tertiary referral center between January 2014 and October 2018. We selected all patients with pathologically confirmed PTC and regional LN metastases at presentation, who underwent TT plus central neck LN dissection with or without lateral neck LN dissection, followed by RAIT (n = 2860). We excluded patients with missing or incomplete data (n = 1673) or those with less than 6 months of follow-up after the initial RAI administration (n = 106). We next excluded patients with detectable thyroglobulin (Tg) antibody (TgAb) or serum thyrotropin (TSH) levels (before the first RAIT) lower than 30 mUI/L (n = 172). The final cohort comprised 909 patients.

Data Collection and Definitions

Patient baseline demographic information, histopathologic characteristics, and treatment details, including age at diagnosis, sex, body mass index, the travel distance between the patient's residence and our facility, American Thyroid Association (ATA) risk stratification, primary tumor stage, extrathyroidal extension (ETE), bilateral tumor, BRAF V600E mutational status, node stage, number of metastatic LNs, number of LNs removed, extranodal extension (ENE), size of metastatic LN foci in the largest diameter, cumulative RAI doses, laboratory test results (TSH, TgAb, and Tg), Charlson-Deyo comorbidity score (0 or ≥1) (24), thyroid surgery–specific postoperative complications, such as hypoparathyroidism/hypocalcemia and vocal cord/fold paralysis (25), surgical hospital type (university affiliated or nonuniversity affiliated), and insurance status, were recorded from the electronic medical record. In addition, as we did not directly collect data on patient socioeconomic status (SES), 2 area-based socioeconomic measures, the annual per capita disposable income (APCDI) and the proportion of students who are receiving a first- to ninth-grade education obtained from the 2019 National Economic and Social Development Statistical Bulletins and Statistical Yearbooks of various regions (26), were used as proxies for SES. The great-circle formula calculated travel distance. The type of insurance comprised self-pay; Publicly Funded Free Medical Care for civil servants; New Cooperative Medical System for rural residents; Urban Resident Basic Medical Insurance for urban unemployed, older individuals, children, and students; and Urban Employee Basic Medical Insurance for urban employed and retired (27). Both income and educational status were stratified by respective medians, creating high income (top 50th percentile) vs low income (bottom 50th percentile) and high educational status (top 50th percentile) vs low educational status (bottom 50th percentile) variables. As defined by the seventh edition of the American Joint Commission on Cancer staging system for differentiated thyroid cancer, ETE can be classified into no, microscopic, and macroscopic ETEs (28). All patients were restaged according to the guidelines of the eighth edition of the American Joint Commission on Cancer (29). The risk of disease recurrence was classified according to the 2015 ATA risk stratification criteria (15). Serum Tg and TgAb tests were assessed using an electrochemiluminescent immunometric assay (Cobas e801 instrument, Roche Diagnostics; TgAb, catalog No. 04738578, RRID: AB_263104) with a lower detection limit of 0.04 ng/mL and 10 IU/mL, respectively. Serum TSH levels were measured using an immunoradiometric assay (Weifang Sanwei Bioengineering Group Co Ltd) with a lower detection limit of 0.01 mIU/L.

Treatment and Follow-up Protocol

Specialized head-neck surgeons performed all TTs and LN dissections in eligible individuals; all surgical specimens were carefully examined by experienced surgeons and pathologists. Routine prophylactic central-compartment neck dissection was performed in patients with TT. For patients with clinically apparent lateral neck disease or biopsy-proven disease, therapeutic lateral neck dissection of the involved compartments was performed. Voice/laryngeal assessment and serum calcium, phosphate, and parathyroid hormone levels were measured at 2 to 8 weeks postoperatively at the discretion of the surgical team in the hospital where the patients underwent the first surgery. Postoperative initial RAIT was typically administered following 3 weeks of thyroid hormone withdrawal. All individuals were advised to adhere to an iodine-restricted diet for approximately 1 to 2 weeks preceding RAIT. A TSH level of 30 mIU/L or greater before RAIT was confirmed in all individuals. The stimulated serum Tg and TgAb levels were measured immediately before the first RAIT. The choice of the initial RAI-administered activities between 30 and 100 mCi, 100 and 150 mCi, and 100 and 200 mCi was at the treating physician's discretion in low-, intermediate-, and high-risk patients, respectively.

Clinical assessments after completion of the first RAIT were scheduled at 1 and 4 months, every 3 to 6 months for the first 2 to 3 years, and then annually thereafter. The routine follow-up protocol consisted of a medical history review and clinical examination; basal serum Tg level on thyroxine therapy, TgAb, and TSH tests; and neck ultrasonography. Six to 18 months after the initial administration of RAI, a diagnostic RAI whole-body scan combined with stimulated Tg measurement in the presence of TSH of 30 mIU/L or greater through thyroid hormone withdrawal was performed in most patients. In patients with suspected persistent or newly identified evidence, additional imaging studies, such as cross-sectional imaging (computed tomography [CT] or magnetic resonance imaging [MRI]) or positron emission tomography (PET)/CT scan, were performed when required for further clarification. During follow-up, additional treatments (ie, repeated surgery or RAIT, or other therapies) for persistent or recurrent disease were administered when required based on the judgment of the attending clinicians and the consensus of the thyroid cancer multidisciplinary team. Repeated RAI treatment with empirically prescribed activity of iodine-131 (100-200 mCi) was administered to patients with persistent nonresectable disease every 6 to 12 months as long as the disease continued to concentrate RAI and an objective benefit was demonstrated (eg, decreases in lesion size and Tg level). Follow-up duration was calculated from the initial surgical treatment date to the most recent clinic visit. The cutoff date for follow-up was October 10, 2022.

Waiting Time for Radioactive Iodine Therapy and Definition of Clinical Outcome

We defined WRAIT as the interval in days from TT for thyroid cancer to the first administration of RAI. The outcome measures of our study were a biochemical or structural persistent disease, recurrence, and disease remission, based on reassessment at each follow-up visit after initial postsurgery RAIT. Patients were considered to have an excellent response or long-term disease remission status at the final known follow-up examination if they had no biochemical (suppressed Tg level of <1 ng/mL, stimulated Tg level of <10 ng/mL, and no detectable TgAb), no evidence of disease on cross-sectional imaging (ultrasonography, CT, or MRI), functional imaging (RAI whole-body scan or PET/CT), or biopsy-proven disease (cytology or histology). Patients were deemed to have persistent biochemical disease only if they displayed a suppressed Tg level of 1 ng/mL or greater, stimulated Tg level of 10 ng/mL or greater, or increasing levels of TgAbs without a structural correlate. Persistent structural disease involving local soft tissues, LNs, and distant areas was defined as evidence of disease on imaging studies or biopsy-proven disease, regardless of Tg or TgAb values. Recurrence was noted when new evidence of cancer's structural or biochemical persistence was detected after a disease-free status for at least 1 year (30). Disease-free survival (DFS) was defined as the time between initial RAI administration and disease persistence/recurrence (31).

Statistical Analysis

Continuous variables are presented as median (interquartile range [IQR]) and compared using the Wilcoxon rank sum test, and categorical variables are presented as numbers (percentages) and compared using the χ2 test. An optimal threshold for WRAIT (in days) was identified using the recursive partitioning analysis (RPA) to determine the greatest differences in the estimated risk of developing the persistent or recurrent disease (32). The cutoff was validated by iterative Wald tests for the odds of persistent or recurrent disease at increasing WRAIT thresholds of 5 days. The WRAIT producing the maximal χ2 scores was designated to dichotomize time and divide patients into receiving either early or delayed RAIT. Univariable and multivariable logistic regression models were performed to evaluate the effect of WRAIT on outcomes and identify risk factors for prolonged WRAIT based on the identified optimal cutoff. Prognostic factors were retained in the multivariable model only if they had P values less than .05 on univariable analysis. Cramer's V correlation coefficient and variance inflation factor (VIF) were calculated to investigate the collinearity between the receipt of reoperation before the first RAIT and surgical hospital type. In general, Cramer's value of 0.3 or greater denotes a strong association (33). Collinearity was considered present when dependent variables were highly correlated, and the VIF of 1 variable was at least 10 (34). In addition, to address selection bias, an inverse probability of treatment weighting (IPTW) via propensity score (PS) methods were performed (35). The PS, defined as the probability that a patient would receive either early RAIT or late RAIT, was calculated through logistic regression, including all baseline covariates. In the IPTW analysis, the PSs were used to calculate the stabilized inverse probability weights. Logistic models that used the stabilized inverse probability weights were reported. The balance of covariates between treatment groups was examined using the SDs, in which differences below 20% indicated successful balance. DFS curves were estimated with the Kaplan-Meier method and compared using the log-rank test.

To assess the robustness of the main results, 4 post hoc sensitivity analyses were performed. First, we conducted a reduced multivariable logistic regression for evaluating variables related to extended WRAIT. Second, we examined the association between longer WRAIT and incident persistent or recurrent disease in a PS-matched cohort, in which the one-to-one nearest-neighbor matching algorithm and caliper width 0.1 were applied to create a subsample controlled for potential confounders. Third, we repeated analyses stratified by age at PTC diagnosis (<55 years vs ≥55 years), ATA risk stratification (low vs intermediate vs high risk), T stage (T1-2 vs T3-4), N stage (N1a vs N1b), and different levels of stimulated Tg before the first RAIT (Tg <1 ng/mL vs 1≤ Tg ≤10 ng/mL vs Tg >10 ng/mL). Fourth, we repeated analyses using previously published WRAIT thresholds of 90-day or 180-day defining early and delayed groups. To explore whether the effect of WRAIT on the outcome was proportional to histopathological characteristics, we used causal mediation analysis to characterize the causality relationship in our retrospective study (36). RPA was performed with the use of JMP Pro 14 software (SAS Institute Inc). All other analyses were conducted with the use of R software (R Foundation, version 4.1.1). Two-sided P less than .05 denoted statistical significance.

Results

Patient Characteristics and Propensity Adjustment

Our cohort comprised 909 patients with node-positive PTC. There were 869 patients with a conventional form of PTC (95.6%), 29 with a follicular variant, 2 with a diffuse sclerosing variant, 1 with tall cell variant, 1 with columnar variant, 2 with insular variant, and 5 with a solid variant. Overall, 808 patients (88.9%) were younger than 55 years, and 63.4% were female (Table 1). A total of 391 (43%) patients had LN metastases involving the central compartment alone, 101 (11.1%) involved the lateral compartment alone, and 417 (45.9%) involved both central and lateral compartments. Of these 909 patients, the maximum sizes of the metastatic foci within an LN were less than 0.2 cm in 148 (16.3%), 0.2 to 3 cm in 756 (83.2%), and greater than 3 cm in 5 (0.6%). The median numbers of dissected and positive LNs were 15 and 5, respectively. ENE was observed in 351 (38.6%) patients. Overall, 95 patients (10.4%) underwent at least 1 reoperation after the initial surgical procedure for PTC. No patient had a history of high-dose iodine exposure (eg, intravenous contrast, amiodarone, or others) within 3 months before RAIT. The median WRAIT was 91 days (IQR, 63-188 days) post surgery. RPA and Wald statistics identified a WRAIT of 88 days for maximizing differences in the risk of persistent or recurrent disease when comparing patients with a longer time to initial RAIT with those with an equal or shorter time to initial RAIT. Delayed RAIT was defined as the use of RAI occurring more than 88 days after surgery.

Table 1.

Baseline characteristics of unweighted sample, inverse probability of treatment-weighted sample, and propensity score-matched sample based on waiting time for radioactive iodine threshold of 88 days

CharacteristicUnweighted sample, No. (%)Weighted sample, %PS-matched sample, No. (%)
WRAIT ≤88 d(n = 445)WRAIT > 88 d(n = 464)PWRAIT ≤88 d(n = 457.7)WRAIT > 88 d(n = 463.9)Standardized differenceWRAIT ≤88 d(n = 306)WRAIT > 88 d(n = 306)Standardized difference
Age at diagnosis, y.606
 <55398 (89.4)410 (88.4)85.588.7−0.097273 (89.2)272 (88.9)0.010
 ≥5547 (10.6)54 (11.6)14.511.30.09733 (10.8)34 (11.1)−0.010
Sex.785
 Female280 (62.9)296 (63.8)59.564.9−0.111190 (62.1)195 (63.7)−0.034
 Male165 (37.1)168 (36.2)40.535.10.111116 (37.9)111 (36.3)0.034
Body mass index22.7 (20.3-24.8)22.2 (19.9-24.7).05322.8 (20.2-25.4)22.6 (20.1-24.8)0.07222.2 (20.0-24.6)22.2 (20.0-24.5)0.054
Insurance type.057
 PFFMC14 (3.1)13 (2.8)2.52.6−0.0065 (1.6)9 (2.9)−0.088
 UEBMI363 (81.6)390 (84.1)81.583.9−0.065262 (85.6)250 (81.7)0.106
 URBMI8 (1.8)8 (1.7)1.61.9−0.0203 (1.0)7 (2.3)−0.103
 NCMS35 (7.9)44 (9.5)10.48.40.06823 (7.5)31 (10.1)−0.092
 Self-pay25 (5.6)9 (1.9)4.03.10.04513 (4.2)9 (2.9)0.070
Education level.847
 Low339 (76.2)356 (76.7)78.177.10.022238 (77.8)239 (78.1)−0.008
 High106 (23.8)108 (23.3)21.922.9−0.02268 (22.2)67 (21.9)0.008
Income level<.001
 Low213 (47.9)273 (58.8)57.654.20.069173 (56.5)175 (57.2)−0.013
 High232 (52.1)191 (41.2)42.445.8−0.069133 (43.5)131 (42.8)0.013
Distance, km.273
 ≤54231 (51.9)224 (48.3)46.150.1−0.079147 (48.0)152 (49.7)−0.033
 >54214 (48.1)240 (51.7)53.949.90.079159 (52.0)154 (50.3)0.033
CCI score.164
 0342 (76.9)338 (72.8)76.974.10.066231 (75.5)230 (75.2)0.008
 ≥1103 (23.1)126 (27.2)23.125.9−0.06675 (24.5)76 (24.8)−0.008
Surgical hospital type<.001
 University affiliated339 (76.2)267 (57.5)62.465.8−0.072216 (70.6)204 (66.7)0.085
 Nonuniversity affiliated106 (23.8)197 (42.5)37.634.20.07290 (29.4)102 (33.3)−0.085
Thyroid surgery–specific complications.262
 No334 (75.1)333 (71.8)75.273.30.042223 (72.9)222 (72.5)0.007
 Yes111 (24.9)131 (28.2)24.826.7−0.04283 (27.1)84 (27.5)−0.007
Receipt of reoperation before first RAIT<.001
 No435 (97.8)379 (81.7)86.189.5−0.104296 (96.7)292 (95.4)0.067
 Yes10 (2.2)85 (18.3)13.910.50.10410 (3.3)14 (4.6)−0.067
Bilateral tumor.531
 No305 (68.5)309 (66.6)70.867.70.068209 (68.3)205 (67.0)0.028
 Yes140 (31.5)155 (33.4)29.232.3−0.06897 (31.7)101 (33.0)−0.028
Extrathyroidal extension.003
 Nil269 (60.4)233 (50.2)52.154.5−0.048165 (53.9)162 (52.9)0.020
 Microscopic152 (34.2)188 (40.5)41.238.10.063119 (38.9)119 (38.9)0.000
 Macroscopic24 (5.4)43 (9.3)6.77.4−0.02622 (7.2)25 (8.2)−0.037
ATA risk stratification.049
 Low46 (10.3)48 (10.3)8.99.7−0.03035 (11.4)33 (10.8)0.021
 Intermediate358 (80.4)349 (75.2)80.278.90.033243 (79.4)235 (76.8)0.063
 High41 (9.2)67 (14.4)10.911.4−0.01528 (9.2)38 (12.4)−0.106
T stage.125
 T1322 (72.4)317 (68.3)67.869.8−0.045209 (68.3)211 (69.0)−0.014
 T285 (19.1)86 (18.5)22.319.30.07564 (20.9)60 (19.6)0.033
 T324 (5.4)33 (7.1)5.46.1−0.02820 (6.5)18 (5.9)0.027
 T414 (3.1)28 (6.0)4.54.8−0.01513 (4.2)17 (5.6)−0.061
N stage<.001
 N1a224 (50.3)167 (36.0)41.042.8−0.035130 (42.5)131 (42.8)−0.007
 N1b221 (49.7)297 (64.0)59.057.20.035176 (57.5)175 (57.2)0.007
M stage.753
 M0424 (95.3)440 (94.8)96.195.80.015297 (97.1)294 (96.1)0.054
 M121 (4.7)24 (5.2)3.94.2−0.0159 (2.9)12 (3.9)−0.054
Multifocality<.001
 No351 (78.9)319 (68.8)76.274.60.038237 (77.5)224 (73.2)0.099
 Yes94 (21.1)145 (31.2)23.825.4−0.03869 (22.5)82 (26.8)−0.099
BRAF V600E mutation.118
 Negative82 (18.4)99 (21.3)19.718.90.01956 (18.3)55 (18.0)0.008
 Positive275 (61.8)296 (63.8)63.264.3−0.022198 (64.7)207 (67.6)−0.062
 Unknown88 (19.8)69 (14.9)17.116.80.00852 (17.0)44 (14.4)0.072
Extranodal extension.7
 No276 (62.0)282 (60.8)56.360.4−0.083191 (62.4)188 (61.4)0.020
 Yes169 (38.0)182 (39.2)43.739.60.083115 (37.6)118 (38.6)−0.020
Size of metastatic LN focia, cm0.4 (0.2–0.9)0.5 (0.2–1.1).0270.4 (0.25-1.1)0.5 (0.2-1.0)0.0050.4 (0.2-1.0)0.5 (0.2-1.0)−0.011
No. of LNs examined13 (6–25)16 (8–29.2).00416 (8-26)15 (7-28)−0.02514 (7-27)15 (7-27)0.021
No. of metastatic LNs5 (3–9)5 (3–10).4085 (3-10)5 (3-10)0.0195 (3-9)5 (3-9)0.009
Cumulative RAI dose, mCi150 (140–180)150 (100–180)<.001150 (135-180)150 (120-248)−0.018150 (136.2-174)150 (102-180)0.066
sTg before first RAIT, ng/mL<.001
  <179 (17.8)125 (26.9)18.221.8−0.08957 (18.6)76 (24.8)−0.151
 1-10208 (46.7)164 (35.3)40.239.70.010145 (47.4)117 (38.2)0.186
  >10158 (35.5)175 (37.7)41.638.50.063104 (34.0)113 (36.9)−0.062
CharacteristicUnweighted sample, No. (%)Weighted sample, %PS-matched sample, No. (%)
WRAIT ≤88 d(n = 445)WRAIT > 88 d(n = 464)PWRAIT ≤88 d(n = 457.7)WRAIT > 88 d(n = 463.9)Standardized differenceWRAIT ≤88 d(n = 306)WRAIT > 88 d(n = 306)Standardized difference
Age at diagnosis, y.606
 <55398 (89.4)410 (88.4)85.588.7−0.097273 (89.2)272 (88.9)0.010
 ≥5547 (10.6)54 (11.6)14.511.30.09733 (10.8)34 (11.1)−0.010
Sex.785
 Female280 (62.9)296 (63.8)59.564.9−0.111190 (62.1)195 (63.7)−0.034
 Male165 (37.1)168 (36.2)40.535.10.111116 (37.9)111 (36.3)0.034
Body mass index22.7 (20.3-24.8)22.2 (19.9-24.7).05322.8 (20.2-25.4)22.6 (20.1-24.8)0.07222.2 (20.0-24.6)22.2 (20.0-24.5)0.054
Insurance type.057
 PFFMC14 (3.1)13 (2.8)2.52.6−0.0065 (1.6)9 (2.9)−0.088
 UEBMI363 (81.6)390 (84.1)81.583.9−0.065262 (85.6)250 (81.7)0.106
 URBMI8 (1.8)8 (1.7)1.61.9−0.0203 (1.0)7 (2.3)−0.103
 NCMS35 (7.9)44 (9.5)10.48.40.06823 (7.5)31 (10.1)−0.092
 Self-pay25 (5.6)9 (1.9)4.03.10.04513 (4.2)9 (2.9)0.070
Education level.847
 Low339 (76.2)356 (76.7)78.177.10.022238 (77.8)239 (78.1)−0.008
 High106 (23.8)108 (23.3)21.922.9−0.02268 (22.2)67 (21.9)0.008
Income level<.001
 Low213 (47.9)273 (58.8)57.654.20.069173 (56.5)175 (57.2)−0.013
 High232 (52.1)191 (41.2)42.445.8−0.069133 (43.5)131 (42.8)0.013
Distance, km.273
 ≤54231 (51.9)224 (48.3)46.150.1−0.079147 (48.0)152 (49.7)−0.033
 >54214 (48.1)240 (51.7)53.949.90.079159 (52.0)154 (50.3)0.033
CCI score.164
 0342 (76.9)338 (72.8)76.974.10.066231 (75.5)230 (75.2)0.008
 ≥1103 (23.1)126 (27.2)23.125.9−0.06675 (24.5)76 (24.8)−0.008
Surgical hospital type<.001
 University affiliated339 (76.2)267 (57.5)62.465.8−0.072216 (70.6)204 (66.7)0.085
 Nonuniversity affiliated106 (23.8)197 (42.5)37.634.20.07290 (29.4)102 (33.3)−0.085
Thyroid surgery–specific complications.262
 No334 (75.1)333 (71.8)75.273.30.042223 (72.9)222 (72.5)0.007
 Yes111 (24.9)131 (28.2)24.826.7−0.04283 (27.1)84 (27.5)−0.007
Receipt of reoperation before first RAIT<.001
 No435 (97.8)379 (81.7)86.189.5−0.104296 (96.7)292 (95.4)0.067
 Yes10 (2.2)85 (18.3)13.910.50.10410 (3.3)14 (4.6)−0.067
Bilateral tumor.531
 No305 (68.5)309 (66.6)70.867.70.068209 (68.3)205 (67.0)0.028
 Yes140 (31.5)155 (33.4)29.232.3−0.06897 (31.7)101 (33.0)−0.028
Extrathyroidal extension.003
 Nil269 (60.4)233 (50.2)52.154.5−0.048165 (53.9)162 (52.9)0.020
 Microscopic152 (34.2)188 (40.5)41.238.10.063119 (38.9)119 (38.9)0.000
 Macroscopic24 (5.4)43 (9.3)6.77.4−0.02622 (7.2)25 (8.2)−0.037
ATA risk stratification.049
 Low46 (10.3)48 (10.3)8.99.7−0.03035 (11.4)33 (10.8)0.021
 Intermediate358 (80.4)349 (75.2)80.278.90.033243 (79.4)235 (76.8)0.063
 High41 (9.2)67 (14.4)10.911.4−0.01528 (9.2)38 (12.4)−0.106
T stage.125
 T1322 (72.4)317 (68.3)67.869.8−0.045209 (68.3)211 (69.0)−0.014
 T285 (19.1)86 (18.5)22.319.30.07564 (20.9)60 (19.6)0.033
 T324 (5.4)33 (7.1)5.46.1−0.02820 (6.5)18 (5.9)0.027
 T414 (3.1)28 (6.0)4.54.8−0.01513 (4.2)17 (5.6)−0.061
N stage<.001
 N1a224 (50.3)167 (36.0)41.042.8−0.035130 (42.5)131 (42.8)−0.007
 N1b221 (49.7)297 (64.0)59.057.20.035176 (57.5)175 (57.2)0.007
M stage.753
 M0424 (95.3)440 (94.8)96.195.80.015297 (97.1)294 (96.1)0.054
 M121 (4.7)24 (5.2)3.94.2−0.0159 (2.9)12 (3.9)−0.054
Multifocality<.001
 No351 (78.9)319 (68.8)76.274.60.038237 (77.5)224 (73.2)0.099
 Yes94 (21.1)145 (31.2)23.825.4−0.03869 (22.5)82 (26.8)−0.099
BRAF V600E mutation.118
 Negative82 (18.4)99 (21.3)19.718.90.01956 (18.3)55 (18.0)0.008
 Positive275 (61.8)296 (63.8)63.264.3−0.022198 (64.7)207 (67.6)−0.062
 Unknown88 (19.8)69 (14.9)17.116.80.00852 (17.0)44 (14.4)0.072
Extranodal extension.7
 No276 (62.0)282 (60.8)56.360.4−0.083191 (62.4)188 (61.4)0.020
 Yes169 (38.0)182 (39.2)43.739.60.083115 (37.6)118 (38.6)−0.020
Size of metastatic LN focia, cm0.4 (0.2–0.9)0.5 (0.2–1.1).0270.4 (0.25-1.1)0.5 (0.2-1.0)0.0050.4 (0.2-1.0)0.5 (0.2-1.0)−0.011
No. of LNs examined13 (6–25)16 (8–29.2).00416 (8-26)15 (7-28)−0.02514 (7-27)15 (7-27)0.021
No. of metastatic LNs5 (3–9)5 (3–10).4085 (3-10)5 (3-10)0.0195 (3-9)5 (3-9)0.009
Cumulative RAI dose, mCi150 (140–180)150 (100–180)<.001150 (135-180)150 (120-248)−0.018150 (136.2-174)150 (102-180)0.066
sTg before first RAIT, ng/mL<.001
  <179 (17.8)125 (26.9)18.221.8−0.08957 (18.6)76 (24.8)−0.151
 1-10208 (46.7)164 (35.3)40.239.70.010145 (47.4)117 (38.2)0.186
  >10158 (35.5)175 (37.7)41.638.50.063104 (34.0)113 (36.9)−0.062

Continuous variables are presented as medians with interquartile ranges.

Bold indicates statistically significant P values.

Abbreviations: ATA, American Thyroid Association; CCI, Charlson Comorbidity Index; LN, lymph node; NCMS, New Cooperative Medical System; PS, propensity score; RAI, radioactive iodine; RAIT, radioactive iodine therapy; sTg, stimulated thyroglobulin; PFFMC, Publicly Funded Free Medical Care; UEBMI, Urban Employee Basic Medical Insurance; URBMI, Urban Resident Basic Medical Insurance; WRAIT, waiting time for radioactive iodine therapy.

The longest diameter of the metastatic LN foci.

Table 1.

Baseline characteristics of unweighted sample, inverse probability of treatment-weighted sample, and propensity score-matched sample based on waiting time for radioactive iodine threshold of 88 days

CharacteristicUnweighted sample, No. (%)Weighted sample, %PS-matched sample, No. (%)
WRAIT ≤88 d(n = 445)WRAIT > 88 d(n = 464)PWRAIT ≤88 d(n = 457.7)WRAIT > 88 d(n = 463.9)Standardized differenceWRAIT ≤88 d(n = 306)WRAIT > 88 d(n = 306)Standardized difference
Age at diagnosis, y.606
 <55398 (89.4)410 (88.4)85.588.7−0.097273 (89.2)272 (88.9)0.010
 ≥5547 (10.6)54 (11.6)14.511.30.09733 (10.8)34 (11.1)−0.010
Sex.785
 Female280 (62.9)296 (63.8)59.564.9−0.111190 (62.1)195 (63.7)−0.034
 Male165 (37.1)168 (36.2)40.535.10.111116 (37.9)111 (36.3)0.034
Body mass index22.7 (20.3-24.8)22.2 (19.9-24.7).05322.8 (20.2-25.4)22.6 (20.1-24.8)0.07222.2 (20.0-24.6)22.2 (20.0-24.5)0.054
Insurance type.057
 PFFMC14 (3.1)13 (2.8)2.52.6−0.0065 (1.6)9 (2.9)−0.088
 UEBMI363 (81.6)390 (84.1)81.583.9−0.065262 (85.6)250 (81.7)0.106
 URBMI8 (1.8)8 (1.7)1.61.9−0.0203 (1.0)7 (2.3)−0.103
 NCMS35 (7.9)44 (9.5)10.48.40.06823 (7.5)31 (10.1)−0.092
 Self-pay25 (5.6)9 (1.9)4.03.10.04513 (4.2)9 (2.9)0.070
Education level.847
 Low339 (76.2)356 (76.7)78.177.10.022238 (77.8)239 (78.1)−0.008
 High106 (23.8)108 (23.3)21.922.9−0.02268 (22.2)67 (21.9)0.008
Income level<.001
 Low213 (47.9)273 (58.8)57.654.20.069173 (56.5)175 (57.2)−0.013
 High232 (52.1)191 (41.2)42.445.8−0.069133 (43.5)131 (42.8)0.013
Distance, km.273
 ≤54231 (51.9)224 (48.3)46.150.1−0.079147 (48.0)152 (49.7)−0.033
 >54214 (48.1)240 (51.7)53.949.90.079159 (52.0)154 (50.3)0.033
CCI score.164
 0342 (76.9)338 (72.8)76.974.10.066231 (75.5)230 (75.2)0.008
 ≥1103 (23.1)126 (27.2)23.125.9−0.06675 (24.5)76 (24.8)−0.008
Surgical hospital type<.001
 University affiliated339 (76.2)267 (57.5)62.465.8−0.072216 (70.6)204 (66.7)0.085
 Nonuniversity affiliated106 (23.8)197 (42.5)37.634.20.07290 (29.4)102 (33.3)−0.085
Thyroid surgery–specific complications.262
 No334 (75.1)333 (71.8)75.273.30.042223 (72.9)222 (72.5)0.007
 Yes111 (24.9)131 (28.2)24.826.7−0.04283 (27.1)84 (27.5)−0.007
Receipt of reoperation before first RAIT<.001
 No435 (97.8)379 (81.7)86.189.5−0.104296 (96.7)292 (95.4)0.067
 Yes10 (2.2)85 (18.3)13.910.50.10410 (3.3)14 (4.6)−0.067
Bilateral tumor.531
 No305 (68.5)309 (66.6)70.867.70.068209 (68.3)205 (67.0)0.028
 Yes140 (31.5)155 (33.4)29.232.3−0.06897 (31.7)101 (33.0)−0.028
Extrathyroidal extension.003
 Nil269 (60.4)233 (50.2)52.154.5−0.048165 (53.9)162 (52.9)0.020
 Microscopic152 (34.2)188 (40.5)41.238.10.063119 (38.9)119 (38.9)0.000
 Macroscopic24 (5.4)43 (9.3)6.77.4−0.02622 (7.2)25 (8.2)−0.037
ATA risk stratification.049
 Low46 (10.3)48 (10.3)8.99.7−0.03035 (11.4)33 (10.8)0.021
 Intermediate358 (80.4)349 (75.2)80.278.90.033243 (79.4)235 (76.8)0.063
 High41 (9.2)67 (14.4)10.911.4−0.01528 (9.2)38 (12.4)−0.106
T stage.125
 T1322 (72.4)317 (68.3)67.869.8−0.045209 (68.3)211 (69.0)−0.014
 T285 (19.1)86 (18.5)22.319.30.07564 (20.9)60 (19.6)0.033
 T324 (5.4)33 (7.1)5.46.1−0.02820 (6.5)18 (5.9)0.027
 T414 (3.1)28 (6.0)4.54.8−0.01513 (4.2)17 (5.6)−0.061
N stage<.001
 N1a224 (50.3)167 (36.0)41.042.8−0.035130 (42.5)131 (42.8)−0.007
 N1b221 (49.7)297 (64.0)59.057.20.035176 (57.5)175 (57.2)0.007
M stage.753
 M0424 (95.3)440 (94.8)96.195.80.015297 (97.1)294 (96.1)0.054
 M121 (4.7)24 (5.2)3.94.2−0.0159 (2.9)12 (3.9)−0.054
Multifocality<.001
 No351 (78.9)319 (68.8)76.274.60.038237 (77.5)224 (73.2)0.099
 Yes94 (21.1)145 (31.2)23.825.4−0.03869 (22.5)82 (26.8)−0.099
BRAF V600E mutation.118
 Negative82 (18.4)99 (21.3)19.718.90.01956 (18.3)55 (18.0)0.008
 Positive275 (61.8)296 (63.8)63.264.3−0.022198 (64.7)207 (67.6)−0.062
 Unknown88 (19.8)69 (14.9)17.116.80.00852 (17.0)44 (14.4)0.072
Extranodal extension.7
 No276 (62.0)282 (60.8)56.360.4−0.083191 (62.4)188 (61.4)0.020
 Yes169 (38.0)182 (39.2)43.739.60.083115 (37.6)118 (38.6)−0.020
Size of metastatic LN focia, cm0.4 (0.2–0.9)0.5 (0.2–1.1).0270.4 (0.25-1.1)0.5 (0.2-1.0)0.0050.4 (0.2-1.0)0.5 (0.2-1.0)−0.011
No. of LNs examined13 (6–25)16 (8–29.2).00416 (8-26)15 (7-28)−0.02514 (7-27)15 (7-27)0.021
No. of metastatic LNs5 (3–9)5 (3–10).4085 (3-10)5 (3-10)0.0195 (3-9)5 (3-9)0.009
Cumulative RAI dose, mCi150 (140–180)150 (100–180)<.001150 (135-180)150 (120-248)−0.018150 (136.2-174)150 (102-180)0.066
sTg before first RAIT, ng/mL<.001
  <179 (17.8)125 (26.9)18.221.8−0.08957 (18.6)76 (24.8)−0.151
 1-10208 (46.7)164 (35.3)40.239.70.010145 (47.4)117 (38.2)0.186
  >10158 (35.5)175 (37.7)41.638.50.063104 (34.0)113 (36.9)−0.062
CharacteristicUnweighted sample, No. (%)Weighted sample, %PS-matched sample, No. (%)
WRAIT ≤88 d(n = 445)WRAIT > 88 d(n = 464)PWRAIT ≤88 d(n = 457.7)WRAIT > 88 d(n = 463.9)Standardized differenceWRAIT ≤88 d(n = 306)WRAIT > 88 d(n = 306)Standardized difference
Age at diagnosis, y.606
 <55398 (89.4)410 (88.4)85.588.7−0.097273 (89.2)272 (88.9)0.010
 ≥5547 (10.6)54 (11.6)14.511.30.09733 (10.8)34 (11.1)−0.010
Sex.785
 Female280 (62.9)296 (63.8)59.564.9−0.111190 (62.1)195 (63.7)−0.034
 Male165 (37.1)168 (36.2)40.535.10.111116 (37.9)111 (36.3)0.034
Body mass index22.7 (20.3-24.8)22.2 (19.9-24.7).05322.8 (20.2-25.4)22.6 (20.1-24.8)0.07222.2 (20.0-24.6)22.2 (20.0-24.5)0.054
Insurance type.057
 PFFMC14 (3.1)13 (2.8)2.52.6−0.0065 (1.6)9 (2.9)−0.088
 UEBMI363 (81.6)390 (84.1)81.583.9−0.065262 (85.6)250 (81.7)0.106
 URBMI8 (1.8)8 (1.7)1.61.9−0.0203 (1.0)7 (2.3)−0.103
 NCMS35 (7.9)44 (9.5)10.48.40.06823 (7.5)31 (10.1)−0.092
 Self-pay25 (5.6)9 (1.9)4.03.10.04513 (4.2)9 (2.9)0.070
Education level.847
 Low339 (76.2)356 (76.7)78.177.10.022238 (77.8)239 (78.1)−0.008
 High106 (23.8)108 (23.3)21.922.9−0.02268 (22.2)67 (21.9)0.008
Income level<.001
 Low213 (47.9)273 (58.8)57.654.20.069173 (56.5)175 (57.2)−0.013
 High232 (52.1)191 (41.2)42.445.8−0.069133 (43.5)131 (42.8)0.013
Distance, km.273
 ≤54231 (51.9)224 (48.3)46.150.1−0.079147 (48.0)152 (49.7)−0.033
 >54214 (48.1)240 (51.7)53.949.90.079159 (52.0)154 (50.3)0.033
CCI score.164
 0342 (76.9)338 (72.8)76.974.10.066231 (75.5)230 (75.2)0.008
 ≥1103 (23.1)126 (27.2)23.125.9−0.06675 (24.5)76 (24.8)−0.008
Surgical hospital type<.001
 University affiliated339 (76.2)267 (57.5)62.465.8−0.072216 (70.6)204 (66.7)0.085
 Nonuniversity affiliated106 (23.8)197 (42.5)37.634.20.07290 (29.4)102 (33.3)−0.085
Thyroid surgery–specific complications.262
 No334 (75.1)333 (71.8)75.273.30.042223 (72.9)222 (72.5)0.007
 Yes111 (24.9)131 (28.2)24.826.7−0.04283 (27.1)84 (27.5)−0.007
Receipt of reoperation before first RAIT<.001
 No435 (97.8)379 (81.7)86.189.5−0.104296 (96.7)292 (95.4)0.067
 Yes10 (2.2)85 (18.3)13.910.50.10410 (3.3)14 (4.6)−0.067
Bilateral tumor.531
 No305 (68.5)309 (66.6)70.867.70.068209 (68.3)205 (67.0)0.028
 Yes140 (31.5)155 (33.4)29.232.3−0.06897 (31.7)101 (33.0)−0.028
Extrathyroidal extension.003
 Nil269 (60.4)233 (50.2)52.154.5−0.048165 (53.9)162 (52.9)0.020
 Microscopic152 (34.2)188 (40.5)41.238.10.063119 (38.9)119 (38.9)0.000
 Macroscopic24 (5.4)43 (9.3)6.77.4−0.02622 (7.2)25 (8.2)−0.037
ATA risk stratification.049
 Low46 (10.3)48 (10.3)8.99.7−0.03035 (11.4)33 (10.8)0.021
 Intermediate358 (80.4)349 (75.2)80.278.90.033243 (79.4)235 (76.8)0.063
 High41 (9.2)67 (14.4)10.911.4−0.01528 (9.2)38 (12.4)−0.106
T stage.125
 T1322 (72.4)317 (68.3)67.869.8−0.045209 (68.3)211 (69.0)−0.014
 T285 (19.1)86 (18.5)22.319.30.07564 (20.9)60 (19.6)0.033
 T324 (5.4)33 (7.1)5.46.1−0.02820 (6.5)18 (5.9)0.027
 T414 (3.1)28 (6.0)4.54.8−0.01513 (4.2)17 (5.6)−0.061
N stage<.001
 N1a224 (50.3)167 (36.0)41.042.8−0.035130 (42.5)131 (42.8)−0.007
 N1b221 (49.7)297 (64.0)59.057.20.035176 (57.5)175 (57.2)0.007
M stage.753
 M0424 (95.3)440 (94.8)96.195.80.015297 (97.1)294 (96.1)0.054
 M121 (4.7)24 (5.2)3.94.2−0.0159 (2.9)12 (3.9)−0.054
Multifocality<.001
 No351 (78.9)319 (68.8)76.274.60.038237 (77.5)224 (73.2)0.099
 Yes94 (21.1)145 (31.2)23.825.4−0.03869 (22.5)82 (26.8)−0.099
BRAF V600E mutation.118
 Negative82 (18.4)99 (21.3)19.718.90.01956 (18.3)55 (18.0)0.008
 Positive275 (61.8)296 (63.8)63.264.3−0.022198 (64.7)207 (67.6)−0.062
 Unknown88 (19.8)69 (14.9)17.116.80.00852 (17.0)44 (14.4)0.072
Extranodal extension.7
 No276 (62.0)282 (60.8)56.360.4−0.083191 (62.4)188 (61.4)0.020
 Yes169 (38.0)182 (39.2)43.739.60.083115 (37.6)118 (38.6)−0.020
Size of metastatic LN focia, cm0.4 (0.2–0.9)0.5 (0.2–1.1).0270.4 (0.25-1.1)0.5 (0.2-1.0)0.0050.4 (0.2-1.0)0.5 (0.2-1.0)−0.011
No. of LNs examined13 (6–25)16 (8–29.2).00416 (8-26)15 (7-28)−0.02514 (7-27)15 (7-27)0.021
No. of metastatic LNs5 (3–9)5 (3–10).4085 (3-10)5 (3-10)0.0195 (3-9)5 (3-9)0.009
Cumulative RAI dose, mCi150 (140–180)150 (100–180)<.001150 (135-180)150 (120-248)−0.018150 (136.2-174)150 (102-180)0.066
sTg before first RAIT, ng/mL<.001
  <179 (17.8)125 (26.9)18.221.8−0.08957 (18.6)76 (24.8)−0.151
 1-10208 (46.7)164 (35.3)40.239.70.010145 (47.4)117 (38.2)0.186
  >10158 (35.5)175 (37.7)41.638.50.063104 (34.0)113 (36.9)−0.062

Continuous variables are presented as medians with interquartile ranges.

Bold indicates statistically significant P values.

Abbreviations: ATA, American Thyroid Association; CCI, Charlson Comorbidity Index; LN, lymph node; NCMS, New Cooperative Medical System; PS, propensity score; RAI, radioactive iodine; RAIT, radioactive iodine therapy; sTg, stimulated thyroglobulin; PFFMC, Publicly Funded Free Medical Care; UEBMI, Urban Employee Basic Medical Insurance; URBMI, Urban Resident Basic Medical Insurance; WRAIT, waiting time for radioactive iodine therapy.

The longest diameter of the metastatic LN foci.

Based on the estimated threshold, 445 patients (49%) received early RAIT, and 464 patients (51%) received delayed RAIT (see Table 1). Before weighting, patients who received delayed initial RAIT were significantly more likely to undergo TT at a nonuniversity-affiliated hospital and reoperation before the initial RAIT; have a lower cumulative radioiodine dose; and present with ETE, N1b disease, multifocality, an ATA high risk of recurrence, a larger size of the largest metastatic focus in metastatic LNs, an increasing number of LNs examined, and a lower stimulated Tg value before the initial RAIT (all P < .05). In terms of urban and rural regions metrics of SES, patients who underwent delayed RAIT tended to live in regions with a lower APCDI level (P < .001). After IPTW adjustment, the 2cohorts were well-balanced on these and the other baseline factors, with all SD less than 20%.

Factors Associated With Delayed Postoperative Administration of Radioactive Iodine

Following multivariable adjustment (Table 2), independent factors predictive of delayed WRAIT longer than 88 days included reoperation before the initial RAIT, TT at a nonuniversity-affiliated hospital, ETE, multifocality, and stage N1b disease. In contrast, patients who lived in areas with a higher income level or presented with stimulated Tg of 1 ng/mL or greater before the first RAIT had a higher likelihood of initial RAI administration before 88 days. In addition, no differences in the WRAIT were found concerning T stage, BRAF V600E mutation, size of the largest metastatic focus in metastatic LNs, or number of LNs examined. Analysis was repeated in a reduced model for robustness with similar results (see Table 2). Furthermore, reoperation before the first RAIT displayed a significantly weak correlation with surgical hospital type (Cramer's V correlation coefficient = 0.1933; P < .0001), and the VIFs were 1.08 for reoperation before the first RAIT and 1.07 for surgical hospital type (Supplementary e-Tables S1 and S2) (37), showing no sign of collinearity between these 2 predictors. Finally, the key clinical characteristics of ETE, multifocality, and N1b disease were compared according to surgical hospital type, receipt of reoperation before initial RAI administration, and income level. The results are presented in Supplementary e-Table S3 (37). Compared with patients who did not undergo reoperation, a higher percentage of patients who underwent reoperation had stage N1b (14.9% vs 4.6%; P < .0001) and multifocal (13.8% vs 9.3%; P = .0482) diseases. The prevalence of reoperation (8.4%) in patients with ETE was similar to the frequency of reoperation (12.2%) in patients without ETE (P = .063). When comparing patients who underwent thyroid surgery in a nonuniversity-affiliated hospital with those in a university-affiliated hospital, the proportions of ETE, multifocality, and N1b disease significantly differed at 37.8% and 29.7%, 39.8% and 31.0%, and 36.7% and 28.9%, respectively (all P < .05). No statistically significant association was noted between the APCDI and any of these 3 tumor features (all P > .1).

Table 2.

Multivariable logistic regression model of factors potentially predictive of prolonged waiting time for radioactive iodine therapy based on the optimal threshold of 88 days

CharacteristicFull modelReduced model
OR (95% CI)POR (95% CI)P
Receipt of reoperation before first RAIT (Ref: no)9.33 (4.65-18.73)<.0018.36 (4.20-16.65)<.001
 Yes
Income level (Ref: low)
 High0.64 (0.48-0.86).0030.66 (0.50-0.87).004
Surgical hospital type (Ref: university affiliated)
 Nonuniversity affiliated1.92 (1.40-2.62)<.0011.90 (1.40-2.59)<.001
Extrathyroidal extension (Ref: none)
 Present1.45 (1.07-1.96).0161.53 (1.15-2.04).004
T stage (Ref: T1)
 T20.97 (0.67-1.40).875NANA
 T31.15 (0.63-2.08).651NANA
 T41.40 (0.68-2.89).366NANA
N stage (Ref: N1a)
 N1b1.63 (1.15-2.32).0061.49 (1.11-1.99).007
Multifocality (Ref: no)
 Yes1.44 (1.03-2.01).0341.42 (1.03-1.98).034
BRAF V600E mutation (Ref: negative)
 Positive1.26 (0.86-1.84).239NANA
Size of metastatic LN focia, cm0.94 (0.71-1.25).674NANA
No. of LNs examined1.00 (0.99-1.01).623NANA
sTg before first RAIT, ng/mL (Ref: <1)
 1-100.44 (0.30-0.64)<.0010.45 (0.31-0.66)<.001
 >100.49 (0.33-0.73)<.0010.48 (0.33-0.71)<.001
CharacteristicFull modelReduced model
OR (95% CI)POR (95% CI)P
Receipt of reoperation before first RAIT (Ref: no)9.33 (4.65-18.73)<.0018.36 (4.20-16.65)<.001
 Yes
Income level (Ref: low)
 High0.64 (0.48-0.86).0030.66 (0.50-0.87).004
Surgical hospital type (Ref: university affiliated)
 Nonuniversity affiliated1.92 (1.40-2.62)<.0011.90 (1.40-2.59)<.001
Extrathyroidal extension (Ref: none)
 Present1.45 (1.07-1.96).0161.53 (1.15-2.04).004
T stage (Ref: T1)
 T20.97 (0.67-1.40).875NANA
 T31.15 (0.63-2.08).651NANA
 T41.40 (0.68-2.89).366NANA
N stage (Ref: N1a)
 N1b1.63 (1.15-2.32).0061.49 (1.11-1.99).007
Multifocality (Ref: no)
 Yes1.44 (1.03-2.01).0341.42 (1.03-1.98).034
BRAF V600E mutation (Ref: negative)
 Positive1.26 (0.86-1.84).239NANA
Size of metastatic LN focia, cm0.94 (0.71-1.25).674NANA
No. of LNs examined1.00 (0.99-1.01).623NANA
sTg before first RAIT, ng/mL (Ref: <1)
 1-100.44 (0.30-0.64)<.0010.45 (0.31-0.66)<.001
 >100.49 (0.33-0.73)<.0010.48 (0.33-0.71)<.001

Bold indicates statistically significant P values.

Abbreviations: LN, lymph node; NA, not applicable; OR, odds ratio; Ref, reference; sTg, stimulated thyroglobulin.

The longest diameter of the metastatic LN foci.

Table 2.

Multivariable logistic regression model of factors potentially predictive of prolonged waiting time for radioactive iodine therapy based on the optimal threshold of 88 days

CharacteristicFull modelReduced model
OR (95% CI)POR (95% CI)P
Receipt of reoperation before first RAIT (Ref: no)9.33 (4.65-18.73)<.0018.36 (4.20-16.65)<.001
 Yes
Income level (Ref: low)
 High0.64 (0.48-0.86).0030.66 (0.50-0.87).004
Surgical hospital type (Ref: university affiliated)
 Nonuniversity affiliated1.92 (1.40-2.62)<.0011.90 (1.40-2.59)<.001
Extrathyroidal extension (Ref: none)
 Present1.45 (1.07-1.96).0161.53 (1.15-2.04).004
T stage (Ref: T1)
 T20.97 (0.67-1.40).875NANA
 T31.15 (0.63-2.08).651NANA
 T41.40 (0.68-2.89).366NANA
N stage (Ref: N1a)
 N1b1.63 (1.15-2.32).0061.49 (1.11-1.99).007
Multifocality (Ref: no)
 Yes1.44 (1.03-2.01).0341.42 (1.03-1.98).034
BRAF V600E mutation (Ref: negative)
 Positive1.26 (0.86-1.84).239NANA
Size of metastatic LN focia, cm0.94 (0.71-1.25).674NANA
No. of LNs examined1.00 (0.99-1.01).623NANA
sTg before first RAIT, ng/mL (Ref: <1)
 1-100.44 (0.30-0.64)<.0010.45 (0.31-0.66)<.001
 >100.49 (0.33-0.73)<.0010.48 (0.33-0.71)<.001
CharacteristicFull modelReduced model
OR (95% CI)POR (95% CI)P
Receipt of reoperation before first RAIT (Ref: no)9.33 (4.65-18.73)<.0018.36 (4.20-16.65)<.001
 Yes
Income level (Ref: low)
 High0.64 (0.48-0.86).0030.66 (0.50-0.87).004
Surgical hospital type (Ref: university affiliated)
 Nonuniversity affiliated1.92 (1.40-2.62)<.0011.90 (1.40-2.59)<.001
Extrathyroidal extension (Ref: none)
 Present1.45 (1.07-1.96).0161.53 (1.15-2.04).004
T stage (Ref: T1)
 T20.97 (0.67-1.40).875NANA
 T31.15 (0.63-2.08).651NANA
 T41.40 (0.68-2.89).366NANA
N stage (Ref: N1a)
 N1b1.63 (1.15-2.32).0061.49 (1.11-1.99).007
Multifocality (Ref: no)
 Yes1.44 (1.03-2.01).0341.42 (1.03-1.98).034
BRAF V600E mutation (Ref: negative)
 Positive1.26 (0.86-1.84).239NANA
Size of metastatic LN focia, cm0.94 (0.71-1.25).674NANA
No. of LNs examined1.00 (0.99-1.01).623NANA
sTg before first RAIT, ng/mL (Ref: <1)
 1-100.44 (0.30-0.64)<.0010.45 (0.31-0.66)<.001
 >100.49 (0.33-0.73)<.0010.48 (0.33-0.71)<.001

Bold indicates statistically significant P values.

Abbreviations: LN, lymph node; NA, not applicable; OR, odds ratio; Ref, reference; sTg, stimulated thyroglobulin.

The longest diameter of the metastatic LN foci.

Prolonged Waiting Time for Radioactive Iodine Therapy and Risk of Persistent or Recurrent Disease

During a median follow-up time of 53.2 months (IQR, 35.6-68.6 months; range, 7.1-193.7 months), 652 patients achieved long-term complete remission or excellent response to therapy. A total of 257 patients (28.3%) showed persistent (n = 218) or recurrent (n = 39) disease at the end of follow-up. Of 257 patients, 65 (7.2%), 153 (16.8%), and 7 (0.8%) patients had a biochemical incomplete response, a structural incomplete response, and a biochemical recurrence, respectively. None of the patients died. When classified based on the 2015 ATA risk stratification system, 12.7% (12/94), 26.3% (186/707), and 54.6% (59/108) of the low-, intermediate-, and high-risk patients (P < .001), respectively, had persistent or recurrent disease detected during follow-up.

Before applying the PS methods, 21.8% (97/445) of patients treated with early RAIT and 34.5% (160/464) of patients treated with delayed RAIT presented with a disease event. The crude 5-year DFS rate was 62.48% for patients treated with delayed RAIT vs 77.03% for patients treated with early RAIT (log-rank P < .0001) (Fig. 1A). On univariable logistic regression analysis (Table 3 and Supplementary e-Table S4) (37), WRAIT was a significant prognostic factor for disease nonremission (odds ratio [OR], 1.89; P < .001). After adjusting for confounding variables (see Table 3 and Supplementary e-Table S4) (37), WRAIT remained an independent prognostic factor for disease nonremission (OR, 2.50; P < .001). Among the patients in the IPTW sample, the crude 5-year DFS rate among patients treated with delayed RAIT was still statistically significantly lower than those treated with early RAIT (62.28% vs 75.53%; log-rank P < .0001) (Fig. 1B). On multivariable analysis with IPTW (see Table 3 and Supplementary e-Table S5) (37), there was a statistically significant positive association between delayed RAI use and a greater risk of persistent or recurrent disease (OR, 2.47; P < .001). As presented in Supplementary e-Table S6 (37), the association of WRAIT with persistent or recurrent disease was not mediated by ETE, multifocality, or stage N1b disease (data not shown).

Kaplan-Meier survival curves for patients treated with early radioactive iodine (RAI) (≤88 days) and delayed RAI (>88 days) following surgery. A, Kaplan-Meier survival curves for patients treated with early RAI and delayed RAI following surgery before propensity score weighting. B, Kaplan-Meier survival curves for patients treated with early RAI and delayed RAI following surgery after propensity score weighting. P value is shown for comparison by using the log-rank test.
Figure 1.

Kaplan-Meier survival curves for patients treated with early radioactive iodine (RAI) (≤88 days) and delayed RAI (>88 days) following surgery. A, Kaplan-Meier survival curves for patients treated with early RAI and delayed RAI following surgery before propensity score weighting. B, Kaplan-Meier survival curves for patients treated with early RAI and delayed RAI following surgery after propensity score weighting. P value is shown for comparison by using the log-rank test.

Table 3.

Effect of waiting time for radioactive iodine therapy on odds of persistent or recurrent disease

Analytic approachNo. of events/No. of patients at risk (%)OR (95% CI)P
Early RAIT groupDelayed RAIT group
Crude analysis97/445 (21.8)160/464 (34.5)1.89 (1.41-2.54)<.001
Multivariable analysis97/445 (21.8)160/464 (34.5)2.50 (1.67-3.75)<.001
PS analyses
 With IPTW107.1/457.7 (23.4)162.8/463.9 (35.1)2.47 (1.60-3.82)a<.001
 With matchingc63/306 (20.6)107/306 (35.0)2.64 (1.68-4.14)b<.001
Analytic approachNo. of events/No. of patients at risk (%)OR (95% CI)P
Early RAIT groupDelayed RAIT group
Crude analysis97/445 (21.8)160/464 (34.5)1.89 (1.41-2.54)<.001
Multivariable analysis97/445 (21.8)160/464 (34.5)2.50 (1.67-3.75)<.001
PS analyses
 With IPTW107.1/457.7 (23.4)162.8/463.9 (35.1)2.47 (1.60-3.82)a<.001
 With matchingc63/306 (20.6)107/306 (35.0)2.64 (1.68-4.14)b<.001

Full regression models can be found in Supplementary e-Tables S4, S5, and S7.

Bold indicates statistically significant P values.

Abbreviations: IPTW, inverse probability of treatment weighting; OR, odds ratio; PS, propensity score; RAIT, radioactive iodine therapy.

Shown is the primary analysis with an OR from the multivariable logistical regression model with inverse probability weighting according to propensity score.

Shown is the OR from a multivariable logistical regression model with matching according to propensity score.

Matching was performed based on nearest neighbor method in a 1:1 ratio without any replacement, with a caliper distance of 0.1 of the SD of the logit of the propensity score.

Table 3.

Effect of waiting time for radioactive iodine therapy on odds of persistent or recurrent disease

Analytic approachNo. of events/No. of patients at risk (%)OR (95% CI)P
Early RAIT groupDelayed RAIT group
Crude analysis97/445 (21.8)160/464 (34.5)1.89 (1.41-2.54)<.001
Multivariable analysis97/445 (21.8)160/464 (34.5)2.50 (1.67-3.75)<.001
PS analyses
 With IPTW107.1/457.7 (23.4)162.8/463.9 (35.1)2.47 (1.60-3.82)a<.001
 With matchingc63/306 (20.6)107/306 (35.0)2.64 (1.68-4.14)b<.001
Analytic approachNo. of events/No. of patients at risk (%)OR (95% CI)P
Early RAIT groupDelayed RAIT group
Crude analysis97/445 (21.8)160/464 (34.5)1.89 (1.41-2.54)<.001
Multivariable analysis97/445 (21.8)160/464 (34.5)2.50 (1.67-3.75)<.001
PS analyses
 With IPTW107.1/457.7 (23.4)162.8/463.9 (35.1)2.47 (1.60-3.82)a<.001
 With matchingc63/306 (20.6)107/306 (35.0)2.64 (1.68-4.14)b<.001

Full regression models can be found in Supplementary e-Tables S4, S5, and S7.

Bold indicates statistically significant P values.

Abbreviations: IPTW, inverse probability of treatment weighting; OR, odds ratio; PS, propensity score; RAIT, radioactive iodine therapy.

Shown is the primary analysis with an OR from the multivariable logistical regression model with inverse probability weighting according to propensity score.

Shown is the OR from a multivariable logistical regression model with matching according to propensity score.

Matching was performed based on nearest neighbor method in a 1:1 ratio without any replacement, with a caliper distance of 0.1 of the SD of the logit of the propensity score.

The robustness of our findings was assessed using post hoc sensitivity analyses. To further confirm the effects of WRAIT, we implemented PS-matching analysis. From the delayed RAIT group, 306 patients (65.9%) were matched with 306 patients who underwent early RAIT. Patients treated with and without timely RAI were balanced on 25 baseline characteristics (see Table 1). In our propensity-matched cohort, WRAIT exceeding 88 days was still associated with an increased risk of PTC persistence or recurrence (OR, 2.64; P < .001) (see Table 3 and Supplementary e-Table S7) (37). Moreover, the magnitude of the associations between RAI administration delay and the persistent or recurrent disease outcomes of interest did not substantially change when we did subgroup analyses regarding age, ATA risk group, T stage, N stage, and stimulated Tg levels before the initial RAIT (Supplementary e-Table S8) (37), and when we repeated analyses by selecting 90 days or 180 days as a cutoff at which delay is defined (Supplementary e-Table S9) (37).

Discussion

This study investigated the effect of and factors for delayed RAIT following TT for an exclusively node-positive PTC population. The results demonstrated that WRAIT exceeding 88 days was associated with a statistically significantly higher rate of persistent or recurrent PTC events. This association remained even when accounting for numerous predetermined potentially confounding patient, socioeconomic, hospital, and tumor characteristics and in the PS-weighting and multiple sensitivity analyses. Although a causal relationship between WRAIT and clinical outcomes cannot be established, our results implied that shorter WRAIT could potentially reduce the incidence of PTC persistence or recurrence. Moreover, using regression analysis, we identified several independent risk factors for delayed RAIT, which included residence in lower-income areas, reoperation before the initial RAIT, TT at a nonuniversity-affiliated hospital, multifocality, ETE, stage N1b disease, and pre–RAIT-stimulated Tg level less than 1 ng/mL.

Previous investigations into the prolonged WRAIT for patients with low-, intermediate-, or high-risk differentiated thyroid cancer differed concerning the patient, the sample size, the treatment characteristics included, the statistical approaches used, the outcome measures used, as well as the arbitrarily selected cutoff for the definition of early vs delayed RAITs (14, 18-23). As a result, while some investigations have found that an increasing WRAIT is associated with an increased probability of disease persistence or recurrence, suggesting that longer WRAIT may risk tumor regrowth (14, 20, 21), other investigations have demonstrated that a reasonable delay between thyroidectomy and RAIT does not carry a worse prognosis on patient outcomes (18, 19, 22, 23). Importantly, no studies have examined the prognostic significance of the WRAIT, specifically in patients with nodal metastases, who represent a common clinical scenario and are associated with increased PTC recurrence and mortality. To our knowledge, the present study is the first to address the optimal timing of postoperative administration of RAI in a sizable cohort of 909 patients with node-positive disease. A strength of this study is that the RPA approach improves on previous efforts by examining the effect of WRAIT as a continuous variable on patient outcomes in an objective fashion without prespecified an arbitrary cutoff time. Another strength of this study was that our sample size is significantly larger than previously published studies focusing on the effect of RAI timing on patients with low-, intermediate-, or high-risk disease. Furthermore, we used the PS-weighting method to minimize inherent bias, and we relied on the conservative analytic approach, IPTW-adjusted models, and sensitivity analyses to base our conclusions.

There has been less research on the predictors of RAI treatment delay in differentiated thyroid cancer. Similar to what others have described regarding head and neck cancers (38, 39), we observed that patients who lived in lower-income areas were more likely to have extended WRAIT, implying that socioeconomic barriers remain considerable in preventing them from receiving timely and appropriate care. Interestingly, patients who underwent TT in nonuniversity-affiliated hospitals reported longer WRAIT than those who underwent TT in university-affiliated hospitals. Indeed, many pathology reports from these nonuniversity-affiliated hospitals did not contain additional clinical factors required for initial risk stratification and postoperative RAI decision-making, such as the presence or absence of ENE, size of the largest metastatic focus within an LN, BRAF mutational status, and specific histological subtypes. In such contexts, surgical pathology specimens are typically requested from referring surgical facilities and reevaluated independently by 2 experienced pathologists at our facility. Additionally, most nonuniversity-affiliated hospitals are located farther away. Thus, obtaining and analyzing additional pathologic materials increase the WRAIT. It is also possible that head and neck surgeons at nonuniversity-affiliated hospitals are more likely to delay referral to nuclear medicine specialists for RAIT. However, we were unable to assess this based on our analyses. Repeated local surgery after the initial surgery for PTC can also affect the timeliness. In our study cohort, patients who underwent reoperation took longer to receive their first RAI administration, potentially because of the additional workup required to evaluate the disease extent and ensure that they were adequate candidates for reoperation. In addition, a postoperative stimulated Tg value of less than 1 ng/mL generally suggests long-term remission and reflects a prolonged event-free survival (40, 41). We cannot exclude that physicians may have selected patients to receive a late RAI administration owing to a postoperative stimulated Tg value of less than 1 ng/mL.

In our study, the clinical characteristics associated with increased WRAIT were ETE, multifocality, and N1b disease. The reason for this is not completely clear; however, the heterogeneity of patients referred to our facility may partially explain this finding. Some patients with adverse tumor features may have experienced more pretreatment analyses that postponed the initial RAIT. This hypothesis is supported by the fact that there was a higher rate of reoperation before the initial RAIT in patients with multifocal disease and lateral cervical LN metastases (see Supplementary e-Table S3) (37). There are also data suggesting that advanced tumor stage is associated with neck reoperation (42). Furthermore, patients with ETE, multifocality, or N1b disease more frequently received TT in a nonuniversity-affiliated hospital (see Supplementary e-Table S3) (37), which may exacerbate delays in RAIT because TT for PTC in a nonuniversity-affiliated hospital has been shown to be a determinant of increased WRAIT (see Table 2). Factors that may influence the variation between hospitals in the WRAIT should be explored in the future.

Although the presence of ETE, multifocality, and N1b disease had emerged as important risk factors for the persistence or recurrence of PTC (15, 43), it was difficult to explore whether these key clinical features played a substantial role in the relationship between WRAIT and disease persistence or recurrence in our retrospective observational study. Causal mediation analysis was applied to address this limitation. In the present study, we used WRAIT as the treatment and ETE, multifocality, or N1b disease as a mediator variable and found that the effect of increased WRAIT was not mediated by each of these tumor characteristics (see Supplementary e-Table S6) (37). This finding may reflect the true effect of WRAIT in real-world clinical practice and further confirm the benefit of reduced WRAIT in patients with node-positive PTC, regardless of whether the patients present with adverse tumor features at the time of diagnosis. Further studies are warranted to identify the main factors involved in the association between delayed RAI administration and poor prognosis.

The present study has some limitations. First, this is a nonrandomized, retrospective observational study and hence suffers from potential selection and ascertainment bias despite robust PS weighting and matching. Although we attempted to include as many potential confounders as possible, it is still possible that some amount of unmeasured confounding remains, which may have affected our results. Second, although the single-center nature of this study may provide insight into the role of RAIT delay in the setting of a highly select group of patients with node-positive PTC, it may not necessarily be generalizable to histologic types other than PTC. Larger prospective studies in more diverse populations must validate our findings for generalizability. Third, we were not able to evaluate several important variables relevant to clinical outcome and prolonged WRAIT in the analyses, such as the diagnosis-to-surgery interval and individual reasons for delayed RAI use (eg, nonadherence, missed appointments, unrelated medical comorbidities, or unknown) because of incomplete clinical data before RAIT at our facility.

In conclusion, this is the first study showing that longer WRAIT beyond 88 days after TT in N1 PTC was associated with greater disease persistence or recurrence. Residence in lower-income areas, reoperation before the initial RAIT, TT at a nonuniversity-affiliated hospital, multifocality, extrathyroidal extension, N1b disease, and pre–RAIT-stimulated Tg level less than 1 ng/mL were associated with delayed RAI administration. Identifying predictive factors for RAIT delay can assist in recognizing at-risk patients with node-positive PTC and pinpoint specific areas to reduce delay.

Disclosures

All authors declare that they do not have any commercial association that might create a conflict of interest in connection with this manuscript.

Data Availability

Some or all data sets generated during and/or analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

References

1

Sung
H
,
Ferlay
J
,
Siegel
RL
, et al.
Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries
.
CA Cancer J Clin
.
2021
;
71
(
3
):
209
249
.

2

Lim
H
,
Devesa
SS
,
Sosa
JA
,
Check
D
,
Kitahara
CM
.
Trends in thyroid cancer incidence and mortality in the United States, 1974-2013
.
JAMA
.
2017
;
317
(
13
):
1338
1348
.

3

LeClair
K
,
Bell
KJL
,
Furuya-Kanamori
L
,
Doi
SA
,
Francis
DO
,
Davies
L
.
Evaluation of gender inequity in thyroid cancer diagnosis: differences by sex in US thyroid cancer incidence compared with a meta-analysis of subclinical thyroid cancer rates at autopsy
.
JAMA Intern Med
.
2021
;
181
(
10
):
1351
1358
.

4

Wada
N
,
Suganuma
N
,
Nakayama
H
, et al.
Microscopic regional lymph node status in papillary thyroid carcinoma with and without lymphadenopathy and its relation to outcomes
.
Langenbecks Arch Surg
.
2007
;
392
(
4
):
417
422
.

5

Randolph
GW
,
Duh
QY
,
Heller
KS
, et al.
The prognostic significance of nodal metastases from papillary thyroid carcinoma can be stratified based on the size and number of metastatic lymph nodes, as well as the presence of extranodal extension
.
Thyroid
.
2012
;
22
(
11
):
1144
1152
.

6

Nixon
IJ
,
Wang
LY
,
Palmer
FL
, et al.
The impact of nodal status on outcome in older patients with papillary thyroid cancer
.
Surgery
.
2014
;
156
(
1
):
137
146
.

7

Adam
MA
,
Pura
J
,
Goffredo
P
, et al.
Presence and number of lymph node metastases are associated with compromised survival for patients younger than age 45 years with papillary thyroid cancer
.
J Clin Oncol
.
2015
;
33
(
21
):
2370
2375
.

8

Mazzaferri
EL
,
Jhiang
SM
.
Long-term impact of initial surgical and medical therapy on papillary and follicular thyroid cancer
.
Am J Med
.
1994
;
97
(
5
):
418
428
.

9

Ruel
E
,
Thomas
S
,
Dinan
M
,
Perkins
JM
,
Roman
SA
,
Sosa
JA
.
Adjuvant radioactive iodine therapy is associated with improved survival for patients with intermediate-risk papillary thyroid cancer
.
J Clin Endocrinol Metab
.
2015
;
100
(
4
):
1529
1536
.

10

Ruel
E
,
Thomas
S
,
Perkins
JM
,
Roman
SA
,
Sosa
JA
.
The impact of pathologically positive lymph nodes in the clinically negative neck: an analysis of 39,301 patients with papillary thyroid cancer
.
Ann Surg Oncol
.
2017
;
24
(
7
):
1935
1942
.

11

Sun
Y
,
Gong
J
,
Guo
B
,
Shang
J
,
Cheng
Y
,
Xu
H
.
Association of adjuvant radioactive iodine therapy with survival in node-positive papillary thyroid cancer
.
Oral Oncol
.
2018
;
87
:
152
157
.

12

Hanna
TP
,
King
WD
,
Thibodeau
S
, et al.
Mortality due to cancer treatment delay: systematic review and meta-analysis
.
BMJ
.
2020
;
371
:
m4087
.

13

Liao
DZ
,
Schlecht
NF
,
Rosenblatt
G
, et al.
Association of delayed time to treatment initiation with overall survival and recurrence among patients with head and neck squamous cell carcinoma in an underserved urban population
.
JAMA Otolaryngol Head Neck Surg
.
2019
;
145
(
11
):
1001
1009
.

14

Higashi
T
,
Nishii
R
,
Yamada
S
, et al.
Delayed initial radioactive iodine therapy resulted in poor survival in patients with metastatic differentiated thyroid carcinoma: a retrospective statistical analysis of 198 cases
.
J Nucl Med
.
2011
;
52
(
5
):
683
689
.

15

Haugen
BR
,
Alexander
EK
,
Bible
KC
, et al.
2015 American Thyroid Association management guidelines for adult patients with thyroid nodules and differentiated thyroid cancer: the American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer
.
Thyroid
.
2016
;
26
(
1
):
1
133
.

16

National Comprehensive Cancer Network Inc 2021 practice guidelines in oncology–thyroid carcinoma v.3. Accessed October 15, 2021. http://www.nccn.org/professionals/physician_gls/PDF/thyroid.pdf

17

Pacini
F
,
Fuhrer
D
,
Elisei
R
, et al.
2022 ETA consensus statement: what are the indications for post-surgical radioiodine therapy in differentiated thyroid cancer?
Eur Thyroid J
.
2022
;
11
(
1
):
e210046
.

18

Tsirona
S
,
Vlassopoulou
V
,
Tzanela
M
, et al.
Impact of early vs late postoperative radioiodine remnant ablation on final outcome in patients with low-risk well-differentiated thyroid cancer
.
Clin Endocrinol (Oxf)
.
2014
;
80
(
3
):
459
463
.

19

Scheffel
RS
,
Zanella
AB
,
Dora
JM
,
Maia
AL
.
Timing of radioactive iodine administration does not influence outcomes in patients with differentiated thyroid carcinoma
.
Thyroid
.
2016
;
26
(
11
):
1623
1629
.

20

Krajewska
J
,
Jarzab
M
,
Kukulska
A
, et al.
Postoperative radioiodine treatment within 9 months from diagnosis significantly reduces the risk of relapse in low-risk differentiated thyroid carcinoma
.
Nucl Med Mol Imaging
.
2019
;
53
(
5
):
320
327
.

21

Li
H
,
Zhang
YQ
,
Wang
C
,
Zhang
X
,
Li
X
,
Lin
YS
.
Delayed initial radioiodine therapy related to incomplete response in low- to intermediate-risk differentiated thyroid cancer
.
Clin Endocrinol (Oxf)
.
2018
;
88
(
4
):
601
606
.

22

Matrone
A
,
Gambale
C
,
Torregrossa
L
, et al.
Delayed 131-I first treatment after surgery has no impact on the median term outcome of patients with intermediate risk differentiated thyroid cancer
.
Endocr Pract
.
2020
;
26
(
1
):
58
71
.

23

Kim
M
,
Han
M
,
Jeon
MJ
, et al.
Impact of delayed radioiodine therapy in intermediate-/high-risk papillary thyroid carcinoma
.
Clin Endocrinol (Oxf)
.
2019
;
91
(
3
):
449
455
.

24

Charlson
ME
,
Pompei
P
,
Ales
KL
,
MacKenzie
CR
.
A new method of classifying prognostic comorbidity in longitudinal studies: development and validation
.
J Chronic Dis
.
1987
;
40
(
5
):
373
383
.

25

Papaleontiou
M
,
Hughes
DT
,
Guo
C
,
Banerjee
M
,
Haymart
MR
.
Population-based assessment of complications following surgery for thyroid cancer
.
J Clin Endocrinol Metab
.
2017
;
102
(
7
):
2543
2551
.

26

Zhan
G
,
Ling
L
,
Ying
L
,
Li
YP
,
Jin
W
,
Li
W
.
Regional economic levels and adverse events linked to intrauterine devices
.
J Evid Based Med
.
2011
;
4
(
1
):
8
14
.

27

Liu
H
,
Song
B
,
Jin
J
, et al.
Length of stay, hospital costs and mortality associated with comorbidity according to the Charlson Comorbidity Index in immobile patients after ischemic stroke in China: a national study
.
Int J Health Policy Manag
.
2021
;
11
(
9
):
1780
1787
.

28

Edge
SB
,
Compton
CC
.
The American Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM
.
Ann Surg Oncol
.
2010
;
17
(
6
):
1471
1474
.

29

Perrier
ND
,
Brierley
JD
,
Tuttle
RM
.
Differentiated and anaplastic thyroid carcinoma: major changes in the American Joint Committee on Cancer eighth edition cancer staging manual
.
CA Cancer J Clin
.
2018
;
68
(
1
):
55
63
.

30

van Velsen
EFS
,
Stegenga
MT
,
van Kemenade
FJ
, et al.
Evaluation of the 2015 ATA guidelines in patients with distant metastatic differentiated thyroid cancer
.
J Clin Endocrinol Metab
.
2020
;
105
(
3
):
e457
e465
.

31

Trimboli
P
,
Zilioli
V
,
Imperiali
M
,
Ceriani
L
,
Giovanella
L
.
High-sensitive basal serum thyroglobulin 6-12 months after thyroid ablation is strongly associated with early response to therapy and event-free survival in patients with low-to-intermediate risk differentiated thyroid carcinomas
.
Eur J Endocrinol
.
2017
;
176
(
5
):
497
504
.

32

Kass
GV
.
An exploratory technique for investigating large quantities of categorical data
.
J R Stat Soc
.
1980
;
29
(
2
):
119
127
.

33

Camy
F
,
Karpathiou
G
,
Dumollard
JM
, et al.
Brain metastasis PD-L1 and CD8 expression is dependent on primary tumor type and its PD-L1 and CD8 status
.
J Immunother Cancer
.
2020
;
8
(
2
):
e000597
.

34

Koch
E
,
Bogado
M
,
Araya
F
, et al.
Impact of parity on anthropometric measures of obesity controlling by multiple confounders: a cross-sectional study in Chilean women
.
J Epidemiol Community Health
.
2008
;
62
(
5
):
461
470
.

35

Cole
SR
,
Hernán
MA
.
Constructing inverse probability weights for marginal structural models
.
Am J Epidemiol
.
2008
;
168
(
6
):
656
664
.

36

Zhang
Z
,
Zheng
C
,
Kim
C
,
Van Poucke
S
,
Lin
S
,
Lan
P
.
Causal mediation analysis in the context of clinical research
.
Ann Transl Med
.
2016
;
4
(
21
):
425
.

37

Sun
YG
,
Sun
QL
,
Tian
JY
,
He
XC
.
Supplementary data for “Impact of waiting time for radioactive iodine therapy on outcome in N1 stage papillary thyroid cancer.” Figshare 2023. Deposited April 23, 2023. https://doi.org/10.6084/m9.figshare.22304140.v5

38

van Harten
MC
,
Hoebers
FJP
,
Kross
KW
,
van Werkhoven
ED
,
van den Brekel
MWM
,
van Dijk
BAC
.
Determinants of treatment waiting times for head and neck cancer in the Netherlands and their relation to survival
.
Oral Oncol
.
2015
;
51
(
3
):
272
278
.

39

Murphy
CT
,
Galloway
TJ
,
Handorf
EA
, et al.
Survival impact of increasing time to treatment initiation for patients with head and neck cancer in the United States
.
J Clin Oncol
.
2016
;
34
(
2
):
169
178
.

40

Rosario
PW
,
Xavier
AC
,
Calsolari
MR
.
Value of postoperative thyroglobulin and ultrasonography for the indication of ablation and 131I activity in patients with thyroid cancer and low risk of recurrence
.
Thyroid
.
2011
;
21
(
1
):
49
53
.

41

Pacini
F
,
Molinaro
E
,
Lippi
F
, et al.
Prediction of disease status by recombinant human TSH-stimulated serum Tg in the postsurgical follow-up of differentiated thyroid carcinoma
.
J Clin Endocrinol Metab
.
2001
;
86
(
12
):
5686
5690
.

42

Semrad
TJ
,
Keegan
THM
,
Semrad
A
,
Brunson
A
,
Farwell
DG
.
Predictors of neck reoperation and mortality after initial total thyroidectomy for differentiated thyroid cancer
.
Thyroid
.
2018
;
28
(
9
):
1143
1152
.

43

Kim
H
,
Kwon
H
,
Moon
BI
.
Association of multifocality with prognosis of papillary thyroid carcinoma: a systematic review and meta-analysis
.
JAMA Otolaryngol Head Neck Surg
.
2021
;
147
(
10
):
847
854
.

Abbreviations

     
  • APCDI

    annual per capita disposable income

  •  
  • ATA

    American Thyroid Association

  •  
  • CT

    computed tomography

  •  
  • DFS

    disease-free survival

  •  
  • ENE

    extranodal extension

  •  
  • ETE

    extrathyroidal extension

  •  
  • IPTW

    inverse probability of treatment weighting

  •  
  • IQR

    interquartile range

  •  
  • LN

    lymph node

  •  
  • MRI

    magnetic resonance imaging

  •  
  • OR

    odds ratio

  •  
  • PET

    positron emission tomography

  •  
  • PS

    propensity score

  •  
  • PTC

    papillary thyroid cancer

  •  
  • RAI

    radioactive iodine

  •  
  • RAIT

    radioactive iodine therapy

  •  
  • RPA

    recursive partitioning analysis

  •  
  • SES

    socioeconomic status

  •  
  • Tg

    thyroglobulin

  •  
  • TgAb

    thyroglobulin antibody

  •  
  • TSH

    thyrotropin

  •  
  • TT

    total thyroidectomy

  •  
  • VIF

    variance inflation factor

  •  
  • WRAIT

    waiting time for radioactive iodine therapy

Author notes

Y.S., Q.S., and J.T. contributed equally to this work.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/pages/standard-publication-reuse-rights)