Sample size requirements to detect gene-environment interactions in genome-wide association studies

Is the gene-environment interaction paradigm relevant to genome-wide studies? The case of education and body mass index

This study uses data from the Framingham Heart Study to examine the relevance of the gene-environment interaction paradigm for genome-wide association studies (GWAS). We use completed college education as our environmental measure and estimate the interactive effect of genotype and education on body mass index (BMI) using 260,402 single-nucleotide polymorphisms (SNPs). Our results highlight the sensitivity of parameter estimates obtained from GWAS models and the difficulty of framing genome-wide results using the existing gene-environment interaction typology. We argue that SNP-environment interactions across the human genome are not likely to provide consistent evidence regarding genetic influences on health that differ by environment. Nevertheless, genome-wide data contain rich information about individual respondents, and we demonstrate the utility of this type of data. We highlight the fact that GWAS is just one use of genome-wide data, and we encourage demographers to develop methods that incorporate this vast amount of information from respondents into their analyses.

Genomewide interaction and enrichment analysis on antidepressant response

BACKGROUND

Genomewide association studies (GWASs) on antidepressant efficacy have yielded modest results. A possible reason is that response is influenced by other factors, which possibly interact with genetic variation. We used a GWAS model to predict antidepressant response, by including predictors previously known to affect response, such as quality of life (QoL). We also evaluated the association between genes, previously implicated in gene-environment (G × E) interactions, and response using an enrichment analysis.

METHOD

We examined a sample of 1426 depressed patients from the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) trial: 774 responders, 652 non-responders and 418,865 single nucleotide polymorphisms (SNPs) were analysed. First, in a GWAS model, we investigated whether genetic variations interact with patients' levels of QoL to predict response, after controlling for demographic characteristics, severity and population stratification. Second, we conducted an enrichment analysis exploring whether candidate genes that have emerged from prior G × E interaction studies on depression are associated with treatment response.

RESULTS

The GWAS model, with QoL as a moderator, yielded one SNP (rs520210) associated with response in the NEDD4L gene (p = 3.64 × 10⁻⁸). In the Caucasian sample only, we observed a drop in significance for this SNP. The enrichment analysis showed that SNPs within serotonergic genes contained more significant markers that predicted response, compared with a random set of genes in the genome.

CONCLUSIONS

Our findings point to possible target genes, which are proposed for further independent replication. Our enrichment analysis provides further support, in a genomewide context, of the role of serotonergic genes in influencing antidepressant response.

Genome-wide association studies and the clinic: a focus on breast cancer

Breast cancer is the most frequently diagnosed cancer among women worldwide, and has long been considered to be a genetic disease. A wide range of genetic variants, both rare mutations and more common variants, have been shown to influence breast cancer risk. In particular, recent studies have identified a number of common genetic variants, or single nucleotide polymorphisms, that are associated with breast cancer risk. In this review, we will briefly present the genetic epidemiology of breast cancer, genome-wide association study technology and how this technology may influence breast cancer screening in the clinic.

Identification of new genetic susceptibility loci for breast cancer through consideration of gene-environment interactions

Genes that alter disease risk only in combination with certain environmental exposures may not be detected in genetic association analysis. By using methods accounting for gene-environment (G × E) interaction, we aimed to identify novel genetic loci associated with breast cancer risk. Up to 34,475 cases and 34,786 controls of European ancestry from up to 23 studies in the Breast Cancer Association Consortium were included. Overall, 71,527 single nucleotide polymorphisms (SNPs), enriched for association with breast cancer, were tested for interaction with 10 environmental risk factors using three recently proposed hybrid methods and a joint test of association and interaction. Analyses were adjusted for age, study, population stratification, and confounding factors as applicable. Three SNPs in two independent loci showed statistically significant association: SNPs rs10483028 and rs2242714 in perfect linkage disequilibrium on chromosome 21 and rs12197388 in ARID1B on chromosome 6. While rs12197388 was identified using the joint test with parity and with age at menarche (P-values = 3 × 10(-07)), the variants on chromosome 21 q22.12, which showed interaction with adult body mass index (BMI) in 8,891 postmenopausal women, were identified by all methods applied. SNP rs10483028 was associated with breast cancer in women with a BMI below 25 kg/m(2) (OR = 1.26, 95% CI 1.15-1.38) but not in women with a BMI of 30 kg/m(2) or higher (OR = 0.89, 95% CI 0.72-1.11, P for interaction = 3.2 × 10(-05)). Our findings confirm comparable power of the recent methods for detecting G × E interaction and the utility of using G × E interaction analyses to identify new susceptibility loci.

