Profile‐likelihood Bayesian model averaging for two‐sample summary data Mendelian randomization in the presence of horizontal pleiotropy

Accounting for the impact of rare variants on causal inference with RARE: a novel multivariable Mendelian randomization method

Mendelian randomization (MR) method utilizes genetic variants as instrumental variables to infer the causal effect of an exposure on an outcome. However, the impact of rare variants on traits is often neglected, and traditional MR assumptions can be violated by correlated horizontal pleiotropy (CHP) and uncorrelated horizontal pleiotropy (UHP). To address these issues, we propose a multivariable MR approach, an extension of the standard MR framework: MVMR incorporating Rare variants Accounting for multiple Risk factors and shared horizontal plEiotropy (RARE). In the simulation studies, we demonstrate that RARE effectively detects the causal effects of exposures on outcome with accounting for the impact of rare variants on causal inference. Additionally, we apply RARE to study the effects of high density lipoprotein and low density lipoprotein on type 2 diabetes and coronary atherosclerosis, respectively, thereby illustrating its robustness and effectiveness in real data analysis.

Cathepsins and their role in gynecological cancers: Evidence from two-sample Mendelian randomization analysis

Prior studies have reported connections between cathepsins (CTS) and gynecological cancers; however, the exact causal links are yet to be fully understood. Leveraging publicly accessible genome-wide association study summary datasets, we performed a two-sample bidirectional Mendelian randomization (MR) and multivariate MR (MVMR) analysis, with the inverse variance weighted (IVW) method as the primary approach. MR analysis demonstrated inverse associations between CTSB and cervical cancer (IVW: odds ratio [OR] = 0.9995, 95% confidence interval [CI] = 0.9991-0.9999, P = .0418), CTSE and ovarian cancer (IVW: OR = 0.9197, 95% CI = 0.8505-0.9944, P = .0358), CTSZ and ovarian cancer (IVW: OR = 0.9449, 95% CI = 0.8938-0.9990, P = .0459), CTSE and high grade serous ovarian cancer (IVW: OR = 0.8939, 95% CI = 0.8248-0.9689, P = .0063), and CTSZ and high grade serous ovarian cancer (IVW: OR = 0.9269, 95% CI = 0.8667-0.9913, P = .0268). A positive correlation was identified between CTSH and clear cell ovarian cancer (IVW: OR = 1.1496, 95% CI = 1.0368-1.2745, P = .0081). Nevertheless, subsequent adjustment for the false discovery rate revealed that none of the P-values retained statistical significance (PFDR > 0.05). MVMR analysis results elucidated that CTSZ was inversely associated with cervical cancer (IVW: OR = 0.9988, 95% CI = 0.9981-0.9996, P = .0022). Moreover, a positive association was noted between CTSF and cervical cancer (IVW: OR = 1.0007, 95% CI = 1.0000-1.0014, P = .0364), and similarly, between CTSS and cervical cancer (IVW: OR = 1.0005, 95% CI = 1.0000-1.0011, P = .0490). CTSO exhibited a positive association with non-endometrioid endometrial cancer (IVW: OR = 1.4405, 95% CI = 1.1864-1.7490, P < .001), and CTSH was positively associated with clear cell ovarian cancer (IVW: OR = 1.1167, 95% CI = 1.0131-1.2310, P = .0263). The MVMR analysis findings reveal that CTSZ emerges as a protective element against cervical cancer, whereas CTSF and CTSS represent risk factors for this disease. CTSO stands out as a risk factor for non-endometrioid endometrial cancer, and CTSH acts as a risk factor for clear cell ovarian cancer. This study elucidates causative connections between CTS and gynecological cancers, providing innovative insights for diagnostic and therapeutic optimization.

The genetic landscape of substance use disorders

Substance use disorders represent a significant public health concern with considerable socioeconomic implications worldwide. Twin and family-based studies have long established a heritable component underlying these disorders. In recent years, genome-wide association studies of large, broadly phenotyped samples have identified regions of the genome that harbour genetic risk variants associated with substance use disorders. These regions have enabled the discovery of putative causal genes and improved our understanding of genetic relationships among substance use disorders and other traits. Furthermore, the integration of these data with clinical information has yielded promising insights into how individuals respond to medications, allowing for the development of personalized treatment approaches based on an individual's genetic profile. This review article provides an overview of recent advances in the genetics of substance use disorders and demonstrates how genetic data may be used to reduce the burden of disease and improve public health outcomes.

Mendelian randomization analysis using multiple biomarkers of an underlying common exposure

Mendelian randomization (MR) analysis is increasingly popular for testing the causal effect of exposures on disease outcomes using data from genome-wide association studies. In some settings, the underlying exposure, such as systematic inflammation, may not be directly observable, but measurements can be available on multiple biomarkers or other types of traits that are co-regulated by the exposure. We propose a method for MR analysis on latent exposures (MRLE), which tests the significance for, and the direction of, the effect of a latent exposure by leveraging information from multiple related traits. The method is developed by constructing a set of estimating functions based on the second-order moments of GWAS summary association statistics for the observable traits, under a structural equation model where genetic variants are assumed to have indirect effects through the latent exposure and potentially direct effects on the traits. Simulation studies show that MRLE has well-controlled type I error rates and enhanced power compared to single-trait MR tests under various types of pleiotropy. Applications of MRLE using genetic association statistics across five inflammatory biomarkers (CRP, IL-6, IL-8, TNF-α, and MCP-1) provide evidence for potential causal effects of inflammation on increasing the risk of coronary artery disease, colorectal cancer, and rheumatoid arthritis, while standard MR analysis for individual biomarkers fails to detect consistent evidence for such effects.

