Dataset Information

How Well Can Multivariate and Univariate GWAS Distinguish Between True and Spurious Pleiotropy?

ABSTRACT: Quantification of the simultaneous contributions of loci to multiple traits, a phenomenon called pleiotropy, is facilitated by the increased availability of high-throughput genotypic and phenotypic data. To understand the prevalence and nature of pleiotropy, the ability of multivariate and univariate genome-wide association study (GWAS) models to distinguish between pleiotropic and non-pleiotropic loci in linkage disequilibrium (LD) first needs to be evaluated. Therefore, we used publicly available maize and soybean genotypic data to simulate multiple pairs of traits that were either (i) controlled by quantitative trait nucleotides (QTNs) on separate chromosomes, (ii) controlled by QTNs in various degrees of LD with each other, or (iii) controlled by a single pleiotropic QTN. We showed that multivariate GWAS could not distinguish between QTNs in LD and a single pleiotropic QTN. In contrast, a unique QTN detection rate pattern was observed for univariate GWAS whenever the simulated QTNs were in high LD or pleiotropic. Collectively, these results suggest that multivariate and univariate GWAS should both be used to infer whether or not causal mutations underlying peak GWAS associations are pleiotropic. Therefore, we recommend that future studies use a combination of multivariate and univariate GWAS models, as both models could be useful for identifying and narrowing down candidate loci with potential pleiotropic effects for downstream biological experiments.

SUBMITTER: Fernandes SB

PROVIDER: S-EPMC7873880 | biostudies-literature | 2020

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

How Well Can Multivariate and Univariate GWAS Distinguish Between True and Spurious Pleiotropy?

Fernandes Samuel B SB Zhang Kevin S KS Jamann Tiffany M TM Lipka Alexander E AE

Frontiers in genetics 20210108

Quantification of the simultaneous contributions of loci to multiple traits, a phenomenon called pleiotropy, is facilitated by the increased availability of high-throughput genotypic and phenotypic data. To understand the prevalence and nature of pleiotropy, the ability of multivariate and univariate genome-wide association study (GWAS) models to distinguish between pleiotropic and non-pleiotropic loci in linkage disequilibrium (LD) first needs to be evaluated. Therefore, we used publicly availa ...[more]

PMID: 33584799

Similar Datasets

Project description:BackgroundObservational studies have shown associations between psychiatric disorders and osteoarthritis (OA). However, the causal impact of different psychiatric disorder types on specific sites of osteoarthritis remains unclear. This study aimed to comprehensively understand the potential causal associations between psychiatric disorders and osteoarthritis using Mendelian randomization (MR) analysis.MethodsWe collected data from genome-wide association studies of knee osteoarthritis (KOA) (n = 403,124), hip osteoarthritis (HOA) (n = 393,873), osteoarthritis of the knee or hip (KHOA) (n = 417,596), as well as three psychiatric disorders: bipolar disorder (n = 41,917), major depressive disorder (n = 170,756), and schizophrenia (n = 76,755) among European populations. We applied bidirectional univariate and multivariate MR analyses, including inverse variance weighted, Mendelian randomization-Egger, weighted median, simple mode, and weighted mode. We considered p < .05 as a criterion for identifying potential evidence of association. Bonferroni correction was used for multiple tests.ResultsOur univariate MR analysis results demonstrated that bipolar disorder is a protective factor for KOA (OR = 0.90, 95% CI = 0.83 to 0.97, p = 0.0048) and may also be protective for KHOA (p = 0.02). Conversely, major depression has a positive causal effect on both KOA (OR = 1.27; 95% CI = 1.08 to 1.49; p = 0.0036) and KHOA (OR = 1.24; 95% CI = 1.12 to 1.37; p = 3.62×10-05 ). Furthermore, our analysis suggested that KHOA may be a risk factor for major depression (OR = 1.06; 95% CI = 1.00 to 1.12; p = 0.0469) in reverse MR. After adjusting smoking (OR = 1.46; 95% CI = 1.19 to 1.65; p = 0.0032) and body mass index (OR = 1.44; 95% CI = 1.09 to 1.81; p = 8.56×10-04 ), the casual association between major depression and KHOA remained.ConclusionOur study indicates that major depression is a great risk factor for KHOA, increasing the likelihood of their occurrence. However, further in-depth studies will be required to validate these results and elucidate the underlying molecular mechanisms.