Gene-environment interactions in cancer epidemiology: a National Cancer Institute Think Tank report

Cancer risk is determined by a complex interplay of genetic and environmental factors. Genome-wide association studies (GWAS) have identified hundreds of common (minor allele frequency [MAF] > 0.05) and less common (0.01 < MAF < 0.05) genetic variants associated with cancer. The marginal effects of most of these variants have been small (odds ratios: 1.1-1.4). There remain unanswered questions on how best to incorporate the joint effects of genes and environment, including gene-environment (G × E) interactions, into epidemiologic studies of cancer. To help address these questions, and to better inform research priorities and allocation of resources, the National Cancer Institute sponsored a "Gene-Environment Think Tank" on January 10-11, 2012. The objective of the Think Tank was to facilitate discussions on (1) the state of the science, (2) the goals of G × E interaction studies in cancer epidemiology, and (3) opportunities for developing novel study designs and analysis tools. This report summarizes the Think Tank discussion, with a focus on contemporary approaches to the analysis of G × E interactions. Selecting the appropriate methods requires first identifying the relevant scientific question and rationale, with an important distinction made between analyses aiming to characterize the joint effects of putative or established genetic and environmental factors and analyses aiming to discover novel risk factors or novel interaction effects. Other discussion items include measurement error, statistical power, significance, and replication. Additional designs, exposure assessments, and analytical approaches need to be considered as we move from the current small number of success stories to a fuller understanding of the interplay of genetic and environmental factors.

Genes-environment interactions in obesity- and diabetes-associated pancreatic cancer: a GWAS data analysis

BACKGROUND

Obesity and diabetes are potentially alterable risk factors for pancreatic cancer. Genetic factors that modify the associations of obesity and diabetes with pancreatic cancer have previously not been examined at the genome-wide level.

METHODS

Using genome-wide association studies (GWAS) genotype and risk factor data from the Pancreatic Cancer Case Control Consortium, we conducted a discovery study of 2,028 cases and 2,109 controls to examine gene-obesity and gene-diabetes interactions in relation to pancreatic cancer risk by using the likelihood-ratio test nested in logistic regression models and Ingenuity Pathway Analysis (IPA).

RESULTS

After adjusting for multiple comparisons, a significant interaction of the chemokine signaling pathway with obesity (P = 3.29 × 10(-6)) and a near significant interaction of calcium signaling pathway with diabetes (P = 1.57 × 10(-4)) in modifying the risk of pancreatic cancer were observed. These findings were supported by results from IPA analysis of the top genes with nominal interactions. The major contributing genes to the two top pathways include GNGT2, RELA, TIAM1, and GNAS. None of the individual genes or single-nucleotide polymorphism (SNP) except one SNP remained significant after adjusting for multiple testing. Notably, SNP rs10818684 of the PTGS1 gene showed an interaction with diabetes (P = 7.91 × 10(-7)) at a false discovery rate of 6%.

CONCLUSIONS

Genetic variations in inflammatory response and insulin resistance may affect the risk of obesity- and diabetes-related pancreatic cancer. These observations should be replicated in additional large datasets.

IMPACT

A gene-environment interaction analysis may provide new insights into the genetic susceptibility and molecular mechanisms of obesity- and diabetes-related pancreatic cancer.

Gene-environment interactions in genome-wide association studies: current approaches and new directions

BACKGROUND

Complex psychiatric traits have long been thought to be the result of a combination of genetic and environmental factors, and gene-environment interactions are thought to play a crucial role in behavioral phenotypes and the susceptibility and progression of psychiatric disorders. Candidate gene studies to investigate hypothesized gene-environment interactions are now fairly common in human genetic research, and with the shift toward genome-wide association studies, genome-wide gene-environment interaction studies are beginning to emerge.

