A practical problem with Egger regression in Mendelian randomization

Causal relationship between serum metabolites and idiopathic pulmonary fibrosis: Insights from a two-sample Mendelian randomization study

Background

Idiopathic pulmonary fibrosis (IPF) is an irreversible lung disease with unclear pathological mechanisms. In this study, we utilized bidirectional Mendelian randomization (MR) to analyze the relationship between serum metabolites and IPF, and conducted metabolic pathway analysis.

Aim

To determine the causal relationship between serum metabolites and IPF using MR analysis.

Methods

A two-sample MR analysis was conducted to evaluate the causal relationship between 824 serum metabolites and IPF. The inverse variance weighted (IVW) method was used to estimate the causal relationship between exposure and results. Sensitivity analysis was conducted using MR Egger, weighted median, and maximum likelihood to eliminate pleiotropy. Additionally, metabolic pathway analysis was conducted to identify potential metabolic pathways.

Results

We identified 12 serum metabolites (6 risks and 6 protective) associated with IPF from 824 metabolites. Among them, 11 were known and 1 was unknown. 1-Eicosatrienoylglycophorophospholine and 1-myristoylglycophorophospholine were bidirectional MR positive factors, with 1-myristoylglycophorophospholine being a risk factor (1.0013, 1.0097) and 1-eicosatrienoylglycophorine being a protective factor (0.9914, 0.9990). The four lipids (1-linoleoylglycerophoethanolamine*, total cholesterol in large high-density lipoprotein [HDL], cholesterol esters in very large HDL, and phospholipids in very large HDL) and one NA metabolite (degree of unsaturation) were included in the known hazardous metabolites. The known protective metabolites included three types of lipids (carnitine, 1-linoleoylglycerophoethanolamine*, and 1-eicosatrienoylglycerophophophorine), one amino acid (hypoxanthine), and two unknown metabolites (the ratio of omega-6 fatty acids to omega-3 fatty acids, and the ratio of photoshopids to total lipids ratio in chylomicrons and extremely large very low-density lipoprotein [VLDL]). Moreover, sn-Glycerol 3-phosphate and 1-Acyl-sn-glycero-3-phosphocline were found to be involved in the pathogenesis of IPF through metabolic pathways such as Glycerolide metabolism and Glycerophospholipid metabolism.

Conclusion

Our study identified 6 causal risks and 6 protective serum metabolites associated with IPF. Additionally, 2 metabolites were found to be involved in the pathogenesis of IPF through metabolic pathways, providing a new perspective for further understanding the metabolic pathway and the pathogenesis of IPF.

The causality between use of glucocorticoids and risk of pancreatitis: a Mendelian randomization study

Background and aim

To date, the association between glucocorticoid use and the risk of pancreatitis remains controversial. The aim of this study was the investigation of this possible relationship.

Methods

We carried out a two-sample Mendelian randomization (MR) analysis using GWAS data from European ancestry, East Asian descendants and the FinnGen Biobank Consortium to evaluate this potential causal relationship. Genetic variants associated with glucocorticoid use were selected based on genome-wide significance (p < 5×10-8).

Results

Our MR analysis of European ancestry data revealed no significant causal relationship between glucocorticoid use and AP (IVW: OR=1.084, 95% CI= 0.945-1.242, P=0.249; MR-Egger: OR=1.049, 95% CI= 0.686-1.603, P=0.828; weighted median: OR=1.026, 95% CI= 0.863-1.219, P=0.775) or CP (IVW: OR=1.027, 95% CI= 0.850-1.240, P=0.785; MR-Egger: OR= 1.625, 95% CI= 0.913-2.890, P= 0.111; weighted median: OR= 1.176, 95% CI= 0.909-1.523, P= 0.218). Sensitivity analyses, including MR-Egger and MR-PRESSO, indicated no evidence of pleiotropy or heterogeneity, confirming the robustness of our findings. Multivariable MR analysis adjusted for alcohol consumption, BMI, cholelithiasis and C-reactive protein levels supported these findings. Replicated analysis was performed on datasets from the FinnGen Biobank Consortium and East Asian descendants, and similar results were obtained.

