TwoStepCisMR: A Novel Method and R Package for Attenuating Bias in cis-Mendelian Randomization Analyses

Associations of lipids and lipid-modifying drug target genes with atrial fibrillation risk based on genomic data

BACKGROUND

The causal associations of lipids and the drug target genes with atrial fibrillation (AF) risk remain obscure. We aimed to investigate the causal associations using genetic evidence.

METHODS

Mendelian randomization (MR) analyses were conducted using summary-level genome-wide association studies (GWASs) in European and East Asian populations. Lipid profiles (low-density lipoprotein cholesterol, triglyceride, and lipoprotein[a]) and lipid-modifying drug target genes (3-hydroxy-3-methylglutaryl-CoA reductase, proprotein convertase subtilisin/kexin type 9, NPC1-like intracellular cholesterol transporter 1, apolipoprotein C3, angiopoietin-like 3, and lipoprotein[a]) were used as exposures. AF was used as an outcome. The inverse variance weighted method was applied as the primary method. Summary-data-based Mendelian randomization analyses were performed for further validation using expression quantitative trait loci data. Mediation analyses were conducted to explore the indirect effect of coronary heart disease.

RESULTS

In the European population, MR analyses demonstrated that elevated levels of lipoprotein(a) increased AF risk. Moreover, analyses focusing on drug targets revealed that the genetically proxied target gene LPA, which simulates the effects of drug intervention by reducing lipoprotein(a), exhibited an association with AF risk. This association was validated in independent datasets. There were no consistent and significant associations observed for other traits when analyzed in different datasets. This finding was also corroborated by Summary-data-based Mendelian randomization analyses between LPA and AF. Mediation analyses revealed that coronary heart disease plays a mediating role in this association. However, in the East Asian population, no statistically significant evidence was observed to support these associations.

CONCLUSIONS

This study provided genetic evidence that Lp(a) may be a causal factor for AF and that LPA may represent a promising pharmacological target for preventing AF in the European population.

Assessing the Role of Cortisol in Anxiety, Major Depression, and Neuroticism: A Mendelian Randomization Study Using SERPINA6/SERPINA1 Variants

Background

Previous evidence informed by the toxic stress model suggests that higher cortisol causes anxiety and major depression, but clinical success is lacking. To clarify the role of cortisol, we used Mendelian randomization to estimate its associations with anxiety, major depression, and neuroticism, leveraging the largest available genome-wide association studies including from the Psychiatric Genomics Consortium, the UK Biobank, and FinnGen.

Methods

After meta-analyzing 2 genome-wide association studies on morning plasma cortisol (n = 32,981), we selected single nucleotide polymorphisms (SNPs) at p < 5 × 10-8 and r2 < 0.3 in the SERPINA6/SERPINA1 gene region encoding proteins that influence cortisol bioavailability. We applied these SNPs to summary genetic associations with the outcomes considered (n = 17,310-449,484), and systolic blood pressure as a positive outcome, using inverse-variance weighted meta-analysis accounting for correlation. Sensitivity analyses addressing SNP correlation and confounding by childhood maltreatment and follow-up analyses using only SNPs that colocalized with SERPINA6 expression were conducted.

Results

Cortisol was associated with anxiety (pooled odds ratio [OR] 1.16 per cortisol z score; 95% CI, 1.04 to 1.31), but not major depression (pooled OR 1.02, 95% CI, 0.95 to 1.10) or neuroticism (β -0.025; 95% CI, -0.071 to 0.022). Sensitivity analyses yielded similar estimates. Cortisol was positively associated with systolic blood pressure, as expected. Using rs9989237 and rs2736898, selected using colocalization, cortisol was associated with anxiety in the UK Biobank (OR 1.32; 95% CI, 1.01 to 1.74) but not with major depression in FinnGen (OR 1.14; 95% CI, 0.95 to 1.37).

Conclusions

Cortisol was associated with anxiety and may be a potential target for prevention. Other targets may be more relevant to major depression and neuroticism.

Association between genetically proxied lipid-lowering drug targets, lipid traits, and amyotrophic lateral sclerosis: a mendelian randomization study

BACKGROUND