METHODS

We summarize the basic ideas behind gene-environment interaction, and provide an overview of possible study designs and traditional analysis methods in the context of genome-wide analysis. We then discuss novel approaches beyond the traditional strategy of analyzing the interaction between the environmental factor and each polymorphism individually.

RESULTS

Two-step filtering approaches that reduce the number of polymorphisms tested for interactions can substantially increase the power of genome-wide gene-environment studies. New analytical methods including data-mining approaches, and gene-level and pathway-level analyses, also have the capacity to improve our understanding of how complex genetic and environmental factors interact to influence psychologic and psychiatric traits. Such methods, however, have not yet been utilized much in behavioral and mental health research.

CONCLUSIONS

Although methods to investigate gene-environment interactions are available, there is a need for further development and extension of these methods to identify gene-environment interactions in the context of genome-wide association studies. These novel approaches need to be applied in studies of psychology and psychiatry.

Comparisons of power of statistical methods for gene-environment interaction analyses

Any genome-wide analysis is hampered by reduced statistical power due to multiple comparisons. This is particularly true for interaction analyses, which have lower statistical power than analyses of associations. To assess gene-environment interactions in population settings we have recently proposed a statistical method based on a modified two-step approach, where first genetic loci are selected by their associations with disease and environment, respectively, and subsequently tested for interactions. We have simulated various data sets resembling real world scenarios and compared single-step and two-step approaches with respect to true positive rate (TPR) in 486 scenarios and (study-wide) false positive rate (FPR) in 252 scenarios. Our simulations confirmed that in all two-step methods the two steps are not correlated. In terms of TPR, two-step approaches combining information on gene-disease association and gene-environment association in the first step were superior to all other methods, while preserving a low FPR in over 250 million simulations under the null hypothesis. Our weighted modification yielded the highest power across various degrees of gene-environment association in the controls. An optimal threshold for step 1 depended on the interacting allele frequency and the disease prevalence. In all scenarios, the least powerful method was to proceed directly to an unbiased full interaction model, applying conventional genome-wide significance thresholds. This simulation study confirms the practical advantage of two-step approaches to interaction testing over more conventional one-step designs, at least in the context of dichotomous disease outcomes and other parameters that might apply in real-world settings.

Test for interactions between a genetic marker set and environment in generalized linear models

We consider in this paper testing for interactions between a genetic marker set and an environmental variable. A common practice in studying gene-environment (GE) interactions is to analyze one single-nucleotide polymorphism (SNP) at a time. It is of significant interest to analyze SNPs in a biologically defined set simultaneously, e.g. gene or pathway. In this paper, we first show that if the main effects of multiple SNPs in a set are associated with a disease/trait, the classical single SNP-GE interaction analysis can be biased. We derive the asymptotic bias and study the conditions under which the classical single SNP-GE interaction analysis is unbiased. We further show that, the simple minimum p-value-based SNP-set GE analysis, can be biased and have an inflated Type 1 error rate. To overcome these difficulties, we propose a computationally efficient and powerful gene-environment set association test (GESAT) in generalized linear models. Our method tests for SNP-set by environment interactions using a variance component test, and estimates the main SNP effects under the null hypothesis using ridge regression. We evaluate the performance of GESAT using simulation studies, and apply GESAT to data from the Harvard lung cancer genetic study to investigate GE interactions between the SNPs in the 15q24-25.1 region and smoking on lung cancer risk.

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.

Empirical hierarchical bayes approach to gene-environment interactions: development and application to genome-wide association studies of lung cancer in TRICL

