Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Variation in CHI3LI in Relation to Type 2 Diabetes and Related Quantitative Traits

  • Camilla Noelle Rathcke ,

    cnr@dadlnet.dk

    Affiliation Department of Endocrinology, Herlev Hospital, University of Copenhagen, Copenhagen, Denmark

  • Johan Holmkvist,

    Affiliation Hagedorn Research Institute, Gentofte, Denmark

  • Torben Jørgensen,

    Affiliation Research Centre for Prevention and Health, Glostrup Hospital, University of Copenhagen, Copenhagen, Denmark

  • Knut Borch-Johnsen,

    Affiliation Steno Diabetes Center, Gentofte, Denmark

  • Torben Hansen,

    Affiliation Hagedorn Research Institute, Gentofte, Denmark

  • Oluf Borbye Pedersen,

    Affiliations Hagedorn Research Institute, Gentofte, Denmark, Faculty of Health Sciences, University of Aarhus, Aarhus, Denmark, Faculty of Health Sciences, University of Copenhagen, Copenhagen, Denmark

  • Henrik Vestergaard

    Affiliations Department of Endocrinology, Herlev Hospital, University of Copenhagen, Copenhagen, Denmark, Faculty of Health Sciences, University of Copenhagen, Copenhagen, Denmark

Abstract

Background

CHI3LI encoding the inflammatory glycoprotein YKL-40 is located on chromosome 1q32.1. YKL-40 is involved in inflammatory processes and patients with Type 2 Diabetes (T2D) have elevated circulating YKL-40 levels which correlate with their level of insulin resistance. Interestingly, it has been reported that rs10399931 (−329 G/A) of CHI3LI contributes to the inter-individual plasma YKL-40 levels in patients with sarcoidosis, and that rs4950928 (−131 C/G) is a susceptibility polymorphism for asthma and a decline in lung function. We hypothesized that single nucleotide polymorphisms (SNPs) or haplotypes thereof the CHI3LI locus might influence risk of T2D. The aim of the present study was to investigate the putative association between SNPs and haplotype blocks of CHI3LI and T2D and T2D related quantitative traits.

Methods/Principal Findings

Eleven SNPs of CHI3LI were genotyped in 6514 individuals from the Inter99 cohort and 2924 individuals from the outpatient clinic at Steno Diabetes Center. In cas-control studies a total of 2345 T2D patients and 5302 individuals with a normal glucose tolerance test were examined.

We found no association between rs10399931 (OR, 0.98 (CI, 0.88–1.10), p = 0.76), rs4950928 (0.98 (0.87–1.10), p = 0.68) or any of the other SNPs with T2D. Similarly, we found no significant association between any of the 11 tgSNPs and T2D related quantitative traits, all p>0.14. None of the identified haplotype blocks of CHI3LI showed any association with T2D, all p>0.16.

Conclusions/Significance

None of the examined SNPs or haplotype blocks of CHI3LI showed any association with T2D or T2D related quantitative traits. Estimates of insulin resistance and dysregulated glucose homeostasis in T2D do not seem to be accounted for by the examined variations of CHI3LI.

Introduction

The primary molecular defects in Type 2 diabetes (T2D) remain largely unknown, but it is widely recognized that both genetic and environmental risk factors have interacting roles in the pathogenesis and the progression of the disease. However, recent discoveries have led to a progress in the understanding of the molecular genetic background of T2D brought forward by genome-wide association (GWA) studies which have dramatically increased the number of validated type 2 diabetes loci with modest impact on relative diabetes risk [1][3]. Several loci of interest may be located in the region of human chromosone 1q21–24 which shows replicated linkage to T2D in at least eight populations including North European Caucasian families, Pima Indians and Old Order Amish [1]. Furthermore, linkage analysis has also provided evidence for linkage of the metabolic syndrome (central obesity, dyslipidemia, hyperglycemia and hypertension) to 1q23–31 [4]. These regions as well as other closely related regions of chromosome 1 are under intense search for identification of T2D susceptibility genes.