Conclusions

This MR analysis suggests that there is no causal association between glucocorticoid use and the risk of pancreatitis.

A robust cis-Mendelian randomization method with application to drug target discovery

Mendelian randomization (MR) uses genetic variants as instrumental variables (IVs) to investigate causal relationships between traits. Unlike conventional MR, cis-MR focuses on a single genomic region using only cis-SNPs. For example, using cis-pQTLs for a protein as exposure for a disease opens a cost-effective path for drug target discovery. However, few methods effectively handle pleiotropy and linkage disequilibrium (LD) of cis-SNPs. Here, we propose cisMR-cML, a method based on constrained maximum likelihood, robust to IV assumption violations with strong theoretical support. We further clarify the severe but largely neglected consequences of the current practice of modeling marginal, instead of conditional genetic effects, and only using exposure-associated SNPs in cis-MR analysis. Numerical studies demonstrated our method's superiority over other existing methods. In a drug-target analysis for coronary artery disease (CAD), including a proteome-wide application, we identified three potential drug targets, PCSK9, COLEC11 and FGFR1 for CAD.

Evaluation of bi-directional causal association between obstructive sleep apnoea syndrome and diabetic microangiopathy: a Mendelian randomization study

Background

The relationship between obstructive sleep apnea syndrome (OSAS) and diabetic microangiopathy remains controversial.

Objective

This study aimed to use bidirectional two-sample Mendelian Randomization (MR) to assess the causal relationship between OSAS and diabetic microangiopathy.

Methods

First, we used the Linkage Disequilibrium Score Regression(LDSC) analysis to assess the genetic correlation. Then, the bidirectional two-sample MR study was conducted in two stages: OSAS and lung function-related indicators (forced vital capacity (FVC) and forced expiratory volume in 1 s (FEV1)) were investigated as exposures, with diabetic microangiopathy as the outcome in the first stage, and genetic tools were used as proxy variables for OSAS and lung function-related measures in the second step. Genome-wide association study data came from the open GWAS database. We used Inverse-Variance Weighted (IVW), MR-Egger regression, Weighted median, Simple mode, and Weighted mode for effect estimation and pleiotropy testing. We also performed sensitivity analyses to test the robustness of the results. Furthermore, we performed multivariate and mediation MR analyses.

Results

In the LDSC analysis, We found a genetic correlation between OSAS, FVC, FEV 1, and diabetic microangiopathy. In the MR analysis, based on IVW analysis, genetically predicted OSAS was positively correlated with the incidence of diabetic retinopathy (DR), diabetic kidney disease (DKD), and diabetic neuropathy (DN). In the subgroup analysis of DR, there was a significant causal relationship between OSAS and background diabetic retinopathy (BDR) and proliferative diabetic retinopathy (PDR). The reverse MR did not show a correlation between the incidence of diabetic microangiopathy and OSAS. Reduced FVC had a potential causal relationship with increased incidence of DR and PDR. Reduced FEV1 had a potential causal relationship with the increased incidence of BDR, PDR, and DKD. Multivariate MR analysis showed that the association between OSAS and diabetic microangiopathy remained significant after adjusting for confounding factors. However, we did not find the significant mediating factors.

Conclusion

Our results suggest that OSAS may be a cause of the development of diabetic microangiopathy, and OSAS may also be associated with a high risk of diabetic microangiopathy, providing a reference for a better understanding of the prevention of diabetic microangiopathy.

Assessment of causal associations between obesity and peripheral artery disease: a bidirectional Mendelian randomization study

Background

Several observational studies have documented a potential link between obesity and peripheral artery disease (PAD), although conflicting findings exist. The causal relationship between obesity and PAD continues to be a subject of ongoing debate in the medical community.

Objectives

In this study, we employed a bidirectional Mendelian randomization (MR) analysis to explore the potential causal relationship between obesity and the risk of PAD.

Methods

To investigate these causal relationships, we conducted bidirectional MR analysis using publicly available genome-wide association study (GWAS) data. Effect estimates were calculated using the random-effects inverse variance-weighted (IVW) method.

Results

