Cigarette smoking and thyroid cancer risk: a cohort study.

BACKGROUND: We evaluated the association between smoking status and thyroid cancer risk and whether this association is mediated by body mass index (BMI) and thyroid-stimulating hormone (TSH).
METHODS: We performed a cohort study of 96,855 Korean adults who were followed annually or biennially for a median of 5.9 years.
RESULTS: During 511,052.9 person-years of follow-up, 1,250 participants developed thyroid cancer. In men, we observed a dose-dependent inverse association between current smoking, pack-years, and thyroid cancer. After adjustment for confounders, adjusted hazard ratios (95% confidence intervals) for thyroid cancer comparing current and former smokers to never smokers were 0.58 (0.45-0.75) and 0.93 (0.73-1.18), respectively. After further adjustment for BMI and TSH as potential mediators, this association was slightly attenuated, but remained significant. For women, current smokers tended to have a lower risk of thyroid cancer, but this association did not reach statistical significance.
CONCLUSIONS: In this cohort study, current smoking was associated with a decreased risk of incident thyroid cancer in men but not in women and this association was observed even after adjusting for TSH and BMI levels as potential mediators. Further mechanistic studies are needed to elucidate the possible effect of smoking on the pathogenesis of thyroid cancer development.

Introduction

The incidence of thyroid cancer has increased substantially over the past few decades in most areas of the world.[1,2] This increase can be attributed to enhanced detection of small, clinically insignificant thyroid tumors through increasingly widespread use of sensitive screening and diagnostic procedures.[3,4] Besides the contribution of over-diagnosis, evidence also indicates that a true increase is also occurring.[4] Some factors such as ionizing radiation exposure, female sex, non-cancerous benign thyroid conditions, and obesity appear to increase the risk of thyroid cancer, but the causes of thyroid cancer are not fully understood.[4]

Smoking is the single greatest avoidable risk factor for various types of cancer, including acute myeloid leukemia, lung, esophagus, larynx, mouth, throat, kidney, bladder, stomach, cervix, and colorectal cancers.[5] The association between cigarette smoking and thyroid cancer has been evaluated in observational studies, but the results remain controversial. Some case-control[6-16] and cohort studies showed an inverse relationship between smoking and thyroid cancer,[17-20] while other studies found no association.[21-23] Recently, in a meta-analysis of 25 case-control studies and six cohort studies, smoking was associated with a lower risk of thyroid cancer in current smokers, but this association was observed only in case-control studies, requiring a well-planned cohort study with a larger sample size.[24] Additionally, in previous studies, the mechanism underlying the inverse association between smoking and thyroid cancer has been proposed to be mediated by lower BMI and lower thyroid-stimulating hormone (TSH) in smokers.[25] Indeed, smoking has been reported to be associated with reduced body weight, leading to a lower BMI[25,26] Smoking appears to be associated with adverse effects on the thyroid, resulting in changes in circulating hormone concentrations, while higher TSH values are associated with a higher frequency and more advanced stages of thyroid cancer.[27,28]Until now, no cohort studies have evaluated the effect of smoking on the incidence of thyroid cancer while taking into account the effects of BMI and TSH.

Therefore, we performed a longitudinal cohort study to evaluate the association between smoking status and the development of thyroid cancer in young and middle aged Korean men and women and examined whether this association was mediated by TSH levels and BMI.

Methods

Study population

The Kangbuk Samsung Health Study is a retrospective cohort study of men and women 18 years of age or older who underwent a comprehensive annual or biennial health examination at the clinics of the Kangbuk Samsung Hospital Total Healthcare Screening Center in Seoul and Suwon, South Korea, from 2002 to present.[29,30] Most examinees (over 80%) are employees of various companies or local governmental organizations and their spouses. In South Korea, the Industrial Safety and Health Law requires annual or biennial health screening examinations of all employees, free of charge. The rest of the examinees voluntarily purchased health screening exams at the healthcare center.

The present study used de-identified retrospective data routinely collected during the health screening process of all study participants with at least one follow-up visit who underwent a comprehensive health examination at clinics between 2006 and 2008 and were followed annually or biennially until December 2014 (N = 98,553). A total of 1,698 participants were excluded due to missing data on smoking status, TSH levels, or BMI (n = 264), or a history of cancer (n = 1,437). Because some individuals met more than one exclusion criterion, the total number of patients eligible for the study was 96,855. This study was approved by the Institutional Review Board of Kangbuk Samsung Hospital (IRB No. KBSMC 2017-10-015), and the requirement for informed consent was waived because we used de-identified retrospective data routinely collected during the health screening process.

