help button home button Endocrine Society JCEM
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Journal of Clinical Endocrinology & Metabolism , doi:10.1210/jc.2004-1606
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wang, H.
Right arrow Articles by Elbein, S. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wang, H.
Right arrow Articles by Elbein, S. C.
Right arrowPubmed/NCBI databases
*Gene*GEO Profiles
*HomoloGene*OMIM
*UniGene
Medline Plus Health Information
*Diabetes
*Diabetic Kidney Problems
Related Collections
Right arrow Diabetes and Insulin
The Journal of Clinical Endocrinology & Metabolism Vol. 90, No. 2 1123-1129
Copyright © 2005 by The Endocrine Society

Molecular Screening and Association Analyses of the Interleukin 6 Receptor Gene Variants with Type 2 Diabetes, Diabetic Nephropathy, and Insulin Sensitivity

Hua Wang, Zhengxian Zhang, Winston Chu, Terri Hale, Judith J. Cooper and Steven C. Elbein

Division of Endocrinology and Metabolism (H.W., Z.Z., W.C., T.H., J.J.C., S.C.E.), Department of Internal Medicine, College of Medicine, University of Arkansas for Medical Sciences, Little Rock, Arkansas 72205; and Endocrinology Section (T.H., J.J.C., S.C.E.), Medicine and Research Services, Central Arkansas Veterans Healthcare System, Little Rock, Arkansas 72205

Address all correspondence and requests for reprints to: Steven C. Elbein, M.D., Division of Endocrinology and Metabolism, Central Arkansas Veterans Healthcare System, Endocrinology 111J-1/LR, Little Rock, Arkansas 72205. E-mail: elbeinstevenc{at}uams.edu.


    Abstract
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
IL-6 levels and polymorphisms have been implicated in type 2 diabetes mellitus (T2DM) and insulin resistance. The IL-6 receptor (IL-6R) comprises two subunits, IL-6R and gp130, of which IL-6R confers specificity to IL-6 action and is located in a region of replicated linkage to T2DM on chromosome 1q21. We screened this gene for variation in Northern European Caucasian and African-American ethnic groups. We identified 11 variants with a minor allele frequency over 5%, including two amino acid changes (D358A and V385I) and four variants in the 3' untranslated region. No variant was associated with obesity or measures of insulin sensitivity, but two single nucleotide polymorphisms in the 3' untranslated region showed a trend to an association with T2DM in all Caucasians, and three single nucleotide polymorphisms, including D358A, showed a trend (P < 0.06) to an association with T2DM among the subset of Northern European Caucasians. Variant V385I was unique to African-Americans and was significantly associated with diabetes and diabetic nephropathy (P < 0.05). Among individuals heterozygous for the four variants in the transcribed sequence, one allele was significantly overrepresented, thus suggesting the existence of a regulatory variant controlling mRNA stability or expression. IL-6R is not likely to explain the linkage to diabetes in this region, but our work supports a minor role of variants in T2DM risk and suggests that sequence variants may alter IL-6R mRNA levels and possibly levels of soluble IL-6R.


    Introduction
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
TYPE 2 DIABETES MELLITUS (T2DM) is characterized by insulin resistance and progressive pancreatic ß-cell failure. Increasing evidence suggests that chronic activation of the innate immune system reduces insulin sensitivity and hence may precede the development of T2DM (1, 2). Molecular markers of inflammation in T2DM have included the acute-phase response proteins sialic acid, {alpha}-1 acid glycoprotein, serum amyloid A, C-reactive protein, cortisol, and the cytokine IL-6 (3).

IL-6, a key mediator of the inflammatory response, is expressed widely, including expression in adipocytes. Although immune cells, fibroblasts, endothelial cells, and monocytes have been considered the major sources, adipocytes contribute significantly to circulating IL-6 (4, 5). Accordingly, IL-6 levels correlate with several measures of obesity, including body mass index (BMI), waist to hip ratio, and percent body fat (4), and are increased in T2DM and insulin-resistant states (6). Hence, elevated IL-6 levels and increased IL-6 action may contribute to the development of T2DM (7). In support of this hypothesis, polymorphisms of the IL-6 gene have been associated with T2DM in several studies (8, 9, 10). Furthermore, IL-6 levels and polymorphisms have been associated with the risk of diabetic nephropathy (11).

IL-6 acts through a receptor comprising two subunits, IL-6 receptor (IL-6R) and gp130. IL-6 first binds to the IL-6R, and this complex in turn binds subsequently to the ubiquitous gp130 subunit (12, 13). Specificity in the IL-6 response is conferred by the IL-6R subunit, which is expressed in its membrane-bound form in hepatocytes and leukocytes. Although most IL-6 action may be through membrane-bound IL-6R, membrane-bound IL-6R is cleaved proteolytically to a circulating, soluble form of IL-6R (sIL-6R), which may activate cells that express the gp130 subunit but not IL-6R (13). Increased levels of IL-6/sIL-6R have been reported in both T2DM (14) and diabetic nephropathy (15).

