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

Non-Cholesterol Sterol Levels Predict Hyperglycemia and Conversion to Type 2 Diabetes in Finnish Men

  • Henna Cederberg,

    Affiliation Department of Medicine, University of Eastern Finland and Kuopio University Hospital, Kuopio, Finland

  • Helena Gylling,

    Affiliations Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland, Department of Medicine, Division of Internal Medicine, University of Helsinki, Helsinki, Finland

  • Tatu A. Miettinen,

    Affiliation Department of Medicine, Division of Internal Medicine, University of Helsinki, Helsinki, Finland

  • Jussi Paananen,

    Affiliation Department of Medicine, University of Eastern Finland, Kuopio, Finland

  • Jagadish Vangipurapu,

    Affiliation Department of Medicine, University of Eastern Finland, Kuopio, Finland

  • Jussi Pihlajamäki,

    Affiliation Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland

  • Teemu Kuulasmaa,

    Affiliation Department of Medicine, University of Eastern Finland, Kuopio, Finland

  • Alena Stančáková,

    Affiliation Department of Medicine, University of Eastern Finland, Kuopio, Finland

  • Ulf Smith,

    Affiliation Department of Molecular and Clinical Medicine, Sahlgrenska Academy at the University of Gothenburg, Gothenburg, Sweden

  • Johanna Kuusisto,

    Affiliation Department of Medicine, University of Eastern Finland and Kuopio University Hospital, Kuopio, Finland

  • Markku Laakso

    markku.laakso@kuh.fi

    Affiliation Department of Medicine, University of Eastern Finland and Kuopio University Hospital, Kuopio, Finland

Abstract

We investigated the levels of non-cholesterol sterols as predictors for the development of hyperglycemia (an increase in the glucose area under the curve in an oral glucose tolerance test) and incident type 2 diabetes in a 5-year follow-up study of a population-based cohort of Finnish men (METSIM Study, N = 1,050) having non-cholesterol sterols measured at baseline. Additionally we determined the association of 538,265 single nucleotide polymorphisms (SNP) with non-cholesterol sterol levels in a cross-sectional cohort of non-diabetic offspring of type 2 diabetes (the Kuopio cohort of the EUGENE2 Study, N = 273). We found that in a cross-sectional METSIM Study the levels of sterols indicating cholesterol absorption were reduced as a function of increasing fasting glucose levels, whereas the levels of sterols indicating cholesterol synthesis were increased as a function of increasing 2-hour glucose levels. A cholesterol synthesis marker desmosterol significantly predicted an increase, and two absorption markers (campesterol and avenasterol) a decrease in the risk of hyperglycemia and incident type 2 diabetes in a 5-year follow-up of the METSIM cohort, mainly attributable to insulin sensitivity. A SNP of ABCG8 was associated with fasting plasma glucose levels in a cross-sectional study but did not predict hyperglycemia or incident type 2 diabetes. In conclusion, the levels of some, but not all non-cholesterol sterols are markers of the worsening of hyperglycemia and type 2 diabetes.

Introduction

Type 2 diabetes is associated with atherogenic dyslipidemia [1] increasing the risk of cardiovascular mortality and morbidity [2]. Abnormalities in cholesterol metabolism in type 2 diabetes include enhanced cholesterol synthesis and reduced cholesterol absorption [3][6]. Also insulin resistance and obesity have been previously associated with elevated cholesterol synthesis and low cholesterol absorption in cross-sectional studies [7][10]. However, the mechanisms underlying the association of non-cholesterol sterol levels with hyperglycemia have remained unclear. Atherogenic dyslipidemia predicts the development of type 2 diabetes [1], and therefore sterols associated with cholesterol synthesis and absorption could also contribute to the risk of the development of hyperglycemia and type 2 diabetes.