The analysis of gene-environment (G × E) interactions remains one of the greatest challenges in the postgenome-wide association studies (GWASs) era. Recent methods constitute a compromise between the robust but underpowered case-control and powerful case-only methods. Inferences of the latter are biased when the assumption of gene-environment (G-E) independence in controls fails. We propose a novel empirical hierarchical Bayes approach to G × E interaction (EHB-GE), which benefits from greater rank power while accounting for population-based G-E correlation. Building on Lewinger et al.'s ([2007] Genet Epidemiol 31:871-882) hierarchical Bayes prioritization approach, the method first obtains posterior G-E correlation estimates in controls for each marker, borrowing strength from G-E information across the genome. These posterior estimates are then subtracted from the corresponding case-only G × E estimates. We compared EHB-GE with rival methods using simulation. EHB-GE has similar or greater rank power to detect G × E interactions in the presence of large numbers of G-E correlations with weak to strong effects or only a low number of such correlations with large effect. When there are no or only a few weak G-E correlations, Murcray et al.'s method ([2009] Am J Epidemiol 169:219-226) identifies markers with low G × E interaction effects better. We applied EHB-GE and competing methods to four lung cancer case-control GWAS from the Interdisciplinary Research in Cancer of the Lung/International Lung Cancer Consortium with smoking as environmental factor. A number of genes worth investigating were identified by the EHB-GE approach.

SBERIA: set-based gene-environment interaction test for rare and common variants in complex diseases

Identification of gene-environment interaction (G × E) is important in understanding the etiology of complex diseases. However, partially due to the lack of power, there have been very few replicated G × E findings compared to the success in marginal association studies. The existing G × E testing methods mainly focus on improving the power for individual markers. In this paper, we took a different strategy and proposed a set-based gene-environment interaction test (SBERIA), which can improve the power by reducing the multiple testing burdens and aggregating signals within a set. The major challenge of the signal aggregation within a set is how to tell signals from noise and how to determine the direction of the signals. SBERIA takes advantage of the established correlation screening for G × E to guide the aggregation of genotypes within a marker set. The correlation screening has been shown to be an efficient way of selecting potential G × E candidate SNPs in case-control studies for complex diseases. Importantly, the correlation screening in case-control combined samples is independent of the interaction test. With this desirable feature, SBERIA maintains the correct type I error level and can be easily implemented in a regular logistic regression setting. We showed that SBERIA had higher power than benchmark methods in various simulation scenarios, both for common and rare variants. We also applied SBERIA to real genome-wide association studies (GWAS) data of 10,729 colorectal cancer cases and 13,328 controls and found evidence of interaction between the set of known colorectal cancer susceptibility loci and smoking.

Efficient two-step testing of gene-gene interactions in genome-wide association studies

Exhaustive testing of all possible SNP pairs in a genome-wide association study (GWAS) generally yields low power to detect gene-gene (G × G) interactions because of small effect sizes and stringent requirements for multiple-testing correction. We introduce a new two-step procedure for testing G × G interactions in case-control GWAS to detect interacting single nucleotide polymorphisms (SNPs) regardless of their marginal effects. In an initial screening step, all SNP pairs are tested for gene-gene association in the combined sample of cases and controls. In the second step, the pairs that pass the screening are followed up with a traditional test for G × G interaction. We show that the two-step method is substantially more powerful to detect G × G interactions than the exhaustive testing approach. For example, with 2,000 cases and 2,000 controls, the two-step method can have more than 90% power to detect an interaction odds ratio of 2.0 compared to less than 50% power for the exhaustive testing approach. Moreover, we show that a hybrid two-step approach that combines our newly proposed two-step test and the two-step test that screens for marginal effects retains the best power properties of both. The two-step procedures we introduce have the potential to uncover genetic signals that have not been previously identified in an initial single-SNP GWAS. We demonstrate the computational feasibility of the two-step G × G procedure by performing a G × G scan in the asthma GWAS of the University of Southern California Children's Health Study.

The case-only test for gene-environment interaction is not uniformly powerful: an empirical example

The case-only test has been proposed as a more powerful approach to detect gene-environment (G × E) interactions. This approach assumes that the genetic and environmental factors are independent. Although it is well known that Type I error rate will increase if this assumption is violated, it is less widely appreciated that G × E correlation can also lead to power loss. We illustrate this phenomenon by comparing the performance of the case-only test to other approaches to detect G × E interactions in a genome-wide association study (GWAS) of esophageal squamous-cell carcinoma (ESCC) in Chinese populations. Some of these approaches do not use information on the correlation between exposure and genotype (standard logistic regression), whereas others seek to use this information in a robust fashion to boost power without increasing Type I error (two-step, empirical Bayes, and cocktail methods). G × E interactions were identified involving drinking status and two regions containing genes in the alcohol metabolism pathway, 4q23 and 12q24. Although the case-only test yielded the most significant tests of G × E interaction in the 4q23 region, the case-only test failed to identify significant interactions in the 12q24 region which were readily identified using other approaches. The low power of the case-only test in the 12q24 region is likely due to the strong inverse association between the single nucleotide polymorphism (SNPs) in this region and drinking status. This example underscores the need to consider multiple approaches to detect G × E interactions, as different tests are more or less sensitive to different alternative hypotheses and violations of the G × E independence assumption.

