INTRODUCTION: Genetic discoveries are validated through the meta-analysis of genome-wide association scans in large international consortia. Because environmental variables may interact with genetic factors, investigation of differing genetic effects for distinct levels of an environmental exposure in these large consortia may yield additional susceptibility loci undetected by main effects analysis. We describe a method of joint meta-analysis (JMA) of SNP and SNP by Environment (SNP × E) regression coefficients for use in gene-environment interaction studies. METHODS: In testing SNP × E interactions, one approach uses a two degree of freedom test to identify genetic variants that influence the trait of interest. This approach detects both main and interaction effects between the trait and the SNP. We propose a method to jointly meta-analyze the SNP and SNP × E coefficients using multivariate generalized least squares. This approach provides confidence intervals of the two estimates, a joint significance test for SNP and SNP × E terms, and a test of homogeneity across samples. RESULTS: We present a simulation study comparing this method to four other methods of meta-analysis and demonstrate that the JMA performs better than the others when both main and interaction effects are present. Additionally, we implemented our methods in a meta-analysis of the association between SNPs from the type 2 diabetes-associated gene PPARG and log-transformed fasting insulin levels and interaction by body mass index in a combined sample of 19,466 individuals from five cohorts.