Knowledge-driven analysis identifies a gene-gene interaction affecting high-density lipoprotein cholesterol levels in multi-ethnic populations

Identification of epistatic SNP combinations in rheumatoid arthritis using LAMPLINK and Japanese cohorts

Genome-wide association studies have enabled the identification of important genetic factors in many trait studies. However, only a fraction of the heritability can be explained by known genetic factors, even in the most common diseases. Genetic loci combinations, or epistatic contributions expressed by combinations of single nucleotide polymorphisms (SNPs), have been argued to be one of the critical factors explaining some of the missing heritability, especially in oligogenic/polygenic diseases. Rheumatoid arthritis (RA) is a complex disease with more than 100 reported SNP associations, as well as various HLA haplotypes and amino acids; however, many associations between RA and inter-chromosomal SNP combinations are unknown. To discover novel associations of epistatic interactions with high odds ratios in RA, we applied the LAMPLINK method, a systematic enumerative procedure for identifying high-order SNP combinations, to a Japanese RA cohort (discovery cohort; 4024 patients with RA and 7731 controls). We validated the identified associations in a different Japanese cohort (validation cohort; 810 RA patients and 6303 controls). In this study, we identified 90 significant genetic associations in the discovery cohort. Among these, 74 (82.2%) associations were replicated in the validation cohort, and eight combinations were inter-chromosomal, all of which comprised rs7765379 or rs35265698 located in the HLA region. These two SNPs exhibited strong correlations with valine at amino acid position 11 in HLA-DRB1 (HLA-DRB1-11-Val). Finally, we discovered that rs9624 showed an association with RA through an epistatic interaction with HLA-DRB1-11-Val. Overall, LAMPLINK showed high reliability for identifying epistatic genetic contributions hidden in complex traits.

BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

A biobank-scale test of marginal epistasis reveals genome-wide signals of polygenic epistasis

The contribution of epistasis (interactions among genes or genetic variants) to human complex trait variation remains poorly understood. Methods that aim to explicitly identify pairs of genetic variants, usually single nucleotide polymorphisms (SNPs), associated with a trait suffer from low power due to the large number of hypotheses tested while also having to deal with the computational problem of searching over a potentially large number of candidate pairs. An alternate approach involves testing whether a single SNP modulates variation in a trait against a polygenic background. While overcoming the limitation of low power, such tests of polygenic or marginal epistasis (ME) are infeasible on Biobank-scale data where hundreds of thousands of individuals are genotyped over millions of SNPs. We present a method to test for ME of a SNP on a trait that is applicable to biobank-scale data. We performed extensive simulations to show that our method provides calibrated tests of ME. We applied our method to test for ME at SNPs that are associated with 53 quantitative traits across ≈ 300 K unrelated white British individuals in the UK Biobank (UKBB). Testing 15, 601 trait-loci associations that were significant in GWAS, we identified 16 trait-loci pairs across 12 traits that demonstrate strong evidence of ME signals (p-value p < 5 × 10 - 8 53 ). We further partitioned the significant ME signals across the genome to identify 6 trait-loci pairs with evidence of local (within-chromosome) ME while 15 show evidence of distal (cross-chromosome) ME. Across the 16 trait-loci pairs, we document that the proportion of trait variance explained by ME is about 12x as large as that explained by the GWAS effects on average (range: 0.59 to 43.89). Our results show, for the first time, evidence of interaction effects between individual genetic variants and overall polygenic background modulating complex trait variation.

Genetic Dissection of Epistatic Interactions Contributing Yield-Related Agronomic Traits in Rice Using the Compressed Mixed Model

Rice (Oryza sativa) is one of the most important cereal crops in the world, and yield-related agronomic traits, including plant height (PH), panicle length (PL), and protein content (PC), are prerequisites for attaining the desired yield and quality in breeding programs. Meanwhile, the main effects and epistatic effects of quantitative trait nucleotides (QTNs) are all important genetic components for yield-related quantitative traits. In this study, we conducted genome-wide association studies (GWAS) for 413 rice germplasm resources, with 36,901 single nucleotide polymorphisms (SNPs), to identify QTNs, QTN-by-QTN interaction (QQI), and their candidate genes, using a multi-locus compressed variance component mixed model, 3VmrMLM. As a result, two significant QTNs and 56 paired QQIs were detected, amongst 5219 genes of these QTNs, and 26 genes were identified as the yield-related confirmed genes, such as LCRN1, OsSPL3, and OsVOZ1 for PH, and LOG and QsBZR1 for PL. To reveal the substantial contributions related to the variation of yield-related agronomic traits in rice, we further implemented an enrichment analysis and expression analysis. As the results showed, 114 genes, nearly all significant QQIs, were involved in 37 GO terms; for example, the macromolecule metabolic process (GO:0043170), intracellular part (GO:0044424), and binding (GO:0005488). It was revealed that most of the QQIs and the candidate genes were significantly involved in the biological process, molecular function, and cellular component of the target traits. The demonstrated genetic interactions play a critical role in yield-related agronomic traits of rice, and such epistatic interactions contributed to large portions of the missing heritability in GWAS. These results help us to understand the genetic basis underlying the inheritance of the three yield-related agronomic traits and provide implications for rice improvement.

Polygenic scores for dyslipidemia: the emerging genomic model of plasma lipoprotein trait inheritance

PURPOSE OF REVIEW

Contemporary polygenic scores, which summarize the cumulative contribution of millions of common single-nucleotide variants to a phenotypic trait, can have effects comparable to monogenic mutations. This review focuses on the emerging use of 'genome-wide' polygenic scores for plasma lipoproteins to define the etiology of clinical dyslipidemia, modify the severity of monogenic disease, and inform therapeutic options.

RECENT FINDINGS

Polygenic scores for low-density lipoprotein cholesterol (LDL-C), triglycerides, and high-density lipoprotein cholesterol are associated with severe hypercholesterolemia, hypertriglyceridemia, or hypoalphalipoproteinemia, respectively. These polygenic scores for LDL-C or triglycerides associate with risk of incident coronary artery disease (CAD) independent of polygenic scores designed specifically for CAD and may identify individuals that benefit most from lipid-lowering medication. Additionally, the severity of hypercholesterolemia and CAD associated with familial hypercholesterolemia-a common monogenic disorder-is modified by these polygenic factors. The current focus of polygenic scores for dyslipidemia is to design predictive polygenic scores for diverse populations and determining how these polygenic scores could be implemented and standardized for use in the clinic.

SUMMARY

Polygenic scores have shown early promise for the management of dyslipidemias, but several challenges need to be addressed before widespread clinical implementation to ensure that potential benefits are robust and reproducible, equitable, and cost-effective.

Functional Haplotype of LIPC Induces Triglyceride-Mediated Suppression of HDL-C Levels According to Genome-Wide Association Studies

Hepatic lipase (encoded by LIPC) is a glycoprotein in the triacylglycerol lipase family and mainly synthesized in and secreted from the liver. Previous studies demonstrated that hepatic lipase is crucial for reverse cholesterol transport and modulating metabolism and the plasma levels of several lipoproteins. This study was conducted to investigate the suppression effect of high-density lipoprotein cholesterol (HDL-C) levels in a genome-wide association study and explore the possible mechanisms linking triglyceride (TG) to LIPC variants and HDL-C. Genome-wide association data for TG and HDL-C were available for 4657 Taiwan-biobank participants. The prevalence of haplotypes in the LIPC promoter region and their effects were calculated. The cloned constructs of the haplotypes were expressed transiently in HepG2 cells and evaluated in a luciferase reporter assay. Genome-wide association analysis revealed that HDL-C was significantly associated with variations in LIPC after adjusting for TG. Three haplotypes (H1: TCG, H2: CTA and H3: CCA) in LIPC were identified. H2: CTA was significantly associated with HDL-C levels and H1: TCG suppressed HDL-C levels when a third factor, TG, was included in mediation analysis. The luciferase reporter assay further showed that the H2: CTA haplotype significantly inhibited luciferase activity compared with the H1: TCG haplotype. In conclusion, we identified a suppressive role for TG in the genome-wide association between LIPC and HDL-C. A functional haplotype of hepatic lipase may reduce HDL-C levels and is suppressed by TG.