Sample size considerations of prediction-validation methods in high-dimensional data for survival outcomes

A variety of prediction methods are used to relate high-dimensional genome data with a clinical outcome using a prediction model. Once a prediction model is developed from a data set, it should be validated using a resampling method or an independent data set. Although the existing prediction methods have been intensively evaluated by many investigators, there has not been a comprehensive study investigating the performance of the validation methods, especially with a survival clinical outcome. Understanding the properties of the various validation methods can allow researchers to perform more powerful validations while controlling for type I error. In addition, sample size calculation strategy based on these validation methods is lacking. We conduct extensive simulations to examine the statistical properties of these validation strategies. In both simulations and a real data example, we have found that 10-fold cross-validation with permutation gave the best power while controlling type I error close to the nominal level. Based on this, we have also developed a sample size calculation method that will be used to design a validation study with a user-chosen combination of prediction. Microarray and genome-wide association studies data are used as illustrations. The power calculation method in this presentation can be used for the design of any biomedical studies involving high-dimensional data and survival outcomes.

Gene-environment interactions in asthma and allergic diseases: challenges and perspectives

The concept of gene-environment (GxE) interactions has dramatically evolved in the last century and has now become a central theme in studies that assess the causes of human disease. Despite the numerous efforts to discover genes associated in asthma and allergy through various approaches, including the recent genome-wide association studies, investigation of GxE interactions has been mainly limited to candidate genes, candidate environmental exposures, or both. This review discusses the various strategies from hypothesis-driven strategies to the full agnostic search of GxE interactions with an illustration from recently published articles. Challenges raised by each piece of the puzzle (ie, phenotype, environment, gene, and analysis of GxE interaction) are put forward, and tentative solutions are proposed. New perspectives to integrate various types of data generated by new sequencing technologies and to progress toward a systems biology approach of disease are outlined. The future of a molecular network-based approach of disease to which GxE interactions are related requires space for innovative and multidisciplinary research. Assembling the various parts of a puzzle in a complex system could well occur in a way that might not necessarily follow the rules of logic.

Genome-wide joint meta-analysis of SNP and SNP-by-smoking interaction identifies novel loci for pulmonary function

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

The TNF-α -308 Promoter Gene Polymorphism and Chronic HBV Infection

Background and Aims. TNF-α -308 allele promoter polymorphism has been known to be a potential prognostic factor in patients with chronic HBV infection. We tried to determine how TNF-α -308 allele promoter polymorphism would affect the prognosis in patients with chronic HBV infection. Methods. We searched MEDLINE, EMBASE, and reference lists of relevant review articles related to the association between “TNF-α G-308A promoter polymorphism” with “chronic HBV infection”. We only focused on searching -308 locus in published studies. We reviewed 21 original articles about TNF-α -308 allele polymorphism and its effect on prognosis in patients with chronic HBV infection and discussed the results. Results. conflicting results were observed. The results were divided into 3 groups including neutral, negative, and positive associations between TNF-α -308 allele polymorphism and prognosis in patients with chronic HBV infection. We summarized the primary data as a table. Conclusions. Authors concluded that although there is an upward trend in evidence to claim that there is a positive relation between TNF-α G-308A promoter polymorphisms and resolution of chronic HBV infection, due to many biases and limitations observed in reviewed studies, an organized well-designed study is needed for clarifying the real association.

Powerful cocktail methods for detecting genome-wide gene-environment interaction

