MR-LDP: a two-sample Mendelian randomization for GWAS summary statistics accounting for linkage disequilibrium and horizontal pleiotropy

The goldmine of GWAS summary statistics: a systematic review of methods and tools

Genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of complex traits and diseases. GWAS summary statistics have become essential tools for various genetic analyses, including meta-analysis, fine-mapping, and risk prediction. However, the increasing number of GWAS summary statistics and the diversity of software tools available for their analysis can make it challenging for researchers to select the most appropriate tools for their specific needs. This systematic review aims to provide a comprehensive overview of the currently available software tools and databases for GWAS summary statistics analysis. We conducted a comprehensive literature search to identify relevant software tools and databases. We categorized the tools and databases by their functionality, including data management, quality control, single-trait analysis, and multiple-trait analysis. We also compared the tools and databases based on their features, limitations, and user-friendliness. Our review identified a total of 305 functioning software tools and databases dedicated to GWAS summary statistics, each with unique strengths and limitations. We provide descriptions of the key features of each tool and database, including their input/output formats, data types, and computational requirements. We also discuss the overall usability and applicability of each tool for different research scenarios. This comprehensive review will serve as a valuable resource for researchers who are interested in using GWAS summary statistics to investigate the genetic basis of complex traits and diseases. By providing a detailed overview of the available tools and databases, we aim to facilitate informed tool selection and maximize the effectiveness of GWAS summary statistics analysis.

Genetic prediction of blood metabolites mediating the relationship between gut microbiota and Alzheimer's disease: a Mendelian randomization study

Background

Observational studies have suggested an association between gut microbiota and Alzheimer's disease (AD); however, the causal relationship remains unclear, and the role of blood metabolites in this association remains elusive.

Purpose

To elucidate the causal relationship between gut microbiota and AD and to investigate whether blood metabolites serve as potential mediators.

Materials and methods

Univariable Mendelian randomization (UVMR) analysis was employed to assess the causal relationship between gut microbiota and AD, while multivariable MR (MVMR) was utilized to mitigate confounding factors. Subsequently, a two-step mediation MR approach was employed to explore the role of blood metabolites as potential mediators. We primarily utilized the inverse variance-weighted method to evaluate the causal relationship between exposure and outcome, and sensitivity analyses including Contamination mixture, Maximum-likelihood, Debiased inverse-variance weighted, MR-Egger, Bayesian Weighted Mendelian randomization, and MR pleiotropy residual sum and outlier were conducted to address pleiotropy.

Results

After adjustment for reverse causality and MVMR correction, class Actinobacteria (OR: 1.03, 95% CI: 1.01-1.06, p = 0.006), family Lactobacillaceae (OR: 1.03, 95% CI: 1.00-1.05, p = 0.017), genus Lachnoclostridium (OR: 1.03, 95% CI: 1.00-1.06, p = 0.019), genus Ruminiclostridium9 (OR: 0.97, 95% CI: 0.94-1.00, p = 0.027) and genus Ruminiclostridium6 (OR: 1.03, 95% CI: 1.01-1.05, p = 0.009) exhibited causal effects on AD. Moreover, 1-ribosyl-imidazoleacetate levels (-6.62%), Metabolonic lactone sulfate levels (2.90%), and Nonadecanoate (19:0) levels (-12.17%) mediated the total genetic predictive effects of class Actinobacteria on AD risk. Similarly, 2-stearoyl-GPE (18:0) levels (-9.87%), Octadecanedioylcarnitine (C18-DC) levels (4.44%), 1-(1-enyl-stearoyl)-2-oleoyl-GPE (p-18:0/18:1) levels (38.66%), and X-23639 levels (13.28%) respectively mediated the total genetic predictive effects of family Lactobacillaceae on AD risk. Furthermore, Hexadecanedioate (C16-DC) levels (5.45%) mediated the total genetic predictive effects of genus Ruminiclostridium 6 on AD risk; Indole-3-carboxylate levels (13.91%), X-13431 levels (7.08%), Alpha-ketoglutarate to succinate ratio (-13.91%), 3-phosphoglycerate to glycerate ratio (15.27%), and Succinate to proline ratio (-14.64%) respectively mediated the total genetic predictive effects of genus Ruminiclostridium 9 on AD risk.

Conclusion

