Ordered subset analysis for case-control studies

Finding associations in a heterogeneous setting: statistical test for aberration enrichment

Most two-group statistical tests find broad patterns such as overall shifts in mean, median, or variance. These tests may not have enough power to detect effects in a small subset of samples, e.g., a drug that works well only on a few patients. We developed a novel statistical test targeting such effects relevant for clinical trials, biomarker discovery, feature selection, etc. We focused on finding meaningful associations in complex genetic diseases in gene expression, miRNA expression, and DNA methylation. Our test outperforms traditional statistical tests in simulated and experimental data and detects potentially disease-relevant genes with heterogeneous effects.

Detecting latent exposure in genome-wide association studies using a breakpoint model for logistic regression

Detecting gene-environment (G × E) interactions in the context of genome-wide association studies (GWAS) is a challenging problem since standard methods generally present a lack of power. An additional difficulty arises from the fact that the causal exposure is seldom observed and only a proxy of this exposure is observed. This leads to an additional drop in terms of power and it explains the failure of standard methods in detecting interactions, even very strong ones. In this article, we consider the latent exposure as a source of heterogeneity and we propose a new powerful method, named "Breakpoint Model for Logistic Regression" (BMLR), based on a breakpoint model, in order to detect G × E interactions when causal exposure is unobserved. First, the BMLR method is compared to the ordered-subset analysis for case-control method, which has been developed for the same purpose, through simulations. This highlights the ability of BMLR to detect the heterogeneity, and therefore, to detect interaction with latent exposure. Finally, the BMLR method is compared to standard methods, such as Plink, to perform a GWAS on a published realistic benchmark.

JBASE: Joint Bayesian Analysis of Subphenotypes and Epistasis

MOTIVATION

Rapid advances in genotyping and genome-wide association studies have enabled the discovery of many new genotype-phenotype associations at the resolution of individual markers. However, these associations explain only a small proportion of theoretically estimated heritability of most diseases. In this work, we propose an integrative mixture model called JBASE: joint Bayesian analysis of subphenotypes and epistasis. JBASE explores two major reasons of missing heritability: interactions between genetic variants, a phenomenon known as epistasis and phenotypic heterogeneity, addressed via subphenotyping.

RESULTS

Our extensive simulations in a wide range of scenarios repeatedly demonstrate that JBASE can identify true underlying subphenotypes, including their associated variants and their interactions, with high precision. In the presence of phenotypic heterogeneity, JBASE has higher Power and lower Type 1 Error than five state-of-the-art approaches. We applied our method to a sample of individuals from Mexico with Type 2 diabetes and discovered two novel epistatic modules, including two loci each, that define two subphenotypes characterized by differences in body mass index and waist-to-hip ratio. We successfully replicated these subphenotypes and epistatic modules in an independent dataset from Mexico genotyped with a different platform.

AVAILABILITY AND IMPLEMENTATION

JBASE is implemented in C++, supported on Linux and is available at http://www.cs.toronto.edu/∼goldenberg/JBASE/jbase.tar.gz. The genotype data underlying this study are available upon approval by the ethics review board of the Medical Centre Siglo XXI. Please contact Dr Miguel Cruz at mcruzl@yahoo.com for assistance with the application.

CONTACT

anna.goldenberg@utoronto.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ordered subset analysis of copy number variation association with age at onset of Alzheimer's disease

Genetic heterogeneity is a common problem for genome-wide association studies of complex human diseases. Ordered-subset analysis (OSA) reduces genetic heterogeneity and optimizes the use of phenotypic information, thus improving power under some disease models. We hypothesized that in a genetically heterogeneous disorder such as Alzheimer's disease (AD), utilizing OSA by age at onset (AAO) of AD may increase the power to detect relevant loci. Using this approach, 8 loci were detected, including the chr15 : 30,44 region harboring CHRFAM7A. The association was replicated in the NIA-LOAD Familial Study dataset. CHRFAM7A is a dominant negative regulator of CHRNA7 function, the receptor that facilitates amyloid-β1-42 internalization through endocytosis and has been implicated in AD. OSA, using AAO as a quantitative trait, optimized power and detected replicable signals suggesting that AD is genetically heterogeneous between AAO subsets.