Systematic proteome-wide Mendelian randomization using the human plasma proteome to identify therapeutic targets for lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the predominant histological subtype of lung cancer and the leading cause of cancer-related mortality. Identifying effective drug targets is crucial for advancing LUAD treatment strategies.

METHODS

This study employed proteome-wide Mendelian randomization (MR) and colocalization analyses. We collected data on 1394 plasma proteins from a protein quantitative trait loci (pQTL) study involving 4907 individuals. Genetic associations with LUAD were derived from the Transdisciplinary Research in Cancer of the Lung (TRICL) study, including 11,245 cases and 54,619 controls. We integrated pQTL and LUAD genome-wide association studies (GWASs) data to identify candidate proteins. MR utilizes single nucleotide polymorphisms (SNPs) as genetic instruments to estimate the causal effect of exposure on outcome, while Bayesian colocalization analysis determines the probability of shared causal genetic variants between traits. Our study applied these methods to assess causality between plasma proteins and LUAD. Furthermore, we employed a two-step MR to quantify the proportion of risk factors mediated by proteins on LUAD. Finally, protein-protein interaction (PPI) analysis elucidated potential links between proteins and current LUAD medications.

RESULTS

We identified nine plasma proteins significantly associated with LUAD. Increased levels of ALAD, FLT1, ICAM5, and VWC2 exhibited protective effects, with odds ratios of 0.79 (95% CI 0.72-0.87), 0.39 (95% CI 0.28-0.55), 0.91 (95% CI 0.72-0.87), and 0.85 (95% CI 0.79-0.92), respectively. Conversely, MDGA2 (OR, 1.13; 95% CI 1.08-1.19), NTM (OR, 1.12; 95% CI 1.09-1.16), PMM2 (OR, 1.35; 95% CI 1.18-1.53), RNASET2 (OR, 1.15; 95% CI 1.08-1.21), and TFPI (OR, 4.58; 95% CI 3.02-6.94) increased LUAD risk. Notably, none of the nine proteins showed evidence of reverse causality. Bayesian colocalization indicated that RNASET2, TFPI, and VWC2 shared the same variant with LUAD. Furthermore, NTM and FLT1 demonstrated interactions with targets of current LUAD medications. Additionally, FLT1 and TFPI are currently under evaluation as therapeutic targets, while NTM, RNASET2, and VWC2 are potentially druggable. These findings shed light on LUAD pathogenesis, highlighting the tumor-promoting effects of RNASET2, TFPI, and NTM, along with the protective effects of VWC2 and FLT1, providing a significant biological foundation for future LUAD therapeutic targets.

CONCLUSIONS

Our proteome-wide MR analysis highlighted RNASET2, TFPI, VWC2, NTM, and FLT1 as potential drug targets for further clinical investigation in LUAD. However, the specific mechanisms by which these proteins influence LUAD remain elusive. Targeting these proteins in drug development holds the potential for successful clinical trials, providing a pathway to prioritize and reduce costs in LUAD therapeutics.

Association between viral infections and glioma risk: a two-sample bidirectional Mendelian randomization analysis

BACKGROUND

Glioma is one of the leading types of brain tumor, but few etiologic factors of primary glioma have been identified. Previous observational research has shown an association between viral infection and glioma risk. In this study, we used Mendelian randomization (MR) analysis to explore the direction and magnitude of the causal relationship between viral infection and glioma.

METHODS

We conducted a two-sample bidirectional MR analysis using genome-wide association study (GWAS) data. Summary statistics data of glioma were collected from the largest meta-analysis GWAS, involving 12,488 cases and 18,169 controls. Single-nucleotide polymorphisms (SNPs) associated with exposures were used as instrumental variables to estimate the causal relationship between glioma and twelve types of viral infections from corresponding GWAS data. In addition, sensitivity analyses were performed.

RESULTS

After correcting for multiple tests and sensitivity analysis, we detected that genetically predicted herpes zoster (caused by Varicella zoster virus (VZV) infection) significantly decreased risk of low-grade glioma (LGG) development (OR = 0.85, 95% CI: 0.76-0.96, P = 0.01, FDR = 0.04). No causal effects of the other eleven viral infections on glioma and reverse causality were detected.

CONCLUSIONS

This is one of the first and largest studies in this field. We show robust evidence supporting that genetically predicted herpes zoster caused by VZV infection reduces risk of LGG. The findings of our research advance understanding of the etiology of glioma.

Mendelian randomization accounting for complex correlated horizontal pleiotropy while elucidating shared genetic etiology

Mendelian randomization (MR) harnesses genetic variants as instrumental variables (IVs) to study the causal effect of exposure on outcome using summary statistics from genome-wide association studies. Classic MR assumptions are violated when IVs are associated with unmeasured confounders, i.e., when correlated horizontal pleiotropy (CHP) arises. Such confounders could be a shared gene or inter-connected pathways underlying exposure and outcome. We propose MR-CUE (MR with Correlated horizontal pleiotropy Unraveling shared Etiology and confounding), for estimating causal effect while identifying IVs with CHP and accounting for estimation uncertainty. For those IVs, we map their cis-associated genes and enriched pathways to inform shared genetic etiology underlying exposure and outcome. We apply MR-CUE to study the effects of interleukin 6 on multiple traits/diseases and identify several S100 genes involved in shared genetic etiology. We assess the effects of multiple exposures on type 2 diabetes across European and East Asian populations.