Meta-regression of gene-environment interaction in genome-wide association studies

Genome-Wide Gene-Environment Interaction Analysis Identifies Novel Candidate Variants for Growth Traits in Beef Cattle

Complex traits are widely considered to be the result of a compound regulation of genes, environmental factors, and genotype-by-environment interaction (G × E). The inclusion of G × E in genome-wide association analyses is essential to understand animal environmental adaptations and improve the efficiency of breeding decisions. Here, we systematically investigated the G × E of growth traits (including weaning weight, yearling weight, 18-month body weight, and 24-month body weight) with environmental factors (farm and temperature) using genome-wide genotype-by-environment interaction association studies (GWEIS) with a dataset of 1350 cattle. We validated the robust estimator's effectiveness in GWEIS and detected 29 independent interacting SNPs with a significance threshold of 1.67 × 10-6, indicating that these SNPs, which do not show main effects in traditional genome-wide association studies (GWAS), may have non-additive effects across genotypes but are obliterated by environmental means. The gene-based analysis using MAGMA identified three genes that overlapped with the GEWIS results exhibiting G × E, namely SMAD2, PALMD, and MECOM. Further, the results of functional exploration in gene-set analysis revealed the bio-mechanisms of how cattle growth responds to environmental changes, such as mitotic or cytokinesis, fatty acid β-oxidation, neurotransmitter activity, gap junction, and keratan sulfate degradation. This study not only reveals novel genetic loci and underlying mechanisms influencing growth traits but also transforms our understanding of environmental adaptation in beef cattle, thereby paving the way for more targeted and efficient breeding strategies.

Meta-Analysis of Joint Test of SNP and SNP-Environment Interaction with Heterogeneity

Many complex diseases are caused by single nucleotide polymorphisms (SNPs), environmental factors, and the interaction between SNPs and environment. Joint tests of the SNP and SNP-environment interaction effects (JMA) and meta-regression (MR) are commonly used to evaluate these SNP-environment interactions. However, these two methods do not consider genetic heterogeneity. We previously presented a random-effect MR, which provided higher power than the MR in datasets with high heterogeneity. However, this method requires group-level data, which sometimes are not available. Given this, we designed this study to evaluate the introduction of the random effects of SNP and SNP-environment interaction into the JMA, and then extended this to the random effect model. Likelihood ratio statistic is applied to test the JMA and the new method we proposed in this paper. We evaluated the null distributions of these tests, and the powers for this method. This method was verified by simulation and was shown to provide similar powers to the random effect meta-regression method (RMR). However, this method only requires study-level data which relaxed the condition of the RMR. Our study suggests that this method is more suitable for finding the association between SNP and diseases in the absence of group-level data.

Meta-analysis of SNP-environment interaction with heterogeneity for overlapping data

Meta-analysis is a popular method used in genome-wide association studies, by which the results of multiple studies are combined to identify associations. This process generates heterogeneity. Recently, we proposed a random effect model meta-regression method (MR) to study the effect of single nucleotide polymorphism (SNP)-environment interactions. This method takes heterogeneity into account and produces high power. We also proposed a fixed effect model overlapping MR in which the overlapping data is taken into account. In the present study, a random effect model overlapping MR that simultaneously considers heterogeneity and overlapping data is proposed. This method is based on the random effect model MR and the fixed effect model overlapping MR. A new way of solving the logarithm of the determinant of covariance matrices in likelihood functions is also provided. Tests for the likelihood ratio statistic of the SNP-environment interaction effect and the SNP and SNP-environment joint effects are given. In our simulations, null distributions and type I error rates were proposed to verify the suitability of our method, and powers were applied to evaluate the superiority of our method. Our findings indicate that this method is effective in cases of overlapping data with a high heterogeneity.

Meta-Analysis of SNP-Environment Interaction With Overlapping Data

