Identifying drug targets for neurological and psychiatric disease via genetics and the brain transcriptome.

Discovering drugs that efficiently treat brain diseases has been challenging. Genetic variants that modulate the expression of potential drug targets can be utilized to assess the efficacy of therapeutic interventions. We therefore employed Mendelian Randomization (MR) on gene expression measured in brain tissue to identify drug targets involved in neurological and psychiatric diseases. We conducted a two-sample MR using cis-acting brain-derived expression quantitative trait loci (eQTLs) from the Accelerating Medicines Partnership for Alzheimer's Disease consortium (AMP-AD) and the CommonMind Consortium (CMC) meta-analysis study (n = 1,286) as genetic instruments to predict the effects of 7,137 genes on 12 neurological and psychiatric disorders. We conducted Bayesian colocalization analysis on the top MR findings (using P<6x10-7 as evidence threshold, Bonferroni-corrected for 80,557 MR tests) to confirm sharing of the same causal variants between gene expression and trait in each genomic region. We then intersected the colocalized genes with known monogenic disease genes recorded in Online Mendelian Inheritance in Man (OMIM) and with genes annotated as drug targets in the Open Targets platform to identify promising drug targets. 80 eQTLs showed MR evidence of a causal effect, from which we prioritised 47 genes based on colocalization with the trait. We causally linked the expression of 23 genes with schizophrenia and a single gene each with anorexia, bipolar disorder and major depressive disorder within the psychiatric diseases and 9 genes with Alzheimer's disease, 6 genes with Parkinson's disease, 4 genes with multiple sclerosis and two genes with amyotrophic lateral sclerosis within the neurological diseases we tested. From these we identified five genes (ACE, GPNMB, KCNQ5, RERE and SUOX) as attractive drug targets that may warrant follow-up in functional studies and clinical trials, demonstrating the value of this study design for discovering drug targets in neuropsychiatric diseases.

Role of circulating polyunsaturated fatty acids on cardiovascular diseases risk: analysis using Mendelian randomization and fatty acid genetic association data from over 114,000 UK Biobank participants.

BACKGROUND: Despite early interest in the health effects of polyunsaturated fatty acids (PUFA), there is still substantial controversy and uncertainty on the evidence linking PUFA to cardiovascular diseases (CVDs). We investigated the effect of plasma concentration of omega-3 PUFA (i.e. docosahexaenoic acid (DHA) and total omega-3 PUFA) and omega-6 PUFA (i.e. linoleic acid and total omega-6 PUFA) on the risk of CVDs using Mendelian randomization. METHODS: We conducted the largest genome-wide association study (GWAS) of circulating PUFA to date including a sample of 114,999 individuals and incorporated these data in a two-sample Mendelian randomization framework to investigate the involvement of circulating PUFA on a wide range of CVDs in up to 1,153,768 individuals of European ancestry (i.e. coronary artery disease, ischemic stroke, haemorrhagic stroke, heart failure, atrial fibrillation, peripheral arterial disease, aortic aneurysm, venous thromboembolism and aortic valve stenosis). RESULTS: GWAS identified between 46 and 64 SNPs for the four PUFA traits, explaining 4.8-7.9% of circulating PUFA variance and with mean F statistics >100. Higher genetically predicted DHA (and total omega-3 fatty acids) concentration was related to higher risk of some cardiovascular endpoints; however, these findings did not pass our criteria for multiple testing correction and were attenuated when accounting for LDL-cholesterol through multivariable Mendelian randomization or excluding SNPs in the vicinity of the FADS locus. Estimates for the relation between higher genetically predicted linoleic acid (and total omega-6) concentration were inconsistent across different cardiovascular endpoints and Mendelian randomization methods. There was weak evidence of higher genetically predicted linoleic acid being related to lower risk of ischemic stroke and peripheral artery disease when accounting by LDL-cholesterol. CONCLUSIONS: We have conducted the largest GWAS of circulating PUFA to date and the most comprehensive Mendelian randomization analyses. Overall, our Mendelian randomization findings do not support a protective role of circulating PUFA concentration on the risk of CVDs. However, horizontal pleiotropy via lipoprotein-related traits could be a key source of bias in our analyses.