The role of phenotype in gene discovery in the whole genome sequencing era

As whole genome sequence becomes a routine component of gene discovery studies in humans, we will have an exhaustive catalog of genetic variation and the challenge becomes understanding the phenotypic consequences of these variants. Statistical genetic methods and analytical approaches that are concerned with optimizing phenotypes for gene discovery for complex traits offer two general categories of advantages. They may increase power to localize genes of interest and also aid in interpreting associations between genetic variants and disease outcomes by suggesting potential mechanisms and pathways through which genes may affect outcomes. Such phenotype optimization approaches include use of allied phenotypes such as symptoms or ages of onset to reduce genetic heterogeneity within a set of cases, study of quantitative risk factors or endophenotypes, joint analyses of related phenotypes, and derivation of new phenotypes designed to extract independent measures underlying the correlations among a set of related phenotypes through approaches such as principal components. New opportunities are also presented by technological advances that permit efficient collection of hundreds or thousands of phenotypes on an individual, including phenotypes more proximal to the level of gene action such as levels of gene expression, microRNAs, or metabolic and proteomic profiles.

TDT-HET: a new transmission disequilibrium test that incorporates locus heterogeneity into the analysis of family-based association data

BACKGROUND

Locus heterogeneity is one of the most documented phenomena in genetics. To date, relatively little work had been done on the development of methods to address locus heterogeneity in genetic association analysis. Motivated by Zhou and Pan's work, we present a mixture model of linked and unlinked trios and develop a statistical method to estimate the probability that a heterozygous parent transmits the disease allele at a di-allelic locus, and the probability that any trio is in the linked group. The purpose here is the development of a test that extends the classic transmission disequilibrium test (TDT) to one that accounts for locus heterogeneity.

RESULTS

Our simulations suggest that, for sufficiently large sample size (1000 trios) our method has good power to detect association even the proportion of unlinked trios is high (75%). While the median difference (TDT-HET empirical power - TDT empirical power) is approximately 0 for all MOI, there are parameter settings for which the power difference can be substantial. Our multi-locus simulations suggest that our method has good power to detect association as long as the markers are reasonably well-correlated and the genotype relative risk are larger. Results of both single-locus and multi-locus simulations suggest our method maintains the correct type I error rate.Finally, the TDT-HET statistic shows highly significant p-values for most of the idiopathic scoliosis candidate loci, and for some loci, the estimated proportion of unlinked trios approaches or exceeds 50%, suggesting the presence of locus heterogeneity.

CONCLUSIONS

We have developed an extension of the TDT statistic (TDT-HET) that allows for locus heterogeneity among coded trios. Benefits of our method include: estimates of parameters in the presence of heterogeneity, and reasonable power even when the proportion of linked trios is small. Also, we have extended multi-locus methods to TDT-HET and have demonstrated that the empirical power may be high to detect linkage. Last, given that we obtain PPBs, we conjecture that the TDT-HET may be a useful method for correctly identifying linked trios. We anticipate that researchers will find this property increasingly useful as they apply next-generation sequencing data in family based studies.

IGF-I gene variability is associated with an increased risk for AD

Insulin-like growth factor I (IGF-I), a neuroprotective factor with a wide spectrum of actions in the adult brain, is involved in the pathogenesis of Alzheimer's disease (AD). Circulating levels of IGF-I change in AD patients and are implicated in the clearance of brain amyloid beta (Aβ) complexes. To investigate this hypothesis, we screened the IGF-I gene for various well known single nucleotide polymorphisms (SNPs) covering % of the gene variability in a population of 2352 individuals. Genetic analysis indicated different distribution of genotypes of 1 single nucleotide polymorphism, and 1 extended haplotype in the AD population compared with healthy control subjects. In particular, the frequency of rs972936 GG genotype was significantly greater in AD patients than in control subjects (63% vs. 55%). The rs972936 GG genotype was associated with an increased risk for disease, independently of apolipoprotein E genotype, and with enhanced circulating levels of IGF-I. These findings suggest that polymorphisms within the IGF-I gene could infer greater risk for AD through their effect on IGF-I levels, and confirm the physiological role IGF-I in the pathogenesis of AD.