Our mediation MR analysis provides genetic evidence suggesting the potential mediating role of blood metabolites in the causal relationship between gut microbiota and AD. Further large-scale randomized controlled trials are warranted to validate the role of blood metabolites in the specific mechanisms by which gut microbiota influence AD.

Using Genetics to Investigate Relationships between Phenotypes: Application to Endometrial Cancer

Genome-wide association studies (GWAS) have accelerated the exploration of genotype-phenotype associations, facilitating the discovery of replicable genetic markers associated with specific traits or complex diseases. This narrative review explores the statistical methodologies developed using GWAS data to investigate relationships between various phenotypes, focusing on endometrial cancer, the most prevalent gynecological malignancy in developed nations. Advancements in analytical techniques such as genetic correlation, colocalization, cross-trait locus identification, and causal inference analyses have enabled deeper exploration of associations between different phenotypes, enhancing statistical power to uncover novel genetic risk regions. These analyses have unveiled shared genetic associations between endometrial cancer and many phenotypes, enabling identification of novel endometrial cancer risk loci and furthering our understanding of risk factors and biological processes underlying this disease. The current status of research in endometrial cancer is robust; however, this review demonstrates that further opportunities exist in statistical genetics that hold promise for advancing the understanding of endometrial cancer and other complex diseases.

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.

simmrd: An open-source tool to perform simulations in Mendelian randomization

Mendelian randomization (MR) has become a popular tool for inferring causality of risk factors on disease. There are currently over 45 different methods available to perform MR, reflecting this extremely active research area. It would be desirable to have a standard simulation environment to objectively evaluate the existing and future methods. We present simmrd, an open-source software for performing simulations to evaluate the performance of MR methods in a range of scenarios encountered in practice. Researchers can directly modify the simmrd source code so that the research community may arrive at a widely accepted framework for researchers to evaluate the performance of different MR methods.

Genetic evidence strengthens the bidirectional connection between gut microbiota and Shigella infection: insights from a two-sample Mendelian randomization study

Background

In recent investigations, substantial strides have been made in the precise modulation of the gut microbiota to prevent and treat a myriad of diseases. Simultaneously, the pressing issue of widespread antibiotic resistance and multidrug resistance resulting from Shigella infections demands urgent attention. Several studies suggest that the antagonistic influence of the gut microbiota could serve as a novel avenue for impeding the colonization of pathogenic microorganisms or treating Shigella infections. However, conventional research methodologies encounter inherent challenges in identifying antagonistic microbial agents against Shigella, necessitating a comprehensive and in-depth analysis of the causal relationship between Shigella infections and the gut microbiota.

Materials and methods

Utilizing the aggregated summary statistics from Genome-Wide Association Studies (GWAS), we conducted Mendelian Randomization (MR) analyses encompassing 18,340 participants to explore the interplay between the gut microbiota and Shigella infections. This investigation also involved 83 cases of Shigella infection patients and 336,396 control subjects. In the positive strand of our findings, we initially performed a preliminary analysis using the Inverse Variance Weighting (IVW) method. Subsequently, we undertook sensitivity analyses to assess the robustness of the results, addressing confounding factors' influence. This involved employing the Leave-One-Out method and scrutinizing funnel plots to ensure the reliability of the MR analysis outcomes. Conclusively, a reverse MR analysis was carried out, employing the Wald ratio method due to the exposure of individual Single Nucleotide Polymorphisms (SNPs). This was undertaken to explore the plausible associations between Shigella infections and genetically predicted compositions of the gut microbiota.

Results

In this study, we employed 2,818 SNPs associated with 211 species of gut microbiota as instrumental variables (IVs). Through IVW analysis, our positive MR findings revealed a significant negative correlation between the occurrence of Shigella infections and the phylum Tenericutes (OR: 0.18, 95% CI: 0.04-0.74, p = 0.02), class Mollicutes (OR: 0.18, 95% CI: 0.04-0.74, p = 0.02), genus Intestinimonas (OR: 0.16, 95% CI: 0.04-0.63, p = 0.01), genus Gordonibacter (OR: 0.39, 95% CI: 0.16-0.93, p = 0.03), and genus Butyrivibrio (OR: 0.44, 95% CI: 0.23-0.87, p = 0.02). Conversely, a positive correlation was observed between the occurrence of Shigella infections and genus Sutterella (OR: 10.16, 95% CI: 1.87-55.13, p = 0.01) and genus Alistipes (OR: 12.24, 95% CI: 1.71-87.34, p = 0.01). In sensitivity analyses, utilizing MR-Egger regression analysis and MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) detection, all outcomes demonstrated robust stability. Simultaneously, in the reverse MR analysis, Shigella infections resulted in an upregulation of four bacterial genera and a downregulation of three bacterial genera.

Conclusion

In summation, the MR analysis outcomes corroborate the presence of bidirectional causal relationships between the gut microbiota and Shigella infections. This study not only unveils novel perspectives for the prevention and treatment of Shigella infections but also furnishes fresh insights into the mechanistic underpinnings of how the gut microbiota contributes to the pathogenesis of Shigella infections. Consequently, the established dual causal association holds promise for advancing our understanding and addressing the complexities inherent in the interplay between the gut microbiota and Shigella infections, thereby paving the way for innovative therapeutic interventions and preventive strategies in the realm of Shigella-related diseases.

Impact of gut microbiota on nonalcoholic fatty liver disease: insights from a leave-one-out cross-validation study

Introduction

Enteric dysbacteriosis is strongly associated with nonalcoholic fatty liver disease (NAFLD). However, the underlying causal relationship remains unknown. Thus, the present study aimed to investigate the relationship between gut microbiota and NAFLD using Mendelian randomization (MR) and analyze the target genes potentially regulated by specific microbiota.

Methods

Bidirectional two-sample MR analysis was performed using inverse variance weighted (IVW) supplemented by MR-Egger, weighted median, simple mode, and weighted mode methods. Data were pooled from gut microbiota and NAFLD association studies. The least absolute shrinkage, selection operator regression, and the Support Vector Machine algorithm were used to identify genes regulated by these intestinal flora in NAFLD. The liver expression of these genes was verified in methionine choline-deficient (MCD) diet-fed mice.

Results

IVW results confirmed a causal relationship between eight specific gut microbes and NAFLD. Notably, the order Actinomycetales, NB1n, the family Actinomycetaceae, Oxalobacteraceae and the genus Ruminococcaceae UCG005 were positively correlated, whereas Lactobacillaceae, the Christensenellaceae R7 group, and Intestinibacter were negatively correlated with NAFLD onset. In NAFLD, these eight bacteria regulated four genes: colony-stimulating factor 2 receptor β, fucosyltransferase 2, 17-beta-hydroxysteroid dehydrogenase 14, and microtubule affinity regulatory kinase 3 (MAPK3). All genes, except MARK3, were differentially expressed in the liver tissues of MCD diet-fed mice.

Discussion

The abundance of eight gut microbiota species and NAFLD progression displayed a causal relationship based on the expression of the four target genes. Our findings contributed to the advancement of intestinal microecology-based diagnostic technologies and targeted therapies for NAFLD.

The causal effects between selenium levels and the brain cortical structure: A two-sample Mendelian randomization study

Extensive research has demonstrated the critical role of selenium (Se) and selenoproteins in brain function and cognition. However, the impact of Se on brain cortical structure remains enigmatic. Therefore, this study used Mendelian randomization (MR) analysis to investigate the causal effect between Se levels and brain cortical structure.

METHODS

This study utilizes 11 genetic variants associated with Se level variations, extracted from a large-scale genome-wide association study (GWAS) encompassed circulating Se levels (n = 5477) and toenail Se levels (n = 4162) in the European population. Outcome data were sourced from the summary statistics of the ENIGMA Consortium, comprising GWAS data from 51,666 individuals. The variables include cortical surface area (SA), thickness (TH) at the global level, and 34 functional cortical regions evaluated by magnetic resonance imaging. The inverse-variance-weighted method was used as the primary estimate. Additionally, sensitivity analyses were conducted to detect potential violations of assumptions underlying MR.

RESULTS

At the global level, Se levels were not correlated with SA but showed a significant negative correlation with TH (β = -0.00485 mm, SE = 0.00192, p = .0115). Heterogeneity was observed across different brain regions, with positive correlations found between Se levels and the TH of the parahippocampal gyrus, superior frontal gyrus, and frontal pole, whereas negative correlations were found with the TH of the inferior parietal lobe and middle temporal lobe. Regarding SA, Se levels exhibit positive correlations with the pars triangularis, caudal anterior cingulate, inferior parietal lobe, and banks of the superior temporal sulcus. Conversely, negative correlations were observed with the medial orbitofrontal cortex, posterior cingulate gyrus, insula, and the middle, superior, and transverse gyrus of the temporal lobe. No pleiotropy was detected.

