Analysis of multiple related phenotypes in genome-wide association studies

Transethnic meta-analysis of exome-wide variants identifies new loci associated with male-specific metabolic syndrome

BACKGROUND

Metabolic syndrome (MetS) is a group of very common human conditions promoting strong understand the impact of rare variants, beyond exome-wide association studies, to potentially discover causative variants, across different ethnic populations.

OBJECTIVE

We performed transethnic, exome-wide MetS association studies on MetS in men.

METHODS

We analyzed genotype data of 5302 European subjects (2658 cases and 2644 controls), in the discovery stage of the European METabolic Syndrome In Men study, generated from exome chips, and 2481 subjects (714 cases and 1767 controls), in the replication stage, across 6 independent cohorts of 5 ancestries (T2D-GENES consortium), using whole-exome sequencing. We therefore evaluated gene-level and variant-level associations, of rare variants for MetS, using logistic regression (LR) and multivariate analyses (MulA).

RESULTS

Gene-based association found the gene for the cholesteryl ester transfer protein (CETP) (from MulA, p value = 4.67 × 10-9; from LR, p value = 0.009) to well associate with MetS. At two missense variants, from 8 rare variants in CETP, Ala390Pro (rs5880) (from MulA, p value = 1.28 × 10-7; from LR, p value = 1.34 × 10-4) and Arg468Gln (rs1800777) (from MulA, p value = 2.40 × 10-5; from LR, p value = 1.49 × 10-3) significantly associated with MetS across five ancestries.

CONCLUSIONS

Our findings highlight novel rare variants of genes that confer MetS susceptibility, in Europeans, that are shared with diverse populations, emphasizing an opportunity to further understand the biological target or genes that underlie MetS, across populations.

Spatial rank-based multifactor dimensionality reduction to detect gene-gene interactions for multivariate phenotypes

Background

Identifying interaction effects between genes is one of the main tasks of genome-wide association studies aiming to shed light on the biological mechanisms underlying complex diseases. Multifactor dimensionality reduction (MDR) is a popular approach for detecting gene–gene interactions that has been extended in various forms to handle binary and continuous phenotypes. However, only few multivariate MDR methods are available for multiple related phenotypes. Current approaches use Hotelling’s T² statistic to evaluate interaction models, but it is well known that Hotelling’s T² statistic is highly sensitive to heavily skewed distributions and outliers.

Results

We propose a robust approach based on nonparametric statistics such as spatial signs and ranks. The new multivariate rank-based MDR (MR-MDR) is mainly suitable for analyzing multiple continuous phenotypes and is less sensitive to skewed distributions and outliers. MR-MDR utilizes fuzzy k-means clustering and classifies multi-locus genotypes into two groups. Then, MR-MDR calculates a spatial rank-sum statistic as an evaluation measure and selects the best interaction model with the largest statistic. Our novel idea lies in adopting nonparametric statistics as an evaluation measure for robust inference. We adopt tenfold cross-validation to avoid overfitting. Intensive simulation studies were conducted to compare the performance of MR-MDR with current methods. Application of MR-MDR to a real dataset from a Korean genome-wide association study demonstrated that it successfully identified genetic interactions associated with four phenotypes related to kidney function. The R code for conducting MR-MDR is available at https://github.com/statpark/MR-MDR.

Conclusions

Intensive simulation studies comparing MR-MDR with several current methods showed that the performance of MR-MDR was outstanding for skewed distributions. Additionally, for symmetric distributions, MR-MDR showed comparable power. Therefore, we conclude that MR-MDR is a useful multivariate non-parametric approach that can be used regardless of the phenotype distribution, the correlations between phenotypes, and sample size.

Enhanced Permutation Tests via Multiple Pruning

Big multi-omics data in bioinformatics often consists of a huge number of features and relatively small numbers of samples. In addition, features from multi-omics data have their own specific characteristics depending on whether they are from genomics, proteomics, metabolomics, etc. Due to these distinct characteristics, standard statistical analyses using parametric-based assumptions may sometimes fail to provide exact asymptotic results. To resolve this issue, permutation tests can be a way to exactly analyze multi-omics data because they are distribution-free and flexible to use. In permutation tests, p-values are evaluated by estimating the locations of test statistics in an empirical null distribution generated by random shuffling. However, the permutation approach can be infeasible when the number of features increases, because more stringent control of type I error is needed for multiple hypothesis testing, and consequently, much larger numbers of permutations are required to reach significance. To address this problem, we propose a well-organized strategy, “ENhanced Permutation tests via multiple Pruning (ENPP).” ENPP prunes the features in every permutation round if they are determined to be non-significant. In other words, if the feature statistics from the permuted datasets exceed the feature statistics from the original dataset, beyond a predetermined threshold, the feature is determined to be non-significant. If so, ENPP removes the feature and iterates the process without the feature in the next permutation round. Our simulation study showed that the ENPP method could remove about 50% of the features at the first permutation round, and, by the 100th permutation round, 98% of the features had been removed and only 7.4% of the computation time with the original unpruned permutation approach had elapsed. In addition, we applied this approach to a real data set (Korea Association REsource: KARE) of 327,872 SNPs to find association with a non-normally distributed phenotype (fasting plasma glucose), interpreted the results, and discussed the feasibility and advantages of the approach.

Multivariate Cluster-Based Multifactor Dimensionality Reduction to Identify Genetic Interactions for Multiple Quantitative Phenotypes

To understand the pathophysiology of complex diseases, including hypertension, diabetes, and autism, deleterious phenotypes are unlikely due to the effects of single genes, but rather, gene-gene interactions (GGIs), which are widely analyzed by multifactor dimensionality reduction (MDR). Early MDR methods mainly focused on binary traits. More recently, several extensions of MDR have been developed for analyzing various traits such as quantitative traits and survival times. Newer technologies, such as genome-wide association studies (GWAS), have now been developed for assessing multiple traits, to simultaneously identify genetic variants associated with various pathological phenotypes. It has also been well demonstrated that analyzing multiple traits has several advantages over single trait analysis. While there remains a need to find GGIs for multiple traits, such studies have become more difficult, due to a lack of novel methods and software. Herein, we propose a novel multi-CMDR method, by combining fuzzy clustering and MDR, to find GGIs for multiple traits. Multi-CMDR showed similar power to existing methods, when phenotypes followed bivariate normal distributions, and showed better power than others for skewed distributions. The validity of multi-CMDR was confirmed by analyzing real-life Korean GWAS data.

New Common and Rare Variants Influencing Metabolic Syndrome and Its Individual Components in a Korean Population

To identify novel loci for susceptibility to MetS, we conducted genome-wide association and exome wide association studies consisting of a discovery stage cohort (KARE, 1946 cases and 6427 controls), and a replication stage cohort (HEXA, 430 cases and 3,264 controls). For finding genetic variants for MetS, with its components, we performed multivariate analysis for common and rare associations, using a standard logistic regression analysis for MetS. From the discovery and replication GWA studies, we confirmed 21 genome-wide signals significantly associated with MetS. Of these 21, four were previously unreported to associate with any MetS components: rs765547 near LPL; rs3782889 in MYL2; and rs11065756 and rs10849915 in CCDC63. Using exome chip variants, gene-based analysis of rare variants revealed three genes, CETP, SH2B1, and ZFP2, in the discovery stage, among which only CETP was confirmed in the replication stage. Finally, CETP D442G (rs2303790) associated, as a less common variant, with decreased risk of MetS. In conclusion, we discovered a total of five new MetS-associated loci, and their overlap with other disease-related components, suggest roles in the various etiologies of MetS, and its possible preventive strategies.