Circulating inflammatory cytokines and risk of five cancers: a Mendelian randomization analysis.

BACKGROUND: Epidemiological and experimental evidence has linked chronic inflammation to cancer aetiology. It is unclear whether associations for specific inflammatory biomarkers are causal or due to bias. In order to examine whether altered genetically predicted concentration of circulating cytokines are associated with cancer development, we performed a two-sample Mendelian randomisation (MR) analysis. METHODS: Up to 31,112 individuals of European descent were included in genome-wide association study (GWAS) meta-analyses of 47 circulating cytokines. Single nucleotide polymorphisms (SNPs) robustly associated with the cytokines, located in or close to their coding gene (cis), were used as instrumental variables. Inverse-variance weighted MR was used as the primary analysis, and the MR assumptions were evaluated in sensitivity and colocalization analyses and a false discovery rate (FDR) correction for multiple comparisons was applied. Corresponding germline GWAS summary data for five cancer outcomes (breast, endometrial, lung, ovarian, and prostate), and their subtypes were selected from the largest cancer-specific GWASs available (cases ranging from 12,906 for endometrial to 133,384 for breast cancer). RESULTS: There was evidence of inverse associations of macrophage migration inhibitory factor with breast cancer (OR per SD = 0.88, 95% CI 0.83 to 0.94), interleukin-1 receptor antagonist with endometrial cancer (0.86, 0.80 to 0.93), interleukin-18 with lung cancer (0.87, 0.81 to 0.93), and beta-chemokine-RANTES with ovarian cancer (0.70, 0.57 to 0.85) and positive associations of monokine induced by gamma interferon with endometrial cancer (3.73, 1.86 to 7.47) and cutaneous T-cell attracting chemokine with lung cancer (1.51, 1.22 to 1.87). These associations were similar in sensitivity analyses and supported in colocalization analyses. CONCLUSIONS: Our study adds to current knowledge on the role of specific inflammatory biomarker pathways in cancer aetiology. Further validation is needed to assess the potential of these cytokines as pharmacological or lifestyle targets for cancer prevention.

Causal relationships between body mass index, smoking and lung cancer: Univariable and multivariable Mendelian randomization.

At the time of cancer diagnosis, body mass index (BMI) is inversely correlated with lung cancer risk, which may reflect reverse causality and confounding due to smoking behavior. We used two-sample univariable and multivariable Mendelian randomization (MR) to estimate causal relationships of BMI and smoking behaviors on lung cancer and histological subtypes based on an aggregated genome-wide association studies (GWASs) analysis of lung cancer in 29 266 cases and 56 450 controls. We observed a positive causal effect for high BMI on occurrence of small-cell lung cancer (odds ratio (OR) = 1.60, 95% confidence interval (CI) = 1.24-2.06, P = 2.70 × 10-4 ). After adjustment of smoking behaviors using multivariable Mendelian randomization (MVMR), a direct causal effect on small cell lung cancer (ORMVMR = 1.28, 95% CI = 1.06-1.55, PMVMR = .011), and an inverse effect on lung adenocarcinoma (ORMVMR = 0.86, 95% CI = 0.77-0.96, PMVMR = .008) were observed. A weak increased risk of lung squamous cell carcinoma was observed for higher BMI in univariable Mendelian randomization (UVMR) analysis (ORUVMR = 1.19, 95% CI = 1.01-1.40, PUVMR = .036), but this effect disappeared after adjustment of smoking (ORMVMR = 1.02, 95% CI = 0.90-1.16, PMVMR = .746). These results highlight the histology-specific impact of BMI on lung carcinogenesis and imply mediator role of smoking behaviors in the association between BMI and lung cancer.