RESULTS

This MR study indicated that Se levels causally influence the brain cortical structure.

Association between the gut microbiota and nonalcoholic fatty liver disease: A two-sample Mendelian randomization study

BACKGROUND

Increasing studies have shown that there is a significant association between gut microbiota and non-alcoholic fatty liver disease.

AIMS

To show the potential association between gut microbiota and non-alcoholic fatty liver disease, we performed a two-sample Mendelian randomization analysis.

METHODS

We analyzed summary statistics from genome-wide association studies of gut microbiota and non-alcoholic fatty liver disease and conducted Mendelian randomization studies to evaluate relationships between these factors.

RESULTS

Of the 211 gut microbiota taxa examined, the inverse variance weighted method identified Lactobacillaceae (OR = 0.83, 95% CI = 0.72 - 0.95, P = 0.007), Christensenellaceae (OR = 0.74, 95% CI = 0.59 - 0.92, P = 0.007), and Intestinibacter (OR = 0.85, 95% CI = 0.73 - 0.99, P = 0.035) were negatively correlated with non-alcoholic fatty liver disease. And Coriobacteriia (OR = 1.22, 95% CI = 1.01 - 1.42, P = 0.038), Actinomycetales (OR = 1.25, 95% CI = 1.02 - 1.53, P = 0.031), Oxalobacteraceae (OR = 1.10, 95% CI = 1.01 - 1.21, P = 0.036), Ruminococcaceae_UCG005 (OR = 1.18, 95% CI = 1.01 - 1.38, P = 0.033) are positively associated with non-alcoholic fatty liver disease.

CONCLUSIONS

Our study found that the abundance of certain strains was associated with the progression of nonalcoholic fatty liver disease.

Modeling of horizontal pleiotropy identifies possible causal gene expression in systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is a chronic autoimmune condition with complex causes involving genetic and environmental factors. While genome-wide association studies (GWASs) have identified genetic loci associated with SLE, the functional genomic elements responsible for disease development remain largely unknown. Mendelian Randomization (MR) is an instrumental variable approach to causal inference based on data from observational studies, where genetic variants are employed as instrumental variables (IVs).

Methods

This study utilized a two-step strategy to identify causal genes for SLE. In the first step, the classical MR method was employed, assuming the absence of horizontal pleiotropy, to estimate the causal effect of gene expression on SLE. In the second step, advanced probabilistic MR methods (PMR-Egger, MRAID, and MR-MtRobin) were applied to the genes identified in the first step, considering horizontal pleiotropy, to filter out false positives. PMR-Egger and MRAID analyses utilized whole blood expression quantitative trait loci (eQTL) and SLE GWAS summary data, while MR-MtRobin analysis used an independent eQTL dataset from multiple immune cell types along with the same SLE GWAS data.

Results

The initial MR analysis identified 142 genes, including 43 outside of chromosome 6. Subsequently, applying the advanced MR methods reduced the number of genes with significant causal effects on SLE to 66. PMR-Egger, MRAID, and MR-MtRobin, respectively, identified 13, 7, and 16 non-chromosome 6 genes with significant causal effects. All methods identified expression of PHRF1 gene as causal for SLE. A comprehensive literature review was conducted to enhance understanding of the functional roles and mechanisms of the identified genes in SLE development.

Conclusions

The findings from the three MR methods exhibited overlapping genes with causal effects on SLE, demonstrating consistent results. However, each method also uncovered unique genes due to different modelling assumptions and technical factors, highlighting the complementary nature of the approaches. Importantly, MRAID demonstrated a reduced percentage of causal genes from the Major Histocompatibility complex (MHC) region on chromosome 6, indicating its potential in minimizing false positive findings. This study contributes to unraveling the mechanisms underlying SLE by employing advanced probabilistic MR methods to identify causal genes, thereby enhancing our understanding of SLE pathogenesis.

Leveraging trans-ethnic genetic risk scores to improve association power for complex traits in underrepresented populations