The use of circulating lipid traits as biomarkers to predict the risk of amyotrophic lateral sclerosis (ALS) is currently controversial, and the evidence-based medical evidence for the use of lipid-lowering agents, especially statins, on ALS risk remains insufficient. Our aim was to apply a Mendelian randomization (MR) approach to assess the causal impact of lipid-lowering agents and circulating lipid traits on ALS risk.

MATERIALS AND METHODS

Our study included primary and secondary analyses, in which the risk associations of lipid-lowering gene inhibitors, lipid traits, and ALS were assessed by the inverse variance weighting method as the primary approach. The robustness of the results was assessed using LDSC assessment, conventional MR sensitivity analysis, and used Mediating MR to explore potential mechanisms of occurrence. In the secondary analysis, the association of lipid-lowering genes with ALS was validated using the Summary data-based Mendelian Randomization (SMR) method.

RESULTS

Our results showed strong evidence between genetic proxies for Apolipoprotein B (ApoB) inhibitor (OR = 0.76, 95% CI = 0.68 - 0.86; P = 5.58 × 10-6) and reduced risk of ALS. Additionally, 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) inhibitor (OR = 1.06, 95% CI = 0. 85-1.33) was not found to increase ALS risk. SMR results suggested that ApoB expression was associated with increased ALS risk, and colocalization analysis did not support a significant common genetic variation between ApoB and ALS. Mediator MR analysis suggested a possible mediating role for interleukin-6 and low-density lipoprotein cholesterol (LDL-C). While elevated LDL-C was significantly associated with increased risk of ALS among lipid traits, total cholesterol (TC) and ApoB were weakly associated with ALS. LDSC results suggested a potential genetic correlation between these lipid traits and ALS.

CONCLUSIONS

Using ApoB inhibitor can lower the risk of ALS, statins do not trigger ALS, and LDL-C, TC, and ApoB levels can predict the risk of ALS.

Unraveling the link between plasma caffeine concentrations and inflammatory bowel disease risk through Mendelian randomization

BACKGROUND

Caffeine is believed to possess anti-inflammatory properties, yet direct population-based evidence regarding its impact on inflammatory bowel disease (IBD) remains scarce.

OBJECTIVES

In this study, we used 2-sample Mendelian randomization (MR) study to investigate the causal relationship between long-term plasma caffeine concentrations and IBD and its subtypes, ulcerative colitis (UC) and Crohn disease (CD).

METHODS

We used single nucleotide polymorphisms (SNPs) associated with plasma caffeine concentrations at genome-wide significance within a ±100-kb range around the CYP1A2 or AHR genes as instrumental variables. Genome-wide association study (GWAS) data for IBD and its subtypes were obtained from FinnGen and International Inflammatory Bowel Disease Genetics Consortium. We conducted a meta-analysis of MR-related SNPs from both sources and used a multiplicative inverse variance-weighted random effects model to combine the effects of each SNP proxy on exposure to outcomes.

RESULTS

In our study, genetically predicted higher plasma caffeine concentrations were associated with a lower risk of IBD, with an odds ratio (OR) of 0.78 (95% confidence interval [CI]: 0.66, 0.91; PFDR = 0.004). This trend was also observed in UC and CD, with ORs of 0.79 (95% CI: 0.66, 0.94; PFDR = 0.014) and 0.78 (95% CI: 0.62, 0.98; PFDR = 0.032), respectively.

CONCLUSION

Our study indicates a potential causal link between genetically predicted higher plasma caffeine concentrations and a reduced risk of IBD, including its subtypes UC and CD.

Genetic association of lipid-lowering drug target genes with erectile dysfunction and male reproductive health

Objective

The effect of hypolipidemic drugs on male erectile function is still controversial. This Mendelian randomization (MR) study aimed to explore the potential impact of lipid-lowering drug targets on ED.

Methods

We collected seven genetic variants encoding lipid-lowering drug targets (LDLR, HMGCR, NPC1L1, PCSK9, APOB, APOC3 and LPL) from published genome-wide association study (GWAS) statistics, and performed drug target MR analysis. The risk of ED was defined as the primary outcome, sex hormone levels and other diseases as the secondary outcomes. Mediation analyses were performed to explore potential mediating factors.

Results

