Mendelian randomization (MR) has been used to estimate the causal effect of body mass index (BMI) on particular traits thought to be affected by BMI. However, BMI may also be a modifiable, causal ...risk factor for outcomes where there is no prior reason to suggest that a causal effect exists. We performed a MR phenome-wide association study (MR-pheWAS) to search for the causal effects of BMI in UK Biobank (n = 334 968), using the PHESANT open-source phenome scan tool. A subset of identified associations were followed up with a formal two-stage instrumental variable analysis in UK Biobank, to estimate the causal effect of BMI on these phenotypes. Of the 22 922 tests performed, our MR-pheWAS identified 587 associations below a stringent P value threshold corresponding to a 5% estimated false discovery rate. These included many previously identified causal effects, for instance, an adverse effect of higher BMI on risk of diabetes and hypertension. We also identified several novel effects, including protective effects of higher BMI on a set of psychosocial traits, identified initially in our preliminary MR-pheWAS in circa 115,000 UK Biobank participants and replicated in a different subset of circa 223,000 UK Biobank participants. Our comprehensive MR-pheWAS identified potential causal effects of BMI on a large and diverse set of phenotypes. This included both previously identified causal effects, and novel effects such as a protective effect of higher BMI on feelings of nervousness.
Spending more time active (and less sedentary) is associated with health benefits such as improved cardiovascular health and lower risk of all-cause mortality. It is unclear whether these ...associations differ depending on whether time spent sedentary or in moderate-vigorous physical activity (MVPA) is accumulated in long or short bouts. In this study, we used a novel method that accounts for substitution (i.e., more time in MVPA means less time sleeping, in light activity or sedentary) to examine whether length of sedentary and MVPA bouts associates with all-cause mortality.
We used data on 79,503 adult participants from the population-based UK Biobank cohort, which recruited participants between 2006 and 2010 (mean age at accelerometer wear 62.1 years SD = 7.9, 54.5% women; mean length of follow-up 5.1 years SD = 0.73). We derived (1) the total time participants spent in activity categories-sleep, sedentary, light activity, and MVPA-on average per day; (2) time spent in sedentary bouts of short (1 to 15 minutes), medium (16 to 40 minutes), and long (41+ minutes) duration; and (3) MVPA bouts of very short (1 to 9 minutes), short (10 to 15 minutes), medium (16 to 40 minutes), and long (41+ minutes) duration. We used Cox proportion hazards regression to estimate the association of spending 10 minutes more average daily time in one activity or bout length category, coupled with 10 minutes less time in another, with all-cause mortality. Those spending more time in MVPA had lower mortality risk, irrespective of whether this replaced time spent sleeping, sedentary, or in light activity, and these associations were of similar magnitude (e.g., hazard ratio HR 0.96 95% CI: 0.94, 0.97; P < 0.001 per 10 minutes more MVPA, coupled with 10 minutes less light activity per day). Those spending more time sedentary had higher mortality risk if this replaced light activity (HR 1.02 95% CI: 1.01, 1.02; P < 0.001 per 10 minutes more sedentary time, with 10 minutes less light activity per day) and an even higher risk if this replaced MVPA (HR 1.06 95% CI: 1.05, 1.08; P < 0.001 per 10 minutes more sedentary time, with 10 minutes less MVPA per day). We found little evidence that mortality risk differed depending on the length of sedentary or MVPA bouts. Key limitations of our study are potential residual confounding, the limited length of follow-up, and use of a select sample of the United Kingdom population.
We have shown that time spent in MVPA was associated with lower mortality, irrespective of whether it replaced time spent sleeping, sedentary, or in light activity. Time spent sedentary was associated with higher mortality risk, particularly if it replaced MVPA. This emphasises the specific importance of MVPA. Our findings suggest that the impact of MVPA does not differ depending on whether it is obtained from several short bouts or fewer longer bouts, supporting the recent removal of the requirement that MVPA should be accumulated in bouts of 10 minutes or more from the UK and the United States policy. Further studies are needed to investigate causality and explore health outcomes beyond mortality.
Psychiatric disorders are highly heritable and associated with a wide variety of social adversity and physical health problems. Using genetic liability (rather than phenotypic measures of disease) as ...a proxy for psychiatric disease risk can be a useful alternative for research questions that would traditionally require large cohort studies with long-term follow up. Here we conducted a hypothesis-free phenome-wide association study in about 330,000 participants from the UK Biobank to examine associations of polygenic risk scores (PRS) for five psychiatric disorders (major depression (MDD), bipolar disorder (BP), schizophrenia (SCZ), attention-deficit/ hyperactivity disorder (ADHD) and autism spectrum disorder (ASD)) with 23,004 outcomes in UK Biobank, using the open-source PHESANT software package. There was evidence after multiple testing (p<2.55x10-06) for associations of PRSs with 294 outcomes, most of them attributed to associations of PRSMDD (n = 167) and PRSSCZ (n = 157) with mental health factors. Among others, we found strong evidence of association of higher PRSADHD with 1.1 months younger age at first sexual intercourse 95% confidence interval CI: -1.25,-0.92 and a history of physical maltreatment; PRSASD with 0.01% lower erythrocyte distribution width 95%CI: -0.013,-0.007; PRSSCZ with 0.95 lower odds of playing computer games 95%CI:0.95,0.96; PRSMDD with a 0.12 points higher neuroticism score 95%CI:0.111,0.135 and PRSBP with 1.03 higher odds of having a university degree 95%CI:1.02,1.03. We were able to show that genetic liabilities for five major psychiatric disorders associate with long-term aspects of adult life, including socio-demographic factors, mental and physical health. This is evident even in individuals from the general population who do not necessarily present with a psychiatric disorder diagnosis.
Mendelian randomization (MR) is an established approach to evaluate the effect of an exposure on an outcome. The gene-by-environment (GxE) study design can be used to determine whether the genetic ...instrument affects the outcome through pathways other than via the exposure of interest (horizontal pleiotropy). MR phenome-wide association studies (MR-pheWAS) search for the effects of an exposure, and can be conducted in UK Biobank using the PHESANT package. In this proof-of-principle study, we introduce the novel GxE MR-pheWAS approach, that combines MR-pheWAS with the use of GxE interactions. This method aims to identify the presence of effects of an exposure while simultaneously investigating horizontal pleiotropy. We systematically test for the presence of 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. We used PHESANT to test for the presence of effects of smoking heaviness, instrumented by genetic variant rs16969968, among never and ever smokers respectively, in UK Biobank. We ranked results by the strength of interaction between ever and never smokers. We replicated previously established effects of smoking heaviness, including detrimental effects on lung function. Novel results included a detrimental effect of heavier smoking on facial aging. We have demonstrated how GxE MR-pheWAS can be used to identify potential effects of an exposure, while simultaneously assessing whether results may be biased by horizontal pleiotropy.
Epidemiological cohorts typically contain a diverse set of phenotypes such that automation of phenome scans is non-trivial, because they require highly heterogeneous models. For this reason, phenome ...scans have to date tended to use a smaller homogeneous set of phenotypes that can be analysed in a consistent fashion. We present PHESANT (PHEnome Scan ANalysis Tool), a software package for performing comprehensive phenome scans in UK Biobank.
PHESANT tests the association of a specified trait with all continuous, integer and categorical variables in UK Biobank, or a specified subset. PHESANT uses a novel rule-based algorithm to determine how to appropriately test each trait, then performs the analyses and produces plots and summary tables.
The PHESANT phenome scan is implemented in R. PHESANT includes a novel Javascript D3.js visualization and accompanying Java code that converts the phenome scan results to the required JavaScript Object Notation (JSON) format.
PHESANT is available on GitHub at https://github.com/MRCIEU/PHESANT. Git tag v0.5 corresponds to the version presented here.
Insomnia symptoms are widespread in the population and might have effects on many chronic conditions and their risk factors but previous research has focused on select hypothesised ...associations/effects rather than taking a systematic hypothesis-free approach across many health outcomes.
We performed a Mendelian randomisation (MR) phenome-wide association study (PheWAS) in 336,975 unrelated white-British UK Biobank participants. Self-reported insomnia symptoms were instrumented by a genetic risk score (GRS) created from 129 single-nucleotide polymorphisms (SNPs). A total of 11,409 outcomes from UK Biobank were extracted and processed by an automated pipeline (PHESANT) for the MR-PheWAS. Potential causal effects (those passing a Bonferroni-corrected significance threshold) were followed up with two-sample MR in MR-Base, where possible.
Four hundred thirty-seven potential causal effects of insomnia symptoms were observed for a diverse range of outcomes, including anxiety, depression, pain, body composition, respiratory, musculoskeletal and cardiovascular traits. We were able to undertake two-sample MR for 71 of these 437 and found evidence of causal effects (with directionally concordant effect estimates across main and sensitivity analyses) for 30 of these. These included novel findings (by which we mean not extensively explored in conventional observational studies and not previously explored using MR based on a systematic search) of an adverse effect on risk of spondylosis (OR 95%CI = 1.55 1.33, 1.81) and bronchitis (OR 95%CI = 1.12 1.03, 1.22), among others.
Insomnia symptoms potentially cause a wide range of adverse health-related outcomes and behaviours. This has implications for developing interventions to prevent and treat a number of diseases in order to reduce multimorbidity and associated polypharmacy.
Alzheimer's disease (AD) has no proven causal and modifiable risk factors, or effective interventions. We report a phenome-wide association study (PheWAS) of genetic liability for AD in 334,968 ...participants of the UK Biobank study, stratified by age. We also examined the effects of AD genetic liability on previously implicated risk factors. We replicated these analyses in the HUNT study. PheWAS hits and previously implicated risk factors were followed up in a Mendelian randomization (MR) framework to identify the causal effect of each risk factor on AD risk. A higher genetic liability for AD was associated with medical history and cognitive, lifestyle, physical and blood-based measures as early as 39 years of age. These effects were largely driven by the APOE gene. The follow-up MR analyses were primarily null, implying that most of these associations are likely to be a consequence of prodromal disease or selection bias, rather than the risk factor causing the disease.
Age at menarche has been associated with various health outcomes. We aimed to identify potential causal effects of age at menarche on health-related traits in a hypothesis-free manner.
We conducted a ...Mendelian randomization phenome-wide association study (MR-pheWAS) of age at menarche with 17,893 health-related traits in UK Biobank (n = 181,318) using PHESANT. The exposure of interest was the genetic risk score for age at menarche. We conducted a second MR-pheWAS after excluding SNPs associated with BMI from the genetic risk score, to examine whether results might be due to the genetic overlap between age at menarche and BMI. We followed up a subset of health-related traits to investigate MR assumptions and seek replication in independent study populations.
Of the 17,893 tests performed in our MR-pheWAS, we identified 619 associations with the genetic risk score for age at menarche at a 5% false discovery rate threshold, of which 295 were below a Bonferroni-corrected P value threshold. These included potential effects of younger age at menarche on lower lung function, higher heel bone-mineral density, greater burden of psychosocial/mental health problems, younger age at first birth, higher risk of childhood sexual abuse, poorer cardiometabolic health, and lower physical activity. After exclusion of variants associated with BMI, the genetic risk score for age at menarche was related to 37 traits at a 5% false discovery rate, of which 29 were below a Bonferroni-corrected P value threshold. We attempted to replicate findings for bone-mineral density, lung function, neuroticism, and childhood sexual abuse using 5 independent cohorts/consortia. While estimates for lung function, higher bone-mineral density, neuroticism, and childhood sexual abuse in replication cohorts were consistent with UK Biobank estimates, confidence intervals were wide and often included the null.
The genetic risk score for age at menarche was related to a broad range of health-related traits. Follow-up analyses indicated imprecise evidence of an effect of younger age at menarche on greater bone-mineral density, lower lung function, higher neuroticism score, and greater risk of childhood sexual abuse in the smaller replication samples available; hence, these findings need further exploration when larger independent samples become available.
Abstract
Background
A lack of genetic data across generations makes transgenerational Mendelian randomization (MR) difficult. We used UK Biobank and a novel proxy gene-by-environment MR to ...investigate effects of maternal smoking heaviness in pregnancy on offspring health, using participants’ (generation one: G1) genotype (rs16969968 in CHRNA5) as a proxy for their mothers’ (G0) genotype.
Methods
We validated this approach by replicating an established effect of maternal smoking heaviness on offspring birthweight. Then we applied this approach to explore effects of maternal (G0) smoking heaviness on offspring (G1) later life outcomes and on birthweight of G1 women’s children (G2).
Results
Each additional smoking-increasing allele in offspring (G1) was associated with a 0.018 95% confidence interval (CI): -0.026, -0.009 kg lower G1 birthweight in maternal (G0) smoking stratum, but no meaningful effect (-0.002 kg; 95% CI: -0.008, 0.003) in maternal non-smoking stratum (interaction P-value = 0.004). The differences in associations of rs16969968 with grandchild’s (G2) birthweight between grandmothers (G0) who did, versus did not, smoke were heterogeneous (interaction P-value = 0.042) among mothers (G1) who did (-0.020 kg/allele; 95% CI: -0.044, 0.003), versus did not (0.007 kg/allele; 95% CI: -0.005, 0.020), smoke in pregnancy.
Conclusions
Our study demonstrated how offspring genotype can be used to proxy for the mother’s genotype in gene-by-environment MR. We confirmed the causal effect of maternal (G0) smoking on offspring (G1) birthweight, but found little evidence of an effect on G1 longer-term health outcomes. For grandchild’s (G2) birthweight, the effect of grandmother’s (G0) smoking heaviness in pregnancy may be modulated by maternal (G1) smoking status in pregnancy.
Mendelian randomization may give biased causal estimates if the instrument affects the outcome not solely via the exposure of interest (violating the exclusion restriction assumption). We demonstrate ...use of a global randomization test as a falsification test for the exclusion restriction assumption. Using simulations, we explored the statistical power of the randomization test to detect an association between a genetic instrument and a covariate set due to (a) selection bias or (b) horizontal pleiotropy, compared to three approaches examining associations with individual covariates: (i) Bonferroni correction for the number of covariates, (ii) correction for the effective number of independent covariates, and (iii) an r
permutation-based approach. We conducted proof-of-principle analyses in UK Biobank, using CRP as the exposure and coronary heart disease (CHD) as the outcome. In simulations, power of the randomization test was higher than the other approaches for detecting selection bias when the correlation between the covariates was low (r
< 0.1), and at least as powerful as the other approaches across all simulated horizontal pleiotropy scenarios. In our applied example, we found strong evidence of selection bias using all approaches (e.g., global randomization test p < 0.002). We identified 51 of the 58 CRP genetic variants as horizontally pleiotropic, and estimated effects of CRP on CHD attenuated somewhat to the null when excluding these from the genetic risk score (OR = 0.96 95% CI: 0.92, 1.00 versus 0.97 95% CI: 0.90, 1.05 per 1-unit higher log CRP levels). The global randomization test can be a useful addition to the MR researcher's toolkit.