Asthma and incident coronary heart disease: an observational and Mendelian randomisation study.

BACKGROUND: Observational studies suggest asthma is a risk factor for coronary heart disease (CHD) and sex modifies the risk, but they may suffer from methodological limitations. To overcome these, we applied a "triangulation approach", where different methodologies, with different potential biases, were leveraged to enhance confidence in findings. METHODS: First, we conducted an observational study using UK medical records to match asthma patients 1:1, by age, sex and general practitioner (GP) practice, to the general population. We measured the association between asthma and incident CHD (myocardial infarction: hospitalisation/death) by applying minimal sufficient adjustment: model 1, smoking, body mass index, oral corticosteroids, atopy and deprivation; model 2, additionally adjusting for healthcare behaviour (GP consultation frequency). Second, we conducted a Mendelian randomisation (MR) study using data from the UK Biobank, Trans-National Asthma Genetic Consortium (TAGC) and Coronary Artery Disease Genome-wide Replication and Meta-analysis consortium (CARDIoGRAM). Using 64 asthma single nucleotide polymorphisms, the effect of asthma on CHD was estimated with inverse variance-weighted meta-analysis and methods that adjust for pleiotropy. RESULTS: In our observational study (n=1 522 910), we found asthma was associated with 6% increased risk of CHD (model 1: HR 1.06, 95% CI 1.01-1.13); after accounting for healthcare behaviour, we found no association (model 2: HR 0.99, 95% CI 0.94-1.05). Asthma severity did not modify the association, but sex did (females: HR 1.11, 95% CI 1.01-1.21; males: HR 0.91, 95% CI 0.84-0.98). Our MR study (n=589 875) found no association between asthma and CHD (OR 1.01, 95% CI 0.98-1.04) and no modification by sex. CONCLUSIONS: Our findings suggest that asthma is not a risk factor for CHD. Previous studies may have suffered from detection bias or residual confounding.

The causal relationship between gastro-oesophageal reflux disease and idiopathic pulmonary fibrosis: a bidirectional two-sample Mendelian randomisation study.

BACKGROUND: Gastro-oesophageal reflux disease (GORD) is associated with idiopathic pulmonary fibrosis (IPF) in observational studies. It is not known if this association arises because GORD causes IPF or because IPF causes GORD, or because of confounding by factors, such as smoking, associated with both GORD and IPF. We used bidirectional Mendelian randomisation (MR), where genetic variants are used as instrumental variables to address issues of confounding and reverse causation, to examine how, if at all, GORD and IPF are causally related. METHODS: A bidirectional two-sample MR was performed to estimate the causal effect of GORD on IPF risk and of IPF on GORD risk, using genetic data from the largest GORD (78 707 cases and 288 734 controls) and IPF (4125 cases and 20 464 controls) genome-wide association meta-analyses currently available. RESULTS: GORD increased the risk of IPF, with an OR of 1.6 (95% CI 1.04-2.49; p=0.032). There was no evidence of a causal effect of IPF on the risk of GORD, with an OR of 0.999 (95% CI 0.997-1.000; p=0.245). CONCLUSIONS: We found that GORD increases the risk of IPF, but found no evidence that IPF increases the risk of GORD. GORD should be considered in future studies of IPF risk and interest in it as a potential therapeutic target should be renewed. The mechanisms underlying the effect of GORD on IPF should also be investigated.

A Mendelian randomization study investigating the causal role of inflammation on Parkinson's disease.

There is increasing evidence for inflammation as a determinant in the pathogenesis of Parkinson's disease, but its role in parkinsonian neurodegeneration remains elusive. It is not clear whether inflammatory cascades are causes or consequences of dopamine neuron death. In the present study, we aim to perform an in-depth statistical investigation of the causal relationship between inflammation and Parkinson's disease using a two-sample Mendelian randomization design. Genetic instruments were selected using summary-level data from the largest genome-wide association studies to date (sample size ranging from 13 955 to 204 402 individuals) conducted on a European population for the following inflammation biomarkers: C-reactive protein, interleukin-6, interleukin 1 receptor antagonist and tumour necrosis factor α. Genetic association data on Parkinson's disease (56 306 cases and 1 417 791 controls) and age at onset of Parkinson's disease (28 568 cases) were obtained from the International Parkinson's Disease Genomics Consortium. On primary analysis, causal associations were estimated on sets of strong (P-value < 5 × 10-8; F-statistic > 10) and independent (linkage disequilibrium r2 < 0.001) genetic instruments using the inverse-variance weighted method. In sensitivity analysis, we estimated causal effects using robust Mendelian randomization methods and after removing pleiotropic genetic variants. Reverse causation was also explored. We repeated the analysis on different data sources for inflammatory biomarkers to check the consistency of the findings. In all the three data sources selected for interleukin-6, we found statistical evidence for an earlier age at onset of Parkinson's disease associated with increased interleukin-6 concentration [years difference per 1 log-unit increase = -2.364, 95% confidence interval (CI) = -4.789-0.060; years difference per 1 log-unit increase = -2.011, 95% CI = -3.706 to -0.317; years difference per 1 log-unit increase = -1.569, 95% CI = -2.891 to -0.247]. We did not observe any statistical evidence for causal effects of C-reactive protein, interleukin 1 receptor antagonist and tumour necrosis factor α on both Parkinson's disease and its age at onset. Results after excluding possible pleiotropic genetic variants were consistent with findings from primary analyses. When investigating reverse causation, we did not find evidence for a causal effect of Parkinson's disease or age at onset on any biomarkers of inflammation. We found evidence for a causal association between the onset of Parkinson's disease and interleukin-6. The findings of this study suggest that the pro-inflammatory activity of the interleukin-6 cytokine could be a determinant of prodromal Parkinson's disease.

Age at menarche and adult body mass index: a Mendelian randomization study.

BACKGROUND: Pubertal timing has psychological and physical sequelae. While observational studies have demonstrated an association between age at menarche and adult body mass index (BMI), confounding makes it difficult to infer causality. METHODS: The Mendelian randomization (MR) technique is not limited by traditional confounding and was used to investigate the presence of a causal effect of age at menarche on adult BMI. MR uses genetic variants as instruments under the assumption that they act on BMI only through age at menarche (no pleiotropy). Using a two-sample MR approach, heterogeneity between the MR estimates from individual instruments was used as a proxy for pleiotropy, with sensitivity analyses performed if detected. Genetic instruments and estimates of their association with age at menarche were obtained from a genome-wide association meta-analysis on 182,416 women. The genetic effects on adult BMI were estimated using data on 80,465 women from the UK Biobank. The presence of a causal effect of age at menarche on adult BMI was further investigated using data on 70,692 women from the GIANT Consortium. RESULTS: There was evidence of pleiotropy among instruments. Using the UK Biobank data, after removing instruments associated with childhood BMI that were likely exerting pleiotropy, fixed-effect meta-analysis across instruments demonstrated that a 1 year increase in age at menarche reduces adult BMI by 0.38 kg/m2 (95% CI 0.25-0.51 kg/m2). However, evidence of pleiotropy remained. MR-Egger regression did not suggest directional bias, and similar estimates to the fixed-effect meta-analysis were obtained in sensitivity analyses when using a random-effect model, multivariable MR, MR-Egger regression, a weighted median estimator and a weighted mode-based estimator. The direction and significance of the causal effect were replicated using GIANT Consortium data. CONCLUSION: MR provides evidence to support the hypothesis that earlier age at menarche causes higher adult BMI. Complex hormonal and psychological factors may be responsible.

The Effect of Iron Status on Risk of Coronary Artery Disease: A Mendelian Randomization Study-Brief Report.

OBJECTIVE: Iron status is a modifiable trait that has been implicated in cardiovascular disease. This study uses the Mendelian randomization technique to investigate whether there is any causal effect of iron status on risk of coronary artery disease (CAD). APPROACH AND RESULTS: A 2-sample Mendelian randomization approach is used to estimate the effect of iron status on CAD risk. Three loci (rs1800562 and rs1799945 in the HFE gene and rs855791 in TMPRSS6) that are each associated with serum iron, transferrin saturation, ferritin, and transferrin in a pattern suggestive of an association with systemic iron status are used as instruments. SNP (single-nucleotide polymorphism)-iron status association estimates are based on a genome-wide association study meta-analysis of 48 972 individuals. SNP-CAD estimates are derived by combining the results of a genome-wide association study meta-analysis of 60 801 CAD cases and 123 504 controls with those of a meta-analysis of 63 746 CAD cases and 130 681 controls obtained from Metabochip and genome-wide association studies. Combined Mendelian randomization estimates are obtained for each marker by pooling results across the 3 instruments. We find evidence of a protective effect of higher iron status on CAD risk (iron odds ratio, 0.94 per SD unit increase; 95% confidence interval, 0.88-1.00; P=0.039; transferrin saturation odds ratio, 0.95 per SD unit increase; 95% confidence interval, 0.91-0.99; P=0.027; log-transformed ferritin odds ratio, 0.85 per SD unit increase; 95% confidence interval, 0.73-0.98; P=0.024; and transferrin odds ratio, 1.08 per SD unit increase; 95% confidence interval, 1.01-1.16; P=0.034). CONCLUSIONS: This Mendelian randomization study supports the hypothesis that higher iron status reduces CAD risk. These findings may highlight a therapeutic target.

Serum iron level and kidney function: a Mendelian randomization study.

Background: Iron depletion is a known consequence of chronic kidney disease (CKD), but there is contradicting epidemiological evidence on whether iron itself affects kidney function and whether its effect is protective or detrimental in the general population. While epidemiological studies tend to be affected by confounding and reverse causation, Mendelian randomization (MR) can provide unconfounded estimates of causal effects by using genes as instruments. Methods: We performed an MR study of the effect of serum iron levels on estimated glomerular filtration rate (eGFR), using genetic variants known to be associated with iron. MR estimates of the effect of iron on eGFR were derived based on the association of each variant with iron and eGFR from two large genome-wide meta-analyses on 48 978 and 74 354 individuals. We performed a similar MR analysis for ferritin, which measures iron stored in the body, using variants associated with ferritin. Results: A combined MR estimate across all variants showed a 1.3% increase in eGFR per standard deviation increase in iron (95% confidence interval 0.4­2.1%; P = 0.004). The results for ferritin were consistent with those for iron. Secondary MR analyses of the effects of iron and ferritin on CKD did not show significant associations but had very low statistical power. Conclusions: Our study suggests a protective effect of iron on kidney function in the general population. Further research is required to confirm this causal association, investigate it in study populations at higher risk of CKD and explore its underlying mechanism of action.

A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization.

Mendelian randomization (MR) uses genetic data to probe questions of causality in epidemiological research, by invoking the Instrumental Variable (IV) assumptions. In recent years, it has become commonplace to attempt MR analyses by synthesising summary data estimates of genetic association gleaned from large and independent study populations. This is referred to as two-sample summary data MR. Unfortunately, due to the sheer number of variants that can be easily included into summary data MR analyses, it is increasingly likely that some do not meet the IV assumptions due to pleiotropy. There is a pressing need to develop methods that can both detect and correct for pleiotropy, in order to preserve the validity of the MR approach in this context. In this paper, we aim to clarify how established methods of meta-regression and random effects modelling from mainstream meta-analysis are being adapted to perform this task. Specifically, we focus on two contrastin g approaches: the Inverse Variance Weighted (IVW) method which assumes in its simplest form that all genetic variants are valid IVs, and the method of MR-Egger regression that allows all variants to violate the IV assumptions, albeit in a specific way. We investigate the ability of two popular random effects models to provide robustness to pleiotropy under the IVW approach, and propose statistics to quantify the relative goodness-of-fit of the IVW approach over MR-Egger regression. © 2017 The Authors. Statistics in Medicine Published by JohnWiley & Sons Ltd.

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.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium.

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Bayesian analysis of censored response data in family-based genetic association studies.

Biomarkers are subject to censoring whenever some measurements are not quantifiable given a laboratory detection limit. Methods for handling censoring have received less attention in genetic epidemiology, and censored data are still often replaced with a fixed value. We compared different strategies for handling a left-censored continuous biomarker in a family-based study, where the biomarker is tested for association with a genetic variant, S, adjusting for a covariate, X. Allowing different correlations between X and S, we compared simple substitution of censored observations with the detection limit followed by a linear mixed effect model (LMM), Bayesian model with noninformative priors, Tobit model with robust standard errors, the multiple imputation (MI) with and without S in the imputation followed by a LMM. Our comparison was based on real and simulated data in which 20% and 40% censoring were artificially induced. The complete data were also analyzed with a LMM. In the MICROS study, the Bayesian model gave results closer to those obtained with the complete data. In the simulations, simple substitution was always the most biased method, the Tobit approach gave the least biased estimates at all censoring levels and correlation values, the Bayesian model and both MI approaches gave slightly biased estimates but smaller root mean square errors. On the basis of these results the Bayesian approach is highly recommended for candidate gene studies; however, the computationally simpler Tobit and the MI without S are both good options for genome-wide studies.

Discovery and fine mapping of serum protein loci through transethnic meta-analysis.

Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p < 5 × 10(-8)) for serum albumin (HPN-SCN1B, GCKR-FNDC4, SERPINF2-WDR81, TNFRSF11A-ZCCHC2, FRMD5-WDR76, and RPS11-FCGRT, in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein (TNFRS13B, 6q21.3, and ELL2, in up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

Coffee Intake, Plasma Caffeine Levels, and Kidney Function: Two-Sample Mendelian Randomization Among East Asian and European Ancestries.

Introduction: Previous Mendelian randomization (MR) studies for the coffee-kidney association have reported inconsistent relationships in European populations and never examined mediators of this association. We aimed to evaluate this causal relationship using two-sample MR among both East Asian and European ancestries and to explore underlying mechanisms using plasma caffeine levels. Methods: Among East Asians, the largest genome-wide association study (GWAS) results for coffee intake, plasma caffeine levels, and kidney outcomes were obtained from 152,634; 8940; and 47,070 Japanese adults. Among Europeans, summary statistics were acquired from European GWAS with 428,860; 7719; and 564,470 adults for each trait. We applied different MR methods (inverse-variance weighted [IVW] with random effects, weighted median, weighted mode, and MR-Egger). Results: After excluding possible pleiotropic variants, among East Asian ancestry, drinking an extra coffee intake per week showed a protective association on serum creatinine-based estimated glomerular filtration rate (eGFRcre) (ß = 0.077; 95% confidence interval [CI] = 0.003 to 0.150). Analysis in European ancestry also showed a causal relationship between drinking an extra coffee intake per day and eGFRcre (ß = 0.052; 95% CI = 0.027 to 0.078). These results were consistent across different MR methods accounting for invalid instruments. Higher plasma caffeine levels were associated with lower eGFRcre among both East Asian (ß = -0.071; 95% CI = -0.137 to -0.006) and European ancestries (ß = -0.048; 95% CI = -0.057 to -0.040). Conclusions: Our cross-ancestry MR study found beneficial effects of coffee intake on eGFRcre. However, given the possible adverse effects of plasma caffeine levels on eGFRcre, interpretation of the results should be carefully considered and further investigations on noncaffeine and biological pathways are needed.

Genetic determinants of complement activation in the general population.

Complement is a fundamental innate immune response component. Its alterations are associated with severe systemic diseases. To illuminate the complement's genetic underpinnings, we conduct genome-wide association studies of the functional activity of the classical (CP), lectin (LP), and alternative (AP) complement pathways in the Cooperative Health Research in South Tyrol study (n = 4,990). We identify seven loci, encompassing 13 independent, pathway-specific variants located in or near complement genes (CFHR4, C7, C2, MBL2) and non-complement genes (PDE3A, TNXB, ABO), explaining up to 74% of complement pathways' genetic heritability and implicating long-range haplotypes associated with LP at MBL2. Two-sample Mendelian randomization analyses, supported by transcriptome- and proteome-wide colocalization, confirm known causal pathways, establish within-complement feedback loops, and implicate causality of ABO on LP and of CFHR2 and C7 on AP. LP causally influences collectin-11 and KAAG1 levels and the risk of mouth ulcers. These results build a comprehensive resource to investigate the role of complement in human health.

Serum iron levels and the risk of Parkinson disease: a Mendelian randomization study.

BACKGROUND: Although levels of iron are known to be increased in the brains of patients with Parkinson disease (PD), epidemiological evidence on a possible effect of iron blood levels on PD risk is inconclusive, with effects reported in opposite directions. Epidemiological studies suffer from problems of confounding and reverse causation, and mendelian randomization (MR) represents an alternative approach to provide unconfounded estimates of the effects of biomarkers on disease. We performed a MR study where genes known to modify iron levels were used as instruments to estimate the effect of iron on PD risk, based on estimates of the genetic effects on both iron and PD obtained from the largest sample meta-analyzed to date. METHODS AND FINDINGS: We used as instrumental variables three genetic variants influencing iron levels, HFE rs1800562, HFE rs1799945, and TMPRSS6 rs855791. Estimates of their effect on serum iron were based on a recent genome-wide meta-analysis of 21,567 individuals, while estimates of their effect on PD risk were obtained through meta-analysis of genome-wide and candidate gene studies with 20,809 PD cases and 88,892 controls. Separate MR estimates of the effect of iron on PD were obtained for each variant and pooled by meta-analysis. We investigated heterogeneity across the three estimates as an indication of possible pleiotropy and found no evidence of it. The combined MR estimate showed a statistically significant protective effect of iron, with a relative risk reduction for PD of 3% (95% CI 1%-6%; p = 0.001) per 10 µg/dl increase in serum iron. CONCLUSIONS: Our study suggests that increased iron levels are causally associated with a decreased risk of developing PD. Further studies are needed to understand the pathophysiological mechanism of action of serum iron on PD risk before recommendations can be made.

Genome-wide association analysis and fine mapping of NT-proBNP level provide novel insight into the role of the MTHFR-CLCN6-NPPA-NPPB gene cluster.

High blood concentration of the N-terminal cleavage product of the B-type natriuretic peptide (NT-proBNP) is strongly associated with cardiac dysfunction and is increasingly used for heart failure diagnosis. To identify genetic variants associated with NT-proBNP level, we performed a genome-wide association analysis in 1325 individuals from South Tyrol, Italy, and followed up the most significant results in 1746 individuals from two German population-based studies. A genome-wide significant signal in the MTHFR-CLCN6-NPPA-NPPB gene cluster was replicated, after correction for multiple testing (replication one-sided P-value = 8.4 × 10(-10)). A conditional regression analysis of 128 single-nucleotide polymorphisms in the region of interest identified novel variants in the CLCN6 gene as independently associated with NT-proBNP. In this locus, four haplotypes were associated with increased NT-proBNP levels (haplotype-specific combined P-values from 8.3 × 10(-03) to 9.3 × 10(-11)). The observed increase in the NT-proBNP level was proportional to the number of haplotype copies present (i.e. dosage effect), with an increase associated with two copies that varied between 20 and 100 pg/ml across populations. The identification of novel variants in the MTHFR-CLCN6-NPPA-NPPB cluster provides new insights into the biological mechanisms of cardiac dysfunction.

Molecular phenotypes of mitochondrial dysfunction in clinically non-manifesting heterozygous PRKN variant carriers.

Homozygous or compound heterozygous (biallelic) variants in PRKN are causal for PD with highly penetrant symptom expression, while the much more common heterozygous variants may predispose to PD with highly reduced penetrance, through altered mitochondrial function. In the presence of pathogenic heterozygous variants, it is therefore important to test for mitochondrial alteration in cells derived from variant carriers to establish potential presymptomatic molecular markers. We generated lymphoblasts (LCLs) and human induced pluripotent stem cell (hiPSC)-derived neurons from non-manifesting heterozygous PRKN variant carriers and tested them for mitochondrial functionality. In LCLs, we detected hyperactive mitochondrial respiration, and, although milder compared to a biallelic PRKN-PD patient, hiPSC-derived neurons of non-manifesting heterozygous variant carriers also displayed several phenotypes of altered mitochondrial function. Overall, we identified molecular phenotypes that might be used to monitor heterozygous PRKN variant carriers during the prodromal phase. Such markers might also be useful to identify individuals at greater risk of eventual disease development and for testing potential mitochondrial function-based neuroprotective therapies before neurodegeneration advances.

Systematic review of Mendelian randomization studies on Parkinson's disease.

Background: Parkinson's disease (PD) is known to be associated with non-genetic factors. To infer causality, Mendelian randomization (MR) studies are increasingly used. Here, genetic variants are used as instrumental variables for the risk factor but have no direct effect on PD themselves. Methods: We performed a systematic literature review on MR studies for PD. Studies were identified searching the PubMed database. Upon data extraction, we evaluated the methodological quality and summarized the evidence. Results: Twelve articles were included. Most studies showed "good" methodological quality, but most did not report proper power estimations. Twelve analyses yielded nominally significant effects. Conclusions: Our systematic review shows that most MR studies were well performed and allow to identify causal exposures, which may inform further studies on the prevention and early intervention of PD.

Whole Exome Sequencing Enhanced Imputation Identifies 85 Metabolite Associations in the Alpine CHRIS Cohort.

Metabolites are intermediates or end products of biochemical processes involved in both health and disease. Here, we take advantage of the well-characterized Cooperative Health Research in South Tyrol (CHRIS) study to perform an exome-wide association study (ExWAS) on absolute concentrations of 175 metabolites in 3294 individuals. To increase power, we imputed the identified variants into an additional 2211 genotyped individuals of CHRIS. In the resulting dataset of 5505 individuals, we identified 85 single-variant genetic associations, of which 39 have not been reported previously. Fifteen associations emerged at ten variants with >5-fold enrichment in CHRIS compared to non-Finnish Europeans reported in the gnomAD database. For example, the CHRIS-enriched ETFDH stop gain variant p.Trp286Ter (rs1235904433-hexanoylcarnitine) and the MCCC2 stop lost variant p.Ter564GlnextTer3 (rs751970792-carnitine) have been found in patients with glutaric acidemia type II and 3-methylcrotonylglycinuria, respectively, but the loci have not been associated with the respective metabolites in a genome-wide association study (GWAS) previously. We further identified three gene-trait associations, where multiple rare variants contribute to the signal. These results not only provide further evidence for previously described associations, but also describe novel genes and mechanisms for diseases and disease-related traits.

A multi-omics study of circulating phospholipid markers of blood pressure.

High-throughput techniques allow us to measure a wide-range of phospholipids which can provide insight into the mechanisms of hypertension. We aimed to conduct an in-depth multi-omics study of various phospholipids with systolic blood pressure (SBP) and diastolic blood pressure (DBP). The associations of blood pressure and 151 plasma phospholipids measured by electrospray ionization tandem mass spectrometry were performed by linear regression in five European cohorts (n = 2786 in discovery and n = 1185 in replication). We further explored the blood pressure-related phospholipids in Erasmus Rucphen Family (ERF) study by associating them with multiple cardiometabolic traits (linear regression) and predicting incident hypertension (Cox regression). Mendelian Randomization (MR) and phenome-wide association study (Phewas) were also explored to further investigate these association results. We identified six phosphatidylethanolamines (PE 38:3, PE 38:4, PE 38:6, PE 40:4, PE 40:5 and PE 40:6) and two phosphatidylcholines (PC 32:1 and PC 40:5) which together predicted incident hypertension with an area under the ROC curve (AUC) of 0.61. The identified eight phospholipids are strongly associated with triglycerides, obesity related traits (e.g. waist, waist-hip ratio, total fat percentage, body mass index, lipid-lowering medication, and leptin), diabetes related traits (e.g. glucose, insulin resistance and insulin) and prevalent type 2 diabetes. The genetic determinants of these phospholipids also associated with many lipoproteins, heart rate, pulse rate and blood cell counts. No significant association was identified by bi-directional MR approach. We identified eight blood pressure-related circulating phospholipids that have a predictive value for incident hypertension. Our cross-omics analyses show that phospholipid metabolites in the circulation may yield insight into blood pressure regulation and raise a number of testable hypothesis for future research.

The use of two-sample methods for Mendelian randomization analyses on single large datasets.

BACKGROUND: With genome-wide association data for many exposures and outcomes now available from large biobanks, one-sample Mendelian randomization (MR) is increasingly used to investigate causal relationships. Many robust MR methods are available to address pleiotropy, but these assume independence between the gene-exposure and gene-outcome association estimates. Unlike in two-sample MR, in one-sample MR the two estimates are obtained from the same individuals, and the assumption of independence does not hold in the presence of confounding. METHODS: With simulations mimicking a typical study in UK Biobank, we assessed the performance, in terms of bias and precision of the MR estimate, of the fixed-effect and (multiplicative) random-effects meta-analysis method, weighted median estimator, weighted mode estimator and MR-Egger regression, used in both one-sample and two-sample data. We considered scenarios differing by the: presence/absence of a true causal effect; amount of confounding; and presence and type of pleiotropy (none, balanced or directional). RESULTS: Even in the presence of substantial correlation due to confounding, all two-sample methods used in one-sample MR performed similarly to when used in two-sample MR, except for MR-Egger which resulted in bias reflecting direction and magnitude of the confounding. Such bias was much reduced in the presence of very high variability in instrument strength across variants (IGX2 of 97%). CONCLUSIONS: Two-sample MR methods can be safely used for one-sample MR performed within large biobanks, expect for MR-Egger. MR-Egger is not recommended for one-sample MR unless the correlation between the gene-exposure and gene-outcome estimates due to confounding can be kept low, or the variability in instrument strength is very high.