Skip to main content
  • More from ADA
    • Diabetes
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care
  • Subscribe
  • Log in
  • My Cart
  • Follow ada on Twitter
  • RSS
  • Visit ada on Facebook
Diabetes Care

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • Special Article Collections
    • ADA Standards of Medical Care
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • Special Article Collections
    • ADA Standards of Medical Care
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes Care
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • Special Article Collections
    • ADA Standards of Medical Care
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • Special Article Collections
    • ADA Standards of Medical Care
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Cardiovascular and Metabolic Risk

Epigenetic Link Between Statin Therapy and Type 2 Diabetes

  1. Carolina Ochoa-Rosales1,2⇑,
  2. Eliana Portilla-Fernandez1,
  3. Jana Nano3,4,
  4. Rory Wilson3,5,
  5. Benjamin Lehne6,
  6. Pashupati P. Mishra7,
  7. Xu Gao8,9,
  8. Mohsen Ghanbari1,10,
  9. Oscar L. Rueda-Ochoa1,11,
  10. Diana Juvinao-Quintero12,
  11. Marie Loh6,13,
  12. Weihua Zhang6,14,
  13. Jaspal S. Kooner14,15,16,17,
  14. Hans J. Grabe18,
  15. Stephan B. Felix19,20,
  16. Ben Schöttker8,21,
  17. Yan Zhang8,
  18. Christian Gieger3,5,
  19. Martina Müller-Nurasyid22,23,24,
  20. Margit Heier3,25,
  21. Annette Peters3,5,
  22. Terho Lehtimäki7,
  23. Alexander Teumer19,26,
  24. Hermann Brenner8,21,
  25. Melanie Waldenberger3,5,
  26. M. Arfan Ikram1,
  27. Joyce B.J. van Meurs27,
  28. Oscar H. Franco28,
  29. Trudy Voortman1,
  30. John Chambers6,13,14,16,17,
  31. Bruno H. Stricker1 and
  32. Taulant Muka28⇑
  1. 1Department of Epidemiology, Erasmus MC University Medical Center, Rotterdam, the Netherlands
  2. 2Centro de Vida Saludable de la Universidad de Concepción, Concepción, Chile
  3. 3Institute of Epidemiology, Helmholtz Zentrum München, Neuherberg, Germany
  4. 4German Center for Diabetes Research, Neuherberg, Germany
  5. 5Research Unit Molecular Epidemiology, Helmholtz Zentrum München, Neuherberg, Germany
  6. 6Department of Epidemiology and Biostatistics, Imperial College London, London, U.K.
  7. 7Department of Clinical Chemistry, Fimlab Laboratories, and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
  8. 8Division of Clinical Epidemiology and Aging Research, German Cancer Research Center (DKFZ), Heidelberg, Germany
  9. 9Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY
  10. 10Department of Genetics, School of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
  11. 11Electrocardiography Research group, School of Medicine, Universidad Industrial de Santander, Bucaramanga, Colombia
  12. 12MRC Integrative Epidemiology Unit, Bristol Medical School, University of Bristol, U.K.
  13. 13Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore, Singapore
  14. 14Department of Cardiology, Ealing Hospital, London North West University Healthcare NHS Trust, Middlesex, U.K.
  15. 15National Heart and Lung Institute, Imperial College London, London, U.K.
  16. 16Imperial College Healthcare NHS Trust, London, U.K.
  17. 17MRC-PHE Centre for Environment and Health, Imperial College London, London, U.K.
  18. 18Department of Psychiatry and Psychotherapy, University Medicine Greifswald, Greifswald, Germany
  19. 19Partner Site Greifswald, German Center for Cardiovascular Research (DZHK), Greifswald, Germany
  20. 20Department of Internal Medicine B, University Medicine Greifswald, Greifswald, Germany
  21. 21Network Aging Research, University of Heidelberg, Heidelberg, Germany
  22. 22Institute for Medical Information Processing, Biometry, and Epidemiology (IBE), Faculty of Medicine, Hospital of the Ludwig-Maximilians-University (LMU) Munich, Munich, Germany
  23. 23Department of Internal Medicine I (Cardiology), Hospital of the Ludwig-Maximilians-University (LMU) Munich, Munich, Germany
  24. 24Institute of Genetic Epidemiology, Helmholtz Zentrum München - German Research Center for Environmental Health, Neuherberg, Germany
  25. 25KORA Study Centre, University Hospital of Augsburg, Augsburg, Germany
  26. 26Institute for Community Medicine, University Medicine Greifswald, Greifswald, Germany
  27. 27Department of Internal Medicine, Erasmus MC University Medical Center, Rotterdam, the Netherlands
  28. 28Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
  1. Corresponding authors: Taulant Muka, taulant.muka{at}ispm.unibe.ch, and Carolina Ochoa-Rosales, c.ochoarosales{at}erasmusmc.nl
Diabetes Care 2020 Apr; 43(4): 875-884. https://doi.org/10.2337/dc19-1828
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

OBJECTIVE To investigate the role of epigenetics in statins’ diabetogenic effect comparing DNA methylation (DNAm) between statin users and nonusers in an epigenome-wide association study in blood.

RESEARCH DESIGN AND METHODS Five cohort studies’ participants (n = 8,270) were classified as statin users when they were on statin therapy at the time of DNAm assessment with Illumina 450K or EPIC array or noncurrent users otherwise. Associations of DNAm with various outcomes like incident type 2 diabetes, plasma glucose, insulin, and insulin resistance (HOMA of insulin resistance [HOMA-IR]) as well as with gene expression were investigated.

RESULTS Discovery (n = 6,820) and replication (n = 1,450) phases associated five DNAm sites with statin use: cg17901584 (1.12 × 10−25 [DHCR24]), cg10177197 (3.94 × 10−08 [DHCR24]), cg06500161 (2.67 × 10−23 [ABCG1]), cg27243685 (6.01 × 10−09 [ABCG1]), and cg05119988 (7.26 × 10−12 [SC4MOL]). Two sites were associated with at least one glycemic trait or type 2 diabetes. Higher cg06500161 methylation was associated with higher fasting glucose, insulin, HOMA-IR, and type 2 diabetes (odds ratio 1.34 [95% CI 1.22, 1.47]). Mediation analyses suggested that ABCG1 methylation partially mediates the effect of statins on high insulin and HOMA-IR. Gene expression analyses showed that statin exposure and ABCG1 methylation were associated with ABCG1 downregulation, suggesting epigenetic regulation of ABCG1 expression. Further, outcomes insulin and HOMA-IR were significantly associated with ABCG1 expression.

CONCLUSIONS This study sheds light on potential mechanisms linking statins with type 2 diabetes risk, providing evidence on DNAm partially mediating statins’ effects on insulin traits. Further efforts shall disentangle the molecular mechanisms through which statins may induce DNAm changes, potentially leading to ABCG1 epigenetic regulation.