The results showed that LDLR, LPL agonists and APOC3 inhibitors were significantly associated with a reduced risk of ED occurrence. APOB inhibitors were associated with an increased risk of ED occurrence. In terms of sex hormone levels, LDLR and LPL agonists were significantly associated with increased TT levels, and HMGCR was associated with decreased TT and BT levels significantly. In terms of male-related disease, MR results showed that LDLR agonists and PCSK9 inhibitors were significantly associated with an elevated risk of PH; HMGCR, NPC1L1 inhibitors were associated with a reduced risk of PCa; and LDLR agonists were significantly associated with a reduced risk of AS and MI; in addition, HMGCR inhibitors were associated with a reduced risk of PCa.

Conclusion

After performing drug-targeted MR analysis, we found that that there was a causal relationship between lipid-lowering drug targets and ED. APOC3, APOB, LDLR and LPL may be new candidate drug targets for the treatment of ED.

A drug target for erectile dysfunction to help improve fertility, sexual activity, and wellbeing: mendelian randomisation study

OBJECTIVE

To investigate the association of genetically proxied (using a surrogate biomarker) inhibition of phosphodiesterase 5 (PDE5), an established drug target for erectile dysfunction, with fertility, sexual behaviour, and subjective wellbeing.

DESIGN

Two sample cis-mendelian randomisation study.

SETTING

Summary data on genetic associations obtained from the International Consortium for Blood Pressure and UK Biobank.

PARTICIPANTS

Individuals of European ancestry from the International Consortium for Blood Pressure (n=757 601) for estimating PDE5 inhibition (using the surrogate biomarker of diastolic blood pressure reduction), and UK Biobank (n=211 840) for estimating the fertility, sexual behaviour, and subjective wellbeing outcomes in male participants.

INTERVENTION

Genetically proxied PDE5 inhibition.

MAIN OUTCOME MEASURES

Number of children fathered, number of sexual partners, probability of never having had sexual intercourse, and subjective wellbeing.

RESULTS

Genetically proxied PDE5 inhibition was associated with male participants having 0.28 (95% confidence interval 0.16 to 0.39) more children (false discovery rate corrected P<0.001). This association was not identified in female participants. No evidence was found of an association between genetically proxied PDE5 inhibition and number of sexual partners, probability of never having had sexual intercourse, or self-reported wellbeing in male participants.

CONCLUSIONS

The findings of this study provide genetic support for PDE5 inhibition potentially increasing the number of children fathered by male individuals. Absence of this association in female participants supports increased propensity for sustained and robust penile erections as a potential underlying mechanism. Further studies are required to confirm this, however, and these findings should not promote indiscriminate use of PDE5 inhibitors, which can also have harmful adverse effects.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Association between testosterone and cancers risk in women: a two-sample Mendelian randomization study

OBJECTIVE

Previous observational studies have explored the correlation between testosterone and cancer risk. However, the causal association between testosterone and various cancer types in women remains inconclusive. The objective of this Mendelian randomization study is to evaluate the causal links between total testosterone (TT) and bioavailable testosterone (BT) with cancer risk in females.

METHODS

Initially, a rigorous quality control process was employed to identify suitable instrumental single nucleotide polymorphisms (SNPs) associated with the exposure under investigation that exhibited a significant association. The genetic causal relationship between female testosterone levels and the risk of developing cancers was examined through a two-sample Mendelian randomization. Various analytical methods, including inverse-variance weighted (IVW), MR-Egger, weighted median, simple mode, and weighted mode, were applied in the investigation. Key findings were primarily based on the results obtained via IVW (random effects), and sensitivity analyses were conducted to assess the reliability of the obtained results. Furthermore, maximum likelihood, penalized weighted median, and IVW (fixed effects) methods were utilized to further validate the robustness of the results.

RESULTS

