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

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • 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
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • 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
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Genetics/Genomes/Proteomics/Metabolomics

Bilirubin as a Potential Causal Factor in Type 2 Diabetes Risk: A Mendelian Randomization Study

  1. Ali Abbasi1,2,3⇑,
  2. Petronella E. Deetman2,
  3. Eva Corpeleijn1,
  4. Ron T. Gansevoort2,
  5. Rijk O.B. Gans2,
  6. Hans L. Hillege1,
  7. Pim van der Harst4,5,6,
  8. Ronald P. Stolk1,
  9. Gerjan Navis2,
  10. Behrooz Z. Alizadeh1 and
  11. Stephan J.L. Bakker2
  1. 1Department of Epidemiology, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands
  2. 2Department of Internal Medicine, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands
  3. 3Medical Research Council Epidemiology Unit, University of Cambridge School of Clinical Medicine, Institute of Metabolic Science, Cambridge Biomedical Campus, Addenbrooke’s Hospital, Cambridge, U.K.
  4. 4Department of Cardiology, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands
  5. 5Department of Genetics, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands
  6. 6Durrer Center for Cardiogenetic Research, ICIN-Netherlands Heart Institute, Utrecht, the Netherlands
  1. Corresponding author: Ali Abbasi, ali.abbasi{at}mrc-epid.cam.ac.uk or a.abbasi{at}umcg.nl.
Diabetes 2015 Apr; 64(4): 1459-1469. https://doi.org/10.2337/db14-0228
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

Circulating bilirubin, a natural antioxidant, is associated with decreased risk of type 2 diabetes (T2D), but the nature of the relationship remains unknown. We performed Mendelian randomization in a prospective cohort of 3,381 participants free of diabetes at baseline (age 28–75 years; women 52.6%). We used rs6742078 located in the uridine diphosphate–glucuronosyltransferase locus as an instrumental variable (IV) to study a potential causal effect of serum total bilirubin level on T2D risk. T2D developed in a total of 210 participants (6.2%) during a median follow-up period of 7.8 years. In adjusted analyses, rs6742078, which explained 19.5% of bilirubin variation, was strongly associated with total bilirubin (a 0.68-SD increase in bilirubin levels per T allele; P < 1 × 10−122) and was also associated with T2D risk (odds ratio [OR] 0.69 [95% CI 0.54–0.90]; P = 0.006). Per 1-SD increase in log-transformed bilirubin levels, we observed a 25% (OR 0.75 [95% CI 0.62–0.92]; P = 0.004) lower risk of T2D. In Mendelian randomization analysis, the causal risk reduction for T2D was estimated to be 42% (causal OR for IV estimation per 1-SD increase in log-transformed bilirubin 0.58 [95% CI 0.39–0.84]; P = 0.005), which was comparable to the observational estimate (Durbin-Wu-Hausman χ2 test, P for difference = 0.19). These novel results provide evidence that an elevated bilirubin level is causally associated with the risk of T2D and support its role as a protective determinant.

Introduction

Experimental and human studies show that elevated bilirubin levels are associated with decreased risk of type 2 diabetes (T2D) and diabetes-related outcomes (1–6). Bilirubin might antagonize oxidative stress by acting as an antioxidant and cytoprotectant, which may have beneficial effects in diseases related to oxidative stress (7–10). For example, we have previously shown that elevated bilirubin levels are associated with a reduction in mortality in the carriers of iron overload genotypes, which are phenotypically associated with increased oxidative stress (11). It is, however, unclear whether the association between bilirubin and diabetes is free of unobserved confounders. Apart from confounding, bilirubin concentrations could be changed as a consequence of T2D pathology (i.e., reverse causality). Therefore, causal inference of the effect of bilirubin on T2D (i.e., bilirubin being an intermediate phenotype in the etiology of T2D) is uncertain in conventional observational studies (12).

One approach to overcome these aforementioned limitations is to implement an instrumental variable (IV) that represents bilirubin levels and allows for assessment of the relation to T2D that is robust to confounding and reverse causality biases (13,14). Recent genome-wide association studies (GWAS) have identified that rs6742078, the single nucleotide polymorphism (SNP) that is mapped to hepatic uridine diphosphate-glucuronyltransferase (UGT1A1), explains 18% of the variation in serum total bilirubin levels (15–17). The UGT1A1 gene encodes hepatic uridine diphosphate-glucuronyltransferase, which is a major enzyme in bilirubin conjugation and in the regulation of bilirubin levels (5,15–17). When UGT1A1 is downregulated, bilirubin levels subsequently rise. Given the fact that the inheritance of either a T or a G allele of rs6742078 is randomly defined prior to disease onset (i.e., at the time of conception, based on the laws of Mendel), the variation of rs6742078 is robust to unobserved confounding and reverse causality. Therefore, rs6742078 is a suitable candidate to be used as an IV representing bilirubin variation between individuals. The approach of using genetic instruments (i.e., rs6742078) to test the causal association of a given intermediate exposure (i.e., bilirubin) with a disease outcome (i.e., T2D) is called Mendelian randomization (14). We hypothesized that implementing rs6742078 as an IV could provide evidence that the previously reported association between serum bilirubin level and T2D is causal or is not.

We aimed to investigate whether elevated total bilirubin levels are causally associated with a lower risk of new-onset T2D. By implementation of Mendelian randomization, we examined whether rs6742078(T) is associated with serum total bilirubin levels and T2D risk in a population-based cohort study.

Research Design and Methods

Study Population and Design

The Prevention of Renal and Vascular End-stage Disease (PREVEND) study is a Dutch population-based cohort for which 8,592 participants (age range 28–75 years) from the city of Groningen, the Netherlands, were recruited between 1997 and 1998. We examined a random sample of 4,016 individuals of the total cohort in which genotyping was performed. Details of the PREVEND study have been published previously (18). Among 4,016 individuals, 367 participants were excluded because of not passing quality control of DNA analyses, 124 participants were excluded because of missing follow-up data or information about medication, and 144 individuals were excluded because of prevalent diabetes at baseline, leaving 3,381 participants eligible for analyses (Fig. 1).

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