Introduction

Statins effectively reduce the risk of cardiovascular disease (1). However, clinical trials and observational studies show that statins lead to insulin resistance and type 2 diabetes (2,3). The underlying mechanisms remain unclear.

Statins are associated with epigenetic changes, including histone acetylation, miRNA regulation (4), and DNA methylation (DNAm), particularly at genes related to lipid and insulin metabolism (5). DNAm is linked to type 2 diabetes pathophysiology (6); thus, it may be a potential mechanism contributing to the increased risk of type 2 diabetes observed in statin therapy. Nevertheless, this hypothesis has not been investigated (4).

We conducted an epigenome-wide association study (EWAS) in blood investigating the association between statin use and changes in DNAm at sites in the genome called CpGs. Further, we sought to replicate the findings and study the associations of the statin-related CpGs with gene expression, glycemic traits, and type 2 diabetes. Finally, we explored the potential role of the statin-related CpGs mediating the association of statins with glycemic traits and type 2 diabetes.

Research Design and Methods

Study Design and Population

Overall, 8,270 participants from five prospective cohort studies were included. In EWAS analyses, the discovery panel included 6,820 Caucasian individuals from Epidemiologische Studie zu Chancen der Verhütung, Früherkennung und optimierten Therapie chronischer Erkrankungen in der älteren Bevölkerung (ESTHER) (7), a population-based study recruited in Saarland, Germany, with participants aged 50–75 years; Kooperative Gesundheitsforschung in der Region Augsburg-F4 (KORA-F4) (8), a follow-up study of the KORA-S4 cohort, established in 1996 in Augsburg, Germany; Study of Health in Pomerania (SHIP)-Trend (9), a population-based study of participants from northeast Germany; and a population-based cohort of South Asians residing in London, U.K., aged 35–75 years: the London Life Sciences Prospective Population (LOLIPOP) study (10). Our study included 998 individuals from ESTHER, 1,726 from KORA-F4, 247 from SHIP-Trend, and 3,848 from LOLIPOP. Hence, the discovery panel was composed of 56% South Asians and 44% Caucasians.

The replication panel was selected based on ethnic homogeneity and availability of data for post-EWAS analyses. It included 1,450 Caucasians from the Rotterdam Study (RS) (subcohorts RS-III-1 [first visit of the third RS cohort] and RS-BIOS [comprising the third visit of the second cohort (RS-II-3) and the second visit of the third RS cohort (RS-III-2)]) (11), a prospective population-based cohort study of participants living in Ommoord, Rotterdam, the Netherlands, aged 45 years and older. Only RS had information on former statin users (n = 74). Additional association analyses between DNAm and glycemic traits (n = 1,165) and gene expression (n = 731) included RS data, whereas association analyses with incident type 2 diabetes included RS, KORA, and LOLIPOP (n = 4,422). See Fig. 1 for study design.

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Study design. An EWAS was conducted to identify DNAm changes at CpG sites related to current statin use compared with noncurrent statin use. Four cohorts were used as a discovery panel, and findings were replicated in an independent cohort. Association analyses of the identified statin-related CpGs with glycemic traits, incident type 2 diabetes (T2D), and gene expression were performed. The large white arrows show the proceeding steps in the study design. The dotted black arrows represent a potential indirect effect of statins through changes in cg06500161 methylation. The black thin arrows indicate the direct effect of statins on glycemic traits (insulin and HOMA-IR).

Statin Use

Information on medication was obtained from pharmacy records in RS; self-declaration in KORA-F4, LOLIPOP, and SHIP-Trend; and general practitioner records in ESTHER.

Statin use status was categorized in current users and noncurrent users. Users were defined as current users were if the statins prescription occurred at the time of blood drawing for DNAm assessment or as noncurrent users otherwise. RS additionally identified former statin users as participants who had previously used statins but were no longer users on the blood draw date.

DNAm Data

DNA was extracted from whole peripheral blood (stored in EDTA tubes) by standardized salting out methods. Genome-wide DNAm was measured in bisulfite-converted genomic DNA using the Illumina Infinium HumanMethylation450 BeadChip or Illumina Infinium MethylationEPIC BeadChip according to the manufacturer’s protocols. Samples showing inadequate hybridization, incomplete bisulfite treatment, and gender swaps were excluded. Details on cohort-specific methods are provided in Supplementary Table 1. CpG methylation proportion was reported as a normalized β-value ranging from 0 to 1, where 1 represents 100% methylation. For avoidance of sex bias, CpGs annotated to genes in sex chromosomes were excluded. After normalization and quality control, overlapping CpGs across studies were meta-analyzed. CpGs missing in one cohort were also included, while CpGs missing in two or more cohorts were excluded, leaving 332,357 CpGs for final analysis.

mRNA Expression Data