Trans-ethnic genome-wide association studies have revealed that many loci identified in European populations can be reproducible in non-European populations, indicating widespread trans-ethnic genetic similarity. However, how to leverage such shared information more efficiently in association analysis is less investigated for traits in underrepresented populations. We here propose a statistical framework, trans-ethnic genetic risk score informed gene-based association mixed model (GAMM), by hierarchically modeling single-nucleotide polymorphism effects in the target population as a function of effects of the same trait in well-studied populations. GAMM powerfully integrates genetic similarity across distinct ancestral groups to enhance power in understudied populations, as confirmed by extensive simulations. We illustrate the usefulness of GAMM via the application to 13 blood cell traits (i.e. basophil count, eosinophil count, hematocrit, hemoglobin concentration, lymphocyte count, mean corpuscular hemoglobin, mean corpuscular hemoglobin concentration, mean corpuscular volume, monocyte count, neutrophil count, platelet count, red blood cell count and total white blood cell count) in Africans of the UK Biobank (n = 3204) while utilizing genetic overlap shared in Europeans (n = 746 667) and East Asians (n = 162 255). We discovered multiple new associated genes, which had otherwise been missed by existing methods, and revealed that the trans-ethnic information indirectly contributed much to the phenotypic variance. Overall, GAMM represents a flexible and powerful statistical framework of association analysis for complex traits in underrepresented populations by integrating trans-ethnic genetic similarity across well-studied populations, and helps attenuate health inequities in current genetics research for people of minority populations.

Deconstructing a Syndrome: Genomic Insights Into PCOS Causal Mechanisms and Classification

Polycystic ovary syndrome (PCOS) is among the most common disorders in women of reproductive age, affecting up to 15% worldwide, depending on the diagnostic criteria. PCOS is characterized by a constellation of interrelated reproductive abnormalities, including disordered gonadotropin secretion, increased androgen production, chronic anovulation, and polycystic ovarian morphology. It is frequently associated with insulin resistance and obesity. These reproductive and metabolic derangements cause major morbidities across the lifespan, including anovulatory infertility and type 2 diabetes (T2D). Despite decades of investigative effort, the etiology of PCOS remains unknown. Familial clustering of PCOS cases has indicated a genetic contribution to PCOS. There are rare Mendelian forms of PCOS associated with extreme phenotypes, but PCOS typically follows a non-Mendelian pattern of inheritance consistent with a complex genetic architecture, analogous to T2D and obesity, that reflects the interaction of susceptibility genes and environmental factors. Genomic studies of PCOS have provided important insights into disease pathways and have indicated that current diagnostic criteria do not capture underlying differences in biology associated with different forms of PCOS. We provide a state-of-the-science review of genetic analyses of PCOS, including an overview of genomic methodologies aimed at a general audience of non-geneticists and clinicians. Applications in PCOS will be discussed, including strengths and limitations of each study. The contributions of environmental factors, including developmental origins, will be reviewed. Insights into the pathogenesis and genetic architecture of PCOS will be summarized. Future directions for PCOS genetic studies will be outlined.

Proteome-wide Mendelian randomization in global biobank meta-analysis reveals multi-ancestry drug targets for common diseases

Proteome-wide Mendelian randomization (MR) shows value in prioritizing drug targets in Europeans but with limited evidence in other ancestries. Here, we present a multi-ancestry proteome-wide MR analysis based on cross-population data from the Global Biobank Meta-analysis Initiative (GBMI). We estimated the putative causal effects of 1,545 proteins on eight diseases in African (32,658) and European (1,219,993) ancestries and identified 45 and 7 protein-disease pairs with MR and genetic colocalization evidence in the two ancestries, respectively. A multi-ancestry MR comparison identified two protein-disease pairs with MR evidence in both ancestries and seven pairs with specific effects in the two ancestries separately. Integrating these MR signals with clinical trial evidence, we prioritized 16 pairs for investigation in future drug trials. Our results highlight the value of proteome-wide MR in informing the generalizability of drug targets for disease prevention across ancestries and illustrate the value of meta-analysis of biobanks in drug development.

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.

Association Between Telomere Length and Skin Cancer and Aging: A Mendelian Randomization Analysis

Background: Telomere shortening is a hallmark of cellular senescence. However, telomere length (TL)-related cellular senescence has varying effects in different cancers, resulting in a paradoxical relationship between senescence and cancer. Therefore, we used observational epidemiological studies to investigate the association between TL and skin cancer and aging, and to explore whether such a paradoxical relationship exists in skin tissue.