The IL-6R gene maps to human chromosome 1q21, a region of well-replicated linkage to T2DM described by our laboratory in both Northern European Caucasians (16, 17) and African-Americans (18) and in multiple other diverse populations (19). Thus, IL-6R must be considered a strong functional and positional candidate gene for T2DM. Based on these arguments, we tested the hypothesis that sequence variation in IL-6R was associated with T2DM, insulin resistance, and diabetic nephropathy in Caucasian and African-American populations. We have extensively screened the IL-6R gene for sequence variation, we have evaluated the role of these variants in T2DM and insulin sensitivity in Caucasians and in T2DM and diabetic nephropathy in African-Americans, and we have tested for possible regulatory polymorphisms by searching for an imbalance in IL-6R mRNA levels between parental alleles.


    Subjects and Methods
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
Experimental subjects

Caucasian individuals were ascertained in Utah for Northern European ancestry, as described previously (20). Additional Caucasian individuals and African-American individuals were ascertained from Arkansas for similar criteria. All diabetic individuals had at least one other first-degree relative with T2DM. Nondiabetic control individuals had no known family history of diabetes and had either a normal 75-g oral glucose tolerance test or a fasting glucose of less than 5.6 mmol/liter (100 mg/dl). Screening for mutations was conducted on DNA from 32 unrelated subjects with T2DM, consisting of 16 Caucasian individuals who were selected from families with the greatest evidence for linkage to chromosome 1q21-q24 and 16 African-American subjects from families ascertained in Arkansas for at least two diabetic siblings.

The association of variants with T2DM was assessed in both Caucasian and African-American populations. Studies of Caucasians were conducted in 192 unrelated, nondiabetic individuals and 192 unrelated individuals with T2DM requiring therapy, of whom 69 were drawn from families used in our linkage studies (16). Analyses of samples only ascertained in Utah included 191 cases and 118 control individuals. Association of single nucleotide polymorphisms (SNPs) among African-American individuals was tested in 129 unrelated nondiabetic individuals who had no family history of T2DM and 275 T2DM subjects who had at least one diabetic first-degree relative. The diabetic cases included 150 individuals who had diabetes and diabetic nephropathy (end-stage renal disease or a spot albumin to creatinine ratio > 300 mg/g creatinine) and 125 individuals with T2DM but no nephropathy (albumin to creatinine ratio < 30 mg/g creatinine). Alternatively, the association of IL-6R variants with diabetes and diabetic nephropathy in African-Americans was tested in samples of DNA pools constructed in triplicate from each of the three African-American subject groups (nondiabetic controls, individuals with T2DM and diabetic nephropathy, and individuals with T2DM but without nephropathy).

Insulin sensitivity was tested in 126 Caucasian nondiabetic members from 26 families for whom measurements of insulin sensitivity were available, as described previously (20). All subjects provided written informed consent under protocols approved by either the University of Utah or the University of Arkansas for Medical Sciences Institutional Review Boards.

Detection of sequence variation

IL-6R genomic sequence was obtained from the public database (accession no. NT_021933) and compared with the mRNA sequence (accession no. NM_000565) to determine the intron-exon boundaries. We designed 11 sets of primers to amplify each of the 10 exons, 100–200 bp of each flanking intron, and 1023 bp of the 5' flanking region using Primer 3.0 software (Whitehead Institute for Biomedical Research, Boston, MA). Each fragment of 480–580 bp was amplified in the presence of 32P-deoxycytidine triphosphate and digested with the appropriate restriction enzyme to generate fragments of 100–280 bp before single-strand conformation polymorphism analysis, as described previously (21). Conformation variants were characterized using bidirectional sequence analysis (22) on LI-COR GR-4200 sequencers (LI-COR Biotech, Lincoln, NE). Primer sequences and detailed methods are available from the authors.

Genetic analysis

Biallelic SNPs were genotyped by pyrosequencing on a PSQ-96 machine (Pyrosequencing, Inc., Uppsala, Sweden) according to the manufacturer methods but with the modification that each reaction was coamplified with a universal biotinylated primer and the complementary sequence was added to each allele-specific primer. SNP2 (rs4845617) could not be assayed accurately using pyrosequencing and was typed using an oligonucleotide ligation assay as described previously (23). In African-Americans, noncoding and synonymous SNPs were typed initially in pools using Pyrosequencing and Allele Quantification software (Pyrosequencing, Inc.). Each pool was constructed in triplicate and analyzed in duplicate. SNPs were typed in individual African-American samples if they altered the coding sequence, showed an association in Caucasians, showed a 5% difference in allele frequency estimated from pooled samples, or if pooled typing failed quality control checks. Finally, DNA from at least 92 unrelated individuals was tested for each ethnic group to determine linkage disequilibrium and Hardy Weinberg equilibrium. All reported SNPs were in Hardy Weinberg equilibrium.

Allele-specific expression

Differential expression of the two IL-6R alleles was examined in individuals heterozygous for SNPs located within the mRNA sequence. Total RNA (1 µg) from Epstein-Barr virus-transformed lymphocytes was reverse transcribed using random primers and Superscript II reverse transcriptase (Invitrogen Life Technologies, Carlsbad, CA). Allele-specific expression of SNPs 11, 12, 13, and 16 was assayed using pyrosequencing assays as described earlier, and expression of each allele was quantified using Allele Quantification software (Pyrosequencing, Inc.). The ratio of the alleles was compared with genomic DNA using the same assay. Primer sequences are available upon request.

Statistical analysis

The significance of differences in allele frequencies at SNPs was assessed by allelic association using Fisher’s exact test. The effect of each SNP on quantitative traits was evaluated using general linear regression models and tested in SPSS for Windows, version 12.0 software (SPSS, Inc., Chicago, IL). Skewed variables (BMI and insulin sensitivity) were natural logarithmically transformed to normality before analysis. All models included diagnosis, genotype, and gender as fixed factors. Family studies included age as a covariate and family membership as a random factor. Nominal significance was considered to be P < 0.05, and results are presented without correction for multiple tests given the correlations between SNPs. Linkage disequilibrium was assessed by allele counting and the expectation maximization algorithm, as implemented in the 2LD program (24) to assess D' and r2 values. Differences in haplotype distributions between cases and controls were tested using Arlequin version 2.0 (Genetics and Biometry Laboratory, Department of Anthropology, University of Geneva. Geneva, Switzerland) (25). Allele-specific expression was tested by comparing the major allele frequency estimated from Allele Quantification software (Pyrosequencing, Inc.) in the cDNA from transformed lymphocytes to that in at least 97 DNA samples from heterozygous individuals using a Mann-Whitney U test. Reported P values are based on the comparison of the ratio in cDNA to that in genomic DNA, but the reported allelic ratios for mRNA were corrected by dividing the mRNA ratio by the mean ratio estimated from at least 100 DNA samples from heterozygous individuals. Distortion of allelic ratios in mRNA was also tested by comparing the number of individuals with 50% increased expression between mRNA and DNA samples using Fisher’s exact test.


    Results
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
Sequence variation in the IL-6R gene

We screened the exons, adjacent intronic sequences, and 5' flanking region of the IL-6R gene spanning more than 60 kb of genomic DNA in both Caucasian and African-American individuals. We identified 17 SNPs and the previously examined microsatellite polymorphism (17), including two nonsynonymous SNPs (D358A, SNP11, or rs8192284; V385I, SNP12, or rs2228146) and two synonymous SNPs (SNP3, A31, or rs8192282; SNP4, H70, or rs2229237; Table 1Go; Fig. 1Go). In total, five SNPs were unique to African-American subjects, including V385I (SNP12), and two SNPs were observed only in Caucasian subjects. We identified one additional SNP from the public database (www.ncbi.nih.gov/SNP) in a region of intron 1 that was not screened (rs1386821; SNP17).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Sequence variation in the IL-6R gene

 


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 1. Map of IL-6R gene. Exons are shown as black boxes; the 5' and 3' untranslated regions are shown as gray boxes. Each SNP is shown in the approximate location, as detailed in Table 1Go.

 
We examined 10 SNPs and one simple tandem repeat polymorphism that were sufficiently polymorphic (minor allele frequency > 5%) in Caucasians to be informative (Table 2Go). Two SNPs in the 3' untranslated region showed a trend to an association with T2DM (SNP13, rs2229238, +59794 A>T; SNP16, rs4072391, +60776 C>T; Table 2Go), but the two SNPs were in strong linkage disequilibrium. When only subjects ascertained in Utah (196 cases and 118 controls) were considered, SNP17 (intron 1, +3943, rs1386822) was nominally significant (P = 0.02), and the only nonsynonymous SNP in Caucasians, D358A (SNP11; +48867, rs8192284), showed a trend to an association (P = 0.053), as did SNP7 (rs4845371; +30237 bp) in intron 5 (P = 0.057). No other SNP showed any association with T2DM in Caucasians. As reported previously (17), the microsatellite in intron 8 showed neither excess transmission to diabetic family members nor any association with T2DM among unrelated individuals (data not shown).


View this table:
[in this window]
[in a new window]
 
TABLE 2. Association of IL-6R sequence variants with T2DM in Caucasian and African-American populations

 
We typed seven SNPs in African-American individuals. We chose SNPs for individual genotyping if they were nonsynonymous, if they showed a trend to an association in Caucasians, or if the estimated allele frequencies differed by over 0.05 between pooled DNA from diabetic and control individuals (Table 2Go). Four SNPs were not tested due to low allele frequencies (SNP4, rs2229237; SNP8, and SNP14) or assay design problems (SNP10). V385I (SNP12) was significantly more common among subjects with T2DM than among controls [odds ratio (OR) = 1.62, P = 0.042; Table 2Go] and was most strongly associated in the subset of T2DM patients without diabetic nephropathy (frequency 0.25 for T2DM patients vs. 0.13 for controls, OR = 2.2, P = 0.0045). No other SNP, including D358A, showed any association with T2DM or with diabetic nephropathy among African-American case and control individuals (Table 2Go).