Single nucleotide polymorphisms (SNPs) in the sitosterolemia genes ABCG5 and ABCG8, encoding the ATP-binding cassette sub-family G member proteins, have been shown to regulate cholesterol absorption [11][13]. These transporters actively efflux plant sterols and, to a lesser extent cholesterol, back into the intestinal lumen. In a previous genome-wide association study SNPs rs4245791 and rs41360247 in/near ABCG8 were associated with serum phytosterol levels [14]. The impact of SNPs in ABCG5/8, and possible other, yet undiscovered SNPs on the association between non-cholesterol sterol levels and glucose metabolism is largely unknown.

To investigate in more detail the association of non-cholesterol sterols with the risk of the worsening of hyperglycemia and incident type 2 diabetes we addressed 1) the relationship of cholesterol synthesis markers (squalene, cholestenol, lathosterol, and desmosterol) and absorption markers (campesterol, sitosterol, avenasterol, and cholestanol) with fasting and 2-hour glucose levels in a cross-sectional setting, 2) the association of non-cholesterol sterol levels as predictors for the worsening of hyperglycemia and the development of incident type 2 diabetes in a 5-year follow-up of the METSIM (METabolic Syndrome In Men) Study, and 3) the role of insulin sensitivity, insulin secretion and genetic factors in this relationship.

Subjects and Methods

Participants and Clinical Measurements

METSIM study.

Non-cholesterol sterol levels were determined in a random sample of 1,050 consequently examined men from the population-based METSIM Study including 10,197 men. This study has been described in detail elsewhere [15]. Glucose tolerance was evaluated according to the American Diabetes Association criteria [16]. Participants with previously diagnosed diabetes mellitus or receiving statin or ezetimibe medication were excluded, leaving 746 men for statistical analysis. A total of 474 METSIM participants (64%) have been re-examined during an ongoing follow-up study [follow-up time 5.4±0.3 years, mean±SD); 25.1% with normal glucose tolerance, 51.3% with isolated impaired fasting (IFG) glucose, 3.6% with isolated impaired glucose tolerance (IGT), 12.1% with IFG and IGT, and 7.9% with incident type 2 diabetes]. Body mass index (BMI) was calculated as weight (kg) divided by the height (m) squared. Ethical Approval The study was approved by the Ethics Committee of the University of Kuopio and Kuopio University Hospital, conducted in accordance with the principles of the Declaration of Helsinki, and all study participants gave written informed consent.

Kuopio cohort of the EUGENE2 study.

The participants were healthy, non-diabetic offspring of patients with type 2 diabetes, as previously described in detail [17]. The probands were randomly selected among type 2 diabetic patients living in the region of Kuopio, Finland. One of the parents had to have type 2 diabetes and the other parent normal glucose tolerance in an OGTT. Altogether 273 offspring (128 males, 145 females) were included in the current study aiming to identify SNPs regulating sterol levels. Ethical Approval The study was approved by the Ethics Committee of the University of Kuopio and Kuopio University Hospital, conducted in accordance with the principles of the Declaration of Helsinki, and all study participants gave written informed consent.

Sterol and Lipid Measurements

The methods to measure sterol and lipid levels have been described in detail elsewhere [7]. Briefly, serum cholesterol, cholesterol precursors (squalene, cholestenol, lathosterol, and desmosterol), campesterol, sitosterol and avenasterol (plant sterols), and cholestanol, a metabolite of cholesterol, were quantified from non-saponifiable serum material by capillary gas chromatography (Agilent 6890N Network GC System, Agilent Technologies, Wilmington, DE) equipped with a 50 m long Ultra 2 capillary column (5% Phenyl-methyl siloxane) (Agilent Technologies, Wilmington, DE) [18]. α-cholestane was used as the internal standard. The inter-assay coefficients of variation for sterols at relevant concentrations were as follows: cholesterol 3.2%, cholestanol 2.7%, desmosterol 6.0%, lathosterol 3.7%, campesterol 1.8%, sitosterol 2.4%, and avenasterol 4.1%, respectively. External quality control samples were run and analyzed with every 25 study samples. The serum values were expressed as 102 µmol/mmol of cholesterol (called ratio in the text) by dividing non-cholesterol sterol concentrations with the cholesterol value of the same GC run in order to eliminate the changing concentrations of sterol transporters, mainly LDL. The ratios to cholesterol of serum cholesterol precursors reflect whole-body cholesterol synthesis and those of plant sterols and cholestanol cholesterol absorption.