Systematic Mendelian randomization framework elucidates hundreds of CpG sites which may mediate the influence of genetic variants on disease.

We have undertaken a systematic Mendelian randomization (MR) study using methylation quantitative trait loci (meQTL) as genetic instruments to assess the relationship between genetic variation, DNA methylation and 139 complex traits. Using two-sample MR, we identified 1148 associations across 61 traits where genetic variants were associated with both proximal DNA methylation (i.e. cis-meQTL) and complex trait variation (P < 1.39 × 10-08). Joint likelihood mapping provided evidence that the genetic variant which influenced DNA methylation levels for 348 of these associations across 47 traits was also responsible for variation in complex traits. These associations showed a high rate of replication in the BIOS QTL and UK Biobank datasets for 14 selected traits, as 101 of the attempted 128 associations survived multiple testing corrections (P < 3.91 × 10-04). Integrating expression quantitative trait loci (eQTL) data suggested that genetic variants responsible for 306 of the 348 refined meQTL associations also influence gene expression, which indicates a coordinated system of effects that are consistent with causality. CpG sites were enriched for histone mark peaks in tissue types relevant to their associated trait and implicated genes were enriched across relevant biological pathways. Though we are unable to distinguish mediation from horizontal pleiotropy in these analyses, our findings should prove valuable in prioritizing candidate loci where DNA methylation may influence traits and help develop mechanistic insight into the aetiology of complex disease.

The influence of obesity-related factors in the etiology of renal cell carcinoma-A mendelian randomization study.

BACKGROUND: Several obesity-related factors have been associated with renal cell carcinoma (RCC), but it is unclear which individual factors directly influence risk. We addressed this question using genetic markers as proxies for putative risk factors and evaluated their relation to RCC risk in a mendelian randomization (MR) framework. This methodology limits bias due to confounding and is not affected by reverse causation. METHODS AND FINDINGS: Genetic markers associated with obesity measures, blood pressure, lipids, type 2 diabetes, insulin, and glucose were initially identified as instrumental variables, and their association with RCC risk was subsequently evaluated in a genome-wide association study (GWAS) of 10,784 RCC patients and 20,406 control participants in a 2-sample MR framework. The effect on RCC risk was estimated by calculating odds ratios (ORSD) for a standard deviation (SD) increment in each risk factor. The MR analysis indicated that higher body mass index increases the risk of RCC (ORSD: 1.56, 95% confidence interval [CI] 1.44-1.70), with comparable results for waist-to-hip ratio (ORSD: 1.63, 95% CI 1.40-1.90) and body fat percentage (ORSD: 1.66, 95% CI 1.44-1.90). This analysis further indicated that higher fasting insulin (ORSD: 1.82, 95% CI 1.30-2.55) and diastolic blood pressure (DBP; ORSD: 1.28, 95% CI 1.11-1.47), but not systolic blood pressure (ORSD: 0.98, 95% CI 0.84-1.14), increase the risk for RCC. No association with RCC risk was seen for lipids, overall type 2 diabetes, or fasting glucose. CONCLUSIONS: This study provides novel evidence for an etiological role of insulin in RCC, as well as confirmatory evidence that obesity and DBP influence RCC risk.

LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis.

MOTIVATION: LD score regression is a reliable and efficient method of using genome-wide association study (GWAS) summary-level results data to estimate the SNP heritability of complex traits and diseases, partition this heritability into functional categories, and estimate the genetic correlation between different phenotypes. Because the method relies on summary level results data, LD score regression is computationally tractable even for very large sample sizes. However, publicly available GWAS summary-level data are typically stored in different databases and have different formats, making it difficult to apply LD score regression to estimate genetic correlations across many different traits simultaneously. RESULTS: In this manuscript, we describe LD Hub - a centralized database of summary-level GWAS results for 173 diseases/traits from different publicly available resources/consortia and a web interface that automates the LD score regression analysis pipeline. To demonstrate functionality and validate our software, we replicated previously reported LD score regression analyses of 49 traits/diseases using LD Hub; and estimated SNP heritability and the genetic correlation across the different phenotypes. We also present new results obtained by uploading a recent atopic dermatitis GWAS meta-analysis to examine the genetic correlation between the condition and other potentially related traits. In response to the growing availability of publicly accessible GWAS summary-level results data, our database and the accompanying web interface will ensure maximal uptake of the LD score regression methodology, provide a useful database for the public dissemination of GWAS results, and provide a method for easily screening hundreds of traits for overlapping genetic aetiologies. AVAILABILITY AND IMPLEMENTATION: The web interface and instructions for using LD Hub are available at http://ldsc.broadinstitute.org/ CONTACT: jie.zheng@bristol.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator.

Developments in genome-wide association studies and the increasing availability of summary genetic association data have made application of Mendelian randomization relatively straightforward. However, obtaining reliable results from a Mendelian randomization investigation remains problematic, as the conventional inverse-variance weighted method only gives consistent estimates if all of the genetic variants in the analysis are valid instrumental variables. We present a novel weighted median estimator for combining data on multiple genetic variants into a single causal estimate. This estimator is consistent even when up to 50% of the information comes from invalid instrumental variables. In a simulation analysis, it is shown to have better finite-sample Type 1 error rates than the inverse-variance weighted method, and is complementary to the recently proposed MR-Egger (Mendelian randomization-Egger) regression method. In analyses of the causal effects of low-density lipoprotein cholesterol and high-density lipoprotein cholesterol on coronary artery disease risk, the inverse-variance weighted method suggests a causal effect of both lipid fractions, whereas the weighted median and MR-Egger regression methods suggest a null effect of high-density lipoprotein cholesterol that corresponds with the experimental evidence. Both median-based and MR-Egger regression methods should be considered as sensitivity analyses for Mendelian randomization investigations with multiple genetic variants.

Estimating the causal influence of body mass index on risk of Parkinson disease: A Mendelian randomisation study.

BACKGROUND: Both positive and negative associations between higher body mass index (BMI) and Parkinson disease (PD) have been reported in observational studies, but it has been difficult to establish causality because of the possibility of residual confounding or reverse causation. To our knowledge, Mendelian randomisation (MR)-the use of genetic instrumental variables (IVs) to explore causal effects-has not previously been used to test the effect of BMI on PD. METHODS AND FINDINGS: Two-sample MR was undertaken using genome-wide association (GWA) study data. The associations between the genetic instruments and BMI were obtained from the GIANT consortium and consisted of the per-allele difference in mean BMI for 77 independent variants that reached genome-wide significance. The per-allele difference in log-odds of PD for each of these variants was estimated from a recent meta-analysis, which included 13,708 cases of PD and 95,282 controls. The inverse-variance weighted method was used to estimate a pooled odds ratio (OR) for the effect of a 5-kg/m2 higher BMI on PD. Evidence of directional pleiotropy averaged across all variants was sought using MR-Egger regression. Frailty simulations were used to assess whether causal associations were affected by mortality selection. A combined genetic IV expected to confer a lifetime exposure of 5-kg/m2 higher BMI was associated with a lower risk of PD (OR 0.82, 95% CI 0.69-0.98). MR-Egger regression gave similar results, suggesting that directional pleiotropy was unlikely to be biasing the result (intercept 0.002; p = 0.654). However, the apparent protective influence of higher BMI could be at least partially induced by survival bias in the PD GWA study, as demonstrated by frailty simulations. Other important limitations of this application of MR include the inability to analyse non-linear associations, to undertake subgroup analyses, and to gain mechanistic insights. CONCLUSIONS: In this large study using two-sample MR, we found that variants known to influence BMI had effects on PD in a manner consistent with higher BMI leading to lower risk of PD. The mechanism underlying this apparent protective effect warrants further study.

Investigating the association between genetically proxied circulating levels of immune checkpoint proteins and cancer survival: protocol for a Mendelian randomisation analysis.

INTRODUCTION: Compared with the traditional drug development pathway, investigating alternative uses for existing drugs (ie, drug repurposing) requires substantially less time, cost and resources. Immune checkpoint inhibitors are licensed for the treatment of certain breast, colorectal, head and neck, lung and melanoma cancers. These drugs target immune checkpoint proteins to reduce the suppression of T cell activation by cancer cells. As T cell suppression is a hallmark of cancer common across anatomical sites, we hypothesise that immune checkpoint inhibitors could be repurposed for the treatment of additional cancers beyond the ones already indicated. METHODS AND ANALYSIS: We will use two-sample Mendelian randomisation to investigate the effect of genetically proxied levels of protein targets of two immune checkpoint inhibitors-programmed cell death protein 1 and programmed death ligand 1-on survival of seven cancer types (breast, colorectal, head and neck, lung, melanoma, ovarian and prostate). Summary genetic association data will be obtained from prior genome-wide association studies of circulating protein levels and cancer survival in populations of European ancestry. Various sensitivity analyses will be performed to examine the robustness of findings to potential violations of Mendelian randomisation assumptions, collider bias and the impact of alternative genetic instrument construction strategies. The impact of treatment history and tumour stage on the findings will also be investigated using summary-level and individual-level genetic data where available. ETHICS AND DISSEMINATION: No separate ethics approval will be required for these analyses as we will be using data from previously published genome-wide association studies which individually gained ethical approval and participant consent. Results from analyses will be submitted as an open-access peer-reviewed publication and statistical code will be made freely available on the completion of the analysis.

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome.

Background: Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours. Methods: We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343). Results: Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity. Conclusions: This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas. Funding: Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17-022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

The association between genetically elevated polyunsaturated fatty acids and risk of cancer.

BACKGROUND: The causal relevance of polyunsaturated fatty acids (PUFAs) for risk of site-specific cancers remains uncertain. METHODS: Using a Mendelian randomization (MR) framework, we assessed the causal relevance of PUFAs for risk of cancer in European and East Asian ancestry individuals. We defined the primary exposure as PUFA desaturase activity, proxied by rs174546 at the FADS locus. Secondary exposures were defined as omega 3 and omega 6 PUFAs that could be proxied by genetic polymorphisms outside the FADS region. Our study used summary genetic data on 10 PUFAs and 67 cancers, corresponding to 562,871 cases and 1,619,465 controls, collected by the Fatty Acids in Cancer Mendelian Randomization Collaboration. We estimated odds ratios (ORs) for cancer per standard deviation increase in genetically proxied PUFA exposures. FINDINGS: Genetically elevated PUFA desaturase activity was associated (P < 0.0007) with higher risk (OR [95% confidence interval]) of colorectal cancer (1.09 [1.07-1.11]), esophageal squamous cell carcinoma (1.16 [1.06-1.26]), lung cancer (1.06 [1.03-1.08]) and basal cell carcinoma (1.05 [1.02-1.07]). There was little evidence for associations with reproductive cancers (OR = 1.00 [95% CI: 0.99-1.01]; Pheterogeneity = 0.25), urinary system cancers (1.03 [0.99-1.06], Pheterogeneity = 0.51), nervous system cancers (0.99 [0.95-1.03], Pheterogeneity = 0.92) or blood cancers (1.01 [0.98-1.04], Pheterogeneity = 0.09). Findings for colorectal cancer and esophageal squamous cell carcinoma remained compatible with causality in sensitivity analyses for violations of assumptions. Secondary MR analyses highlighted higher omega 6 PUFAs (arachidonic acid, gamma-linolenic acid and dihomo-gamma-linolenic acid) as potential mediators. PUFA biosynthesis is known to interact with aspirin, which increases risk of bleeding and inflammatory bowel disease. In a phenome-wide MR study of non-neoplastic diseases, we found that genetic lowering of PUFA desaturase activity, mimicking a hypothetical intervention to reduce cancer risk, was associated (P < 0.0006) with increased risk of inflammatory bowel disease but not bleeding. INTERPRETATION: The PUFA biosynthesis pathway may be an intervention target for prevention of colorectal cancer and esophageal squamous cell carcinoma but with potential for increased risk of inflammatory bowel disease. FUNDING: Cancer Resesrch UK (C52724/A20138, C18281/A19169). UK Medical Research Council (MR/P014054/1). National Institute for Health Research (NIHR202411). UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4). National Cancer Institute (R00 CA215360). National Institutes of Health (U01 CA164973, R01 CA60987, R01 CA72520, U01 CA74806, R01 CA55874, U01 CA164973 and U01 CA164973).