Project description:Orienting the causal relationship between pairs of traits is a fundamental task in scientific research with significant implications in practice, such as in prioritizing molecular targets and modifiable risk factors for developing therapeutic and interventional strategies for complex diseases. A recent method, called Steiger's method, using a single SNP as an instrument variable (IV) in the framework of Mendelian randomization (MR), has since been widely applied. We report the following new contributions. First, we propose a single SNP-based alternative, overcoming a severe limitation of Steiger's method in simply assuming, instead of inferring, the existence of a causal relationship. We also clarify a condition necessary for the validity of the methods in the presence of hidden confounding. Second, to improve statistical power, we propose combining the results from multiple, and possibly correlated, SNPs as multiple instruments. Third, we develop three goodness-of-fit tests to check modeling assumptions, including those required for valid IVs. Fourth, by relaxing one of the three IV assumptions in MR, we propose several methods, including an Egger regression-like approach and its multivariable version (analogous to multivariable MR), to account for horizontal pleiotropy of the SNPs/IVs, which is often unavoidable in practice. All our methods can simultaneously infer both the existence and (if so) the direction of a causal relationship, largely expanding their applicability over that of Steiger's method. Although we focus on uni-directional causal relationships, we also briefly discuss an extension to bi-directional relationships. Through extensive simulations and an application to infer the causal directions between low density lipoprotein (LDL) cholesterol, or high density lipoprotein (HDL) cholesterol, and coronary artery disease (CAD), we demonstrate the superior performance and advantage of our proposed methods over Steiger's method and bi-directional MR. In particular, after accounting for horizontal pleiotropy, our method confirmed the well known causal direction from LDL to CAD, while other methods, including bi-directional MR, might fail.

Project description:ObjectiveTo explore the impact of inflammatory factors on the incidence of cerebral small vessel disease (CSVD), we performed a mendelian randomization (MR) study to analyze the causal relationship between multiple inflammatory factors and CSVD imaging markers and utilized summary-data-based mendelian randomization (SMR) analysis to infer whether the impact of instrumental variables (IVs) on disease is mediated by gene expression or DNA methylation.MethodsUsing public databases such as UKB and IEU, and original genome-wide association studies, we obtained IVs related to exposure (inflammatory factors) and outcome (CSVD imaging markers). We performed the inverse variance weighted, weighted median, and MR-Egger methods to assess causal effects between exposure and outcome in univariate MR analysis. To evaluate their heterogeneity, a series of sensitivity analyses were conducted, including the Cochrane Q test, MR-Egger intercept test, MR-Presso, and leave-one-out analysis. We also applied mediation and multivariate MR analysis to explore the interactions between positive exposures on the same outcome. Additionally, we conducted the SMR, which utilizes instruments within or near relevant genes in blood or brain tissues, to elucidate the causal associations with CSVD markers.ResultsABO Univariate MR of multiple cohorts revealed that the risk of small vessel stroke (SVS) increases with elevated levels of TNF-related apoptosis-inducing ligand (TRAIL, OR, 1.23, 95% CI, 1.08-1.39) and interleukin-1 receptor-like 2, (IL-1RL2, OR, 1.29, 95% CI, 1.04-1.61). IL-18 was a potential risk factor for extensive basal ganglia perivascular space burden (BGPVS, OR, 1.02, 95% CI, 1.00-1.05). Moreover, the risk of extensive white matter perivascular space burden (WMPVS) decreased with rising levels of E-selectin (OR, .98, 95% CI, .97-1.00), IL-1RL2 (OR, .97, 95% CI, .95-1.00), IL-3 receptor subunit alpha (IL-3Ra, OR, .98, 95% CI, .97-1.00), and IL-5 receptor subunit alpha (IL-5Ra, OR, .98, 95% CI, .97-1.00). Mediation and multivariate MR analysis indicated that E-selectin and IL-3Ra might interact during the pathogenesis of WMPVS. SMR estimates showed that TRAIL-related IVs rs5030044 and rs2304456 increased the risk of SVS by increasing the expression of gene Kininogen-1 (KNG1) in the cerebral cortex, particularly in the frontal cortex (βsmr = .10, Psmr = .003, FDR = .04). Instruments (rs507666 and rs2519093) related to E-selectin and IL-3Ra could increase the risk of WMPVS by enhancing DNA methylation of the gene ABO in blood tissue (βsmr = .01-.02, Psmr = .001, FDR = .01-.03).ConclusionAccording to MR and SMR analysis, higher levels of TRAIL increased the risk of SVS by upregulating gene expression of KNG1 in brain cortex tissues. In addition, protective effects of E-selectin and IL-3a levels on WMPVS were regulated by increased DNA methylation of gene ABO in blood tissue.

Dataset Information

How Well Can Multivariate and Univariate GWAS Distinguish Between True and Spurious Pleiotropy?

Publications

How Well Can Multivariate and Univariate GWAS Distinguish Between True and Spurious Pleiotropy?

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets