Hostname: page-component-7bb8b95d7b-l4ctd Total loading time: 0 Render date: 2024-10-03T07:39:45.992Z Has data issue: false hasContentIssue false

Evidence for causal effects of lifetime smoking on risk for depression and schizophrenia: a Mendelian randomisation study

Published online by Cambridge University Press:  06 November 2019

Robyn E. Wootton*
Affiliation:
School of Experimental Psychology, University of Bristol, BristolBS8 1TU, UK MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK NIHR Biomedical Research Centre at the University Hospitals Bristol NHS Foundation Trust and the University of Bristol, BristolBS8 2BN, UK
Rebecca C. Richmond
Affiliation:
MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK
Bobby G. Stuijfzand
Affiliation:
Jean Golding Institute, Royal Fort House, University of Bristol, BristolBS8 1UH, UK
Rebecca B. Lawn
Affiliation:
School of Experimental Psychology, University of Bristol, BristolBS8 1TU, UK MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK
Hannah M. Sallis
Affiliation:
School of Experimental Psychology, University of Bristol, BristolBS8 1TU, UK MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK
Gemma M. J. Taylor
Affiliation:
Department of Psychology, Addiction and Mental Health Group (AIM), University of Bath, BathBA2 7AY, UK
Gibran Hemani
Affiliation:
MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK
Hannah J. Jones
Affiliation:
MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK
Stanley Zammit
Affiliation:
Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK Division of Psychological Medicine and Clinical Neurosciences, MRC Centre for Neuropsychiatric Genetics and Genomics, Cardiff University, CardiffCF24 4HQ, UK
George Davey Smith
Affiliation:
MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK Department of Population Health Sciences, Bristol Medical School, University of Bristol, BristolBS8 2PR, UK
Marcus R. Munafò
Affiliation:
School of Experimental Psychology, University of Bristol, BristolBS8 1TU, UK MRC Integrative Epidemiology Unit, University of Bristol, BristolBS8 2PR, UK UK Centre for Tobacco and Alcohol Studies, University of Bristol, BristolBS8 1TU, UK
*
Author for correspondence: Robyn E. Wootton, E-mail: robyn.wootton@bristol.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Background

Smoking prevalence is higher amongst individuals with schizophrenia and depression compared with the general population. Mendelian randomisation (MR) can examine whether this association is causal using genetic variants identified in genome-wide association studies (GWAS).

Methods

We conducted two-sample MR to explore the bi-directional effects of smoking on schizophrenia and depression. For smoking behaviour, we used (1) smoking initiation GWAS from the GSCAN consortium and (2) we conducted our own GWAS of lifetime smoking behaviour (which captures smoking duration, heaviness and cessation) in a sample of 462690 individuals from the UK Biobank. We validated this instrument using positive control outcomes (e.g. lung cancer). For schizophrenia and depression we used GWAS from the PGC consortium.

Results

There was strong evidence to suggest smoking is a risk factor for both schizophrenia (odds ratio (OR) 2.27, 95% confidence interval (CI) 1.67–3.08, p < 0.001) and depression (OR 1.99, 95% CI 1.71–2.32, p < 0.001). Results were consistent across both lifetime smoking and smoking initiation. We found some evidence that genetic liability to depression increases smoking (β = 0.091, 95% CI 0.027–0.155, p = 0.005) but evidence was mixed for schizophrenia (β = 0.022, 95% CI 0.005–0.038, p = 0.009) with very weak evidence for an effect on smoking initiation.

Conclusions

These findings suggest that the association between smoking, schizophrenia and depression is due, at least in part, to a causal effect of smoking, providing further evidence for the detrimental consequences of smoking on mental health.

Type
Original Articles
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2019

Introduction

Smoking is a major risk factor for lung cancer, cardiovascular and respiratory diseases making it the leading cause of preventable death worldwide (World Health Organization, 2011). In developed nations, smoking is more common amongst individuals with mental health conditions (Coulthard, Farrell, Singleton, & Meltzer, Reference Coulthard, Farrell, Singleton and Meltzer2002; Lasser et al., Reference Lasser, Boyd, Woolhandler, Himmelstein, McCormick and Bor2000; Lawrence, Mitrou, & Zubrick, Reference Lawrence, Mitrou and Zubrick2009; McClave, McKnight-Eily, Davis, & Dube, Reference McClave, McKnight-Eily, Davis and Dube2010), in particular schizophrenia (Royal College of Physicians, 2013) and depression (Byers et al., Reference Byers, Vittinghoff, Lui, Hoang, Blazer, Covinsky and Fredman2012; Leung, Gartner, Dobson, Lucke, & Hall, Reference Leung, Gartner, Dobson, Lucke and Hall2011; Tjora et al., Reference Tjora, Hetland, Aarø, Wold, Wiium and Øverland2014). In the UK, estimates suggest that up to 45% of individuals with schizophrenia, and 31% of individuals with depression smoke (Royal College of Physicians, 2013), compared with around 15% of the general population (Office for National Statistics, 2019). Individuals with mental health conditions smoke more heavily (Coulthard et al., Reference Coulthard, Farrell, Singleton and Meltzer2002) and experience up to 18 years reduced life expectancy compared with the general population (Chang et al., Reference Chang, Hayes, Perera, Broadbent, Fernandes, Lee and Stewart2011; Royal College of Physicians, 2013). Much of this reduction can be explained by smoking-related diseases (Royal College of Physicians, 2013), making it important to understand the relationship between smoking and mental health.