Chitinase 3-like 1 (CHI3LI) encoding the inflammatory glycoprotein YKL-40 is located in a phylogenetically highly conserved area on chromosome 1q32.1. YKL-40 is involved in inflammatory processes and extracellular matrix remodelling including angiogenesis and functions as a growth factor for several cell types [5]. Studies have shown, that patients with T2D have elevated plasma YKL-40 levels [6], [7] which correlate with insulin resistance and with their diabetic lipid profile, but not with estimates of glycemic regulation [6]. Also patients with type 1 diabetes (T1D) have elevated YKL-40 levels which increase with increasing levels of albuminuria [8].

The CHI3LI promotor variant −329 G/A (rs10399931) associates with the inter-individual plasma YKL-40 levels in patients with sarcoidosis [9]. Furthermore, the functional CHI3LI promotor SNP −131C→G (rs4950928) increases the mRNA transcript levels coding for the YKL-40 protein in a genome wide study of gene expression in Epstein-Barr virus-transformed lymphoblastoid cells from children with asthma [10] and is shown to be predictive of asthma, bronchial hyperresponsiveness and reduced lung function [11]. Since CHI3LI is a susceptibility gene for asthma with elevated circulating YKL-40 levels as a biomarker, it is possible that one of these single nucleotide polymorphisms (SNPs) or another gene variation with which they could be in tight linkage disequilibrium with may be functional in relation to T2D.

Applying a HapMap–based tgSNP approach the objective of the present study was to investigate the putative association of common variation in the CHI3LI locus with T2D, obesity and several T2D associated quantitative traits.

Methods

Ethics Statement

Informed written consent was obtained from all participants before participation. The study was approved by the Ethical Committees of Copenhagen and Aarhus and was in accordance with the principle of the Helsinki Declaration II.

Study population

A total of 11 tgSNPs (including the rs10399931 and rs4950928) covering all haplotype blocks in CHI3LI locus was genotyped. Genotyping was performed in 9438 individuals including Inter99, a population-based sample of unrelated, middle-aged Danish individuals living in the greater Copenhagen area, and 2924 individuals recruited from the outpatient clinic at Steno Diabetes Center.

The Inter99 study population has been described in detail previously [12]. In brief, both participants from the Inter99 cohort and individuals recruited from the outpatient clinic were characterized by an oral glucose tolerance test (OGTT). Impaired glucose homeostasis (IFG, IGT and diabetes) was defined according to WHO criteria after a 75 g OGTT [13] giving a classification of normal glucose tolerance (NGT; n = 5302), impaired fasting glycemia (IFG; n = 540), impaired glucose tolerance (IGT; n = 830) or screen-detected untreated T2D (SDU DM; n = 294). Missing data was seen in 421 individuals (4.5%) and a total of 2091 individuals had known T2D. The participants were investigated for an association between genotype and quantitative metabolic traits; individuals with known type 2 diabetes were excluded from these analyses. All participants were of Danish nationality.

Biochemical and anthropometric measurements

Height and body weight were measured in light indoor clothes and without shoes, and BMI was calculated as weight (kg) divided by (height (m))2. Waist circumference was measured in the standing position midway between the iliac crest and the lower costal margin and hip circumference at its maximum.

In both Inter99- and Steno Diabetes Center participants, blood samples were drawn in the fasting state after a 12-h overnight fast. Plasma glucose was analysed by a glucose oxidase method (Granutest, Merck, Darmstadt, Germany). HbA1C was measured by ion-exchange high performance liquid chromatography (normal reference range: 4.1–6.4%) and serum insulin (excluding des(31, 32) and intact proinsulin) was measured using the AutoDELFIA insulin kit (Perkin-Elmer/Wallac, Turku, Finland). Serum C-peptide concentrations were measured by a time-resolved fluoroimmunoassay (AutoDELFIA C-peptide kit; Perkin-Elmer/Wallac, Turku, Finland).

Serum triglyceride and total and HDL-cholesterol were analysed using enzymatic colorimetric methods (GPO-PAP and CHOD-PAP, Roche Molecular Biochemicals, Mannheim, Germany). HOMA-IR was calculated as fasting plasma glucose [mmol/l] multiplied by fasting serum insulin [pmol/l] and divided by 22.5. Insulinogenic index at 30 min was calculated as fasting serum insulin subtracted from plasma insulin at 30 min [pmol/l] divided by plasma glucose at 30 min [mmol/l].

