Background: Recent genome-wide association studies (GWAS) have shown that single nucleotide polymorphisms (SNPs) in the Chr9p21 region are associated with coronary artery disease (CAD). Most of the SNPs identified in this region are non-coding SNPs, suggesting that they may influence gene expression by cis or trans mechanisms to affect disease susceptibility. Since all cells from an individual have the same DNA sequence variations, levels of gene expression in immortalized cell lines can reflect the functional effects of DNA sequence variations that influence or regulate gene expression. The objective of this study is to evaluate the functional consequences of the risk variants in the Chr9p21 region on gene expression.
Methods: We examined the association between the variants in the Chr9p21 region and the transcript-level mRNA expression of the adjacent genes (cis) as well as all other genes across the whole genome (trans) from transformed beta-lymphocytes in 801 non-Hispanic white participants from The Genetic Epidemiology Network of Arteriopathy (GENOA) study.
Results: We found that the CAD risk variants in the Chr9p21 region were significantly associated with the mRNA expression of the ANRIL transcript ENST00000428597 (p = 8.58e-06). Importantly, a few distant transcripts were also found to be associated with the variants in this region, including the well-known CAD risk gene ABCA1 (p = 1.01e-05). Gene enrichment testing suggests that retinol metabolism, N-Glycan biosynthesis, and TGF signaling pathways may be involved.
Conclusion: These results suggest that the effect of risk variants in the Chr9p21 region on susceptibility to CAD is likely to be mediated through both cis and trans mechanisms.