We identified eight independent single nucleotide polymorphisms (SNPs) associated with obesity in 218,735 samples involving 16,380,465 SNPs, all of which met the genome-wide significance threshold (p < 5 × 10-⁸). The IVW analysis indicates a significant positive association between genetic obesity and multiple datasets with PAD as the outcome: Queue-1 (GWAS ID: finn-b-I9_PAD) (OR = 1.138, 95% CI: 1.027-1.261, p = 0.013), Queue-2 (GWAS ID: bbj-a-144) (OR = 1.190, 95% CI: 1.019-1.390, p = 0.028), Queue-3 (GWAS ID: ebi-a-GCST90018670) (OR = 1.174, 95% CI: 1.014-1.360, p = 0.032), and Queue-4 (GWAS ID: ebi-a-GCST90018890) (OR = 1.194, 95% CI: 1.099-1.296, p < 0.001). However, we did not observe a significant genetic-level association between obesity and PAD for Queue-5 (GWAS ID: ukb-d-I9_PAD) (OR = 1.001, 95% CI: 1.000-1.002, p = 0.071). Furthermore, we conducted a reverse causal MR analysis to explore the potential reverse causal relationship between obesity and PAD. This comprehensive analysis did not provide evidence of a reverse causal association between these two factors.

Conclusions

In summary, our study offers genetic evidence suggesting a possible causal link between obesity and PAD. While we did not find evidence supporting the "obesity paradox", prudent weight management remains crucial, as lower weight does not necessarily guarantee better outcomes. As with any study, caution is required in interpreting the findings. Further research is essential to assess the clinical relevance of weight in preventing PAD, which could inform the development of more precise intervention strategies.

Mendelian randomisation and mediation analysis of self-reported walking pace and coronary artery disease

The aim of this study was to assess the causal relationship between habitual walking pace and cardiovascular disease risk using a Mendelian randomisation approach. We performed both one- and two-sample Mendelian randomisation analyses in a sample of 340,000 European ancestry participants from UK Biobank, applying a range of sensitivity analyses to assess pleiotropy and reverse causality. We used a latent variable framework throughout to model walking pace as a continuous exposure, despite being measured in discrete categories, which provided more robust and interpretable causal effect estimates. Using one-sample Mendelian randomisation, we estimated that a 1 mph (i.e., 1.6 kph) increase in self-reported habitual walking pace corresponds to a 63% (hazard ratio (HR) = 0.37, 95% confidence interval (CI), 0.25-0.55, P = 2.0 × 10-6) reduction in coronary artery disease risk. Using conditional analyses, we also estimated that the proportion of the total effect on coronary artery disease mediated through BMI was 45% (95% CI 16-70%). We further validated findings from UK Biobank using two-sample Mendelian randomisation with outcome data from the CARDIoGRAMplusC4D consortium. Our findings suggest that interventions that seek to encourage individuals to walk more briskly should lead to protective effects on cardiovascular disease risk.

Collider bias correction for multiple covariates in GWAS using robust multivariable Mendelian randomization

Genome-wide association studies (GWAS) have identified many genetic loci associated with complex traits and diseases in the past 20 years. Multiple heritable covariates may be added into GWAS regression models to estimate direct effects of genetic variants on a focal trait, or to improve the power by accounting for environmental effects and other sources of trait variations. When one or more covariates are causally affected by both genetic variants and hidden confounders, adjusting for them in GWAS will produce biased estimation of SNP effects, known as collider bias. Several approaches have been developed to correct collider bias through estimating the bias by Mendelian randomization (MR). However, these methods work for only one covariate, some of which utilize MR methods with relatively strong assumptions, both of which may not hold in practice. In this paper, we extend the bias-correction approaches in two aspects: first we derive an analytical expression for the collider bias in the presence of multiple covariates, then we propose estimating the bias using a robust multivariable MR (MVMR) method based on constrained maximum likelihood (called MVMR-cML), allowing the presence of invalid instrumental variables (IVs) and correlated pleiotropy. We also established the estimation consistency and asymptotic normality of the new bias-corrected estimator. We conducted simulations to show that all methods mitigated collider bias under various scenarios. In real data analyses, we applied the methods to two GWAS examples, the first a GWAS of waist-hip ratio with adjustment for only one covariate, body-mass index (BMI), and the second a GWAS of BMI adjusting metabolomic principle components as multiple covariates, illustrating the effectiveness of bias correction.