Data collection

Baseline and follow-up examinations were conducted at Kangbuk Samsung Hospital Health Screening Center clinics in Seoul and Suwon. All participants were asked to complete a standardized, self-administered questionnaire about their socioeconomic status, smoking status, alcohol intake, regular exercise habits, medical history and medication use. Current alcohol use was assessed regarding as the frequency of alcohol drinking per week and amount of alcohol consumed per drinking day.[31]Alcohol consumption was categorized into none, <20 g of ethanol/day, and ≥20 g of ethanol/day for men and none, <10 g/day, and ≥10 g of ethanol/day for women.[32] Regular exercise habits were evaluated the weekly frequency of moderate- or vigorous-intensity physical activity. Smoking status was categorized into never, former, and current smokers. The pack-years was calculated as the number of cigarettes smoked per day multiplied by the number of years of previous or current smoking.

Weight, height, and sitting blood pressure were assessed by trained nurses. Weight was measured in light clothing and no shoes to the nearest 0.1 kg and height was measured to the nearest 0.1 cm using a digital scale. Body mass index (BMI) was calculated as weight in kg divided square of height and obesity was defined as BMI ≥25 kg/m2 according to Asian-specific criteria.[33,34] Hypertension was defined as a systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg, or current use of antihypertensive medication.

Incident thyroid cancer cases

Thyroid cancer cases were primarily ascertained based on self-report of physician-diagnosed thyroid cancer by standardized, self-administered questionnaires at the health screening visits. Thyroid cancer cases were defined as participants who reported physician-diagnosed thyroid cancer during follow-up. In a subsample of 57,741 participants who underwent health examination and provided informed consent for linkage to the national registries between 2011 and 2014, we evaluated the accuracy of self-reported physician-diagnosed thyroid cancer. When self-reported thyroid cancer was compared with confirmed thyroid cancer through national cancer registry data among 57,741 participants, the sensitivity and specificity for self-report of physician-diagnosed thyroid cancer were 98.1 % and 99.8 %, respectively. Specific histologic types of thyroid cancer included papillary (97.4 %), follicular (1.4 %), carcinoid (0.1 %), medullary (0.1 %), mixed medullary-papillary carcinoma (0.9 %), and unspecified cancer (0.2 %).

Laboratory analyses

Blood samples were taken from the antecubital vein after at least a 10-h fast. Measurement of serum uric acid, glucose, total cholesterol, triglyceride low-density lipoprotein cholesterol (LDL-C), alanine aminotransferase, aspartate aminotransferase levels, γ-glutamyl transferase, insulin, and high sensitivity C-reactive protein (hsCRP) has been described in detail elsewhere.[35] Insulin resistance was assessed with HOMA-IR according to the following equation: (fasting blood insulin (μU/ml) × fasting blood glucose (mmol/l))/22.5.[36] To assess thyroid function, serum TSH was measured by radioimmunoassay using a commercial kit (RIA-gnost VRhTSH, Schering-Cis Bio International, Gif-sur-Yvette, France), with a lower detection limit of 0.025 μIU/ml. The normal range was 0.25–5.0 μIU/ml, and the intra- and inter-assay coefficients of variation for quality control specimens were 1.2–5.7 and 2.4–5.4%, respectively, for TSH. Diabetes was defined as a fasting serum glucose ≥126 mg/dl, or current use of insulin or antidiabetic medications. The Laboratory Medicine Department at Kangbuk Samsung Hospital, Korea has been accredited by the Korean Society of Laboratory Medicine (KSLM) and the Korean Association of Quality Assurance for Clinical Laboratories (KAQACL). The laboratory participates in the College of American Pathologists (CAP) survey proficiency testing.

Statistical analysis

Subject characteristics were presented according to smoking status (never, former and current smokers) and were compared according to the smoking status using a linear regression test for continuous variables or a χ2 test for categorical variables. To test for linear trends, category numbers were used as continuous variables in regression models. All analyses were performed separately in men and women, as the association between smoking and thyroid cancer differs significantly by sex.

