The challenges of genome-wide interaction studies: Lessons to learn from the analysis of HDL blood levels

Genome-wide association studies (GWAS) have revealed 74 single nucleotide polymorphisms (SNPs) associated with high-density lipoprotein cholesterol (HDL) blood levels. This study is, to our knowledge, the first genome-wide interaction study (GWIS) to identify SNP6SNP interactions associated with HDL levels. We performed a GWIS in the Rotterdam Study (RS) cohort I (RS-I) using the GLIDE tool which leverages the massively parallel computing power of Graphics Processing Units (GPUs) to perform linear regression on all genome-wide pairs of SNPs. By performing a meta-analysis together with Rotterdam Study cohorts II and III (RS-II and RS-III), we were able to filter 181 interaction terms with a p-value, 1 · 1028 that replicated in the two independent cohorts. We were not able to replicate any of these interaction term in the AGES, ARIC, CHS, ERF, FHS and NFBC-66 cohorts (Ntotal = 30, 011) when adjusting for multiple testing. Our GWIS resulted in the consistent finding of a possible interaction between rs774801 in ARMC8 (ENSG00000114098) and rs12442098 in SPATA8 (ENSG00000185594) being associated with HDL levels. However, p-values do not reach the preset Bonferroni correction of the p-values. Our study suggest that even for highly genetically determined traits such as HDL the sample sizes needed to detect SNP6SNP interactions are large and the 2-step filtering approaches do not yield a solution. Here we present our analysis plan and our reservations concerning GWIS

Best Practices and Joint Calling of the HumanExome BeadChip: The CHARGE Consortium

Genotyping arrays are a cost effective approach when typing previously-identified genetic polymorphisms in large numbers of samples. One limitation of genotyping arrays with rare variants (e.g., minor allele frequency [MAF] <0.01) is the difficulty that automated clustering algorithms have to accurately detect and assign genotype calls. Combining intensity data from large numbers of samples may increase the ability to accurately call the genotypes of rare variants. Approximately 62,000 ethnically diverse samples from eleve

Gene-gene Interaction Analyses for Atrial Fibrillation

Atrial fibrillation (AF) is a heritable disease that affects more than thirty million individuals worldwide. Extensive efforts have been devoted to the study of genetic determinants of AF. The objective of our study is to examine the effect of gene-gene interaction on AF susceptibility. We performed a large-scale association analysis of gene-gene interactions with AF in 8,173 AF cases, and 65,237 AF-free referents collected from 15 studies for discovery. We examined putative interactions between genome-wide SNPs and 17 known AF-related SNPs. The top interactions were then tested for association in a

GWAS analysis of handgrip and lower body strength in older adults in the CHARGE consortium

Decline in muscle strength with aging is an important predictor of health trajectory in the elderly. Several factors, including genetics, are proposed contributors to variability in muscle strength. To identify genetic contributors to muscle strength, a meta-analysis of genomewide association studies of handgrip was conducted. Grip strength was measured using a handheld dynamometer in 27 581 individuals of European descent over 65 years of age from 14 cohort studies. Genomewide association analysis was conducted on ~2.7 million imputed and genotyped variants (SNPs). Replication of the most significant findings was conducted using data from 6393 individuals from three cohorts. GWAS of lower body strength was also characterized in a subset of cohorts. Two genomewide significant (P-value< 5 × 10−8) and 39 suggestive (P-value< 5 × 10−5) associations were observed from meta-analysis of the discovery cohorts. After meta-analysis with replication cohorts, genomewide significant association was observed for rs752045 on chromosome 8 (β = 0.47, SE = 0.08, P-value = 5.20 × 10−10). This SNP is mapped to an intergenic region and is located within an accessible chromatin region (DNase hypersensitivity site) in skeletal muscle myotubes differentiated from the human skeletal muscle myoblasts cell line. This locus alters a binding motif of the CCAAT/enhancer-binding protein-β (CEBPB) that is implicated in muscle repair mechanisms. GWAS of lower body strength did not yield significant results. A common genetic variant in a chromosomal region that regulates myotube differentiation and muscle repair may contribute to variability in grip strength in the elderly. Further studies are needed to uncover the mechanisms that link this genetic variant with muscle strength

Genetic Interactions with Age, Sex, Body Mass Index, and Hypertension in Relation to Atrial Fibrillation: The AFGen Consortium

It is unclear whether genetic markers interact with risk factors to influence atrial fibrillation (AF) risk. We performed genome-wide interaction analyses between genetic variants and age, sex, hypertension, and body mass index in the AFGen Consortium. Study-specific results were combined using meta-analysis (88,383 individuals of European descent, including 7,292 with AF). Variants with nominal interaction associations in the discovery analysis were tested for association in four independent studies (131,441 individuals, including 5,722 with AF). In the discovery analysis, the AF risk associated with the minor rs6817105 allele (at the PITX2 locus) was greater among subjects ≤ 65 years of age than among those > 65 years (interaction p-value = 4.0 × 10-5). The interaction p-value exceeded genome-wide significance in combined discovery and replication analyses (interaction p-value = 1.7 × 10-8). We observed one genome-wide significant interaction with body mass index and several suggestive interactions with age, sex, and body mass index in the discovery analysis. However, none was replicated in the independent sample. Our findings suggest that the pathogenesis of AF may differ according to age in individuals of European descent, but we did not observe evidence of statistically significant genetic interactions with sex, body mass index, or hypertension on AF risk

Genome-Wide Association Study for Incident Myocardial Infarction and Coronary Heart Disease in Prospective Cohort Studies: The CHARGE Consortium

Background Data are limited on genome-wide association studies (GWAS) for incident coronary heart disease (CHD). Moreover, it is not known whether genetic variants identified to date also associate with risk of CHD in a prospective setting. Methods We performed a two-stageGWAS analysis of incident myocardial infarction (MI) and CHD in a total of 64,297 individuals (including 3898MI cases, 5465 CHD cases). SNPs that passed an arbitrary threshold of 5×10-6 in Stage I were taken to Stage II for further discovery. Furthermore, in an analysis of prognosis, we studied whether known SNPs from former GWAS were associated with totalmortality in individuals who experienced MI during follow-up. Results In Stage I 15 loci passed the threshold of 5×10-6; 8 loci for MI and 8 loci for CHD, for which one locus overlapped and none were reported in previous GWAS meta-analyses. We took 60 SNPs representing these 15 loci to Stage II of discovery. Four SNPs near QKI showed nominally significant association with MI (p-value<8.8×10-3) and three exceeded the genome-wide significance threshold when Stage I and Stage II results were combined (top SNP rs6941513: p = 6.2×10-9). Despite excellent power, the 9p21 locus SNP (rs1333049) was only modestly associated with MI (HR = 1.09, p-value = 0.02) and marginally with CHD (HR = 1.06, p-value = 0.08). Among an inception cohort of those who experienced MI during follow-up, the risk allele of rs1333049 was associated with a decreased risk of subsequent mortality (HR = 0.90, p-value = 3.2×10-3). Conclusions QKI represents a novel locus that may serve as a predictor of incident CHD in prospective studies. The association of the 9p21 locus both with increased risk of first myocardial infarction and longer survival after MI highlights the importance of study design in investigating genetic determinants of complex disorders

RANTES/CCL5 and risk for coronary events: Results from the MONICA/KORA Augsburg case-cohort, Athero-express and CARDIoGRAM studies

Background: The chemokine RANTES (regulated on activation, normal T-cell expressed and secreted)/CCL5 is involved in the pathogenesis of cardiovascular disease in mice, whereas less is known in humans. We hypothesised that its relevance for atherosclerosis should be reflected by associations between CCL5 gene variants, RANTES serum concentrations and protein levels in atherosclerotic plaques and risk for coronary events. Methods and Findings: We conducted a case-cohort study within the population-based MONICA/KORA Augsburg studies. Baseline RANTES serum levels were measured in 363 individuals with incident coronary events and 1,908 non-cases (mean follow-up: 10.2±

Meta-analysis of 49 549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels

Background So far, more than 170 loci have been associated with circulating lipid levels through genomewide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels. Methods We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ~60 000 individuals in the discovery stage and ~90 000 samples in the replication stage. Results Our study resu