Genotyping

SNPs in a region 22 kb upstream and 10 kb downstream of CHI3LI were chosen from the HapMap project (www.hapmap.org). Haplotype tgSNPs were selected from the HapMap data using TAGGER [14]. TAGGER was used with a 5% MAF cut off, agressive tagging (r2>0.8). Two haplotype blocks with 2 and 5 SNPs, respectively, and with similar structures in both patients with T2D and in controls were identified (Figure 1). The other tgSNPs were not included in the haplotype block analyses since they were not part of any block. Genotyping was performed using TaqMan allelic discrimination (Kbioscience, Herts, UK) with a success rate >97.0%. Discordance was 0% as judged from re-genotyping of 965 random duplicate samples. Genotype distribution obeyed Hardy Weinberg Equilibrium (HWE), all p>0.25.

thumbnail
Figure 1. Linkage disequilibrium (LD) pattern of the investigated variations of CHI3L1.

https://doi.org/10.1371/journal.pone.0005469.g001

Statistical analyses

Logistic regression analyses were used to analyse differences in allele and genotype frequencies between case and control subjects; analyses were adjusted for age, sex and BMI where appropriate. The statistical analyses were performed using RGui version 2.5.0 (available at http://www.r-project.org). The haplotype block analyses were analysed for an association with T2D after adjustment for age, sex and BMI. Haplotype analyses were performed using PLINK [15]. P-values are not adjusted for multiple testing.

Results

Clinical characteristics for the Inter99 cohort and the individuals recruited at Steno Diabetes Center are presented in Table 1. A total of 2345 T2D patients and 5302 individuals with a normal glucose tolerance test were identified.

Overall we found no association of rs10399931, rs4950928 or any of the other tgSNPs of CHI3LI with T2D after adjustment for age, sex, BMI and glucose tolerance, all p>0.25 (Table 2). In similar analyses, none of the tgSNPs were associated with impaired fasting glucose (IFG) or impaired glucose tolerance (IGT), although a borderline protective association with IGT was detected for rs12123883 (Table S1 and S2).

thumbnail
Table 2. Case control association studies of 2345 T2D patients and 5302 glucose tolerant control participants in relation to the 11 tgSNPs of CHI3LI.

https://doi.org/10.1371/journal.pone.0005469.t002

Table 3 and 4 show the multiadjusted association of rs10399931 and rs4950928 with BMI and other T2D related quantitative traits among NGT, IFG, IGT and screen-detected untreated T2D in Inter99 and Steno Diabetes Center collected individuals. We found a significant association between rs10399931 and the serum C-peptide concentration after 2 hours, p = 0.03, but in all we found no significant association between rs10399931 or rs4950928 and any of the T2D quantitative traits. Similarly, none of the other SNPs showed any associations (data not shown).

thumbnail
Table 3. Association of rs10399931 of CHI3LI and T2D related quantitative traits among NGT, IFG, IGT and screen-detected untreated T2D in Inter99 and Steno Diabetes Center collected individuals.

https://doi.org/10.1371/journal.pone.0005469.t003

thumbnail
Table 4. Association of rs49509028 of CHI3LI and T2D related quantitative traits among NGT, IFG, IGT and screen-detected untreated T2D in Inter99 and Steno Diabetes Center collected individuals.

https://doi.org/10.1371/journal.pone.0005469.t004

The two identified haplotype blocks of CHI3LI did not show any association with T2D (Table 5).

thumbnail
Table 5. CHI3LI haplotype association studies in relation to T2D.

https://doi.org/10.1371/journal.pone.0005469.t005

Discussion

This is a large-scale and to our knowledge the first study of variation in CHI3LI in relation to T2D and pre-diabetic quantitative traits. CHI3LI encodes the inflammatory protein YKL-40 which previously has been found to be elevated in plasma from patients with T2D [6], [7] and from patients with T1D where it also correlates with increasing levels of albuminuria [8]. In our study, we have exmained 11 tgSNPs covering the common variation (MAF>5%) in CHI3LI including the SNPs rs10399931 and rs4950928 which previously have been found to contribute to interindividual variations in circulating YKL-40 levels in patients with sarcoidosis [9] and asthma [11], respectively. Furthermore, CHI3LI has been found to be a susceptibility gene for asthma, since rs4950928 is predictive of asthma, bronchial hyperresponsiveness and reduced lung function.

A total of 9438 individuals have been examined, but no significant association of rs10399931, rs4950928 or any of the other tgSNPs of CHI3LI were found with T2D or any of the intermediate T2D related quantitative traits. A borderline protective association with IGT was detected for rs12123883, but rs12123883 is a rare SNP with only few homozygotes and the protective association is most likely a chance finding. Furthermore, the p-value would not withstand correction for multiple testing considering the number of tests performed. On the same line, none of the two identified haplotype blocks showed any association with T2D.

In association analyses of rs10399931 and SNPs analysed in the DGI genome wide association scan study, rs10399931 is in LD with rs4950929 (r2 = 0.74) and rs946263 (r2 = 0.78). Other correlations had r2-values<0.05. In the DGI genome wide association scan study, no significant associations of rs4950929 with any of the type 2 diabetic phenotypic traits were found, but the rs946263 G-allele was nominally found to be associated with T2D (p = 0.027) [16]. rs946263, which is a synonymous SNP located −9639 base pairs from the translational start site in CHI3LI, is also in perfect LD (r2 = 1.0) with rs4950928, which accounts for approximately 10% of the variance in plasma YKL-40 levels in asthma patients [11]. Beside the LD of rs10399931 and rs4950928 with rs946263, none of the examined SNPs or haplotype blocks of CHI3LI seem to be of any significance for the T2D diagnosis or any of the T2D related quantitative traits. Since the encoded inflammatory protein, YKL-40, of CHI3LI seems to be both a marker of inflammation and endothelial dysfunction [5], [17] it is possible that it is the endothelial dysfunction and the later micro- and macrovascular complications of T2D that accounts for the elevated YKL-40 levels in T2D patients. In endothelial dysfunction, elevated YKL-40 levels seem to play a role in relation to cell migration, reorganization and tissue remodelling as a response to endothelial damage [18][20].

The possible assocation between YKL-40 and vascular damage is supported by the finding of an independent association between elevated circulating YKL-40 levels and increasing levels of albuminuria in T1D patients [8]. Elevated plasma YKL-40 is also found to be associated with the presence and extent of coronary artery disease as assessed by coronary angiography [21] and just recently, YKL-40 levels have been found to be elevated in patients with myocardial infarction [22]. Furthermore, increasing YKL-40 levels do predict cardiovascular mortality in individuals without known DM or CHD after adjustment for known CV risk factors and markers [7].

Substantial evidence supports the notion that endothelial dysfunction and low-grade inflammation play a major pathophysiological role in the development of insulin resistance and the progression to manifest T2D [23][29]. We have previously shown, that YKL-40 levels are elevated in patients with T2D where it correlates with insulin resistance as assessed by the HOMA model and with levels of triglycerides and non-esterified fatty acids [6]. Since we found no association of any of the tgSNPs and surrogate markers of insulin resistance or the diabetic lipid profile in the present study, it does not seem to be any of these tgSNPs that explain the previously found correlation with these variables [6].

Multiple genes on different chromosomes may influence inflammatory biomarker levels and may have a potential role in the development of insulin resistance through a low-grade inflammatory state [30][33]. In proximity of the YKL-40 encoding gene, locus 1q32.2 has been documented to influence the plasma levels of ICAM-1[30], an adhesion molecule well-documented to participate in endothelial dysfunction and the development of T2D [23], [25]. However, association studies between YKL-40 and ICAM-1, either at protein or transcript levels, have not been conducted. Neither have association studies of variations of CHI3L1 in relation to micro- or macrovascular conditions never been made.

In conclusion, both rs10399931 and rs4950928 are in high LD with rs946263 which is reported to associate with phenotypic traits of T2D in DGI [16]. However, neither rs1039991 nor rs4950928 nor any other tgSNPs of CHI3LI which were investigated in the present study are associated with T2D or T2D related quantitative traits including estimates of insulin resistance and dysregulated glucose homeostasis in the examined approximately 9500 Danes. Further studies are required to determine YKL-40 levels as a marker of inflammation in the prediabetic state and the role of CHI3LI in the pathogenesis of T2D.

Supporting Information

Table S1.

Case control association studies of 540 individuals with impaired fasting glucose (IFG) and 5302 glucose tolerant control participants in relation to the 11 tgSNPs of CHI3LI.

https://doi.org/10.1371/journal.pone.0005469.s001

(0.08 MB DOC)

Table S2.

Case control association studies of 830 individuals with impaired glucose tolerance (IGT) and 5302 glucose tolerant control participants in relation to the 11 tgSNPs of CHI3LI

https://doi.org/10.1371/journal.pone.0005469.s002

(0.08 MB DOC)

Acknowledgments

We wish to thank A. Forman, I.-L. Wantzin and M. Stendal for technical assistance and Michel Kristensen, Arne Nielsen and Grete Lademann for management and secretarial support.

Author Contributions

Conceived and designed the experiments: CNR TJ KBJ TH OBP HV. Performed the experiments: CNR TJ KBJ TH OBP HV. Analyzed the data: CNR JH. Contributed reagents/materials/analysis tools: CNR JH TH HV. Wrote the paper: CNR. Critical revision of the manuscript: JH TJ KBJ TH OBP HV.

References

  1. 1. McCarthy MI (2003) Growing evidence for diabetes susceptibility genes from genome scan data. Curr Diab Rep 3: 159–67.
  2. 2. McCarthy MI, Zeggini E (2006) Genetics of type 2 diabetes. Curr Diab Rep 6: 147–54.
  3. 3. McCarthy MI, Zeggini E (2007) Genome-wide association scans for Type 2 diabetes: new insights into biology and therapy. Trends Pharmacol Sci 28: 598–601.
  4. 4. Langefeld CD, Wagenknecht LE, Rotter JI, Williams AH, Hokanson JE, et al. (2004) Linkage of the metabolic syndrome to 1q23–q31 in Hispanic families: the Insulin Resistance Atherosclerosis Study Family Study. Diabetes 53: 1170–4.
  5. 5. Rathcke CN, Vestergaard H (2006) YKL-40, a new inflammatory marker with relation to insulin resistance and with a role in endothelial dysfunction and atherosclerosis. Inflamm Res 55: 221–7.
  6. 6. Rathcke CN, Johansen JS, Vestergaard H (2006) YKL-40, a biomarker of inflammation, is elevated in patients with type 2 diabetes and is related to insulin resistance. Inflamm Res 55: 53–9.
  7. 7. Rathcke CN, Raymond I, Kistorp C, Hildebrandt P, Faber J, et al. (2009) Low grade inflammation as measured by levels of YKL-40: Association with an increased overall and cardiovascular mortality rate in an elderly population. Int J Cardiol.. In Press.
  8. 8. Rathcke CN, Persson F, Tarnow L, Rossing P, Vestergaard H (2009) YKL-40, a marker of inflammation and endothelial dysfunction, is elevated in patients with Type 1 Diabetes and increases with levels of albuminuria. Diabetes Care.. In Press.
  9. 9. Kruit A, Grutters JC, Ruven HJ, van Moorsel CC, van den Bosch JM (2007) A CHI3L1 gene polymorphism is associated with serum levels of YKL-40, a novel sarcoidosis marker. Respir Med 101: 1563–71.
  10. 10. Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, et al. (2007) A genome-wide association study of global gene expression. Nat Genet 39: 1202–7.
  11. 11. Ober C, Tan Z, Sun Y, Possicj JD, Pan L, et al. (2008) Effect of variation in CHI3L1 on serum YKL-40 level, risk of asthma, and lung function. N Engl J Med 358: 1682–91.
  12. 12. Jorgensen T, Borch-Johnsen K, Thomsen TF, Ibsen H, Glumer C, et al. (2003) A randomized non-pharmacological intervention study for prevention of ischaemic heart disease: baseline results Inter99. Eur J Cardiovasc Prev Rehabil 10: 377–86.
  13. 13. Alberti KG, Zimmet PZ (1998) Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med 15: 539–53.
  14. 14. de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, et al. (2005) Efficiency and power in genetic association studies. Nat Genet 37: 1217–23.
  15. 15. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–75.
  16. 16. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, et al. (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316: 1331–6.
  17. 17. Johansen JS (2006) Studies on serum YKL-40 as a biomarker in diseases with inflammation, tissue remodelling, fibroses and cancer. Dan Med Bull 53: 172–209.
  18. 18. Shackelton LM, Mann DM, Millis AJ (1995) Identification of a 38-kDa heparin-binding glycoprotein (gp38k) in differentiating vascular smooth muscle cells as a member of a group of proteins associated with tissue remodeling. J Biol Chem 270: 13076–83.
  19. 19. Malinda KM, Ponce L, Kleinman HK, Shackelton LM, Millis AJ (1999) Gp38k, a protein synthesized by vascular smooth muscle cells, stimulates directional migration of human umbilical vein endothelial cells. Exp Cell Res 250: 168–73.
  20. 20. Millis AJ, Hoyle M, Kent L (1986) In vitro expression of a 38,000 dalton heparin-binding glycoprotein by morphologically differentiated smooth muscle cells. J Cell Physiol 127: 366–72.
  21. 21. Kucur M, Isman FK, Karadag B, Vural VA, Tavsanoglu S (2007) Serum YKL-40 levels in patients with coronary artery disease. Coron Artery Dis 18: 391–6.
  22. 22. Nojgaard C, Host NB, Christensen IJ, Poulsen SH, Egstrup K, et al. (2008) Serum levels of YKL-40 increases in patients with acute myocardial infarction. Coron Artery Dis 19: 257–63.
  23. 23. Meigs JB, Hu FB, Rifai N, Manson JE (2004) Biomarkers of endothelial dysfunction and risk of type 2 diabetes mellitus. JAMA 291: 1978–86.
  24. 24. Pickup JC (2004) Inflammation and activated innate immunity in the pathogenesis of type 2 diabetes. Diabetes Care 27: 813–23.
  25. 25. Kolb H, Mandrup-Poulsen T (2005) An immune origin of type 2 diabetes? Diabetologia 48: 1038–50.
  26. 26. Eckel RH, Grundy SM, Zimmet PZ (2005) The metabolic syndrome. Lancet 365: 1415–28.
  27. 27. Festa A, Hanley AJ, Tracy RP, D'Agostino R Jr, Haffner SM (2003) Inflammation in the prediabetic state is related to increased insulin resistance rather than decreased insulin secretion. Circulation 108: 1822–30.
  28. 28. Rutter MK, Meigs JB, Sullivan LM, D'Agostino RB Sr, Wilson PW (2004) C-reactive protein, the metabolic syndrome, and prediction of cardiovascular events in the Framingham Offspring Study. Circulation 110: 380–5.
  29. 29. Ingelsson E, Hulthe J, Lind L (2008) Inflammatory markers in relation to insulin resistance and the metabolic syndrome. Eur J Clin Invest 38: 502–9.
  30. 30. Dupuis J, Larson MG, Vasan RS, Massaro JM, Wilson PW, et al. (2005) Genome scan of systemic biomarkers of vascular inflammation in the Framingham Heart Study: evidence for susceptibility loci on 1q. Atherosclerosis 182: 307–14.
  31. 31. Pantsulaia I, Trofimov S, Kobyliansky E, Livshits G (2002) Genetic and environmental influences on IL-6 and TNF-alpha plasma levels in apparently healthy general population. Cytokine 19: 138–46.
  32. 32. de Maat MP, Bladbjerg EM, Hjelmborg JB, Bathum L, Jespersen J, et al. (2004) Genetic influence on inflammation variables in the elderly. Arterioscler Thromb Vasc Biol 24: 2168–73.
  33. 33. Vickers MA, Green FR, Terry C, Mayosi BM, Julier C, et al. (2002) Genotype at a promoter polymorphism of the interleukin-6 gene is associated with baseline levels of plasma C-reactive protein. Cardiovasc Res 53: 1029–34.