DOI: 10.2337/db06-0600 © 2006 by the American Diabetes Association
IL6 Gene Promoter Polymorphisms and Type 2 DiabetesJoint Analysis of Individual Participants Data From 21 Studies
1 GSF-Institute of Epidemiology, Neuherberg, Germany Address correspondence and reprint requests to Dr. Thomas Illig, Institute of Epidemiology, GSF-National Research Center for Environment and Health, Ingolstaedter Landstrasse 1, 85764 Neuherberg, Germany. E-mail: illig{at}gsf.de
Abbreviations:
FHS, Framingham Heart Study; HWE, Hardy-Weinberg equilibrium; IL, interleukin; IPD, individual participants data
Several lines of evidence indicate a causal role of the cytokine interleukin (IL)-6 in the development of type 2 diabetes in humans. Two common polymorphisms in the promoter of the IL-6 encoding gene IL6, –174G>C (rs1800795) and –573G>C (rs1800796), have been investigated for association with type 2 diabetes in numerous studies but with results that have been largely equivocal. To clarify the relationship between the two IL6 variants and type 2 diabetes, we analyzed individual data on >20,000 participants from 21 published and unpublished studies. Collected data represent eight different countries, making this the largest association analysis for type 2 diabetes reported to date. The GC and CC genotypes of IL6 –174G>C were associated with a decreased risk of type 2 diabetes (odds ratio 0.91, P = 0.037), corresponding to a risk modification of nearly 9%. No evidence for association was found between IL6 –573G>C and type 2 diabetes. The observed association of the IL6 –174 C-allele with a reduced risk of type 2 diabetes provides further evidence for the hypothesis that immune mediators are causally related to type 2 diabetes; however, because the association is borderline significant, additional data are still needed to confirm this finding. Recent studies have investigated the role of variants within genes encoding immune-related markers in mediating increased type 2 diabetes risk. One of the most widely studied immune genes is the interleukin (IL)-6 encoding gene IL6, which maps to chromosome 7p21. IL-6 exerts crucial effects not only in inflammation and infection but also within the nervous and endocrine systems (1). A vast number of epidemiological, genetic, rodent, and human in vivo and in vitro studies have investigated the putative role of IL-6 in the pathogenesis underlying type 2 diabetes. The impact of IL-6 on hepatocytes, skeletal muscle cells, ß-cells, and the central nervous system has been described, and both protective and pathogenic activity of IL-6 in type 2 diabetes was suggested (2,3). Functional relevance has been ascribed to several IL6 variants located in the promoter region, including –174G>C (rs1800795) and –573G>C (rs1800796, previously denoted as –572G>C), with in vitro data demonstrating unequivocally that the IL6 –174G>C sequence affects promoter strength (4,5). The relation between –174G>C and circulating IL-6 is not completely consistent in the literature. Whereas several studies indicate that –174G>C is associated with plasma levels of IL-6, particularly in inflammatory situations (6,7), no association between –174G>C and IL-6 was found within 718 nondiabetic women of the Nurses Health Study (8). Association between –174G>C and type 2 diabetes was first reported in U.S. Pima Indians and Spanish Caucasians (9), the C-allele being statistically significantly associated with a decreased risk of type 2 diabetes. One study subsequently replicated these initial findings (10), although most did not (11–14). The only major study on –573G>C was performed in Danish Caucasians and showed a significantly increased risk of type 2 diabetes by the C-allele, but the –573G>C control genotypes were not in Hardy-Weinberg equilibrium (HWE) (14). Because of the ambiguity in interpreting the role of IL6 polymorphisms in type 2 diabetes susceptibility based on these disparate reports, we assembled an international IL6–type 2 diabetes consortium in order to perform a joint analysis. The consortium utilized individual participants data (IPD) and recruited all published and unpublished data on the association of the IL6 –174G>C or –573G>C polymorphisms and type 2 diabetes. This approach overcomes many of the problems associated with meta-analyses of published estimates such as variability in study design, poor data quality, insufficient or heterogeneous confounder adjustment, and publication bias (15). As of late 2005, investigators from the U.S., Greece, Spain, Germany, U.K., Denmark, Sweden, and Finland participated in the consortium and contributed raw data on >30,000, mostly Caucasian, subjects. As such, this study is one of the largest genetic epidemiologic association studies on IPD ever conducted. The aim of this joint analysis is to provide conclusive evidence whether the two IL6 variants, –174G>C and –573G>C, are associated with risk of type 2 diabetes.
All available published and unpublished studies fulfilling the following criteria were included in this joint analysis: 1) association study conducted in humans, 2) polymorphic genotype data for IL6 –174G>C or –573G>C, 3) type 2 diabetic case and nondiabetic control subjects, 4) published before September 2005 or unpublished, and 5) availability of IPD. Studies were excluded if the control group consisted only of individuals with pre-diabetes (16) or if ethnic admixture of unrelated study subjects was reported in the original publication (Pima Indian case-control study, 9). Information on search strategy, study recruitment, data collection and cleaning, and genotyping methods is provided in the online appendix (available at http://diabetes.diabetesjournals.org).
Definition of analyzed samples.
Statistical analyses. Study-specific odds ratios (ORs) with SEs for association between IL6 variants and type 2 diabetes were estimated from the IPD by logistic regression for case-control comparisons (SAS Proc Logistic) and by conditional logistic regression for discordant-sib comparisons (SAS Proc Phreg). The correlation due to linkage between disease status and investigated variants among sibs sharing the same marker alleles was accounted for by a jackknife variance estimate (18). All analyses were adjusted for age, sex, and BMI. Effect modification by BMI (quantitative and dichotomized at 28 kg/m2) and sex was tested.
As the CC genotype of IL6 –573G>C was rare (<1.5% in all studies), C-allele carriers (CC and GC genotypes) were compared with GG subjects. For –174G>C, ORs comparing either CC or GC with the wild-type GG were calculated, according to which the appropriate genetic model was chosen for the main analysis. Between-study heterogeneity was tested by the For the summary OR, study-specific ORs were combined by using the inverse-variance fixed-effect and the DerSimonian and Laird random-effects models. As the heterogeneity between study-specific ORs was low in all main analyses, the two models provided identical or very similar results. Thus, only the fixed-effect results are reported. The summary ORs of all studies where the control group was in HWE are reported as main results. Publication bias was investigated by visual inspection of funnel plots and formally tested using Eggers regression method (20). Funnel and forest plots were prepared using Review Manager software version 4.2 (Cochrane Collaboration, Copenhagen, Denmark).
For the IL6 –174G>C polymorphism, 10 published studies met the inclusion criteria. All of these studies, with the exception of the Framingham Heart Study (FHS) (21), provided IPD and were included in the joint analysis. Additionally, 12 unpublished studies were available for –174G>C and included in our analyses. For –573G>C, only one published study was available. However, data from eight unpublished studies met our inclusion criteria and were additionally used in our analyses. Data from 30,636 (–174G>C) and 21,352 (–573G>C) individuals were initially compiled in the central database; 22,626 and 17,305 subjects met the requirements for the analyzed samples, respectively. Except for one discordant-sib study on admixed Pima Indians, all studies consisted of Caucasian subjects.
Study-specific descriptive statistics.
IL6 –174G>C polymorphism and risk of type 2 diabetes. Figure 1A shows the ORs and 95% CIs for 18 individual studies for the association between IL6 –174 C-allele dominant and type 2 diabetes, adjusted for age, sex, and BMI. The pooled OR for 4,746 case and 16,230 control subjects was 0.91 (P = 0.037); I2, the impact of heterogeneity, was 0% (95% CI 0–50). The dominant genetic model appeared most consistent with the data, the pooled model-free ORGCvsGG and ORCCvsGG being 0.92 (0.83–1.01) and 0.90 (0.80–1.01), respectively, and was thus chosen for the main analysis. Visual inspection of the funnel plot of all 18 studies showed that studies with high, as well as low, precision of the OR estimate were distributed symmetrically around the pooled OR (online appendix Figure A1). Thus, no publication bias is suggested, which was further supported by the nonsignificant Eggers regression test (P = 0.71).
Two studies were not included in the main analysis due to HWE violation in the control groups and one (GREEK_CC study) because BMI adjustment was not possible. Their study-specific ORs were 1.9 (95% CI 1.2–3.0) for the RMIFAM_DS study, 0.6 (0.3–1.3) for the TGN_CC study (adjusted for age, sex, and BMI), and 2.2 (0.7–7.1) for the GREEK_CC study (adjusted for age and sex). Sensitivity analyses, including these studies or showing the impact of BMI adjustment, are presented in Table 2. Further sensitivity analyses are presented in online appendix Table A4; no major difference was found between case-control/discordant-sib studies, between studies which originally were designed/not designed as type 2 diabetes studies, between studies which used/did not use an oral glucose tolerance test to exclude subjects with impaired glucose tolerance from the control subjects, between studies enriched/not enriched for myocardial infarction patients, and between published/unpublished studies, respectively. Analyzing men and women separately also did not appreciably affect the size of the pooled OR. Likewise, there was no major change when excluding each study at a turn, with the pooled ORs ranging between 0.89 (0.81–0.99) and 0.92 (0.84–1.01) (online appendix Figure A2 A). There was no evidence that BMI (P > 0.4, no evidence for heterogeneity between studies) or sex (P = 0.93, no heterogeneity) significantly modified the relationship between IL6 –174G>C and type 2 diabetes.
IL6 –573G>C polymorphism and risk of type 2 diabetes. Figure 1B shows the ORs and 95% CIs for eight individual studies for the association between IL6 –573 C-allele dominant and type 2 diabetes, adjusted for age, sex, and BMI. The pooled OR for 2,392 case and 9,265 control subjects was 1.05 (P = 0.65); I2 was estimated as 0% (95% CI 0–68). The DANISH_CC study (OR 1.7 [95% CI 1.3–2.2]) was not included in this main analysis because control genotypes for –573G>C were not in HWE. Sensitivity analyses, presented in Table 2, show that heterogeneity between studies was substantially reduced by eliminating the DANISH_CC study and by adjusting for BMI (reduction from I2 = 44.8% [P = 0.07] to I2 = 0.0% [P = 0.46]). Further sensitivity analyses for subgroups of studies and stratification for sex show no remarkable change in the pooled result (online appendix Table A5). Removing each study at a turn yielded pooled ORs ranging between 0.98 and 1.15 with 95% CIs that always included unity, indicating that the pooled OR was not unduly influenced by any single study (online appendix Figure A2 B). There was no effect modification of BMI (P > 0.6) or sex (P = 0.28) on the relationship between –573G>C and type 2 diabetes.
The results presented here, based on IPD from 5,601 type 2 diabetic case and 17,019 control subjects and representing 21 association studies, provide evidence that the IL6 –174G>C polymorphism is associated with type 2 diabetes and that individuals carrying the C-allele have a 9% lower odds of suffering from type 2 diabetes compared with individuals with the GG genotype (P = 0.037). We did not find a statistically significant relationship between IL6 –573G>C and type 2 diabetes. It is plausible that the shown association of –174G>C with type 2 diabetes reflects a true modulating effect of –174G>C or another variant in linkage disequilibrium with –174G>C. The closest known gene (TOMM7) is situated about 100 kb from IL6 and is located within a different linkage disequilibrium block (http://www.hapmap.org).
Putative impact of unincluded studies. Since the deadline for inclusion of newly published studies has elapsed until today (June 2006), only two large studies (>500 participants) fulfilling the inclusion criteria for our joint analysis have been published (8). Their study-specific ORs for association between IL6 –174 C-allele dominant and type 2 diabetes, adjusted for age and BMI, were 0.95 (95% CI 0.82–1.10) for the Nurses Health Study (1,315 female case and 2,265 female control subjects) and 0.95 (0.77–1.17) for the Health Professional Follow-up Study (885 male case and 894 male control subjects) (Dr. Lu Qi, personal communication). The pooled OR for the joint analysis, including these studies, was 0.92 (0.86–0.99) and had a slightly lower P value of 0.030 than our main analysis.
Analysis strategy.
Strengths and limitations of this joint analysis. Joint analyses based on IPD have several advantages compared with meta-analyses that are based on published estimates or summary data (15). Here, standardized methods were applied, incoming data were checked and cleaned, genotypes were tested for HWE violation, putative confounders for type 2 diabetes were uniformly adjusted for, stratified and interaction analyses were performed, and a consistent genetic model was applied. The observed low heterogeneity among studies may have resulted from these standardized procedures. The greatest limitation of any meta-analysis is the risk of publication bias. To avoid this bias, we have strived to include all existing data involving the IL6 –174G>C and –573G>C variants and type 2 diabetes susceptibility and managed to include predominantly unpublished data. Nevertheless, we conducted analyses to assess the effect of publication bias on our results for –174G>C. Utilizing the funnel plot and Eggers regression test, there was no evidence for publication bias, suggesting that our study sample is comprised of a representative dataset. Although this study including >20,000 subjects is among the largest genetic association studies performed to date on IPD, the observed inverse association of the IL6 –174 C-allele with type 2 diabetes, showing a P value of 0.037, is borderline significant. Bonferroni correction for the two analyzed single nucleotide polymorphisms would turn the result to statistical nonsignificance. However, the ORs of the recently published Nurses Health Study and Health Professional Follow-up Study point in the same direction as our joint analysis, thus adding strength to the reported association of –174G>C with type 2 diabetes. The weak OR of 0.91 is plausible, as type 2 diabetes is a complex disease whose etiology is dependent upon multiple genetic and environmental factors and consistent with estimates obtained in other genes that affect susceptibility to type 2 diabetes (24). In conclusion, this joint analysis is the largest association study on the genetics of type 2 diabetes published to date. We have assessed the role of two widely studied polymorphisms in the IL6 gene, using IPD from published and unpublished studies, and did not find evidence for an association between IL6 –573G>C and type 2 diabetes. In contrast, we determined that the GC and CC genotypes of IL6 –174G>C show an OR of 0.91 for association with type 2 diabetes, which corresponds to a risk reduction of nearly 9%. However, because the association between the IL6 –174G>C polymorphism and type 2 diabetes is borderline significant, a secondary analysis including additional data is critical. Thus, the present work represents a crucial first step toward elucidating the extent to which the IL6 –174G>C plays a role in type 2 diabetes susceptibility and provides additional evidence supporting a direct relationship between chronic subclinical inflammation and type 2 diabetes etiology.
Parts of this work were supported by the GSF Research Center for Environment and Health, the German Diabetes Center, the German Federal Ministry of Education, Science, Research and Technology (01 ER 9701/4)/National Genome Research Network, the German Federal Ministry of Health and Social Security, the Ministry of Science and Research of the State North Rhine-Westphalia, the German Research Foundation (HO1073/8-1, Wi621/12-1, and TH 784/2-1), the European Foundation for the Study of Diabetes, and the intramural research program of the U.S. National Institute of Diabetes and Digestive and Kidney Diseases. We gratefully acknowledge the participation in the original studies of all individuals used in this joint analysis. We thank Lu Qi for supplying summary data on the Nurses Health and the Health Professional Follow-up Studies. We thank Hubert Kolb for invaluable advice and discussions throughout the project and Andrea Schneider for excellent handling of a large part of the original data. We also thank Arsin Sabunchi and Anika Luze for expert technical support.
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication May 3, 2006 and accepted in revised form July 12, 2006
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||