Design and quality control of large-scale two-sample Mendelian randomization studies.

Background: Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods: We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results: We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions: In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Trans-ethnic Mendelian-randomization study reveals causal relationships between cardiometabolic factors and chronic kidney disease.

BACKGROUND: This study was to systematically test whether previously reported risk factors for chronic kidney disease (CKD) are causally related to CKD in European and East Asian ancestries using Mendelian randomization. METHODS: A total of 45 risk factors with genetic data in European ancestry and 17 risk factors in East Asian participants were identified as exposures from PubMed. We defined the CKD by clinical diagnosis or by estimated glomerular filtration rate of <60 ml/min/1.73 m2. Ultimately, 51 672 CKD cases and 958 102 controls of European ancestry from CKDGen, UK Biobank and HUNT, and 13 093 CKD cases and 238 118 controls of East Asian ancestry from Biobank Japan, China Kadoorie Biobank and Japan-Kidney-Biobank/ToMMo were included. RESULTS: Eight risk factors showed reliable evidence of causal effects on CKD in Europeans, including genetically predicted body mass index (BMI), hypertension, systolic blood pressure, high-density lipoprotein cholesterol, apolipoprotein A-I, lipoprotein(a), type 2 diabetes (T2D) and nephrolithiasis. In East Asians, BMI, T2D and nephrolithiasis showed evidence of causality on CKD. In two independent replication analyses, we observed that increased hypertension risk showed reliable evidence of a causal effect on increasing CKD risk in Europeans but in contrast showed a null effect in East Asians. Although liability to T2D showed consistent effects on CKD, the effects of glycaemic phenotypes on CKD were weak. Non-linear Mendelian randomization indicated a threshold relationship between genetically predicted BMI and CKD, with increased risk at BMI of >25 kg/m2. CONCLUSIONS: Eight cardiometabolic risk factors showed causal effects on CKD in Europeans and three of them showed causality in East Asians, providing insights into the design of future interventions to reduce the burden of CKD.

Exploiting horizontal pleiotropy to search for causal pathways within a Mendelian randomization framework.

In Mendelian randomization (MR) analysis, variants that exert horizontal pleiotropy are typically treated as a nuisance. However, they could be valuable in identifying alternative pathways to the traits under investigation. Here, we develop MR-TRYX, a framework that exploits horizontal pleiotropy to discover putative risk factors for disease. We begin by detecting outliers in a single exposure-outcome MR analysis, hypothesising they are due to horizontal pleiotropy. We search across hundreds of complete GWAS summary datasets to systematically identify other (candidate) traits that associate with the outliers. We develop a multi-trait pleiotropy model of the heterogeneity in the exposure-outcome analysis due to pathways through candidate traits. Through detailed investigation of several causal relationships, many pleiotropic pathways are uncovered with already established causal effects, validating the approach, but also alternative putative causal pathways. Adjustment for pleiotropic pathways reduces the heterogeneity across the analyses.

Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases.