Linkage disequilibrium and haplotype analyses

Among Caucasians, SNPs from exon 2 (SNP3) to the 3' untranslated region (SNP16) were in strong linkage disequilibrium when measured by D' values (D' > 0.88; Table 3Go). We selected seven tag SNPs using r2 > 0.8 for the 10 SNPs that were over 10% minor allele frequency using the LDSelect program (26) (Table 3Go). Among Caucasians, only eight haplotypes comprising SNPs 2, 17, 3, 6, 7, 15, and 16 were predicted with a frequency of over 5%, thus accounting for 70% of all haplotypes (Table 4Go). Global comparison showed no significant difference in the distribution of the haplotypes between case and control populations (P = 0.35). Among African-American individuals, the block structure was similar with SNPs 3 through 12 showing strong linkage disequilibrium by D' (D' > 0.83), but no SNP pair exceeded the r2 threshold of 0.8 (Table 5Go). Only SNPs 5 and 7 were associated (r2 = 0.73). Although 10 tag SNPs are predicted based on the lower levels of linkage disequilibrium, our primary strategy in African-Americans was to test pooled DNA samples. Hence, we tested haplotypes comprising the seven SNPs that were typed in individual samples (Table 5Go). Again, we observed only seven of 128 possible haplotypes with a predicted frequency of more than 5%, and these seven haplotypes accounted for 72% of all predicted haplotypes (Table 6Go). In contrast to individual SNP data, the global haplotype distribution between cases and controls was not significantly different (P = 0.15; Table 6Go), although two rare haplotypes not shown in Table 5Go were significantly different between cases and controls (GTCAACT; frequency 0.041 vs. 0.009, P = 0.035; GGCAGCT; frequency 0.008 vs. 0.029, P = 0.03). Notably, the rare A allele of SNP12 (V385I) was present on only one common haplotype (GTCAACC; Table 6Go), where a trend to an association was noted (case frequency, 0.107; control frequency, 0.165), but the same allele was found on two rare haplotypes, including GTCAACT noted earlier, that were present in 3.3% and 1.3% of the full population. The presence of this allele on three haplotypes is consistent with low levels of linkage disequilibrium among the seven SNPs typed in African-American individuals (Table 5Go).


View this table:
[in this window]
[in a new window]
 
TABLE 3. Linkage disequilibrium in Caucasian subjects

 

View this table:
[in this window]
[in a new window]
 
TABLE 4. Haplotypes estimated from tag SNPs in Caucasians

 

View this table:
[in this window]
[in a new window]
 
TABLE 5. Linkage disequilibrium for tested SNPs in African-American subjects

 

View this table:
[in this window]
[in a new window]
 
TABLE 6. Estimated haplotype frequencies in African-Americans

 
IL-6R and quantitative traits

IL-6 levels are reported to correlate well with insulin resistance, and data support an association of IL-6 polymorphisms with insulin sensitivity. To address the hypothesis that IL-6R polymorphisms might alter insulin sensitivity, we tested SNPs 3, 6, 11, 16, and 17, which were both common and not in strong linkage disequilibrium, in 126 Caucasian individuals for whom insulin sensitivity measures were available. No SNP showed any trend to influence insulin sensitivity when age, gender, BMI, and family membership were included as covariates in a general linear model (P > 0.5 for all five SNPs). A comparable African-American population was not available for this study.

A previous study suggested that IL-6R variants might influence BMI (27). To test this hypothesis, we examined the following three populations: the 126 subjects who had undergone careful phenotyping, the Caucasian case control population, and our African-American case control population. Among the carefully phenotyped samples, we again found no trend for SNPs 3, 6, 11, or 16 to influence BMI (P > 0.5 for all analyses) when age, gender, glucose tolerance status, and family membership were included in the model. Similarly, among 360 Caucasian individuals with or without diabetes, we found no trend for SNPs 17, 3, 5, 6, 9, 11, 13, or 16 to influence BMI, neither when tested using a general linear model with gender, diagnosis, and IL-6R genotype as covariates nor when BMI was compared by genotype without considering covariates (data not shown). Among African-American subjects, SNP7 was associated with BMI whether genotypes were compared by ANOVA without consideration of diabetes status (P = 0.030) or when age, gender, and diagnosis were included in the model (P = 0.015). However, no other SNP, including V385I and D358A, showed any association with obesity (data not shown). Notably, most subjects with end-stage renal disease were not included in these analyses because accurate height and weight could not be obtained. Among those individuals for whom BMI was available, we would anticipate that advanced nephropathy would result in weight loss, but this was not included in our analysis.

Allele-specific expression

Several recent studies have demonstrated that differences in allele-specific expression for polymorphisms in the coding region of genes are common and inherited in transformed lymphocytes (28). SNPs in coding regions offer the possibility of comparing the ratio of the levels between two alleles among heterozygous individuals. Hence, we tested whether sequence variation in IL-6R might alter the mRNA ratio between the alleles by comparing the ratio of the two alleles for SNPs 11 (D358A), 12 (V1385I), 13 (3' untranslated region), and 16 (3' untranslated region) (Table 7Go). To correct for possible unequal amplification of the alleles, we used the same assay for both DNA and cDNA, and we corrected the raw measured ratio in cDNA using the ratio in genomic DNA. Significance was judged by comparing the uncorrected ratio in cDNA prepared from transformed lymphocyte to that observed in genomic DNA using a nonparametric Mann-Whitney U test. The results (Table 4Go) suggested that some but not all heterozygous individuals had increased expression of one allele by 50–250%, for mean levels that showed 30–60% excess expression from one allele.


View this table:
[in this window]
[in a new window]
 
TABLE 7. Allele-specific mRNA levels in transformed lymphocytes

 
Because of the large variation in expression and the possibility that the heterozygous coding region SNPs might be tagging a functional SNP at some distance, we analyzed the data to determine the number of samples in DNA and RNA showing extreme values of at a least 50% increase in one allele. Despite modestly increased mean ratios, the number of extreme values for expression ratios was markedly increased for transformed-lymphocyte mRNA for all four SNPs (P < 0.000001).


    Discussion
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 
Considerable data support a role for adipocyte-derived cytokines (adipokines), such as TNF-{alpha}, leptin, adiponectin, resistin, and IL-6, in insulin sensitivity and possibly T2DM. The metabolic pathways for these factors are thus potential candidates for diabetes pathogenesis. We examined the IL-6R based on considerable data supporting IL-6 in the pathogenesis of insulin resistance in obesity and T2DM (14, 29), the key role of IL-6R in IL-6 action, and the location of the IL-6R gene on chromosome 1q21 in a region of widely replicated linkage to T2DM. We identified 17 SNPs, including two that alter amino acids in the transmembrane domain. D358A is predicted to alter the splicing of IL-6R to sIL-6R and thus might alter the response to IL-6 in nonhepatic tissues (30). Despite the putative functional consequences, we only found a trend to an association of D358A among the subset of our Caucasian population that was ascertained in Utah for Northern European ancestry. Among all Caucasians, we could not show an association of this variant with T2DM. Similarly, SNPs 13 and 16, which were associated with unequal allelic mRNA levels, showed a trend to an association with T2DM among all Caucasians, but no other sequence variation in the IL-6R gene or haplotypes comprising tag SNPs were associated with T2DM among Caucasians. Additionally, no SNP was associated with measures of insulin sensitivity or obesity among our Caucasian population.

Our findings differ from the recent report of Wolford et al. (27) in Pima Indians, in which only two SNPs were identified from screening and only five SNPs were evaluated. These authors also reported no association of IL-6R SNPs with T2DM or prediabetic measures, but they reported an association of genotypes for all five SNPs, including D358A, with BMI. Escobar-Morreale et al. (31) examined the association of a microsatellite polymorphism of the IL-6R gene with obesity in Spanish women, but they did not test other variation in the gene. They reported an increased frequency of the 149-bp allele with obesity defined as a BMI of more than 25 kg/m2. Pedersen et al. (32) recently typed the single D358A variant in a large Danish cohort of 1393 subjects with T2DM and 4596 glucose-tolerant subjects. In contrast to Wolford et al. (27) and in support of our results in Caucasian Americans of Northern European ancestry, they found a modest association of D358A with T2DM (P = 0.02) but not with obesity or surrogate measures of insulin sensitivity (Pedersen, O., personal communication). We have 80% power to detect a 10% difference in allele frequencies for common variants, which translates into an OR of 1.5. Hence, the modest risk of the D358A allele (OR = 1.14) in the Danish study would likely not have reached statistical significance in our smaller study, and thus, our finding of a trend to an association is consistent with the Danish study.

Among African-American subjects, we identified a previously unreported amino acid change, V385I, that was associated with both T2DM and diabetic nephropathy. This variant, which was not seen in Caucasians in our study, not observed in Pima Indians in the study by Wolford et al. (27), and not reported among the seven SNPs detected in Koreans (33), is predicted to lie within the transmembrane domain of the IL-6R gene. However, the functional consequences of this conservative substitution are unknown. The apparent association of different alleles of V385I with diabetes and diabetic nephropathy was unexpected. The I385 variant was more common among diabetic individuals without nephropathy than among either controls or diabetic individuals with nephropathy, whereas the frequency was similar among control individuals and individuals with diabetic nephropathy. This finding may represent a spurious association. Alternatively, I385 may predispose to diabetes but protect against diabetic nephropathy, whereas the V385 allele may impart a modest protection against diabetes but increase the risk of nephropathy among those with diabetes and hyperglycemia. Given the complicated association and the modest P values, most of which would not be significant with a Bonferroni correction, this finding must be confirmed in other African-American populations. Nonetheless, a previous study supported a role for IL-6 polymorphisms in the progression of diabetic nephropathy (11). We also found a nominally significant association of SNP7 (intron 5) with BMI among African-American individuals. Again, considering the seven SNPs examined, the significance of this finding is uncertain.

In addition to the evidence for an association of V385I with diabetes and diabetic nephropathy in African-Americans and the trend to an association of SNPs 11 (D358A), 13, and 16 in Caucasians, we found evidence for possible cis-acting regulatory SNPs in the IL-6R gene through altered ratios of mRNA between two alleles from heterozygous individuals. We anticipated that either regulatory SNPs that altered mRNA expression or SNPs that stabilized the mRNA, such as 3' untranslated region SNPs, would alter the ratio between two alleles. Two scenarios for cis-acting SNPs are possible; there is either a direct effect from the SNP responsible for the difference or an indirect effect through partial linkage disequilibrium with an upstream regulatory SNP. The former scenario should give uniform ratios of one allele to the other, whereas incomplete linkage disequilibrium would result in some individuals with an altered ratio, and other individuals may be homozygous for the regulatory variant and would not show an altered mRNA ratio. Our data are most compatible with the latter scenario. Hence, modest increases in the mean ratio, in fact, reflect a large number of individuals in whom the mRNA ratio between alleles is 1.5 or greater. Furthermore, this observation applied to all four SNPs tested and to both ethnic groups. We propose that the altered allelic expression reflects an upstream polymorphism that altered IL-6R mRNA expression in individuals who carried the regulatory SNP. The apparently incomplete linkage disequilibrium between this undetected SNP and our transcribed SNPs suggests that the variant is likely in the 5' region of the gene, beyond the haplotype block and possibly outside of the region screened. Alternatively, the altered ratio might reflect epigenetic regulation, such as imprinting, although no imprinting in this region has been reported. Were the SNP in an unknown upstream regulatory region, low levels of linkage disequilibrium between the SNPs we tested and that putative regulatory SNP would likely mask any association with diabetes or metabolic traits.

In summary, we have identified a large number of SNPs in two populations in the IL-6R gene. Many of these SNPs were not identified in previous screens of the gene in Pima Indians or Koreans. We find the best evidence for an association of an amino acid change (V385I) with T2DM and diabetic nephropathy in African-Americans, and we find a trend to an association of several SNPs, including D358A, with T2DM in Caucasians. The latter finding corroborates unpublished data from Pedersen et al. (32). Finally, we demonstrate altered mRNA levels between the alleles of approximately 50% of individuals heterozygous for coding and 3' untranslated region SNPs, suggesting that, in addition to the amino acid variants (nonsynonymous SNPs), SNPs altering mRNA expression or perhaps stability are likely. Further studies are required to determine the identity of the regulatory SNP and to confirm the associations of D358A and V385I with T2DM.


    Acknowledgments
 
We thank Yiwen Jia and Xiaoqin Wang for RNA preparation and technical support, and the General Clinical Research Center nursing and laboratory staff for invaluable assistance.


    Footnotes
 
This work was supported by Grants DK039311 and DK054636 from the National Institutes of Health/National Institute of Diabetes and Digestive and Kidney Diseases, by the Research Service of the Department of Veterans Affairs, by grant support from the American Diabetes Association, and by General Clinical Research Center Grant M01RR14288 from National Center for Research Resources (National Institutes of Health) to the University of Arkansas for Medical Sciences.

First Published Online November 23, 2004

Abbreviations: BMI, Body mass index; IL-6R, IL-6 receptor; sIL-6R, soluble IL-6 receptor; SNP, single nucleotide polymorphism; T2DM, type 2 diabetes mellitus.

Received August 12, 2004.

Accepted November 10, 2004.


    References
 Top
 Abstract
 Introduction
 Subjects and Methods
 Results
 Discussion
 References
 

  1. Festa A, D'Agostino Jr R, Howard G, Mykkanen L, Tracy RP, Haffner SM 2000 Chronic subclinical inflammation as part of the insulin resistance syndrome: the Insulin Resistance Atherosclerosis Study (IRAS). Circulation 102:42–47[Abstract/Free Full Text]
  2. Vozarova B, Weyer C, Lindsay RS, Pratley RE, Bogardus C, Tataranni PA 2002 High white blood cell count is associated with a worsening of insulin sensitivity and predicts the development of type 2 diabetes. Diabetes 51:455–461[Abstract/Free Full Text]
  3. Pickup JC, Mattock MB, Chusney GD, Burt D 1997 NIDDM as a disease of the innate immune system: association of acute-phase reactants and interleukin-6 with metabolic syndrome X. Diabetologia 40:1286–1292[CrossRef][Medline]
  4. Fried SK, Bunkin DA, Greenberg AS 1998 Omental and subcutaneous adipose tissues of obese subjects release interleukin-6: depot difference and regulation by glucocorticoid. J Clin Endocrinol Metab 83:847–850[Abstract/Free Full Text]
  5. Mohamed-Ali V, Pinkney JH, Coppack SW 1998 Adipose tissue as an endocrine and paracrine organ. Int J Obes Relat Metab Disord 22:1145–1158[CrossRef][Medline]
  6. Vozarova B, Weyer C, Hanson K, Tataranni PA, Bogardus C, Pratley RE 2001 Circulating interleukin-6 in relation to adiposity, insulin action, and insulin secretion. Obes Res 9:414–417[Medline]
  7. Pradhan AD, Manson JE, Rifai N, Buring JE, Ridker PM 2001 C-reactive protein, interleukin 6, and risk of developing type 2 diabetes mellitus. JAMA 286:327–334[Abstract/Free Full Text]
  8. Vozarova B, Fernandez-Real JM, Knowler WC, Gallart L, Hanson RL, Gruber JD, Ricart W, Vendrell J, Richart C, Tataranni PA, Wolford JK 2003 The interleukin-6 (-174) G/C promoter polymorphism is associated with type-2 diabetes mellitus in Native Americans and Caucasians. Hum Genet 112:409–413[Medline]
  9. Fernandez-Real JM, Broch M, Vendrell J, Gutierrez C, Casamitjana R, Pugeat M, Richart C, Ricart W 2000 Interleukin-6 gene polymorphism and insulin sensitivity. Diabetes 49:517–520[Abstract]
  10. Kubaszek A, Pihlajamaki J, Punnonen K, Karhapaa P, Vauhkonen I, Laakso M 2003 The C-174G promoter polymorphism of the IL-6 gene affects energy expenditure and insulin sensitivity. Diabetes 52:558–561[Abstract/Free Full Text]
  11. Kitamura A, Hasegawa G, Obayashi H, Kamiuchi K, Ishii M, Yano M, Tanaka T, Yamaguchi M, Shigeta H, Ogata M, Nakamura N, Yoshikawa T 2002 Interleukin-6 polymorphism (-634C/G) in the promotor region and the progression of diabetic nephropathy in type 2 diabetes. Diabet Med 19:1000–1005[CrossRef][Medline]
  12. Kishimoto T, Tanaka T, Yoshida K, Akira S, Taga T 1995 Cytokine signal transduction through a homo- or heterodimer of gp130. Ann NY Acad Sci 766:224–234[Medline]
  13. Jones SA, Horiuchi S, Topley N, Yamamoto N, Fuller GM 2001 The soluble interleukin 6 receptor: mechanisms of production and implications in disease. FASEB J 15:43–58[Abstract/Free Full Text]
  14. Kado S, Nagase T, Nagata N 1999 Circulating levels of interleukin-6, its soluble receptor and interleukin-6/interleukin-6 receptor complexes in patients with type 2 diabetes mellitus. Acta Diabetol 36:67–72[CrossRef][Medline]
  15. Pecoits-Filho R, Barany P, Lindholm B, Heimburger O, Stenvinkel P 2002 Interleukin-6 is an independent predictor of mortality in patients starting dialysis treatment. Nephrol Dial Transplant 17:1684–1688[Abstract/Free Full Text]
  16. Elbein SC, Hoffman MD, Teng K, Leppert MF, Hasstedt SJ 1999 A genome-wide search for type 2 diabetes susceptibility genes in Utah Caucasians. Diabetes 48:1175–1182[Abstract]
  17. Das SK, Hasstedt SJ, Zhang Z, Elbein SC 2004 Linkage and association mapping of a chromosome 1q21–q24 type 2 diabetes susceptibility locus in Northern European Caucasians. Diabetes 53:492–499[Abstract/Free Full Text]
  18. Elbein SC, Zheng W, Wang X, Zhang Z, Hasstedt SJ 2004 Novel T2DM locus on chromosome 4 and confirmation of chromosome 1q21 T2DM linkage in multiplex African American families. Diabetes 53(Suppl 2):A270 (Abstract)
  19. McCarthy MI 2003 Growing evidence for diabetes susceptibility genes from genome scan data. Curr Diab Rep 3:159–167[Medline]
  20. Elbein SC, Wegner K, Kahn SE 2000 Reduced ß-cell compensation to the insulin resistance associated with obesity in members of Caucasian familial type 2 diabetic kindreds. Diabetes Care 23:221–227[Abstract]
  21. Wang H, Chu W, Hemphill C, Hasstedt SJ, Elbein SC 2002 Mutation screening and association of human retinoid X receptor {gamma} variation with lipid levels in familial type 2 diabetes. Mol Genet Metab 76:3–12
  22. Shevchenko YO, Bale SJ, Compton JG 2000 Mutation screening using automated bidirectional dideoxy fingerprinting. Biotechniques 28:134–138[Medline]
  23. Wang H, Chu W, Das SK, Ren Q, Hasstedt SJ, Elbein SC 2002 Liver pyruvate kinase polymorphisms are associated with type 2 diabetes in Northern European Caucasians. Diabetes 51:2861–2865[Abstract/Free Full Text]
  24. Ott J 1999 Analysis of human genetic linkage. 3rd ed. Baltimore: Johns Hopkins University Press
  25. University of Geneva 2000 Arlequin: a software for population genetics data analysis. Geneva: Genetics and Biometry Laboratory, Department of Anthropology, University of Geneva
  26. Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA 2004 Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 74:106–120[CrossRef][Medline]
  27. Wolford JK, Colligan PB, Gruber JD, Bogardus C 2003 Variants in the interleukin 6 receptor gene are associated with obesity in Pima Indians. Mol Genet Metab 80:338–343[CrossRef][Medline]
  28. Cheung VG, Conlin LK, Weber TM, Arcaro M, Jen KY, Morley M, Spielman RS 2003 Natural variation in human gene expression assessed in lymphoblastoid cells. Nat Genet 33:422–425[CrossRef][Medline]
  29. Kern PA, Ranganathan S, Li C, Wood L, Ranganathan G 2001 Adipose tissue tumor necrosis factor and interleukin-6 expression in human obesity-associated insulin resistance. Am J Physiol Endocrinol Metab 280:E745–E751
  30. Mullberg J, Oberthur W, Lottspeich F, Mehl E, Dittrich E, Graeve L, Heinrich PC, Rose-John S 1994 The soluble human IL-6 receptor. Mutational characterization of the proteolytic cleavage site. J Immunol 152:4958–4968[Abstract]
  31. Escobar-Morreale HF, Calvo RM, Villuendas G, Sancho J, San Millan JL 2003 Association of polymorphisms in the interleukin 6 receptor complex with obesity and hyperandrogenism. Obes Res 11:987–996[Medline]
  32. Hamid YH, Urhammer SA, Jensen DP, Glumer C, Borch-Johnsen K, Jorgensen T, Hansen T, Pedersen O 2004 Variation in the interleukin-6 receptor gene associates with type 2 diabetes in Danish whites. Diabetes 53:3342–3345[Abstract/Free Full Text]
  33. Kim LH, Lee HS, Kim YJ, Jung JH, Kim JY, Park BL, Shin HD 2003 Identification of novel SNPs in the interleukin 6 receptor gene (IL6R). Hum Mutat 21:450–451



This article has been cited by other articles:


Home page
Physiol. Rev.Home page
B. K. Pedersen and M. A. Febbraio
Muscle as an Endocrine Organ: Focus on Muscle-Derived Interleukin-6
Physiol Rev, October 1, 2008; 88(4): 1379 - 1406.
[Abstract] [Full Text] [PDF]


Home page
J. Am. Soc. Nephrol.Home page
J. F. Navarro-Gonzalez and C. Mora-Fernandez
The Role of Inflammatory Cytokines in Diabetic Nephropathy
J. Am. Soc. Nephrol., March 1, 2008; 19(3): 433 - 442.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
L. Qi, N. Rifai, and F. B. Hu
Interleukin-6 Receptor Gene Variations, Plasma Interleukin-6 Levels, and Type 2 Diabetes in U.S. Women
Diabetes, December 1, 2007; 56(12): 3075 - 3081.
[Abstract] [Full Text] [PDF]


Home page
Eur J EndocrinolHome page
M Bustamante, X Nogues, L Mellibovsky, L Agueda, S Jurado, E Caceres, J Blanch, R Carreras, A Diez-Perez, D Grinberg, et al.
Polymorphisms in the interleukin-6 receptor gene are associated with bone mineral density and body mass index in Spanish postmenopausal women
Eur. J. Endocrinol., November 1, 2007; 157(5): 677 - 684.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
W. S. Chu, S. K. Das, H. Wang, J. C. Chan, P. Deloukas, P. Froguel, L. J. Baier, W. Jia, M. I. McCarthy, M. C.Y. Ng, et al.
Activating Transcription Factor 6 (ATF6) Sequence Polymorphisms in Type 2 Diabetes and Pre-Diabetic Traits
Diabetes, March 1, 2007; 56(3): 856 - 862.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
S. C. Elbein, X. Wang, M. A. Karim, B. I. Freedman, D. W. Bowden, A. R. Shuldiner, F. L. Brancati, and W. H. L. Kao
Role of a Proline Insertion in the Insulin Promoter Factor 1 (IPF1) Gene in African Americans With Type 2 Diabetes.
Diabetes, October 1, 2006; 55(10): 2909 - 2914.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
S. K. Das, W. S. Chu, T. C. Hale, X. Wang, R. L. Craig, H. Wang, A. R. Shuldiner, P. Froguel, P. Deloukas, M. I. McCarthy, et al.
Polymorphisms in the Glucokinase-Associated, Dual-Specificity Phosphatase 12 (DUSP12) Gene Under Chromosome 1q21 Linkage Peak Are Associated With Type 2 Diabetes
Diabetes, September 1, 2006; 55(9): 2631 - 2639.
[Abstract] [Full Text] [PDF]


Home page
DiabetesHome page
O. P. Kristiansen and T. Mandrup-Poulsen
Interleukin-6 and Diabetes: The Good, the Bad, or the Indifferent?
Diabetes, December 1, 2005; 54(suppl_2): S114 - S124.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wang, H.
Right arrow Articles by Elbein, S. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wang, H.
Right arrow Articles by Elbein, S. C.
Right arrowPubmed/NCBI databases
*Gene*GEO Profiles
*HomoloGene*OMIM
*UniGene
Medline Plus Health Information
*Diabetes
*Diabetic Kidney Problems
Related Collections
Right arrow Diabetes and Insulin


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals