A Bayesian approach to Mendelian randomization using summary statistics in the univariable and multivariable settings with correlated pleiotropy.

Mendelian randomization uses genetic variants as instrumental variables to make causal inferences on the effect of an exposure on an outcome. Due to the recent abundance of high-powered genome-wide association studies, many putative causal exposures of interest have large numbers of independent genetic variants with which they associate, each representing a potential instrument for use in a Mendelian randomization analysis. Such polygenic analyses increase the power of the study design to detect causal effects; however, they also increase the potential for bias due to instrument invalidity. Recent attention has been given to dealing with bias caused by correlated pleiotropy, which results from violation of the "instrument strength independent of direct effect" assumption. Although methods have been proposed that can account for this bias, a number of restrictive conditions remain in many commonly used techniques. In this paper, we propose a Bayesian framework for Mendelian randomization that provides valid causal inference under very general settings. We propose the methods MR-Horse and MVMR-Horse, which can be performed without access to individual-level data, using only summary statistics of the type commonly published by genome-wide association studies, and can account for both correlated and uncorrelated pleiotropy. In simulation studies, we show that the approach retains type I error rates below nominal levels even in high-pleiotropy scenarios. We demonstrate the proposed approaches in applied examples in both univariable and multivariable settings, some with very weak instruments.

Multi-response Mendelian randomization: Identification of shared and distinct exposures for multimorbidity and multiple related disease outcomes.

The existing framework of Mendelian randomization (MR) infers the causal effect of one or multiple exposures on one single outcome. It is not designed to jointly model multiple outcomes, as would be necessary to detect causes of more than one outcome and would be relevant to model multimorbidity or other related disease outcomes. Here, we introduce multi-response Mendelian randomization (MR2), an MR method specifically designed for multiple outcomes to identify exposures that cause more than one outcome or, conversely, exposures that exert their effect on distinct responses. MR2 uses a sparse Bayesian Gaussian copula regression framework to detect causal effects while estimating the residual correlation between summary-level outcomes, i.e., the correlation that cannot be explained by the exposures, and vice versa. We show both theoretically and in a comprehensive simulation study how unmeasured shared pleiotropy induces residual correlation between outcomes irrespective of sample overlap. We also reveal how non-genetic factors that affect more than one outcome contribute to their correlation. We demonstrate that by accounting for residual correlation, MR2 has higher power to detect shared exposures causing more than one outcome. It also provides more accurate causal effect estimates than existing methods that ignore the dependence between related responses. Finally, we illustrate how MR2 detects shared and distinct causal exposures for five cardiovascular diseases in two applications considering cardiometabolic and lipidomic exposures and uncovers residual correlation between summary-level outcomes reflecting known relationships between cardiovascular diseases.

Using genetic association data to guide drug discovery and development: Review of methods and applications.

Evidence on the validity of drug targets from randomized trials is reliable but typically expensive and slow to obtain. In contrast, evidence from conventional observational epidemiological studies is less reliable because of the potential for bias from confounding and reverse causation. Mendelian randomization is a quasi-experimental approach analogous to a randomized trial that exploits naturally occurring randomization in the transmission of genetic variants. In Mendelian randomization, genetic variants that can be regarded as proxies for an intervention on the proposed drug target are leveraged as instrumental variables to investigate potential effects on biomarkers and disease outcomes in large-scale observational datasets. This approach can be implemented rapidly for a range of drug targets to provide evidence on their effects and thus inform on their priority for further investigation. In this review, we present statistical methods and their applications to showcase the diverse opportunities for applying Mendelian randomization in guiding clinical development efforts, thus enabling interventions to target the right mechanism in the right population group at the right time. These methods can inform investigators on the mechanisms underlying drug effects, their related biomarkers, implications for the timing of interventions, and the population subgroups that stand to gain the most benefit. Most methods can be implemented with publicly available data on summarized genetic associations with traits and diseases, meaning that the only major limitations to their usage are the availability of appropriately powered studies for the exposure and outcome and the existence of a suitable genetic proxy for the proposed intervention.

Relaxing parametric assumptions for non-linear Mendelian randomization using a doubly-ranked stratification method.

Non-linear Mendelian randomization is an extension to standard Mendelian randomization to explore the shape of the causal relationship between an exposure and outcome using an instrumental variable. A stratification approach to non-linear Mendelian randomization divides the population into strata and calculates separate instrumental variable estimates in each stratum. However, the standard implementation of stratification, referred to as the residual method, relies on strong parametric assumptions of linearity and homogeneity between the instrument and the exposure to form the strata. If these stratification assumptions are violated, the instrumental variable assumptions may be violated in the strata even if they are satisfied in the population, resulting in misleading estimates. We propose a new stratification method, referred to as the doubly-ranked method, that does not require strict parametric assumptions to create strata with different average levels of the exposure such that the instrumental variable assumptions are satisfied within the strata. Our simulation study indicates that the doubly-ranked method can obtain unbiased stratum-specific estimates and appropriate coverage rates even when the effect of the instrument on the exposure is non-linear or heterogeneous. Moreover, it can also provide unbiased estimates when the exposure is coarsened (that is, rounded, binned into categories, or truncated), a scenario that is common in applied practice and leads to substantial bias in the residual method. We applied the proposed doubly-ranked method to investigate the effect of alcohol intake on systolic blood pressure, and found evidence of a positive effect of alcohol intake, particularly at higher levels of alcohol consumption.

Uncovering associations between pre-existing conditions and COVID-19 Severity: A polygenic risk score approach across three large biobanks.

OBJECTIVE: To overcome the limitations associated with the collection and curation of COVID-19 outcome data in biobanks, this study proposes the use of polygenic risk scores (PRS) as reliable proxies of COVID-19 severity across three large biobanks: the Michigan Genomics Initiative (MGI), UK Biobank (UKB), and NIH All of Us. The goal is to identify associations between pre-existing conditions and COVID-19 severity. METHODS: Drawing on a sample of more than 500,000 individuals from the three biobanks, we conducted a phenome-wide association study (PheWAS) to identify associations between a PRS for COVID-19 severity, derived from a genome-wide association study on COVID-19 hospitalization, and clinical pre-existing, pre-pandemic phenotypes. We performed cohort-specific PRS PheWAS and a subsequent fixed-effects meta-analysis. RESULTS: The current study uncovered 23 pre-existing conditions significantly associated with the COVID-19 severity PRS in cohort-specific analyses, of which 21 were observed in the UKB cohort and two in the MGI cohort. The meta-analysis yielded 27 significant phenotypes predominantly related to obesity, metabolic disorders, and cardiovascular conditions. After adjusting for body mass index, several clinical phenotypes, such as hypercholesterolemia and gastrointestinal disorders, remained associated with an increased risk of hospitalization following COVID-19 infection. CONCLUSION: By employing PRS as a proxy for COVID-19 severity, we corroborated known risk factors and identified novel associations between pre-existing clinical phenotypes and COVID-19 severity. Our study highlights the potential value of using PRS when actual outcome data may be limited or inadequate for robust analyses.

Robust use of phenotypic heterogeneity at drug target genes for mechanistic insights: Application of cis-multivariable Mendelian randomization to GLP1R gene region.

Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.

Combining evidence from Mendelian randomization and colocalization: Review and comparison of approaches.

Mendelian randomization and colocalization are two statistical approaches that can be applied to summarized data from genome-wide association studies (GWASs) to understand relationships between traits and diseases. However, despite similarities in scope, they are different in their objectives, implementation, and interpretation, in part because they were developed to serve different scientific communities. Mendelian randomization assesses whether genetic predictors of an exposure are associated with the outcome and interprets an association as evidence that the exposure has a causal effect on the outcome, whereas colocalization assesses whether two traits are affected by the same or distinct causal variants. When considering genetic variants in a single genetic region, both approaches can be performed. While a positive colocalization finding typically implies a non-zero Mendelian randomization estimate, the reverse is not generally true: there are several scenarios which would lead to a non-zero Mendelian randomization estimate but lack evidence for colocalization. These include the existence of distinct but correlated causal variants for the exposure and outcome, which would violate the Mendelian randomization assumptions, and a lack of strong associations with the outcome. As colocalization was developed in the GWAS tradition, typically evidence for colocalization is concluded only when there is strong evidence for associations with both traits. In contrast, a non-zero estimate from Mendelian randomization can be obtained despite only nominally significant genetic associations with the outcome at the locus. In this review, we discuss how the two approaches can provide complementary information on potential therapeutic targets.

Noise-augmented directional clustering of genetic association data identifies distinct mechanisms underlying obesity.

Clustering genetic variants based on their associations with different traits can provide insight into their underlying biological mechanisms. Existing clustering approaches typically group variants based on the similarity of their association estimates for various traits. We present a new procedure for clustering variants based on their proportional associations with different traits, which is more reflective of the underlying mechanisms to which they relate. The method is based on a mixture model approach for directional clustering and includes a noise cluster that provides robustness to outliers. The procedure performs well across a range of simulation scenarios. In an applied setting, clustering genetic variants associated with body mass index generates groups reflective of distinct biological pathways. Mendelian randomization analyses support that the clusters vary in their effect on coronary heart disease, including one cluster that represents elevated body mass index with a favourable metabolic profile and reduced coronary heart disease risk. Analysis of the biological pathways underlying this cluster identifies inflammation as potentially explaining differences in the effects of increased body mass index on coronary heart disease.

Birth weight influences cardiac structure, function, and disease risk: evidence of a causal association.

BACKGROUND AND AIMS: Low birth weight is a common pregnancy complication, which has been associated with higher risk of cardiometabolic disease in later life. Prior Mendelian randomization (MR) studies exploring this question do not distinguish the mechanistic contributions of variants that directly influence birth weight through the foetal genome (direct foetal effects), vs. variants influencing birth weight indirectly by causing an adverse intrauterine environment (indirect maternal effects). In this study, MR was used to assess whether birth weight, independent of intrauterine influences, is associated with cardiovascular disease risk and measures of adverse cardiac structure and function. METHODS: Uncorrelated (r2 < .001), genome-wide significant (P < 5 × 10-8) single nucleotide polymorphisms were extracted from genome-wide association studies summary statistics for birth weight overall, and after isolating direct foetal effects only. Inverse-variance weighted MR was utilized for analyses on outcomes of atrial fibrillation, coronary artery disease, heart failure, ischaemic stroke, and 16 measures of cardiac structure and function. Multiple comparisons were accounted for by Benjamini-Hochberg correction. RESULTS: Lower genetically-predicted birth weight, isolating direct foetal effects only, was associated with an increased risk of coronary artery disease (odds ratio 1.21, 95% confidence interval 1.06-1.37; P = .031), smaller chamber volumes, and lower stroke volume, but higher contractility. CONCLUSIONS: The results of this study support a causal role of low birth weight in cardiovascular disease, even after accounting for the influence of the intrauterine environment. This suggests that individuals with a low birth weight may benefit from early targeted cardiovascular disease prevention strategies, independent of whether this was linked to an adverse intrauterine environment during gestation.

Statistical methods for cis-Mendelian randomization with two-sample summary-level data.

Mendelian randomization (MR) is the use of genetic variants to assess the existence of a causal relationship between a risk factor and an outcome of interest. Here, we focus on two-sample summary-data MR analyses with many correlated variants from a single gene region, particularly on cis-MR studies which use protein expression as a risk factor. Such studies must rely on a small, curated set of variants from the studied region; using all variants in the region requires inverting an ill-conditioned genetic correlation matrix and results in numerically unstable causal effect estimates. We review methods for variable selection and estimation in cis-MR with summary-level data, ranging from stepwise pruning and conditional analysis to principal components analysis, factor analysis, and Bayesian variable selection. In a simulation study, we show that the various methods have comparable performance in analyses with large sample sizes and strong genetic instruments. However, when weak instrument bias is suspected, factor analysis and Bayesian variable selection produce more reliable inferences than simple pruning approaches, which are often used in practice. We conclude by examining two case studies, assessing the effects of low-density lipoprotein-cholesterol and serum testosterone on coronary heart disease risk using variants in the HMGCR and SHBG gene regions, respectively.

Risk Factors for Intracerebral Hemorrhage: Genome-Wide Association Study and Mendelian Randomization Analyses.

BACKGROUND: The genetic and nongenetic causes of intracerebral hemorrhage (ICH) remain obscure. The present study aimed to uncover the genetic and modifiable risk factors for ICH. METHODS: We meta-analyzed genome-wide association study data from 3 European biobanks, involving 7605 ICH cases and 711 818 noncases, to identify the genomic loci linked to ICH. To uncover the potential causal associations of cardiometabolic and lifestyle factors with ICH, we performed Mendelian randomization analyses using genetic instruments identified in previous genome-wide association studies of the exposures and ICH data from the present genome-wide association study meta-analysis. We performed multivariable Mendelian randomization analyses to examine the independent associations of the identified risk factors with ICH and evaluate potential mediating pathways. RESULTS: We identified 1 ICH risk locus, located at the APOE genomic region. The lead variant in this locus was rs429358 (chr19:45411941), which was associated with an odds ratio of ICH of 1.17 (95% CI, 1.11-1.20; P=6.01×10-11) per C allele. Genetically predicted higher levels of body mass index, visceral adiposity, diastolic blood pressure, systolic blood pressure, and lifetime smoking index, as well as genetic liability to type 2 diabetes, were associated with higher odds of ICH after multiple testing corrections. Additionally, a genetic increase in waist-to-hip ratio and liability to smoking initiation were consistently associated with ICH, albeit at the nominal significance level (P<0.05). Multivariable Mendelian randomization analysis showed that the association between body mass index and ICH was attenuated on adjustment for type 2 diabetes and further that type 2 diabetes may be a mediator of the body mass index-ICH relationship. CONCLUSIONS: Our findings indicate that the APOE locus contributes to ICH genetic susceptibility in European populations. Excess adiposity, elevated blood pressure, type 2 diabetes, and smoking were identified as the chief modifiable cardiometabolic and lifestyle factors for ICH.

LDL-c Lowering, Ischemic Stroke, and Small Vessel Disease Brain Imaging Biomarkers: A Mendelian Randomization Study.

BACKGROUND: The effects of lipid-lowering drug targets on different ischemic stroke subtypes are not fully understood. We aimed to explore the mechanisms by which lipid-lowering drug targets differentially affect the risk of ischemic stroke subtypes and their underlying pathophysiology. METHODS: Using a 2-sample Mendelian randomization approach, we assessed the effects of genetically proxied low-density lipoprotein cholesterol (LDL-c) and 3 clinically approved LDL-lowering drugs (HMGCR [3-hydroxy-3-methylglutaryl-CoA reductase], PCSK9 [proprotein convertase subtilisin/kexin type 9], and NPC1L1 [Niemann-Pick C1-Like 1]) on stroke subtypes and brain imaging biomarkers associated with small vessel stroke (SVS), including white matter hyperintensity volume and perivascular spaces. RESULTS: In genome-wide Mendelian randomization analyses, lower genetically predicted LDL-c was significantly associated with a reduced risk of any stroke, ischemic stroke, and large artery stroke, supporting previous findings. Significant associations between genetically predicted LDL-c and cardioembolic stroke, SVS, and biomarkers, perivascular space and white matter hyperintensity volume, were not identified in this study. In drug-target Mendelian randomization analysis, genetically proxied reduced LDL-c through NPC1L1 inhibition was associated with lower odds of perivascular space (odds ratio per 1-mg/dL decrease, 0.79 [95% CI, 0.67-0.93]) and with lower odds of SVS (odds ratio, 0.29 [95% CI, 0.10-0.85]). CONCLUSIONS: This study provides supporting evidence of a potentially protective effect of LDL-c lowering through NPC1L1 inhibition on perivascular space and SVS risk, highlighting novel therapeutic targets for SVS.

Comparison of caffeine consumption behavior with plasma caffeine levels as exposure measures in drug-target Mendelian randomization.

Mendelian randomization is an epidemiological technique that can explore the potential effect of perturbing a pharmacological target. Plasma caffeine levels can be used as a biomarker to measure the pharmacological effects of caffeine. Alternatively, this can be assessed using a behavioral proxy, such as average number of caffeinated drinks consumed per day. Either variable can be used as the exposure in a Mendelian randomization investigation, and to select which genetic variants to use as instrumental variables. Another possibility is to choose variants in gene regions with known biological relevance to caffeine level regulation. These choices affect the causal question that is being addressed by the analysis, and the validity of the analysis assumptions. Further, even when using the same genetic variants, the sign of Mendelian randomization estimates (positive or negative) can change depending on the choice of exposure. Some genetic variants that decrease caffeine metabolism associate with higher levels of plasma caffeine, but lower levels of caffeine consumption, as individuals with these variants require less caffeine consumption for the same physiological effect. We explore Mendelian randomization estimates for the effect of caffeine on body mass index, and discuss implications for variant and exposure choice in drug target Mendelian randomization investigations.

Genetic investigation into the broad health implications of caffeine: evidence from phenome-wide, proteome-wide and metabolome-wide Mendelian randomization.

BACKGROUND: Caffeine is one of the most utilized drugs in the world, yet its clinical effects are not fully understood. Circulating caffeine levels are influenced by the interplay between consumption behaviour and metabolism. This study aimed to investigate the effects of circulating caffeine levels by considering genetically predicted variation in caffeine metabolism. METHODS: Leveraging genetic variants related to caffeine metabolism that affect its circulating levels, we investigated the clinical effects of plasma caffeine in a phenome-wide association study (PheWAS). We validated novel findings using a two-sample Mendelian randomization framework and explored the potential mechanisms underlying these effects in proteome-wide and metabolome-wide Mendelian randomization. RESULTS: Higher levels of genetically predicted circulating caffeine among caffeine consumers were associated with a lower risk of obesity (odds ratio (OR) per standard deviation increase in caffeine = 0.97, 95% confidence interval (CI) CI: 0.95-0.98, p = 2.47 × 10-4), osteoarthrosis (OR = 0.97, 95% CI: 0.96-0.98, P=1.10 × 10-8) and osteoarthritis (OR: 0.97, 95% CI: 0.96 to 0.98, P = 1.09 × 10-6). Approximately one third of the protective effect of plasma caffeine on osteoarthritis risk was estimated to be mediated through lower bodyweight. Proteomic and metabolomic perturbations indicated lower chronic inflammation, improved lipid profiles, and altered protein and glycogen metabolism as potential biological mechanisms underlying these effects. CONCLUSIONS: We report novel evidence suggesting that long-term increases in circulating caffeine may reduce bodyweight and the risk of osteoarthrosis and osteoarthritis. We confirm prior genetic evidence of a protective effect of plasma caffeine on risk of overweight and obesity. Further clinical study is warranted to understand the translational relevance of these findings before clinical practice or lifestyle interventions related to caffeine consumption are introduced.

Modifiable risk factors and inflammation-related proteins in polymyalgia rheumatica: genome-wide meta-analysis and Mendelian randomisation.

OBJECTIVE: Polymyalgia rheumatica (PMR) is an age-related inflammatory disease of unknown cause. We aimed to identify potentially modifiable risk factors and therapeutic targets for preventing or treating PMR. METHODS: We meta-analysed genetic association data from 8,156 cases of PMR (defined using diagnostic codes and self-report) and 416,495 controls of European ancestry from the UK Biobank and FinnGen. We then performed Mendelian randomization analyses to estimate the association between eight modifiable risk factors (using data from up to 1.2 million individuals) and 65 inflammation-related circulating proteins (up to 55,792 individuals), using the inverse variance weighted and pleiotropy robust methods. RESULTS: We identified three novel genome-wide significant loci in the IL1R1, NEK6 and CCDC88B genes and confirmation of previously described associations with HLA-DRB1 and ANKRD55. Genetically predicted smoking intensity (OR 1.32; 95%CI 1.08-1.60; p = 0.006) and visceral adiposity (OR 1.22; 95%CI 1.10-1.37; p = 3.10x10-4) were associated with PMR susceptibility. Multiple circulating proteins related to IL-1 family signaling were associated with PMR. IL-1 receptor-like 2, also known as IL-36 receptor (OR 1.25; p = 1.89x10-32), serum amyloid A2 (OR 1.06, 9.91x10-10) and CXCL6 (OR 1.09, p = 4.85x10-7) retained significance after correction for multiple testing. CONCLUSION: Reducing smoking and visceral adiposity at a population level might reduce incidence of PMR. We identified proteins that may play causal roles in PMR, potentially suggesting new therapeutic opportunities. Further research is needed before these findings are applied to clinical practice.

Genetically predicted plasma cortisol and common chronic diseases: A Mendelian randomization study.

OBJECTIVE: Cushing's syndrome is characterized by hypercortisolaemia and is frequently accompanied by comorbidities such as type 2 diabetes, hypertension, osteoporosis, depression and schizophrenia. It is unclear whether moderate but lifelong hypercortisolaemia is causally associated with these diseases in the general population. We aimed to address this research gap using a Mendelian randomization approach. METHODS: We used three cortisol-associated genetic variants in the SERPINA6/SERPINA1 region as genetic instruments in a two-sample, inverse-variance-weighted Mendelian randomization analysis. We obtained summary-level statistics for cortisol and disease outcomes from publicly available genetic consortia, and meta-analysed them as appropriate. We conducted a multivariable Mendelian randomization analysis to assess potential mediating effects. RESULTS: A 1 standard deviation higher genetically predicted plasma cortisol was associated with greater odds of hypertension (odds ratio: 1.12; 95% confidence interval [CI]: 1.05-1.18) as well as higher systolic blood pressure (mean difference [MD]: 0.03 SD change; 95% CI: 0.01-0.05) and diastolic blood pressure (MD: 0.03 SD change; 95% CI: 0.01-0.04). There was no evidence of association with type 2 diabetes, osteoporosis, depression and schizophrenia. The association with hypertension was attenuated upon adjustment for waist circumference, suggesting potential mediation through central obesity. CONCLUSION: There is strong evidence for a causal association between plasma cortisol and greater risk for hypertension, potentially mediated by obesity.

Impaired glucose tolerance and cardiovascular risk factors in relation to infertility: a Mendelian randomization analysis in the Norwegian Mother, Father, and Child Cohort Study.

STUDY QUESTION: Are impaired glucose tolerance (as measured by fasting glucose, glycated hemoglobin, and fasting insulin) and cardiovascular disease risk (as measured by low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, triglycerides, systolic blood pressure, and diastolic blood pressure) causally related to infertility? SUMMARY ANSWER: Genetic instruments suggest that higher fasting insulin may increase infertility in women. WHAT IS KNOWN ALREADY: Observational evidence suggests a shared etiology between impaired glucose tolerance, cardiovascular risk, and fertility problems. STUDY DESIGN, SIZE, DURATION: This study included two-sample Mendelian randomization (MR) analyses, in which we used genome-wide association summary data that were publicly available for the biomarkers of impaired glucose tolerance and cardiovascular disease, and sex-specific genome-wide association studies (GWASs) of infertility conducted in the Norwegian Mother, Father, and Child Cohort Study. PARTICIPANTS/MATERIALS, SETTING, METHODS: There were 68 882 women (average age 30, involved in 81 682 pregnancies) and 47 474 of their male partners (average age 33, 55 744 pregnancies) who had available genotype data and who provided self-reported information on time-to-pregnancy and use of ARTs. Of couples, 12% were infertile (having tried to conceive for ≥12 months or used ARTs to conceive). We applied the inverse variance weighted method with random effects to pool data across variants and a series of sensitivity analyses to explore genetic instrument validity. (We checked the robustness of genetic instruments and the lack of unbalanced horizontal pleiotropy, and we used methods that are robust to population stratification.) Findings were corrected for multiple comparisons by the Bonferroni method (eight exposures: P-value < 0.00625). MAIN RESULTS AND THE ROLE OF CHANCE: In women, increases in genetically determined fasting insulin levels were associated with greater odds of infertility (+1 log(pmol/l): odds ratio 1.60, 95% CI 1.17 to 2.18, P-value = 0.003). The results were robust in the sensitivity analyses exploring the validity of MR assumptions and the role of pleiotropy of other cardiometabolic risk factors. There was also evidence of higher glucose and glycated hemoglobin causing infertility in women, but the findings were imprecise and did not pass our P-value threshold for multiple testing. The results for lipids and blood pressure were close to the null, suggesting that these did not cause infertility. LIMITATIONS, REASONS FOR CAUTION: We did not know if underlying causes of infertility were in the woman, man, or both. Our analyses only involved couples who had conceived. We did not have data on circulating levels of cardiometabolic risk factors, and we opted to conduct an MR analysis using GWAS summary statistics. No sex-specific genetic instruments on cardiometabolic risk factors were available. Our results may be affected by selection and misclassification bias. Finally, the characteristics of our study sample limit the generalizability of our results to populations of non-European ancestry. WIDER IMPLICATIONS OF THE FINDINGS: Treatments for lower fasting insulin levels may reduce the risk of infertility in women. STUDY FUNDING/COMPETING INTEREST(S): The MoBa Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Norwegian Ministry of Education and Research. This work was supported by the European Research Council [grant numbers 947684, 101071773, 293574, 101021566], the Research Council of Norway [grant numbers 262700, 320656, 274611], the South-Eastern Norway Regional Health Authority [grant numbers 2020022, 2021045], and the British Heart Foundation [grant numbers CH/F/20/90003, AA/18/1/34219]. Open Access funding was provided by the Norwegian Institute of Public Health. The funders had no role in the study design; the collection, analysis, and interpretation of data; the writing of the report; or the decision to submit the article for publication. D.A.L. has received research support from National and International government and charitable bodies, Roche Diagnostics and Medtronic for research unrelated to the current work. O.A.A. has been a consultant to HealthLytix. The rest of the authors declare that no competing interests exist. TRIAL REGISTRATION NUMBER: N/A.

Arterial Stiffness and Diabetes Risk in Framingham Heart Study and UK Biobank.

BACKGROUND: Microvascular damage from large artery stiffness (LAS) in pancreatic, hepatic, and skeletal muscles may affect glucose homeostasis. Our goal was to evaluate the association between LAS and the risk of type 2 diabetes using prospectively collected, carefully phenotyped measurements of LAS as well as Mendelian randomization analyses. METHODS: Carotid-femoral pulse wave velocity (CF-PWV) and brachial and central pulse pressure were measured in 5676 participants of the FHS (Framingham Heart Study) without diabetes. We used Cox proportional hazards regression to evaluate the association of CF-PWV and pulse pressure with incident diabetes. We subsequently performed 2-sample Mendelian randomization analyses evaluating the associations of genetically predicted brachial pulse pressure with type 2 diabetes in the UKBB (United Kingdom Biobank). RESULTS: In FHS, individuals with higher CF-PWV were older, more often male, and had higher body mass index and mean arterial pressure compared to those with lower CF-PWV. After a median follow-up of 7 years, CF-PWV and central pulse pressure were associated with an increased risk of new-onset diabetes (per SD increase, multivariable-adjusted CF-PWV hazard ratio, 1.36 [95% CI, 1.03-1.76]; P=0.030; central pulse pressure multivariable-adjusted CF-PWV hazard ratio, 1.26 [95% CI, 1.08-1.48]; P=0.004). In United Kingdom Biobank, genetically predicted brachial pulse pressure was associated with type 2 diabetes, independent of mean arterial pressure (adjusted odds ratio, 1.16 [95% CI, 1.00-1.35]; P=0.049). CONCLUSIONS: Using prospective cohort data coupled with Mendelian randomization analyses, we found evidence supporting that greater LAS is associated with increased risk of developing diabetes. LAS may play an important role in glucose homeostasis and may serve as a useful marker of future diabetes risk.

Relationship Between Ascending Thoracic Aortic Diameter and Blood Pressure: A Mendelian Randomization Study.

BACKGROUND: Observational studies identified elevated blood pressure (BP) as a strong risk factor for thoracic aortic dilation, and BP reduction is the primary medical intervention recommended to prevent progression of aortic aneurysms. However, although BP may impact aortic dilation, aortic size may also impact BP. The causal relationship between BP and thoracic aortic size has not been reliably established. METHODS: Genome-wide association studies summary statistics were obtained for BP and ascending thoracic aortic diameter (AscAoD). Causal effects of BP on AscAoD were estimated using 2-sample Mendelian randomization using a range of pleiotropy-robust methods. RESULTS: Genetically predicted increased systolic BP, diastolic BP, and mean arterial pressure all significantly associate with higher AscAoD (systolic BP: ß estimate, 0.0041 mm/mm Hg [95% CI, 0.0008-0.0074]; P=0.02, diastolic BP: ß estimate, 0.0272 mm/mm Hg [95% CI, 0.0224-0.0320]; P<0.001, and mean arterial pressure: ß estimate, 0.0168 mm/mm Hg [95% CI, 0.0130-0.0206]; P<0.001). Genetically predicted pulse pressure, meanwhile, had an inverse association with AscAoD (ß estimate, -0.0155 mm/mm Hg [95% CI, -0.0213 to -0.0096]; P<0.001). Multivariable Mendelian randomization analyses showed that genetically predicted increased mean arterial pressure and reduced pulse pressure were independently associated with AscAoD. Bidirectional Mendelian randomization demonstrated that genetically predicted AscAoD was inversely associated with pulse pressure (ß estimate, -2.0721 mm Hg/mm [95% CI, -3.1137 to -1.0306]; P<0.001) and systolic BP (ß estimate, -1.2878 mm Hg/mm [95% CI, -2.3533 to -0.2224]; P=0.02), while directly associated with diastolic BP (0.8203 mm Hg/mm [95% CI, 0.2735-1.3672]; P=0.004). CONCLUSIONS: BP likely contributes causally to ascending thoracic aortic dilation. Increased AscAoD likely contributes to lower systolic BP and pulse pressure, but not diastolic BP, consistent with the hemodynamic consequences of a reduced aortic diameter.

A data-adaptive method for investigating effect heterogeneity with high-dimensional covariates in Mendelian randomization.

BACKGROUND: Mendelian randomization is a popular method for causal inference with observational data that uses genetic variants as instrumental variables. Similarly to a randomized trial, a standard Mendelian randomization analysis estimates the population-averaged effect of an exposure on an outcome. Dividing the population into subgroups can reveal effect heterogeneity to inform who would most benefit from intervention on the exposure. However, as covariates are measured post-"randomization", naive stratification typically induces collider bias in stratum-specific estimates. METHOD: We extend a previously proposed stratification method (the "doubly-ranked method") to form strata based on a single covariate, and introduce a data-adaptive random forest method to calculate stratum-specific estimates that are robust to collider bias based on a high-dimensional covariate set. We also propose measures based on the Q statistic to assess heterogeneity between stratum-specific estimates (to understand whether estimates are more variable than expected due to chance alone) and variable importance (to identify the key drivers of effect heterogeneity). RESULT: We show that the effect of body mass index (BMI) on lung function is heterogeneous, depending most strongly on hip circumference and weight. While for most individuals, the predicted effect of increasing BMI on lung function is negative, it is positive for some individuals and strongly negative for others. CONCLUSION: Our data-adaptive approach allows for the exploration of effect heterogeneity in the relationship between an exposure and an outcome within a Mendelian randomization framework. This can yield valuable insights into disease aetiology and help identify specific groups of individuals who would derive the greatest benefit from targeted interventions on the exposure.