Meta-analysis, which combines the results of multiple studies, is an important analytical method in genome-wide association studies. In genome-wide association studies practice, studies employing meta-analysis may have overlapping data, which could yield false positive results. Recent studies have proposed models to handle the issue of overlapping data when testing the genetic main effect of single nucleotide polymorphism. However, there is still no meta-analysis method for testing gene-environment interaction when overlapping data exist. Inspired by the methods of testing the main effect of gene with overlapping data, we proposed an overlapping meta-regulation method to address the issue in testing the gene-environment interaction. We generalized the covariance matrices of the regular meta-regression model by employing Lin’s and Han’s correlation structures to incorporate the correlations introduced by the overlapping data. Based on our proposed models, we further provided statistical significance tests of the gene-environment interaction as well as joint effects of the gene main effect and the interaction. Through simulations, we examined type I errors and statistical powers of our proposed methods at different levels of data overlap among studies. We demonstrated that our method well controls the type I error and simultaneously achieves statistical power comparable with the method that removes overlapping samples a priori before the meta-analysis, i.e., the splitting method. On the other hand, ignoring overlapping data will inflate the type I error. Unlike the splitting method that requires individual-level genotype and phenotype data, our proposed method for testing gene-environment interaction handles the issue of overlapping data effectively and statistically efficiently at the meta-analysis level.

Meta-Analysis of SNP-Environment Interaction with Heterogeneity

Meta-analyses are widely used in genome-wide association studies to combine the results obtained from multiple studies. Classical random-effects methods treat genetic heterogeneity as a random effect and consider it as a portion of the variance associated with a fixed effect of the variant. Recent work suggests performing hypothesis testing with the null hypothesis under which neither fixed nor random effects exist for a variant. This method has been shown to perform better than classical random-effects methods. In this work, we propose a meta-analysis of testing single nucleotide polymorphism (SNP)-environment interaction in the presence of genetic heterogeneity. We introduced the random effects of the SNP and SNP-environment interaction under test into a meta-regression model to account for heterogeneity. A test for the SNP-environment interaction was formulated to test for fixed and random effects of the interaction simultaneously. Similarly, a test for total genetic effects was formulated to test for fixed effects of the SNP and the SNP-environment interaction together with their random effects. We performed simulations to study the null distribution and statistical power of the proposed tests. We show that the new methods have higher power than classical random-effects and fixed-effects meta-regression methods when heterogeneity effects are large.

Robustness of meta-analyses in finding gene × environment interactions

Meta-analyses that synthesize statistical evidence across studies have become important analytical tools for genetic studies. Inspired by the success of genome-wide association studies of the genetic main effect, researchers are searching for gene × environment interactions. Confounders are routinely included in the genome-wide gene × environment interaction analysis as covariates; however, this does not control for any confounding effects on the results if covariate × environment interactions are present. We carried out simulation studies to evaluate the robustness to the covariate × environment confounder for meta-regression and joint meta-analysis, which are two commonly used meta-analysis methods for testing the gene × environment interaction or the genetic main effect and interaction jointly. Here we show that meta-regression is robust to the covariate × environment confounder while joint meta-analysis is subject to the confounding effect with inflated type I error rates. Given vast sample sizes employed in genome-wide gene × environment interaction studies, non-significant covariate × environment interactions at the study level could substantially elevate the type I error rate at the consortium level. When covariate × environment confounders are present, type I errors can be controlled in joint meta-analysis by including the covariate × environment terms in the analysis at the study level. Alternatively, meta-regression can be applied, which is robust to potential covariate × environment confounders.

Kernel Approach for Modeling Interaction Effects in Genetic Association Studies of Complex Quantitative Traits

The etiology of complex traits likely involves the effects of genetic and environmental factors, along with complicated interaction effects between them. Consequently, there has been interest in applying genetic association tests of complex traits that account for potential modification of the genetic effect in the presence of an environmental factor. One can perform such an analysis using a joint test of gene and gene-environment interaction. An optimal joint test would be one that remains powerful under a variety of models ranging from those of strong gene-environment interaction effect to those of little or no gene-environment interaction effect. To fill this demand, we have extended a kernel machine based approach for association mapping of multiple SNPs to consider joint tests of gene and gene-environment interaction. The kernel-based approach for joint testing is promising, because it incorporates linkage disequilibrium information from multiple SNPs simultaneously in analysis and permits flexible modeling of interaction effects. Using simulated data, we show that our kernel machine approach typically outperforms the traditional joint test under strong gene-environment interaction models and further outperforms the traditional main-effect association test under models of weak or no gene-environment interaction effects. We illustrate our test using genome-wide association data from the Grady Trauma Project, a cohort of highly traumatized, at-risk individuals, which has previously been investigated for interaction effects.