Insulin Sensitivity and Insulin Secretion Indices

Matsuda index of insulin sensitivity (ISI) was calculated as previously reported [15]. An index of early-phase insulin secretion during an OGTT, InsAUC0–30/GluAUC0–30, was calculated as (fasting insulin +30 min insulin)/(fasting glucose +30 min glucose) (pmol/mmol) [19].

Genetic Analyses

Genotyping.

Genotyping for the participants of the EUGENE2 Study was carried out at the Finnish Genome Center in Helsinki, Finland. Genotyping was performed with the commercial release of the Infinium HumanHap 550 k version 3 SNP microarrays (Illumina, San Diego, CA). Briefly, 750 ng of DNA was used in genotyping according to the manufacturer’s protocol (Illumina). Genotyping for the METSIM Study participants was performed using the TaqMan Allelic Discrimination Assays (Applied Biosystems) or the Sequenom MassARRAY platform at the University of Eastern Finland.

Genotyping quality assessment and control in the EUGENE2 study.

Software package PLINK v1.07 [20] was used for quality assessment and control for the genotyping data. Gender calls from X chromosome genotype data was verified to be in concordance with the reported gender of each individual. To verify known familial relationships and to detect not reported first-degree cryptic relationships, pairwise identity-by-descent analysis was performed. Individuals with genotyping call rates <99% were excluded. All of the 273 samples passed this criterion. Total genotyping call rate in remaining individuals was ≥99.2%. Genotypes advanced to the actual association analysis if they passed the following quality control criteria 1) had a >95% genotype call rate (4,201 markers failed), 2) had a minor allelic frequency (MAF) >0.1% (18,720 markers failed) and 3) demonstrated Hardy-Weinberg Equilibrium (HWE) with a P>1×10−5 (575 markers failed). In addition, 317 heterozygous haploid genotypes were set to be missing. Out of the 561,466 genotypes, a total of 538,265 markers passed these quality control criteria. For each genotype, summary information about MAF, HWE and call rates were recorded, and this information was used when evaluating marker quality for replication.

Genotype replication.

Most significant SNPs from the genome-wide association analysis were replicated in the METSIM cohort. SNPs were selected for replication based on statistical significance, effect size, linkage disequilibrium structure and allele frequencies. Priority was given to genotypes that showed strong associations with several sterol measurements. Total of 177 SNPs were genotyped. Genotypes advanced to the association analysis if they passed the following quality control criteria 1) had a >99% genotype call rate (20 markers failed), 2) had a MAF >5% (30 markers failed), and 3) demonstrated HWE with a P>0.001 (13 markers failed). Out of the 177 genotypes, a total of 134 markers passed these quality control criteria.

Genome-wide association analysis.

Software package GenABEL [21] was used for genome-wide association analysis of the EUGENE2 data. Sterol measurements were transformed to log10-scale to adjust for skewed distributions. To account for the family structures in the data, Identity-by-state kinship matrix between the subjects was calculated and then used within mixed model score test for association in related people [22]. Gender specific differences in the levels of non-cholesterol levels have been previously reported [23] and higher levels of cholesterol synthesis markers squalene, lathosterol and desmosterol were also observed in men as compared women in the EUGENE2 Kuopio cohort (data not shown). Therefore, gender, together with age and BMI were included in the analysis as covariates and the additive genetic model was applied.

Replication analysis.

PLINK was used to analyze the data from the METSIM replication study. Sterol measurements were transformed to log10-scale to adjust for skewed distributions. Linear regression including age and BMI as a covariate and assuming an additive genetic model was used.

Meta-analysis.

Joint fixed-effects meta-analysis combining results from both Stage 1 (genome-wide association analysis) and Stage 2 (replication study) was performed using PLINK. Minor alleles from Stage 1 analysis were used as reference alleles.

Replication of ABCG8 rs4299376 in the entire METSIM cohort.

Genotyping of rs4299376 of ABCG8 in the entire METSIM cohort was performed with TaqMan Allelic Discrimination Assay (Allied Biosystems). Call rate was 100% and the error rate 0%.

Statistical Analysis

Statistical analyses were conducted using SPSS version 19 (SPSS, Chicago, IL). All sterol to cholesterol ratios, glucose, BMI, insulin sensitivity index (Matsuda ISI) and insulin secretion index (InsAUC0–30/GluAUC0–30) were log-transformed to correct for their skewed distribution. Sterols were compared across the fasting and 2-hour glucose categories using the general linear model. P≤0.006 (corrected for 8 tests by Bonferroni method) was considered statistically significant. Correlations between sterol levels and indices of insulin sensitivity and insulin secretion were evaluated with partial correlations. Associations between baseline sterol levels and glucose area under the curve (AUC) at follow-up were evaluated using linear regression analysis. Associations between baseline sterol levels and incident type 2 diabetes were evaluated with Cox regression analysis (participants having incident type 2 diabetes at follow-up vs. participants having NGT/NFG and HbA1c <6.5% at baseline and follow-up examinations, N = 97). Incident type 2 diabetes at follow-up was defined by the American Diabetes Association criteria [24] as follows: fasting plasma glucose (FPG) ≥7.0 mmol/l or 2-h plasma glucose (2hPG) ≥11.1 mmol/l or HbA1c ≥6.5%, or by a physician-based diagnosis between baseline and follow-up examinations and the use of anti-diabetic medication. For statistical analyses of non-cholesterol sterols as predictors for hyperglycemia or incident type 2 diabetes, P value adjusted for multiple comparisons (P≤0.006 corrected for 8 traits) by Bonferroni method for considered as statistically significant. For GWAS analysis, P value adjusted for multiple comparisons (P<6×10−9 corrected for 8 traits) by Bonferroni method was considered to be statistically significant. For genetic association analysis, unstandardized effect sizes [B (SE)] per copy of the risk allele were estimated by linear regression analysis adjusted for age, BMI, smoking and physical activity, using untransformed dependent variables. Variables not normally distributed were log-transformed for statistical analyses.

Results

Associations of Sterols with Glucose Levels, Insulin Sensitivity and Insulin Secretion

After the exclusion of participants with previously diagnosed diabetes mellitus and participants on statin medication, a total of 746 men were included in statistical analyses. Their mean age was 57.6 years (SD 5.8), BMI 27.1 kg/m2 (SD 3.9), 18.8% of them were current smokers, and 63.8% were physically active (undertaking regular physical exercise at least once a week for ≥30 min). Mean FPG at baseline was 5.6 mmol/l (SD 0.6) and 2hPG 6.1 mmol/l (SD 2.20). Levels of sterols at baseline are given in Table S1.

As shown in Figure 1 the levels of three absorption markers, campesterol, sitosterol and avenasterol were significantly reduced (P<1×105) as a function of increasing FPG levels, but no significant changes were observed for cholesterol synthesis markers. Desmosterol significantly (P<1×10−5) and lathosterol and cholestenol nominally (P<0.01) increased as a function of increasing 2hPG levels, whereas the levels of absorption markers decreased (campesterol, P<0.01).

thumbnail
Figure 1. Non-cholesterol sterols indicating cholesterol synthesis (diamonds) and absorption (circles) in relation to hyperglycemia, according to the categories of fasting (FPG) and 2-hour plasma glucose (2hPG) levels (N = 746).

Data points represent unadjusted means of sterol levels in each glucose category. P values were calculated using general linear model adjusted for age, BMI, smoking, and physical activity: **P<0.01, ***P<10−5. Individuals with previously diagnosed diabetes and on statin or ezetimibe treatment were excluded.

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

