The challenge of detecting genotype-by-methylation interaction: GAW20.

BACKGROUND: GAW20 working group 5 brought together researchers who contributed 7 papers with the aim of evaluating methods to detect genetic by epigenetic interactions. GAW20 distributed real data from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study, including single-nucleotide polymorphism (SNP) markers, methylation (cytosine-phosphate-guanine [CpG]) markers, and phenotype information on up to 995 individuals. In addition, a simulated data set based on the real data was provided. RESULTS: The 7 contributed papers analyzed these data sets with a number of different statistical methods, including generalized linear mixed models, mediation analysis, machine learning, W-test, and sparsity-inducing regularized regression. These methods generally appeared to perform well. Several papers confirmed a number of causative SNPs in either the large number of simulation sets or the real data on chromosome 11. Findings were also reported for different SNPs, CpG sites, and SNP-CpG site interaction pairs. CONCLUSIONS: In the simulation (200 replications), power appeared generally good for large interaction effects, but smaller effects will require larger studies or consortium collaboration for realizing a sufficient power.

A Novel Gene ARHGAP44 for Longitudinal Changes in Glycated Hemoglobin (HbA1c) in Subjects without Type 2 Diabetes: Evidence from the Long Life Family Study (LLFS) and the Framingham Offspring Study (FOS).

Glycated hemoglobin (HbA1c) indicates average glucose levels over three months and is associated with insulin resistance and type 2 diabetes (T2D). Longitudinal changes in HbA1c (ΔHbA1c) are also associated with aging processes, cognitive performance, and mortality. We analyzed ΔHbA1c in 1,886 non-diabetic Europeans from the Long Life Family Study to uncover gene variants influencing ΔHbA1c. Using growth curve modeling adjusted for multiple covariates, we derived ΔHbA1c and conducted linkage-guided sequence analysis. Our genome-wide linkage scan identified a significant locus on 17p12. In-depth analysis of this locus revealed a variant rs56340929 (explaining 27% of the linkage peak) in the ARHGAP44 gene that was significantly associated with ΔHbA1c. RNA transcription of ARHGAP44 was associated with ΔHbA1c. The Framingham Offspring Study data further supported these findings on the gene level. Together, we found a novel gene ARHGAP44 for ΔHbA1c in family members without T2D. Follow-up studies using longitudinal omics data in large independent cohorts are warranted.

Simulation of a medication and methylation effects on triglycerides in the Genetic Analysis Workshop 20.

The GAW20 simulation data set is based upon the companion Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study fenofibrate clinical trial data set that forms the real data example for GAW20. The simulated data problem consists of 200 simulated replications of what might happen if we were to repeat the GOLDN clinical trial 200 independent times, for these exact same subjects, but using a new fictitious drug (called "genomethate") that has a pharmaco-epigenetic effect on triglyceride response. For each replication, the pre-genomethate values at visits 1 and 2 are constant (ie, pedigree structures, age, sex, all phenotypes, covariates, genome-wide association study (GWAS) genotypes, and visit 2 methylation values), the same as the real GOLDN data across all 200 replications. Only the post-genomethate treatment data (ie, methylation and triglyceride levels for visits 3 and 4) change across the 200 replications. We postulate a growth curve pharmaco-epigenetic response model, in which each patient's response to genomethate treatment is individualized, and is dependent upon their genotype as well as the methylation state for key genes.

Ascertainment correction for Markov chain Monte Carlo segregation and linkage analysis of a quantitative trait.

Although extended pedigrees are often sampled through probands with extreme levels of a quantitative trait, Markov chain Monte Carlo (MCMC) methods for segregation and linkage analysis have not been able to perform ascertainment corrections. Further, the extent to which ascertainment of pedigrees leads to biases in the estimation of segregation and linkage parameters has not been previously studied for MCMC procedures. In this paper, we studied these issues with a Bayesian MCMC approach for joint segregation and linkage analysis, as implemented in the package Loki. We first simulated pedigrees ascertained through individuals with extreme values of a quantitative trait in spirit of the sequential sampling theory of Cannings and Thompson [Cannings and Thompson [1977] Clin. Genet. 12:208-212]. Using our simulated data, we detected no bias in estimates of the trait locus location. However, in addition to allele frequencies, when the ascertainment threshold was higher than or close to the true value of the highest genotypic mean, bias was also found in the estimation of this parameter. When there were multiple trait loci, this bias destroyed the additivity of the effects of the trait loci, and caused biases in the estimation all genotypic means when a purely additive model was used for analyzing the data. To account for pedigree ascertainment with sequential sampling, we developed a Bayesian ascertainment approach and implemented Metropolis-Hastings updates in the MCMC samplers used in Loki. Ascertainment correction greatly reduced biases in parameter estimates. Our method is designed for multiple, but a fixed number of trait loci.