Based on the results of IVW analysis, our study indicated a positive causal relationship between BT and breast cancer (OR = 1.1407, 95%CI: 1.0627-1.2244, P = 0.0015) and endometrial cancer (OR = 1.4610, 95%CI: 1.2695-1.6813, P = 1.22E-06). Moreover, our findings also showed a positive causal association between TT and breast cancer (OR = 1.1764, 95%CI: 1.0846-1.2761, P = 0.0005), cervical cancer(OR = 1.0020, 95%CI: 1.0007-1.0032, P = 0.0077), and endometrial cancer(OR = 1.4124, 95%CI: 1.2083-1.6511, P = 0.0001). Additionally, our results demonstrated a negative causal relationship between BT and ovarian cancer (OR = 0.8649, 95%CI: 0.7750-0.9653, P = 0.0320). However, no causal relationship was found between BT, TT and other types of cancer (corrected P > 0.05).

CONCLUSIONS

This study elucidates the role of testosterone on the development of breast cancer, endometrial cancer, ovarian cancer, and cervical cancer. It also hints at a potential but fragile link between testosterone and bladder cancer, as well as thyroid cancer. Nonetheless, it's worth noting that no statistically significant relationship between testosterone and various other types of cancer in females was identified.

Association between genetically proxied HMGCR inhibition and male reproductive health: A Mendelian randomization study

BACKGROUND

The causal associations between statin use and male sex hormone levels and related disorders have not been fully understood. In this study, we employed Mendelian randomization for the first time to investigate these associations.

METHODS

In two-sample Mendelian randomization framework, genetic proxies for hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) inhibition were identified as variants in the HMGCR gene that were associated with both levels of gene expression and low density lipoprotein cholesterol (LDL-C). We assessed the causal relationship between HMGCR inhibitor and 5 sex hormone levels/2 male-related diseases. Additionally, we investigated the association between 4 circulating lipid traits and outcomes. The "inverse variance weighting" method was used as the primary approach, and we assessed for potential heterogeneity and pleiotropy. In a secondary analysis, we revalidated the impact of HMGCR on 7 major outcomes using the summary-data-based Mendelian randomization method.

RESULTS

Our study found a significant causal association between genetic proxies for HMGCR inhibitor and decreased levels of total testosterone (TT) (LDL-C per standard deviation = 38.7mg/dL, effect = -0.20, 95% confidence interval [CI] = -0.25 to -0.15) and bioavailable testosterone (BT) (effect = -0.15, 95% CI = -0.21 to -0.10). Obesity-related factors were found to mediate this association. Furthermore, the inhibitor were found to mediate a reduced risk of prostate cancer (odds ratio = 0.81, 95%CI = 0.7-0.93) by lowering bioavailable testosterone levels, without increasing the risk of erectile dysfunction (P = .17). On the other hand, there was a causal association between increased levels of LDL-C, total cholesterol, triglycerides (TG) and decreased levels of TT, sex hormone-binding globulin, and estradiol.

CONCLUSIONS

The use of HMGCR inhibitor will reduce testosterone levels and the risk of prostate cancer without the side effect of erectile dysfunction. LDL-C, total cholesterol, and TG levels were protective factors for TT, sex hormone-binding globulin, and estradiol.

Genetic association of lipids and lipid-lowering drug target genes with non-alcoholic fatty liver disease

BACKGROUND

Some observational studies found that dyslipidaemia is a risk factor for non-alcoholic fatty liver disease (NAFLD), and lipid-lowering drugs may lower NAFLD risk. However, it remains unclear whether dyslipidaemia is causative for NAFLD. This Mendelian randomisation (MR) study aimed to explore the causal role of lipid traits in NAFLD and evaluate the potential effect of lipid-lowering drug targets on NAFLD.

METHODS

Genetic variants associated with lipid traits and variants of genes encoding lipid-lowering drug targets were extracted from the Global Lipids Genetics Consortium genome-wide association study (GWAS). Summary statistics for NAFLD were obtained from two independent GWAS datasets. Lipid-lowering drug targets that reached significance were further tested using expression quantitative trait loci data in relevant tissues. Colocalisation and mediation analyses were performed to validate the robustness of the results and explore potential mediators.

FINDINGS