The associations of sterol levels with insulin secretion and insulin sensitivity are shown in Table 1. Synthesis markers, of which desmosterol particularly, correlated inversely with Matsuda ISI, and positively with insulin secretion. In contrast, all markers of cholesterol absorption correlated positively with Matsuda ISI, and cholestanol correlated inversely with insulin secretion. Partial correlations of synthesis markers and cholestanol with insulin secretion lost their statistical significance after the adjustment for Matsuda, suggesting that low insulin sensitivity is the primary mechanism explaining the association of sterols indicating cholesterol synthesis with glucose levels. Similarly, the correlations of the three absorption markers, campesterol, sitosterol and avenasterol were also more significant with insulin sensitivity than with insulin secretion.

thumbnail
Table 1. Partial correlations of sterol levels indicating cholesterol synthesis and cholesterol absorption with insulin sensitivity and insulin secretion adjusted for age, BMI, smoking and physical activity.

https://doi.org/10.1371/journal.pone.0067406.t001

The correlations of sterol levels with FPG and 2hPG were considerably weaker than those with insulin sensitivity as shown in Table 2. Synthesis markers did not correlate significantly with FPG levels but campesterol, sitosterol and avenasterol had significant inverse correlations with FPG. These correlations considerably weakened after the adjustment for both insulin sensitivity and insulin secretion, whereas no major changes were observed when adjusted for insulin secretion only. Results for 2hPG were quite similar with the exception of higher correlations of synthesis markers with 2hPG than with FPG. Adjustment for all confounding factors and additionally for insulin sensitivity and insulin secretion abolished all statistically significant associations. Taken together, these results provide evidence that the association of glucose metabolism with sterols is mainly mediated by low insulin sensitivity.

thumbnail
Table 2. Partial correlations between sterol levels and fasting plasma glucose (FPG), and 2-hour plasma glucose (2hPG), adjusted for age, BMI, smoking, physical activity, and additionally for insulin sensitivity and insulin secretion.

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

Effects of Sterol Levels on the Worsening of Hyperglycemia and the Development of Incident Type 2 Diabetes

We first evaluated in a 5-year follow-up of the METSIM Study the effects of baseline sterol levels as predictors of the worsening of hyperglycemia, defined by a glucose area under the curve (Glucose AUC) in an OGTT (measurement of glucose levels at 0, 30 and 120 min) in 474 individuals without type 2 diabetes at baseline (Table 3). The levels of synthesis markers, cholestenol (P = 6.5×10−4) and desmosterol (P = 1.6×10−4) significantly predicted (P≤0.006) an increase in Glucose AUC after the adjustment for age, BMI, smoking and physical activity. In contrast, absorption markers campesterol, sitosterol and avenasterol significantly predicted a decrease in Glucose AUC at follow-up, after the adjustment for age, BMI, smoking and physical activity. The associations were largely attributable to low insulin sensitivity since additional adjustment for Matsuda ISI weakened all associations (Table 3), and abolished the association of desmosterol and campesterol with Glucose AUC at follow-up. Additional adjustment for insulin secretion, glucose AUC at baseline and total triglycerides did not substantially alter the results.

thumbnail
Table 3. Association between baseline sterol levels and glucose area under the curve (AUC) at the 5-year follow-up study.

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

A cholesterol synthesis marker desmosterol significantly predicted an increase, and absorption markers campesterol and avenasterol predicted a decrease in the risk of incident type 2 diabetes (Table 4). Associations of desmosterol, campesterol and avenasterol with incident type 2 diabetes were abolished after the adjustment for age, BMI, smoking and physical activity.

thumbnail
Table 4. Association between baseline sterol levels and newly developed type 2 diabetes in a 5-year follow-up study.

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

SNPs Regulating Sterol Levels