Study population sample from the PREVEND study cohort. For quality control, samples were excluded (n = 367) due to call rates <0.95, sex mismatch, duplicate discordance, contamination, and genetic similarity (identity by state >0.2). Population stratification was assessed by principal component analysis over the sample correlation matrix, based on 16,842 independent (LD-pruned) SNPs. We excluded samples when they diverged from the mean with at least 3 SDs (z score >3) for the first five principal components.

The PREVEND study has been approved by the medical ethics committee, University Medical Center Groningen, and conformed to the principles outlined in the Declaration of Helsinki. All participants gave written informed consent.

Baseline Study Measurements

At baseline and three rounds of screening from 1997 to 1998 (i.e., baseline examination) until 1 January 2007 (i.e., the third examination), study participants underwent two outpatient visits to obtain clinical characteristics and to collect two 24-h urine samples on 2 consecutive days (18). Furthermore, information on medication use was substantiated using pharmacy-based data from all community pharmacies in the city of Groningen (19). Blood samples were obtained after overnight fasting and were stored at −80°C until analyses without previous thawing and refreezing (Supplementary Table 1). Total bilirubin level was measured by a colorimetric assay (2,4-dichloroaniline reaction; MEGA; Merck, Darmstadt, Germany).

Definition of Outcome

During follow-up, participants were classified as new-onset T2D case patients if the fasting plasma glucose level was ≥7.0 mmol/L (126 mg/dL), the random sample plasma glucose concentration was ≥11.1 mmol/L (200 mg/dL), they reported a physician diagnosis of T2D, or they received insulin or oral hypoglycemic agents based on a central pharmacy registration. We included case patients in whom T2D was diagnosed from the first 3 months after the baseline screening visits (1997–1998) until January 2007 (20,21).

Genotyping and Quality Control

Genotyping was performed using Illumina HumanCytoSNP-12 arrays. Genotypes of 3,649 individuals were available preimputation (see Fig. 1 for quality control of samples). SNPs were called using Illumina Genome Studio software. SNPs were excluded with a minor allele frequency <0.01, a call rate <0.95, or a deviation from Hardy-Weinberg equilibrium (P < 1 × 10−5). Genome-wide genotype imputation was performed using Beagle version 3.1.0. We used the National Center for Biotechnology Information build 36 of Phase II HapMap CEU data (release 22) as the reference panel (22). Given the expected strong relationship of rs6742078 (imputation accuracy score >0.9) with bilirubin levels (16,17) and the fact that rs6742078 is in high linkage disequilibrium (LD) (r2 = 0.88) with the functional UGT1A1*28 polymorphism (17), we decided that rs6742078 is a suitable IV for bilirubin in our study.

Statistical Analysis

We performed analyses in four steps. First, we examined the associations of bilirubin with participants’ characteristics, including age, sex, smoking status, family history of diabetes, and associated clinical variables. This was performed using linear regression. Next, we used logistic regression to examine the association of bilirubin with new-onset T2D (model A). This model estimates an observed effect of bilirubin on T2D (βbilirubin-T2D coefficient). The odds ratio (OR) and its corresponding 95% CIs for the risk of T2D was calculated per 1-SD increase in log-transformed bilirubin levels. Model A was adjusted for age, sex, family history of diabetes, BMI, hypertension, smoking status, and fasting glucose level (21). Second, we tested whether rs6742078 was associated with log-transformed bilirubin levels. We used linear regression under the assumption of an additive genetic model (model B), which is in agreement with previously published studies (17,23). Model B was adjusted for the same covariates as those included in model A (14). This analysis provided the observed effect estimate, which was expressed as the change in the number of SDs in log-transformed bilirubin per each copy increase in the number of T allele (range 0–2). To estimate the proportion of bilirubin variation explained by the IV of the rs6742078 genotypes, we calculated the difference between r2 values (i.e., explained variance) for model B without the genotype and for model B with the genotype. Third, logistic regression was used to evaluate the association between rs6742078 and T2D (model C). We calculated the OR for T2D per each copy increase in the rs6742078(T) allele. Model C was adjusted for the same covariates as those of model A. Fourth, we used the Wald-type method of IV estimation to test a potential causal relationship between elevated total bilirubin levels and decreased T2D risk. For the Wald-type method, the estimated risk is calculated using the following equation:Embedded ImageIn this equation, ORgenotype-T2D is the OR resulting from logistic regression between rs6742078 and T2D, representing the risk of T2D per each copy of the rs6742078(T) allele (model C). The βgenotype-bilirubin coefficient is obtained from the regression of bilirubin on the SNP per each copy of the T allele of rs6742078 (model B). To validate the causal estimate derived from the Wald-type method, we used two alternative methods, called the two-stage least squares (2SLS) method and the multiplicative generalized method of moments (MGMM) (Table 3). The first stage of the 2SLS method is to examine the observational association between rs6742078 and bilirubin by means of linear regression (model 1), saving the predicted values and the residuals. In the second stage, the predicted values of bilirubin from model 1 are used as covariates, with T2D as the dependent variable in a logistic regression model, providing the βIV2SLS estimate (14). To obtain the MGMM βIV, the Stata command “ivpois” was used (24). The MGMM does not require distributional assumptions about bilirubin and thus can provide consistent estimation of causal OR (24).

To examine whether the association between bilirubin and T2D has “endogeneity,” we compared the IV estimate with the observed estimate from standard logistic regression analysis (i.e., βbilirubin-T2D coefficient) using the Durbin-Wu-Hausman χ2 test, which assesses whether there is a difference between these two estimates (14,24). We calculated the power of our study to detect a causal effect for bilirubin, given our sample size and bilirubin variation explained by rs6742078 (17,25).

To account for potential confounding bias, we tested whether rs6742078 is associated with common diabetes risk factors, other variables (Table 2), and 24-h urine albumin excretion (UAE). The association with UAE was tested because of the enrichment of the PREVEND cohort with microalbuminuric individuals at baseline (18).