Genetics of Familial Hypercholesterolemia: New Insights

Familial hypercholesterolemia (FH) is one of the most common monogenic diseases, leading to an increased risk of premature atherosclerosis and its cardiovascular complications due to its effect on plasma cholesterol levels. Variants of three genes (LDL-R, APOB and PCSK9) are the major causes of FH, but in some probands, the FH phenotype is associated with variants of other genes. Alternatively, the typical clinical picture of FH can result from the accumulation of common cholesterol-increasing alleles (polygenic FH). Although the Czech Republic is one of the most successful countries with respect to FH detection, approximately 80% of FH patients remain undiagnosed. The opportunities for international collaboration and experience sharing within international programs (e.g., EAS FHSC, ScreenPro FH, etc.) will improve the detection of FH patients in the future and enable even more accessible and accurate genetic diagnostics.

Discovering genetic interactions bridging pathways in genome-wide association studies

Genetic interactions have been reported to underlie phenotypes in a variety of systems, but the extent to which they contribute to complex disease in humans remains unclear. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions, but existing methods for identifying them from GWAS data tend to focus on testing individual locus pairs, which undermines statistical power. Importantly, a global genetic network mapped for a model eukaryotic organism revealed that genetic interactions often connect genes between compensatory functional modules in a highly coherent manner. Taking advantage of this expected structure, we developed a computational approach called BridGE that identifies pathways connected by genetic interactions from GWAS data. Applying BridGE broadly, we discover significant interactions in Parkinson's disease, schizophrenia, hypertension, prostate cancer, breast cancer, and type 2 diabetes. Our novel approach provides a general framework for mapping complex genetic networks underlying human disease from genome-wide genotype data.

Pleiotropic association of LIPC variants with lipid and urinary 8-hydroxy deoxyguanosine levels in a Taiwanese population

BACKGROUND

Hepatic lipase (HL, encoded by LIPC) is a glycoprotein primarily synthesized and secreted by hepatocytes. Previous studies had demonstrated that HL is crucial for reverse cholesterol transport and affects the metabolism, composition, and level of several lipoproteins. In current study, we investigated the association of LIPC (Lipase C, Hepatic Type) variants with circulating and urinary biomarker levels by using subgroup and mediation analyses.

METHODS

A total of 572 participants from Taiwan were genotyped for three LIPC single nucleotide polymorphisms (SNPs) by using TaqMan assay. Fasting levels of glucose, lipid profile, inflammation markers, urine creatinine and 8-hydroxy deoxyguanosine (8-OHdG) were measured. The chi-square test, 2-sample t test and Analysis of variance (ANOVA) were used to examine differences among variables and genotype frequencies.

RESULTS

SNPs rs2043085 and rs1532085 were significantly associated with urinary 8-OHdG levels, whereas all three SNPs were more significantly associated with Triglycerides (TG) or HDL-cholesterol (HDL-C) levels after additional adjustment for HDL-C or TG levels, respectively. Subgroup analyses revealed that the association of the LIPC SNPs with the levels of serum TG, HDL-C, and urinary 8-OHdG were predominantly observed in the men but not in the women. Differential associations of the LIPC SNPs with various lipid levels were observed in participants with different adiposity statuses. Mediation analyses indicated that TG levels acted as a suppressor masking the association of the LIPC genotypes with HDL-C levels, particularly in the men (Sobel test, all P < 0.01).

CONCLUSION

Our data revealed that interaction and suppression effects mediated the pleiotropic association of the LIPC variants. The effects of the LIPC SNPs depended on sex, adiposity status, and TG levels. Thus, our findings can provide a method for identifying high-risk populations of cardiovascular diseases for clinical diagnosis.

Candidate Variants for Additive and Interactive Effects on Bioenergy Traits in Switchgrass (Panicum virgatum L.) Identified by Genome-Wide Association Analyses

Switchgrass ( L.) is a promising herbaceous energy crop, but further gains in biomass yield and quality must be achieved to enable a viable bioenergy industry. Developing DNA markers can contribute to such progress, but depiction of genetic bases should be reliable, involving simple additive marker effects and also interactions with genetic backgrounds (e.g., ecotypes) or synergies with other markers. We analyzed plant height, C content, N content, and mineral concentration in a diverse panel consisting of 512 genotypes of upland and lowland ecotypes. We performed association analyses based on exome capture sequencing and tested 439,170 markers for marginal effects, 83,290 markers for marker × ecotype interactions, and up to 311,445 marker pairs for pairwise interactions. Analyses of pairwise interactions focused on subsets of marker pairs preselected on the basis of marginal marker effects, gene ontology annotation, and pairwise marker associations. Our tests identified 12 significant effects. Homology and gene expression information corroborated seven effects and indicated plausible causal pathways: flowering time and lignin synthesis for plant height; plant growth and senescence for C content and mineral concentration. Four pairwise interactions were detected, including three interactions preselected on the basis of pairwise marker correlations. Furthermore, a marker × ecotype interaction and a pairwise interaction were confirmed in an independent switchgrass panel. Our analyses identified reliable candidate variants for important bioenergy traits. Moreover, they exemplified the importance of interactive effects for depicting genetic bases and illustrated the usefulness of preselecting marker pairs for identifying pairwise marker interactions in association studies.

Detecting epistasis within chromatin regulatory circuitry reveals CAND2 as a novel susceptibility gene for obesity

BACKGROUND

Genome-wide association studies have identified many susceptibility loci for obesity. However, missing heritability problem is still challenging and ignorance of genetic interactions is believed to be an important cause. Current methods for detecting interactions usually do not consider regulatory elements in non-coding regions. Interaction analyses within chromatin regulatory circuitry may identify new susceptibility loci.

METHODS

We developed a pipeline named interaction analyses within chromatin regulatory circuitry (IACRC), to identify genetic interactions impacting body mass index (BMI). Potential interacting SNP pairs were obtained based on Hi-C datasets, PreSTIGE (Predicting Specific Tissue Interactions of Genes and Enhancers) algorithm, and super enhancer regions. SNP × SNP analyses were next performed in three GWAS datasets, including 2286 unrelated Caucasians from Kansas City, 3062 healthy Caucasians from the Gene Environment Association Studies initiative, and 3164 Hispanic subjects from the Women's Health Initiative.

RESULTS

A total of 16,643,227 SNP × SNP analyses were performed. Meta-analyses showed that two SNP pairs, rs6808450-rs9813534 (combined P = 2.39 × 10^-9) and rs6808450-rs3773306 (combined P = 2.89 × 10^-9) were associated with BMI after multiple testing corrections. Single-SNP analyses did not detect significant association signals for these three SNPs. In obesity relevant cells, rs6808450 is located in intergenic enhancers, while rs9813534 and rs3773306 are located in the region of strong transcription regions of CAND2 and RPL32, respectively. The expression of CAND2 was significantly downregulated after the differentiation of human Simpson-Golabi-Behmel syndrome (SGBS) preadipocyte cells (P = 0.0241). Functional validation in the International Mouse Phenotyping Consortium database showed that CAND2 was associated with increased lean body mass and decreased total body fat amount.

CONCLUSIONS

Detecting epistasis within chromatin regulatory circuitry identified CAND2 as a novel obesity susceptibility gene. We hope IACRC could facilitate the interaction analyses for complex diseases and offer new insights into solving the missing heritability problem.