We performed a genome-wide association study of 538,265 SNPs and non-cholesterol sterol levels in the Kuopio cohort of the EUGENE2 study, and a targeted replication of 134 SNPs in the METSIM cohort. Table 5 presents genome-wide significant and Figure 2 nominally significant results (P<0.001) of a meta-analysis of the discovery sample and the replication sample. The three genome-wide significant SNPs were all associated with the cholesterol absorption marker levels, rs4299376 (from P = 4.7×10−15 to P = 4.3×10−19) and rs6544713 (P = 2.5×10−15 to 3.7×10−19, linkage disequilibrium in European populations with rs429976, r2 = 0.99) of ABCG8, and rs6756629 of ABCG5/ABCG8 (from P = 4.8×10−12 to P = 1.8×10−17). Effect sizes, calculated as per minor allele, were positive for the ABCG8 SNPs and negative for rs6756629 of ABCG5/ABCG8. Genome-wide significant SNPs (rs4299376, rs6544713, and rs6756629) were all associated with the cholesterol absorption marker levels even after the exclusion of participants with previously diagnosed diabetes and those taking statin medication, and adjustment for age, gender and BMI (Table 5). The major (sterol-decreasing) allele of the most significant SNP, rs4299376, was associated with elevated levels of FPG (P = 0.022) in the additive model, but not with the 2hPG (P = 0.425) or Glucose AUC (P = 0.144) in the cross-sectional METSIM Study cohort (5,639 men without previously diagnosed diabetes and without statin treatment). ABCG8 rs4299376 did not predict the worsening of hyperglycemia in a 5-year follow-up study (P≥0.121 for association with FPG, 2hPG and Glucose AUC at follow-up in the additive model).

thumbnail
Figure 2. SNPs associated with sterol levels (left panel) and the effect size evaluated by the standardized β coefficient from linear regression model (right panel).

Genotyping results represent a meta-analysis of 1st stage (GWAS, EUGENE2 Study) and 2nd stage of targeted replication (METSIM Study), total N = 1,013. SNPs having an association P<0.001 are shown. Individuals with previously diagnosed diabetes and on statin or ezetimibe treatment were excluded. P-values are adjusted for age, gender and BMI.

https://doi.org/10.1371/journal.pone.0067406.g002

thumbnail
Table 5. Genotyping results: 1st stage results obtained from genome-wide association study (EUGENE2 Study), 2nd stage results obtained from replication study (METSIM Study), and meta-analysis of both cohorts.

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

Discussion

We demonstrated for the first time that non-cholesterol sterol levels predicted the worsening of hyperglycemia and the development of incident type 2 diabetes in a 5-year prospective follow-up of a Finnish population-based study including 1,050 men (the METSIM cohort). We also demonstrated that rs4299376 and rs6544713 of ABCG8 and rs6756629 of ABCG5/8 were significantly associated with all cholesterol absorption markers (P for association from 4.8×10−12 to 3.7×10−19), but not with cholesterol synthesis markers. The major allele of the most significant SNP, rs4299376 of ABCG8, was nominally associated with elevated levels of FPG at baseline, but did not predict the development of hyperglycemia in a 5-year follow-up study.

Elevated cholesterol synthesis and reduced absorption have been reported in individuals with type 2 diabetes or hyperglycemia 36,89, and have been associated with insulin resistance [7] in cross-sectional studies. However, none of the previous studies has investigated the association of sterol levels in the entire range of hyperglycemia. We showed that in a cross-sectional setting, cholesterol synthesis markers were elevated and cholesterol absorption markers decreased as a function of 2hPG and FPG levels, respectively.

In our 5-year follow-up of the METSIM cohort we found that the cholesterol synthesis markers cholestenol and desmosterol significantly (P≤0.006) predicted increases in Glucose AUC, and desmosterol also incident type 2 diabetes. Adjustment for Matsuda ISI alone largely abolished these associations demonstrating that low insulin sensitivity was the major mechanism explaining the worsening of hyperglycemia and the development of incident diabetes. Absorption markers campesterol, sitosterol, and avenasterol significantly predicted a decrease in Glucose AUC, and campesterol and avenasterol also a lower risk of incident type 2 diabetes. Statistically significant associations were abolished when adjusted for Matsuda ISI and other confounding factors. To the best of our knowledge, this is the first prospective study demonstrating that non-cholesterol sterols predict the worsening of hyperglycemia and incident type 2 diabetes.