To test for pleiotropy, we performed logistic regression with bilirubin as an independent variable and T2D as a dependent variable, and saved the residuals. Next, we tested the correlation between the genotype and the residuals that are independent of total bilirubin levels (26). To further dissect potential pleiotropy, we also calculated IV estimates from two additional bilirubin SNPs, namely, rs4149056 in SLCO1B1 and rs16928809 in SLC22A18 (17), in the PREVEND study. Additionally, we extracted effect estimates for rs6742078, rs4149056, and rs16928809 from publically available GWAS data on bilirubin, T2D (DIAbetes Genetics Replication And Meta-analysis [DIAGRAM]), or glycemic traits (Meta-Analyses of Glucose and Insulin-related traits Consortium [MAGIC]) (17,27–29) (Supplementary Table 2), and we calculated a summary statistics genetic risk score for bilirubin, as described previously (27). We then evaluated the presence of heterogeneity (derived from the likelihood-based method [27]) on the causal effects of these three SNPs on T2D, glucose, insulin, and HOMA of insulin resistance (HOMA-IR) (27–29) (Table 3). The presence of heterogeneity indicates potential pleiotropy (30).

We also performed a series of secondary analyses. First, we repeated regression models for bilirubin levels and rs6742078 by using a weighted method to compensate for the enrichment of the PREVEND cohort with microalbuminuric individuals (18). Second, we added the interaction terms of bilirubin × sex and bilirubin × age to model A or rs6742078 × sex and rs6742078 × age to model C, to evaluate potential sex or age differences in the association of bilirubin or rs6742078 with T2D, respectively. Third, we examined whether the effect estimates for bilirubin or rs6742078 were consistent across age strata of people <50 or ≥50 years old. Fourth, we accounted for the time between obtaining the sample for the measurement of bilirubin (and determining the genotype) and the development of T2D by using a Cox regression in which we applied the approximate time to the incidence of T2D as the time scale in our analysis. The time to the incidence of T2D cannot be ascertained exactly, because we screened for incident T2D periodically. We tested the proportional hazards assumption indicating that the effect estimates do not vary over time. Fifth, a potential correlation between the dropout rates and the genotype was tested. Sixth, we examined the association of rs6742078 with prevalent or a combination of both prevalent and incident T2D. A single imputation and predictive mean matching method was applied for missing data (21). For most baseline clinical variables, <1% of data was missing, whereas this was up to 8% for self-reported variables such as family history of diabetes. All the statistical analyses were carried out using IBM SPSS version 19.0, Stata/IC version 11.2 (StataCorp, College Station, TX), or R version 2.15.1 (Vienna, Austria) for Windows (http://cran.r-project.org/).

Results

Baseline characteristics of the population are shown in Table 1. During a median follow-up period of 7.7 years (Q1–Q3 7.4–8.0 years), T2D developed in a total of 210 individuals (6.2%) (from 1997 to 1998 until 2007).

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

Baseline characteristics of the study participants

Association Between Bilirubin and T2D

We observed that total bilirubin levels were lower in individuals in whom T2D developed versus those in whom T2D did not develop (7.1 ± 3.1 vs. 7.7 ± 3.6 μmol/L, respectively; P < 0.001). Figure 2 depicts a graphical representation of the distribution of bilirubin with an incremental decrease in ORs for T2D over the range of bilirubin levels (2–35 μmol/L). The multivariable adjusted OR for T2D was 0.75 (95% CI 0.62–0.92), per 1-SD increase in log-transformed bilirubin levels (model A). The associations of bilirubin with baseline variables are shown in Table 2. When adding separately the variables, which were associated with bilirubin, to model A, corresponding ORs for T2D ranged from 0.78 (95% CI 0.67–0.92) to 0.83 (95% CI 0.69–0.94).

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

Histogram of serum total bilirubin levels with ORs for new-onset (i.e., incident) T2D over the range of bilirubin levels. The number of participants are represented by the white bars, which correspond to the right axis. The solid line denotes the ORs, which were centered on the median value of bilirubin.

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

Associations of baseline variables with serum bilirubin and rs6742078

Association Between rs6742078 and Bilirubin

We examined whether rs6742078, which has been previously discovered as the top-ranked SNP in a GWAS of total bilirubin (16,17), was associated with bilirubin levels in our cohort. In model B, the rs6742078*T allele was strongly associated with bilirubin levels (a 0.68-SD increase in log-transformed bilirubin levels per allele, P < 1 × 10−122). The proportion of bilirubin variation explained by rs6742078 was 19.5% (the difference between r2 = 0.070 for model B without the genotype and r2 = 0.265 for model B with the genotype). Using linear regression analysis, the rs6742078*T allele was associated with a 2.607 μmol/L higher total bilirubin level (Fig. 3, left panel).

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

Association of the rs6742078 genotypes (TT, GT, and GG) with levels of total bilirubin and incidence of T2D in the study population. The mean (SE) level of bilirubin for each genotype is represented by the gray bars in the left panel. The ORs (95% CIs) for the association of the genotypes with new-onset T2D are shown in the right panel.

Association Between rs6742078 and T2D

We observed that a one-copy increase in the number of the rs6742078*T allele was associated with T2D risk with an OR of 0.69 (95% CI 0.54–0.90; P = 0.006), which corresponds to a 31% decreased risk of T2D (Fig. 3, right panel). Adjustment for bilirubin level attenuated the association between rs6742078 and T2D risk (OR 0.80 [95% CI 0.60–1.06]; P = 0.13).

Causal Estimates for the Effect of Bilirubin on T2D risk

To estimate the causal effect of bilirubin on T2D, we used the Wald-type method, which assesses the ratio of the association of rs6742078 with T2D and that of rs6742078 with bilirubin levels (Fig. 4). Based on this, we calculated that the estimated causal effect of total bilirubin (per 1-SD increase in log-transformed bilirubin levels) corresponds to a 42% lower risk of T2D (causal ORIVestimation 0.58 [95% CI 0.39–0.84]; P = 0.005). Given a sample size of 3,381 individuals (including 210 incident cases of T2D) and 19.5% of bilirubin variation explained by rs6742078 in our study, we had 80% power at a significance level of 0.05 to estimate an OR of 0.64 per 1-SD increase in log-transformed bilirubin levels. We observed that the direction of the IV and observed estimates were similar with a slightly larger effect size for the IV estimate (Pendogeneity = 0.19). The use of two alternative methods for IV estimation provided evidence of causal effect, which was also comparable to the observed OR for T2D (Pendogeneity = 0.18 or 0.22; Table 3).

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

Mendelian randomization analysis for the association of levels of serum total bilirubin and the incidence of T2D. The observed effect of genotype on bilirubin levels (βgenotype-bilirubin) and the risk of T2D (βgenotype-T2D) was estimated per each copy of the T allele of rs6742078. The observed effect of bilirubin on T2D was estimated per 1-SD increase in log-transformed bilirubin levels. The estimated causal risk of T2D was calculated by the Wald-type method, which is a ratio of βgenotype-T2D to βgenotype-bilirubin, where risk = (1 − exp(β)) × 100.

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

Causal estimates for T2D and glycemic traits using Mendelian randomization analysis

Confounding and Pleiotropy

We found no evidence that rs6742078 is associated with potential confounders like various demographic, lifestyle, and clinical variables (Table 2). Thus, we assumed that the association between the IV and T2D was mainly free of confounding. To test for pleiotropy, we assessed the correlation between rs6742078 and the residuals from the logistic regression for T2D (as dependent variable) and bilirubin (as independent variable). There was no evidence for pleiotropy (P for correlation = 0.11). We consistently observed no heterogeneity for the causal effects of bilirubin on T2D risk or glycemic traits when we used a summary statistics genetic score including rs6742078, rs4149056, and rs16928809 (Table 3). Details of the association of individual SNPs with bilirubin, T2D, and glycemic traits are given in Supplementary Table 2. In the PREVEND study, while rs6742078 explained 19.5% of bilirubin variation, rs4149056 and rs16928809 explained only 0.18% and 0.15% of the variance, respectively. This limits the use of these SNPs as a single IV. For rs4149056 and rs16928809, we estimated the causal ORs of 0.58 (95% CI 0.05–6.08) and 0.17 (95% CI 0.02–2.16) in the PREVEND study. These estimates lack power for the detection of associations as these two SNPs explained only a very small amount of bilirubin variation and therefore are to be considered weak IVs.

Secondary Analyses

In subsequent analyses, we repeated the analyses in models A and C, when weighted for individuals with a mean urine albumin level of >10 mg/L. The results of these complex design analyses were similar to those of the main analyses, with an OR of 0.79 (95% CI 0.66–0.95; P = 0.01) per 1-SD increase in log-transformed bilirubin, and an OR of 0.59 (95% CI 0.47–0.97; P < 0.001) for the rs6742078*T allele. Second, we found no evidence for interactions of sex and age with bilirubin and rs6742078 (P for sex interaction = 0.49 and P for age interaction = 0.38 with bilirubin; P for sex interaction = 0.96 and P for age interaction = 0.25 with rs6742078) for the risk of T2D. Third, we observed that both observational and IV estimates were age consistent, with ORs of 0.69 (95% CI 0.51–1.04) and 0.45 (95% CI 0.21–0.94) in people <50 years of age, and 0.86 (95% CI 0.70–1.04) and 0.55 (95% CI 0.34–0.90) in people ≥50 years of age, respectively. Fourth, the assumption of proportional hazards was met for the Cox regression model, indicating there is no evidence of time-varying associations of bilirubin and rs6742078 with T2D (P = 0.38 and P = 0.54, respectively). Fifth, we observed no evidence of an association between the dropout rates (82.9% and 71.3% of participants, respectively, who underwent the second and the third round) and the genotype (P for correlation = 0.15). Sixth, we found an OR of 0.93 (95% CI 0.69–1.26) per each T allele of rs6742078 for prevalent T2D; and when we performed a combined analysis of both prevalent and incident T2D, our findings remained materially unchanged (OR 0.77 [95% CI 0.64–0.93]). In further analysis, we tested for another UGT1A1 SNP, rs887829 (in perfect LD with rs6742078; r2 = 1.0), which has been shown to be strongly associated with bilirubin levels in GWAS data (31). In our cohort, the rs887829*T allele was associated with a decreased risk of T2D (OR 0.70 [95% CI 0.54–0.89]), which was comparable to that of the rs6742078*T allele. Finally, we observed no evidence of systematic differences between the participants who were included in the main analysis compared with the rest of the cohort participants without diabetes at baseline (data not shown).

Discussion

In this prospective cohort study, we investigated whether we could provide evidence that bilirubin is causally associated with the risk of new-onset T2D. We demonstrated that the rs6742078*T allele in the UGT1A1 locus is strongly associated with elevated levels of endogenous bilirubin and is also associated with new-onset T2D. Several studies (1–6) have previously demonstrated associations between bilirubin and T2D in animal models and in humans. However, confounding and reverse causality might have contributed to the implication of causal relationships in observational studies. To our knowledge, this is the first study to report a potential causal association between elevated total bilirubin levels and decreased T2D risk that is likely free of potential confounders using Mendelian randomization.

Accumulating evidence shows that bilirubin has powerful antioxidant properties (7–10). Bilirubin may compensate for the oxidative stress that might be an important factor in the pathophysiology of diabetes (32–36). Oxidative stress is consistently related to increased plasma glucose levels and may contribute to the development of diabetes (37), and to microvascular and macrovascular complications of diabetes (38,39). Bilirubin has been shown to prevent oxidative stress in a number of diseases, including atherosclerosis, cancer, and diabetic nephropathy (40,41). In mice, biliverdin, a precursor of bilirubin, protects against the deterioration of glucose tolerance (3), while increased levels of bilirubin reduced streptozotocin-induced pancreatic β-cell damage through attenuating oxidative stress (42) or increasing insulin sensitivity in mice (43). In line with experimental studies, previous observational studies (2,44) have shown that elevated bilirubin level is associated with a 26–31% decreased risk of diabetes in humans.

In previous studies (5,15–17), variations in the UGT1A1 and SLCO1B1 loci seemed to explain part of interindividual differences in bilirubin levels. Johnson et al. (17) found that rs6742078, very close to the UGT1A1*28 TATA box polymorphism, accounted for ∼18% of the variation in total bilirubin levels (P < 5.0 × 10−324). Less powerful, but still present, was the association with rs4149056, a functional SNP of the SLCO1B1 locus, which accounted for 0.6% of bilirubin variation (17). In our cohort study, we observed that 19.5% of bilirubin variation was explained by the UGT1A1 rs6742078 SNP. It is likely that the association of rs6742078 with bilirubin can be explained to a large extent by the functional UGT1A1*28 genetic variant, which has been found to be in high LD (r2 = 0.88) with rs6742078 (17). Therefore, the contribution of genetic variation in UGT1A1 to the risk of T2D could be greater than estimated from the associations observed for rs6742078. The genetically determined component of protective effects may also vary globally as the TA repeat in the UGT1A1 promoter has other, less common functional variants (e.g., the UGT1A1*6 variant, which is more often observed in Asians with Japanese or Chinese ancestry) (45).

From phenotypic perspective, the Framingham Heart Study (5) has found that the UGT1A1*28 variant was associated with both bilirubin level and cardiovascular disease (CVD). However, the relationships of the UGT1A1 SNPs with CVD, glycemic traits, and diabetes have been inconsistent, possibly due to differences in design and populations across studies (23). In the current study, we found that the rs6742078 SNP was associated with T2D risk. We tested whether rs6742078 was associated with glycemic traits that were recently evaluated in meta-analysis of 21 GWA studies (29). rs6742078 was associated with fasting glucose level (β = −0.009, P = 0.02), insulin level (β = −0.009, P = 0.03), and HOMA-IR (β = −0.009, P = 0.03) based on a priori selected single-test association analysis in published GWAS data (29). In a recent GWAS meta-analysis (28) (12,171 case subjects and 56,862 control subjects), no evidence of an association between rs6742078 and T2D has been reported (P = 0.1). The lack of association between rs6742078 and T2D might be due to the heterogeneity in the T2D phenotype (28,46), as T2D cases have been ascertained using different sources (28). Of note, we identified a similar pattern of absence of an association in the GWAS meta-analysis, while another prospective Mendelian randomization study (47) reported the presence of such an association. In the latter study, the presence of evidence for the association between two SNPs in the gene encoding sex hormone-binding globulin and T2D could not be confirmed in the case-control GWAS meta-analysis. Case-control studies typically included both prevalent and incident cases, which differ with respect to the disease stage, diagnostic criteria, and duration of diabetes (46). Also, these discrepancies may be partly explained by differences in data quality control and preanalysis preparations (28,47). A prospective cohort design and its accompanying statistical methods will reduce potential biases such as incidence-prevalence bias and nonresponse bias (46). Although we estimated the risk of future T2D, there still might be a chance of a type I error in our analysis. The current results, together with the associations of rs6742078 with glycemic traits in the meta-analysis of GWAS, support the idea that elevated bilirubin levels may lead to decreased risk of insulin resistance and T2D. Our finding remains to be confirmed for the potential role of bilirubin in the trajectories of glycemic traits and the development of T2D by prospective cohort studies.

In our study, we used a strong IV, rs6742078, which satisfied the requirements for Mendelian randomization (48), to examine the potential causal relationship between bilirubin and T2D in a large population-based study. Using this approach, we directly investigated a relationship between genotype and phenotype independent of potential confounders. Other bilirubin SNPs, rs4149056 and rs16928809 (17), may not be selected as a single IV, given the fact that the strength of these instruments is limited (25) because of the small amount of bilirubin variation that is explained by these two SNPs (0.18% and 0.15% in the PREVEND study, respectively).

Our findings may have implications for the prevention or prediction of T2D and its complications. For example, Gilbert syndrome, which is characterized by moderate hyperbilirubinemia, exists in 5–10% of the general population (5,16,49). It has been shown that diabetic patients with Gilbert syndrome have reduced markers of oxidative stress and decreased risk of diabetic nephropathy and CVD (5,49). One may speculate that the modulation of bilirubin levels (e.g., using medical interventions that inhibit UGT1A1) might have protective effects against the risk of T2D (45). A search for the sources of increase in bilirubin levels or its antioxidative capacity (e.g., effective therapeutics with minimal adverse effects) merits further investigations.

One major concern remaining in the context of Mendelian randomization is that we cannot completely control for the possibility of reverse causality, unobserved confounding, and pleiotropy (14,50). The functional UGT1A1*28 variant (in high LD with rs6742078), which is associated with a 10–35% reduction in UGT1A1 activity, is mainly involved in glucuronidation of bilirubin, but also of other compounds (e.g., estrogen) and drugs (e.g., irinotecan) (45,51). Yet, the effects of the genotype on phenotypic traits other than bilirubin and their consequences on disease outcomes (e.g., cancer) remain inconclusive. An extensive knowledge of gene function and associated biological processes is needed to better understand the extent to which such genetic variation in the human genome contributes to lifelong protection against T2D. Our cohort is predominantly comprised of white adults of Dutch descent, and it is therefore unknown whether the findings of our study can be generalized to other populations (18). The PREVEND cohort was enriched for individuals with microalbuminuria at baseline. However, in a secondary analysis, we used a weighted method to compensate for this, and this did not affect our results.

In summary, endogenous bilirubin and the most highly associated SNP, rs6742078 in the UGT1A1 locus, are associated with the risk of T2D. We observed consistent effect estimates derived from alternative genetic instruments, rs4149056 and rs16928809. However, these two SNPs explain far less of bilirubin variation, and therefore the estimates derived from these individual instruments should be interpreted cautiously. Our findings provide evidence that lifelong genetically elevated bilirubin level are likely to be causally protective against the development of T2D. Further studies are warranted to validate this finding and to elucidate the exact underlying mechanisms of action.

Article Information

Acknowledgments. The authors thank Professor Dr. L.T.W. de Jong-van den Berg and Dr. S.T. Visser from the Department of Social Pharmacy, Pharmacoepidemiology, and Pharmacotherapy, Groningen University Institute for Drug Exploration, University of Groningen, University Medical Center Groningen, for providing the data on pharmacy-registered use of insulin or oral hypoglycemic agents.

Funding. This work was supported by the Netherlands Heart Foundation, Dutch Diabetes Research Foundation, Dutch Kidney Foundation, the Netherlands Organization for Scientific Research project (NWO), and the Medical Research Council UK (grant no. MC_U106179471). A.A. is supported by a Rubicon grant from the NWO (project no. 825.13.004).

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

None of the study sponsors had a role in the study design; data collection, analysis, and interpretation; report writing; or the decision to submit the report for publication.

Author Contributions. A.A., P.E.D., and S.J.L.B. helped to conceive and design the study, analyze the data, write the first draft of the manuscript, and agreed with the manuscript results and conclusions. E.C., R.T.G., R.O.B.G., H.L.H., P.v.d.H., and G.N. contributed to the writing of the manuscript, and agreed with the manuscript results and conclusions. R.P.S. helped to conceive and design the study, contributed to the writing of the manuscript, and agreed with the manuscript results and conclusions. B.Z.A. helped to conceive and design the study, analyzed the data, contributed to the writing of the manuscript, and agreed with the manuscript results and conclusions. A.A., B.Z.A., and S.J.L.B. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Footnotes

  • This article contains Supplementary Data online at http://diabetes.diabetesjournals.org/lookup/suppl/doi:10.2337/db14-0228/-/DC1.

  • Received February 12, 2014.
  • Accepted October 25, 2014.
  • © 2015 by the American Diabetes Association. 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.

References

  1. ↵
    1. Lin R,
    2. Wang Y,
    3. Wang Y, et al
    . Common variants of four bilirubin metabolism genes and their association with serum bilirubin and coronary artery disease in Chinese Han population. Pharmacogenet Genomics 2009;19:310–318pmid:19238116
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Cheriyath P,
    2. Gorrepati VS,
    3. Peters I, et al
    . High total bilirubin as a protective factor for diabetes mellitus: an analysis of NHANES data from 1999 - 2006. J Clin Med Res 2010;2:201–206pmid:21629541
    OpenUrlPubMed
  3. ↵
    1. Ikeda N,
    2. Inoguchi T,
    3. Sonoda N, et al
    . Biliverdin protects against the deterioration of glucose tolerance in db/db mice. Diabetologia 2011;54:2183–2191pmid:21614569
    OpenUrlCrossRefPubMed
    1. Horsfall LJ,
    2. Nazareth I,
    3. Petersen I
    . Cardiovascular events as a function of serum bilirubin levels in a large, statin-treated cohort. Circulation 2012;126:2556–2564pmid:23110860
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Lin JP,
    2. O’Donnell CJ,
    3. Schwaiger JP, et al
    . Association between the UGT1A1*28 allele, bilirubin levels, and coronary heart disease in the Framingham Heart Study. Circulation 2006;114:1476–1481pmid:17000907
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Jenko-Pražnikar Z,
    2. Petelin A,
    3. Jurdana M,
    4. Žiberna L
    . Serum bilirubin levels are lower in overweight asymptomatic middle-aged adults: an early indicator of metabolic syndrome? Metabolism 2013;62:976–985pmid:23414908
    OpenUrlCrossRefPubMed
  6. ↵
    1. Duann P,
    2. Lianos EA
    . GEC-targeted HO-1 expression reduces proteinuria in glomerular immune injury. Am J Physiol Renal Physiol 2009;297:F629–F638pmid:19587144
    OpenUrlAbstract/FREE Full Text
    1. Abraham NG,
    2. Asija A,
    3. Drummond G,
    4. Peterson S
    . Heme oxygenase -1 gene therapy: recent advances and therapeutic applications. Curr Gene Ther 2007;7:89–108pmid:17430129
    OpenUrlCrossRefPubMedWeb of Science
    1. Baranano DE,
    2. Rao M,
    3. Ferris CD,
    4. Snyder SH
    . Biliverdin reductase: a major physiologic cytoprotectant. Proc Natl Acad Sci USA 2002;99:16093–16098pmid:12456881
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Stocker R,
    2. Yamamoto Y,
    3. McDonagh AF,
    4. Glazer AN,
    5. Ames BN
    . Bilirubin is an antioxidant of possible physiological importance. Science 1987;235:1043–1046pmid:3029864
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Alizadeh BZ,
    2. Njajou OT,
    3. Houwing-Duistermaat JJ, et al
    . Does bilirubin protect against hemochromatosis gene (HFE) related mortality? Am J Med Genet A 2004;129A:39–43pmid:15266614
    OpenUrlPubMed
  9. ↵
    1. Iribarren C,
    2. Reed DM,
    3. Chen R,
    4. Yano K,
    5. Dwyer JH
    . Low serum cholesterol and mortality. Which is the cause and which is the effect? Circulation 1995;92:2396–2403pmid:7586337
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Thanassoulis G,
    2. O’Donnell CJ
    . Mendelian randomization: nature’s randomized trial in the post-genome era. JAMA 2009;301:2386–2388pmid:19509388
    OpenUrlCrossRefPubMedWeb of Science
  11. ↵
    1. Lawlor DA,
    2. Harbord RM,
    3. Sterne JA,
    4. Timpson N,
    5. Davey Smith G
    . Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med 2008;27:1133–1163pmid:17886233
    OpenUrlCrossRefPubMed
  12. ↵
    1. Buch S,
    2. Schafmayer C,
    3. Völzke H, et al
    . Loci from a genome-wide analysis of bilirubin levels are associated with gallstone risk and composition. Gastroenterology 2010;139:1942–1951
    OpenUrlCrossRefPubMed
  13. ↵
    1. Lin JP,
    2. Vitek L,
    3. Schwertner HA
    . Serum bilirubin and genes controlling bilirubin concentrations as biomarkers for cardiovascular disease. Clin Chem 2010;56:1535–1543pmid:20693308
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Johnson AD,
    2. Kavousi M,
    3. Smith AV, et al
    . Genome-wide association meta-analysis for total serum bilirubin levels. Hum Mol Genet 2009;18:2700–2710pmid:19414484
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Lambers Heerspink HJ,
    2. Brantsma AH,
    3. de Zeeuw D,
    4. Bakker SJ,
    5. de Jong PE,
    6. Gansevoort RT; PREVEND Study Group
    . Albuminuria assessed from first-morning-void urine samples versus 24-hour urine collections as a predictor of cardiovascular morbidity and mortality. Am J Epidemiol 2008;168:897–905pmid:18775924
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Monster TB,
    2. Janssen WM,
    3. de Jong PE,
    4. de Jong-van den Berg LT; PREVEND Study Group Prevention of REnal and Vascular ENT Stage Disease
    . Pharmacy data in epidemiological studies: an easy to obtain and reliable tool. Pharmacoepidemiol Drug Saf 2002;11:379–384pmid:12271879
    OpenUrlCrossRefPubMedWeb of Science
  17. ↵
    1. Abbasi A,
    2. Corpeleijn E,
    3. Meijer E, et al
    . Sex differences in the association between plasma copeptin and incident type 2 diabetes: the Prevention of Renal and Vascular Endstage Disease (PREVEND) study. Diabetologia 2012;55:1963–1970pmid:22526609
    OpenUrlCrossRefPubMed
  18. ↵
    1. Abbasi A,
    2. Peelen LM,
    3. Corpeleijn E, et al
    . Prediction models for risk of developing type 2 diabetes: systematic literature search and independent external validation study. BMJ 2012;345:e5900pmid:22990994
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Verweij N,
    2. Mahmud H,
    3. Mateo Leach I, et al
    . Genome-wide association study on plasma levels of midregional-proadrenomedullin and C-terminal-pro-endothelin-1. Hypertension 2013;61:602–608pmid:23381795
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Stender S,
    2. Frikke-Schmidt R,
    3. Nordestgaard BG,
    4. Grande P,
    5. Tybjaerg-Hansen A
    . Genetically elevated bilirubin and risk of ischaemic heart disease: three Mendelian randomization studies and a meta-analysis. J Intern Med 2013;273:59–68pmid:22805420
    OpenUrlCrossRefPubMed
  21. ↵
    1. Palmer TM,
    2. Sterne JA,
    3. Harbord RM, et al
    . Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses. Am J Epidemiol 2011;173:1392–1403pmid:21555716
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Burgess S
    . Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol 2014;43:922–929pmid:24608958
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Li N,
    2. van der Sijde MR,
    3. Bakker SJ, et al, .; LifeLines Cohort Study Group
    . Pleiotropic effects of lipid genes on plasma glucose, HbA1c, and HOMA-IR levels. Diabetes 2014;63:3149–3158pmid:24722249
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Burgess S,
    2. Butterworth A,
    3. Thompson SG
    . Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol 2013;37:658–665pmid:24114802
    OpenUrlCrossRefPubMed
  25. ↵
    1. Morris AP,
    2. Voight BF,
    3. Teslovich TM, et al, .; Wellcome Trust Case Control Consortium, ; Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) Investigators, ; Genetic Investigation of ANthropometric Traits (GIANT) Consortium, ; Asian Genetic Epidemiology Network–Type 2 Diabetes (AGEN-T2D) Consortium, ; South Asian Type 2 Diabetes (SAT2D) Consortium, ; DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium
    . Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet 2012;44:981–990pmid:22885922
    OpenUrlCrossRefPubMed
  26. ↵
    1. Dupuis J,
    2. Langenberg C,
    3. Prokopenko I, et al, .; DIAGRAM Consortium, ; GIANT Consortium, ; Global BPgen Consortium, ; Anders Hamsten on behalf of Procardis Consortium, ; MAGIC Investigators
    . New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 2010;42:105–116pmid:20081858
    OpenUrlCrossRefPubMedWeb of Science
  27. ↵
    1. Burgess S,
    2. Thompson SG
    . Use of allele scores as instrumental variables for Mendelian randomization. Int J Epidemiol 2013;42:1134–1144pmid:24062299
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Sanna S,
    2. Busonero F,
    3. Maschio A, et al
    . Common variants in the SLCO1B3 locus are associated with bilirubin levels and unconjugated hyperbilirubinemia. Hum Mol Genet 2009;18:2711–2718pmid:19419973
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Van Campenhout A,
    2. Van Campenhout C,
    3. Lagrou AR, et al
    . Impact of diabetes mellitus on the relationships between iron-, inflammatory- and oxidative stress status. Diabetes Metab Res Rev 2006;22:444–454pmid:16506275
    OpenUrlCrossRefPubMed
    1. Robertson RP,
    2. Harmon JS
    . Diabetes, glucose toxicity, and oxidative stress: a case of double jeopardy for the pancreatic islet beta cell. Free Radic Biol Med 2006;41:177–184pmid:16814095
    OpenUrlCrossRefPubMedWeb of Science
    1. Seppen J,
    2. Bosma P
    . Bilirubin, the gold within. Circulation 2012;126:2547–2549pmid:23110861
    OpenUrlFREE Full Text
    1. Wellen KE,
    2. Hotamisligil GS
    . Inflammation, stress, and diabetes. J Clin Invest 2005;115:1111–1119pmid:15864338
    OpenUrlCrossRefPubMedWeb of Science
  30. ↵
    1. Li N,
    2. Frigerio F,
    3. Maechler P
    . The sensitivity of pancreatic beta-cells to mitochondrial injuries triggered by lipotoxicity and oxidative stress. Biochem Soc Trans 2008;36:930–934pmid:18793163
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    1. Yamagishi S
    . [Role of advanced glycation end products (AGE) and soluble receptor for AGE (sRAGE) in vascular complications in diabetes]. Nihon Rinsho 2012;70(Suppl. 5):243–247pmid:23156401
    OpenUrlPubMed
  32. ↵
    1. Baynes JW,
    2. Thorpe SR
    . Role of oxidative stress in diabetic complications: a new perspective on an old paradigm. Diabetes 1999;48:1–9pmid:9892215
    OpenUrlAbstract
  33. ↵
    1. Finkel T,
    2. Holbrook NJ
    . Oxidants, oxidative stress and the biology of ageing. Nature 2000;408:239–247pmid:11089981
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    1. Vítek L
    . The role of bilirubin in diabetes, metabolic syndrome, and cardiovascular diseases. Front Pharmacol 2012;3:55pmid:22493581
    OpenUrlCrossRefPubMed
  35. ↵
    1. Riphagen IJ,
    2. Deetman PE,
    3. Bakker SJ, et al
    . Bilirubin and progression of nephropathy in type 2 diabetes: a post hoc analysis of RENAAL with independent replication in IDNT. Diabetes 2014;63:2845–2853pmid:24677717
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Fu YY,
    2. Kang KJ,
    3. Ahn JM, et al
    . Hyperbilirubinemia reduces the streptozotocin-induced pancreatic damage through attenuating the oxidative stress in the Gunn rat. Tohoku J Exp Med 2010;222:265–273pmid:21139377
    OpenUrlCrossRefPubMed
  37. ↵
    1. Dong H,
    2. Huang H,
    3. Yun X, et al
    . Bilirubin increases insulin sensitivity in leptin-receptor deficient and diet-induced obese mice through suppression of ER stress and chronic inflammation. Endocrinology 2014;155:818–828pmid:24424052
    OpenUrlCrossRefPubMedWeb of Science
  38. ↵
    1. Jung CH,
    2. Lee MJ,
    3. Kang YM, et al
    . Higher serum bilirubin level as a protective factor for the development of diabetes in healthy Korean men: a 4 year retrospective longitudinal study. Metabolism 2014;63:87–93pmid:24140099
    OpenUrlCrossRefPubMed
  39. ↵
    1. Erlinger S,
    2. Arias IM,
    3. Dhumeaux D
    . Inherited disorders of bilirubin transport and conjugation: new insights into molecular mechanisms and consequences. Gastroenterology 2014;146:1625–1638pmid:24704527
    OpenUrlCrossRefPubMed
  40. ↵
    1. Manolio TA,
    2. Bailey-Wilson JE,
    3. Collins FS
    . Genes, environment and the value of prospective cohort studies. Nat Rev Genet 2006;7:812–820pmid:16983377
    OpenUrlCrossRefPubMedWeb of Science
  41. ↵
    1. Ding EL,
    2. Song Y,
    3. Manson JE, et al
    . Sex hormone-binding globulin and risk of type 2 diabetes in women and men. N Engl J Med 2009;361:1152–1163pmid:19657112
    OpenUrlCrossRefPubMedWeb of Science
  42. ↵
    1. Pierce BL,
    2. Ahsan H,
    3. Vanderweele TJ
    . Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol 2011;40:740–752pmid:20813862
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Maruhashi T,
    2. Soga J,
    3. Fujimura N, et al
    . Hyperbilirubinemia, augmentation of endothelial function, and decrease in oxidative stress in Gilbert syndrome. Circulation 2012;126:598–603pmid:22773454
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Pfister R,
    2. Sharp S,
    3. Luben R, et al
    . Mendelian randomization study of B-type natriuretic peptide and type 2 diabetes: evidence of causal association from population studies. PLoS Med 2011;8:e1001112pmid:22039354
    OpenUrlCrossRefPubMed
  45. ↵
    1. Stender S,
    2. Frikke-Schmidt R,
    3. Nordestgaard BG,
    4. Tybjærg-Hansen A
    . Extreme bilirubin levels as a causal risk factor for symptomatic gallstone disease. JAMA Intern Med 2013;173:1222–1228pmid:23753274
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Diabetes: 64 (4)