No significant effect of lipid traits and eight lipid-lowering drug targets on NAFLD risk was found. Genetic mimicry of lipoprotein lipase (LPL) enhancement was associated with lower NAFLD risks in two independent datasets (OR₁ = 0.60 [95% CI 0.50-0.72], p₁ = 2.07 × 10^-8; OR₂ = 0.57 [95% CI 0.39-0.82], p₂ = 3.00 × 10^-3). A significant MR association (OR = 0.71 [95% CI, 0.58-0.87], p = 1.20 × 10^-3) and strong colocalisation association (PP.H₄ = 0.85) with NAFLD were observed for LPL expression in subcutaneous adipose tissue. Fasting insulin and type 2 diabetes mediated 7.40% and 9.15%, respectively, of the total effect of LPL on NAFLD risk.

INTERPRETATION

Our findings do not support dyslipidaemia as a causal factor for NAFLD. Among nine lipid-lowering drug targets, LPL is a promising candidate drug target in NAFLD. The mechanism of action of LPL in NAFLD may be independent of its lipid-lowering effects.

FUNDING

Capital's Funds for Health Improvement and Research (2022-4-4037). CAMS Innovation Fund for Medical Sciences (CIFMS, grant number: 2021-I2M-C&T-A-010).

Proteome-wide Mendelian randomization implicates nephronectin as an actionable mediator of the effect of obesity on COVID-19 severity

Obesity is a major risk factor for Coronavirus disease (COVID-19) severity; however, the mechanisms underlying this relationship are not fully understood. As obesity influences the plasma proteome, we sought to identify circulating proteins mediating the effects of obesity on COVID-19 severity in humans. Here, we screened 4,907 plasma proteins to identify proteins influenced by body mass index using Mendelian randomization. This yielded 1,216 proteins, whose effect on COVID-19 severity was assessed, again using Mendelian randomization. We found that an s.d. increase in nephronectin (NPNT) was associated with increased odds of critically ill COVID-19 (OR = 1.71, P = 1.63 × 10^-10). The effect was driven by an NPNT splice isoform. Mediation analyses supported NPNT as a mediator. In single-cell RNA-sequencing, NPNT was expressed in alveolar cells and fibroblasts of the lung in individuals who died of COVID-19. Finally, decreasing body fat mass and increasing fat-free mass were found to lower NPNT levels. These findings provide actionable insights into how obesity influences COVID-19 severity.

Appraisal of the causal effect of plasma caffeine on adiposity, type 2 diabetes, and cardiovascular disease: two sample mendelian randomisation study

Objective

To investigate the potential causal effects of long term plasma caffeine concentrations on adiposity, type 2 diabetes, and major cardiovascular diseases.

Design

Two sample mendelian randomisation study.

Setting

Genome-wide association study summary data for associations of two single nucleotide polymorphisms associated with plasma caffeine at the genome-wide significance threshold (rs2472297 near the CYP1A2 gene and rs4410790 near the AHR gene) and their association with the outcomes.

Participants

Primarily individuals of European ancestry participating in cohorts contributing to genome-wide association study consortia.

Main outcome measures

Outcomes studied were body mass index, whole body fat mass, whole body fat-free mass, type 2 diabetes, ischaemic heart disease, atrial fibrillation, heart failure, and stroke.

Results

Higher genetically predicted plasma caffeine concentrations were associated with lower body mass index (beta -0.08 standard deviation (SD) (95% confidence interval -0.10 to -0.06), where 1 SD equals about 4.8 kg/m² in body mass index, for every standard deviation increase in plasma caffeine) and whole body fat mass (beta -0.06 SD (-0.08 to -0.04), 1 SD equals about 9.5 kg; P<0.001) but not fat-free mass (beta -0.01 SD (-0.02 to -0.00), 1 SD equals about 11.5 kg; P=0.17). Higher genetically predicted plasma caffeine concentrations were associated with a lower risk of type 2 diabetes in two consortia (FinnGen and DIAMANTE), with a combined odds ratio of 0.81 ((95% confidence interval 0.74 to 0.89); P<0.001). Approximately half (43%; 95% confidence interval 30% to 61%) of the effect of caffeine on type 2 diabetes was estimated to be mediated through body mass index reduction. No strong associations were reported between genetically predicted plasma caffeine concentrations and a risk of any of the studied cardiovascular diseases.

Conclusions

Higher plasma caffeine concentrations might reduce adiposity and risk of type 2 diabetes. Further clinical study is warranted to investigate the translational potential of these findings towards reducing the burden of metabolic disease.