We also investigated the role of genetic factors in the association of sterol levels with the risk of type 2 diabetes. We showed that rs4299376 and rs6544713 (in high linkage disequilibrium with rs4299376) of ABCG8 and rs6756629 of ABCG5/8 were very significantly associated with plant sterols (campesterol, sitosterol and avenasterol) and cholestanol (β from 0.038 to 0.087, P<1×10−14 for rs4299376, and β from −0.046 to −0.108, P<1×10−11 for rs6756629) after adjustment for age, gender and BMI. The ATP hemitransporter proteins ABCG5 and ABCG8 are key regulators of dietary cholesterol absorption in the intestinal lumen. SNP rs4299376 of ABCG8 has previously been associated with LDL cholesterol level and rs6756629 of ABCG5/8 with total cholesterol level [25], [26]. Previously published GWAS identified two SNPs of ABCG8 (rs41360247 and rs4245791) as regulators of plant sterol (campesterol and sitosterol) levels [14]. Our results are in agreement with these findings as SNPs identified in our study are in close linkage equilibrium (r2 = 0.99–1.0) with the SNPs reported earlier [14]. However, the major allele of rs4299376 of ABCG8, which was very significantly associated with surrogate markers of cholesterol absorption, was nominally associated with elevated fasting glucose in 5,639 men in cross-sectional analysis of the METSIM cohort (P = 0.022), but not with the worsening of hyperglycemia or the development of incident type 2 diabetes in our follow-up study, probably due to a relatively small sample size in our follow-up study.

Our study has limitations. Only Finnish men were included in the METSIM Study, and therefore, we do not know whether our results are applicable to women and to different ethnic or racial groups. We had only modest power to demonstrate statistically significant associations of gene variants with non-cholesterol sterol levels.

In conclusion, our population-based study is the first to show that desmosterol, marker of cholesterol synthesis, and campesterol and avenasterol, markers of cholesterol absorption, significantly predict the worsening of hyperglycemia and incident type 2 diabetes. The predominant mechanism linking sterol metabolism to hyperglycemia is likely to be low insulin sensitivity.

Supporting Information

Table S1.

Mean level and the range of sterol levels at baseline (N = 746).

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

(DOCX)

Author Contributions

Conceived and designed the experiments: ML HG J. Pihlajamäki. Performed the experiments: TAM TK. Analyzed the data: HC J. Paananen JV AS. Contributed reagents/materials/analysis tools: TAM. Wrote the paper: HC ML. Collected the data: US JK ML Designed the study: ML HG. Reviewed/edited and approved the final version of the manuscript: HC HG TAM J. Paananen JV J. Pihlajamäki TK AS US JK ML.