mRNA data included a subset of the replication panel (n = 731 [RS-III-1]). Total RNA was isolated (PAXgene Blood RNA kits; QIAGEN) from whole blood (PAXgene Tubes; Becton Dickinson). Samples were analyzed using the LabChip GX (Caliper Life Sciences) according to the manufacturer’s instructions. Samples with an RNA quality score >7 were amplified, labeled (TotalPrep RNA; Ambion), and hybridized to the Illumina HumanHT-12 v4 Expression BeadChips according to the manufacturer’s protocol. Expression data were quantile normalized to the median distribution and log (2) transformed. Probe and sample means were centered to zero. Genes were considered significantly expressed when detection P values calculated by Genome Studio were <0.05 in >10% of all discovery samples. Thus, 21,238 probes were available to study. Quality control used eQTL (expression quantitative trait loci) mapping pipeline (https://github.com/molgenis/systemsgenetics/tree/master/eqtl-mapping-pipeline) (12). Only probes that uniquely mapped to the human genome build 37 and represented gene mRNA expression were analyzed (13). Samples were processed at the Genetic Laboratory of Internal Medicine, Erasmus Medical Center. Expression data are available for 881 samples at the Gene Expression Omnibus (GEO) public repository (accession no. GSE338828).

Covariates, Glycemic Traits, and Type 2 Diabetes

Covariates were selected based on previous literature on DNAm studies, including age, sex, BMI, smoking status, history of coronary heart disease (CHD), antihypertensive medication, systolic blood pressure (SBP), total cholesterol, HDL, triglycerides, and LDL. BMI was calculated by dividing weight in kilograms by the square of height in meters. Smoking history was self-reported and classified as current, former, and never smokers. SBP was the mean of two consecutive measurements at the right brachial artery with a random-zero sphygmomanometer with the participant in a sitting position. Antihypertensive medication use was defined as the use of diuretics, β-blockers, ACE inhibitors, and calcium channel blockers. CHD was defined as the occurrence of myocardial infarction, revascularization, coronary artery bypass graft surgery, or percutaneous coronary intervention. Triglycerides, HDL, and total cholesterol concentrations were measured in blood using an automated enzymatic method.

Glycemic traits examined were plasma fasting glucose, insulin, HOMA of insulin resistance (HOMA-IR), and type 2 diabetes incidence. Type 2 diabetes cases were defined as the presence of either serum glucose level ≥7.0 mmol/L or glucose-lowering medication use in RS; self-report in KORA; physician diagnosis or use of glucose-lowering medication in ESTHER; diabetes medication use, HbA1c ≥6.5%, serum glucose level (nonfasting) ≥11.1 mmol/L, or self-report in SHIP-Trend; and diabetes history or fasting glucose ≥7.0 mmol/L or HbA1c >6.5% in LOLIPOP. The interassay coefficients of variations are insulin <8%, glucose <1.4%, and lipids <2.1%. LDL was calculated using the Friedewald formula (total cholesterol − HDL − triglycerides/5).

Statistical Analyses

EWAS

The association between CpG methylation and current versus noncurrent statin users was investigated in a two-step EWAS approach: discovery and replication. Normalized DNAm β-values at each CpG were the dependent variable, and current statin use (yes = 1, no = 0) was the predictor. Per individual CpG, participants with methylation levels higher than three times the interquartiles range were excluded. The regression analysis used the lme4 package in R, version 3.4.2 (https://cran.r-project.org/web/packages/lme4/index.html), to perform linear mixed effects models. The random part of the model included adjustment for array number and position on array as random terms to account for batch effects in the DNAm measurement. The fixed part of the model included leukocyte proportions (to account for cell mixture), sex, age, smoking status, BMI, SBP, antihypertensive medication, and prevalent CHD and type 2 diabetes. Leukocyte proportions (B cells, CD4 T cells, CD8 T cells, granulocytes, monocytes, and NK cells) were estimated as previously described by Houseman (14), using the R minfi package (15). Statistical significance for discovery was determined by a Bonferroni-corrected threshold (1.0 × 10−7). Heterogeneity (I2) of effect estimates was assessed to account for differences between cohorts. Results across the discovery cohorts were combined in a meta-analysis by inverse variance–weighted method in METAL, version 2011-03-25.

The identified CpGs were replicated in an independent panel using identical models. Bonferroni correction for replication was 0.05 divided by the number of discovery findings. Results of discovery and replication panels were combined in a meta-analysis as aforementioned.

Pearson test was used to assess the correlation among the replicated CpGs in the replication panel (RS).

Associations Between Statin-Related CpGs and Glycemic Traits and Incident Type 2 Diabetes

Standardized methylation values at the replicated CpGs as predictors of fasting glucose, insulin, and HOMA-IR (RS [n = 1,165]) were tested in linear mixed effects models. For exploration of reverse causation, prevalent type 2 diabetes case subjects and former statin users were excluded. Logistic regression was used to assess the association with incident type 2 diabetes (RS, KORA, and LOLIPOP [n = 4,422, including 1,173 cases]). Results across cohorts were meta-analyzed as aforementioned. The adjusting covariates included leukocyte proportions, batch effects, sex, age, smoking status, SBP, antihypertensive medication, CHD, and statin use. Bonferroni correction was applied (<0.0025). Variables with right-skewed distributions were transformed to the natural logarithmic scale (glucose, insulin, HOMA-IR, and HDL).

Causal Mediation Analysis

To investigate whether DNAm changes are part of the pathway though which statins exert their diabetogenic effect, we explored the significant associations in the previous stage in a nonparametric causal mediation analysis using the R mediation package (16) as previously described (17). The mediation analysis dissects the total effect of the exposure (statin treatment) on the outcome (glycemic traits and incident type 2 diabetes) into the 1) indirect effect acting via the mediator of interest (DNAm) and 2) direct effect acting directly or via a mediator other than the one under study. The mediation package compares the effect estimates of each association between statin use and traits in subjects with different methylation levels. Nonparametric bootstrap with 20,000 resamples was performed. Data from the RS were used (n = 1,180), with exclusion of prevalent type 2 diabetes case subjects and former statin users. The normalized CpG methylation β-value was modeled as mediator and the CpG-associated trait as outcome. Supplementary Fig. 1 shows a schematic representation. Sex, age, and cohort were considered confounders. For investigation of potential unmeasured confounding, the sequential ignorability assumption (no confounder assumption) was tested in a sensitivity analysis using the correlated residuals method (18), which estimates the correlation that quantifies how strong the confounder would have to be to change the conclusion.

Gene Expression-Association Analysis

For investigation of possible biological pathways, relations between the identified CpGs and transcriptome-wide gene expression were explored (RS) (n = 731). Statin-related CpGs were regressed against methylation batch variables as random terms; a similar approach was applied for 21,238 transcription probes and the expression batch variables. Next, linear regressions between the residuals from the CpGs regression (predictors) and the expression probe residuals (dependent variable) were performed. The Bonferroni correction was 0.05/number of probes × number of CpGs tested.

Additionally, the significant expression probes in the aforementioned analysis were tested in association with statin use exposure (n = 631) and as predictors of the significant glycemic traits in the mediation analyses (n = 616). Type 2 diabetes case subjects and former statin users were excluded.

Dose Effects