The search for gene-gene interactions in genome-wide association studies: challenges in abundance of methods, practical considerations, and biological interpretation

One of the primary goals in this era of precision medicine is to understand the biology of human diseases and their treatment, such that each individual patient receives the best possible treatment for their disease based on their genetic and environmental exposures. One way to work towards achieving this goal is to identify the environmental exposures and genetic variants that are relevant to each disease in question, as well as the complex interplay between genes and environment. Genome-wide association studies (GWAS) have allowed for a greater understanding of the genetic component of many complex traits. However, these genetic effects are largely small and thus, our ability to use these GWAS finding for precision medicine is limited. As more and more GWAS have been performed, rather than focusing only on common single nucleotide polymorphisms (SNPs) and additive genetic models, many researchers have begun to explore alternative heritable components of complex traits including rare variants, structural variants, epigenetics, and genetic interactions. While genetic interactions are a plausible reality that could explain some of the heritabliy that has not yet been identified, especially when one considers the identification of genetic interactions in model organisms as well as our understanding of biological complexity, still there are significant challenges and considerations in identifying these genetic interactions. Broadly, these can be summarized in three categories: abundance of methods, practical considerations, and biological interpretation. In this review, we will discuss these important elements in the search for genetic interactions along with some potential solutions. While genetic interactions are theoretically understood to be important for complex human disease, the body of evidence is still building to support this component of the underlying genetic architecture of complex human traits. Our hope is that more sophisticated modeling approaches and more robust computational techniques will enable the community to identify these important genetic interactions and improve our ability to implement precision medicine in the future.

Another Round of "Clue" to Uncover the Mystery of Complex Traits

A plethora of genetic association analyses have identified several genetic risk loci. Technological and statistical advancements have now led to the identification of not only common genetic variants, but also low-frequency variants, structural variants, and environmental factors, as well as multi-omics variations that affect the phenotypic variance of complex traits in a population, thus referred to as complex trait architecture. The concept of heritability, or the proportion of phenotypic variance due to genetic inheritance, has been studied for several decades, but its application is mainly in addressing the narrow sense heritability (or additive genetic component) from Genome-Wide Association Studies (GWAS). In this commentary, we reflect on our perspective on the complexity of understanding heritability for human traits in comparison to model organisms, highlighting another round of clues beyond GWAS and an alternative approach, investigating these clues comprehensively to help in elucidating the genetic architecture of complex traits.

The role of gene polymorphisms in endometriosis

Endometriosis is a benign gynecologic disorder, affecting up to 10% of women, characterized by the presence of functional endometrial tissue at ectopic positions generally within the peritoneum. It is a heritable condition influenced by multiple genetic and environmental factors, with an overall heritability estimated at approximately 50%. In this study, we investigated whether single nucleotide polymorphisms (SNPs) rs7521902, rs10859871 and rs11031006, mapping to WNT4, VEZT and FSHB genetic loci, respectively, are associated with risk for endometriosis in a Greek population. This study included 166 women with histologically confirmed endometriosis diagnosed through surgery and 150 normal controls. Genotyping of the rs7521902, rs10859871 and rs11031006 SNPs was performed with Taqman primer/probe sets. A significant association was detected with the AC genotype of rs7521902 (WNT4) in patients with stage III and IV disease only. Evidence for association with endometriosis was also found for the AC genotype of the rs10859871 of VEZT. Notably, a significant difference in the distribution of the AG genotype and the minor allele A of FSHB rs11031006 SNP was found between the endometriosis patients and controls. We found a genetic association between rs11031006 (FSHB) SNP and endometriosis. WNT4 and VEZT genes constitute the most consistently associated genes with endometriosis. In the present study, an association of rs7521902 (WNT4) and rs10859871 (VEZT) was confirmed in women with endometriosis at the genotypic but not the allelic level.

Genetic associations with lipoprotein subfraction measures differ by ethnicity in the multi-ethnic study of atherosclerosis (MESA)

A recent genome-wide association study associated 62 single nucleotide polymorphisms (SNPs) from 43 genomic loci, with fasting lipoprotein subfractions in European-Americans (EAs) at genome-wide levels of significance across three independent samples. Whether these associations are consistent across ethnicities with a non-European ancestry is unknown. We analyzed 15 lipoprotein subfraction measures, on 1677 African-Americans (AAs), 1450 Hispanic-Americans (HAs), and 775 Chinese-Americans (CHN) participating in the multi-ethnic study of atherosclerosis (MESA). Genome-wide data were obtained using the Affymetrix 6.0 and Illumina HumanOmni chips. Linear regression models between genetic variables and lipoprotein subfractions were adjusted for age, gender, body mass index, smoking, study center, and genetic ancestry (based on principal components), and additionally adjusted for Mexican/Non-Mexican status in HAs. A false discovery rate correction was applied separately within the results for each ethnicity to correct for multiple testing. Power calculations revealed that we did not have the power for SNP-based measures of association, so we analyzed phenotype-specific genetic risk scores (GRSs), constructed as in the original genome-wide analysis. We successfully replicated all 15 GRS-lipoprotein associations in 2527 EAs. Among the 15 significant GRS-lipoprotein associations in EAs, 11 were significant in AAs, 13 in HAs, and 1 in CHNs. Further analyses revealed that ethnicity differences could not be explained by differences in linkage disequilibrium, lipid lowering drugs, diabetes, or gender. Our study emphasizes the importance of ethnicity (here indexing genetic ancestry) in genetic risk for CVD and highlights the need to identify ethnicity-specific genetic variants associated with CVD risk.

Dissection of additive, dominance, and imprinting effects for production and reproduction traits in Holstein cattle

BACKGROUND

Although genome-wide association and genomic selection studies have primarily focused on additive effects, dominance and imprinting effects play an important role in mammalian biology and development. The degree to which these non-additive genetic effects contribute to phenotypic variation and whether QTL acting in a non-additive manner can be detected in genetic association studies remain controversial.

RESULTS

To empirically answer these questions, we analyzed a large cattle dataset that consisted of 42,701 genotyped Holstein cows with genotyped parents and phenotypic records for eight production and reproduction traits. SNP genotypes were phased in pedigree to determine the parent-of-origin of alleles, and a three-component GREML was applied to obtain variance decomposition for additive, dominance, and imprinting effects. The results showed a significant non-zero contribution from dominance to production traits but not to reproduction traits. Imprinting effects significantly contributed to both production and reproduction traits. Interestingly, imprinting effects contributed more to reproduction traits than to production traits. Using GWAS and imputation-based fine-mapping analyses, we identified and validated a dominance association signal with milk yield near RUNX2, a candidate gene that has been associated with milk production in mice. When adding non-additive effects into the prediction models, however, we observed little or no increase in prediction accuracy for the eight traits analyzed.

CONCLUSIONS

Collectively, our results suggested that non-additive effects contributed a non-negligible amount (more for reproduction traits) to the total genetic variance of complex traits in cattle, and detection of QTLs with non-additive effect is possible in GWAS using a large dataset.

Comparing GWAS Results of Complex Traits Using Full Genetic Model and Additive Models for Revealing Genetic Architecture

Most of the genome-wide association studies (GWASs) for human complex diseases have ignored dominance, epistasis and ethnic interactions. We conducted comparative GWASs for total cholesterol using full model and additive models, which illustrate the impacts of the ignoring genetic variants on analysis results and demonstrate how genetic effects of multiple loci could differ across different ethnic groups. There were 15 quantitative trait loci with 13 individual loci and 3 pairs of epistasis loci identified by full model, whereas only 14 loci (9 common loci and 5 different loci) identified by multi-loci additive model. Again, 4 full model detected loci were not detected using multi-loci additive model. PLINK-analysis identified two loci and GCTA-analysis detected only one locus with genome-wide significance. Full model identified three previously reported genes as well as several new genes. Bioinformatics analysis showed some new genes are related with cholesterol related chemicals and/or diseases. Analyses of cholesterol data and simulation studies revealed that the full model performs were better than the additive-model performs in terms of detecting power and unbiased estimations of genetic variants of complex traits.