References

  1. 1. Adiels M, Olofsson SO, Taskinen MR, Borén J (2006) Diabetic dyslipidaemia. Curr Opin Lipidol 17: 238–246.
  2. 2. Laakso M, Cederberg H (2012) Glucose control in diabetes: which target level to aim for? J Intern Med 272: 1–12.
  3. 3. Sutherland WH, Scott RS, Lintott CJ, Robertson MC, Stapely SA, et al. (1992) Plasma non-cholesterol sterols in patients with non-insulin dependent diabetes mellitus. Horm Metab Res 24: 172–175.
  4. 4. Simonen PP, Gylling HK, Miettinen TA (2002) Diabetes contributes to cholesterol metabolism regardless of obesity. Diabetes Care 25: 1511–1515.
  5. 5. Simonen P, Gylling H, Miettinen TA (2008) The validity of serum squalene and non-cholesterol sterols as surrogate markers of cholesterol synthesis and absorption in type 2 diabetes. Atherosclerosis 197: 883–888.
  6. 6. Briones ER, Steiger DL, Palumbo PJ, O’Fallon WM, Langworthy AL, et al. (1986) Sterol excretion and cholesterol absorption in diabetics and nondiabetics with and without hyperlipidemia. Am J Clin Nutr 44: 353–361.
  7. 7. Gylling H, Hallikainen M, Pihlajamäki J, Simonen P, Kuusisto J, et al. (2010) Insulin sensitivity regulates cholesterol metabolism to a greater extent than obesity: lessons from the METSIM Study. 51: 2422–2427.
  8. 8. Stranberg TE, Salomaa V, Vanhanen H, Miettinen TA (1996) Associations of fasting blood glucose with cholesterol absorption and synthesis in nondiabetic middle-aged men. Diabetes 45: 755–761.
  9. 9. Pihlajamäki J, Gylling H, Miettinen TA, Laakso M (2004) Insulin resistance is associated with increased cholesterol synthesis and decreased cholesterol absorption in normoglycemic men. J Lipid Res 45: 507–512.
  10. 10. Miettinen TA, Gylling H (2000) Cholesterol absorption efficiency and sterol metabolism in obesity. Atherosclerosis 153: 241–248.
  11. 11. Gylling H, Hallikainen M, Pihlajamäki J, Agren J, Laakso M, et al. (2004) Polymorphisms in the ABCG5 and ABCG8 genes associate with cholesterol absoption and insulin sensitivity. J Lipid Res 45: 1660–1665.
  12. 12. Berge KE, von Bergmann K, Lutjohann D, Guerra R, Grundy SM, et al. (2002) Heritability of plasma noncholesterol sterols and relationship to DNA sequence polymorphism in ABCG5 and ABCG8. J Lipid Res 43: 486–494.
  13. 13. Yu L, Li-Hawkins J, Hammer RE, Berge KE, Horton JD, et al. (2002) Overexpression of ABCG5 and ABCG8 promotes biliary cholesterol secretion and reduces fractional absorption of dietary cholesterol. J Clin Invest 110: 671–680.
  14. 14. Teupser D, Baber R, Ceglarek U, Scholz M, Illig T, et al. (2010) Genetic regulation of serum phytosterol levels and risk of coronary artery disease. Circ Cardiovasc Genet 3: 331–339.
  15. 15. Stancáková A, Javorský M, Kuulasmaa T, Haffner SM, Kuusisto J, et al. (2009) Changes in insulin sensitivity and insulin release in relation to glycemia and glucose tolerance in 6,414 Finnish men. Diabetes 58: 1212–1221.
  16. 16. Genuth S, Alberti KG, Bennett P, Buse J, Defronzo R, et al. (2003) Follow-up report on the diagnosis of diabetes mellitus. Diabetes Care 26: 3160–3167.
  17. 17. Laakso M, Zilinskaite J, Hansen T, Boesgaard TW, Vänttinen M, et al. (2008) Insulin sensitivity, insulin release and glucagon-like peptide-1 levels in persons with impaired fasting glucose and/or impaired glucose tolerance in the EUGENE2 study. Diabetologia 51: 502–511.
  18. 18. Miettinen TA (1988) Cholesterol metabolism during ketoconazole treatment in man. J Lipid Res 29: 43–51.
  19. 19. Abdul-Ghani MA, Matsuda M, Jani R, Jenkinson CP, Coletta DK, et al. (2008) The relationship between fasting hyperglycemia and insulin secretion in subjects with normal or impaired glucose tolerance. Am J Physiol Endocrinol Metab 295: E401–E406.
  20. 20. 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–575.
  21. 21. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM (2007) GenABEL: an R library for genome-wide association analysis. Bioinformatics 23: 1294–1296.
  22. 22. Chen WM, Abecasis GR (2007) Family-based association tests for genomewide association scans. Am J Hum Genet 81: 913–926.
  23. 23. Matthan NR, Zhu L, Pencina M, D'Agostino RB, Schaefer EJ, et al. (2013) Sex-specific differences in the predictive value of cholesterol homeostasis markers and 10-year cardiovascular disease event rate in framingham offspring study participants. JAHA 2: e005066.
  24. 24. International Expert Committee (2009) International Expert Committee report on the role of the A1C assay in the diagnosis of diabetes. Diabetes Care 32: 1327–1334.
  25. 25. Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, et al. (2009) Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet 41: 56–65.
  26. 26. Lu Y, Feskens EJ, Boer JM, Imholz S, Verschuren WM, et al. (2010) Exploring genetic determinants of plasma total cholesterol levels and their predictive value in a longitudinal study. Atherosclerosis 213: 200–205.