The human proteome is a major source of therapeutic targets. Recent genetic association analyses of the plasma proteome enable systematic evaluation of the causal consequences of variation in plasma protein levels. Here we estimated the effects of 1,002 proteins on 225 phenotypes using two-sample Mendelian randomization (MR) and colocalization. Of 413 associations supported by evidence from MR, 130 (31.5%) were not supported by results of colocalization analyses, suggesting that genetic confounding due to linkage disequilibrium is widespread in naïve phenome-wide association studies of proteins. Combining MR and colocalization evidence in cis-only analyses, we identified 111 putatively causal effects between 65 proteins and 52 disease-related phenotypes ( https://www.epigraphdb.org/pqtl/ ). Evaluation of data from historic drug development programs showed that target-indication pairs with MR and colocalization support were more likely to be approved, evidencing the value of this approach in identifying and prioritizing potential therapeutic targets.

Immune-mediated genetic pathways resulting in pulmonary function impairment increase lung cancer susceptibility.

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Fetal alcohol spectrum disorders: the epigenetic perspective.

Ethanol is a classic teratogen capable of inducing a wide range of developmental abnormalities. Studies in animal models suggest that differences in timing and dosage underlie this variability, with three particularly important developmental periods: preconception, preimplantation, and gastrulation. These periods of teratogenesis correlate with peak periods of epigenetic reprogramming which, together with the evidence that ethanol interferes with one-carbon metabolism, DNA methylation, histone modifications, and noncoding RNA, suggests an important role for epigenetic mechanisms in the etiology of fetal alcohol spectrum disorders (FASDs). In addition to a number of testable hypotheses, an epigenetic model suggests that the concept of a "fetal alcohol spectrum" should be expanded to include "preconceptional effects." This proposal has important public health implications, highlighting the urgency of research into the epigenetic basis of FASDs.

Exposure of mouse embryos to ethanol during preimplantation development: effect on DNA methylation in the h19 imprinting control region.

In the present study, it was hypothesized that disruption of imprinting control in the H19/Igf2 domain may be a mechanism of ethanol-induced growth retardation-a key clinical feature of the fetal alcohol spectrum disorders (FASD). To test this prediction, genomic bisulphite sequencing was carried out on 473 bp of the H19 imprinting control region in DNA obtained from midgestation F(1) hybrid mouse embryos (C57BL/6 x Mus musculus castaneus) exposed to ethanol during preimplantation development. Although ethanol-exposed placentae and embryos were severely growth retarded in comparison with saline-treated controls, DNA methylation at paternal and maternal alleles was unaffected in embryos. However, paternal alleles were significantly less methylated in ethanol-treated placentae in comparison with saline-treated controls. Partial correlations suggested that the relationship between ethanol and placental weight partly depended on DNA methylation at a CCCTC-binding factor site on the paternal allele in placentae, suggesting a novel mechanism of ethanol-induced growth retardation. In contrast, partial correlations suggested that embryo growth retardation was independent of placental growth retardation. Relaxation of allele-specific DNA methylation in control placentae in comparison with control embryos was also observed, consistent with a model of imprinting in which 1) regulation of allele-specific DNA methylation in the placenta depends on a stochastic interplay between silencer and enhancer chromatin assembly factors and 2) imprinting control mechanisms in the embryo are more robust to environmental perturbations.

Using the MR-Base platform to investigate risk factors and drug targets for thousands of phenotypes.

Mendelian randomization (MR) estimates the causal effect of exposures on outcomes by exploiting genetic variation to address confounding and reverse causation. This method has a broad range of applications, including investigating risk factors and appraising potential targets for intervention. MR-Base has become established as a freely accessible, online platform, which combines a database of complete genome-wide association study results with an interface for performing Mendelian randomization and sensitivity analyses. This allows the user to explore millions of potentially causal associations. MR-Base is available as a web application or as an R package. The technical aspects of the tool have previously been documented in the literature. The present article is complementary to this as it focuses on the applied aspects. Specifically, we describe how MR-Base can be used in several ways, including to perform novel causal analyses, replicate results and enable transparency, amongst others. We also present three use cases, which demonstrate important applications of Mendelian randomization and highlight the benefits of using MR-Base for these types of analyses.