In 303 current statin users from SHIP-Trend and RS with available data on statin dose, we examined the association between the replicated CpGs and the defined daily dose (DDD) (WHO ATC/DDD classification [https://www.whocc.no/atc_ddd_index/]), using the EWAS model. Next, the association between gene expression and statin dose was studied in the RS.

Confounding and Ethnicity

To investigate whether the associations of statin use with the replicated CpGs were independent of blood lipids, we additionally adjusted the EWAS for serum total cholesterol or for individual blood lipids instead (HDL, LDL, and triglycerides) in both discovery and replication.

In an effort to rule out confounding by indication, we restricted the EWAS to subjects with LDL ≥70 mg/dL or ≥100 mg/dL—cutoffs used to advise statin prescription (19,20). Further, we reran the EWAS excluding prevalent type 2 diabetes and prediabetes case subjects.

To account for ethnicity heterogeneity in the discovery panel (South Asians and Caucasians), we reran the EWAS, taking Caucasians (n = 4,349) for discovery and South Asians (n = 3,849) for replication.

Post Hoc Power Analysis

To calculate the range of power for replication of the findings, we used the GPower 3.1 tool (21) with six predictors representing all discovered CpGs except for cg03467813, which was not included due to its high heterogeneity (I2 = 99%). CpGs’ smallest β-estimate (0.0039 [cg10177197]) and probability error 0.008 or the largest estimate (0.0219 [cg17901584]) was used.

Results

Association Between Statin Use and DNA Methylation

The mean (SD) age in the discovery panel ranged from 51.5 (13.8) years (SHIP-Trend) to 62.1 (6.5) years (KORA-F4). Of the total discovery panel, 2,969 (43.5%) participants were women. Supplementary Table 2A and B shows other baseline characteristics.

In the discovery panel, seven CpGs passed the Bonferroni threshold for significance (P < 1.0 × 10−7), being differentially methylated in current statins users compared with noncurrent users, of which current users had lower methylation at four CpGs, annotated to genes DHCR24 (cg17901584, located at TSS1500), FAM50B (cg03467813, located at the gene body region), SC4MOL (cg05119988, located at 5′ untranslated region [5′UTR]) and AHRR (cg05575921, located at gene body). Three CpGs annotated to ABCG1 (cg06500161 and cg27243685, located at gene body) and DHCR24 (cg10177197, located at gene body) showed higher methylation among current users (Table 1 and Supplementary Fig. 1). All associations were independently replicated, except for two CpGs annotated to AHRR and FAM50B, which failed to reach statistical significance after Bonferroni correction for replication (P < 7.14 × 10−3). The post hoc power analysis showed a power ranging from 39% to 99% for our replication study. In the replication panel, correlations among CpGs ranged from r = −0.011 between cg06500161 (ABCG1) and cg17901584 (DHCR24) to r = 0.477 for cg06500161 and cg27243685 (both annotated to ABCG1) (Supplementary Table 3).

View this table:
  • View inline
  • View popup
Table 1

DNAm sites associated with current statin use (compared with noncurrent use)

The meta-analysis across the discovery and replication panels revealed one new CpGs at DHCR24 (cg17475467). The heterogeneity in the associations between statin use and the replicated CpGs across cohorts varied from I2 = 0 for cg10177197 (DHCR24) to I2 = 74.8 for cg17901584 (DHCR24) (Table 1). Unreplicated CpGs were not further investigated.

The sensitivity analysis using RS data with never statin users as the reference group (former users excluded) revealed similar results, replicating the five CpGs (Supplementary Table 4).

Associations of the Statin-Related CpGs With Glycemic Traits and Type 2 Diabetes

After a Bonferroni correction (2.5 × 10−3) and comprehensive assessment of potential confounders, increase in 1 SD of methylation at cg06500161 (ABCG1) was associated with incident type 2 diabetes (odds ratio 1.34 [95% CI 1.22, 1.47]), fasting glucose (1.03 × 10−3), insulin (4.63 × 10−8), and HOMA-IR (1.05 × 10−8). Increase in 1 SD of methylation at cg27243685 (ABCG1) was associated with incident type 2 diabetes (1.24 [1.13, 1.38]) and approaching borderline statistical significance with augmented insulin (4.62 × 10−3) but not statistically significantly associated with glucose or HOMA-IR. Increase in 1 SD of methylation at cg17901584 was approaching borderline statistical significance with incident type 2 diabetes (6.2 × 10−3) but was not significantly associated with glucose, insulin, or HOMA-IR. No significant findings were observed for cg05119988 (SC4MOL) or cg10177197 (DHCR24) (Table 2).

View this table:
  • View inline
  • View popup
Table 2

Association of DNAm of the replicated CpGs associated with statin use with glycemic traits and type 2 diabetes

Causal Mediation

Significant results were obtained for the models with outcome fasting insulin and HOMA-IR. For both cases, the analyses showed that 1) current statin use has a significant overall effect of β = 0.275 (0.190, 0.360) on insulin and β = 0.291 (0.199, 0.380) on HOMA-IR and 2) part of that effect goes directly or via mediator(s) other than cg06500161, called average direct effect (ADE). We found an ADE of β = 0.233 (0.149, 0.320) on insulin and β = 0.242 (0.151, 0.330) on HOMA-IR. 3) Another part of such effect may operate via an indirect path (indirect effect), possibly through cg06500161 methylation, with an average causal mediator effect of β = 0.043 (0.024, 0.060) on insulin and β = 0.049 (0.028, 0.070) on HOMA-IR; 4) consequently, of the total effect of statins on insulin and HOMA-IR, 15.5% (0.086, 0.260) and 16.8% (0.095, 0.270) are suggested to act via cg06500161 methylation, respectively. According to these results, cg06500161 methylation does not explain 100% of the effect of statins on insulin/HOMA-IR; thus, we may call cg06500161 methylation a partial mediator instead of a total mediator.

No significant findings were observed for the mediation model with outcome glucose (Table 3). We also investigated the potential mediation of cg27243685 on incident type 2 diabetes, with no significant findings. The sensitivity analysis showed that residual correlation due to unmeasured confounding of r ≥ 0.2 would be needed to violate the sequential ignorability assumption.

View this table:
  • View inline
  • View popup
Table 3

Causal mediation analysis on the significant associations between the statin-related CpGs (cg06500161 and cg27243685) as mediators and glycemic traits and incident type 2 diabetes as outcomes

Further mediation models were performed to account for potential confounders. In addition to age, sex, and cohort, we further adjusted for 1) BMI, 2) BMI and HDL, and 3) smoking status, SBP, antihypertensive medication, and presence of CHD. Results did not materially change with these further adjustments, although the effect sizes attenuated (Supplementary Tables 5–7).

Gene Expression

For exploration of potential epigenetic regulations of gene expression by DNAm, the identified CpGs (predictors) were tested in association with transcriptome-wide gene expression in the RS. cg06500161 was inversely associated with two expression probes at ABCG1 (ILMN_1794782 and ILMN_2329927) after Bonferroni correction (4.71 × 10−7). An inverse association was observed for cg27243685 and one probe (ILMN_1794782). These results indicate a lower expression of ABCG1 with increasing cg06500161 and cg27243685 methylation. Furthermore, cg17901584 methylation was directly associated with ABCG1 expression (ILMN_1794782) and suggestively associated with another probe at the same gene (ILMN_2329927). There were no significant findings for the other two CpGs after Bonferroni correction (Supplementary Table 8). Moreover, exposure to statin treatment was associated with lower ABCG1 expression (Supplementary Table 9). When ABCG1 expression was studied as a predictor of insulin traits, we found an inverse association between ILMN_1794782 and fasting insulin and HOMA-IR (Supplementary Table 10).

Dose Effect

Statin dose exposure was significantly associated with the five identified CpGs after Bonferroni correction (P < 0.01) (Supplementary Table 11). Exposure to increasing statin dose was associated nominally with lower levels of ABCG1 expression at ILMN_1794782 and suggestively with ILMN_2329927 (Supplementary Table 12).

Sensitivity Analyses