Identifying gene and environment interaction (G × E) can provide insights into biological networks of complex diseases, identify novel genes that act synergistically with environmental factors, and inform risk prediction. However, despite the fact that hundreds of novel disease-associated loci have been identified from genome-wide association studies (GWAS), few G × Es have been discovered. One reason is that most studies are underpowered for detecting these interactions. Several new methods have been proposed to improve power for G × E analysis, but performance varies with scenario. In this article, we present a module-based approach to integrating various methods that exploits each method's most appealing aspects. There are three modules in our approach: (1) a screening module for prioritizing Single Nucleotide Polymorphisms (SNPs); (2) a multiple comparison module for testing G × E; and (3) a G × E testing module. We combine all three of these modules and develop two novel "cocktail" methods. We demonstrate that the proposed cocktail methods maintain the type I error, and that the power tracks well with the best existing methods, despite that the best methods may be different under various scenarios and interaction models. For GWAS, where the true interaction models are unknown, methods like our "cocktail" methods that are powerful under a wide range of situations are particularly appealing. Broadly speaking, the modular approach is conceptually straightforward and computationally simple. It builds on common test statistics and is easily implemented without additional computational efforts. It also allows for an easy incorporation of new methods as they are developed. Our work provides a comprehensive and powerful tool for devising effective strategies for genome-wide detection of gene-environment interactions.

Two-stage testing procedures with independent filtering for genome-wide gene-environment interaction

Several two-stage multiple testing procedures have been proposed to detect gene-environment interaction in genome-wide association studies. In this article, we elucidate general conditions that are required for validity and power of these procedures, and we propose extensions of two-stage procedures using the case-only estimator of gene-treatment interaction in randomized clinical trials. We develop a unified estimating equation approach to proving asymptotic independence between a filtering statistic and an interaction test statistic in a range of situations, including marginal association and interaction in a generalized linear model with a canonical link. We assess the performance of various two-stage procedures in simulations and in genetic studies from Women's Health Initiative clinical trials.

Age of onset and effect size in genome-wide association studies

BACKGROUND

Genome-wide association studies (GWAS) have identified many susceptibility loci for complex traits, but have not identified the majority of the genetic contribution to common diseases. We explored whether the magnitude of associations detected in GWAS and, therefore, the likelihood of detecting a significant association for a given sample size, is generally greater for childhood-onset traits (e.g., birth defects) than for traits with onset in adulthood.

METHODS

Data were obtained from the National Human Genome Research Institute Catalog of Published GWAS. Traits were categorized as having an average age of onset in childhood (<18 years, n = 15 traits), early adulthood (18-54 years, n = 32 traits), or late adulthood (≥55 years, n = 31 traits). The relationship between age of onset category and the magnitude of significant associations detected by GWAS was assessed using logistic regression.

RESULTS

Associations characterized by an odds ratio (OR) ≥ 1.5 were significantly more common for GWAS of childhood traits than for late adulthood-onset traits after adjustment for several covariates (adjusted OR, 2.55; 95% confidence interval, 1.37-4.73). Results in subgroup analyses using more stringent inclusion criteria (based on sample size, effect size, p value threshold for inclusion, and novel variant-trait associations) were similar.

CONCLUSIONS

These findings suggest that, on average, marker-trait associations detected in GWAS for traits with young onset may have a larger magnitude of effect than those for traits with adult onset. Therefore, GWAS for young-onset traits, such as birth defects, may be more likely than those for adult-onset traits to identify major genetic risk factors. Birth Defects Research (Part A) 2012. © 2012 Wiley Periodicals, Inc.

Genomics and the respiratory effects of air pollution exposure

