MR-PheWAS: hypothesis prioritization among potential causal effects of body mass index on many outcomes, using Mendelian randomization

The MR-Base platform supports systematic causal inference across the human phenome

Results from genome-wide association studies (GWAS) can be used to infer causal relationships between phenotypes, using a strategy known as 2-sample Mendelian randomization (2SMR) and bypassing the need for individual-level data. However, 2SMR methods are evolving rapidly and GWAS results are often insufficiently curated, undermining efficient implementation of the approach. We therefore developed MR-Base (<ext-link ext-link-type="uri" xlink:href="http://www.mrbase.org">http://www.mrbase.org</ext-link>): a platform that integrates a curated database of complete GWAS results (no restrictions according to statistical significance) with an application programming interface, web app and R packages that automate 2SMR. The software includes several sensitivity analyses for assessing the impact of horizontal pleiotropy and other violations of assumptions. The database currently comprises 11 billion single nucleotide polymorphism-trait associations from 1673 GWAS and is updated on a regular basis. Integrating data with software ensures more rigorous application of hypothesis-driven analyses and allows millions of potential causal relationships to be efficiently evaluated in phenome-wide association studies.

The MR-Base platform supports systematic causal inference across the human phenome

Results from genome-wide association studies (GWAS) can be used to infer causal relationships between phenotypes, using a strategy known as 2-sample Mendelian randomization (2SMR) and bypassing the need for individual-level data. However, 2SMR methods are evolving rapidly and GWAS results are often insufficiently curated, undermining efficient implementation of the approach. We therefore developed MR-Base (http://www.mrbase.org): a platform that integrates a curated database of complete GWAS results (no restrictions according to statistical significance) with an application programming interface, web app and R packages that automate 2SMR. The software includes several sensitivity analyses for assessing the impact of horizontal pleiotropy and other violations of assumptions. The database currently comprises 11 billion single nucleotide polymorphism-trait associations from 1673 GWAS and is updated on a regular basis. Integrating data with software ensures more rigorous application of hypothesis-driven analyses and allows millions of potential causal relationships to be efficiently evaluated in phenome-wide association studies.

Genetics of obesity: what genetic association studies have taught us about the biology of obesity and its complications

Genome-wide association studies (GWAS) for BMI, waist-to-hip ratio, and other adiposity traits have identified more than 300 single-nucleotide polymorphisms (SNPs). Although there is reason to hope that these discoveries will eventually lead to new preventive and therapeutic agents for obesity, this will take time because such developments require detailed mechanistic understanding of how an SNP influences phenotype (and this information is largely unavailable). Fortunately, absence of functional information has not prevented GWAS findings from providing insights into the biology of obesity. Genes near loci regulating total body mass are enriched for expression in the CNS, whereas genes for fat distribution are enriched in adipose tissue itself. Gene by environment and lifestyle interaction analyses have revealed that our increasingly obesogenic environment might be amplifying genetic risk for obesity, yet those at highest risk could mitigate this risk by increasing physical activity and possibly by avoiding specific dietary components. GWAS findings have also been used in mendelian randomisation analyses probing the causal association between obesity and its many putative complications. In supporting a causal association of obesity with diabetes, coronary heart disease, specific cancers, and other conditions, these analyses have clinical relevance in identifying which outcomes could be preventable through weight loss interventions.

MR-PheWAS: exploring the causal effect of SUA level on multiple disease outcomes by using genetic instruments in UK Biobank

Objectives

We aimed to investigate the role of serum uric acid (SUA) level in a broad spectrum of disease outcomes using data for 120 091 individuals from UK Biobank.

Methods

We performed a phenome-wide association study (PheWAS) to identify disease outcomes associated with SUA genetic risk loci. We then implemented conventional Mendelianrandomisation (MR) analysis to investigate the causal relevance between SUA level and disease outcomes identified from PheWAS. We next applied MR Egger analysis to detect and account for potential pleiotropy, which conventional MR analysis might mistake for causality, and used the HEIDI (heterogeneity in dependent instruments) test to remove cross-phenotype associations that were likely due to genetic linkage.

Results

Our PheWAS identified 25 disease groups/outcomes associated with SUA genetic risk loci after multiple testing correction (P<8.57e-05). Our conventional MR analysis implicated a causal role of SUA level in three disease groups: inflammatory polyarthropathies (OR=1.22, 95% CI 1.11 to 1.34), hypertensive disease (OR=1.08, 95% CI 1.03 to 1.14) and disorders of metabolism (OR=1.07, 95% CI 1.01 to 1.14); and four disease outcomes: gout (OR=4.88, 95% CI 3.91 to 6.09), essential hypertension (OR=1.08, 95% CI 1.03 to 1.14), myocardial infarction (OR=1.16, 95% CI 1.03 to 1.30) and coeliac disease (OR=1.41, 95% CI 1.05 to 1.89). After balancing pleiotropic effects in MR Egger analysis, only gout and its encompassing disease group of inflammatory polyarthropathies were considered to be causally associated with SUA level. Our analysis highlighted a locus (ATXN2/S2HB3) that may influence SUA level and multiple cardiovascular and autoimmune diseases via pleiotropy.

Conclusions

Elevated SUA level is convincing to cause gout and inflammatory polyarthropathies, and might act as a marker for the wider range of diseases with which it associates. Our findings support further investigation on the clinical relevance of SUA level with cardiovascular, metabolic, autoimmune and respiratory diseases.

Current Scope and Challenges in Phenome-Wide Association Studies

PURPOSE OF REVIEW

Over many decades, researchers have been designing studies to investigate the relationship between genotypes and phenotypes to gain an understanding about the effect of genetics on disease. Recently, a high-throughput approach called phenome-wide associations studies (PheWAS) have been extensively used to identify associations between genetic variants and many diseases and traits simultaneously. In this review, we describe the value of PheWAS along with methodological issues and challenges in interpretation for current applications of PheWAS.

RECENT FINDINGS

PheWAS have uncovered a paradigm to identify new associations for genetic loci across many diseases. The application of PheWAS have been effective with phenotype data from electronic health records, epidemiological studies, and clinical trials data.

SUMMARY

The key strength of a PheWAS is to identify the association of one or more genetic variants with multiple phenotypes, which can showcase interconnections among the phenotypes due to shared genetic associations. While the PheWAS approach appears promising, there are a number of challenges that need to be addressed to provide additional robustness to PheWAS findings.

Recent Developments in Mendelian Randomization Studies

PURPOSE OF REVIEW

Mendelian randomization (MR) is a strategy for evaluating causality in observational epidemiological studies. MR exploits the fact that genotypes are not generally susceptible to reverse causation and confounding, due to their fixed nature and Mendel's First and Second Laws of Inheritance. MR has the potential to provide information on causality in many situations where randomized controlled trials are not possible, but the results of MR studies must be interpreted carefully to avoid drawing erroneous conclusions.

RECENT FINDINGS

In this review, we outline the principles behind MR, as well as assumptions and limitations of the method. Extensions to the basic approach are discussed, including two-sample MR, bidirectional MR, two-step MR, multivariable MR, and factorial MR. We also consider some new applications and recent developments in the methodology, including its ability to inform drug development, automation of the method using tools such as MR-Base, and phenome-wide and hypothesis-free MR.

SUMMARY

In conjunction with the growing availability of large-scale genomic databases, higher level of automation and increased robustness of the methods, MR promises to be a valuable strategy to examine causality in complex biological/omics networks, inform drug development and prioritize intervention targets for disease prevention in the future.

The role of glycaemic and lipid risk factors in mediating the effect of BMI on coronary heart disease: a two-step, two-sample Mendelian randomisation study

AIMS/HYPOTHESIS

The extent to which effects of BMI on CHD are mediated by glycaemic and lipid risk factors is unclear. In this study we examined the effects of these traits using genetic evidence.

METHODS

We used two-sample Mendelian randomisation to determine: (1) the causal effect of BMI on CHD (60,801 case vs 123,504 control participants), type 2 diabetes (34,840 case vs 114,981 control participants), fasting glucose (n = 46,186), insulin (n = 38,238), HbA_1c (n = 46,368) and LDL-cholesterol, HDL-cholesterol and triacylglycerols (n = 188,577); (2) the causal effects of glycaemic and lipids traits on CHD; and (3) the extent to which these traits mediate any effect of BMI on CHD.

RESULTS

One SD higher BMI (~ 4.5 kg/m²) was associated with higher risk of CHD (OR 1.45 [95% CI 1.27, 1.66]) and type 2 diabetes (1.96 [95% CI 1.35, 2.83]), higher levels of fasting glucose (0.07 mmol/l [95% CI 0.03, 0.11]), HbA_1c (0.05% [95% CI 0.01, 0.08]), fasting insulin (0.18 log pmol/l [95% CI 0.14, 0.22]) and triacylglycerols (0.20 SD [95% CI 0.14, 0.26]) and lower levels of HDL-cholesterol (-0.23 SD [95% CI -0.32, -0.15]). There was no evidence for a causal relation between BMI and LDL-cholesterol. The causal associations of higher triacylglycerols, HbA_1c and diabetes risk with CHD risk remained after performing sensitivity analyses that considered different models of horizontal pleiotropy. The BMI-CHD effect reduced from 1.45 to 1.16 (95% CI 0.99, 1.36) and to 1.36 (95% CI 1.19, 1.57) with genetic adjustment for triacylglycerols or HbA_1c, respectively, and to 1.09 (95% CI 0.94, 1.27) with adjustment for both.

CONCLUSIONS/INTERPRETATION

Increased triacylglycerol levels and poor glycaemic control appear to mediate much of the effect of BMI on CHD.

Software Application Profile: PHESANT: a tool for performing automated phenome scans in UK Biobank

Motivation

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.

General features

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.

Implementation

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.

Availability

PHESANT is available on GitHub at [https://github.com/MRCIEU/PHESANT]. Git tag v0.5 corresponds to the version presented here.

Integrative variants, haplotypes and diplotypes of the CAPN3 and FRMD5 genes and several environmental exposures associate with serum lipid variables

To determine whether the integrative variants, haplotypes and diplotypes of the calpain 3 (CAPN3) and the FERM domain containing 5 genes (FRMD5) and several environmental exposures are associated with an implication in lipid homeostasis, which are associated with cardiovascular risk. Genotyping of the CAPN3 rs4344713 and FRMD5 rs524908 was performed by Sanger sequencing in 1,640 subjects (Jing, 819 and Han, 821). Multivariate analyses of covariance models that adjusted by age, gender, body mass index (BMI), blood pressure and lifestyle (smoking and drinking), were constructed using variants, haplotypes and diplotypes of the CAPN3 rs4344713 and FRMD5 rs524908 as predictors and changes in lipid variables. Significant associations with low-density lipoprotein cholesterol and apolipoprotein (Apo) B were found. Linkage disequilibrium with each other showed the haplotype-phenotype associations with triglyceride and ApoA1. This study also suggested pleiotropic associations of the CAPN3-FRMD5 diplotypes with lipid variables. As potential confounders, diastolic blood pressure (DBP) and BMI were significantly associated with lipid variables. We conclude that integrative variants, haplotypes and diplotypes of the CAPN3 rs4344713 and FRMD5 rs524908, as well as DBP and BMI are associated with serum lipid variables in the Jing and Han populations.

Genetic Obesity and the Risk of Atrial Fibrillation: Causal Estimates from Mendelian Randomization

BACKGROUND

Observational studies have identified an association between body mass index (BMI) and incident atrial fibrillation (AF). Inferring causality from observational studies, however, is subject to residual confounding, reverse causation, and bias. The primary objective of this study was to evaluate the causal association between BMI and AF by using genetic predictors of BMI.

METHODS

We identified 51 646 individuals of European ancestry without AF at baseline from 7 prospective population-based cohorts initiated between 1987 and 2002 in the United States, Iceland, and the Netherlands with incident AF ascertained between 1987 and 2012. Cohort-specific mean follow-up ranged from 7.4 to 19.2 years, over which period there was a total of 4178 cases of incident AF. We performed a Mendelian randomization with instrumental variable analysis to estimate a cohort-specific causal hazard ratio for the association between BMI and AF. Two genetic instruments for BMI were used: FTO genotype (rs1558902) and a BMI gene score comprising 39 single-nucleotide polymorphisms identified by genome-wide association studies to be associated with BMI. Cohort-specific estimates were combined by random-effects, inverse variance-weighted meta-analysis.

RESULTS

In age- and sex-adjusted meta-analysis, both genetic instruments were significantly associated with BMI (FTO: 0.43 [95% confidence interval, 0.32-0.54] kg/m2 per A-allele, P<0.001; BMI gene score: 1.05 [95% confidence interval, 0.90-1.20] kg/m2 per 1-U increase, P<0.001) and incident AF (FTO, hazard ratio, 1.07 [1.02-1.11] per A-allele, P=0.004; BMI gene score, hazard ratio, 1.11 [1.05-1.18] per 1-U increase, P<0.001). Age- and sex-adjusted instrumental variable estimates for the causal association between BMI and incident AF were hazard ratio, 1.15 (1.04-1.26) per kg/m2, P=0.005 (FTO) and 1.11 (1.05-1.17) per kg/m2, P<0.001 (BMI gene score). Both of these estimates were consistent with the meta-analyzed estimate between observed BMI and AF (age- and sex-adjusted hazard ratio 1.05 [1.04-1.06] per kg/m2, P<0.001). Multivariable adjustment did not significantly change findings.

CONCLUSIONS

Our data are consistent with a causal relationship between BMI and incident AF. These data support the possibility that public health initiatives targeting primordial prevention of obesity may reduce the incidence of AF.

Phenome-Wide Association Studies as a Tool to Advance Precision Medicine

Beginning in the early 2000s, the accumulation of biospecimens linked to electronic health records (EHRs) made possible genome-phenome studies (i.e., comparative analyses of genetic variants and phenotypes) using only data collected as a by-product of typical health care. In addition to disease and trait genetics, EHRs proved a valuable resource for analyzing pharmacogenetic traits and developing reverse genetics approaches such as phenome-wide association studies (PheWASs). PheWASs are designed to survey which of many phenotypes may be associated with a given genetic variant. PheWAS methods have been validated through replication of hundreds of known genotype-phenotype associations, and their use has differentiated between true pleiotropy and clinical comorbidity, added context to genetic discoveries, and helped define disease subtypes, and may also help repurpose medications. PheWAS methods have also proven to be useful with research-collected data. Future efforts that integrate broad, robust collection of phenotype data (e.g., EHR data) with purpose-collected research data in combination with a greater understanding of EHR data will create a rich resource for increasingly more efficient and detailed genome-phenome analysis to usher in new discoveries in precision medicine.

Mendelian Randomization as an Approach to Assess Causality Using Observational Data

Mendelian randomization refers to an analytic approach to assess the causality of an observed association between a modifiable exposure or risk factor and a clinically relevant outcome. It presents a valuable tool, especially when randomized controlled trials to examine causality are not feasible and observational studies provide biased associations because of confounding or reverse causality. These issues are addressed by using genetic variants as instrumental variables for the tested exposure: the alleles of this exposure-associated genetic variant are randomly allocated and not subject to reverse causation. This, together with the wide availability of published genetic associations to screen for suitable genetic instrumental variables make Mendelian randomization a time- and cost-efficient approach and contribute to its increasing popularity for assessing and screening for potentially causal associations. An observed association between the genetic instrumental variable and the outcome supports the hypothesis that the exposure in question is causally related to the outcome. This review provides an overview of the Mendelian randomization method, addresses assumptions and implications, and includes illustrative examples. We also discuss special issues in nephrology, such as inverse risk factor associations in advanced disease, and outline opportunities to design Mendelian randomization studies around kidney function and disease.

Body mass index: Has epidemiology started to break down causal contributions to health and disease?

OBJECTIVES

To review progress in understanding the methods and results concerning the causal contribution of body mass index (BMI) to health and disease.

METHODS

In the context of conventional evidence focused on the relationship between BMI and health, this review considers current literature on the common, population-based, genetic contribution to BMI and how this has fed into the developing field of applied epidemiology.

RESULTS

Technological and analytical developments have driven considerable success in identifying genetic variants relevant to BMI. This has enabled the implementation of Mendelian randomization to address questions of causality. The product of this work has been the implication of BMI as a causal agent in a host of health outcomes. Further breakdown of causal pathways by integration with other "omics" technologies promises to deliver additional benefit.

CONCLUSIONS

Gaps remain in our understanding of BMI as a risk factor for health and disease, and while promising, applied genetic epidemiology should be considered alongside alternative methods for assessing the impact of BMI on health. Potential limitations, relating to inappropriate or nonspecific measures of obesity and the improper use of genetic instruments, will need to be explored and incorporated into future research aiming to dissect BMI as a risk factor.

Unravelling the human genome-phenome relationship using phenome-wide association studies

Advances in genotyping technology have, over the past decade, enabled the focused search for common genetic variation associated with human diseases and traits. With the recently increased availability of detailed phenotypic data from electronic health records and epidemiological studies, the impact of one or more genetic variants on the phenome is starting to be characterized both in clinical and population-based settings using phenome-wide association studies (PheWAS). These studies reveal a number of challenges that will need to be overcome to unlock the full potential of PheWAS for the characterization of the complex human genome-phenome relationship.