Residual confounding and confounding by indication were as follows: 1) additional adjustment of the EWAS for serum total cholesterol or blood lipids HDL, LDL, and triglycerides instead (Supplementary Table 13); 2) exclusion of prevalent type 2 diabetes and pre–type 2 diabetes case subjects (Supplementary Table 14); or 3) restriction of the EWAS to participants with LDL ≥70 mg/dL or ≥100 mg/dL (Supplementary Table 15) did not change the associations of the five CpGs. Further, exclusion of type 2 diabetes case subjects from the other analyses did not affect the results (data not shown).

With regard to transethnic replication, the new EWAS taking Caucasians only for discovery and South Asians for replication showed that all five CpGs passed the Bonferroni correction threshold (Supplementary Table 16).

Conclusions

The current study sheds light on potential mechanisms linking statin use and type 2 risk diabetes by, firstly, identifying and replicating associations between statin therapy and methylation at five CpGs (cg06500161, cg27243685, cg17901584, cg10177197, and cg05119988) independent of blood lipids and, secondly, providing evidence on the partial mediation via ABCG1 methylation in the effect of statins on increased levels of insulin and insulin resistance under the sequential ignorability assumption.

Our results on the association between statin use and DNAm agree with a report from the Framingham Study, which found epigenome-wide significance at only two sites (cg17901584 and cg06500161) (5). The Framingham Study examined this association in a smaller sample (n = 1,545 vs. 6,820 for our discovery panel), and contrary to our investigation, they did not proceed with replication or examine associations among gene expression, metabolic markers, and type 2 diabetes. Moreover, our work reports novel associations and adds to current knowledge that statin dosage might be implicated in the degree of methylation and ABCG1 expression. The latter finding goes in line with an experimental study where macrophages treated with various types of statins showed lower ABCG1 expression as dose increased (22).

Sites cg06500161 and cg27243685 are annotated to ABCG1, cg17901584 and cg10177197 to DHCR24, and site cg05119988 is annotated to SC4MOL. ABCG1 (ATP-binding cassette member-1 subfamily-G) encodes a protein that mediates the transport of different molecules, such as cholesterol efflux to the HDL, oxysterols, and phospholipid transport in macrophages (23,24). It is also involved in insulin secretion and sensitivity (25). ABCG1 expression was reduced in statin users compared with nonusers (26) and in patients with type 2 diabetes (27,28). Genes DHCR24 (24-dehydrocholesterol reductase) and SC4MOL (sterol-C4-methyl oxidase-like) code enzymes catalyzing different steps during cholesterol biosynthesis. DHCR24 mutations are related to desmosterolosis (29,30) and Alzheimer disease (31). Deficiency of the protein coded by SC4MOL produces congenital cataracts, microcephaly, growth delay, skin conditions, and immune dysfunction (32). An observational study found that DHCR24 and SC4MOL were upregulated among statin users (26), while SC4MOL upregulation increased type 2 diabetes risk (27,28).

Methylation at four of the identified CpGs (cg06500161, cg27243685, cg17901584, and cg05119988) was previously associated with glycemic traits, type 2 diabetes, and blood lipids, while cg06500161 (ABCG1) was associated with increased glucose, insulin, HbA1c, HDL, and triglyceride levels (33–35). Our causal mediation analyses provided evidence of cg06500161 methylation partially mediating the association between statin use and higher fasting insulin and HOMA-IR. However, given the cross-sectional nature of these analyses, our results must be interpreted with caution. Since the mediator was measured at the same time as the outcome, reverse causation cannot be completely ruled out. However, in an effort to overcome this issue, we excluded type 2 diabetes case subjects. The mediation effect observed may be the result of more complex interactions of combined effects across the epigenome. Moreover, the lack of significant findings in mediation by the other CpG at ABCG1 (cg27243685) could reflect a power issue. Hence, future studies shall explore how the CpGs within the same gene interact and mediate the effects.

Based on our findings and the available evidence, we hypothesize that cg06500161 hypermethylation may be a consequence of statin use, possibly inducing a decrease in ABCG1 transcription in blood. Further research is needed to investigate to what extent this downregulation could in turn compromise downstream signals, resulting in impaired insulin metabolism. In this line, impaired insulin sensitivity and secretion as a consequence of statin treatment were recently observed in a longitudinal study (36). Furthermore, a functional study suggested an epigenetic regulation of ABCG1 mediated by methylation-dependent transcription factor binding (37). Our hypothesized model is displayed in Fig. 2. Nevertheless, further investigations should test this hypothesis and assess the effect of statins in target tissues and to what extent these findings are generalizable to populations with ethnic backgrounds other than Caucasians and South Asians.

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

Scheme of the hypothesized mechanism linking statin therapy and risk of type 2 diabetes. The solid black arrow represents the effect of statins on plasmatic insulin and HOMA-IR levels that operates directly or through a pathway different from the mediator analyzed in the current study (methylation at cg06500161). The dotted black arrows represent the suggested alternative pathway, where an indirect effect of statins on insulin and HOMA-IR is mediated by cg06500161 methylation. The black thin arrows indicate decreased or increased level of the parameter; in this figure they are representing decrease in gene expression of ABCG1 and increase in serum levels of fasting insulin and HOMA-IR.

We highlight the large sample size with DNAm data in our study, the use of a replication panel, and the comprehensive assessment of confounding factors. Addition of complementary data like gene expression, transethnic replication, and use of a causal inference method as mediation analysis also add value to our research.

One limitation is the lack of data on DNAm and gene expression from specific tissues of interest for type 2 diabetes and drug metabolism (pancreas, liver, adipose tissue), as both DNAm and gene expression may be tissue specific (38). Instead, we used whole peripheral blood, allowing for cell mixture. Consequently, we included adjustments for estimated blood leukocytes proportions.

Most cohorts had limited information on former statin users; thus, it is likely that they were included in the reference group. However, the replication panel showed that the same results are valid for the group of current statin users only, although we had limited numbers to investigate such associations. Furthermore, it was not possible to explore potential effects of cessation time or treatment duration on DNAm. A Mendelian randomization approach based on genetic instrumental variables to explore a potential causal effect of statins on DNAm was not possible due to the lack of adequate genetic variants that may mimic the randomization of the hypothesized adverse effect of statins on type 2 diabetes risk. The use of variants at the gene encoding 3-hydroxy-3-methylglutaryl CoA reductase (39), statins’ target protein, would mimic the main effect of statins, rather than an alternative pathway leading to the hypothesized effect. Moreover, these variants were associated with LDL (39), reflecting cholesterol metabolism. This violates the Mendelian randomization assumption that the instrumental variable is associated with the outcome only through the exposure, since high cholesterol is a risk factor for type 2 diabetes.

All cohorts used Illumina 450K array to measure DNAm, except one, which used Illumina EPIC array (n = 247). Despite differences between the two, EPIC array covers ∼90% of the same sites on the 450K chip; both use the same technology and showed overall per-sample correlations of r > 0.99 (40), thus encouraging the application of combined data. Moreover, cohorts used different methods to diagnose type 2 diabetes, which might introduce error measurement.