Role of rs3846662 and HMGCR alternative splicing in statin efficacy and baseline lipid levels in familial hypercholesterolemia

OBJECTIVES

To assess the contribution of the rs3846662 polymorphism of HMGCR on serum lipid levels and statin efficacy, we measured in vivo HMGCR mRNA and lipid levels in French Canadian individuals affected by heterozygous familial hypercholesterolemia due to the deletion of more than 15 kb of the LDLR gene.

RESULTS

Men and women carrying the AA genotype at rs3846662, and no APOE4 allele, had higher levels of total cholesterol (5.43 vs. 4.58 mmol/l, P<0.05) and LDL-cholesterol (5.20 vs. 4.39 mmol/l, P<0.05) at baseline. However, with regard to statin efficacy, the penetrance of the AA genotype was sex dependent. Indeed, the percentage reduction in LDL-cholesterol upon statin treatment was significantly decreased in women with the AA genotype compared with women without it (38.4 vs. 46.2%, P<0.05), whereas this was not observed in men. Although both men and women bearing the AA genotype showed a higher ratio of full-length HMGCR mRNA/total HMGCR mRNA compared with individuals without it (n=37, P<0.05), overall transcription of HMGCR was decreased and increased in men and women carrying this genotype, respectively (n=37, P<0.01 and P<0.05). Finally, in our familial hypercholesterolemia cohort, HMGCR alternative splicing explained between 22 and 55% of the variance in statin response.

CONCLUSION

rs3846662 polymorphism and the alternative splicing of HMGCR mRNA significantly impact women's response to statin therapy.

Polymorphisms involving gain or loss of CpG sites are significantly enriched in trait-associated SNPs

Some single nucleotide polymorphisms (SNPs) influence the existence of CpG sites, the basis of DNA modification such as methylation and hydroxymethylation. These polymorphisms can lead to gain or loss of CpG sites and were defined as CpG site related SNPs (cgSNPs) in this study. The cgSNPs change DNA sequence and might potentially affect DNA modification such as methylation. However, the functional consequence of cgSNPs is poorly understood. We observed that a considerable proportion (23.0%) of common variants were cgSNPs in human genome. Mutations involving loss of CpG sites were associated with reduced levels of methylation (~20.2%) using The Cancer Genome Atlas (TCGA) data. Using public databases (SCAN and seeQTL) of expression quantitative trait loci (eQTLs), we found that the cgSNPs were significantly enriched in eQTLs via logistic regression and simulation test. Furthermore, we observed that cgSNPs were more likely to be trait-associated loci especially cancers using a catalog of published genome-wide association studies (GWAS) recorded by National Human Genome Research Institute (NHGRI). Our results indicated that cgSNP might be meaningful as annotation either in SNP functional prediction or in screening for trait-associated SNPs.

Pathway Analysis Incorporating Protein-Protein Interaction Networks Identified Candidate Pathways for the Seven Common Diseases

Pathway analysis has become popular as a secondary analysis strategy for genome-wide association studies (GWAS). Most of the current pathway analysis methods aggregate signals from the main effects of single nucleotide polymorphisms (SNPs) in genes within a pathway without considering the effects of gene-gene interactions. However, gene-gene interactions can also have critical effects on complex diseases. Protein-protein interaction (PPI) networks have been used to define gene pairs for the gene-gene interaction tests. Incorporating the PPI information to define gene pairs for interaction tests within pathways can increase the power for pathway-based association tests. We propose a pathway association test, which aggregates the interaction signals in PPI networks within a pathway, for GWAS with case-control samples. Gene size is properly considered in the test so that genes do not contribute more to the test statistic simply due to their size. Simulation studies were performed to verify that the method is a valid test and can have more power than other pathway association tests in the presence of gene-gene interactions within a pathway under different scenarios. We applied the test to the Wellcome Trust Case Control Consortium GWAS datasets for seven common diseases. The most significant pathway is the chaperones modulate interferon signaling pathway for Crohn’s disease (p-value = 0.0003). The pathway modulates interferon gamma, which induces the JAK/STAT pathway that is involved in Crohn’s disease. Several other pathways that have functional implications for the seven diseases were also identified. The proposed test based on gene-gene interaction signals in PPI networks can be used as a complementary tool to the current existing pathway analysis methods focusing on main effects of genes. An efficient software implementing the method is freely available at http://puppi.sourceforge.net.

Epistatic Gene-Based Interaction Analyses for Glaucoma in eMERGE and NEIGHBOR Consortium

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.

The complex nature of primary-open angle glaucoma (POAG) has left researchers exploring the genetic architecture and searching for the missing heritability using a number of different study designs. Over the past decade, many studies have been conducted to explain the etiology of POAG; however, a high proportion of estimated heritability still remains unexplained. GWA studies for POAG have identified significant associations but these associations have only explained a small proportion of the genetic risk (odds ratios range between 1–3). In this paper, we sought to confirm the primary genome-wide significant associations that have been discovered so far for glaucoma in phenotypes developed from EMR data in an effort to show that EMR data can be a powerful resource for finding genetic variants influencing POAG susceptibility. Next, we tested for statistical interactions, which can be presented as an important tool in an attempt to explain POAG heritability. We used a reduced list of variants filtered by marginal main effect analysis to look for epistatic interactions. We present our results from replication of gene-based interaction analyses performed in eMERGE and the NEIGHBOR consortium data. Using expression data and annotations from various publicly available databases, the most significant genes that replicated in our analyses show expression in the eye and trabecular meshwork. Analysis for estimation of genetic variance explained by significant associations from previous GWAS and replicated variants from gene-based interactions suggest that these explain 5.6% of variance in eMERGE dataset and also explain 3.4% variance in NEIGHBOR dataset.

A Novel Test for Detecting SNP-SNP Interactions in Case-Only Trio Studies

Epistasis plays a significant role in the genetic architecture of many complex phenotypes in model organisms. To date, there have been very few interactions replicated in human studies due in part to the multiple-hypothesis burden implicit in genome-wide tests of epistasis. Therefore, it is of paramount importance to develop the most powerful tests possible for detecting interactions. In this work we develop a new SNP-SNP interaction test for use in case-only trio studies called the trio correlation (TC) test. The TC test computes the expected joint distribution of marker pairs in offspring conditional on parental genotypes. This distribution is then incorporated into a standard 1 d.f. correlation test of interaction. We show via extensive simulations under a variety of disease models that our test substantially outperforms existing tests of interaction in case-only trio studies. We also demonstrate a bias in a previous case-only trio interaction test and identify its origin. Finally, we show that a previously proposed permutation scheme in trio studies mitigates the known biases of case-only tests in the presence of population stratification. We conclude that the TC test shows improved power to identify interactions in existing, as well as emerging, trio association studies. The method is publicly available at www.github.com/BrunildaBalliu/TrioEpi.

PHENOME-WIDE INTERACTION STUDY (PheWIS) IN AIDS CLINICAL TRIALS GROUP DATA (ACTG)