In this Issue

April 2015, 64(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
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

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

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.
Bilirubin as a Potential Causal Factor in Type 2 Diabetes Risk: A Mendelian Randomization Study
(Your Name) has forwarded a page to you from Diabetes
(Your Name) thought you would like to see this page from the Diabetes web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Bilirubin as a Potential Causal Factor in Type 2 Diabetes Risk: A Mendelian Randomization Study
Ali Abbasi, Petronella E. Deetman, Eva Corpeleijn, Ron T. Gansevoort, Rijk O.B. Gans, Hans L. Hillege, Pim van der Harst, Ronald P. Stolk, Gerjan Navis, Behrooz Z. Alizadeh, Stephan J.L. Bakker
Diabetes Apr 2015, 64 (4) 1459-1469; DOI: 10.2337/db14-0228

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

Bilirubin as a Potential Causal Factor in Type 2 Diabetes Risk: A Mendelian Randomization Study
Ali Abbasi, Petronella E. Deetman, Eva Corpeleijn, Ron T. Gansevoort, Rijk O.B. Gans, Hans L. Hillege, Pim van der Harst, Ronald P. Stolk, Gerjan Navis, Behrooz Z. Alizadeh, Stephan J.L. Bakker
Diabetes Apr 2015, 64 (4) 1459-1469; DOI: 10.2337/db14-0228
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
    • Discussion
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Low-Frequency Genetic Variant in the Hepatic Glucokinase Gene Is Associated With Type 2 Diabetes and Insulin Resistance in Chinese Population
  • Baseline Assessment of Circulating MicroRNAs Near Diagnosis of Type 1 Diabetes Predicts Future Stimulated Insulin Secretion
  • Polycystic Ovary Syndrome and Risk of Type 2 Diabetes, Coronary Heart Disease, and Stroke
Show more Genetics/Genomes/Proteomics/Metabolomics

Similar Articles

Navigate

  • Current Issue
  • Online Ahead of Print
  • Scientific Sessions Abstracts
  • Collections
  • 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 Care
  • 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 Print ISSN: 0012-1797, Online ISSN: 1939-327X.