The study endpoint was the development of thyroid cancer. Follow-up for each participant extended from the baseline exam until either the development of thyroid cancer or the last health exam conducted prior to 31 December 2014. Incidence rate was calculated as the number of incident cases divided by person-years of follow-up. Since we knew that thyroid cancer had occurred at some point between the two visits but did not know the precise timing of development, we used a parametric proportional hazards model to take this type of interval censoring into account (stpm command in Stata).[37] In these models, the baseline hazard function was parameterized with restricted cubic splines in log time with four degrees of freedom. The adjusted hazard ratio (aHR) with a 95% confidence interval (CI) for incident thyroid cancer was estimated separately in men and women. The model was initially adjusted for age and then further adjusted for center (Seoul and Suwon), year of screening exam, alcohol intake (≥20 vs. <20 g of ethanol/day), regular exercise (≥3 vs. <3 times/week), and educational level (≥college graduate vs. college graduate or below) in model 1. To evaluate whether the association between smoking and incident thyroid cancer was mediated by BMI and TSH on a priori grounds, we included these variables in a multivariable-adjusted model (model 2). The following criteria for mediation analysis were used: (1) the predictor of interest (smoking status) was significantly related to the mediator (BMI and TSH), (2) the mediator (BMI and TSH) was significantly related to the outcome (thyroid cancer), and (3) the addition of the mediator (BMI and TSH) to the model both attenuated the coefficient of the predictor and had a statistically significant mediation effect. Moreover, to evaluate the impact of change of smoking status, pack- years and covariates over follow-up, we conducted additional analyses introducing smoking status and other covariates as a time-varying covariate in the models.

We assessed the proportional hazards assumption by examining graphs of estimated log (-log) survival. The number of categories was used as a continuous variable and tested on each model to determine linear trends of incidence.

In addition, we performed stratified analysis in pre-specified subgroups defined by age ( ≥ 50 vs. < 50 years), alcohol intake (none, < 20 vs. ≥ 20 g of ethanol/day for men; none, < 10 vs. ≥ 10 g for women), vigorous exercise (none, < 3 vs. ≥ 3 times/week), BMI ( ≥ 25 vs. < 25 kg/m2), HOMA-IR ( ≥ 2.5 vs. < 2.5), and hsCRP ( ≥ 1.0 vs. < 1.0 mg/L).

We conducted a series of sensitivity analyses to test the robustness of our primary results. First, Sensitivity analysis was conducted after implementing a 3-year wash-out period, during which any thyroid cancer cases identified were excluded. This approach was taken to assess if the original results may have been influenced by reverse causality in which thyroid cancer was present, but undiagnosed. Second, we performed the sensitivity analysis among 57,741 participants with available information on the cancer registry data. For this analysis, thyroid cancer cases were defined as participants who were diagnosed with a first primary thyroid malignant neoplasm during follow-up. Third, we performed the sensitivity analysis among 57,741 participants with available information on the cancer registry data using only papillary thyroid carcinoma as an endpoint because papillary thyroid carcinomas and medullary thyroid carcinomas have been considered different in terms of their incident rates, cell origins, risk factors and histopathological features.[38]

All statistical analyses were performed using STATA version 13.0 software (Stata Corp LP, College Station, TX).

Results

Out of 96,855 subjects, 57.9% were male. The mean (SD) age of male and female participants were 40.4 (7.7) and 39.4 (7.7), respectively. Baseline characteristics of the subjects are shown by sex (Tables 1 and 2). Among male participants, current smokers were more likely to be obese, to drink alcohol, and have higher levels of LDL-C, triglycerides, and hsCRP and lower levels of systolic BP than never or ex-smokers (Table 1). Among female participants, the mean (SD) age and BMI was 39.4 (7.7) and 22.0 (2.9) kg/m2 respectively. Female current smokers were more likely to drink alcohol, to have higher levels of HDL-C and triglycerides, and to have lower levels of systolic BP, glucose, and LDL-C (Table 2).

Table 1

Baseline characteristics of study participants by smoking status among men

Characteristics		Smoking status			P for trend
	Overall	Never smoker	Ex-smoker	Current smoker