Association studies have shown and continue to show a substantial amount of success in identifying links between multiple single nucleotide polymorphisms (SNPs) and phenotypes. These studies are also believed to provide insights toward identification of new drug targets and therapies. Albeit of all the success, challenges still remain for applying and prioritizing these associations based on available biological knowledge. Along with single variant association analysis, genetic interactions also play an important role in uncovering the etiology and progression of complex traits. For gene-gene interaction analysis, selection of the variants to test for associations still poses a challenge in identifying epistatic interactions among the large list of variants available in high-throughput, genome-wide datasets. Therefore in this study, we propose a pipeline to identify interactions among genetic variants that are associated with multiple phenotypes by prioritizing previously published results from main effect association analysis (genome-wide and phenome-wide association analysis) based on a-priori biological knowledge in AIDS Clinical Trials Group (ACTG) data. We approached the prioritization and filtration of variants by using the results of a previously published single variant PheWAS and then utilizing biological information from the Roadmap Epigenome project. We removed variants in low functional activity regions based on chromatin states annotation and then conducted an exhaustive pairwise interaction search using linear regression analysis. We performed this analysis in two independent pre-treatment clinical trial datasets from ACTG to allow for both discovery and replication. Using a regression framework, we observed 50,798 associations that replicate at p-value 0.01 for 26 phenotypes, among which 2,176 associations for 212 unique SNPs for fasting blood glucose phenotype reach Bonferroni significance and an additional 9,970 interactions for high-density lipoprotein (HDL) phenotype and fasting blood glucose (total of 12,146 associations) reach FDR significance. We conclude that this method of prioritizing variants to look for epistatic interactions can be used extensively for generating hypotheses for genomewide and phenome-wide interaction analyses. This original Phenome-wide Interaction study (PheWIS) can be applied further to patients enrolled in randomized clinical trials to establish the relationship between patient's response to a particular drug therapy and non-linear combination of variants that might be affecting the outcome.

Revisiting Human Cholesterol Synthesis and Absorption: The Reciprocity Paradigm and its Key Regulators

Hypercholesterolemia is a major risk factor for cardiovascular disease. Cholesterol homeostasis in the body is governed by the interplay between absorption, synthesis, and excretion or conversion of cholesterol into bile acids. A reciprocal relationship between cholesterol synthesis and absorption is known to regulate circulating cholesterol in response to dietary or therapeutic interventions. However, the degree to which these factors affect synthesis and absorption and the extent to which one vector shifts in response to the other are not thoroughly understood. Also, huge inter-individual variability exists in the manner in which the two systems act in response to any cholesterol-lowering treatment. Various factors are known to account for this variability and in light of recent experimental advances new players such as gene-gene interactions, gene-environmental effects, and gut microbiome hold immense potential in offering an explanation to the complex traits of inter-individual variability in human cholesterol metabolism. In this context, the objective of the present review is to provide an overview on cholesterol metabolism and discuss the role of potential factors such as genetics, epigenetics, epistasis, and gut microbiome, as well as other regulators in modulating cholesterol metabolism, especially emphasizing the reciprocal relationship between cholesterol synthesis and absorption. Furthermore, an evaluation of the implications of this push-pull mechanism on cholesterol-lowering strategies is presented.

Revisiting Human Cholesterol Synthesis and Absorption: The Reciprocity Paradigm and its Key Regulators

Hypercholesterolemia is a major risk factor for cardiovascular disease. Cholesterol homeostasis in the body is governed by the interplay between absorption, synthesis, and excretion or conversion of cholesterol into bile acids. A reciprocal relationship between cholesterol synthesis and absorption is known to regulate circulating cholesterol in response to dietary or therapeutic interventions. However, the degree to which these factors affect synthesis and absorption and the extent to which one vector shifts in response to the other are not thoroughly understood. Also, huge inter-individual variability exists in the manner in which the two systems act in response to any cholesterol-lowering treatment. Various factors are known to account for this variability and in light of recent experimental advances new players such as gene-gene interactions, gene-environmental effects, and gut microbiome hold immense potential in offering an explanation to the complex traits of inter-individual variability in human cholesterol metabolism. In this context, the objective of the present review is to provide an overview on cholesterol metabolism and discuss the role of potential factors such as genetics, epigenetics, epistasis, and gut microbiome, as well as other regulators in modulating cholesterol metabolism, especially emphasizing the reciprocal relationship between cholesterol synthesis and absorption. Furthermore, an evaluation of the implications of this push-pull mechanism on cholesterol-lowering strategies is presented.

Genome-Wide Analysis of Gene-Gene and Gene-Environment Interactions Using Closed-Form Wald Tests

Despite the successful discovery of hundreds of variants for complex human traits using genome-wide association studies, the degree to which genes and environmental risk factors jointly affect disease risk is largely unknown. One obstacle toward this goal is that the computational effort required for testing gene-gene and gene-environment interactions is enormous. As a result, numerous computationally efficient tests were recently proposed. However, the validity of these methods often relies on unrealistic assumptions such as additive main effects, main effects at only one variable, no linkage disequilibrium between the two single-nucleotide polymorphisms (SNPs) in a pair or gene-environment independence. Here, we derive closed-form and consistent estimates for interaction parameters and propose to use Wald tests for testing interactions. The Wald tests are asymptotically equivalent to the likelihood ratio tests (LRTs), largely considered to be the gold standard tests but generally too computationally demanding for genome-wide interaction analysis. Simulation studies show that the proposed Wald tests have very similar performances with the LRTs but are much more computationally efficient. Applying the proposed tests to a genome-wide study of multiple sclerosis, we identify interactions within the major histocompatibility complex region. In this application, we find that (1) focusing on pairs where both SNPs are marginally significant leads to more significant interactions when compared to focusing on pairs where at least one SNP is marginally significant; and (2) parsimonious parameterization of interaction effects might decrease, rather than increase, statistical power.

Confronting the missing epistasis problem: on the reproducibility of gene-gene interactions

Epistasis (gene-gene interaction) is thought to play an integral role in the genetic basis of complex traits, and a significant amount of research has been invested into identifying this phenomenon in human disease. However, the overall success of empirical studies of epistasis in humans is unclear, as such studies are rarely systematically evaluated. Here, we have selected asthma as an example of a well-studied, complex human disease, and provide a critical analysis and replication attempt of nearly all prior reports of epistasis for this disease. Of 191 previously reported interactions, we find that 39.8% were not originally identified using an explicit test for interaction and thus may not have been true epistatic effects to begin with. Moreover, directions of effect were not described for 46.1% of the interactions, which prevents their rigorous replication. In the original studies, attempts at replication were made for 15.2% of the interactions, and 7.3% were actually replicated. In the current study, we were able to evaluate 85.9% of the interactions using a large asthma dataset from the GABRIEL Consortium. None of these interactions could be replicated based on strict criteria. However, we found nominally significant (p < 0.05) evidence in support of 23.8% of the evaluated interactions. Although many reports of epistasis are not robustly supported in the published literature, our results suggest that at least some of these reports may have been true-positive examples of epistasis. In general, improvements in empirical studies of epistasis are called for, in order to better understand the importance of this phenomenon in human disease.

Biology-Driven Gene-Gene Interaction Analysis of Age-Related Cataract in the eMERGE Network

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)-rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

Investigating the Role of Gene-Gene Interactions in TB Susceptibility

Tuberculosis (TB) is the second leading cause of mortality from infectious disease worldwide. One of the factors involved in developing disease is the genetics of the host, yet the field of TB susceptibility genetics has not yielded the answers that were expected. A commonly posited explanation for the missing heritability of complex disease is gene-gene interactions, also referred to as epistasis. In this study we investigate the role of gene-gene interactions in genetic susceptibility to TB using a cohort recruited from a high TB incidence community from Cape Town, South Africa. Our discovery data set incorporates genotypes from a large a number of candidate gene studies as well as genome-wide data. After limiting our search space to pairs of putative TB susceptibility genes, as well as pairs of genes that have been curated in online databases as potential interactors, we use statistical modelling to identify pairs of interacting SNPs. We attempt to validate the top models identified in our discovery data set using an independent genome-wide TB case-control data set from The Gambia. A number of models were successfully validated, indicating that interplay between the NRG1 - NRG3, GRIK1 - GRIK3 and IL23R - ATG4C gene pairs may modify susceptibility to TB. Gene pairs involved in the NF-κB pathway were also identified in the discovery data set (SFTPD - NOD2, ISG15 - TLR8 and NLRC5 - IL12RB1), but could not be tested in the Gambian study group due to lack of overlapping data.

