PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosmedPLoS MedicineSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Med. May 2012; 9(5): e1001212.
Published online May 1, 2012. doi:  10.1371/journal.pmed.1001212
PMCID: PMC3341326
The Effect of Elevated Body Mass Index on Ischemic Heart Disease Risk: Causal Estimates from a Mendelian Randomisation Approach
Børge G. Nordestgaard,#1,2,3,4* Tom M. Palmer,#5,6 Marianne Benn,1,2,4 Jeppe Zacho,1,2,4 Anne Tybjærg-Hansen,2,3,4,7 George Davey Smith,5,6 and Nicholas J. Timpson5,6*
1Department of Clinical Biochemistry, Herlev Hospital, Copenhagen University Hospital, Copenhagen, Denmark
2The Copenhagen General Population Study, Herlev Hospital, Copenhagen University Hospital, Copenhagen, Denmark
3The Copenhagen City Heart Study, Bispebjerg Hospital, Copenhagen University Hospital, Copenhagen, Denmark
4Faculty of Health Sciences, University of Copenhagen, Copenhagen, Denmark
5MRC Centre for Causal Analyses in Translational Epidemiology, University of Bristol, Bristol, United Kingdom
6School of Social and Community Medicine, University of Bristol, Bristol, United Kingdom
7Department of Clinical Biochemistry, Rigshospitalet, Copenhagen University Hospital, Copenhagen, Denmark
Cosetta Minelli, Academic Editor
European Academy, Italy
#Contributed equally.
* E-mail: Boerge.Nordestgaard/at/regionh.dk (BGN); n.j.timpson/at/bris.ac.uk (NJT)
Conceived and designed the experiments: BGN NJT GDS. Analyzed the data: NJT TMP MB. Wrote the first draft of the manuscript: NJT. Contributed to the writing of the manuscript: NJT TMP BGN GDS JZ ATH MB. ICMJE criteria for authorship read and met: NJT TMP BGN GDS JZ ATH MB. Agree with manuscript results and conclusions: NJT TMP BGN GDS JZ ATH MB. Enrolled patients: BGN JZ ATH.
Received September 2, 2011; Accepted March 22, 2012.
Background
Adiposity, assessed as elevated body mass index (BMI), is associated with increased risk of ischemic heart disease (IHD); however, whether this is causal is unknown. We tested the hypothesis that positive observational associations between BMI and IHD are causal.
Methods and Findings
In 75,627 individuals taken from two population-based and one case-control study in Copenhagen, we measured BMI, ascertained 11,056 IHD events, and genotyped FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238). Using genotypes as a combined allele score in instrumental variable analyses, the causal odds ratio (OR) between BMI and IHD was estimated and compared with observational estimates. The allele score-BMI and the allele score-IHD associations used to estimate the causal OR were also calculated individually. In observational analyses the OR for IHD was 1.26 (95% CI 1.19–1.34) for every 4 kg/m2 increase in BMI. A one-unit allele score increase associated with a 0.28 kg/m2 (95 CI% 0.20–0.36) increase in BMI and an OR for IHD of 1.03 (95% CI 1.01–1.05) (corresponding to an average 1.68 kg/m2 BMI increase and 18% increase in the odds of IHD for those carrying all six BMI increasing alleles). In instrumental variable analysis using the same allele score the causal IHD OR for a 4 kg/m2 increase in BMI was 1.52 (95% CI 1.12–2.05).
Conclusions
For every 4 kg/m2 increase in BMI, observational estimates suggested a 26% increase in odds for IHD while causal estimates suggested a 52% increase. These data add evidence to support a causal link between increased BMI and IHD risk, though the mechanism may ultimately be through intermediate factors like hypertension, dyslipidemia, and type 2 diabetes. This work has important policy implications for public health, given the continuous nature of the BMI-IHD association and the modifiable nature of BMI. This analysis demonstrates the value of observational studies and their ability to provide unbiased results through inclusion of genetic data avoiding confounding, reverse causation, and bias.
Please see later in the article for the Editors' Summary
Background
Ischemic heart disease (IHD; also known as coronary heart disease) is the leading cause of death among adults in developed countries. In the US alone, IHD kills nearly half a million people every year. With age, fatty deposits (atherosclerotic plaques) build up in the walls of the coronary arteries, the blood vessels that supply the heart with oxygen and nutrients. The resultant reduction in the heart's blood supply causes shortness of breath, angina (chest pains that are usually relieved by rest), and potentially fatal heart attacks (myocardial infarctions). Risk factors for IHD include smoking, high blood pressure (hypertension), abnormal amounts of cholesterol and other fat in the blood (dyslipidemia), type 2 diabetes, and being overweight or obese (having excess body fat). Treatments for IHD include lifestyle changes (for example, losing weight) and medications that lower blood pressure and blood cholesterol levels. The narrowed arteries can also be widened using a device called a stent or surgically bypassed.
Why Was This Study Done?
Prospective observational studies have shown an association between a high body mass index (BMI, a measure of body fat that is calculated by dividing a person's weight in kilograms by their height in meters squared; a BMI greater than 30 kg/m2 indicates obesity) and an increased risk of IHD. Observational studies, which ask whether people who are exposed to a suspected risk factor develop a specific disease more often than people who are not exposed to the risk factor, cannot prove, however, that changes in BMI/adiposity cause IHD. Obese individuals may share other characteristics that cause both IHD and obesity (confounding) or, rather than obesity causing IHD, IHD may cause obesity (reverse causation). Here, the researchers use “Mendelian randomization” to examine whether elevations in BMI across the lifecourse have a causal impact on IHD risk. Three common genetic variants—FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238)—which have the largest single genetic variant associations with BMI were used in this study. Given that gene variants are inherited essentially randomly with respect to conventional confounding factors and are not subject reverse causation, use of these as instruments (or proxy measures) for variation in BMI as a risk factor (as opposed to measuring BMI directly) allows researchers to comment on whether obesity is causally involved in IHD.
What Did the Researchers Do and Find?
The researchers analyzed data from two population-based studies in which adults were physically examined and answered a lifestyle questionnaire before being followed to see how many developed IDH. They also analyzed data from a case-control study on IDH (in a case-control study, people with a disease are matched with similar people without the disease and the occurrence of risk factors in the patients and controls is compared). Overall, the researchers measured the BMI of 75,627 white individuals, among whom 11,056 already had IDH or developed it, and determined which of the BMI-increasing genetic variants each participant carried. On the basis of the observational data, every 4 kg/m2 increase in BMI increased the odds of IDH by 26% (an odds ratio of 1.26). Using a score derived from the combination of the three genetic variants, the researchers confirmed an association between each BMI increasing allele and both BMI (as expected) and IHD (0.28 kg/m2 and an odds ratio for IHD of 1.03, respectively). On average, compared to people carrying no BMI-increasing gene variants, people carrying six BMI-increasing gene variants had a 1.68 kg/m2 increase in BMI and an 18% increase in IHD risk. To extend this and to essentially reassess the original, observational, relationship between BMI and IHD risk, an “instrumental variable analysis” was used to examine the causal effect of a lifetime change in BMI on the risk of IDH. In this, it was found that for every 4 kg/m2 increase in BMI increased the odds of IDH by 52%.
What Do These Findings Mean?
These findings support a causal link between increased BMI and IDH risk, although it may be that BMI affects IDH through intermediate factors such as hypertension, dyslipidemia, and diabetes. The findings also show that observational studies into the impact of elevated BMI on IHD risk were consistent with this, but also that the inclusion of genetic data increases the value of observational studies by making it possible to avoid issues such as confounding and reverse causation. Finally, these findings and those of recent, observational studies have important implications for public-health policy because they show that the association between BMI (which is modifiable by lifestyle changes) and IHD is continuous. That is, any increase in BMI increases the risk of IHD; there is no threshold below which a BMI increase has no effect on IDH risk. Thus, public-health policies that aim to reduce BMI by even moderate levels could substantially reduce the occurrence of IDH in populations.
Additional Information
Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1001212.
Observational examination of the prospective association between body mass index (BMI) and ischemic heart disease (IHD) has been undertaken in a range of populations [1][3] and has consistently shown a positive relationship between BMI and the risk of IHD and other vascular endpoints [4][9]; however, causality in this relationship has not been convincingly demonstrated.
Observational association of BMI with IHD is impaired by confounding [10], reverse causation [11],[12], and bias [12], making it difficult to infer causality. An alternative approach free of confounding, reverse causation, and bias is that of Mendelian randomisation [13]. Mendelian randomisation uses analyses analogous to those in a randomized trial, but where randomization to risk-factor–related genotypes takes place at conception. In the case of IHD, genotypes with the largest known effects on BMI [14], (FTO[rs9939609], MC4R[rs17782313], and TMEM18[rs6548238]), can be used [15] for the reassessment of BMI as a causal risk factor for disease risk.
Genetic variation at the fat mass and obesity related locus (FTO) is thought to have a role in the hypothalamic regulation of appetite and food intake or metabolic rate [16][21] and has already been used to interrogate the relationship between BMI and acute coronary syndrome, lipid profile in myocardial infarction patients, and mortality [22][24]. Furthermore, variation at the same locus has been associated with atherogenic lipid profile and myocardial infarction risk in a manner that suggests an important role for BMI, but that remains unclear [23]. Although, the mechanism of the association between common variation at MC4R and TMEM18 loci and adiposity is currently unclear, these genotypes may still be used to provide evidence of the role of BMI in IHD risk. Indeed, the use of multiple independent genetic instruments for a single risk factor in the undertaking of Mendelian randomisation experiments is recognised as a favourable approach [25].
The aim here was to test the hypothesis that known positive observational relationships between BMI and IHD are causal in two well-sized general population studies (the Copenhagen General Population Study [CGPS]; n = 54,613 [3,780 patients with IHD] and the Copenhagen City Heart Study [CCHS]; n = 10,474 [2,006 patients with IHD]) and in a large collection of IHD patients with matched controls (the Copenhagen Ischaemic Heart Disease Study [CIHDS]; n = 10,540 [5,270 patients with IHD]). Analyses were designed to use an allele score as an unconfounded marker of BMI, and to take advantage of the use of loci FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238) as multiple independent instruments for BMI in efforts to address the potential complicating issue of biological confounding (or pleiotropy). In these analyses, allele score is used as an instrument for assessing the causal relationship between BMI increase and increased IHD risk, which is then directly compared with the observational BMI-IHD association. To evaluate “raw” data, we also show the different parts of the analyses separately, that is, first the BMI-IHD association, second the allele score-BMI association, and third the allele score-IHD association.
Participants
All participants were white and of Danish descent; this information is available through the national Danish Central Person Registry. No participants appeared in more than one of the three studies. The studies were approved by Danish ethical committees and Herlev Hospital.
Copenhagen General Population Study
This general population study was initiated in 2003 with ongoing enrolment [26][29]. IHD endpoints have been collected from 1976 to May 2009. Individuals were selected on the basis of the national Danish Civil Registration System to reflect the adult Danish population aged 20–100 y. Data were obtained from a questionnaire, a physical examination, blood samples, and from DNA. At the time of genotyping 59,883 participants had been included; of these, 5,270 were used as controls in the CIHDS (see below), leaving 54,613 for analyses in the CGPS.
Copenhagen City Heart Study
This prospective general population study was initiated in 1976–1978 with follow-up examinations in 1981–1983, 1991–1994, and 2001–2003 [26]. Participants were recruited and examined exactly as in the CGPS. DNA was available on 10,474 participants attending the 1991–1994 and/or 2001–2003 examinations.
Copenhagen Ischemic Heart Disease Study
This case-control study comprises 5,270 patients from the greater Copenhagen area referred for coronary angiography to Copenhagen University Hospital during the period 1991–2009 and 5,270 unmatched controls without IHD randomly sampled from the CGPS. Beside a diagnosis of IHD as described below, these patients also had stenosis/atherosclerosis on coronary angiography and/or a positive exercise electrocardiography test.
Ischemic Heart Disease
In all three studies, information on diagnosis of IHD (World Health Organization International Classification of Diseases: ICD8 410–414; ICD10 I20–I25) was collected and verified from existing data from 1976 until May 2009 by reviewing all hospital admissions and diagnoses entered in the national Danish Patient Registry and all causes of death entered in the national Danish Causes of Death Registry. Even though some individuals entered into our studies after 1976, we have complete information on all participants on any hospitalisation or death from IHD from 1976 through 2009 through these registries. IHD was angina pectoris and/or myocardial infarction (ICD8 410; ICD10 I21–I22), based on characteristic chest pain, electrocardiographic changes, and/or elevated cardiac enzymes. Follow-up was 100% complete, that is, no individual was lost to follow-up in any of the studies.
Genotyping
Genotyping was conducted blind to phenotypic data. In the absence of genomewide data, the ABI PRISM 7900HT Sequence Detection system (Applied Biosystems Inc.) was used to genotype the BMI instrument loci: FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238) using TaqMan assays. Polymorphisms were selected as those with the largest known common effect sizes for association with BMI in European populations [14],[30]. Genotyping was verified by DNA sequencing, quality control genotyping re-runs were performed twice, and 99.96% of all available participants were genotyped. To act as an aggregate instrument for BMI, a simple score of 0–6 was constructed as the sum of the number of higher BMI-associated alleles across the three genotypes in each study. To check results from this simplest form of allelic score, weighted allele score was also calculated. This score was generated by taking the weighted sum of the number of BMI “elevating” alleles at each locus scaled per standard deviation increase in BMI. Independent gene variant/BMI effect estimates used were taken from the largest available GWAS for BMI to date [14].
Body Mass Index
BMI was calculated as measured weight (kg) divided by measured height squared (m2). To exclude influence of age and sex on our results, BMI was standardised into age- and sex-adjusted Z-scores within each study separately (Table S1). One Z-score corresponds to a BMI standard deviation of 4 kg/m2; thus, for easy interpretation of data, results in Z-scores were converted to BMI values in kg/m2. Participants with BMI greater than 50 kg/m2 were excluded from observational and instrumental variable estimates. BMI measurements were not available for the patients in the CIHDS, hence observational estimates of the IHD-BMI association are not given for this study.
Other Covariates
Smoking was categorized from self-reported data as ever versus never smoked, alcohol consumption was categorized as <14/21 or ≥14/21 units per week for women/men (1 unit = 12 g), education as schooling for <10, ≥10 to <13, and ≥13 y, and annual income as <100,000DKK, 100,000DKK–400,000DKK, 400,000DKK–600,000DKK, and >600,000DKK. Use of statins was recorded at examination. Measured systolic and diastolic blood pressure was recorded as described previously [31] and adjusted for antihypertensive medication by adding a constant value of 10 mmHg and systolic blood pressure and 5 mmHg for diastolic blood pressure [32]. Plasma levels of triglycerides, glucose, high-density lipoprotein cholesterol, and low-density lipoprotein cholesterol were measured using standard hospital assays. Time between IHD event and BMI measurement (“event time”) was assessed as a confounding factor in supplementary analyses.
Statistical Analysis
All analyses were performed in Stata version 11 (StataCorp). For genotypes a deviation from Hardy-Weinberg equilibrium was investigated using a Pearson chi-squared test. Observational estimates of odds ratios (ORs) of IHD per unit increase in standardised BMI, allele score, and potential confounders were estimated using logistic regression. Estimates of the association of standardised BMI with allele score and potential confounders were performed using linear regression. Associations of potential confounders by genotype and allele score were estimated using logistic regression, linear regression, and Pearson's chi-squared test.
Instrumental variable estimates of causal ORs were derived using the Wald-type estimator [33], which involves taking the ratio of the IHD-allele score log OR to the standardised BMI-allele score coefficient and then exponentiating to express as an OR. The standardised BMI-allele score coefficient comes from all individuals from the general population with BMI available, that is, all participants other than those in the CIHDS. Standard errors of Wald-type instrumental variable log ORs were derived using the delta method [34]. Instrumental variable estimates derived from this approach (which is the most efficient use of data) were verified by comparing them to analyses in the CGPS and CCHS using a logistic structural mean model fitted using the three genotypes as separate, weighted, instruments [35]; these models were fitted in the CGPS and CCHS as they require observations of genotypes, BMI, and IHD status for each individual. In these models the joint validity of the multiple instruments was investigated using Hansen's over-identification test [36]. If such a test is rejected it indicates one or more of the instruments may not be valid instruments. As a comparison to using a non-weighted allele score, analyses were also undertaken using a weighted allele score.
In a case-control study like the CIHDS, the allele score-IHD association is valid as in cohort studies like the CGPS and the CCHS because it is not affected by IHD status and represents lifecourse BMI-associated IHD risk. BMI may be affected by IHD, because some patients loose weight after an IHD diagnosis. Therefore, the allele score-BMI association entering into the instrumental variable analyses is best derived from people in the general population as done in the present study. The advantage of including the CIHDS together with the CGPS and the CCHS is that it adds considerable statistical power to the combined analyses, as done by ourselves in previous studies [27][29]. Importantly, however, when the CGPS and the CCHS are combined and CIHDS is excluded, the results are similar.
Meta-analysis pooled estimates were obtained using the random effects meta-analysis model implemented in the user-written Stata command “metan” [37]. In each meta-analysis, between-study heterogeneity was assessed using the I2 statistic and the Cochran heterogeneity Q test. Meta-analysis of ORs was performed on the log scale. Analysing the three studies as a single dataset was undertaken, but results are not reported as they were similar to those of the meta-analyses throughout.
The numbers of IHD patients were 3,780 of 54,613 participants in the CGPS (79% prevalent and 21% incident), 2,006 of 10,474 participants in the CCHS (22% prevalent and 78% incident), and 5,270 of 10,540 participants in the CIHDS (100% prevalent). Baseline characteristics of participants in the three studies are shown in Table 1. Genotypes scored for FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238) were in Hardy-Weinberg equilibrium in the three studies and the distribution of the allele score was approximately normal (Table 2).
Table 1
Table 1
Characteristics of the participants in the CGPS, CCHS, and CIHDS by BMI tertile.
Table 2
Table 2
Distribution of FTO, MC4R, and TMEM18 genotypes in the CGPS, CCHS, and CIHDS.
BMI and IHD Risk: Observational Estimate
In the general population, for one standard deviation (4 kg/m2) increase in BMI, the observational ORs for IHD were 1.23 (95% CI 1.19–1.28) for the CGPS and 1.31 (95% CI 1.23–1.39) for the CCHS (Figure 1). The meta-analysis pooled OR was 1.26 (95% CI 1.19–1.34), and there was some evidence of between-study heterogeneity (I2 = 63; p = 0.10), but this was incorporated in the random effects model (Figure 2). These observational associations were unaffected on adjustment for the potential confounders of sex, smoking status, drinking status, years of education, and income level and attenuated on additional adjustment for event time (Table 3). If ever smokers were separated into current and former smokers and drinkers into moderate and heavy drinkers, observational associations were unaffected. Where available (CPGS), analysis including adjustment for statin use showed evidence for an attenuation in estimates of the relationship between BMI and IHD risk (adjusted OR for IHD per 4 kg/m2 increase in BMI 1.12 [1.08–1.16]).
Figure 1
Figure 1
ORs and 95% CIs for IHD by BMI categories in the CGPS and the CCHS.
Figure 2
Figure 2
Meta-analysis forest plots of observational and instrumental variable causal estimates using allele score of the relationship between IHD and BMI.
Table 3
Table 3
Observational ORs for IHD per 4 kg/m2 increase in BMI.
Genotype and Allele Score Associations with BMI
Pooled across the three studies, each additional adiposity-related allele from the allele score was associated with a 0.28 kg/m2 (95% CI 0.22–0.34) increase in BMI (Figure 3), corresponding to a 1.68 kg/m2 (95% CI 1.31–2.06) BMI increase for a comparison of the minimum (0) to maximum (6) carriage of adiposity increasing alleles. These results were similar across cases and controls in the three separate studies.
Figure 3
Figure 3
Meta-analysis forest plots of the relationships between FTO rs9939609, MC4R rs17782313, and TMEM18 rs6548238 allele score and BMI.
Associations with Confounding Factors
In an analysis such as this one, it is important to examine whether potential confounding factors could be part of the explanation behind an observational BMI-IHD association, or behind a causal BMI-IHD association estimated using allele score in instrumental variable analysis. For a potential confounder to be part of the explanation behind an association, the factor in question needs to associate both with the exposure and the outcome.
We therefore examined both the association between potentially confounding factors (sex, age, smoking status, drinking status, education, income, and event time) and our primary exposure (BMI), our primary outcome (IHD), and our instrumental variables (genotypes and allele score). For both our primary outcome and exposure, there was strong evidence for association between both BMI, IHD, and drinking, education, income, and event time (Tables S2 and S3). In contrast, genotypes at FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238) and the allele score demonstrated no reliable evidence of association with potential confounders (Table 4).
Table 4
Table 4
Associations of potential confounders with genotypes in the three studies.
In aggregate, these data suggest that several factors likely could confound the observational BMI-IHD association. In contrast, it is unlikely that these same factors should confound the instrumental variable analyses assessing the causal BMI-IHD association, as these factors were not associated with genotype or allele score.
Allele Score and IHD Risk
In the CGPS the OR for IHD per risk allele, that is per 0.28 kg/m2 increase, was 1.02 (95% CI 0.98–1.05), in the CCHS 1.05 (95% CI 1.01–1.11), and in the CIHDS 1.03 (95% CI 1.00–1.07). The meta-analysed OR per risk allele was 1.03 (95% CI 1.01–1.05), with no evidence of between-study heterogeneity (I2 = 0, p = 0.43) (Figure 4). This finding corresponds to an OR for IHD of 1.20 (95% CI 1.05–1.36) for a maximum allele score comparison of 0–6 adiposity increasing alleles corresponding to a 1.68 kg/m2 increase. The increasing trend in BMI with allele score, distribution of allele score, and the relative differences in BMI by IHD status can be seen simultaneously in Figure 5. The relatively high value of BMI at 0 allele score is likely due to the play of chance, given that there are relatively few individuals with this allele score. The absolute BMI values are higher in the CGPS because participants were recruited approximately 15 y after the CCHS, during which time Danish people on average become heavier for height. Where available (CPGS), analysis including adjustment for statin use showed no strong evidence of attenuation, (adjusted OR for IHD per risk allele 1.01 [95% CI 0.98–1.05]).
Figure 4
Figure 4
Meta-analysis forest plots of the relationships between FTO rs9939609, MC4R rs17782313, and TMEM18 rs6548238 allele score and IHD.
Figure 5
Figure 5
Mean BMI and 95% CIs by allele score and IHD status and distribution of allele score in the CGPS and CCHS.
BMI and IHD Risk: Instrumental Variable Causal Estimate Using Allele Score
An instrumental variable analysis using allele score examines the causal effect of a lifecourse change in BMI on the risk of IHD. The instrumental variable estimate of the causal relationship between a 4 kg/m2 increase in BMI and IHD showed an OR of 1.31 (95% CI 0.76–2.26) in the CGPS, 2.11 (95% CI 1.05–4.24) in the CCHS, and 1.46 (95% CI 0.96–2.24) in the CIHDS. The meta-analysed causal OR was 1.52 (95% CI 1.12–2.05), with no evidence of between-study heterogeneity (I2 = 0, p = 0.56) (Figure 2). These results were similar using a weighted allele score (Figure S1), rather than the non-weighted model used above. There was no strong evidence of a difference between instrumental variable and observational ORs (p = 0.25 for difference). Where available (CPGS), analysis including adjustment for statin use showed that, in contrast to observational analysis, there was no evidence for an attenuation of the estimate (adjusted OR for IHD per 4 kg/m2 increase in BMI 1.28 [0.69–2.37]).
Instrumental variable estimates based on the use of individual SNPs gave broadly similar results, but with reduced statistical power (Figure S2). Also instrumental variable estimates from a logistic structural mean model gave similar estimated causal ORs (Figure S3). In these models there was no evidence against the joint validity of using these SNPs as multiple instruments (CGPS, p = 0.97; CCHS, p = 0.91).
As the relationships between BMI and health outcomes can be attenuated with age, we performed subgroup analyses in individuals above and below age 60 y. There was evidence of smaller observational ORs for IHD in those above versus below 60 y (p = 0.02 for difference). This finding was reflected in evidence for interaction between age group and BMI when included in the logistic regression of IHD on BMI and age group (p = 0.01 [CGPS], p<0.001 [CCHS] for interaction). Instrumental variable estimates gave no evidence of difference by age (Figure 6).
Figure 6
Figure 6
Meta-analysis forest plots of observational and instrumental variable estimates of the relationship between IHD and BMI stratified by age group.
BMI and IHD Risk: Observational Versus Instrumental Variable Causal Estimate Using Allele Score
In summary, in observational analyses for every 4 kg/m2 increase in BMI, the OR for IHD was 1.26 (95% CI 1.19–1.34), corresponding to a 26% increased IHD risk (Figure 2). For comparison, in instrumental variable analysis using the allele score the causal IHD ORs for a 4 kg/m2 increase in BMI was 1.52 (95% CI 1.12–2.05), corresponding to a 52% increased IHD risk.
In two large studies of the general population, observational estimates suggested a 26% increase in risk of IHD for every 4 kg/m2 increase in BMI. A unit change of an allelic score combining genotypes from three established genetic associates of BMI was associated with a 0.28 kg/m2 change in BMI, a change neither correlated with classic confounding features nor affected by reverse causation. Using this as an instrument for lifecourse BMI difference in 75,627 Danish individuals from the same general population studies and a further case-control series, instrumental variable analysis was employed to re-estimate the causal relationship between BMI and IHD risk. In doing this, estimates suggest that the same increase in BMI is causally related to an increased risk of IHD consistent with observational estimates, if not greater. Whilst features such as statin use appear to impact observational estimates, those based on genetic instruments for BMI appear consistent across sub-analyses. Importantly, we not only qualify the likely causal role of BMI (rather than just an observational associate), we are able to quantify this effect. It should be mentioned that the causal relationship between increased BMI and increased risk of IHD may be realised through intermediate factors like hypertension, dyslipidemia, and type 2 diabetes.
The design of this study and the use of genetic variation as a proxy marker for elevated BMI lead to additional discussion points. Firstly, the use of effective instruments for BMI has allowed for forms of study bias to be effectively accounted for. For example, observational relationships between BMI and IHD risk tend to increase with time from BMI measurement largely as a result of tracking and the natural progression of BMI with age [2]. In contrast to this, where IHD events are recorded before BMI, reverse causation can reduce the strength of relationships between BMI and IHD, though illness-induced weight loss adds further complication to the assessment of relationships between body weight and mortality [11]. Genetic instruments are not related to this time difference and the use of these as proxy markers for BMI removes the limitations that can be brought about by the existence of these biases. Use of genetic variation as a proxy measure for lifecourse BMI change also helps to avoid regression dilution bias, or the artefactual underestimation of epidemiological associations owing to a failure of baseline measurements to reliably collect the “usual” levels of an exposure [38]. Existing studies have made efforts to account for this and other “confounding at baseline” issues [2],[3],[6]; however, these effects are difficult to quantify. Owing to the consistent nature of genotype effects through adult life [14],[39], this type of loss in measurement precision is minimised. Indeed, as a result of Mendelian randomisation-derived estimates modelling lifecourse differences in BMI, we would anticipate effect estimates that can be larger than those from more conventional observational approaches, as seen in this case.
Despite these benefits, there are a number of potential limitations in Mendelian randomisation studies like the present one [13],[40],[41]. However, the complicating effect of population stratification in Mendelian randomisation studies is likely to have been avoided through the use of an ethnically homogenous white population in the present study. Also, pleiotropy and the potentially confounding effects of linkage disequilibrium are likely avoided owing to the use of multiple genetic polymorphisms, each associated with increased BMI and each acting on BMI independently and via different pathways. Nevertheless, canalisation cannot be completely excluded as a limitation of the present study; however, as canalisation theoretically acts to buffer the effect of genetic deviations, it can obscure associations between genes and IHD, but it is unlikely to explain away the results in this study. Also, although the analyses provide insight into the likely causal effect of lifecourse elevations in BMI, we are not well placed to comment on the impact of acute changes. Lastly, there are known limitations to the use of BMI as a marker for adiposity or other possible anthropometrically related aetiological contributions to IHD risk, and using alternative measures may prove more informative [42]. Nevertheless, these potential limitations to BMI cannot explain the results in this study.
In the context of available evidence concerning the causal role of BMI as an intermediate risk factor for IHD, we can speculate that the explanation for the causal association is straightforward: increased BMI contributes causally to well-known cardiovascular risk factors including hypertension, dyslipidemia, and type 2 diabetes, factors that may then go on to cause the observed increased risk of IHD. Similar evidence supporting the role of elevated BMI in the generation of a common risk profile is emerging [23],[24],[31],[43],[44] and in combination with the agreement of these findings with those of observational studies, an important role for BMI is evident from work using other techniques to avoid the problems of confounding and reverse causation [22],[45],[46].
In conclusion, for every 4 kg/m2 increase in BMI, observational estimates suggested a 26% increase in IHD risk with instrumental variable analysis suggesting a causal 52% increase in IHD risk. These data add novel evidence to support a causal link between increased BMI and increased IHD risk, while the mechanism of this effect is likely to be operating through intermediate factors. In the context of recent, high impact, observational findings, this work has important policy implications for public health given the continuous nature of the BMI-IHD association, the modifiable nature of BMI, and the likely benefits of reducing BMI even by moderate levels. Finally, this analysis demonstrates the value of observational studies and their ability to provide essentially unbiased results because of inclusion of genetic data avoiding confounding, reverse causation, and bias.
Figure S1
Meta-analysis forest plots of observational and instrumental variable estimates using a weighted allele score of the relationship between IHD and standardised BMI.
(TIF)
Figure S2
Meta-analysis forest plots of instrumental variable causal estimates of the relationship between IHD and BMI stratified by genotype. FTO rs9939609, MC4R rs17782313, and TMEM18 rs6548238.
(TIF)
Figure S3
Meta-analysis of logistic structural mean model causal OR estimates of IHD risk per 4 kg/m2 increase in BMI in the CGPS and CCHS. Logistic structural mean models fitted using FTO, MC4R, and TMEM18 genotypes as multiple instruments, with each genotype under an additive model. The first stage association model was fitted with an intercept and the main effects of standardised BMI and each of the three genotypes.
(TIF)
Table S1
Means and standard deviations of BMI by 5-y age band and sex used to generate standardised BMI in the CGPS and CCHS.
(DOCX)
Table S2
Associations of potential confounders with standardised BMI in the three studies.
(DOCX)
Table S3
Associations of potential confounders with IHD in the three studies.
(DOCX)
Acknowledgments
We thank Dorthe Uldall Andersen and Anne Bank for technical assistance.
Abbreviations
BMIbody mass index
CCHSCopenhagen City Heart Study
CGPSCopenhagen General Population Study
CIHDSCopenhagen Ischaemic Heart Disease Study
IHDischemic heart disease