Adverse health effects from air pollutants remain important, despite improvement in air quality in the past few decades. The exact mechanisms of lung injury from exposure to air pollutants are not yet fully understood. Studying the genome (e.g. single-nucleotide polymorphisms (SNP) ), epigenome (e.g. methylation of genes), transcriptome (mRNA expression) and microRNAome (microRNA expression) has the potential to improve our understanding of the adverse effects of air pollutants. Genome-wide association studies of SNP have detected SNP associated with respiratory phenotypes; however, to date, only candidate gene studies of air pollution exposure have been performed. Changes in epigenetic processes, such DNA methylation that leads to gene silencing without altering the DNA sequence, occur with air pollutant exposure, especially global and gene-specific methylation changes. Respiratory cell line and animal models demonstrate distinct gene expression signatures in the transcriptome, arising from exposure to particulate matter or ozone. Particulate matter and other environmental toxins alter expression of microRNA, which are short non-coding RNA that regulate gene expression. While it is clearly important to contain rising levels of air pollution, strategies also need to be developed to minimize the damaging effects of air pollutant exposure on the lung, especially for patients with chronic lung disease and for people at risk of future lung disease. Careful study of genomic responses will improve our understanding of mechanisms of lung injury from air pollution and enable future clinical testing of interventions against the toxic effects of air pollutants.

Beyond the fourth wave of genome-wide obesity association studies

Obesity and related complications are major health burdens. Almost 700 million adults are currently obese globally and the prevalence is predicted to rise towards 2030. The sudden change of lifestyle with physical inactivity and excessive calorie intake undoubtedly have a major part of the epidemic development; however, some individuals seem to be more prone to be affected by an unhealthy lifestyle than others. Hence, genetic predisposition also has an essential role in determining disease susceptibility and response to lifestyle factors. Since the introduction of genome-wide association studies (GWAS), the success of identifying obesity susceptibility variants have increased, and a total of 32 variants have been identified associating genome-wide significantly with body mass index (BMI) and 18 with measures of fat distribution during four overall obesity GWAS waves. However, the immediate success of the GWAS approach has eased off, but the proportion of explained variance for BMI by the identified obesity variants remains low. This review suggests and discusses new initiatives to take GWAS of obesity to the next level, including gene–environment interactions as modulating/masking factors, low-frequent or rare variants and ways to address such analyses, and finally reflections about the applicability of epigenetic modifications when elucidating the genetic background of obesity.

Simultaneously testing for marginal genetic association and gene-environment interaction

In this article, the authors propose to simultaneously test for marginal genetic association and gene-environment interaction to discover single nucleotide polymorphisms that may be involved in gene-environment or gene-treatment interaction. The asymptotic independence of the marginal association estimator and various interaction estimators leads to a simple and flexible way of combining the 2 tests, allowing for exploitation of gene-environment independence in estimating gene-environment interaction. The proposed test differs from the 2-df test proposed by Kraft et al. (Hum Hered. 2007;63(2):111-119) in two respects. First, for the genetic association component, it tests for marginal association, which is often the primary objective in inference, rather than the main effect in a model with gene-environment interaction. Second, the gene-environment testing component can easily exploit putative gene-environment independence using either the case-only estimator or the empirical Bayes estimator, depending on whether the goal is gene-treatment interaction in a randomized trial or gene-environment interaction in an observational study. The use of the proposed joint test is illustrated through simulations and a genetic study (1993-2005) from the Women's Health Initiative.

Challenges and opportunities in genome-wide environmental interaction (GWEI) studies

The interest in performing gene-environment interaction studies has seen a significant increase with the increase of advanced molecular genetics techniques. Practically, it became possible to investigate the role of environmental factors in disease risk and hence to investigate their role as genetic effect modifiers. The understanding that genetics is important in the uptake and metabolism of toxic substances is an example of how genetic profiles can modify important environmental risk factors to disease. Several rationales exist to set up gene-environment interaction studies and the technical challenges related to these studies-when the number of environmental or genetic risk factors is relatively small-has been described before. In the post-genomic era, it is now possible to study thousands of genes and their interaction with the environment. This brings along a whole range of new challenges and opportunities. Despite a continuing effort in developing efficient methods and optimal bioinformatics infrastructures to deal with the available wealth of data, the challenge remains how to best present and analyze genome-wide environmental interaction (GWEI) studies involving multiple genetic and environmental factors. Since GWEIs are performed at the intersection of statistical genetics, bioinformatics and epidemiology, usually similar problems need to be dealt with as for genome-wide association gene-gene interaction studies. However, additional complexities need to be considered which are typical for large-scale epidemiological studies, but are also related to "joining" two heterogeneous types of data in explaining complex disease trait variation or for prediction purposes.

A general framework for two-stage analysis of genome-wide association studies and its application to case-control studies

