Problems in interpreting and using GWAS of conditional phenotypes illustrated by 'alcohol GWAS'.

Dear Editor,

We read with interest findings of a genome-wide association study (GWAS) in the UK Biobank reported by Clarke et al1 in Molecular Psychiatry. The authors aimed to conduct a GWAS of alcohol consumption. Prior to conducting the GWAS, Clarke and colleagues regressed alcohol intake on body weight and age, and used the residuals from this regression analysis for the GWAS analysis. The authors identified 14 loci harbouring variants associated with the residuals at levels indicative of GWAS significance and took these forward to characterise the relationship of these SNPs with various phenotypes, such as adiposity.

By adjusting alcohol for body weight prior to the GWAS, Clarke and colleagues in effect generate a phenotype which is automatically associated with body weight and its correlates.2 When the SNPs from the GWAS are combined (as the authors do), this would have the effect of generating a genetic instrument that has a positive association with alcohol and an inverse association with body weight. This can lead to spurious findings when the SNPs are then related to other traits. For example, in figure 3, the authors note a negative genetic correlation of alcohol consumption (indexed by SNPs associating with alcohol adjusted for body weight) with various markers of adiposity including waist circumference, body mass index, overweight and obesity. In Table 2, these relationships are more formally characterised: a polygenic risk score associating with 0.08 SD higher alcohol per week was associated with a 0.024 SD lower body mass index (BMI). A naïve interpretation would be that more alcohol leads to less adiposity, however this would be incorrect. Large-scale Mendelian randomization meta-analysis using rs1229984 in ADH1B has identified the genetic variant associated with more alcohol intake to be strongly associated with a higher (not lower) body mass index.3 A recent instrumental variable study that used flushing as a proxy for ALDH2 genotype also suggested strong associations of more alcohol with higher BMI 4. Thus, a large body of evidence using reliable Mendelian randomization approaches indicates that consuming more alcohol, all things being equal, leads to higher BMI.

So why the discordance between these Mendelian randomization studies and the current report by Clarke and colleagues? As described above, by regressing alcohol for body weight prior to the GWAS, this means that SNPs that were identified by the GWAS had a pattern of association indicative of higher alcohol but lower body weight. In effect, the associations that Clarke and colleagues report with adiposity have been induced by the very nature of their statistical model.

An analogy is conducting a GWAS for waist hip ratio (WHR) adjusted for BMI (WHRadjBMI) – on combining the 49 WHRadjBMI SNPs together, the SNPs collectively associate with higher WHR but also associate with lower BMI5, 6. This means that while the aim of deriving the WHRadjBMI phenotype was to proxy central adiposity independently of general adiposity, the net consequence is to derive a genetic instrument that has opposing effects of contrasting adiposity elements, which can introduce confounding7 and make challenging the interpretation of findings when taking these variants forward into the Mendelian randomization setting.

The findings by Clarke et al1 highlight the complexities that can arise in Mendelian randomization analysis8. We caution against interpreting these data as suggesting that alcohol consumption is correlated with “many positive health and behavioural traits”. The burden of evidence indicates that higher alcohol consumption leads to higher levels of many deleterious cardiovascular risk factors including adiposity, blood pressure, inflammation markers, coronary calcification and vascular disease.