Number	56,060	17,005	16,392	22,663	<0.001
Age (years)a	40.4 (7.7)	39.8 (8.1)	42.2 (8.0)	39.7 (6.8)	<0.001
BMI (kg/m2)	24.5 (2.8)	24.2 (2.8)	24.6 (2.7)	24.6 (2.9)	<0.001
Obesity (%)	40.0	36.7	41.1	41.6	<0.001
Alcohol intake (%)b	19.7	10.3	19.3	27.0	<0.001
Vigorous exercise (%)c	20.4	21.2	27.0	15.1	<0.001
High education level (%)d	97.9	97.1	96.8	99.3	0.212
Diabetes (%)	4.7	3.6	5.4	5.0	<0.001
Hypertension (%)	20.4	19.6	24.2	18.4	<0.001
Systolic BP (mmHg)a	118.5 (12.9)	118.7 (13.0)	119.5 (13.1)	117.6 (12.6)	<0.001
Diastolic BP (mmHg)a	76.8 (8.8)	76.7 (8.8)	77.5 (8.8)	76.4 (8.8)	<0.001
Glucose (mg/dL)a	97.8 (15.9)	97.1 (15.0)	98.6 (15.4)	97.7 (17.0)	0.001
Total cholesterol (mg/dL)a	195.2 (33.0)	193.1 (32.3)	196.5 (33.1)	195.8 (33.4)	<0.001
LDL-C (mg/dL)a	115.7 (29.1)	114.7 (28.4)	116.0 (29.0)	116.3 (29.7)	<0.001
HDL-C (mg/dL)a	49.4 (10.5)	49.8 (10.5)	50.1 (10.6)	48.7 (10.3)	<0.001
Triglycerides (mg/dL)e	124 (88–177)	111 (80–157)	121 (87–172)	137 (97–195)	<0.001
HOMA-IRe	1.93 (1.47–2.55)	1.92 (1.47–2.53)	1.97 (1.51–2.58)	1.91 (1.44–2.55)	0.054
hsCRP (mg/L)e	0.6 (0.3–1.1)	0.5 (0.3–1.0)	0.5 (0.3–1.1)	0.6 (0.3–1.2)	<0.001

Data are presented as the mean (standard deviation), median (interquartile range), or percentage.

BMI body mass index, BP blood pressure, LDL-C low-density lipoprotein-cholesterol, HDL-C high-density lipoprotein-cholesterol, HOMA-IR homeostasis model assessment of insulin resistance, hsCRP high sensitivity C-reactive protein.

aMean (standard deviation)

b≥20 g of ethanol per day

c≥3 times per week

d≥college graduate.

eMedian (interquartile range)

Table 2

Baseline characteristics of study participants by smoking status among women

Characteristics		Smoking status			P for trend
	Overall	Never smoker	Ex-smoker	Current smoker
Number	40,795	38,712	1278	805	<0.001
Age (years)a	39.4 (7.7)	39.4 (7.7)	38.4 (8.1)	39.0 (7.9)	<0.001
BMI (kg/m2)	22.0 (2.9)	22.0 (2.9)	22.0 (3.2)	22.0 (3.2)	0.666
Obesity (%)	13.6	13.6	15.2	14.7	0.103
Alcohol intake (%)b	1.3	1.0	3.3	13.6	<0.001
Vigorous exercise (%)c	20.5	20.6	17.3	20.3	0.093
High education level (%)d	98.2	98.0	100.0	100.0	<0.001
Diabetes (%)	1.7	1.6	2.6	1.5	0.287
Hypertension (%)	6.8	6.8	6.0	8.1	0.572
Systolic BP (mmHg)a	107.0 (12.7)	107.1 (12.7)	105.8 (12.3)	105.5 (13.2)	<0.001
Diastolic BP (mmHg)a	68.4 (8.9)	68.5 (8.9)	67.8 (8.8)	67.8 (8.9)	0.002
Glucose (mg/dL)a	92.1 (11.7)	92.2 (11.7)	91.7 (13.3)	91.4 (11.6)	0.015
Total cholesterol (mg/dL)a	184.5 (32.1)	184.5 (32.1)	184.4 (32.8)	183.7 (31.0)	0.493
LDL-C (mg/dL)a	101.7 (27.5)	101.8 (27.5)	100.0 (27.6)	99.0 (26.8)	<0.001
HDL-C (mg/dL)a	58.9 (12.9)	58.8 (12.8)	60.4 (13.6)	61.1 (14.4)	<0.001
Triglycerides (mg/dL)e	77 (58–106)	77 (58–106)	75.5 (58–108)	83 (63–116)	<0.001
HOMA-IRe	1.73 (1.26–2.26)	1.73 (1.26–2.26)	1.68 (1.21–2.23)	1.71 (1.37–2.23)	0.941
hsCRP (mg/L)e	0.3 (0.1–0.7)	0.3 (0.1–0.7)	0.3 (0.1–0.8)	0.3 (0.1–0.7	0.067

Data are presented as the mean (standard deviation), median (interquartile range), or percentage

BMI body mass index, BP blood pressure, LDL-C low-density lipoprotein-cholesterol, HDL-C high-density lipoprotein-cholesterol, HOMA-IR homeostasis model assessment of insulin resistance, hsCRP high sensitivity C-reactive protein.

aMean (standard deviation)

b≥10 g of ethanol per day

c≥3 times per week

d≥College graduate

eMedian (interquartile range)

Table 3 shows the development of thyroid cancer by smoking status. During 511,052.9 person-years of follow-up, 1,250 participants developed thyroid cancer (incidence rate, 3.4 per 1,000 person-years). The median follow-up period for participants was 5.9 years (interquartile range 4.0–6.8). For men, smoking status was significantly associated with a decreased risk of thyroid cancer. In the model adjusted for age, center, year of screening exam, alcohol intake, regular exercise, and educational level, the multivariable-adjusted HRs (95% CIs) for incident thyroid cancer comparing current and former smokers to never smokers were 0.58 (0.45–0.75) and 0.93 (0.73–1.18), respectively. To explore whether the association between smoking status and development of thyroid cancer was mediated by BMI and TSH, we performed additional analyses adjusting for BMI and TSH. The association between smoking and incident thyroid cancer was slightly attenuated, but remained statistically significant in men (p for trend <0.001). In women, current smokers were less likely to develop thyroid cancer, but this association did not reach statistical significance; the multivariable-adjusted HRs (95% CIs) for incident thyroid cancer comparing current and former smokers to never smokers were 0.64 (0.33–1.23) and 0.72 (0.46–1.14), respectively. When the association between smoking status and incident thyroid cancer was evaluated after introducing smoking status and confounders during follow-up as time-varying covariates, the inverse association between current smoking and incident thyroid cancer was still observed in time-dependent models in men but not in women.

Table 3

Development of thyroid cancer by smoking status

Smoking status	Person-years	Incident Cases	Incidence density (per 1,000 person-years)	Age adjusted HR (95% CI)	Multivariable-adjusted HRa (95% CI)		HR (95% CI)b in model using time-dependent variables
					Model 1	Model 2
Men
Never smoker	91,699.7	146	1.6	1.00 (reference)	1.00 (reference)	1.00 (reference)	1.00 (reference)
Ex-smoker	86,462.7	132	1.5	0.94 (0.74–1.19)	0.93 (0.73–1.18)	0.90 (0.71–1.14)	0.95 (0.74–1.21)
Current smoker	119,971.1	112	0.9	0.59 (0.46–0.75)	0.58 (0.45–0.75)	0.56 (0.43–0.72)	0.60 (0.46–0.78)
P for trend				<0.001	<0.001	<0.001	<0.001
Women
Never smoker	202,568.8	832	4.1	1.00 (reference)	1.00 (reference)	1.00 (reference)	1.00 (reference)
Ex-smoker	6542.4	19	2.9	0.73 (0.46–1.15)	0.72 (0.46–1.14)	0.72 (0.46–1.14)	1.01 (0.64–1.62)
Current smoker	3808.3	9	2.4	0.60 (0.31–1.15)	0.64 (0.33–1.23)	0.63 (0.33–1.23)	0.85 (0.60–1.19)
P for trend				0.043	0.060	0.059	0.370

BMI body mass index, HR hazard ratio, CI confidence interval.

aEstimated from parametric Cox models. Multivariable model 1 was adjusted for age, center, year of screening exam, alcohol intake, regular exercise, educational level; model 2: model 1 plus adjustment for BMI and TSH

bEstimated from parametric proportional hazard models with smoking status, alcohol intake, and regular exercise as time-dependent variables and baseline age, center, year of screening exam and education level as time-fixed variables.

We also examined the association between pack-years and the risk of thyroid cancer (Table 4). For men, there was a clear dose-response relationship between smoking amount and thyroid cancer incidence (p for trend <0.001). In model 1, the multivariable-adjusted HRs (95% CIs) for incident thyroid cancer comparing <10, 10–19.9, and ≥20 pack-years with zero pack-years were 0.60 (0.41–0.88), 0.73 (0.54–0.99), and 0.44 (0.29–0.65), respectively. For women, the multivariable-adjusted HRs (95% CIs) for incident thyroid cancer comparing <5, 5–9.9, and ≥10 pack-years with zero pack-years were 0.33 (0.13–0.88), 0.91 (0.30–2.04), and 1.39 (0.57–3.38), respectively. In mediation analysis after adjustment for BMI and TSH, these associations were similarly observed. When pack-years and confounders during follow-up were updated as time-varying covariates, the inverse association between pack-years and the risk of thyroid cancer was still observed in time-dependent models in men. In sensitivity analysis after excluding subjects who developed thyroid cancer within the first 3 years of follow-up, the associations of smoking status and pack-years with incident thyroid cancer were essentially unchanged (Supplementary Tables 3 and 4).

Table 4

Development of thyroid cancer by pack-years

Pack-years	Person-years	Incident cases	Incidence density (per 1000 person-years)	Age-adjusted HR (95% CI)	Multivariable-adjusted HRa (95% CI)		HR (95% CI)b in model using time-dependent variables
					Model 1	Model 2
Men
0	91,699.7	146	1.6	1.00 (reference)	1.00 (reference)	1.00 (reference)	1.00 (reference)
<10	36,712.5	34	0.9	0.61 (0.42–0.89)	0.60 (0.41–0.88)	0.59 (0.40–0.86)	0.64 (0.48–0.85)
10–19.9	50,912.8	60	1.2	0.76 (0.56–1.02)	0.73 (0.54–0.99)	0.69 (0.51–0.94)	0.79 (0.59–1.04)
≥20	42,056.5	32	0.8	0.45 (0.31–0.67)	0.44 (0.29–0.65)	0.40 (0.27–0.60)	0.42 (0.27–0.62)
P for trend				<0.001	<0.001	<0.001	<0.001
Women
0	202,568.8	832	4.1	1.00 (reference)	1.00 (reference)	1.00 (reference)	1.00 (reference)
<5	3067.8	4	1.3	0.33 (0.12–0.88)	0.33 (0.13–0.90)	0.33 (0.12–0.88)	0.83 (0.59–1.15)
5–9.9	1331.7	4	3.0	0.76 (0.29–2.04)	0.91 (0.30–2.17)	0.81 (0.30–2.18)	0.49 (0.15–1.53)
≥10	905.1	5	5.5	1.32 (0.55–3.17)	1.39 (0.57–3.38)	1.40 (0.57–3.40)	1.56 (0.64–3.79)
P for trend				0.410	0.504	0.500	0.421

BMI body mass index, HR hazard ratio, CI confidence interval.

aEstimated from parametric Cox models. Multivariable model 1 was adjusted for age, center, year of screening exam, alcohol intake, regular exercise, educational level; model 2: model 1 plus adjustment for BMI and TSH

bEstimated from parametric proportional hazard models with pack-years, alcohol intake, and regular exercise as time-dependent variables and baseline age, center, year of screening exam, and education level as time-fixed variables.

We performed additional analysis among 57,741 participants with available information on the cancer registry data (Supplementary Tables 3 and 4). In analysis using thyroid cancer based on the national cancer registry data among 57,741 participants only who gave informed consent for the linkage, we found an inverse association between current smoking and thyroid carcinoma and dose-dependent relationship between pack-years and thyroid carcinoma in men. When study endpoint was defined as papillary thyroid carcinoma, these associations were essentially unchanged. For women, the inverse association between current smoking and thyroid carcinoma was observed but did not reach statistical significance.

In subgroup analyses (Supplementary Tables 1 and 2), the association between smoking status and incident thyroid cancer was similar across pre-specified subgroups with no significant interactions with age ( < 50 vs. ≥ 50 years), alcohol intake (none, < 20 vs. ≥ 20 g of ethanol/day for men; none, < 10 vs. ≥ 10 g for women), vigorous exercise (none, < 3 vs. ≥ 3 times/week), BMI ( < 25 kg/m2 vs. ≥ 25 kg/m2), HOMA-IR ( < 2.5 vs. ≥ 2.5), or hsCRP ( < 1.0 mg/L vs. ≥ 1.0 mg/L).

Discussion

In this large cohort of 96,855 young and middle-aged population, current smoking at baseline was significantly associated with a decreased risk of incident thyroid cancer in men. We also observed a non-statistically significant inverse association between current smoking and incidence of thyroid cancer in women. A dose-response relationship was also found between pack-years and incident thyroid cancer in men. These associations persisted after adjustment for potential confounders, indicating that smoking is an independent risk factor for incident thyroid cancer. When we evaluated whether the association between smoking status and development of thyroid cancer was mediated by BMI and TSH, which are proposed mediators of thyroid cancer development in previous studies, the inverse association was slightly attenuated but still observed. Thus, levels of TSH and BMI by smoking could not fully explain the association between current smoking and thyroid cancer.

A pooled analysis of case-control studies of thyroid cancer in 2003 found a reduced risk of thyroid cancer in ever smokers.[11] The relationship was more pronounced in current smokers than former smokers, with odds ratios (95% CIs) of 0.6 (0.6–0.7) in current smokers and 0.9 (0.8–1.1) in former smokers However, a significant association between current smokers and decreased risk of thyroid cancer was observed in women but not in men. In a prospective cohort study of U.S. radiologic technologists followed from 1983 through 2006, current smoking was inversely associated with thyroid cancer risk in women (HR 0.54, 95% CI 0.35–0.82 in women and HR 0.31, 95% CI 0.09–1.04 in men).[17] In another prospective cohort study of postmenopausal women enrolled in the Women’s Health Initiative with a follow-up of 12.7 years, current smokers had a reduced risk for papillary thyroid cancer (HR 0.34, 95% CI 0.15–0.78) compared with never smokers.[18] However, in the same study, because the number of current smokers among cases was small, no associations or trends were seen for amount smoked. In a nine-year prospective study carried out on Koreans in the national insurance program, the inverse association between current smoking and thyroid cancer risk was observed in men (HR 0.6, 95% CI 0.5–1.0), but did not reach statistical significance due to relatively small number of thyroid cancer case.[22] A pooled analysis of five prospective studies in the United States demonstrated that current smokers had a reduced risk of thyroid cancer (HR 0.68, 95% CI 0.55–0.85) compared to never smokers, but this association was not significant in men.[19] In the same meta-analysis, greater smoking intensity, duration, and pack-years were associated with further risk reduction. The inconsistent findings across studies could be attributable to differences in study design, sample size, study population (age, ethnicity, and sex composition), control for potential confounding factors, measures of smoking status, smoking prevalence or incidence of thyroid cancer. Most studies reported that the inverse association between current smoking and thyroid cancer was statistically significant in women[17,18] but not in men due to the relatively lower incidence in men than in women.[22] The present large-scale cohort study demonstrated that current smoking decreased thyroid cancer risk among men in a dose-response manner. The effect of smoking on thyroid cancer might confounded by BMI, and smokers may have lower BMI.[25] However, the association was still observed after adjustment for measured BMI and other confounders. Furthermore, this association was observed even when changes in smoking status and confounders during follow-up were treated as time-varying covariates. For women, we found non-statistically significant inverse association between current smoking and incidence of thyroid cancer. Even though 42 % of the study population was female, the prevalence of current smoking was much lower in females than in males (2 vs. 40%), resulting in lack of power to detect an association between current smoking and thyroid cancer risk in women.

The mechanisms that link smoking and thyroid cancer are not fully understood, but several possibilities have been proposed. First, higher TSH values are associated with a higher frequency and more advanced stage of thyroid cancer,[27] while smoking seems to induce changes in thyroid function tests, such as decreases in TSH and increases in thyroid hormones. Thus, previous studies explained that lower TSH and lower BMI in smokers might influence the risk and protective effect of tobacco on thyroid cancer, which might be mediated by decreasing the body weight and TSH levels of the smoker.[8,19,25,39] However, in the present study, the inverse association between smoking and thyroid cancer was still observed even after adjustment for TSH and BMI. Another possible biological mechanism underlying the association between current smoking and thyroid cancer is the anti-estrogenic effect.[40] Estrogen has been suggested as an associated etiology of thyroid cancer, to explain the higher incidence of this cancer in females than in males.[7,11] Vladan et al. also reported that smoking was negatively associated with thyroid cancer in men, due to its ability to decrease TSH secretion and its innate anti-estrogenic effect.[15] Additionally, nicotine might be involved in the protective effects of smoking. The activation of nicotine receptors on immune cells reduced the risk of Hashimoto’s thyroiditis by shifting the autoimmune profile from pathogenic Th1 and Th17 responses to protective Th2 responses.[25] The tobacco alkaloid anatabine, which has a structure similar to nicotine, also has been shown to decrease the incidence and severity of experimental autoimmune thyroiditis induced by thyroglobulin.[41] Malignant transformation in the thyroid gland might be caused by cellular mediators produced by immune cells in states of chronic inflammation,[42] and autoimmune thyroiditis such as Hashimoto’s disease was reported to be associated with thyroid cancer including papillary carcinoma and thyroid lymphoma.[43,44] Further mechanistic studies are warranted to elucidate the mechanism underlying the protective effect of smoking in thyroid cancer while considering chemical composition of cigarettes such as nicotine and anatabine.

Major strengths of our study include the large sample size, which allowed us to evaluate associations among stratified subgroups, and its prospective cohort study design, which eliminates the possibility of recall bias. Data collection in the Kangbuk Samsung Health Study is standardized carefully and subject to quality control. Our study population is a relatively young, middle aged, and healthy group who are less likely to be affected by biases related to comorbidities, compared to prior studies in older and higher risk populations.

We note a few limitations that should be considered. First, thyroid cancer was primarily ascertained on self-report of physician-diagnosed thyroid cancer. In our study, the accuracy of self-reported thyroid cancer was high when it was compared with confirmed thyroid cancer from national cancer registry data, with a sensitivity of 98.1 % and a specificity of 99.8 %. Additionally, in sensitivity analysis using thyroid cancer based on the national cancer registry data, the results were similar. However, we cannot exclude the potential for measurement error in these variables, which might have resulted in some degree of residual confounding. Second, the number of current smokers in women was small and insufficient to reach a conclusion. We found that current smokers had a lower risk of thyroid cancer than never smokers in women and that the risk was gradually decreased by smoking status as a trend, but it was not statistically significant. It also may be due to the relatively low rate of female smokers in Korea or because smoking status in women was underestimated by the self-reported questionnaire. Generally, women tend to under-report smoking status, and there was a sex difference in self-reported and cotinine-verified smoking prevalence rates in Korea.[45] While this is a limitation, to the best of our knowledge, this cohort study is one of the few studies to evaluate the association between current smoking and thyroid cancer in men with a large sample size, whereas most studies only have female subjects and a relatively lower number of men. Third, lifestyle variables (i.e., smoking, alcohol intake, and physical history) were collected based on a self-administered structured questionnaire used in the health checkup program in Korea as a part of the National Health Insurance plan. We cannot exclude that measurement error in these variables, which might have resulted in some degree of residual confounding. Fourth, number of years since quitting in former smokers was not available, limiting our ability to examine the role of former smoking on thyroid cancer. Finally, our study data included apparently healthy young and middle-aged, educated Koreans. We compared our subjects with a representative sample of the general Korean population (the Korea National Health and Nutrition Examination Survey, KNHANES). The age and sex standardization was performed using the direct method to the age structure of the Korean population aged 20–80 in the year 2010. The age and sex-standardized prevalence of type 2 diabetes (defined as fasting serum glucose level of 126 mg/dL or greater or the use of blood glucose lowering agents), hypertension, obesity (BMI ≥ 25 kg/m2), and current smoker was lower than that of the general population (6.4 vs. 10.5%, 23.7 vs. 29.1%, 28.4 vs. 31.5%, and 18.0 vs. 26.5%, respectively).[46] Therefore, it is highly probable that our study subjects were healthier than the general Korean population. Our findings might not be generalizable to other populations with different characteristics of age and race/ethnicity.

In conclusion, in this prospective cohort study, current smoking was associated with a reduced risk of thyroid cancer in Korean men. Current smoking and/or having high pack-years of exposure was associated with a reduced risk of thyroid cancer even after adjusting for TSH and BMI levels as potential mediators, indicating that the association between current smoking and development of thyroid cancer may not fully elucidated by TSH and BMI levels. Further mechanistic studies underlying the inverse association between smoking and thyroid cancer will help to understand the pathogenesis of thyroid cancer. However, our findings should not be interpreted as indicating that smoking can be recommended to prevent the development of thyroid cancer because smoking is still a major risk factor for various types of cancer and poor prognosis.

Supplementary Table 1 to Table 4