Two-stage analyses of genome-wide association studies have been proposed as a means to improving power for designs including family-based association and gene-environment interaction testing. In these analyses, all markers are first screened via a statistic that may not be robust to an underlying assumption, and the markers thus selected are then analyzed in a second stage with a test that is independent from the first stage and is robust to the assumption in question. We give a general formulation of two-stage designs and show how one can use this formulation both to derive existing methods and to improve upon them, opening up a range of possible further applications. We show how using simple regression models in conjunction with external data such as average trait values can improve the power of genome-wide association studies. We focus on case-control studies and show how it is possible to use allele frequencies derived from an external reference to derive a powerful two-stage analysis. An illustration involving the Wellcome Trust Case-Control Consortium data shows several genome-wide-significant associations, subsequently validated, that were not significant in the standard analysis. We give some analytic properties of the methods and discuss some underlying principles.

Invited commentary: GE-Whiz! Ratcheting gene-environment studies up to the whole genome and the whole exposome

One goal in the post-genome-wide association study era is characterizing gene-environment interactions, including scanning for interactions with all available polymorphisms, not just those showing significant main effects. In recent years, several approaches to such "gene-environment-wide interaction studies" have been proposed. Two contributions in this issue of the American Journal of Epidemiology provide systematic comparisons of the performance of these various approaches, one based on simulation and one based on application to 2 real genome-wide association study scans for type 2 diabetes. The authors discuss some of the broader issues raised by these contributions, including the plausibility of the gene-environment independence assumption that some of these approaches rely upon, the need for replication, and various generalizations of these approaches.

Gene-environment interactions in genome-wide association studies: a comparative study of tests applied to empirical studies of type 2 diabetes

The question of which statistical approach is the most effective for investigating gene-environment (G-E) interactions in the context of genome-wide association studies (GWAS) remains unresolved. By using 2 case-control GWAS (the Nurses' Health Study, 1976-2006, and the Health Professionals Follow-up Study, 1986-2006) of type 2 diabetes, the authors compared 5 tests for interactions: standard logistic regression-based case-control; case-only; semiparametric maximum-likelihood estimation of an empirical-Bayes shrinkage estimator; and 2-stage tests. The authors also compared 2 joint tests of genetic main effects and G-E interaction. Elevated body mass index was the exposure of interest and was modeled as a binary trait to avoid an inflated type I error rate that the authors observed when the main effect of continuous body mass index was misspecified. Although both the case-only and the semiparametric maximum-likelihood estimation approaches assume that the tested markers are independent of exposure in the general population, the authors did not observe any evidence of inflated type I error for these tests in their studies with 2,199 cases and 3,044 controls. Both joint tests detected markers with known marginal effects. Loci with the most significant G-E interactions using the standard, empirical-Bayes, and 2-stage tests were strongly correlated with the exposure among controls. Study findings suggest that methods exploiting G-E independence can be efficient and valid options for investigating G-E interactions in GWAS.

Testing gene-environment interaction in large-scale case-control association studies: possible choices and comparisons

Several methods for screening gene-environment interaction have recently been proposed that address the issue of using gene-environment independence in a data-adaptive way. In this report, the authors present a comparative simulation study of power and type I error properties of 3 classes of procedures: 1) the standard 1-step case-control method; 2) the case-only method that requires an assumption of gene-environment independence for the underlying population; and 3) a variety of hybrid methods, including empirical-Bayes, 2-step, and model averaging, that aim at gaining power by exploiting the assumption of gene-environment independence and yet can protect against false positives when the independence assumption is violated. These studies suggest that, although the case-only method generally has maximum power, it has the potential to create substantial false positives in large-scale studies even when a small fraction of markers are associated with the exposure under study in the underlying population. All the hybrid methods perform well in protecting against such false positives and yet can retain substantial power advantages over standard case-control tests. The authors conclude that, for future genome-wide scans for gene-environment interactions, major power gain is possible by using alternatives to standard case-control analysis. Whether a case-only type scan or one of the hybrid methods should be used depends on the strength and direction of gene-environment interaction and association, the level of tolerance for false positives, and the nature of replication strategies.