Positive association between omega-3/6 polyunsaturated fatty acids and idiopathic normal pressure hydrocephalus: a mendelian randomization study

Background: Idiopathic normal pressure hydrocephalus (iNPH) is a common disease among the elderly, which brings great harm to the health of patients and imposes a huge economic burden on the healthcare system. Research has shown that it is possible to prevent iNPH through nutritional and dietary interventions. Intake of omega-3 and omega-6 polyunsaturated fatty acids (PUFAs) can reduce the risk of many diseases. In this study, we aimed to explore the association between omega-3/6 PUFAs and iNPH. Methods: We conducted a two-sample Mendelian randomization (MR) study using summary data from publicly available genome-wide association studies (GWAS) to evaluate the potential impact of omega-3 and omega-6 PUFAs on the risk of iNPH in European populations. Inverse variance weighting was used as the main method for MR analysis, with Wald ratio, weighted median, MR-Egger, simple mode, and weighted mode as supplementary methods. In addition, we performed a series of instrument variable strength evaluations and sensitivity analyses to test the reliability of the study results. Finally, we also conducted a linkage disequilibrium score regression (LDSC) analysis to assess the genetic correlation and distinguish between causal associations and shared genetic variants between PUFAs and iNPH. Results: One SD increase in genetically predicted levels of total omega-3 PUFAs (OR: 0.748; 95% CI: 0.597-0.937; p = 0.012; IVW), DHA (OR: 0.709; 95% CI: 0.532-0.945; p = 0.019; IVW), ALA (OR: 0.001; 95% CI: 1.17E-06-0.423; p = 0.026; Wald ratio), and DHA (OR: 0.709; 95% CI: 0.532-0.945; p = 0.019; IVW) were associated with a reduced iNPH risk. LDSC did not reveal any significant genetic correlations. Conclusion: Higher genetically predicted levels of total omega-3 PUFAs, ALA, DHA, and DPA are associated with a reduced risk of iNPH. These findings have important implications for preventing iNPH and future nutritional guidance.

No evident causal association between Helicobacter pylori infection and colorectal cancer: a bidirectional mendelian randomization study

Observational studies have reported a correlation between Helicobacter pylori infection and colorectal cancer (CRC); however, the underlying cause has remained unclear. This research was aimed at determining whether there is a correlation between H. pylori infection and CRC by measuring the prevalence of H. pylori CagA antibodies and VacA antibodies. Using data from many genome-wide association studies (GWAS), we conducted a Mendelian randomization (MR) study with two sample GWAS. Then, we used bidirectional MR to evaluate the association between H. pylori infection and CRC for identifying causation. The most common method of analysis was the inverse variance-weighted technique. In addition, we performed supplementary analyses using the weighted median technique and MR-Egger regression. Horizontal pleiotropic outliers were identified and corrected using the MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method. Genetically predicted anti-H. pylori IgG seropositivity was not causally associated with CRC [odds ratio (OR): 1.12; 95% confidence interval (CI): 0.98-1.27, P = 0.08] and neither were H. pylori VacA antibody levels (OR = 0.96, 95% CI: 0.90-1.02, P = 0.25) or H. pylori CagA antibody levels (OR = 1.00, 95% CI: 0.93-1.07, P = 0.92). Furthermore, reverse MR analysis did not reveal evidence for a causal effect of CRC on H. pylori infection. The weighted median, the MR-Egger method, and MR-PRESSO yielded identical results. Using genetic data, MR analysis showed there was no evidence for a causal association between seroprevalence of H. pylori infection and CRC. The relationship between H. pylori infection and CRC requires further research.

Significant association between high neutrophil-lymphocyte ratio and poor prognosis in patients with hepatocellular carcinoma: a systematic review and meta-analysis

Objective