Methods: This study employed two-sample Mendelian randomization (MR) to analyze the causal relationship between TL and skin cancer [melanoma and non-melanoma skin cancers (NMSCs)] and aging. We studied single nucleotide polymorphisms (SNPs) obtained from pooled data belonging to genome-wide association studies (GWAS) in the literature and biobanks. Quality control was performed using pleiotropy, heterogeneity, and sensitivity analyses.

Results: We used five algorithms to analyze the causal relationship between TL and skin aging, melanoma, and NMSCs, and obtained consistent results. TL shortening reduced NMSC and melanoma susceptibility risk with specific odds ratios (ORs) of 1.0344 [95% confidence interval (CI): 1.0168–1.0524, p = 0.01] and 1.0127 (95% CI: 1.0046–1.0209, p = 6.36E-07), respectively. Conversely, TL shortening was validated to increase the odds of skin aging (OR = 0.96, 95% CI: 0.9332–0.9956, p = 0.03). Moreover, the MR-Egger, maximum likelihood, and inverse variance weighted (IVW) methods found significant heterogeneity among instrumental variable (IV) estimates (identified as MR-Egger skin aging Q = 76.72, p = 1.36E-04; melanoma Q = 97.10, p = 1.62E-07; NMSCsQ = 82.02, p = 1.90E-05). The leave-one-out analysis also showed that the SNP sensitivity was robust to each result.

Conclusion: This study found that TL shortening may promote skin aging development and reduce the risk of cutaneous melanoma and NMSCs. The results provide a reference for future research on the causal relationship between skin aging and cancer in clinical practice.

Statistical methods for Mendelian randomization in genome-wide association studies: A review

Genome-wide association studies have yielded thousands of associations for many common diseases and disease-related complex traits. The identified associations made it possible to identify the causal risk factors underlying diseases and investigate the causal relationships among complex traits through Mendelian randomization. Mendelian randomization is a form of instrumental variable analysis that uses SNP associations from genome-wide association studies as instruments to study and uncover causal relationships between complex traits. By leveraging SNP genotypes as instrumental variables, or proxies, for the exposure complex trait, investigators can tease out causal effects from observational data, provided that necessary assumptions are satisfied. We discuss below the development of Mendelian randomization methods in parallel with the growth of genome-wide association studies. We argue that the recent availability of GWAS summary statistics for diverse complex traits has motivated new Mendelian randomization methods with relaxed causality assumptions and that this area continues to offer opportunities for robust biological discoveries.

A two-sample robust Bayesian Mendelian Randomization method accounting for linkage disequilibrium and idiosyncratic pleiotropy with applications to the COVID-19 outcomes