The heterogeneity observed in the meta-analysis for some of the CpGs might be a consequence of differences in statin use prevalence (ranged from 9.2% to 18.1%) and types and dose of statin used across cohorts. Additionally, differences in ethnicity could have influenced heterogeneity (the association of statin use and cg06500161 in the overall meta-analysis showed I2 = 52.30, while I2 = 5.20 across Caucasians only), as could different methods applied to collect data on medication and disease (pharmacy and general practitioner records or self-report).

Finally, some analyses were performed using data from the replication panel only, which may introduce a type II error.

Conclusion

We identified and replicated five DNAm sites associated with statin therapy, of which cg06500161 may be partially mediating the effect of statins on increasing insulin levels and HOMA-IR. This could be one potential mechanism linking statin therapy and type 2 diabetes onset observed in clinical trials and observational studies. Nevertheless, other biological pathways should not be discarded. From our study, it is not clear what might be the mechanism by which statins may induce DNAm changes; this shall be a challenge for future research.

Meanwhile, close monitoring of risk factors for type 2 diabetes in these patients is imperative to prevent further complications.

Article Information

Funding. C.O.-R. reports receiving grant support from Agencia Nacional de Investigación y Desarrollo (ANID)/Scholarship Program/Doctorado Becas Chile/2016 - 72170524. P.P.M. and T.L. reported receiving financial support from the Academy of Finland [grants 286284, 134309 (Eye), 126925, 121584, 124282, 129378 (Salve), 117787 (Gendi), and 41071 (Skidi)], the Social Insurance Institution of Finland, Competitive State Research Financing of the Expert Responsibility area of Kuopio, Tampere and Turku University Hospitals (grant X51001), Juho Vainio Foundation, Paavo Nurmi Foundation, Finnish Foundation for Cardiovascular Research, Finnish Cultural Foundation, Tampere Tuberculosis Foundation, Emil Aaltonen Foundation, Yrjö Jahnsson Foundation, Signe and Ane Gyllenberg Foundation, Diabetes Research Foundation of Finnish Diabetes Association, and EU Horizon 2020 (grant 755320 for TAXINOMISIS Project). H.J.G., S.B.F., and A.T. report receiving financial support from the Federal Ministry of Education and Research (grants 01ZZ9603, 01ZZ0103, and 01ZZ0403), the Ministry of Cultural Affairs as well as the Social Ministry of the Federal State of Mecklenburg-West Pomerania, and the network “Greifswald Approach to Individualized Medicine (GANI_MED)” funded by the Federal Ministry of Education and Research (grant 03IS2061A) and the DZHK (grant 81X34000104). J.C. is supported by the Singapore Ministry of Health’s National Medical Research Council under its Singapore Translational Research Investigator (STaR) Award (NMRC/STaR/0028/2017). KORA was initiated and financed by Helmholtz Zentrum München (German Research Center for Environmental Health), which is funded by the German Federal Ministry of Education and Research (BMBF) and by the State of Bavaria. Furthermore, KORA research has been supported within the Munich Center of Health Sciences (MC-Health), Ludwig-Maximilians-Universität, as part of LMUinnovativ. This work was supported by a grant (WA 4081/1-1) from the German Research Foundation.

Duality of Interest. H.J.G. has received travel grants and speakers honoraria from Fresenius Medical Care, Neuraxpharm, Servier, and Janssen Cilag, as well as research funding from Fresenius Medical Care. O.H.F. reports receiving grants or research support from Metagenics, Inc., and Nestlé Nutrition (Nestec Ltd.). No other potential conflicts of interest relevant to this article were reported.

Author Contributions. T.M. and C.O.-R. contributed to the study concept and design. C.O.-R., R.W., B.L., P.P.M., X.G., M.L., D.J.-Q., J.B.J.v.M., W.Z., J.S.K., H.J.G., S.B.F., B.S., Y.Z., C.G., M.M.-N., M.H., A.P., and A.T. contributed to data generation and analyses. C.O.-R., T.M., R.W., B.L., P.P.M., X.G., D.J.-Q., M.L., E.P.-F., J.N., O.H.F., M.G., O.L.R.-O., A.T., H.B., M.W., M.A.I., J.B.J.v.M., T.V., J.C., and B.H.S. contributed to drafting of the manuscript. O.H.F., M.G., O.L.R.-O., T.L., A.T., H.B., M.W., M.A.I., J.B.J.v.M., T.V., J.C., B.H.S., and T.M. contributed to critical revision of the manuscript for important intellectual content. All authors approved the final version of the manuscript. C.O.-R. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Prior Presentation. Parts of this study were presented in abstract form at the 78th Scientific Sessions of the American Diabetes Association, Orlando, FL, 22–26 June 2018, and at EuroPrevent 2019, Lisbon, Portugal, 11–13 April 2019.

Footnotes

  • This article contains Supplementary Data online at https://care.diabetesjournals.org/lookup/suppl/doi:10.2337/dc19-1828/-/DC1.

  • Received September 13, 2019.
  • Accepted January 14, 2020.
  • © 2020 by the American Diabetes Association.
https://www.diabetesjournals.org/content/license

Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at https://www.diabetesjournals.org/content/license.