Genome-wide epistatic expression quantitative trait loci discovery in four human tissues reveals the importance of local chromosomal interactions governing gene expression

Background

Epistasis (synergistic interaction) among SNPs governing gene expression is likely to arise within transcriptional networks. However, the power to detect it is limited by the large number of combinations to be tested and the modest sample sizes of most datasets. By limiting the interaction search space firstly to cis-trans and then cis-cis SNP pairs where both SNPs had an independent effect on the expression of the most variable transcripts in the liver and brain, we greatly reduced the size of the search space.

Results

Within the cis-trans search space we discovered three transcripts with significant epistasis. Surprisingly, all interacting SNP pairs were located nearby each other on the chromosome (within 290 kb-2.16 Mb). Despite their proximity, the interacting SNPs were outside the range of linkage disequilibrium (LD), which was absent between the pairs (r² < 0.01). Accordingly, we redefined the search space to detect cis-cis interactions, where a cis-SNP was located within 10 Mb of the target transcript. The results of this show evidence for the epistatic regulation of 50 transcripts across the tissues studied. Three transcripts, namely, HLA-G, PSORS1C1 and HLA-DRB5 share common regulatory SNPs in the pre-frontal cortex and their expression is significantly correlated. This pattern of epistasis is consistent with mediation via long-range chromatin structures rather than the binding of transcription factors in trans. Accordingly, some of the interactions map to regions of the genome known to physically interact in lymphoblastoid cell lines while others map to known promoter and enhancer elements. SNPs involved in interactions appear to be enriched for promoter markers.

Conclusions

In the context of gene expression and its regulation, our analysis indicates that the study of cis-cis or local epistatic interactions may have a more important role than interchromosomal interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1300-3) contains supplementary material, which is available to authorized users.

Analysis of gene-gene interactions among common variants in candidate cardiovascular genes in coronary artery disease

Objective

Only a small fraction of coronary artery disease (CAD) heritability has been explained by common variants identified to date. Interactions between genes of importance to cardiovascular regulation may account for some of the missing heritability of CAD. This study aimed to investigate the role of gene-gene interactions in common variants in candidate cardiovascular genes in CAD.

Approach and Results

2,101 patients with CAD from the British Heart Foundation Family Heart Study and 2,426 CAD-free controls were included in the discovery cohort. All subjects were genotyped with the Illumina HumanCVD BeadChip enriched for genes and pathways relevant to the cardiovascular system and disease. The primary analysis in the discovery cohort examined pairwise interactions among 913 common (minor allele frequency >0.1) independent single nucleotide polymorphisms (SNPs) with at least nominal association with CAD in single locus analysis. A secondary exploratory interaction analysis was performed among all 11,332 independent common SNPs surviving quality control criteria. Replication analyses were conducted in 2,967 patients and 3,075 controls from the Myocardial Infarction Genetics Consortium. None of the interactions amongst 913 SNPs analysed in the primary analysis was statistically significant after correction for multiple testing (required P<1.2x10^-7). Similarly, none of the pairwise gene-gene interactions in the secondary analysis reached statistical significance after correction for multiple testing (required P = 7.8x10^-10). None of 36 suggestive interactions from the primary analysis or 31 interactions from the secondary analysis was significant in the replication cohort. Our study had 80% power to detect odds ratios > 1.7 for common variants in the primary analysis.

Conclusions

Moderately large additive interactions between common SNPs in genes relevant to cardiovascular disease do not appear to play a major role in genetic predisposition to CAD. The role of genetic interactions amongst less common SNPs and with medium and small magnitude effects remain to be investigated.

Biological knowledge-driven analysis of epistasis in human GWAS with application to lipid traits

While the importance of epistasis is well established, specific gene-gene interactions have rarely been identified in human genome-wide association studies (GWAS), mainly due to low power associated with such interaction tests. In this chapter, we integrate biological knowledge and human GWAS data to reveal epistatic interactions underlying quantitative lipid traits, which are major risk factors for coronary artery disease. To increase power to detect interactions, we only tested pairs of SNPs filtered by prior biological knowledge, including GWAS results, protein-protein interactions (PPIs), and pathway information. Using published GWAS and 9,713 European Americans (EA) from the Atherosclerosis Risk in Communities (ARIC) study, we identified an interaction between HMGCR and LIPC affecting high-density lipoprotein cholesterol (HDL-C) levels. We then validated this interaction in additional multiethnic cohorts from ARIC, the Framingham Heart Study, and the Multi-Ethnic Study of Atherosclerosis. Both HMGCR and LIPC are involved in the metabolism of lipids and lipoproteins, and LIPC itself has been marginally associated with HDL-C. Furthermore, no significant interaction was detected using PPI and pathway information, mainly due to the stringent significance level required after correcting for the large number of tests conducted. These results suggest the potential of biological knowledge-driven approaches to detect epistatic interactions in human GWAS, which may hold the key to exploring the role gene-gene interactions play in connecting genotypes and complex phenotypes in future GWAS.

CardioGxE, a catalog of gene-environment interactions for cardiometabolic traits

Background

Genetic understanding of complex traits has developed immensely over the past decade but remains hampered by incomplete descriptions of contribution to phenotypic variance. Gene-environment (GxE) interactions are one of these contributors and in the guise of diet and physical activity are important modulators of cardiometabolic phenotypes and ensuing diseases.

Results

We mined the scientific literature to collect GxE interactions from 386 publications for blood lipids, glycemic traits, obesity anthropometrics, vascular measures, inflammation and metabolic syndrome, and introduce CardioGxE, a gene-environment interaction resource. We then analyzed the genes and SNPs supporting cardiometabolic GxEs in order to demonstrate utility of GxE SNPs and to discern characteristics of these important genetic variants. We were able to draw many observations from our extensive analysis of GxEs. 1) The CardioGxE SNPs showed little overlap with variants identified by main effect GWAS, indicating the importance of environmental interactions with genetic factors on cardiometabolic traits. 2) These GxE SNPs were enriched in adaptation to climatic and geographical features, with implications on energy homeostasis and response to physical activity. 3) Comparison to gene networks responding to plasma cholesterol-lowering or regression of atherosclerotic plaques showed that GxE genes have a greater role in those responses, particularly through high-energy diets and fat intake, than do GWAS-identified genes for the same traits. Other aspects of the CardioGxE dataset were explored.

Conclusions

Overall, we demonstrate that SNPs supporting cardiometabolic GxE interactions often exhibit transcriptional effects or are under positive selection. Still, not all such SNPs can be assigned potential functional or regulatory roles often because data are lacking in specific cell types or from treatments that approximate the environmental factor of the GxE. With research on metabolic related complex disease risk embarking on genome-wide GxE interaction tests, CardioGxE will be a useful resource.

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 SNP×SNP 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 · 10⁻⁸ 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 (N_total = 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 SNP×SNP 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.

Detecting epistasis in human complex traits

Genome-wide association studies (GWASs) have become the focus of the statistical analysis of complex traits in humans, successfully shedding light on several aspects of genetic architecture and biological aetiology. Single-nucleotide polymorphisms (SNPs) are usually modelled as having additive, cumulative and independent effects on the phenotype. Although evidently a useful approach, it is often argued that this is not a realistic biological model and that epistasis (that is, the statistical interaction between SNPs) should be included. The purpose of this Review is to summarize recent directions in methodology for detecting epistasis and to discuss evidence of the role of epistasis in human complex trait variation. We also discuss the relevance of epistasis in the context of GWASs and potential hazards in the interpretation of statistical interaction terms.

The Effect of Residing Altitude on Levels of High-Density Lipoprotein Cholesterol: A Pilot Study From the Omani Arab Population

Lower mortality rates from coronary heart disease and higher levels of serum high-density lipoprotein cholesterol (HDL-C) have been observed in populations residing at high altitude. However, this effect has not been investigated in Arab populations, which exhibit considerable genetic homogeneity. We assessed the relationship between residing altitude and HDL-C in 2 genetically similar Omani Arab populations residing at different altitudes. The association between the levels of HDL-C and other metabolic parameters was also investigated. The levels of HDL-C were significantly higher in the high-altitude group compared with the low-altitude group. Stepwise regression analysis showed that altitude was the most significant factor affecting HDL-C, followed by gender, serum triglycerides, and finally the 2-hour postprandial plasma glucose. This finding is consistent with previously published studies from other populations and should be taken into consideration when comparing cardiovascular risk factors in populations residing at different altitudes.

Genetic distance as an alternative to physical distance for definition of gene units in association studies

BACKGROUND

Some association studies, as the implemented in VEGAS, ALIGATOR, i-GSEA4GWAS, GSA-SNP and other software tools, use genes as the unit of analysis. These genes include the coding sequence plus flanking sequences. Polymorphisms in the flanking sequences are of interest because they involve cis-regulatory elements or they inform on untyped genetic variants trough linkage disequilibrium. Gene extensions have customarily been defined as ±50 Kb. This approach is not fully satisfactory because genetic relationships between neighbouring sequences are a function of genetic distances, which are only poorly replaced by physical distances.

RESULTS

Standardized recombination rates (SRR) from the deCODE recombination map were used as units of genetic distances. We searched for a SRR producing flanking sequences near the ±50 Kb offset that has been common in previous studies. A SRR≥2 was selected because it led to gene extensions with median length=45.3 Kb and the simplicity of an integer value. As expected, boundaries of the genes defined with the ±50 Kb and with the SRR≥2 rules were rarely concordant. The impact of these differences was illustrated with the interpretation of top association signals from two large studies including many hits and their detailed analysis based in different criteria. The definition based in genetic distance was more concordant with the results of these studies than the based in physical distance. In the analysis of 18 top disease associated loci form the first study, the SRR≥2 genes led to a fully concordant interpretation in 17 loci; the ±50 Kb genes only in 6. Interpretation of the 43 putative functional genes of the second study based in the SRR≥2 definition only missed 4 of the genes, whereas the based in the ±50 Kb definition missed 10 genes.

CONCLUSIONS

A gene definition based on genetic distance led to results more concordant with expert detailed analyses than the commonly used based in physical distance. The genome coordinates for each gene are provided to maintain a simple use of the new definitions.

Abundant local interactions in the 4p16.1 region suggest functional mechanisms underlying SLC2A9 associations with human serum uric acid

Human serum uric acid concentration (SUA) is a complex trait. A recent meta-analysis of multiple genome-wide association studies (GWAS) identified 28 loci associated with SUA jointly explaining only 7.7% of the SUA variance, with 3.4% explained by two major loci (SLC2A9 and ABCG2). Here we examined whether gene-gene interactions had any roles in regulating SUA using two large GWAS cohorts included in the meta-analysis [the Atherosclerosis Risk in Communities study cohort (ARIC) and the Framingham Heart Study cohort (FHS)]. We found abundant genome-wide significant local interactions in ARIC in the 4p16.1 region located mostly in an intergenic area near SLC2A9 that were not driven by linkage disequilibrium and were replicated in FHS. Taking the forward selection approach, we constructed a model of five SNPs with marginal effects and three epistatic SNP pairs in ARIC-three marginal SNPs were located within SLC2A9 and the remaining SNPs were all located in the nearby intergenic area. The full model explained 1.5% more SUA variance than that explained by the lead SNP alone, but only 0.3% was contributed by the marginal and epistatic effects of the SNPs in the intergenic area. Functional analysis revealed strong evidence that the epistatically interacting SNPs in the intergenic area were unusually enriched at enhancers active in ENCODE hepatic (HepG2, P = 4.7E-05) and precursor red blood (K562, P = 5.0E-06) cells, putatively regulating transcription of WDR1 and SLC2A9. These results suggest that exploring epistatic interactions is valuable in uncovering the complex functional mechanisms underlying the 4p16.1 region.

Analysis pipeline for the epistasis search - statistical versus biological filtering

Gene-gene interactions may contribute to the genetic variation underlying complex traits but have not always been taken fully into account. Statistical analyses that consider gene-gene interaction may increase the power of detecting associations, especially for low-marginal-effect markers, and may explain in part the "missing heritability." Detecting pair-wise and higher-order interactions genome-wide requires enormous computational power. Filtering pipelines increase the computational speed by limiting the number of tests performed. We summarize existing filtering approaches to detect epistasis, after distinguishing the purposes that lead us to search for epistasis. Statistical filtering includes quality control on the basis of single marker statistics to avoid the analysis of bad and least informative data, and limits the search space for finding interactions. Biological filtering includes targeting specific pathways, integrating various databases based on known biological and metabolic pathways, gene function ontology and protein-protein interactions. It is increasingly possible to target single-nucleotide polymorphisms that have defined functions on gene expression, though not belonging to protein-coding genes. Filtering can improve the power of an interaction association study, but also increases the chance of missing important findings.

Human dopamine transporter gene: differential regulation of 18-kb haplotypes

AIM

Since previous functional studies of short haplotypes and polymorphic sites of SLC6A3 have shown variant-dependent and drug-sensitive promoter activity, this study aimed to understand whether a large SLC6A3 regulatory region, containing these small haplotypes and polymorphic sites, can display haplotype-dependent promoter activity in a drug-sensitive and pathway-related manner.

MATERIALS & METHODS

By creating and using a single copy number luciferase-reporter vector, we examined regulation of two different SLC6A3 haplotypes (A and B) of the 5´ 18-kb promoter and two known downstream regulatory variable number tandem repeats by 17 drugs in four different cellular models.

RESULTS

The two regulatory haplotypes displayed up to 3.2-fold difference in promoter activity. The regulations were drug selective (37.5% of the drugs showed effects), and both haplotype and cell type dependent. Pathway analysis revealed at least 13 main signaling hubs targeting SLC6A3, including histone deacetylation, AKT, PKC and CK2 α-chains.

CONCLUSION

SLC6A3 may be regulated via either its promoter or the variable number tandem repeats independently by specific signaling pathways and in a haplotype-dependent manner. Furthermore, we have developed the first pathway map for SLC6A3 regulation. These findings provide a framework for understanding complex and variant-dependent regulations of SLC6A3.

Analysis of multiple association studies provides evidence of an expression QTL hub in gene-gene interaction network affecting HDL cholesterol levels

Epistasis has been suggested to underlie part of the missing heritability in genome-wide association studies. In this study, we first report an analysis of gene-gene interactions affecting HDL cholesterol (HDL-C) levels in a candidate gene study of 2,091 individuals with mixed dyslipidemia from a clinical trial. Two additional studies, the Atherosclerosis Risk in Communities study (ARIC; n = 9,713) and the Multi-Ethnic Study of Atherosclerosis (MESA; n = 2,685), were considered for replication. We identified a gene-gene interaction between rs1532085 and rs12980554 (P = 7.1×10⁻⁷) in their effect on HDL-C levels, which is significant after Bonferroni correction (P_c = 0.017) for the number of SNP pairs tested. The interaction successfully replicated in the ARIC study (P = 7.0×10⁻⁴; P_c = 0.02). Rs1532085, an expression QTL (eQTL) of LIPC, is one of the two SNPs involved in another, well-replicated gene-gene interaction underlying HDL-C levels. To further investigate the role of this eQTL SNP in gene-gene interactions affecting HDL-C, we tested in the ARIC study for interaction between this SNP and any other SNP genome-wide. We found the eQTL to be involved in a few suggestive interactions, one of which significantly replicated in MESA. Importantly, these gene-gene interactions, involving only rs1532085, explain an additional 1.4% variation of HDL-C, on top of the 0.65% explained by rs1532085 alone. LIPC plays a key role in the lipid metabolism pathway and it, and rs1532085 in particular, has been associated with HDL-C and other lipid levels. Collectively, we discovered several novel gene-gene interactions, all involving an eQTL of LIPC, thus suggesting a hub role of LIPC in the gene-gene interaction network that regulates HDL-C levels, which in turn raises the hypothesis that LIPC's contribution is largely via interactions with other lipid metabolism related genes.

How Gene Networks Can Uncover Novel CVD Players

Cardiovascular diseases (CVD) are complex, involving numerous biological entities from genes and small molecules to organ function. Placing these entities in networks where the functional relationships among the constituents are drawn can aid in our understanding of disease onset, progression and prevention. While networks, or interactomes, are often classified by a general term, say lipids or inflammation, it is a more encompassing class of network that is more informative in showing connections among the active entities and allowing better hypotheses of novel CVD players to be formulated. A range of networks will be presented whereby the potential to bring new objects into the CVD milieu will be exemplified.

Genetic interactions found between calcium channel genes modulate amyloid load measured by positron emission tomography

Late-onset Alzheimer's disease (LOAD) is known to have a complex, oligogenic etiology, with considerable genetic heterogeneity. We investigated the influence of genetic interactions between genes in the Alzheimer's disease (AD) pathway on amyloid-beta (Aβ) deposition as measured by PiB or AV-45 ligand positron emission tomography (PET) to aid in understanding LOAD's genetic etiology. Subsets of the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohorts were used for discovery and for two independent validation analyses. A significant interaction between RYR3 and CACNA1C was confirmed in all three of the independent ADNI datasets. Both genes encode calcium channels expressed in the brain. The results shown here support previous animal studies implicating interactions between these calcium channels in amyloidogenesis and suggest that the pathological cascade of this disease may be modified by interactions in the amyloid-calcium axis. Future work focusing on the mechanisms of such relationships may inform targets for clinical intervention.

How genome-wide SNP-SNP interactions relate to nasopharyngeal carcinoma susceptibility

This study is the first to use genome-wide association study (GWAS) data to evaluate the multidimensional genetic architecture underlying nasopharyngeal cancer. Since analysis of data from GWAS confirms a close and consistent association between elevated risk for nasopharyngeal carcinoma (NPC) and major histocompatibility complex class 1 genes, our goal here was to explore lesser effects of gene-gene interactions. We conducted an exhaustive genome-wide analysis of GWAS data of NPC, revealing two-locus interactions occurring between single nucleotide polymorphisms (SNPs), and identified a number of suggestive interaction loci which were missed by traditional GWAS analyses. Although none of the interaction pairs we identified passed the genome-wide Bonferroni-adjusted threshold for significance, using independent GWAS data from the same population (Stage 2), we selected 66 SNP pairs in 39 clusters with P<0.01. We identified that in several chromosome regions, multiple suggestive interactions group to form a block-like signal, effectively reducing the rate of false discovery. The strongest cluster of interactions involved the CREB5 gene and a SNP rs1607979 on chromosome 17q22 (P = 9.86×10(-11)) which also show trans-expression quantitative loci (eQTL) association in Chinese population. We then detected a complicated cis-interaction pattern around the NPC-associated HLA-B locus, which is immediately adjacent to copy-number variations implicated in male susceptibility for NPC. While it remains to be seen exactly how and to what degree SNP-SNP interactions such as these affect susceptibility for nasopharyngeal cancer, future research on these questions holds great promise for increasing our understanding of this disease's genetic etiology, and possibly also that of other gene-related cancers.

A whole-genome simulator capable of modeling high-order epistasis for complex disease

Genome-wide association studies (GWAS) have been successful in finding numerous new risk variants for complex diseases, but the results almost exclusively rely on single-marker scans. Methods that can analyze joint effects of many variants in GWAS data are still being developed and trialed. To evaluate the performance of such methods it is essential to have a GWAS data simulator that can rapidly simulate a large number of samples, and capture key features of real GWAS data such as linkage disequilibrium (LD) among single-nucleotide polymorphisms (SNPs) and joint effects of multiple loci (multilocus epistasis). In the current study, we combine techniques for specifying high-order epistasis among risk SNPs with an existing program GWAsimulator [Li and Li, 2008] to achieve rapid whole-genome simulation with accurate modeling of complex interactions. We considered various approaches to specifying interaction models including the following: departure from product of marginal effects for pairwise interactions, product terms in logistic regression models for low-order interactions, and penetrance tables conforming to marginal effect constraints for high-order interactions or prescribing known biological interactions. Methods for conversion among different model specifications are developed using penetrance table as the fundamental characterization of disease models. The new program, called simGWA, is capable of efficiently generating large samples of GWAS data with high precision. We show that data simulated by simGWA are faithful to template LD structures, and conform to prespecified diseases models with (or without) interactions.

Properties of local interactions and their potential value in complementing genome-wide association studies

Local interactions between neighbouring SNPs are hypothesized to be able to capture variants missing from genome-wide association studies (GWAS) via haplotype effects but have not been thoroughly explored. We have used a new high-throughput analysis tool to probe this underexplored area through full pair-wise genome scans and conventional GWAS in diastolic and systolic blood pressure and six metabolic traits in the Northern Finland Birth Cohort 1966 (NFBC1966) and the Atherosclerosis Risk in Communities study cohort (ARIC). Genome-wide significant interactions were detected in ARIC for systolic blood pressure between PLEKHA7 (a known GWAS locus for blood pressure) and GPR180 (which plays a role in vascular remodelling), and also for triglycerides as local interactions within the 11q23.3 region (replicated significantly in NFBC1966), which notably harbours several loci (BUD13, ZNF259 and APOA5) contributing to triglyceride levels. Tests of the local interactions within the 11q23.3 region conditional on the top GWAS signal suggested the presence of two independent functional variants, each with supportive evidence for their roles in gene regulation. Local interactions captured 9 additional GWAS loci identified in this study (3 significantly replicated) and 73 from previous GWAS (24 in the eight traits and 49 in related traits). We conclude that the detection of local interactions requires adequate SNP coverage of the genome and that such interactions are only likely to be detectable between SNPs in low linkage disequilibrium. Analysing local interactions is a potentially valuable complement to GWAS and can provide new insights into the biology underlying variation in complex traits.

A novel permutation test for case-only analysis identifies epistatic effects on human longevity in the FOXO gene family

Genetic interactions or epistasis could make a substantial contribution to variation in human complex traits including longevity. However, detecting epistatic interactions in high dimensional datasets is difficult due to various reasons including multiple testing of correlated tests. We introduce a novel permutation strategy to the case-only analysis of gene-by-gene interaction using multiple SNPs. The method is applied to genes coding for Forkhead box O transcription factors which recently have been associated with human longevity across different populations hypothesizing that epistatic interaction in the regulation and expression of the FOXO gene family could contribute to the human longevity phenotype. Genotype data were collected from 1088 individuals from the Danish 1905 birth cohort aged over 92-93 years with 12 SNPs in the FOXO1a and 15 SNPs in the FOXO3a genes. Our analysis detected a joint effect between rs9486902 in FOXO3a and rs2701858 in FOXO1a that highly significantly contributes to human longevity (OR = 3.23, 95% CI: 2.93-3.53) which is consistent in both males and females. Our results were compared with published studies, and importance of our novel method and findings was discussed.