Testing and estimation in marker-set association study using semiparametric quantile regression kernel machine

Statistical inference for high-dimensional pathway analysis with multiple responses

Pathway analysis, i.e., grouping analysis, has important applications in genomic studies. Existing pathway analysis approaches are mostly focused on a single response and are not suitable for analyzing complex diseases that are often related with multiple response variables. Although a handful of approaches have been developed for multiple responses, these methods are mainly designed for pathways with a moderate number of features. A multi-response pathway analysis approach that is able to conduct statistical inference when the dimension is potentially higher than sample size is introduced. Asymptotical properties of the test statistic are established and theoretical investigation of the statistical power is conducted. Simulation studies and real data analysis show that the proposed approach performs well in identifying important pathways that influence multiple expression quantitative trait loci (eQTL).

Testing microbiome association using integrated quantile regression models

MOTIVATION

Most existing microbiome association analyses focus on the association between microbiome and conditional mean of health or disease-related outcomes, and within this vein, vast computational tools and methods have been devised for standard binary or continuous outcomes. However, these methods tend to be limited either when the underlying microbiome-outcome association occurs somewhere other than the mean level, or when distribution of the outcome variable is irregular (e.g. zero-inflated or mixtures) such that conditional outcome mean is less meaningful. We address this gap by investigating association analysis between microbiome compositions and conditional outcome quantiles.

RESULTS

We introduce a new association analysis tool named MiRKAT-IQ within the Microbiome Regression-based Kernel Association Test framework using Integrated Quantile regression models to examine the association between microbiome and the distribution of outcome. For an individual quantile, we utilize the existing kernel machine regression framework to examine the association between that conditional outcome quantile and a group of microbial features (e.g. microbiome community compositions). Then, the goal of examining microbiome association with the whole outcome distribution is achieved by integrating all outcome conditional quantiles over a process, and thus our new MiRKAT-IQ test is robust to both the location of association signals (e.g. mean, variance, median) and the heterogeneous distribution of the outcome. Extensive numerical simulation studies have been conducted to show the validity of the new MiRKAT-IQ test. We demonstrate the potential usefulness of MiRKAT-IQ with applications to actual biological data collected from a previous microbiome study.

AVAILABILITY AND IMPLEMENTATION

R codes to implement the proposed methodology is provided in the MiRKAT package, which is available on CRAN.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

kTWAS: integrating kernel machine with transcriptome-wide association studies improves statistical power and reveals novel genes

The power of genotype-phenotype association mapping studies increases greatly when contributions from multiple variants in a focal region are meaningfully aggregated. Currently, there are two popular categories of variant aggregation methods. Transcriptome-wide association studies (TWAS) represent a set of emerging methods that select variants based on their effect on gene expressions, providing pretrained linear combinations of variants for downstream association mapping. In contrast to this, kernel methods such as sequence kernel association test (SKAT) model genotypic and phenotypic variance use various kernel functions that capture genetic similarity between subjects, allowing nonlinear effects to be included. From the perspective of machine learning, these two methods cover two complementary aspects of feature engineering: feature selection/pruning and feature aggregation. Thus far, no thorough comparison has been made between these categories, and no methods exist which incorporate the advantages of TWAS- and kernel-based methods. In this work, we developed a novel method called kernel-based TWAS (kTWAS) that applies TWAS-like feature selection to a SKAT-like kernel association test, combining the strengths of both approaches. Through extensive simulations, we demonstrate that kTWAS has higher power than TWAS and multiple SKAT-based protocols, and we identify novel disease-associated genes in Wellcome Trust Case Control Consortium genotyping array data and MSSNG (Autism) sequence data. The source code for kTWAS and our simulations are available in our GitHub repository (https://github.com/theLongLab/kTWAS).

Robust kernel association testing (RobKAT)

Testing the association between single-nucleotide polymorphism (SNP) effects and a response is often carried out through kernel machine methods based on least squares, such as the sequence kernel association test (SKAT). However, these least-squares procedures are designed for a normally distributed conditional response, which may not apply. Other robust procedures such as the quantile regression kernel machine (QRKM) restrict the choice of the loss function and only allow inference on conditional quantiles. We propose a general and robust kernel association test with a flexible choice of the loss function, no distributional assumptions, and has SKAT and QRKM as special cases. We evaluate our proposed robust association test (RobKAT) across various data distributions through a simulation study. When errors are normally distributed, RobKAT controls type I error and shows comparable power with SKAT. In all other distributional settings investigated, our robust test has similar or greater power than SKAT. Finally, we apply our robust testing method to data from the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) clinical trial to detect associations between selected genes including the major histocompatibility complex (MHC) region on chromosome six and neurotropic herpesvirus antibody levels in schizophrenia patients. RobKAT detected significant association with four SNP sets (HST1H2BJ, MHC, POM12L2, and SLC17A1), three of which were undetected by SKAT.

Reader Reaction: A note on testing and estimation in marker-set association study using semiparametric quantile regression kernel machine

Kong et al. (2016, Biometrics 72, 364-371) presented a quantile regression kernel machine (QRKM) test for robust analysis of genetic marker-set association studies. A potential limitation of QRKM is the permutation-based test design may be unscalable for the massive sizes of modern datasets. In this article, we present an alternative strategy for p-value calculation of QRKM, which is capable of speeding up the QRKM testing procedure dramatically while maintaining the same testing performance as QRKM. The effectiveness of our approach is demonstrated via simulation studies.