Smoking and Swedish smokeless tobacco (snus) are associated with latent autoimmune diabetes in adults (LADA) and type 2 diabetes (T2D). Our aim was to investigate whether genetic susceptibility to T2D, insulin resistance (IR), and insulin secretion (IS) aggravate these associations.
We used data from two population-based Scandinavian studies with case subjects with LADA (n = 839) and T2D (n = 5,771), matched control subjects (n = 3,068), and 1,696,503 person-years at risk. Pooled, multivariate relative risks (RR) with 95% CI were estimated for smoking/genetic risk scores (T2D-GRS, IS-GRS, and IR-GRS), and ORs for snus or tobacco/GRS (case-control data). We estimated additive (proportion attributable to interaction [AP]) and multiplicative interaction between tobacco use and GRS.
The RR of LADA was elevated in high IR-GRS heavy smokers (≥15 pack-years; RR 2.01 [CI 1.30, 3.10]) and tobacco users (≥15 box/pack-years; RR 2.59 [CI 1.54, 4.35]) compared with low IR-GRS individuals without heavy use, with evidence of additive (AP 0.67 [CI 0.46, 0.89]; AP 0.52 [CI 0.21, 0.83]) and multiplicative (P = 0.003; P = 0.034) interaction. In heavy users, there was additive interaction between T2D-GRS and smoking, snus, and total tobacco use. The excess risk conferred by tobacco use did not differ across GRS categories in T2D.
Tobacco use may confer a higher risk of LADA in individuals with genetic susceptibility to T2D and insulin resistance, whereas genetic susceptibility does not seem to influence the increased T2D incidence associated with tobacco use.
Introduction
Latent autoimmune diabetes in adults (LADA) is a common, heterogenous form of diabetes distinguished by a type 1 diabetes–like autoimmune reaction against the insulin-producing β-cells, and a type 2–like phenotype often including obesity and insulin resistance (1). Genetically, LADA is closely related to type 1 diabetes, with a high risk conferred by the HLA gene complex. A genetic overlap with type 2 diabetes exists but is less evident (2). For example, associations with variants of the type 2–associated TCF7L2 gene have been observed, primarily in individuals with glutamic acid decarboxylase antibody (GADA) titers below median (3). The hybrid nature of LADA suggests that environmental factors having either a triggering or promoting effect on autoimmunity, and/or factors known to increase insulin resistance, could enhance the risk. To date, lifestyle risk factors are scarcely investigated in LADA, and neither is their potential interaction with genetic susceptibility (1). However, we recently observed an increased risk of LADA in smokers and users of Swedish smokeless tobacco (snus), using data from two observational studies and a Mendelian randomization study (4), and the excess risk was comparable to that seen in type 2 diabetes (5). Moreover, smoking was associated with insulin resistance in LADA, which is in line with observations in type 2 diabetes (6).
Previous studies on smoking and genetic susceptibility in type 2 diabetes have primarily been carried out in Asian populations, where some have revealed interaction between smoking and type 2–associated single nucleotide polymorphisms (SNPs) (7,8). Because the genetic landscape of type 2 diabetes includes many common genetic variants, each with minor effects (9), susceptibility to type 2 diabetes may be better captured by combining several variants into a genetic risk score (GRS). Interaction between such GRS and smoking has been observed in Asian studies (10,11), but not in a study of European-ancestry populations (12).
Our aim is to investigate whether GRS of type 2 diabetes (T2D-GRS), impaired insulin secretion (IS-GRS), and insulin resistance (IR-GRS) are associated with LADA and to explore potential interaction between tobacco use and GRS in relation to LADA and type 2 diabetes.
Research Design and Methods
The ESTRID Study
Study Population
The Swedish case-control study Epidemiological Study of Risk Factors for LADA and Type 2 Diabetes (ESTRID) collects data on case subjects with LADA and type 2 diabetes from a diabetes register and biobank (All New Diabetics in Scania [ANDIS]) (13) covering >90% of the health care system in the county of Scania, Sweden. In addition, ESTRID enrolls incidence-density sampled control individuals without diabetes through the Swedish Population Register. All participants fill out a health and lifestyle questionnaire, administered close to diagnosis of case subjects. Genetic information was available for case subjects but not control subjects. We therefore include population-based control subjects from the Epidemiological Investigation of Rheumatoid Arthritis (EIRA) study (14), for whom we have both genetic and questionnaire information. These control subjects are free from rheumatoid arthritis and diabetes at inclusion.
We include all case subjects with LADA (n = 594 [n = 345 with genetic information]) and type 2 diabetes (n = 2,045 [n = 1,122]) enrolled in ESTRID 2010–2019, and all control subjects aged ≥35 years in EIRA 2006–2014 (n = 3,068 [n = 1,537]) with information on tobacco use. All participants gave informed consent, and ethical approvals were obtained from the regional ethical review board in Stockholm. The study design is illustrated in Supplementary Figure 1.
Identification and Classification of Cases
Case subjects were diagnosed within the health care system of Scania, and fasting blood samples were drawn at registration. GADA concentrations were measured using an enzyme-linked immunosorbent assay (ELISA; RSR Limited, Cardiff, U.K.). At the chosen cutoff (≥10.7 U/mL), sensitivity and specificity were 0.84/0.98. Values above 250 U/mL were censored. LADA case subjects were aged ≥35 years at diagnosis, GADA positive (≥10 U/mL), and with C-peptide of ≥0.20 nmol/L (IMMULITE 2000; Siemens Healthcare Diagnostics Product Ltd., Llanberis, U.K.) or ≥0.30 nmol/L (Cobas e601; Roche Diagnostics, Mannheim, Germany). Except for C-peptide, used to distinguish LADA from type 1 diabetes, this is consistent with the criteria put forward by the Immunology of Diabetes Society (15), who proposes a minimum age cutoff of 30 years to exclude type 1 diabetes. We arbitrarily chose to use age 35 years to minimize such misclassification. Case subjects with type 2 diabetes were aged ≥35 years, GADA negative, with C-peptide of >0.60 nmol/L (IMMULITE) or >0.72 nmol/L (Cobas).
Genotyping Procedures
Genotyping in ANDIS was done on frozen DNA samples using iPlex (Sequenom, San Diego, CA) or TaqMan assays (Thermo Fisher Scientific, Carlsbad, CA) at the Clinical Research Centre in Malmö, Sweden (13). DNA samples from the EIRA biobank were analyzed with the Illumina Global Screening array or an Infinium Illumina 300 K immunochip custom array (Illumina, San Diego, CA).
The HUNT Study
Study Population
The Trøndelag Health Study (HUNT) (16) invited the entire population ≥20 years of age in the Nord-Trøndelag region, Norway, to an extensive health survey on up to four occasions: HUNT1 (1984–1986), HUNT2 (1995–1997), HUNT3 (2006–2008), and HUNT4 (2017–2019). Participants responded to questionnaires and underwent clinical examination and blood sampling. Eligible for inclusion in the current study were all participants with baseline data in HUNT1 to HUNT3 (n = 94,489). Excluding those with diabetes at baseline (n = 2,481) and those without smoking information (n = 13,466), the analytical sample consisted of 78,542 unique individuals and 1,696,503 person-years (Supplementary Figure 1). The study was approved by the Norwegian Data Protection Authority, Oslo, and the Regional Committee for Medical and Health Research Ethics, Trondheim, and all participants provided informed consent.
Identification and Classification of Case Subjects
Identification of case subjects with diabetes was done through self-report (17). GADA was assessed in blood samples collected at follow-up (median 5 years after diagnosis), analyzed by either immunoprecipitation radioligand assay (Novo Nordisk, Bagsværd, Denmark; HUNT2 and HUNT3) (18) or ELISA (RSR Limited, Cardiff, U.K.; HUNT4) with sensitivity and specificity of 0.64/1.00 and 0.84/0.98, respectively.
To harmonize with ESTRID, case subjects aged ≥35 years were classified as LADA (n = 245 [n = 219 with genetic information]) if they were GADA positive (≥0.08 antibody index [World Health Organization; ≥43 U/mL] in HUNT2 and HUNT3, ≥10 U/mL in HUNT4) and as type 2 diabetes (n = 3,726 [n = 3,421]) if GADA negative. The classification excluded C-peptide, since measurements were not from the time of diagnosis.
Genotyping Procedures
DNA was extracted and genotyped at the Norwegian University of Science and Technology Genomic Core Facility, Trondheim, by HumanCoreExome, Illumina Inc. (San Diego, CA).
GRS
We reproduced previously published GRS of type 2 diabetes (64 SNPs), insulin secretion (16 SNPs), and insulin resistance (5 SNPs) (13). Conferred genotypes were obtained as hard-call dosages with the threshold set to 0.5. The GRS were calculated based on the number of risk alleles in each score weighted by their effect sizes, as reported in previous genome-wide association studies (GWAS) (13) using Plink software (19). Unavailable SNPs (n = 11 [17.2%] in T2D-GRS and n = 1 [6.3%] in IS-GRS in ESTRID; n = 1 [1.6%] in T2D-GRS in HUNT) were replaced by proxy SNPs with r2 >0.8, except one with r2 of 0.61 in the T2D-GRS (ESTRID only). The resulting GRS were z-score standardized when used as continuous exposures.
Definition of Exposures
GRS were treated as continuous variables and measured as increase per SD and as categorical or binary variables based on tertiles. Smoking and snus use were self-reported and analyzed both categorically (smoking/snus use habit; heavy use [≥15 pack-/box-years]) and for continuity (number of pack-/box-years in ever users; current intensity). One pack/box-year was defined as smoking one pack of cigarettes (20 cigarettes) or using seven boxes of snus per week for a year. Total tobacco use was defined as being either a smoker or snus user, and the total number of tobacco-years was the sum of pack-years and box-years. We set the index date for tobacco use to one year prior to diagnosis/participation in ESTRID.
Self-reported educational level (primary school, upper secondary school, or university) and alcohol intake were included as covariates. In ESTRID, beverage-specific responses about the amount and frequency consumed during the preceding year were available. Corresponding questions in HUNT related to frequency of intake during the past 2 weeks up to the past year. Alcohol intake was divided into six and four categories, respectively, in ESTRID and HUNT, ranging from abstainers to high consumers. We also included BMI as a continuous covariate and added a quadratic term. In ESTRID, self-reported weight and height were used for the calculation of BMI (weight [kg]/height [m2]), which correlates well (r = 0.92) with BMI measured at diagnosis (20). In HUNT, BMI was based on weight and height measured at the clinical examination. Adjustment for waist-to-hip ratio (waist [cm]/hip [cm]), available for a subgroup in HUNT, physical activity, and family history of diabetes did not markedly alter the estimates and were not included in the final models. Exposure variables and covariates were updated, if possible, with the information at each survey for HUNT participants who participated in more than one survey.
Statistical Analysis
Differences in baseline characteristics were calculated by Student’s t test for means (± SD), Kruskal-Wallis test for medians (interquartile range), and χ2 test for proportions, and expressed as two-sided P values.
We conducted logistic regression conditioned on age and sex to estimate odds ratios (ORs) with 95% CIs for case-control data. On the cohort data, we performed Cox regression with age as the underlying timescale to estimate hazard ratios with 95% CIs. Participants contributed with person-years up to diabetes diagnosis, death, emigration out of the county, or end of follow-up, whichever came first. Associations between GRS and LADA/type 2 diabetes were adjusted for age and sex in the main model (model 1), and separate estimates were given for LADAhigh (GADA ≥ median) and LADAlow (GADA < median). In sensitivity analyses, GRS were adjusted for HLA (model 2) and additionally, in HUNT, for the first 10 principal components and batch number (model 3). The association with tobacco use, as well as interaction analyses, were adjusted for BMI, educational level, and alcohol intake in addition to age and sex, and analyses of snus use were further adjusted for smoking. To account for the association between BMI and insulin resistance, the interaction between smoking and IR-GRS was further stratified by BMI (<30 or ≥30). For those with missing information on any of the covariates, we performed imputation. In ESTRID, this was based on the median value, and, in HUNT, we first used the value from a previous survey, if available, and, secondly, the median value. We included indicator variables for missing covariate information.
The exposure in the models of additive interaction was the combination of binary variables of tobacco use and GRS, resulting in four combinations, with those with the lowest category of tobacco use and GRS serving as the reference group. Additive interaction was defined as departure from additivity of effects and estimated as attributable proportion (AP) due to interaction with 95% CI, representing the proportion in the doubly exposed case subjects that can be attributed to the interaction between the two exposures.
Multiplicative interaction was examined by including product terms for binary tobacco use and GRS. By fitting these models, we derived estimates contrasting the risk of LADA and type 2 diabetes associated with tobacco use in those with high versus low/intermediate genetic susceptibility. A P value of <0.05 for the product term coefficient indicates a statistically significant difference in the association across strata of GRS. The strongest type of interaction is one which is present on both the additive and multiplicative scales (21). When possible, estimates from both studies were meta-analyzed using a fixed-effects model, resulting in relative risks (RRpooled) for the combined risk. Corresponding APs were obtained by pooling the study-specific APs (APpooled) based on their respective weights, and we calculated P values for the pooled coefficients of the product terms using the log-transformed estimates and CIs. All analyses were performed in SAS 9.4 (SAS Institute, Cary, NC).
Data and Resource Availability
The analyzed data sets are available from the corresponding author upon reasonable request (ESTRID) and with permission from the HUNT Study by applying to the HUNT Study data access committee.
Results
Characteristics
The two study populations consisted of 839 individuals with LADA, 5,771 with type 2 diabetes, and 3,068 control subjects (ESTRID), along with 1,696,503 person-years at risk (HUNT). Those with LADA were, on average, less insulin resistant (HOMA-IR) and had worse β-cell function (HOMA-B) and lower C-peptide concentrations than individuals with type 2 diabetes. They were also more likely to have a high-risk HLA genotype and to be treated with insulin (Table 1). From here onward, we refer to the fully adjusted models, unless otherwise stated.
Characteristics of the study populations
Characteristic . | ESTRID . | HUNT . | ||||||
---|---|---|---|---|---|---|---|---|
Control subjects . | LADA . | Type 2 diabetes . | P . | No diabetes . | LADA . | Type 2 diabetes . | P . | |
Individuals, n | 3,068 | 594 | 2,045 | 78,542 | 245 | 3,726 | ||
Men, % | 27.8 | 53.0 | 60.1 | 0.002 | 48.1 | 47.8 | 53.2 | 0.093 |
Age at diagnosis* | 56.1 (10.3) | 59.1 (12.3) | 63.2 (10.4) | <0.0001 | — | 59.1 (11.5) | 60.3 (11.0) | 0.104 |
Age at baseline | — | — | — | 45.3 (17.0) | 44.7 (13.8) | 44.7 (12.8) | 0.996 | |
BMI at baseline | 25.4 (4.1) | 28.5 (5.6) | 31.2 (5.4) | <0.0001 | 25.2 (3.9) | 28.0 (4.7) | 27.9 (4.5) | 0.865 |
With insulin treatment†, % | — | 40.0 | 6.0 | <0.0001 | — | 42.0 | 13.6 | <0.0001 |
C-peptide, nmol/L | — | 0.72 (0.46, 1.20) | 1.20 (0.97, 1.60) | <0.0001 | — | 0.57 (0.20, 0.97) | 0.86 (0.60, 1.20) | <0.0001 |
HOMA-IR‖ | — | 2.82 (1.83, 4.44) | 3.56 (2.73, 4.78) | <0.0001 | — | 2.04 (1.18, 2.90) | 2.29 (1.65, 3.21) | 0.022 |
HOMA-β‖ | — | 40.5 (15.0, 69.0) | 71.1 (44.3, 95.9) | <0.0001 | — | 59.0 (36.7, 89.4) | 69.7 (45.5, 99.7) | 0.084 |
High T2D-GRS‡, % | 29.3 | 34.8 | 41.7 | 0.022 | 33.5 | 33.3 | 43.2 | 0.004 |
High IS-GRS‡, % | 29.3 | 34.8 | 39.8 | 0.097 | 33.6 | 37.4 | 39.5 | 0.553 |
High IR-GRS‡, % | 31.4 | 34.2 | 35.3 | 0.710 | 33.3 | 37.0 | 37.9 | 0.797 |
HLA high-risk§, % | 33.7 | 60.5 | 31.3 | <0.0001 | 29.0 | 48.9 | 27.2 | <0.0001 |
Characteristic . | ESTRID . | HUNT . | ||||||
---|---|---|---|---|---|---|---|---|
Control subjects . | LADA . | Type 2 diabetes . | P . | No diabetes . | LADA . | Type 2 diabetes . | P . | |
Individuals, n | 3,068 | 594 | 2,045 | 78,542 | 245 | 3,726 | ||
Men, % | 27.8 | 53.0 | 60.1 | 0.002 | 48.1 | 47.8 | 53.2 | 0.093 |
Age at diagnosis* | 56.1 (10.3) | 59.1 (12.3) | 63.2 (10.4) | <0.0001 | — | 59.1 (11.5) | 60.3 (11.0) | 0.104 |
Age at baseline | — | — | — | 45.3 (17.0) | 44.7 (13.8) | 44.7 (12.8) | 0.996 | |
BMI at baseline | 25.4 (4.1) | 28.5 (5.6) | 31.2 (5.4) | <0.0001 | 25.2 (3.9) | 28.0 (4.7) | 27.9 (4.5) | 0.865 |
With insulin treatment†, % | — | 40.0 | 6.0 | <0.0001 | — | 42.0 | 13.6 | <0.0001 |
C-peptide, nmol/L | — | 0.72 (0.46, 1.20) | 1.20 (0.97, 1.60) | <0.0001 | — | 0.57 (0.20, 0.97) | 0.86 (0.60, 1.20) | <0.0001 |
HOMA-IR‖ | — | 2.82 (1.83, 4.44) | 3.56 (2.73, 4.78) | <0.0001 | — | 2.04 (1.18, 2.90) | 2.29 (1.65, 3.21) | 0.022 |
HOMA-β‖ | — | 40.5 (15.0, 69.0) | 71.1 (44.3, 95.9) | <0.0001 | — | 59.0 (36.7, 89.4) | 69.7 (45.5, 99.7) | 0.084 |
High T2D-GRS‡, % | 29.3 | 34.8 | 41.7 | 0.022 | 33.5 | 33.3 | 43.2 | 0.004 |
High IS-GRS‡, % | 29.3 | 34.8 | 39.8 | 0.097 | 33.6 | 37.4 | 39.5 | 0.553 |
High IR-GRS‡, % | 31.4 | 34.2 | 35.3 | 0.710 | 33.3 | 37.0 | 37.9 | 0.797 |
HLA high-risk§, % | 33.7 | 60.5 | 31.3 | <0.0001 | 29.0 | 48.9 | 27.2 | <0.0001 |
Data are shown as mean ± SD or median (interquartile range) unless stated otherwise. Clinical information (C-peptide) was available for 98% of case subjects in ESTRID (LADA = 584, type 2 diabetes n = 1,997), and HOMA was available for 85% (LADA n = 477, type 2 diabetes = 1,758). C-peptide and HOMA (not from diagnosis) were available in HUNT1 through HUNT3 for 41% (LADA = 122, type 2 diabetes = 1,513) and 33% (LADA = 83; type 2 diabetes = 1,226), respectively. In ESTRID, 56% of case subjects and 50% of control subjects had genetic information. In HUNT, genetic information was available for 92% of case subjects and 75% of diabetes-free individuals. The P value is shown for the difference between LADA and type 2 diabetes.
Age at participation for control subjects.
Current use of insulin.
Homeostasis model assessment of insulin resistance (HOMA-IR) and β-cell function (HOMA-β) were calculated based on fasting plasma glucose and serum C-peptide from diagnosis (ESTRID) and from follow-up (median 5 years after diagnosis, HUNT).
Proportion of individuals with a GRS equal to or above the upper tertile.
High-risk genotypes defined as DR3/3, DR3/4, DR4/4, or haplotypes of DR4-DQ8 or DR3-DQ2 (ESTRID) or presence of at least one of the risk variants inferring either DR3-DQ2 or DR4-DQ8 (HUNT).
Genetic Susceptibility, Tobacco Use, and the Risk of LADA and Type 2 Diabetes
The T2D-GRS, IS-GRS, and IR-GRS were approximately normally distributed in the two study populations (Supplementary Figure 2). The T2D-GRS was associated with type 2 diabetes (RRpooled 1.35 [CI 1.31, 1.40] per SD increase) and LADA (RRpooled 1.15 [CI 1.05, 1.26] per SD increase) (Fig. 1). The IS-GRS (RRpooled 1.20 [CI 1.16, 1.23] per SD increase) and IR-GRS (RRpooled 1.13 [CI 1.10, 1.17] per SD increase) were associated with type 2 diabetes, but not with LADA (Fig. 1). Study-specific estimates were similar and did not change markedly with adjustment for HLA, principal components, and genotyping batch (Supplementary Tables 1 and 2). When LADA was stratified by median GADA concentrations, the T2D-GRS and IS-GRS were associated with LADAlow but not with LADAhigh, and there was no association with IR-GRS (Supplementary Table 3). Study-specific estimates pointed in the same directions (Supplementary Tables 4 and 5). Smoking, snus, or tobacco use were associated with higher risk of LADA and type 2 diabetes (Supplementary Figure 3).
RR with 95% CI for the associations between GRS and LADA and type 2 diabetes, adjusted for age and sex. High, upper tertile of GRS; intermediate, middle tertile of GRS; low, lower tertile of GRS.
RR with 95% CI for the associations between GRS and LADA and type 2 diabetes, adjusted for age and sex. High, upper tertile of GRS; intermediate, middle tertile of GRS; low, lower tertile of GRS.
Interaction Between Tobacco Use and Genetic Susceptibility in LADA
There was interaction between tobacco use and IR-GRS in the risk of LADA, which was evident in current and heavy smokers and current and heavy total tobacco users (Table 2, Fig. 2C, and Supplementary Table 6). The combination of current heavy smoking and a high IR-GRS conferred an RRpooled of 2.01 (CI 1.30, 3.10) with significant additive (APpooled 0.67 [CI 0.46, 0.89]) and multiplicative interaction (p = 0.003) (Table 2). Current, heavy total tobacco users with a high IR-GRS had an OR of 2.59 (CI 1.54, 4.35) with AP estimated at 0.52 (CI 0.21, 0.83) and a p for multiplicative interaction of 0.034 (Supplementary Table 6). Study-specific analyses revealed additive interaction between heavy smoking and IR-GRS in the risk of LADA (Supplementary Tables 7 and 8). The results did not differ markedly when stratified by BMI (<30 or ≥30) (Supplementary Table 9). Furthermore, there was no significant interaction between snus use and IR-GRS (Supplementary Table 6).
RR with 95% CI for LADA (A–C) and type 2 diabetes (D–F) by combinations of ever heavy smoking (≥15 vs <15 pack-years) and GRS (upper tertile versus lower and middle tertile), adjusted for age, sex, BMI, educational level, and alcohol intake. AP with 95% CI for additive interaction and P for multiplicative interaction. Pooled data are from the ESTRID (2010–2019) and HUNT (1984–2019) studies.
RR with 95% CI for LADA (A–C) and type 2 diabetes (D–F) by combinations of ever heavy smoking (≥15 vs <15 pack-years) and GRS (upper tertile versus lower and middle tertile), adjusted for age, sex, BMI, educational level, and alcohol intake. AP with 95% CI for additive interaction and P for multiplicative interaction. Pooled data are from the ESTRID (2010–2019) and HUNT (1984–2019) studies.
Interaction between smoking and GRS of type 2 diabetes, insulin secretion, and insulin resistance and risk of LADA
GRS . | Combination of different categories of smoking and GRS . | . | Within strata of GRS . | Multiplicative interaction . | P . | |||
---|---|---|---|---|---|---|---|---|
N case subjects/control subjects/person-years . | RRpooled (95% CI) . | N case subjects/control subjects/person-years . | RRpooled (95% CI) . | Additive interaction . | RRpooled (95% CI) . | RRpooled (95% CI)* . | ||
Nonsmoker . | Current smoker . | APpooled (95% CI) . | Current vs. no smoking . | |||||
T2D-GRSlow/intermediate | 282/879/669,231 | 1.00 (Ref.) | 89/208/270,160 | 1.30 (0.99, 1.71) | 1.27 (0.96, 1.67) | |||
T2D-GRShigh | 151/370/327,629 | 1.22 (0.97, 1.53) | 42/80/138,907 | 1.33 (0.91, 1.93) | −0.08 (−0.57, 0.41) | 1.12 (0.74, 1.69) | 0.84 (0.52, 1.35) | 0.163 |
IS-GRSlow/intermediate | 271/874/662,103 | 1.00 (Ref.) | 91/209/272,347 | 1.38 (1.05, 1.82) | 1.36 (1.03, 1.81) | |||
IS-GRShigh | 162/375/334,757 | 1.29 (1.04, 1.62) | 40/79/136,721 | 1.25 (0.85, 1.83) | −0.33 (−0.91, 0.25) | 0.98 (0.65, 1.47) | 0.71 (0.44, 1.15) | 0.227 |
IR-GRSlow/intermediate | 288/851/666,029 | 1.00 (Ref.) | 77/203/271,698 | 1.08 (0.81, 1.44) | 1.11 (0.83, 1.48) | |||
IR-GRShigh | 145/398/330,832 | 1.10 (0.88, 1.38) | 54/85/137,369 | 1.71 (1.20, 2.43) | 0.51 (0.26, 0.76) | 1.54 (1.04, 2.27) | 1.39 (0.87, 2.22) | 0.169 |
Smokes <15 pack-years (ever) | Smokes ≥15 pack-years (ever) | ≥15 vs. <15 pack-years (ever) | ||||||
T2D-GRSlow/intermediate | 230/845/707,571 | 1.00 (Ref.) | 81/242/147,373 | 1.05 (0.79, 1.40) | 1.04 (0.77, 1.39) | |||
T2D-GRShigh | 112/369/350,282 | 1.06 (0.84, 1.62) | 50/81/76,254 | 1.66 (1.16, 2.38) | 0.38 (0.09, 0.67) | 1.51 (1.01, 2.27) | 1.56 (0.97, 2.51) | 0.067 |
IS-GRSlow/intermediate | 222/849/700,539 | 1.00 (Ref.) | 92/234/150,121 | 1.23 (0.92, 1.63) | 1.21 (0.91, 1.62) | |||
IS-GRShigh | 120/365/357,313 | 1.21 (0.97, 1.53) | 39/89/73,506 | 1.44 (0.99, 2.11) | 0.07 (−0.35, 0.50) | 1.15 (0.76, 1.74) | 1.01 (0.63, 1.64) | 0.971 |
IR-GRSlow/intermediate | 233/831/706,874 | 1.00 (Ref.) | 76/223/148,979 | 1.00 (0.74, 1.35) | 1.05 (0.77, 1.42) | |||
IR-GRShigh | 109/383/350,979 | 1.02 (0.81, 1.30) | 55/100/74,648 | 1.64 (1.17, 2.29) | 0.40 (0.12, 0.68) | 1.51 (1.03, 2.21) | 1.65 (1.04, 2.62) | 0.033 |
Smokes <15 pack-years (current) | Smokes ≥15 pack-years (current) | ≥15 vs. <15 pack-years (current) | ||||||
T2D-GRSlow/intermediate | 316/960/758,169 | 1.00 (Ref.) | 43/127/96,774 | 1.03 (0.71, 1.50) | 1.01 (0.69, 1.47) | |||
T2D-GRShigh | 159/406/375,712 | 1.12 (0.90, 1.39) | 28/44/50,824 | 1.54 (0.97, 2.43) | 0.33 (−0.07, 0.74) | 1.36 (0.84, 2.22) | 1.42 (0.78, 2.60) | 0.256 |
IS-GRSlow/intermediate | 303/958/751,926 | 1.00 (Ref.) | 47/125/98,734 | 1.15 (0.80, 1.65) | 1.12 (0.77, 1.62) | |||
IS-GRShigh | 172/408/381,955 | 1.20 (0.97, 1.49) | 24/46/48,864 | 1.37 (0.85, 2.21) | 0.03 (−0.52, 0.59) | 1.13 (0.68, 1.88) | 1.01 (0.55, 1.86) | 0.977 |
IR-GRSlow/intermediate | 320/924/757,405 | 1.00 (Ref.) | 35/130/98,447 | 0.80 (0.54, 1.19) | 0.83 (0.56, 1.24) | |||
IR-GRShigh | 155/442/376,476 | 1.04 (0.84, 1.29) | 36/41/49,151 | 2.01 (1.30, 3.10) | 0.67 (0.46, 0.89) | 1.88 (1.17, 3.01) | 2.42 (1.34, 4.38) | 0.003 |
GRS . | Combination of different categories of smoking and GRS . | . | Within strata of GRS . | Multiplicative interaction . | P . | |||
---|---|---|---|---|---|---|---|---|
N case subjects/control subjects/person-years . | RRpooled (95% CI) . | N case subjects/control subjects/person-years . | RRpooled (95% CI) . | Additive interaction . | RRpooled (95% CI) . | RRpooled (95% CI)* . | ||
Nonsmoker . | Current smoker . | APpooled (95% CI) . | Current vs. no smoking . | |||||
T2D-GRSlow/intermediate | 282/879/669,231 | 1.00 (Ref.) | 89/208/270,160 | 1.30 (0.99, 1.71) | 1.27 (0.96, 1.67) | |||
T2D-GRShigh | 151/370/327,629 | 1.22 (0.97, 1.53) | 42/80/138,907 | 1.33 (0.91, 1.93) | −0.08 (−0.57, 0.41) | 1.12 (0.74, 1.69) | 0.84 (0.52, 1.35) | 0.163 |
IS-GRSlow/intermediate | 271/874/662,103 | 1.00 (Ref.) | 91/209/272,347 | 1.38 (1.05, 1.82) | 1.36 (1.03, 1.81) | |||
IS-GRShigh | 162/375/334,757 | 1.29 (1.04, 1.62) | 40/79/136,721 | 1.25 (0.85, 1.83) | −0.33 (−0.91, 0.25) | 0.98 (0.65, 1.47) | 0.71 (0.44, 1.15) | 0.227 |
IR-GRSlow/intermediate | 288/851/666,029 | 1.00 (Ref.) | 77/203/271,698 | 1.08 (0.81, 1.44) | 1.11 (0.83, 1.48) | |||
IR-GRShigh | 145/398/330,832 | 1.10 (0.88, 1.38) | 54/85/137,369 | 1.71 (1.20, 2.43) | 0.51 (0.26, 0.76) | 1.54 (1.04, 2.27) | 1.39 (0.87, 2.22) | 0.169 |
Smokes <15 pack-years (ever) | Smokes ≥15 pack-years (ever) | ≥15 vs. <15 pack-years (ever) | ||||||
T2D-GRSlow/intermediate | 230/845/707,571 | 1.00 (Ref.) | 81/242/147,373 | 1.05 (0.79, 1.40) | 1.04 (0.77, 1.39) | |||
T2D-GRShigh | 112/369/350,282 | 1.06 (0.84, 1.62) | 50/81/76,254 | 1.66 (1.16, 2.38) | 0.38 (0.09, 0.67) | 1.51 (1.01, 2.27) | 1.56 (0.97, 2.51) | 0.067 |
IS-GRSlow/intermediate | 222/849/700,539 | 1.00 (Ref.) | 92/234/150,121 | 1.23 (0.92, 1.63) | 1.21 (0.91, 1.62) | |||
IS-GRShigh | 120/365/357,313 | 1.21 (0.97, 1.53) | 39/89/73,506 | 1.44 (0.99, 2.11) | 0.07 (−0.35, 0.50) | 1.15 (0.76, 1.74) | 1.01 (0.63, 1.64) | 0.971 |
IR-GRSlow/intermediate | 233/831/706,874 | 1.00 (Ref.) | 76/223/148,979 | 1.00 (0.74, 1.35) | 1.05 (0.77, 1.42) | |||
IR-GRShigh | 109/383/350,979 | 1.02 (0.81, 1.30) | 55/100/74,648 | 1.64 (1.17, 2.29) | 0.40 (0.12, 0.68) | 1.51 (1.03, 2.21) | 1.65 (1.04, 2.62) | 0.033 |
Smokes <15 pack-years (current) | Smokes ≥15 pack-years (current) | ≥15 vs. <15 pack-years (current) | ||||||
T2D-GRSlow/intermediate | 316/960/758,169 | 1.00 (Ref.) | 43/127/96,774 | 1.03 (0.71, 1.50) | 1.01 (0.69, 1.47) | |||
T2D-GRShigh | 159/406/375,712 | 1.12 (0.90, 1.39) | 28/44/50,824 | 1.54 (0.97, 2.43) | 0.33 (−0.07, 0.74) | 1.36 (0.84, 2.22) | 1.42 (0.78, 2.60) | 0.256 |
IS-GRSlow/intermediate | 303/958/751,926 | 1.00 (Ref.) | 47/125/98,734 | 1.15 (0.80, 1.65) | 1.12 (0.77, 1.62) | |||
IS-GRShigh | 172/408/381,955 | 1.20 (0.97, 1.49) | 24/46/48,864 | 1.37 (0.85, 2.21) | 0.03 (−0.52, 0.59) | 1.13 (0.68, 1.88) | 1.01 (0.55, 1.86) | 0.977 |
IR-GRSlow/intermediate | 320/924/757,405 | 1.00 (Ref.) | 35/130/98,447 | 0.80 (0.54, 1.19) | 0.83 (0.56, 1.24) | |||
IR-GRShigh | 155/442/376,476 | 1.04 (0.84, 1.29) | 36/41/49,151 | 2.01 (1.30, 3.10) | 0.67 (0.46, 0.89) | 1.88 (1.17, 3.01) | 2.42 (1.34, 4.38) | 0.003 |
Pooled data from the ESTRID (2010–2019) and HUNT (1984–2019) studies. GRSlow/intermediate, lower and middle tertiles of GRS; GRShigh, upper tertile of GRS; Ref., reference category. RR adjusted for age, sex, BMI, educational level, and alcohol intake. Smokes <15 pack-years (ever) includes never-smokers. Smokes <15 pack-years (current) includes former smokers and never-smokers.
RR for the interaction term smoking × GRS.
Additive interaction was observed between the T2D-GRS and ever heavy smoking (APpooled 0.38 [CI 0.09, 0.67]), ever (AP 0.72 [CI 0.28, 1.15]) and current (AP 0.67 [CI 0.14, 1.21]) heavy snus use, and heavy total tobacco use (AP 0.49 [CI 0.20, 0.78]) in ever users (Table 2, Fig. 2A, and Supplementary Table 6). A joint exposure to heavy tobacco use and a high T2D-GRS compared with no heavy tobacco use and a low T2D-GRS in ever users resulted in an OR of 2.60 (CI 1.62, 4.17) for LADA, with evidence of additive (AP 0.49 [CI 0.20, 0.78]) and multiplicative (P = 0.038) interaction (Supplementary Table 6). Study-specific analyses of heavy smoking and T2D-GRS yielded similar results, although only significant in the Swedish data (Supplementary Tables 7 and 8). Smoking, snus use, and total tobacco use did not interact with IS-GRS in LADA (Table 2, Fig. 2B, and Supplementary Table 6).
Interaction Between Tobacco Use and Genetic Susceptibility in Type 2 Diabetes
In type 2 diabetes, there was no evidence of interaction between smoking, snus use, or total tobacco use and any of the GRS, neither on the additive nor the multiplicative scale, except for moderate additive interaction between current, heavy smoking and a high IR-GRS (APpooled 0.18 [CI 0.03, 0.33]) (Fig. 2D–F and Supplementary Tables 10 and 11). Study-specific results indicated interaction with IR-GRS in ESTRID but not in HUNT (Supplementary Tables 12 and 13). The interaction between smoking and IR-GRS did not differ when stratified by BMI (Supplementary Table 14).
Conclusions
We used data from two Scandinavian population-based studies and showed that a GRS of type 2 diabetes, but not impaired insulin secretion or insulin resistance, was associated with incidence of LADA. We also found that genetic susceptibility to type 2 diabetes or insulin resistance strengthened the association between tobacco use and LADA risk. The strongest evidence of such interaction was observed between smoking and genetic predisposition to insulin resistance, which was seen in both cohorts and visible on both the additive and multiplicative scale. Genetic susceptibility to type 2 diabetes, insulin resistance, or insulin secretion had limited influence on the association between tobacco use and incidence of type 2 diabetes.
We confirmed findings from GWAS (22) and candidate gene studies (23,24) linking genetic predisposition to type 2 diabetes to an elevated risk of developing LADA, and this was primarily seen in individuals with GADA titers below median. LADA was not associated with genetic susceptibility to impaired insulin secretion or insulin resistance, which is consistent with previous findings from the ANDIS study, where the GRS were examined in relation to severe autoimmune diabetes, including LADA and adult-onset type 1 diabetes (13).
This study is the first investigating the interaction between tobacco use and genetic susceptibility to type 2 diabetes and related traits in LADA. Genetically, LADA is a heterogeneous condition with an etiology that encompasses susceptibility to both type 1 and type 2 diabetes. Our findings indicate that smoking primarily increases the risk of LADA in individuals with a genetic predisposition to type 2 diabetes or insulin resistance. Interestingly, we recently also observed interaction between smoking and HLA genotypes associated with type 1 diabetes in relation to LADA (4). This supports the hypothesis that in LADA, the influence of environmental factors increasing insulin resistance is contingent on genetic predisposition. Whether smoking also interacts with LADA- and type 1 diabetes–associated loci outside the HLA region remains to be investigated.
In contrast to our findings in LADA, the association between smoking and type 2 diabetes seen across strata of genetic susceptibility and the limited evidence of interaction indicate that tobacco use and genetic susceptibility are independently associated with type 2 diabetes risk. Our results correspond with previous findings in populations of primarily European descent (12) but are in contrast to findings in Asian populations (10,11). These discrepancies may be due to differences in genetic risk variants, pathways, and minor allele frequencies across populations (25,26). In general, observational studies have, so far, provided inconsistent evidence regarding interaction between lifestyle factors and type 2 diabetes polygenic risk scores in type 2 diabetes (27).
Adverse effects of smoking and nicotine exposure on insulin sensitivity have been documented in experimental studies (28,29). In line with this, we recently showed that smoking, snus, and total tobacco use are positively associated with HOMA-IR in both LADA and type 2 diabetes (4). Against this background, it is interesting that the strongest evidence of interaction was observed between smoking and IR-GRS in LADA. Genetic susceptibility to insulin resistance also seems to increase the risk of LADA only in individuals who are simultaneously exposed to tobacco. In addition, there was some indication of additive interaction between heavy smoking and IR-GRS in type 2 diabetes. This suggests that inherited and acquired insulin resistance may synergistically accelerate the onset of LADA and possibly also type 2 diabetes. The stronger interactions observed with IR-GRS in LADA, and, to a lesser extent, in type 2 diabetes, indicate that specific pathways involved in insulin resistance may play a significant role in these interactions. To the best of our knowledge, this study is the first to address interaction between a trait-specific GRS and tobacco use in type 2 diabetes.
Strengths and Limitations
The strengths of this study include the large number of incident case subjects with LADA and type 2 diabetes, the population-based design, detailed lifestyle and clinical information, the use of GRS that capture overall type 2 diabetes genetic risk as well as specific diabetes-related traits, and a follow-up time of up to 35.5 years in the cohort study. We also had the opportunity to replicate the results in two independent cohorts.
Recall bias is a concern in ESTRID, since case subjects may retrospectively report tobacco habits differently from control subjects. However, the association between smoking and type 2 diabetes in ESTRID was similar to that observed in HUNT and numerous previous prospective studies (5). Moreover, self-reported smoking correlates well with gold standard measures of blood cotinine concentrations (30).
We could not separate LADA from adult-onset type 1 diabetes in HUNT, because of the lack of C-peptide measurements from diagnosis and limited information on treatment. This implies that some individuals with adult-onset type 1 diabetes are included in the LADA group. We based the LADA classification on the presence of GADA, for which >90% of case subjects are positive (31). Individuals positive only for other autoantibodies will be misclassified as having type 2 diabetes, but this number is likely small and will have minor influence on the results. However, with the specificity of the GADA assay, it follows that some case subjects with type 2 diabetes will be incorrectly categorized as having LADA. Such misclassification could explain why the T2D-GRS was associated with LADAlow and not with LADAhigh. It is, however, unlikely to be driving the observed interactions in LADA, given that these interactions were almost absent in type 2 diabetes.
The GRS used in this study were developed based on published GWAS from European cohorts, which justifies their application on our two study populations that are of primarily European ancestry. More recent GWAS have revealed previously unknown type 2 diabetes–associated loci, not captured by the current T2D-GRS, but we included the loci with the largest effect sizes. Moreover, there are methods for calculating a GRS that may better predict the outcome than a weighted sum approach. However, by replicating these GRS and the method used to calculate them, we can compare our findings with previous findings in subgroups of adult-onset diabetes (13). Our results in LADA were on par with estimates in severe autoimmune diabetes, reassuring the validity of the data and the use of proxy SNPs.
The T2D-GRS includes genetic variants associated with diverse facets of the disease, whereas IS- and IR-GRS pinpoint more specific traits (for further information on included SNPs, we refer to the GWAS catalog [https://www.ebi.ac.uk] and previous association studies). Thus, by incorporating two trait-specific GRS, we could approach potential mechanisms underlying the observed gene-environment interactions. In fact, it has been proposed that the application of GRS targeting specific biological mechanisms may be better aimed at revealing interaction effects due to the etiological complexity of type 2 diabetes (27). It can also be noted that genetic variants in the IR-GRS have been linked to smoke-induced insulin resistance (IRS1) (32) and significantly higher HOMA-IR in obese carriers (PPARG) (33).
A limitation is that we could not adjust the GRS for population genetic structures or genotyping array in the case-control data. However, when adjusting for these factors in the cohort data, it did not change the estimates for the associations between the GRS and LADA and type 2 diabetes. Moreover, estimates were similar for the associations in both data sets, strengthening reliability in the genetic data. Finally, it is important to acknowledge that these results were retrieved in European populations from relatively homogenous environments. They cannot be generalized to non-European populations because there may be differences in genetic makeup and risk factors.
Conclusion
Our results suggest that smoking prevention could potentially inhibit or prolong the development of LADA, especially in individuals with genetic predisposition to insulin resistance or type 2 diabetes. In contrast, tobacco use seems to be a risk factor for type 2 diabetes independent of genetic predisposition. Our findings highlight the importance of tobacco use prevention to reduce the individual and societal burden of LADA, a chronic disease. Replications of the findings are warranted, also in type 1 diabetes, to further promote the discussion on possible underlying mechanisms. Nevertheless, our findings contribute to the sparse body of knowledge about the interplay between environmental and genetic factors in the etiology of LADA.
This article contains supplementary material online at https://doi.org/10.2337/figshare.22087244.
Article Information
Acknowledgments. The authors thank the participants in ESTRID, EIRA, ANDIS, and HUNT as well as administrative personnel, nurses, and research team members from all the studies.
Funding. The ESTRID study was funded by grants from the Swedish Research Council (2018-03035); FORTE (2018-00337); the Novo Nordisk Foundation (NNF19OC0057274); and the Swedish Diabetes Foundation (DIA2022-735). EIRA was funded by the Swedish Research Council; FORTE; the Swedish Rheumatic Foundation; the AFA Insurance Company; and Stockholm County Council. ANDIS was financially supported by the Swedish Research Council (project grants 2020-02191 and 2015-2558 and infrastructure grants 2010-5983, 2012-5538, and 2014-6395), Linnaeus grant 349-2006-237, a strategic research grant 2009-1039 (EXODIAB), Swedish governmental funding of clinical research (ALF), and The Swedish Foundation for Strategic Research (IRC15–0067). B.R. is the recipient of a post-doctoral fellowship from the Novo Nordisk Foundation (grant NNF17OC0027580). The Trøndelag Health Study (HUNT) Study is a collaboration between HUNT Research Centre (Faculty of Medicine and Health Sciences, Norwegian University of Science and Technology [NTNU]), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health. The genetic investigation of the HUNT Study is a collaboration between investigators from the HUNT Study and University of Michigan Medical School and the University of Michigan School of Public Health. The K.G. Jebsen Center for Genetic Epidemiology is financed by Stiftelsen Kristian Gerhard Jebsen, Faculty of Medicine and Health Sciences, NTNU, and Central Norway Regional Health Authority. The genotyping was financed by the National Institutes of Health, University of Michigan, The Norwegian Research Council, and Central Norway Regional Health Authority and the Faculty of Medicine and Health Sciences, NTNU. The genotype quality control and imputation have been conducted by the K.G. Jebsen Center for Genetic Epidemiology, Department of Public Health and Nursing, Faculty of Medicine and Health Sciences, NTNU.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. Contributions to the data collection were made by J.E., B.R., and S.C. (ESTRID); E.A. and T.T. (ANDIS); L.A. (EIRA); and V.G., E.P.S., and B.O.Å. (HUNT). S.C. was responsible for conceptualizing the research objectives, designing the study, and thoroughly revising the manuscript. D.M.A. developed the GRS. J.E. and D.M.A. calculated the GRS in the study. J.E. developed the objectives of the study and was responsible for drafting the manuscript and analyzing the data. All authors approved the final version of the manuscript. J.E. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Prior Presentation. Parts of this study were presented at the following meetings: 56th Annual Meeting of the European Diabetes Epidemiology Group (EDEG), Hersonissos, Greece, 2–5 April 2022; 56th Annual Meeting of the Scandinavian Society for the Study of Diabetes, Reykjavik, Iceland, 12–13 May 2022; the 58th Annual Meeting of the European Association for the Study of Diabetes (EASD), Stockholm, Sweden, 19–23 September 2022; and the 17th International Diabetes Epidemiology Group (IDEG) Symposium, Porto, Portugal, 2–5 December 2022.