It is often assumed that the association between mental health and smoking can be explained by a self-medication model – that is, symptoms of mental illness, or side effects of psychiatric medications, are alleviated by the chemical properties of tobacco (Desai, Seabolt, & Jann, Reference Desai, Seabolt and Jann2001; Khantzian, Reference Khantzian1997; Lerman et al., Reference Lerman, Audrain, Orleans, Boyd, Gold, Main and Caporaso1996; Levin, Wilson, Rose, & McEvoy, Reference Levin, Wilson, Rose and McEvoy1996). However, observational evidence cannot determine whether the association between smoking and mental health is causal or the result of confounding (Lawlor, Harbord, Sterne, Timpson, & Davey Smith, Reference Lawlor, Harbord, Sterne, Timpson and Davey Smith2008). Furthermore, traditional observational evidence cannot robustly identify the direction of causation (Lawlor et al., Reference Lawlor, Harbord, Sterne, Timpson and Davey Smith2008), and there is growing evidence to suggest that smoking may be a causal risk factor for poor mental health. The genome-wide association study (GWAS) of schizophrenia conducted by the Psychiatric Genomics Consortium (PGC) found that variants in the gene cluster CHRNA5-A3-B4 were associated with increased schizophrenia risk (Ripke et al., Reference Ripke, Neale, Corvin, Walters, Farh, Holmans and O'Donovan2014). These variants are known to be strongly associated with heaviness of smoking (Munafò et al., Reference Munafò, Timofeeva, Morris, Prieto-Merino, Sattar, Brennan and Davey Smith2012; Thorgeirsson et al., Reference Thorgeirsson, Geller, Sulem, Rafnar, Wiste, Magnusson and Stefansson2008; Tobacco Consortium, 2010; Ware, van den Bree, & Munafò, Reference Ware, van den Bree and Munafò2011). Therefore, one interpretation is a possible causal effect of smoking on schizophrenia (Gage & Munafò, Reference Gage and Munafò2015). Furthermore, there is evidence of genetic correlations between smoking, schizophrenia and depression (Hartz et al., Reference Hartz, Horton, Hancock, Baker, Caporaso, Chen and McNeil2018; Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Tian2019) warranting further investigation of possible causal effects. Prospective observational studies using related individuals to control for genetic and environmental confounding have suggested a dose–response effect of smoking on schizophrenia (Kendler, Lönn, Sundquist, & Sundquist, Reference Kendler, Lönn, Sundquist and Sundquist2015) and depression (Kendler et al., Reference Kendler, Neale, MacLean, Heath, Eaves and Kessler1993). Meta-analyses show further evidence for an increased relative risk of schizophrenia in smokers over non-smokers (Gurillo, Jauhar, Murray, & MacCabe, Reference Gurillo, Jauhar, Murray and MacCabe2015; Scott et al., Reference Scott, Matuschka, Niemelä, Miettunen, Emmerson and Mustonen2018) and a reduction in depressive symptoms following smoking cessation (Taylor, McNeill et al., Reference Taylor, McNeill, Girling, Farley, Lindson-Hawley and Aveyard2014). Although these studies suggest a potential causal effect, more robust methods are required to triangulate evidence and allow for stronger causal inference (Munafò & Davey Smith, Reference Munafò and Davey Smith2018).

One possible way to overcome bias from residual confounding and reverse causation is Mendelian randomisation (MR) (Davey Smith & Ebrahim, Reference Davey Smith and Ebrahim2003). This method uses genetic variants to proxy for an exposure in an instrumental variable analysis to estimate the causal effect on an outcome (Lawlor et al., Reference Lawlor, Harbord, Sterne, Timpson and Davey Smith2008). Previous MR studies have failed to show any clear evidence for an effect of smoking on depression (Bjørngaard et al., Reference Bjørngaard, Gunnell, Elvestad, Smith, Skorpen, Krokan and Romundstad2013; Taylor, Fluharty et al., Reference Taylor, Fluharty, Bjørngaard, Gabrielsen, Skorpen, Marioni and Munafò2014; Wium-Andersen, Ørsted, & Nordestgaard, Reference Wium-Andersen, Ørsted and Nordestgaard2015) and show suggestive but inconclusive evidence for an effect of smoking on schizophrenia (Gage et al., Reference Gage, Jones, Taylor, Burgess, Zammit and Munafò2017; Wium-Andersen et al., Reference Wium-Andersen, Ørsted and Nordestgaard2015). However, the genetic instruments for smoking used in these MR studies are limited, only capturing individual aspects of smoking behaviour and only having identified limited numbers of suitable genetic variants (Bjørngaard et al., Reference Bjørngaard, Gunnell, Elvestad, Smith, Skorpen, Krokan and Romundstad2013; Gage et al., Reference Gage, Jones, Taylor, Burgess, Zammit and Munafò2017; Taylor, Fluharty et al., Reference Taylor, Fluharty, Bjørngaard, Gabrielsen, Skorpen, Marioni and Munafò2014; Wium-Andersen et al., Reference Wium-Andersen, Ørsted and Nordestgaard2015). Furthermore, any instrument for smoking heaviness requires stratifying samples into smokers and non-smokers. Stratification is not possible using the most common MR method, two-sample MR. In two-sample MR, single-nucleotide polymorphism (SNP)-exposure and SNP-outcome effects are estimated in two independent samples and the effect sizes are obtained from GWAS summary statistics. Therefore, in this context stratification is not possible because GWAS summary data do not provide individual level data regarding smoking status.

Here we improve upon previous MR studies of smoking and mental health in two ways: (1) by using an updated instrument for smoking initiation which comprises 378 SNPs (Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Tian2019) and (2) we develop a novel genetic instrument for lifetime smoking exposure that takes into account smoking status (i.e. ever and never smokers), and among ever smokers takes into account smoking duration, heaviness and cessation. This instrument can therefore capture smoking behaviours beyond initiation and can be used in unstratified samples of smokers and non-smokers. By using the results of these different MR analyses together and triangulating across multiple methods, we aim to determine whether the observational associations between smoking, schizophrenia and depression are likely causal, and the directionality of these relationships.

Methods

Data sources and genetic instruments

Smoking initiation

The most recent GWAS of smoking initiation from the GSCAN consortium identified 378 conditionally independent genome-wide significant SNPs in a sample of 1 232 091 individuals of European ancestry (Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Tian2019). The 378 SNPs explain 2% of the variance in smoking initiation (Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Tian2019). When smoking initiation was the outcome, summary statistics without 23andMe (N = 599 289) were used.

Lifetime smoking

To capture smoking heaviness and duration as well as smoking initiation, we generated a lifetime smoking measure using data from the UK Biobank. The UK Biobank is a national health research resource of 502 647 participants aged 40–69 years, recruited from across the UK between 2006 and 2010 (http://www.ukbiobank.ac.uk). Our sample consisted of 462 690 individuals of European ancestry who had phenotype data and passed genotype inclusion criteria (54% female; mean age = 56.7 years; s.d. = 8.0 years). Overall, 30% of the sample had ever smoked (8% current smokers and 22% former smokers).

Measures of smoking. Smoking measures available in the UK Biobank were self-reported and collected at initial assessment. They included: smoking status (current, former, never – field 20116), age at initiation in years (fields 3436/2867), age at cessation in years (field 2897) and number of cigarettes smoked per day (fields 3456/2887). Anyone self-reporting to smoke more than 100 cigarettes per day was contacted for confirmation. Hand-rolled cigarette smokers were told 1 g of tobacco equates to one cigarette. We calculated duration of smoking and time since cessation. UK Biobank removed individuals smoking fewer than 1 or more than 150 cigarettes a day.

Construction of the lifetime smoking index. Following the method outlined by Leffondré, Abrahamowicz, Xiao, and Siemiatycki (Reference Leffondré, Abrahamowicz, Xiao and Siemiatycki2006), we combined the smoking measures into a lifetime smoking index along with a simulated half-life (τ) constant. Half-life captures the exponentially decreasing effect of smoking at a given time on health outcomes. The value of half-life was determined by simulating the effects of lifetime smoking on lung cancer and overall mortality in the UK Biobank. Both suggested the best fitting value as 18. For full details on construction of the lifetime smoking index see online Supplementary Note.

Genome-wide association study of lifetime smoking index. For full details of genotyping and exclusion procedures see online Supplementary Note. After excluding individuals who did not pass genotype exclusions and who had missing phenotype data, 462 690 individuals remained for the GWAS. Of these individuals, 249 318 were never smokers (54%), 164 649 were former smokers (36%) and 48 723 (11%) were current smokers. The mean value of lifetime smoking score was 0.359 (s.d. = 0.694). GWAS was conducted using the UK Biobank GWAS pipeline set-up for the MRC Integrative Epidemiology Unit (Elsworth et al., Reference Elsworth, Mitchell, Raistrick, Paternoster, Hemani and Gaunt2017). BOLT-LMM was used to conduct the analysis (Loh et al., Reference Loh, Tucker, Bulik-Sullivan, Vilhjálmsson, Finucane, Salem and Price2015), which accounts for population stratification and relatedness using linear mixed modelling. Genotyping chip and sex were included as covariates. As a sensitivity analysis, we reran the GWAS without controlling for genotype chip because the UK BiLEVE sub-sample (which used a different genotyping chip) was selected on the basis of smoking status. Genome-wide significant SNPs were selected at p < 5 × 10−8 and were clumped to ensure independence at linkage disequilibrium (LD) r 2 = 0.001 and a distance of 10 000 kb using the TwoSampleMR package (Hemani et al., Reference Hemani, Zheng, Elsworth, Wade, Haberland, Baird and Langdon2018).

Instrument generated. Our GWAS of lifetime smoking identified 126 independent, genome-wide significant SNPs (see Fig. 1 and online Supplementary Table S1). A standard deviation increase in the lifetime smoking score is equivalent to an individual smoking 20 cigarettes a day for 15 years and stopping 17 years ago or an individual smoking 60 cigarettes a day for 13 years and stopping 22 years ago. We validated our instrument in an independent sample and using positive control outcomes of lung cancer, coronary heart disease and hypomethylation at the aryl-hydrocarbon receptor repressor site cg05575921 because smoking has been robustly shown to predict these outcomes (WHO, 2011; Zeilinger et al., Reference Zeilinger, Kühnel, Klopp, Baurecht, Kleinschmidt, Gieger and Illig2013). For further details of these validations see the online Supplementary Note.

Fig. 1. Manhattan plot of genome-wide association study of lifetime smoking index (N = 462 690). The x-axis represents chromosomal position and the y-axis represents −log10 p value for the association of each SNP with the lifetime smoking index using an additive model and linear regression. The dashed line indicates the genome-wide level of significance (p < 5 × 10−8) and genome-wide significant SNPs are indicated in red.

Schizophrenia

We used summary data from the PGC consortium GWAS, which comprises 36 989 cases and 113 075 controls of European and East Asian ancestry (Ripke et al., Reference Ripke, Neale, Corvin, Walters, Farh, Holmans and O'Donovan2014). Cases were a combination of individuals with schizophrenia and schizoaffective disorder, mostly diagnosed by clinicians, but some samples used research-based assessment. The ascertainment method did not affect GWAS results (Ripke et al., Reference Ripke, Neale, Corvin, Walters, Farh, Holmans and O'Donovan2014). A sensitivity analysis was performed using the GWAS summary data meta-analysed with a further 11 260 cases and 24 542 controls (Pardiñas et al., Reference Pardiñas, Holmans, Pocklington, Escott-Price, Ripke, Carrera and Walters2018). The genetic instrument for schizophrenia came from the PGC GWAS, which identified 114 independent SNPs at 108 loci explaining around 3.4% of the variance in schizophrenia liability (Ripke et al., Reference Ripke, Neale, Corvin, Walters, Farh, Holmans and O'Donovan2014).

Depression

For depression, we used GWAS summary data from the PGC for major depression, which comprises 130 664 major depression cases and 330 470 controls of European ancestry (Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018). Cases were either diagnosed with major depressive disorder (MDD) on inpatient or medical health records or self-reported having a diagnosis or treatment for depression. Therefore, the authors use the term major depression over diagnosed MDD (Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018). 23andMe data (75 607 cases and 231 747 controls) were excluded when major depression was the outcome because genome-wide summary statistics are not available with 23andMe data included. The genetic instrument for major depression from the PGC GWAS was 40 genome-wide significant SNPs which explain 1.9% of the variance in liability (Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018).

Statistical analysis

MR analyses were run using the MR Base R package (Hemani et al., Reference Hemani, Zheng, Elsworth, Wade, Haberland, Baird and Langdon2018; R. Core Team, 2014) and compared across five different methods: inverse-variance weighted, MR Egger (Bowden, Davey Smith, & Burgess, Reference Bowden, Davey Smith and Burgess2015), weighted median (Bowden, Davey Smith, Haycock, & Burgess, Reference Bowden, Davey Smith, Haycock and Burgess2016), weighted mode (Hartwig, Smith, & Bowden, Reference Hartwig, Smith and Bowden2017) and MR RAPS (Zhao, Wang, Hemani, Bowden, & Small, Reference Zhao, Wang, Hemani, Bowden and Small2018). Each of these methods makes slightly different assumptions about the nature of pleiotropy and therefore a roughly consistent point estimate across the multiple methods provides the strongest evidence of causal inference (Lawlor, Tilling, & Davey Smith, Reference Lawlor, Tilling and Davey Smith2016), with the Inverse variance weighted (IVW) method being the main analysis and each other method providing a sensitivity analysis. For additional details on each of the methods see online Supplementary Table S12. Two-sample MR analysis was run bi-directionally, first with smoking as the exposure and then as the outcome. To test the suitability of the MR Egger method, the I 2GX statistic was calculated to quantify the degree of regression dilution bias due to measurement error of SNP-exposure effects (Bowden et al., Reference Bowden, Del Greco, Minelli, Davey Smith, Sheehan and Thompson2016). We also calculated the mean F statistic as an indicator of instrument strength. Steiger filtering was conducted to confirm the direction of effect (Hemani, Tilling, & Davey Smith, Reference Hemani, Tilling and Davey Smith2017). If a SNP from the instrument was unavailable in the outcome, an attempt to find proxies was made with a minimum LD r 2 = 0.8 and palindromic SNPs were aligned with Minor allele frequency < 0.3.

Results

Bi-directional MR analyses provided strong evidence that higher lifetime smoking increases risk of both schizophrenia [IVW: odds ratio (OR) 2.27, 95% confidence interval (CI) 1.67–3.08, p < 0.001] and depression (IVW: OR 1.99, 95% CI 1.71–2.32, p < 0.001) (see Table 1), with consistent direction of effect across all five MR methods. The same was seen for smoking initiation as the instrument on schizophrenia (IVW: OR 1.53, 95% CI 1.35–1.74, p < 0.001) and depression (IVW: OR 1.54, 95% CI 1.44–1.64, p < 0.001) (see Table 1). MR Egger results are the least reliable due to low I 2GX (see online Supplementary Table S3). There was also evidence of a consistent but smaller effect of higher genetic liability for schizophrenia on lifetime smoking (IVW: β = 0.022, 95% CI 0.005–0.038, p = 0.009) and of genetic liability for depression on lifetime smoking (IVW: β = 0.091, 95% CI 0.027–0.155, p = 0.005) (see Table 2). These effects remained for depression on smoking initiation (IVW: β = 0.083, 95% CI 0.039–0.127, p < 0.001) but became even weaker and inconsistent across the different methods for schizophrenia on smoking initiation (IVW: β = 0.010, 95% CI 0.000–0.021, p = 0.04) (see Table 2). Overall, the effect was larger and more consistent for depression than for schizophrenia. We saw similar effects using the more recent meta-analysed GWAS for schizophrenia with an additional 11 260 cases (see online Supplementary Table S9). There was evidence of significant heterogeneity (see online Supplementary Table S4) but MR Egger intercepts suggest directional pleiotropy is not biasing the estimates (see online Supplementary Table S5) and Steiger filtering supported the conclusion that the effects operate in both directions (see online Supplementary Table S6). Given that variants from the CHRNA5-A3-B4 gene complex were identified in both GWAS of lifetime smoking and schizophrenia, we conducted a sensitivity analysis removing these variants and the effects remained consistent (see online Supplementary Table S10).

Table 1. Two-sample MR analyses of the effect of smoking exposure on schizophrenia and depression

SIMEX, simulation extrapolation, MR RAPS, robust adjusted profile score.

Note: Due to low regression dilution I 2GX for lifetime smoking and smoking initiation (see online Supplementary Table S3), MR Egger SIMEX estimates should be interpreted with caution. SIMEX-corrected estimates are unweighted.

Table 2. Two-sample MR analyses of the effect of schizophrenia and depression on smoking

SIMEX, simulation extrapolation, MR RAPS, robust adjusted profile score.

Note: Due to low regression dilution I 2GX (see online Supplementary Table S3), MR Egger estimates could not be conducted apart from for the effect of schizophrenia on smoking initiation, where a weighted MR Egger SIMEX was conducted. Smoking initiation scores are given in βs by GSCAN calculated from the meta-analysed z-score statistic by assuming a prevalence for the binary trait (Liu et al., Reference Liu, Jiang, Wedow, Li, Brazel, Chen and Tian2019).

A total of 14 260 cases and 15 480 controls from the UK Biobank were included in the latest GWAS for major depression (Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018) which could lead to bias from sample overlap (Burgess et al., Reference Burgess, Scott, Timpson, Smith, Thompson and Consortium2015). Therefore, this analysis was repeated using summary data from an earlier GWAS of major depression (Ripke et al., Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen and Sullivan2013), which showed a consistent direction of effect with weaker statistical evidence, possibly due to reduced sample size (N = 18 759) (see online Supplementary Table S7). Second, we repeated this analysis with the recent meta-analysed GWAS which includes broader definitions of depression including self-declared (Howard et al., Reference Howard, Adams, Clarke, Hafferty, Gibson, Shirali and McIntosh2019). These summary data were available with UK Biobank and 23andMe removed leaving a sample size of 500 199 (170 756 cases and 329 443 controls). The results of this analysis are presented in online Supplementary Table S11. They are relatively consistent with the main analysis, showing slightly larger effects of depression on smoking and slightly reduced effects of smoking on depression. Bi-directional MR analyses were repeated using the GWAS of lifetime smoking without controlling for genotyping chip and the same pattern of results was observed (see online Supplementary Table S8). Further sensitivity tests were conducted and are presented in online Supplementary Figs S9–S20.

Discussion

We conducted a GWAS of lifetime smoking exposure which provides a novel genetic instrument that can be used in two-sample MR of summary data without the need to stratify on smoking status. We used this novel genetic instrument along with a new instrument of smoking initiation to explore possible causal pathways between smoking, schizophrenia and depression. The two-sample MR results provide strong evidence that smoking is a risk factor for both schizophrenia and depression. This supports prospective observational evidence controlling for genetic confounding (Kendler et al., Reference Kendler, Neale, MacLean, Heath, Eaves and Kessler1993, Reference Kendler, Lönn, Sundquist and Sundquist2015), as well as meta-analyses of observational studies (Gurillo et al., Reference Gurillo, Jauhar, Murray and MacCabe2015; Taylor, McNeill et al., Reference Taylor, McNeill, Girling, Farley, Lindson-Hawley and Aveyard2014) (although it should be noted that these meta-analyses include estimates not adjusted for known confounders e.g. cannabis use). Effect sizes were similar to a more recent meta-analysis that did adjust for multiple confounders and found a two-fold increased risk of schizophrenia in smokers compared with non-smokers (Scott et al., Reference Scott, Matuschka, Niemelä, Miettunen, Emmerson and Mustonen2018). Some studies which adjust for potential confounders find the effect of smoking attenuates to the null (Jones et al., Reference Jones, Gage, Heron, Hickman, Lewis, Munafò and Zammit2018) or even becomes protective (Zammit et al., Reference Zammit, Allebeck, Dalman, Lundberg, Hemmingsson and Lewis2003), demonstrating that there are likely to be substantial confounding effects in observational studies. Previous MR studies have not found clear evidence to support smoking as a risk factor for either schizophrenia or depression (Bjørngaard et al., Reference Bjørngaard, Gunnell, Elvestad, Smith, Skorpen, Krokan and Romundstad2013; Gage et al., Reference Gage, Jones, Taylor, Burgess, Zammit and Munafò2017; Taylor, Fluharty et al., Reference Taylor, Fluharty, Bjørngaard, Gabrielsen, Skorpen, Marioni and Munafò2014; Wium-Andersen et al., Reference Wium-Andersen, Ørsted and Nordestgaard2015), but our approach offers greater power, captures multiple aspects of smoking behaviour and enables two-sample MR analysis using summary data in unstratified samples. However, it is not possible to precisely estimate from our results what proportion of the observational association between smoking, schizophrenia and depression is causal, or the population attributional fraction of these disorders due to smoking.

In support of the self-medication hypothesis (Khantzian, Reference Khantzian1997), we found evidence to suggest that genetic liability for schizophrenia and depression increases lifetime smoking. This supports previous observational evidence (Desai et al., Reference Desai, Seabolt and Jann2001; Lerman et al., Reference Lerman, Audrain, Orleans, Boyd, Gold, Main and Caporaso1996; Levin et al., Reference Levin, Wilson, Rose and McEvoy1996) and might explain why smoking rates remain so high amongst individuals with schizophrenia and depression compared with the general population (Cook et al., Reference Cook, Wayne, Kafali, Liu, Shu and Flores2014). However, evidence was stronger for self-medication effects in depression than schizophrenia and when using smoking initiation as the outcome rather than lifetime smoking, effects attenuated to the null. Therefore, maybe any self-medication effects of schizophrenia are only on heaviness and duration of smoking (captured by the lifetime smoking index) rather than initiation. However, it is important to note that the effects might be weaker because MR methods typically capture the long-term effects of exposures (Labrecque & Swanson, Reference Labrecque and Swanson2018) with self-medication potentially being more acute.

Alternatively the effects of schizophrenia and depression on lifetime smoking might be explained by the misattribution hypothesis. This proposes that smokers misattribute the ability of cigarettes to relieve withdrawal, to their ability to relieve symptoms of psychological distress (Parrott, Reference Parrott1999, Reference Parrott2003). For example, withdrawal symptoms include depressed mood, anxiety and irritability, and smoking a cigarette alleviates those symptoms (Hughes, Reference Hughes2007). Since many withdrawal symptoms are similar to the negative symptoms of schizophrenia and mood symptoms in depression, their alleviation by smoking could give rise to the strong belief that smoking helps to alleviate mental health symptoms.

A potential biological mechanism for the bi-directional causal effects of smoking, schizophrenia and depression could be neuroadaptations in the dopaminergic and serotonin systems. Nicotine acts on nicotinic cholinergic receptors in the brain stimulating the release of neurotransmitters including dopamine and serotonin (Benowitz, Reference Benowitz2010). Dopamine and serotonin dysfunction have been implicated in the aetiology of schizophrenia and depression, respectively (Howes, McCutcheon, & Stone, Reference Howes, McCutcheon and Stone2015; Jakubovski, Varigonda, Freemantle, Taylor, & Bloch, Reference Jakubovski, Varigonda, Freemantle, Taylor and Bloch2015). It is plausible, therefore, that disruption of these pathways has a causal effect on these disorders. Alternatively, cannabis use could be a mediating mechanism for the effects of smoking on schizophrenia and depression. In prospective studies, cigarette smoking has been shown to increase risk of cannabis dependence (Hindocha et al., Reference Hindocha, Shaban, Freeman, Das, Gale, Schafer and Curran2015). There is strong evidence suggesting that cannabis use increases the risk of psychosis and affective disorders (Moore et al., Reference Moore, Zammit, Lingford-Hughes, Barnes, Jones, Burke and Lewis2007). This vertical pleiotropy does not violate the assumptions of MR, but simply means there are intermediate mechanisms operating between the exposure and the outcome. However, the strong effects we observe for lifetime smoking suggest that any mediating influence of cannabis use is likely to only partially account for these effects, given the relatively low prevalence of cannabis use [e.g. annual prevalence in the UK of ~7% in 2010 (United Nations Office on Drugs and Crime (UNODC), 2011)]. Multivariable MR analysis of tobacco and cannabis use would help resolve this question.

Finally, although we are using the method of MR to provide stronger causal inference, the best evidence of a causal effect comes from many corroborating lines of evidence from study designs with diverse assumptions. Looking at our current results alongside previous literature, we conclude this strengthens evidence for an effect of smoking on increased risk of depression and schizophrenia. Future work should attempt to elucidate the underlying mechanisms with a hope to intervene, inform public health messages or further advance our knowledge on the aetiology of mental illness. In particular, it will be important to consider other constituents of tobacco smoke to determine whether it is exposure to nicotine or some other constituent that increases risk of schizophrenia and depression. This is particularly important in the context of the recent growth in electronic cigarette use.

Strengths and limitations

Our study is the first to generate a genetic instrument for lifetime smoking behaviour in a large sample, which allows the use of two-sample MR with summary data from unstratified samples to answer questions about the association between smoking and other health outcomes. However, there are some limitations which should be noted. First, there is evidence to suggest that even after seemingly controlling for population structure in GWAS of samples as large as the UK Biobank, coincident apparent structure remains (Haworth et al., Reference Haworth, Mitchell, Corbin, Wade, Dudding, Budu-Aggrey and Timpson2019). This might confound the association between smoking and mental health, increasing the risk of false positives. As independent samples with adequate sample size become available, the influence of structure should be further explored. However, it is reassuring that our instrument predicted lifetime smoking in an independent replication sample, where such issues would not arise in the same manner.

Second, sample overlap in two-sample MR can bias results towards the observational estimate (Burgess et al., Reference Burgess, Scott, Timpson, Smith, Thompson and Consortium2015). There was some sample overlap between the major depression GWAS (Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018) and the UK Biobank (used to derive the lifetime smoking instrument and included in the smoking initiation GWAS) meaning that the effects could be inflated. Therefore, a sensitivity analysis was conducted using a previous GWAS of major depression (Ripke et al., Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen and Sullivan2013) which showed the same direction of effect despite lower power. This gives us confidence in the bi-directional effects of smoking and depression, despite sample overlap. This sensitivity analysis also addresses recent concerns over the specificity of the most recent GWAS for major depression (Cai et al., Reference Cai, Revez, Adams, Andlauer, Breen, Byrne and Flint2019). Comparing self-reported ‘seeking help for depression’ with DSM-diagnosed MDD yielded different results (Cai et al., Reference Cai, Revez, Adams, Andlauer, Breen, Byrne and Flint2019). However, the earlier GWAS of depression did use Diagnostic and statistical manual diagnosed cases only (Ripke et al., Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen and Sullivan2013) and showed the same direction of effect despite lower power.

Third, including multiple aspects of smoking behaviour could have introduced more potential for horizontal pleiotropy. The more diffuse the definition of smoking, the more lifestyle factors might be correlated, making it especially important to test for horizontal pleiotropy. We conducted multiple sensitivity analyses (which all make different and largely uncorrelated assumptions), and all demonstrated the same direction of effect. This increases our confidence that the results are not biased by pleiotropy. Furthermore, MR Egger intercepts did not show evidence of directional pleiotropy for schizophrenia or depression. However, further work should still attempt to understand the biological mechanisms underpinning the association to reduce the likelihood of pleiotropic effects.

Fourth, schizophrenia and depression are disorders with an average age of onset around early to mid-adulthood (Office for National Statistics, 2018; World Health Organization, 2018). Our measure of lifetime smoking was generated using participants in the UK Biobank aged over 40 years. Therefore, the causal pathway from smoking to schizophrenia and depression risk might initially seem unclear. However, we were not using participants' smoking behaviour at age 40, but rather retrospective lifetime smoking behaviour from age at initiation. It is plausible that smoking behaviour in earlier life could increase risk of later mental health outcomes and exacerbate symptoms, consequently causing more smokers than non-smokers to seek a diagnosis. Furthermore, we saw consistent effects when using smoking initiation as the exposure. Individuals are more likely to initiate smoking prior to the average age at onset of schizophrenia and depression, with 90% of individuals initiating before 18 years of age (Centers for Disease Control and Prevention, 2018). There has also been recent debate in the field about the interpretation of time varying exposures in MR and one way to minimise bias is to use average SNP effects on phenotype across time, as we have done here with our lifetime smoking instrument (Labrecque & Swanson, Reference Labrecque and Swanson2018). We recommend that future studies wishing to examine the effects of smoking on an outcome use multiple instruments for smoking behaviour with consistent evidence across multiple instruments providing the strongest evidence of a causal effect.

Fifth, there is a high degree of zero inflation in the distribution of our lifetime smoking index scores with 54% of the sample being never smokers and therefore receiving a score of zero. We decided not to transform the variable given the desire to have interpretable effect sizes for MR and we decided not to exclude never smokers because our instrument is designed to be used in two-sample MR without the need to stratify into smokers and non-smokers. Despite the zero inflation, we see similar effects for lifetime smoking and smoking initiation suggesting that this has not impaired the score. Sixth, the lifetime smoking score was simulated using all-cause mortality and lung cancer as outcomes. The pattern of association between smoking and lung cancer risk compared with risk for schizophrenia and depression is likely to be different. However, increased mortality amongst individuals with schizophrenia and depression is in large part due to smoking-related mortality (Royal College of Physicians, 2013). The effects were modelled on all-cause mortality and lung cancer but no difference to the best fitting value of half-life was observed. We hope that using all-cause mortality as an outcome makes the lifetime smoking instrument broadly applicable to exploring multiple outcomes. Furthermore, the same effects are observed using smoking initiation as the exposure, which does not include the simulated variable.

Finally, there is known selection bias in the UK Biobank sample, with participants being more highly educated, less likely to be a smoker and overall healthier than the general UK population (Munafo, Tilling, Taylor, Evans, & Davey Smith, Reference Munafo, Tilling, Taylor, Evans and Davey Smith2017). Of the 9 million individuals contacted, only ~5% consented to take part (Munafo et al., Reference Munafo, Tilling, Taylor, Evans and Davey Smith2017). Due to the lack of representativeness in the UK Biobank sample, prevalence and incidence rates will not reflect underlying population levels and there is potential for collider bias. If both smoking and liability for schizophrenia and depression reduce the likelihood of participating in the UK Biobank, then this would induce a negative correlation between schizophrenia or depression and smoking. That is the opposite of the effects observed, suggesting our estimates may, if anything, be conservative.

Conclusion

In conclusion, we improved upon previous MR studies of smoking and mental health by using an updated instrument of smoking initiation and by developing a novel genetic instrument of lifetime smoking that can be used in two-sample MR of summary data without stratifying on smoking status or reducing power. In two-sample MR analysis, there was evidence of an effect between smoking and mental health in both directions; however, the however, the effects for depression on smoking were stronger than for schizophrenia on smoking. Strong effects of smoking as a risk factor emphasises the detrimental public health consequences of smoking, especially for mental health, and the need to reduce smoking prevalence not only to reduce the burden of physical illness, but also to reduce the burden of mental illness.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291719002678.

Data

Lifetime smoking GWAS summary data is available for download at https://doi.org/10.5523/bris.10i96zb8gm0j81yz0q6ztei23d.

Acknowledgements

We are grateful to the participants of the UK Biobank and the individuals who contributed to each of the previous GWAS analyses conducted as well as all the research staff who worked on the data collection. We are extremely grateful to all the families who took part in the Avon Longitudinal Study of Parents and Children (ALSPAC) study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses.

Financial support

REW, RCR, RBL, HMS, SZ, GH, GDS and MRM are all members of the MRC Integrative Epidemiology Unit at the University of Bristol funded by the MRC: http://www.mrc.ac.uk (MC_UU_00011/1, MC_UU_00011/5, MC_UU_00011/7). This study was supported by the NIHR Biomedical Research Centre at the University Hospitals Bristol NHS Foundation Trust and the University of Bristol. The views expressed in this publication are those of the authors and not necessarily those of the NHS, the National Institute for Health Research or the Department of Health and Social Care. GMJT is funded by Cancer Research UK Population Researcher Postdoctoral Fellowship award (reference: C56067/A21330). This work was also supported by CRUK (grant number C18281/A19169) and the ESRC (grant number ES/N000498/1). GH is supported by the Wellcome Trust (208806/Z/17/Z). The UK Medical Research Council and the Wellcome Trust (Grant ref: 102215/2/13/2) and the University of Bristol provide core support for ALSPAC. The Accessible Resource for Integrated Epigenomics Studies (ARIES) was funded by the UK Biotechnology and Biological Sciences Research Council (BB/I025751/1 and BB/I025263/1). This publication is the work of the authors and REW will serve as guarantor for the contents of this paper. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). This research has been conducted using the UK Biobank Resource under Application Number 9142.

Conflict of interest

REW declares no conflicts of interest on behalf of all co-authors.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. Ethical approval for the replication analysis in the ALSPAC cohort was obtained from the ALSPAC Ethics and Law Committee and the Local Research Ethics Committees. Consent for ALSPAC biological samples has been collected in accordance with the Human Tissue Act (2004).

References

Benowitz, N. L. (2010). Nicotine addiction. The New England Journal of Medicine, 362, 22952303.CrossRefGoogle ScholarPubMed
Bjørngaard, J. H., Gunnell, D., Elvestad, M. B., Smith, G. D., Skorpen, F., Krokan, H., … Romundstad, P. (2013). The causal role of smoking in anxiety and depression: A Mendelian randomization analysis of the HUNT study. Psychological Medicine, 43, 711719.CrossRefGoogle ScholarPubMed
Bowden, J., Davey Smith, G., & Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. International Journal of Epidemiology, 44, 512525.CrossRefGoogle ScholarPubMed
Bowden, J., Davey Smith, G., Haycock, P. C., & Burgess, S. (2016). Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genetic Epidemiology, 40, 304314.CrossRefGoogle ScholarPubMed
Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N. A., & Thompson, J. R. (2016). Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. International Journal of Epidemiology, 45, 19611974.Google ScholarPubMed
Burgess, S., Scott, R. A., Timpson, N. J., Smith, G. D., Thompson, S. G., & Consortium, E.-I. (2015). Using published data in Mendelian randomization: A blueprint for efficient identification of causal risk factors. European Journal of Epidemiology, 30, 543552.CrossRefGoogle ScholarPubMed
Byers, A. L., Vittinghoff, E., Lui, L-Y, Hoang, T., Blazer, D. G., Covinsky, K. E., … Fredman, L. (2012). Twenty-year depressive trajectories among older women. Archives of General Psychiatry, 69, 10731079.CrossRefGoogle ScholarPubMed
Cai, N., Revez, J. A., Adams, M. J., Andlauer, T. F. M., Breen, G., Byrne, E. M., … Flint, J. (2019). Minimal phenotyping yields GWAS hits of low specificity for major depression. bioRxiv, article id 440735, 34.Google Scholar
Centers for Disease Control and Prevention (2018). Smoking and Tobacco Use; Tobacco-Related Disparities; Tobacco Use Among Adults with Mental Illness and Substance Use Disorders. Centers for Disease Control and Prevention. Retrieved from https://www.cdc.gov/tobacco/disparities/mental-illness-substance-use/index.htm (Accessed 24 April 2018).Google Scholar
Chang, C-K, Hayes, R. D., Perera, G., Broadbent, M. T. M., Fernandes, A. C., Lee, W. E., … Stewart, R. (2011). Life expectancy at birth for people with serious mental illness and other major disorders from a secondary mental health care case register in London. PLoS ONE, 6, e19590.CrossRefGoogle Scholar
Cook, B. L., Wayne, G. F., Kafali, E. N., Liu, Z., Shu, C., & Flores, M. (2014). Trends in smoking among adults with mental illness and association between mental health treatment and smoking cessation. JAMA, 311, 172182.CrossRefGoogle ScholarPubMed
Coulthard, M., Farrell, M., Singleton, N., & Meltzer, H. (2002). Tobacco, alcohol and drug use and mental health. London: The Stationery Office.Google Scholar
Davey Smith, G., & Ebrahim, S. (2003). ‘Mendelian randomization’: Can genetic epidemiology contribute to understanding environmental determinants of disease? International Journal of Epidemiology, 32, 122.CrossRefGoogle Scholar
Desai, H. D., Seabolt, J., & Jann, M. W. (2001). Smoking in patients receiving psychotropic medications. CNS Drugs, 15, 469494.CrossRefGoogle ScholarPubMed
Elsworth, B. L., Mitchell, R., Raistrick, C. A., Paternoster, L., Hemani, G., & Gaunt, T. (2017). MRC IEU UK Biobank GWAS pipeline version 1. data.bris. Retrieved from https://doi.org/10.5523/bris.2fahpksont1zi26xosyamqo8rr.CrossRefGoogle Scholar
Gage, S. H., & Munafò, M. R. (2015). Rethinking the association between smoking and schizophrenia. The Lancet Psychiatry, 2, 118.CrossRefGoogle Scholar
Gage, S. H., Jones, H. J., Taylor, A. E., Burgess, S., Zammit, S., & Munafò, M. R. (2017). Investigating causality in associations between smoking initiation and schizophrenia using Mendelian randomization. Scientific Reports, 7, 40653.CrossRefGoogle ScholarPubMed
Gurillo, P., Jauhar, S., Murray, R. M., & MacCabe, J. H. (2015). Does tobacco use cause psychosis? Systematic review and meta-analysis. The Lancet Psychiatry, 2, 718725.CrossRefGoogle ScholarPubMed
Hartwig, F. P., Smith, G. D., & Bowden, J. (2017). Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. International Journal of Epidemiology, 46(6), 19851998.CrossRefGoogle ScholarPubMed
Hartz, S. M., Horton, A. C., Hancock, D. B., Baker, T. B., Caporaso, N. E., Chen, L-S, … McNeil, D. W. (2018). Genetic correlation between smoking behaviors and schizophrenia. Schizophrenia Research, 194, 8690.CrossRefGoogle Scholar
Haworth, S., Mitchell, R., Corbin, L., Wade, K. H., Dudding, T., Budu-Aggrey, A., & Timpson, N. (2019). Apparent latent structure within the UK Biobank sample has implications for epidemiological analysis. Nature Communications, 10(1), article= 333, 1–9.CrossRefGoogle ScholarPubMed
Hemani, G., Tilling, K., & Davey Smith, G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genetics, 13, e1007081.CrossRefGoogle ScholarPubMed
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., & Langdon, R. (2018). The MR-base platform supports systematic causal inference across the human phenome. eLife, 7, e34408.CrossRefGoogle ScholarPubMed
Hindocha, C., Shaban, N. D., Freeman, T. P., Das, R. K., Gale, G., Schafer, G., … Curran, H. V. (2015). Associations between cigarette smoking and cannabis dependence: A longitudinal study of young cannabis users in the United Kingdom. Drug and Alcohol Dependence, 148, 165171.CrossRefGoogle ScholarPubMed
Howard, D. M., Adams, M. J., Clarke, T-K, Hafferty, J. D., Gibson, J., Shirali, M., … McIntosh, A. M. (2019). Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience, 22, 343.CrossRefGoogle ScholarPubMed
Howes, O., McCutcheon, R., & Stone, J. (2015). Glutamate and dopamine in schizophrenia: An update for the 21st century. Journal of Psychopharmacology, 29, 97115.CrossRefGoogle ScholarPubMed
Hughes, J. R. (2007). Effects of abstinence from tobacco: Valid symptoms and time course. Nicotine & Tobacco Research, 9, 315327.CrossRefGoogle ScholarPubMed
Jakubovski, E., Varigonda, A. L., Freemantle, N., Taylor, M. J., & Bloch, M. H. (2015). Systematic review and meta-analysis: Dose-response relationship of selective serotonin reuptake inhibitors in major depressive disorder. American Journal of Psychiatry, 173, 174183.CrossRefGoogle ScholarPubMed
Jones, H. J., Gage, S. H., Heron, J., Hickman, M., Lewis, G., Munafò, M. R., & Zammit, S. (2018). Association of combined patterns of tobacco and cannabis use in adolescence with psychotic experiences. JAMA Psychiatry, 75, 240246.CrossRefGoogle ScholarPubMed
Kendler, K. S., Neale, M. C., MacLean, C. J., Heath, A. C., Eaves, L. J., & Kessler, R. C. (1993). Smoking and major depression: A causal analysis. Archives of General Psychiatry, 50, 3643.CrossRefGoogle ScholarPubMed
Kendler, K. S., Lönn, S. L., Sundquist, J., & Sundquist, K. (2015). Smoking and schizophrenia in population cohorts of Swedish women and men: A prospective co-relative control study. American Journal of Psychiatry, 172, 10921100.CrossRefGoogle ScholarPubMed
Khantzian, E. (1997). The self-medication hypothesis of substance use disorders: A reconsideration and recent applications. Harvard Review of Psychiatry, 4, 231244.CrossRefGoogle ScholarPubMed
Labrecque, J. A., & Swanson, S. A. (2018). Interpretation and potential biases of Mendelian randomization estimates with time-varying exposures. American Journal of Epidemiology, 188, 231238.CrossRefGoogle Scholar
Lasser, K., Boyd, J. W., Woolhandler, S., Himmelstein, D. U., McCormick, D., & Bor, D. H. (2000). Smoking and mental illness: A population-based prevalence study. JAMA, 284, 26062610.CrossRefGoogle ScholarPubMed
Lawlor, D. A., Harbord, R. M., Sterne, J. A., Timpson, N., & Davey Smith, G. (2008). Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Statistics in Medicine, 27, 11331163.CrossRefGoogle ScholarPubMed
Lawlor, D. A., Tilling, K., & Davey Smith, G. (2016). Triangulation in aetiological epidemiology. International Journal of Epidemiology, 45, 18661886.Google ScholarPubMed
Lawrence, D., Mitrou, F., & Zubrick, S. R. (2009). Smoking and mental illness: results from population surveys in Australia and the United States. BMC Public Health, 9, 285.CrossRefGoogle ScholarPubMed
Leffondré, K., Abrahamowicz, M., Xiao, Y., & Siemiatycki, J. (2006). Modelling smoking history using a comprehensive smoking index: application to lung cancer. Statistics in Medicine, 25, 41324146.CrossRefGoogle ScholarPubMed
Lerman, C., Audrain, J., Orleans, C. T., Boyd, R., Gold, K., Main, D., & Caporaso, N. (1996). Investigation of mechanisms linking depressed mood to nicotine dependence. Addictive Behaviors, 21, 919.CrossRefGoogle ScholarPubMed
Leung, J., Gartner, C., Dobson, A., Lucke, J., & Hall, W. (2011). Psychological distress is associated with tobacco smoking and quitting behaviour in the Australian population: Evidence from national cross-sectional surveys. Australian & New Zealand Journal of Psychiatry, 45, 170178.CrossRefGoogle ScholarPubMed
Levin, E. D., Wilson, W., Rose, J. E., & McEvoy, J. (1996). Nicotine–haloperidol interactions and cognitive performance in schizophrenics. Neuropsychopharmacology, 15, 429.CrossRefGoogle ScholarPubMed
Liu, M., Jiang, Y., Wedow, R., Li, Y., Brazel, D. M., Chen, F., … Tian, C. (2019). Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nature Genetics, 51, 237244.CrossRefGoogle ScholarPubMed
Loh, P-R, Tucker, G., Bulik-Sullivan, B. K., Vilhjálmsson, B. J., Finucane, H. K., Salem, R. M., … Price, A. L. (2015). Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nature Genetics, 47, 284290.CrossRefGoogle ScholarPubMed
McClave, A. K., McKnight-Eily, L. R., Davis, S. P., & Dube, S. R. (2010). Smoking characteristics of adults with selected lifetime mental illnesses: Results from the 2007 National Health Interview Survey. American Journal of Public Health, 100, 24642472.CrossRefGoogle ScholarPubMed
Moore, T. H., Zammit, S., Lingford-Hughes, A., Barnes, T. R., Jones, P. B., Burke, M., & Lewis, G. (2007). Cannabis use and risk of psychotic or affective mental health outcomes: A systematic review. The Lancet, 370, 319328.CrossRefGoogle ScholarPubMed
Munafò, M. R., & Davey Smith, G. (2018). Robust research needs many lines of evidence. Nature, 553, 399401.CrossRefGoogle ScholarPubMed
Munafò, M. R., Timofeeva, M. N., Morris, R. W., Prieto-Merino, D., Sattar, N., Brennan, P., & Davey Smith, G. (2012). Association between genetic variants on chromosome 15q25 locus and objective measures of tobacco exposure. JNCI: Journal of the National Cancer Institute, 104, 740748.CrossRefGoogle ScholarPubMed
Munafo, M. R., Tilling, K., Taylor, A. E., Evans, D. M., & Davey Smith, G. (2017). Collider scope: When selection bias can substantially influence observed associations. International Journal of Epidemiology, 47(1), 226235.CrossRefGoogle Scholar
Office for National Statistics (2018). Measuring national well-being: Domains and measures. Office for National Statistics. Retrieved from https://www.ons.gov.uk/peoplepopulationandcommunity/wellbeing/datasets/measuringnationalwellbeingdomainsandmeasures (Accessed 19 December 2018).Google Scholar
Office for National Statistics (2019). Adult smoking habits in the UK: 2018. Office for National Statistics. Retrieved from https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/healthandlifeexpectancies/bulletins/adultsmokinghabitsingreatbritain/2018 (Accessed 27 August 2019).Google Scholar
Pardiñas, A. F., Holmans, P., Pocklington, A. J., Escott-Price, V., Ripke, S., Carrera, N., & Walters, J. T. R. (2018). Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nature Genetics, 50, 381389.CrossRefGoogle ScholarPubMed
Parrott, A. C. (1999). Does cigarette smoking cause stress? American Psychologist, 54, 817.CrossRefGoogle ScholarPubMed
Parrott, A. C. (2003). Cigarette-derived nicotine is not a medicine. The World Journal of Biological Psychiatry, 4, 4955.CrossRefGoogle Scholar
R. Core Team (2014). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. 2013. ISBN 3-900051-07-0.Google Scholar
Ripke, S., Wray, N. R., Lewis, C. M., Hamilton, S. P., Weissman, M. M., Breen, G., … Sullivan, P. F. (2013). A mega-analysis of genome-wide association studies for major depressive disorder. Molecular Psychiatry, 18, 497511.Google ScholarPubMed
Ripke, S., Neale, B. M., Corvin, A., Walters, J. T., Farh, K-H, Holmans, P. A., … O'Donovan, M. C. (2014). Biological insights from 108 schizophrenia-associated genetic loci. Nature, 511, 421427.Google Scholar
Royal College of Physicians (2013). Smoking and mental health. Royal College of Physicians. Retrieved from http://bit.ly/1N2GheS.Google Scholar
Scott, J. G., Matuschka, L., Niemelä, S., Miettunen, J., Emmerson, B., & Mustonen, A. (2018). Evidence of a causal relationship between smoking tobacco and schizophrenia spectrum disorders. Frontiers in Psychiatry, 9, 19.CrossRefGoogle ScholarPubMed
Taylor, A. E., Fluharty, M. E., Bjørngaard, J. H., Gabrielsen, M. E., Skorpen, F., Marioni, R. E., … Munafò, M. R. (2014). Investigating the possible causal association of smoking with depression and anxiety using Mendelian randomisation meta-analysis: the CARTA consortium. BMJ Open, 4, e006141.CrossRefGoogle ScholarPubMed
Taylor, G., McNeill, A., Girling, A., Farley, A., Lindson-Hawley, N., & Aveyard, P. (2014). Change in mental health after smoking cessation: systematic review and meta-analysis. BMJ, 348, g1151.CrossRefGoogle ScholarPubMed
Thorgeirsson, T. E., Geller, F., Sulem, P., Rafnar, T., Wiste, A., Magnusson, K. P., … Stefansson, K. (2008). A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature, 452, 638642.CrossRefGoogle ScholarPubMed
Tjora, T., Hetland, J., Aarø, L. E., Wold, B., Wiium, N., & Øverland, S. (2014). The association between smoking and depression from adolescence to adulthood. Addiction, 109, 10221030.CrossRefGoogle Scholar
Tobacco Consortium (2010). Genome-wide meta-analyses identify multiple loci associated with smoking behavior. Nature Genetics, 42, 441447.CrossRefGoogle Scholar
United Nations Office on Drugs and Crime (UNODC) (2011). World drug report 2011. Vienna: UNODC.Google Scholar
Ware, J. J., van den Bree, M. B., & Munafò, M. R. (2011). Association of the CHRNA5-A3-B4 gene cluster with heaviness of smoking: a meta-analysis. Nicotine & Tobacco Research, 13, 11671175.CrossRefGoogle ScholarPubMed
Wium-Andersen, M. K., Ørsted, D. D., & Nordestgaard, B. G. (2015). Tobacco smoking is causally associated with antipsychotic medication use and schizophrenia, but not with antidepressant medication use or depression. International Journal of Epidemiology, 44, 566577.CrossRefGoogle ScholarPubMed
World Health Organization (2011). WHO report on the global tobacco epidemic. World Health Organization. Retrieved from https://apps.who.int/iris/bitstream/handle/10665/44616/9789240687813_eng.pdf;jsessionid=875DB2CB25F1FB9A84A7581CDAB1647E?sequence=1.Google Scholar
World Health Organization (2018). Schizophrenia. World Health Organization. Retrieved from https://www.who.int/news-room/fact-sheets/detail/schizophrenia (Accessed 19 December 2018).Google Scholar
Wray, N.R., Ripke, S., Mattheisen, M., Trzaskowski, M., Byrne, E. M., Abdellaoui, A., … the Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium (2018). Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature Genetics, 50, 668681.CrossRefGoogle ScholarPubMed
Zammit, S., Allebeck, P., Dalman, C., Lundberg, I., Hemmingsson, T., & Lewis, G. (2003). Investigating the association between cigarette smoking and schizophrenia in a cohort study. American Journal of Psychiatry, 160, 22162221.CrossRefGoogle ScholarPubMed
Zeilinger, S., Kühnel, B., Klopp, N., Baurecht, H., Kleinschmidt, A., Gieger, C., … Illig, T. (2013). Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS ONE, 8, e63812.CrossRefGoogle ScholarPubMed
Zhao, Q., Wang, J., Hemani, G., Bowden, J., & Small, D. S. (2018). Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. arXiv. arXiv:1801.09652, 159.Google Scholar
Figure 0

Fig. 1. Manhattan plot of genome-wide association study of lifetime smoking index (N = 462 690). The x-axis represents chromosomal position and the y-axis represents −log10 p value for the association of each SNP with the lifetime smoking index using an additive model and linear regression. The dashed line indicates the genome-wide level of significance (p < 5 × 10−8) and genome-wide significant SNPs are indicated in red.

Figure 1

Table 1. Two-sample MR analyses of the effect of smoking exposure on schizophrenia and depression

Figure 2

Table 2. Two-sample MR analyses of the effect of schizophrenia and depression on smoking

Supplementary material: File

Wootton et al. supplementary material

Wootton et al. supplementary material

Download Wootton et al. supplementary material(File)
File 6 MB