References

  1. ↵
    1. Cholesterol Treatment Trialists’ (CTT) Collaboration;
    2. Baigent C,
    3. Blackwell L,
    4. Emberson J, et al
    . Efficacy and safety of more intensive lowering of LDL cholesterol: a meta-analysis of data from 170,000 participants in 26 randomised trials. Lancet 2010;376:1670–1681
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Casula M,
    2. Mozzanica F,
    3. Scotti L, et al
    . Statin use and risk of new-onset diabetes: a meta-analysis of observational studies. Nutr Metab Cardiovasc Dis 2017;27:396–406
    OpenUrl
  3. ↵
    1. Thakker D,
    2. Nair S,
    3. Pagada A,
    4. Jamdade V,
    5. Malik A
    . Statin use and the risk of developing diabetes: a network meta-analysis. Pharmacoepidemiol Drug Saf 2016;25:1131–1149
    OpenUrl
  4. ↵
    1. Allen SC,
    2. Mamotte CDS
    . Pleiotropic and adverse effects of statins-do epigenetics play a role? J Pharmacol Exp Ther 2017;362:319–326
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Dogan MV,
    2. Grumbach IM,
    3. Michaelson JJ,
    4. Philibert RA
    . Integrated genetic and epigenetic prediction of coronary heart disease in the Framingham Heart Study. PLoS One 2018;13:e0190549
    OpenUrl
  6. ↵
    1. Muka T,
    2. Nano J,
    3. Voortman T, et al
    . The role of global and regional DNA methylation and histone modifications in glycemic traits and type 2 diabetes: a systematic review. Nutr Metab Cardiovasc Dis 2016;26:553–566
    OpenUrlPubMed
  7. ↵
    1. Raum E,
    2. Rothenbacher D,
    3. Löw M,
    4. Stegmaier C,
    5. Ziegler H,
    6. Brenner H
    . Changes of cardiovascular risk factors and their implications in subsequent birth cohorts of older adults in Germany: a life course approach. Eur J Cardiovasc Prev Rehabil 2007;14:809–814
    OpenUrlCrossRefPubMedWeb of Science
  8. ↵
    1. Rathmann W,
    2. Haastert B,
    3. Icks A, et al
    . High prevalence of undiagnosed diabetes mellitus in Southern Germany: target populations for efficient screening. The KORA survey 2000. Diabetologia 2003;46:182–189
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    1. Völzke H,
    2. Alte D,
    3. Schmidt CO, et al
    . Cohort profile: the study of health in Pomerania. Int J Epidemiol 2011;40:294–307
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    1. Chahal NS,
    2. Lim TK,
    3. Jain P,
    4. Chambers JC,
    5. Kooner JS,
    6. Senior R
    . Does subclinical atherosclerosis burden identify the increased risk of cardiovascular disease mortality among United Kingdom Indian Asians? A population study. Am Heart J 2011;162:460–466
    OpenUrlPubMed
  11. ↵
    1. Ikram MA,
    2. Brusselle GGO,
    3. Murad SD, et al
    . The Rotterdam Study: 2018 update on objectives, design and main results. Eur J Epidemiol 2017;32:807–850
    OpenUrlCrossRef
  12. ↵
    1. Westra HJ,
    2. Peters MJ,
    3. Esko T, et al
    . Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat Genet 2013;45:1238–1243
    OpenUrlCrossRefPubMed
  13. ↵
    1. Schurmann C,
    2. Heim K,
    3. Schillert A, et al
    . Analyzing illumina gene expression microarray data from different tissues: methodological aspects of data analysis in the metaxpress consortium. PLoS One 2012;7:e50938
    OpenUrlCrossRefPubMed
  14. ↵
    1. Houseman EA,
    2. Accomando WP,
    3. Koestler DC, et al
    . DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 2012;13:86
    OpenUrlCrossRefPubMed
  15. ↵
    1. Aryee MJ,
    2. Jaffe AE,
    3. Corrada-Bravo H, et al
    . Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014;30:1363–1369
    OpenUrlCrossRefPubMedWeb of Science
  16. ↵
    1. Tingley D,
    2. Yamamoto T,
    3. Hirose K,
    4. Keele L,
    5. Imai K
    . mediation: R Package for Causal Mediation Analysis. J Stat Softw 2014;59
  17. ↵
    1. Steenaard RV,
    2. Ligthart S,
    3. Stolk L, et al
    . Tobacco smoking is associated with methylation of genes related to coronary artery disease. Clin Epigenetics 2015;7:54
    OpenUrlCrossRef
  18. ↵
    1. Imai K,
    2. Keele L,
    3. Yamamoto T
    . Identification, inference and sensitivity analysis for causal mediation effects. Stat Sci 2010;25:51–71
    OpenUrlCrossRefWeb of Science
  19. ↵
    1. Piepoli MF,
    2. Hoes AW,
    3. Agewall S, et al
    . 2016 European guidelines on cardiovascular disease prevention in clinical practice. Kardiol Pol 2016;74:821–936[in Polish]
    OpenUrl
  20. ↵
    1. Perk J,
    2. De Backer G,
    3. Gohlke H, et al.; European Association for Cardiovascular Prevention & Rehabilitation (EACPR)
    . European guidelines on cardiovascular disease prevention in clinical practice (version 2012): the fifth joint task force of the European society of cardiology and other societies on cardiovascular disease prevention in clinical practice (constituted by representatives of nine societies and by invited experts). Int J Behav Med 2012;19:403–488
    OpenUrlCrossRefPubMed
  21. ↵
    1. Faul F,
    2. Erdfelder E,
    3. Buchner A,
    4. Lang AG
    . Statistical power analyses using G*Power 3.1: tests for correlation and regression analyses. Behav Res Methods 2009;41:1149–1160
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. Wong J,
    2. Quinn CM,
    3. Gelissen IC,
    4. Jessup W,
    5. Brown AJ
    . The effect of statins on ABCA1 and ABCG1 expression in human macrophages is influenced by cellular cholesterol levels and extent of differentiation. Atherosclerosis 2008;196:180–189
    OpenUrl
  23. ↵
    1. Vaughan AM,
    2. Oram JF
    . ABCA1 and ABCG1 or ABCG4 act sequentially to remove cellular cholesterol and generate cholesterol-rich HDL. J Lipid Res 2006;47:2433–2443
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Engel T,
    2. Fobker M,
    3. Buchmann J, et al
    . 3β,5α,6β-cholestanetriol and 25-hydroxycholesterol accumulate in ATP-binding cassette transporter G1 (ABCG1)-deficiency. Atherosclerosis 2014;235:122–129
    OpenUrlCrossRefPubMed
  25. ↵
    1. Olivier M,
    2. Tanck MW,
    3. Out R, et al
    . Human ATP-binding cassette G1 controls macrophage lipoprotein lipase bioavailability and promotes foam cell formation. Arterioscler Thromb Vasc Biol 2012;32:2223–2231
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Obeidat M,
    2. Fishbane N,
    3. Nie Y, et al
    . The effect of statins on blood gene expression in COPD. PLoS One 2015;10:e0140022
    OpenUrlCrossRefPubMed
  27. ↵
    1. Ding J,
    2. Reynolds LM,
    3. Zeller T, et al
    . Alterations of a cellular cholesterol metabolism network are a molecular feature of obesity-related type 2 diabetes and cardiovascular disease. Diabetes 2015;64:3464–3474
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Mauldin JP,
    2. Nagelin MH,
    3. Wojcik AJ, et al
    . Reduced expression of ATP-binding cassette transporter G1 increases cholesterol accumulation in macrophages of patients with type 2 diabetes mellitus. Circulation 2008;117:2785–2792
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Waterham HR,
    2. Koster J,
    3. Romeijn GJ, et al
    . Mutations in the 3beta-hydroxysterol Delta24-reductase gene cause desmosterolosis, an autosomal recessive disorder of cholesterol biosynthesis. Am J Hum Genet 2001;69:685–694
    OpenUrlCrossRefPubMedWeb of Science
  30. ↵
    1. Andersson HC,
    2. Kratz L,
    3. Kelley R
    . Desmosterolosis presenting with multiple congenital anomalies and profound developmental delay. Am J Med Genet 2002;113:315–319
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    1. Lämsä R,
    2. Helisalmi S,
    3. Hiltunen M, et al
    . The association study between DHCR24 polymorphisms and Alzheimer’s disease. Am J Med Genet B Neuropsychiatr Genet 2007;144B:906–910
    OpenUrl
  32. ↵
    1. He M,
    2. Kratz LE,
    3. Michel JJ, et al
    . Mutations in the human SC4MOL gene encoding a methyl sterol oxidase cause psoriasiform dermatitis, microcephaly, and developmental delay. J Clin Invest 2011;121:976–984
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    1. Dekkers KF,
    2. van Iterson M,
    3. Slieker RC, et al.; BIOS Consortium
    . Blood lipids influence DNA methylation in circulating cells. Genome Biol 2016;17:138
    OpenUrlCrossRef
    1. Hidalgo B,
    2. Irvin MR,
    3. Sha J, et al
    . Epigenome-wide association study of fasting measures of glucose, insulin, and HOMA-IR in the Genetics of Lipid Lowering Drugs and Diet Network study. Diabetes 2014;63:801–807
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Dayeh T,
    2. Tuomi T,
    3. Almgren P, et al
    . DNA methylation of loci within ABCG1 and PHOSPHO1 in blood DNA is associated with future type 2 diabetes risk. Epigenetics 2016;11:482–488
    OpenUrl
  35. ↵
    1. Cederberg H,
    2. Stančáková A,
    3. Yaluri N,
    4. Modi S,
    5. Kuusisto J,
    6. Laakso M
    . Increased risk of diabetes with statin treatment is associated with impaired insulin sensitivity and insulin secretion: a 6 year follow-up study of the METSIM cohort. Diabetologia 2015;58:1109–1117
    OpenUrlCrossRefPubMed
  36. ↵
    1. Pfeiffer L,
    2. Wahl S,
    3. Pilling LC, et al
    . DNA methylation of lipid-related genes affects blood lipid levels. Circ Cardiovasc Genet 2015;8:334–342
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Lokk K,
    2. Modhukur V,
    3. Rajashekar B, et al
    . DNA methylome profiling of human tissues identifies global and tissue-specific methylation patterns [published correction appears in Genome Biol 2016;17:224]. Genome Biol 2014;15:r54
    OpenUrlCrossRefPubMed
  38. ↵
    1. Burkhardt R,
    2. Kenny EE,
    3. Lowe JK, et al
    . Common SNPs in HMGCR in micronesians and whites associated with LDL-cholesterol levels affect alternative splicing of exon13. Arterioscler Thromb Vasc Biol 2008;28:2078–2084
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Solomon O,
    2. MacIsaac J,
    3. Quach H, et al
    . Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children. Epigenetics 2018;13:655–664
    OpenUrl