Whether neutrophil-lymphocyte ratio (NLR) is an applicative predictor of poor prognosis in patients with hepatocellular carcinoma (HCC) remains controversial. In response to the current conflicting data, this meta-analysis was conducted to gain a comprehensive and systematic understanding of prognostic value of NLR in HCC.

Methods

Several English databases, including PubMed, EMBASE, and the Cochrane Library, with an update date of February 25, 2023, were systematically searched. We set the inclusion criteria to include randomized controlled trial (RCT) studies that reported the prognostic value of serum NLR levels in patients with HCC receiving treatment. Both the combined ratio (OR) and the diagnosis ratio (DOR) were used to assess the prognostic performance of NLR. Additionally, we completed the risk of bias assessment by Cochrane Risk of Bias Assessment Tool.

Results

This meta-analysis ultimately included 16 studies with a total of 4654 patients with HCC. The results showed that high baseline NLR was significantly associated with poor prognosis or recurrence of HCC. The sensitivity of 0.67 (95% confidence interval [CI]. 0.59-0.73); specificity of 0.723 (95% CI: 0.64-0.78) and DOR of 5.0 (95% CI: 4.0-7.0) were pooled estimated from patient-based analyses. Subsequently, the combined positive likelihood ratio (PLR) and negative likelihood ratio (NLHR) were calculated with the results of 2.4 (95% CI: 1.9-3.0) and 0.46 (95% CI: 0.39-0.56), respectively. In addition, area under the curve (AUC) of the summary receiver operating characteristic (SROC) reflecting prognostic accuracy was calculated to be 0.75 (95% CI: 0.71-0.78). The results of subgroup analysis suggested that high NLR was an effective predictive factor of poor prognosis in HCC in mainland China as well as in the northern region.

Conclusion

Our findings suggest that high baseline NLR is an excellent predictor of poor prognosis or relapse in patients with HCC, especially those from high-incidence East Asian populations.

Systematic review registration

https://www.crd.york.ac.uk/prospero/#recordDetails, identifier CRD42023440640.

Association of the gut microbiota with coronary artery disease and myocardial infarction: A Mendelian randomization study

Background: Previous studies have indicated that the gut microbiota (GM) is associated with coronary artery disease (CAD), but the causality of these associations remains unestablished due to confounding factors and reverse causality. We conducted Mendelian randomization study (MR) to determine the causal effect of the specific bacterial taxa on CAD/myocardial infarction (MI) and identify the mediating factors involved. Methods: Two-sample MR, multivariable MR (MVMR) and mediation analysis were performed. Inverse-variance weighting (IVW) was the main method used to analyze causality, and sensitivity analysis was used to verify the reliability of the study. Causal estimates from CARDIoGRAMplusC4D and FinnGen databases were combined using the meta-analysis method, and repeated validation was conducted based on the UK Biobank (UKB) database. Confounders that may affect the causal estimates were corrected by MVMP and the potential mediation effects were investigated by using mediation analysis. Results: The study suggested that increased abundance of the RuminococcusUCG010 genus leads to a lower risk of CAD (OR, 0.88; 95% CI, 0.78, 1.00; p = 2.88 × 10^-2) and MI (OR, 0.88; 95% CI, 0.79, 0.97; p = 1.08 × 10^-2), with consistent results in both meta-analysis (CAD: OR, 0.86; 95% CI, 0.78, 0.96; p = 4.71 × 10^-3; MI: OR, 0.82; 95% CI, 0.73, 0.92; p = 8.25 × 10^-4) and repeated analysis of the UKB dataset (CAD: OR, 0.99; 95% CI, 0.99, 1.00, p = 2.53 × 10^-4; MI: OR, 0.99; 95% CI, 0.99, 1.00, p = 1.85 × 10-11). Based on multiple databases, T2DM was proved as a mediating factor in the causal effect of RuminococcusUCG010 and CAD/MI, with an average mediation effect proportion of 20% on CAD and 17% on MI, respectively. Conclusion: This MR study provided suggestive genetic evidence that the higher the RuminococcusUCG010 abundance is, the lower the risk of CAD and MI, with T2DM playing a mediating effect. This genus may become a novel target in strategies for treating and preventing CAD and MI.