Footnotes
George Davey Smith is on the PLoS Medicine Editorial Board. The authors have declared that no other competing interests exist.
Supported by the Copenhagen County Foundation and by Herlev Hospital, Copenhagen University Hospital. These sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; or in preparation, review, or approval of the manuscript. TMP, NJT, and GDS work within and are funded by the MRC Centre for Causal Analyses in Translational Epidemiology, which is capacity funded by grant G0600705. BGN, TMP, MB, and NJT 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.
1. Manson JE, Colditz GA, Stampfer MJ, Willett WC, Rosner B, et al. A prospective study of obesity and risk of coronary heart disease in women. N Engl J Med. 1990;322:882–889. [PubMed]
2. Prospective Studies C, Whitlock G, Lewington S, Sherliker P, Clarke R, et al. Body-mass index and cause-specific mortality in 900 000 adults: collaborative analyses of 57 prospective studies. Lancet. 2009;373:1083–1096. [PMC free article] [PubMed]
3. Ni Mhurchu C, Rodgers A, Pan WH, Gu DF, Woodward M, et al. Body mass index and cardiovascular disease in the Asia-Pacific Region: an overview of 33 cohorts involving 310 000 participants. Int J Epidemiol. 2004;33:751–758. [PubMed]
4. Kurth T, Gaziano JM, Rexrode KM, Kase CS, Cook NR, et al. Prospective study of body mass index and risk of stroke in apparently healthy women. Circulation. 2005;111:1992–1998. [PubMed]
5. Zheng W, McLerran DF, Rolland B, Zhang X, Inoue M, et al. Association between body-mass index and risk of death in more than 1 million Asians. N Engl J Med. 2011;364:719–729. [PubMed]
6. Song YM, Sung J, Davey Smith G, Ebrahim S, Song Y-M, et al. Body mass index and ischemic and hemorrhagic stroke: a prospective study in Korean men. Stroke. 2004;35:831–836. [PubMed]
7. Bogers RP, Bemelmans WJ, Hoogenveen RT, Boshuizen HC, Woodward M, et al. Association of overweight with increased risk of coronary heart disease partly independent of blood pressure and cholesterol levels: a meta-analysis of 21 cohort studies including more than 300 000 persons. Arch Intern Med. 2007;167:1720–1728. [PubMed]
8. Rexrode KM, Hennekens CH, Willett WC, Colditz GA, Stampfer MJ, et al. A prospective study of body mass index, weight change, and risk of stroke in women. JAMA. 1997;277:1539–1545. [PubMed]
9. Zhou M, Offer A, Yang G, Smith M, Hui G, et al. Body mass index, blood pressure, and mortality from stroke: a nationally representative prospective study of 212,000 Chinese men. Stroke. 2008;39:753–759. [PubMed]
10. Shahar E. The association of body mass index with health outcomes: causal, inconsistent, or confounded? Am J Epidemiol. 2009;170:957–958. [PubMed]
11. Flegal KM, Graubard BI, Williamson DF, Cooper RS. Reverse causation and illness-related weight loss in observational studies of body weight and mortality. Am J Epidemiol. 2011;173:1–9. [PubMed]
12. Greenberg JA. Correcting biases in estimates of mortality attributable to obesity. Obesity. 2006;14:2071–2079. [PubMed]
13. Davey Smith G, Ebrahim S. ‘Mendelian Randomisation’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22. [PubMed]
14. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, et al. Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. 2010;42:937–948. [PMC free article] [PubMed]
15. Greenland S. An introduction to instrumental variables for epidemiologists. Int J Epidemiol. 2000;29:722–729. [PubMed]
16. Church C, Moir L, McMurray F, Christophe G, Banks GT, et al. Overexpression of Fto leads to increased food intake and results in obesity. Nat Genet. 2010;42:1086–1093. [PMC free article] [PubMed]
17. Church C, Lee S, Bagg EA, McTaggart JS, Deacon R, et al. A mouse model for the metabolic effects of the human fat mass and obesity associated FTO gene. PLoS Genetics. 2009;5:e1000599. doi: 10.1371/journal.pgen.1000599. [PMC free article] [PubMed]
18. Tung YC, Ayuso E, Shan X, Bosch F, O'Rahilly S, et al. Hypothalamic-specific manipulation of Fto, the ortholog of the human obesity gene FTO, affects food intake in rats. PLoS ONE. 2010;5:e8771. doi: 10.1371/journal.pone.0008771. [PMC free article] [PubMed]
19. Gerken T, Girard CA, Tung YL, Webby CJ, Saudek V, et al. The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Science. 2007;318:1469–1472. [PMC free article] [PubMed]
20. Timpson NJ, Emmett P, Frayling TM, Rogers I, Hattersley AT, et al. The FTO/obesity associated locus and dietary intake in children. Am J Clin Nutr. 2008;88:971–978. [PubMed]
21. Fischer J, Koch L, Emmerling C, Vierkotten J, Peters T, et al. Inactivation of the Fto gene protects from obesity. Nature. 2009;458:894–898. [PubMed]
22. Hubacek JA, Stanek V, Gebauerova M, Pilipcincova A, Dlouha D, et al. A FTO variant and risk of acute coronary syndrome. Clinica Chimica Acta. 2010;411:1069–1072. [PubMed]
23. Doney ASF, Dannfald J, Kimber CH, Donnelly LA, Pearson E, et al. The FTO gene is associated with an atherogenic lipid profile and myocardial infarction in patients with type 2 diabetes: a Genetics of Diabetes Audit and Research Study in Tayside Scotland (Go-DARTS) study. Circ Cardiovasc Genet. 2009;2:255–259. [PMC free article] [PubMed]
24. Freathy RM, Timpson NJ, Lawlor DA, Pouta A, Ben-Shlomo Y, et al. Common variation in the FTO gene has a metabolic impact consistent with its effect on BMI. Diabetes. 2008;57:1419–1426. [PMC free article] [PubMed]
25. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2011 In press. [PubMed]
26. Nordestgaard BG, Benn M, Schnohr P, Tybjaerg-Hansen A. Nonfasting triglycerides and risk of myocardial infarction, ischemic heart disease, and death in men and women. JAMA. 2007;298:299–308. [PubMed]
27. Frikke-Schmidt R, Nordestgaard BG, Stene MCA, Sethi AA, Remaley AT, et al. Association of loss-of-function mutations in the ABCA1 gene with high-density lipoprotein cholesterol levels and risk of ischemic heart disease. JAMA. 2008;299:2524–2532. [PubMed]
28. Zacho J, Tybjaerg-Hansen A, Jensen JS, Grande P, Sillensen H, et al. Genetically elevated C-reactive protein and ischaemic vascular disease. N Engl J Med. 2008;359:1897–1908. [PubMed]
29. Kamstrup PR, Tybjaerg-Hansen A, Steffensen R, Nordestgaard BG. Genetically elevated Lipoprotein(a) and increased risk of myocardial infarction. JAMA. 2009;301:2331–2339. [PubMed]
30. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, et al. Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008;40:161–169. [PubMed]
31. Timpson NJ, Harbord R, Davey Smith G, Zacho J, Tybjaerg-Hansen A, et al. Does greater adiposity increase blood pressure and hypertension risk? Mendelian randomisation using FTO/MC4R genotype. Hypertension. 2009;54:84–90. [PubMed]
32. Tobin MD, Sheehan NA, Scurrah KJ, Burton PR. Adjusting for treatment effects in studies of quantitative traits: antihypertensive therapy and systolic blood pressure. Stat Med. 2005;24:2911–2935. [PubMed]
33. Didelez V, Meng S. Assumptions of IV methods for observational epidemiology. Stat Sci. 2010;25:22–40.
34. Thomas DC, Lawlor DA, Thompson JR, Thomas DC, Lawlor DA, et al. re: Estimation of bias in nongenetic observational studies using “Mendelian triangulation” by Bautista et al. Ann Epidemiol. 2007;17:511–513. [PubMed]
35. Vansteelandt S, Goetghebeur E. Causal inference with generalised structural mean models. J R Stat Soc Series B Stat Methodol. 2003;65:817–835.
36. Hansen LP. Large sample properties of generalized method of moments estimators. Econometrica. 1982;50:1029–1054.
37. Harris R. metan: fixed- and random-effects meta-analysis. Stata J. 2008;8:3–28.
38. Clarke R, Shipley M, Lewington S, Youngman L, Collins R, et al. Underestimation of risk associations due to regression dilution in long-term follow-up of prospective studies. Am J Epidemiol. 1999;150:341–353. [PubMed]
39. Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007;316:889–894. [PMC free article] [PubMed]
40. Lawlor DA, Harbord RM, Sterne JAC, Timpson NJ, Davey Smith G. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27:1133–1163. [PubMed]
41. Waddington CH. Canalisation of development and the inheritance of acquireed characters. Nature. 1942;150:563–565.
42. Yusuf S, Hawken S, Ounpuu S, Bautista L, Franzosi MG, et al. Obesity and the risk of myocardial infarction in 27,000 participants from 52 countries: a case-control study. Lancet. 2005;366:1640–1649. [PubMed]
43. Timpson NJ, Nordestgaard BG, Harbord RM, Zacho J, Frayling TM, et al. C-reactive protein levels and body mass index: elucidating direction of causation through reciprocal Mendelian randomization. Int J Obes (Lond) 2010;35:300–308. [PubMed]
44. Timpson NJ, Lawlor DA, Harbord R, Gaunt TR, Day IN, et al. C-reactive protein and metabolic syndrome: Mendelian randomisation suggests associations are non-causal. Lancet. 2005;366:1954–1959. [PubMed]
45. Bjorge T, Engeland A, Tverdal A, Davey Smith G. Body mass index in adolescence in relation to cause-specific mortality: a follow-up of 230,000 Norwegian adolescents. Am J Epidemiol. 2008;168:30–37. [PubMed]
46. Davey Smith G, Sterne JA, Fraser A, Tynelius P, Lawlor DA, et al. The association between BMI and mortality using offspring BMI as an indicator of own BMI: large intergenerational mortality study. BMJ. 2009;339:b5043. [PMC free article] [PubMed]
Articles from PLoS Medicine are provided here courtesy of
Public Library of Science