MR-pheWAS with stratification and interaction: Searching for the causal effects of smoking heaviness identified an effect on facial aging
Authors: Millard Louise A C , Munafò Marcus R , Tilling Kate , Wootton Robyn E , Smith George Davey 
Published: October 2018

Mendelian randomization (MR) is an established approach for estimating the causal effect of an environmental exposure on a downstream outcome. The gene x environment (GxE) study design can be used within an MR framework to determine whether MR estimates may be biased if the genetic instrument affects the outcome through pathways other than via the exposure of interest (known as horizontal pleiotropy). MR phenome-wide association studies (MR-pheWAS) search for the effects of an exposure, and a recently published tool (PHESANT) means that it is now possible to do this comprehensively, across thousands of traits in UK Biobank. In this study, we introduce the GxE MR-pheWAS approach, and search for the causal effects of smoking heaviness – stratifying on smoking status (ever versus never) – as an exemplar. If a genetic variant is associated with smoking heaviness (but not smoking initiation), and this variant affects an outcome (at least partially) via tobacco intake, we would expect the effect of the variant on the outcome to differ in ever versus never smokers. If this effect is entirely mediated by tobacco intake, we would expect to see an effect in ever smokers but not never smokers. We used PHESANT to search for the causal effects of smoking heaviness, instrumented by genetic variant rs16969968, among never and ever smokers respectively, in UK Biobank. We ranked results by: 1) strength of effect of rs16969968 among ever smokers, and 2) strength of interaction between ever and never smokers. We replicated previously established causal effects of smoking heaviness, including a detrimental effect on lung function and pulse rate. Novel results included a detrimental effect of heavier smoking on facial aging. We have demonstrated how GxE MR-pheWAS can be used to identify causal effects of an exposure, while simultaneously assessing the extent that results may be biased by horizontal pleiotropy.

Author summary

Mendelian randomization uses genetic variants associated with an exposure to investigate causality. For instance, a genetic variant that relates to how heavily a person smokes has been used to test whether smoking causally affects health outcomes. Mendelian randomization is biased if the genetic variant also affects the outcome via other pathways. We exploit additional information – that the effect of heavy smoking only occurs in people who actually smoke – to overcome this problem. By testing associations in ever and never smokers separately we can assess whether the genetic variant affects an outcome via smoking or another pathway. If the effect is entirely via smoking heaviness, we would expect to see an effect in ever but not never smokers, and this would suggest that smoking causally influences the outcome. Previous Mendelian randomization studies of smoking heaviness focused on specific outcomes – here we searched for the causal effects of smoking heaviness across over 18,000 traits. We identified previously established effects (e.g. a detrimental effect on lung function) and novel results including a detrimental effect of heavier smoking on facial aging. Our approach can be used to search for the causal effects of other exposures, where the exposure only occurs in known subsets of the population.