Mendelian randomization (MR) is a statistical method exploiting genetic variants as instrumental variables to estimate the causal effect of modifiable risk factors on an outcome of interest. Despite wide uses of various popular two-sample MR methods based on genome-wide association study summary level data, however, those methods could suffer from potential power loss or/and biased inference when the chosen genetic variants are in linkage disequilibrium (LD), and also have relatively large direct effects on the outcome whose distribution might be heavy-tailed which is commonly referred to as the idiosyncratic pleiotropy phenomenon. To resolve those two issues, we propose a novel Robust Bayesian Mendelian Randomization (RBMR) model that uses the more robust multivariate generalized t $t$ -distribution to model such direct effects in a probabilistic model framework which can also incorporate the LD structure explicitly. The generalized t $t$ -distribution can be represented as a Gaussian scaled mixture so that our model parameters can be estimated by the expectation maximization (EM)-type algorithms. We compute the standard errors by calibrating the evidence lower bound using the likelihood ratio test. Through extensive simulation studies, we show that our RBMR has robust performance compared with other competing methods. We further apply our RBMR method to two benchmark data sets and find that RBMR has smaller bias and standard errors. Using our proposed RBMR method, we find that coronary artery disease is associated with increased risk of critically ill coronavirus disease 2019. We also develop a user-friendly R package RBMR (https://github.com/AnqiWang2021/RBMR) for public use.

Mendelian randomization

Mendelian randomization (MR) is a term that applies to the use of genetic variation to address causal questions about how modifiable exposures influence different outcomes. The principles of MR are based on Mendel's laws of inheritance and instrumental variable estimation methods, which enable the inference of causal effects in the presence of unobserved confounding. In this Primer, we outline the principles of MR, the instrumental variable conditions underlying MR estimation and some of the methods used for estimation. We go on to discuss how the assumptions underlying an MR study can be assessed and give methods of estimation that are robust to certain violations of these assumptions. We give examples of a range of studies in which MR has been applied, the limitations of current methods of analysis and the outlook for MR in the future. The difference between the assumptions required for MR analysis and other forms of non-interventional epidemiological studies means that MR can be used as part of a triangulation across multiple sources of evidence for causal inference.

Addressing dyslipidemic risk beyond LDL-cholesterol

Despite the success of LDL-lowering drugs in reducing cardiovascular disease (CVD), there remains a large burden of residual disease due in part to persistent dyslipidemia characterized by elevated levels of triglyceride-rich lipoproteins (TRLs) and reduced levels of HDL. This form of dyslipidemia is increasing globally as a result of the rising prevalence of obesity and metabolic syndrome. Accumulating evidence suggests that impaired hepatic clearance of cholesterol-rich TRL remnants leads to their accumulation in arteries, promoting foam cell formation and inflammation. Low levels of HDL may associate with reduced cholesterol efflux from foam cells, aggravating atherosclerosis. While fibrates and fish oils reduce TRL, they have not been uniformly successful in reducing CVD, and there is a large unmet need for new approaches to reduce remnants and CVD. Rare genetic variants that lower triglyceride levels via activation of lipolysis and associate with reduced CVD suggest new approaches to treating dyslipidemia. Apolipoprotein C3 (APOC3) and angiopoietin-like 3 (ANGPTL3) have emerged as targets for inhibition by antibody, antisense, or RNAi approaches. Inhibition of either molecule lowers TRL but respectively raises or lowers HDL levels. Large clinical trials of such agents in patients with high CVD risk and elevated levels of TRL will be required to demonstrate efficacy of these approaches.

MR-Corr2: a two-sample Mendelian randomization method that accounts for correlated horizontal pleiotropy using correlated instrumental variants

MOTIVATION

Mendelian randomization (MR) is a valuable tool to examine the causal relationships between health risk factors and outcomes from observational studies. Along with the proliferation of genome-wide association studies, a variety of two-sample MR methods for summary data have been developed to account for horizontal pleiotropy (HP), primarily based on the assumption that the effects of variants on exposure (γ) and HP (α) are independent. In practice, this assumption is too strict and can be easily violated because of the correlated HP.

RESULTS

To account for this correlated HP, we propose a Bayesian approach, MR-Corr2, that uses the orthogonal projection to reparameterize the bivariate normal distribution for γ and α, and a spike-slab prior to mitigate the impact of correlated HP. We have also developed an efficient algorithm with paralleled Gibbs sampling. To demonstrate the advantages of MR-Corr2 over existing methods, we conducted comprehensive simulation studies to compare for both type-I error control and point estimates in various scenarios. By applying MR-Corr2 to study the relationships between exposure-outcome pairs in complex traits, we did not identify the contradictory causal relationship between HDL-c and CAD. Moreover, the results provide a new perspective of the causal network among complex traits.

AVAILABILITY AND IMPLEMENTATION

The developed R package and code to reproduce all the results are available at https://github.com/QingCheng0218/MR.Corr2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Performing post-genome-wide association study analysis: overview, challenges and recommendations

Genome-wide association studies (GWAS) provide  huge information on statistically significant single-nucleotide polymorphisms (SNPs) associated with various human complex traits and diseases. By performing GWAS studies, scientists have successfully identified the association of hundreds of thousands to  millions of SNPs to a single phenotype. Moreover, the association of some SNPs with rare diseases has been intensively tested. However, classic GWAS studies have not yet provided solid, knowledgeable insight into functional and biological mechanisms underlying phenotypes or mechanisms of diseases. Therefore, several post-GWAS (pGWAS) methods have been recommended. Currently, there is no simple scientific document to provide a quick guide for performing pGWAS analysis. pGWAS is a crucial step for a better understanding of the biological machinery beyond the SNPs. Here, we provide an overview to performing pGWAS analysis and demonstrate the challenges behind each method. Furthermore, we direct readers to key articles for each pGWAS method and present the overall issues in pGWAS analysis.  Finally, we include a custom pGWAS pipeline to guide new users when performing their research.

A Novel Method for Mendelian Randomization Analyses With Pleiotropy and Linkage Disequilibrium in Genetic Variants From Individual Data

Mendelian randomization makes use of genetic variants as instrumental variables to eliminate the influence induced by unknown confounders on causal estimation in epidemiology studies. However, with the soaring genetic variants identified in genome-wide association studies, the pleiotropy, and linkage disequilibrium in genetic variants are unavoidable and may produce severe bias in causal inference. In this study, by modeling the pleiotropic effect as a normally distributed random effect, we propose a novel mixed-effects regression model-based method PLDMR, pleiotropy and linkage disequilibrium adaptive Mendelian randomization, which takes linkage disequilibrium into account and also corrects for the pleiotropic effect in causal effect estimation and statistical inference. We conduct voluminous simulation studies to evaluate the performance of the proposed and existing methods. Simulation results illustrate the validity and advantage of the novel method, especially in the case of linkage disequilibrium and directional pleiotropic effects, compared with other methods. In addition, by applying this novel method to the data on Atherosclerosis Risk in Communications Study, we conclude that body mass index has a significant causal effect on and thus might be a potential risk factor of systolic blood pressure. The novel method is implemented in R and the corresponding R code is provided for free download.

A robust two-sample transcriptome-wide Mendelian randomization method integrating GWAS with multi-tissue eQTL summary statistics

By treating genetic variants as instrumental variables (IVs), two-sample Mendelian randomization (MR) methods detect genetically regulated risk exposures for complex diseases using only summary statistics. When considering gene expression as exposure in transcriptome-wide MR (TWMR) analyses, the eQTLs (expression-quantitative-trait-loci) may have pleiotropic effects or be correlated with variants that have effects on disease not via expression, and the presence of those invalid IVs would lead to biased inference. Moreover, the number of eQTLs as IVs for a gene is generally limited, making the detection of invalid IVs challenging. We propose a method, "MR-MtRobin," for accurate TWMR inference in the presence of invalid IVs. By leveraging multi-tissue eQTL data in a mixed model, the proposed method makes identifiable the IV-specific random effects due to pleiotropy from estimation errors of eQTL summary statistics, and can provide accurate inference on the dependence (fixed effects) between eQTL and GWAS (genome-wide association study) effects in the presence of invalid IVs. Moreover, our method can improve power and precision in inference by selecting cross-tissue eQTLs as IVs that have improved consistency of effects across eQTL and GWAS data. We applied MR-MtRobin to detect genes associated with schizophrenia risk by integrating summary-level data from the Psychiatric Genomics Consortium and the Genotype-Tissue Expression project (V8).

Aggregating multiple expression prediction models improves the power of transcriptome-wide association studies

Transcriptome-wide association study (TWAS) is an important integrative method for identifying genes that are causally associated with phenotypes. A key step of TWAS involves the construction of expression prediction models for every gene in turn using its cis-SNPs as predictors. Different TWAS methods rely on different models for gene expression prediction, and each such model makes a distinct modeling assumption that is often suitable for a particular genetic architecture underlying expression. However, the genetic architectures underlying gene expression vary across genes throughout the transcriptome. Consequently, different TWAS methods may be beneficial in detecting genes with distinct genetic architectures. Here, we develop a new method, HMAT, which aggregates TWAS association evidence obtained across multiple gene expression prediction models by leveraging the harmonic mean P-value combination strategy. Because each expression prediction model is suited to capture a particular genetic architecture, aggregating TWAS associations across prediction models as in HMAT improves accurate expression prediction and enables subsequent powerful TWAS analysis across the transcriptome. A key feature of HMAT is its ability to accommodate the correlations among different TWAS test statistics and produce calibrated P-values after aggregation. Through numerical simulations, we illustrated the advantage of HMAT over commonly used TWAS methods as well as ad hoc P-value combination rules such as Fisher's method. We also applied HMAT to analyze summary statistics of nine common diseases. In the real data applications, HMAT was on average 30.6% more powerful compared to the next best method, detecting many new disease-associated genes that were otherwise not identified by existing TWAS approaches. In conclusion, HMAT represents a flexible and powerful TWAS method that enjoys robust performance across a range of genetic architectures underlying gene expression.