Finding novel genes by testing G × E interactions in a genome-wide association study

In a genome-wide association study (GWAS), investigators typically focus their primary analysis on the direct (marginal) associations of each single nucleotide polymorphism (SNP) with the trait. Some SNPs that are truly associated with the trait may not be identified in this scan if they have a weak marginal effect and thus low power to be detected. However, these SNPs may be quite important in subgroups of the population defined by an environmental or personal factor, and may be detectable if such a factor is carefully considered in a gene-environment (G × E) interaction analysis. We address the question "Using a genome wide interaction scan (GWIS), can we find new genes that were not found in the primary GWAS scan?" We review commonly used approaches for conducting a GWIS in case-control studies, and propose a new two-step screening and testing method (EDG×E) that is optimized to find genes with a weak marginal effect. We simulate several scenarios in which our two-step method provides 70-80% power to detect a disease locus while a marginal scan provides less than 5% power. We also provide simulations demonstrating that the EDG×E method outperforms other GWIS approaches (including case only and previously proposed two-step methods) for finding genes with a weak marginal effect. Application of this method to a G × Sex scan for childhood asthma reveals two potentially interesting SNPs that were not identified in the marginal-association scan. We distribute a new software program (G×Escan, available at http://biostats.usc.edu/software) that implements this new method as well as several other GWIS approaches.