Statistical and Functional Studies Identify Epistasis of Cardiovascular Risk Genomic Variants From Genome‐Wide Association Studies

BACKGROUND: Epistasis describes how gene‐gene interactions affect phenotypes, and could have a profound impact on human diseases such as coronary artery disease (CAD). The goal of this study was to identify gene‐gene interactions in CAD using an easily generalizable multi‐stage approach. METHODS AND RESULTS: Our forward genetic approach consists of multiple steps that combine statistical and functional approaches, and analyze information from global gene expression profiling, functional interactions, and genetic interactions to robustly identify gene‐gene interactions. Global gene expression profiling shows that knockdown of ANRIL (DQ485454) at 9p21.3 GWAS (genome‐wide association studies) CAD locus upregulates TMEM100 and TMEM106B. Functional studies indicate that the increased monocyte adhesion to endothelial cells and transendothelial migration of monocytes, 2 critical processes in the initiation of CAD, by ANRIL knockdown are reversed by knockdown of TMEM106B, but not of TMEM100. Furthermore, the decreased monocyte adhesion to endothelial cells and transendothelial migration of monocytes induced by ANRIL overexpression was reversed by overexpressing TMEM106B. TMEM106B expression was upregulated by >2‐fold in CAD coronary arteries. A significant association was found between variants in TMEM106B (but not in TMEM100) and CAD (P=1.9×10(−8)). Significant gene‐gene interaction was detected between ANRIL variant rs2383207 and TMEM106B variant rs3807865 (P=0.009). A similar approach also identifies significant interaction between rs6903956 in ADTRP and rs17465637 in MIA3 (P=0.005). CONCLUSIONS: We demonstrate 2 pairs of epistatic interactions between GWAS loci for CAD and offer important insights into the genetic architecture and molecular mechanisms for the pathogenesis of CAD. Our strategy has broad applicability to the identification of epistasis in other human diseases.