PreviousNext
Back to top
Diabetes Care: 43 (4)

In this Issue

April 2020, 43(4)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by Author
  • Masthead (PDF)
Sign up to receive current issue alerts
View Selected Citations (0)
Print
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about Diabetes Care.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Epigenetic Link Between Statin Therapy and Type 2 Diabetes
(Your Name) has forwarded a page to you from Diabetes Care
(Your Name) thought you would like to see this page from the Diabetes Care web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Epigenetic Link Between Statin Therapy and Type 2 Diabetes
Carolina Ochoa-Rosales, Eliana Portilla-Fernandez, Jana Nano, Rory Wilson, Benjamin Lehne, Pashupati P. Mishra, Xu Gao, Mohsen Ghanbari, Oscar L. Rueda-Ochoa, Diana Juvinao-Quintero, Marie Loh, Weihua Zhang, Jaspal S. Kooner, Hans J. Grabe, Stephan B. Felix, Ben Schöttker, Yan Zhang, Christian Gieger, Martina Müller-Nurasyid, Margit Heier, Annette Peters, Terho Lehtimäki, Alexander Teumer, Hermann Brenner, Melanie Waldenberger, M. Arfan Ikram, Joyce B.J. van Meurs, Oscar H. Franco, Trudy Voortman, John Chambers, Bruno H. Stricker, Taulant Muka
Diabetes Care Apr 2020, 43 (4) 875-884; DOI: 10.2337/dc19-1828

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Add to Selected Citations
Share

Epigenetic Link Between Statin Therapy and Type 2 Diabetes
Carolina Ochoa-Rosales, Eliana Portilla-Fernandez, Jana Nano, Rory Wilson, Benjamin Lehne, Pashupati P. Mishra, Xu Gao, Mohsen Ghanbari, Oscar L. Rueda-Ochoa, Diana Juvinao-Quintero, Marie Loh, Weihua Zhang, Jaspal S. Kooner, Hans J. Grabe, Stephan B. Felix, Ben Schöttker, Yan Zhang, Christian Gieger, Martina Müller-Nurasyid, Margit Heier, Annette Peters, Terho Lehtimäki, Alexander Teumer, Hermann Brenner, Melanie Waldenberger, M. Arfan Ikram, Joyce B.J. van Meurs, Oscar H. Franco, Trudy Voortman, John Chambers, Bruno H. Stricker, Taulant Muka
Diabetes Care Apr 2020, 43 (4) 875-884; DOI: 10.2337/dc19-1828
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Research Design and Methods
    • Results
    • Conclusions
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Impact of Obesity on Measures of Cardiovascular and Kidney Health in Youth With Type 1 Diabetes as Compared With Youth With Type 2 Diabetes
  • Serum Galectin-3 and Subsequent Risk of Coronary Heart Disease in Subjects With Childhood-Onset Type 1 Diabetes: A Cohort Study
  • OGTT Glucose Response Curves, Insulin Sensitivity, and β-Cell Function in RISE: Comparison Between Youth and Adults at Randomization and in Response to Interventions to Preserve β-Cell Function
Show more Cardiovascular and Metabolic Risk

Similar Articles

Navigate

  • Current Issue
  • Standards of Care Guidelines
  • Online Ahead of Print
  • Archives
  • Submit
  • Subscribe
  • Email Alerts
  • RSS Feeds

More Information

  • About the Journal
  • Instructions for Authors
  • Journal Policies
  • Reprints and Permissions
  • Advertising
  • Privacy Policy: ADA Journals
  • Copyright Notice/Public Access Policy
  • Contact Us

Other ADA Resources

  • Diabetes
  • Clinical Diabetes
  • Diabetes Spectrum
  • Scientific Sessions Abstracts
  • Standards of Medical Care in Diabetes
  • BMJ Open - Diabetes Research & Care
  • Professional Books
  • Diabetes Forecast

 

  • DiabetesJournals.org
  • Diabetes Core Update
  • ADA's DiabetesPro
  • ADA Member Directory
  • Diabetes.org

© 2021 by the American Diabetes Association. Diabetes Care Print ISSN: 0149-5992, Online ISSN: 1935-5548.