Effect of 9p21.3 coronary artery disease locus neighboring genes on atherosclerosis in mice

Advances and Challenges for GWAS Analysis in Cardiac Diseases: A Focus on Coronary Artery Disease (CAD)

The achievement of genome-wide association studies (GWAS) has rapidly progressed our understanding of the etiology of coronary artery disease (CAD). It unlocks new strategies to strengthen the stalling of CAD drug development. In this review, we highlighted the recent drawbacks, mainly pointing out those involved in identifying causal genes and interpreting the connections between disease pathology and risk variants. We also benchmark the novel insights into the biological mechanism behind the disease primarily based on outcomes of GWAS. Furthermore, we also shed light on the successful discovery of novel treatment targets by introducing various layers of "omics" data and applying systems genetics strategies. Lastly, we discuss in-depth the significance of precision medicine that is helpful to improve through GWAS analysis in cardiovascular research.

Quercetin modulates the liver metabolic profile in a chronic unpredictable mild stress rat model based on metabolomics technology

Depression is the most prevalent psychiatric disease, and its pathogenesis is still unclear. Currently, studies on the pathogenesis of depression are mainly focused on the brain. The liver can modulate brain function via the liver-brain axis, indicating that the liver plays an important role in the development of depression. This study aims to explore the protective effect of quercetin against chronic unpredictable mild stress (CUMS)-induced metabolic changes and the corresponding mechanisms in the rat liver based on untargeted metabolomics technology. In this study, 96 male rats were divided into six groups: control, different doses of quercetin (10 mg per kg bw or 50 mg per kg bw), CUMS, and CUMS + different doses of quercetin. After 8 weeks of CUMS modeling, the liver samples were collected for metabolomics analysis. A total of 17 altered metabolites were identified, including D-glutamic acid, S-adenosylmethionine, lithocholylglycine, L-homocystine, prostaglandin PGE2, leukotriene E4, cholic acid, 5-methyltetrahydrofolic acid, taurochenodeoxycholic acid, S-adenosylhomocysteine, deoxycholic acid, folic acid, L-methionine, leukotriene C5, estriol-17-glucuronide, PE, and PC, indicating that methionine metabolism, bile acid metabolism, and phosphatidylcholine biosynthesis are the major pathways involved in CUMS-induced hepatic metabolic disorders. Hepatic methylation damage may play a role in the pathophysiology of depression, as evidenced by the first discovery of the abnormality of hepatic methionine metabolism. Abnormal changes in hepatic bile acids may provide stronger evidence for depression pathogenesis involving the microbiota-gut-brain axis, suggesting that the liver is involved in depression development and may be a treatment target. The quercetin treatment alleviated the CUMS-induced liver metabolism disorder, suggesting that quercetin may protect against depression by regulating liver metabolism.

Effect of 9p21.3 (lncRNA and CDKN2A/2B) variant on lipid profile

Background

Several 9p21.3 variants, such as rs1333049, rs4977574, rs10757274, rs10757278, and rs10811661, identified from recent genome-wide association studies (GWASs) are reported to be associated with coronary artery disease (CAD) susceptibility but independent of dyslipidemia. This study investigated whether these 9p21.3 variants influenced lipid profiles.

Methods and results

By searching the PubMed and Cochrane databases, 101,099 individuals were included in the analysis. The consistent finding for the rs1333049 C allele on lipid profiles increased the triglyceride (TG) levels. Moreover, the rs4977574 G allele and the rs10757274 G allele, respectively, increased low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) levels. However, the rs10811661 C allele largely reduced LDL-C levels. Subgroup analyses indicated that the effects of the rs1333049 C allele, rs4977574 G allele, and rs10757274 G allele on lipid profiles were stronger in Whites compared with Asians. In contrast, the effect of the rs10811661 C allele on lipid profiles was stronger in Asians compared with Whites.

Conclusion

The rs1333049 C allele, rs4977574 G allele, and rs10757274 G allele of lncRNA, and the rs10811661 G allele of CDKN2A/2B had a significant influence on lipid levels, which may help the understanding of the underlying mechanisms between 9p21.3 variants and CAD.

High shear stress enhances endothelial permeability in the presence of the risk haplotype at 9p21.3

Single nucleotide polymorphisms (SNPs) are exceedingly common in non-coding loci, and while they are significantly associated with a myriad of diseases, their specific impact on cellular dysfunction remains unclear. Here, we show that when exposed to external stressors, the presence of risk SNPs in the 9p21.3 coronary artery disease (CAD) risk locus increases endothelial monolayer and microvessel dysfunction. Endothelial cells (ECs) derived from induced pluripotent stem cells of patients carrying the risk haplotype (R/R WT) differentiated similarly to their non-risk and isogenic knockout (R/R KO) counterparts. Monolayers exhibited greater permeability and reactive oxygen species signaling when the risk haplotype was present. Addition of the inflammatory cytokine TNFα further enhanced EC monolayer permeability but independent of risk haplotype; TNFα also did not substantially alter haplotype transcriptomes. Conversely, when wall shear stress was applied to ECs in a microfluidic vessel, R/R WT vessels were more permeable at lower shear stresses than R/R KO vessels. Transcriptomes of sheared cells clustered more by risk haplotype than by patient or clone, resulting in significant differential regulation of EC adhesion and extracellular matrix genes vs static conditions. A subset of previously identified CAD risk genes invert expression patterns in the presence of high shear concomitant with altered cell adhesion genes, vessel permeability, and endothelial erosion in the presence of the risk haplotype, suggesting that shear stress could be a regulator of non-coding loci with a key impact on CAD.

Genetic Risk Stratification: A Paradigm Shift in Prevention of Coronary Artery Disease

•

CAD is a pandemic that can be prevented. Conventional risk factors are inadequate to detect who is at risk early in the asymptomatic stage.

•

Genetic risk for CAD can be determined at birth, and those at highest genetic risk have been shown to respond to lifestyle changes and statin therapy with a 40% to 50% reduction in cardiac events.

•

Genetic risk stratification for CAD should be brought to the bedside in an attempt to prevent this pandemic disease.

Coronary artery disease (CAD) is a pandemic disease that is highly preventable as shown by secondary prevention. Primary prevention is preferred knowing that 50% of the population can expect a cardiac event in their lifetime. Risk stratification for primary prevention using the American Heart Association/American College of Cardiology predicted 10-year risk based on conventional risk factors for CAD is less than optimal. Conventional risk factors such as hypertension, cholesterol, and age are age-dependent and not present until the sixth or seventh decade of life. The genetic risk score (GRS), which is estimated from the recently discovered genetic variants predisposed to CAD, offers a potential solution to this dilemma. The GRS, which is derived from genotyping the population with a microarray containing these genetic risk variants, has indicated that genetic risk stratification based on the GRS is superior to that of conventional risk factors in detecting those at high risk and who would benefit most from statin therapy. Studies performed in >1 million individuals confirmed genetic risk stratification is superior and primarily independent of conventional risk factors. Prospective clinical trials based on risk stratification for CAD using the GRS have shown lifestyle changes, physical activity, and statin therapy are associated with 40% to 50% reduction in cardiac events in the high genetic risk group (20%). Genetic risk stratification has the advantage of being innate to an individual’s DNA, and because DNA does not change in a lifetime, it is independent of age. Genetic risk stratification is inexpensive and can be performed worldwide, providing risk analysis at any age and thus has the potential to revolutionize primary prevention.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

The SNP rs7865618 of 9p21.3 locus emerges as the most promising marker of coronary artery disease in the southern Indian population

Development of coronary artery disease (CAD) is primarily due to the process of atherosclerosis, however the prognosis of CAD depends on pleiotropic effects of the genes located at 9p21.3 region. Genome wide association studies revealed association of variants in this region with CAD pathology. However, specific marker in predicting CAD development or progression is not yet identified. In the present study, 35 SNPs at 9p21.3 region, located in the cyclin dependent kinase inhibitor (CDKN2A/CDKN2B) genes, were genotyped among 350 CAD cases and 480 controls from the southern Indian population of Hyderabad using fluidigm nanofluidic SNP genotyping system and the data were analyzed using PLINK and R softwares. Of the 35 SNPs analysed, only one SNP, rs7865618, was found to be highly significantly associated with CAD, even after correction for multiple testing (p = 0.008). The AG and GG genotypes of this SNP conferred 3.08 and 1.93 folds increased risk for CAD respectively. In particular, this SNP was significantly associated with severe anatomic (triple vessel disease p = 0.023) and phenotypic (acute coronary syndrome p = 0.007) categories of CAD. Pair wise SNP interaction analysis between the SNPs of 9p21.3 and 11q23.3 regions revealed significantly increased risk of three SNPs of 11q23.3 region that were not associated individually, in conjunction with rs7865618 of 9p21.3.

Over 1000 genetic loci influencing blood pressure with multiple systems and tissues implicated

High blood pressure (BP) remains the major heritable and modifiable risk factor for cardiovascular disease. Persistent high BP, or hypertension, is a complex trait with both genetic and environmental interactions. Despite swift advances in genomics, translating new discoveries to further our understanding of the underlying molecular mechanisms remains a challenge. More than 500 loci implicated in the regulation of BP have been revealed by genome-wide association studies (GWAS) in 2018 alone, taking the total number of BP genetic loci to over 1000. Even with the large number of loci now associated to BP, the genetic variance explained by all loci together remains low (~5.7%). These genetic associations have elucidated mechanisms and pathways regulating BP, highlighting potential new therapeutic and drug repurposing targets. A large proportion of the BP loci were discovered and reported simultaneously by multiple research groups, creating a knowledge gap, where the reported loci to date have not been investigated in a harmonious way. Here, we review the BP-associated genetic variants reported across GWAS studies and investigate their potential impact on the biological systems using in silico enrichment analyses for pathways, tissues, gene ontology and genetic pleiotropy.

Mechanical influences on cardiovascular differentiation and disease modeling

Tissues are continuously exposed to forces in vivo, whether from fluid pressure in an artery from our blood or compressive forces on joints from our body weight. The forces that cells are exposed to arise almost immediately after conception; it is therefore important to understand how forces shape stem cell differentiation into lineage committed cells, how they help organize cells into tissues, and how forces can cause or exacerbate disease. No tissue is exempt, but cardiovascular tissues in particular are exposed to these forces. While animal models have been used extensively in the past, there is growing recognition of their limitations when modeling disease complexity or human genetics. In this mini review, we summarize current understanding of the mechanical influences on the differentiation of cardiovascular progeny, how the transduction of forces influence the onset of disease, and how engineering approaches applied to this problem have yielded systems that create mature-like human tissues in vitro in which to assess the impact of disease on cell function.

Mechanical activation of noncoding-RNA-mediated regulation of disease-associated phenotypes in human cardiomyocytes

How common polymorphisms in noncoding genome regions can regulate cellular function remains largely unknown. Here we show that cardiac fibrosis, mimicked using a hydrogel with controllable stiffness, affects the regulation of the phenotypes of human cardiomyocytes by a portion of the long noncoding RNA ANRIL, the gene of which is located in the disease-associated 9p21 locus. In a physiological environment, cultured cardiomyocytes derived from induced pluripotent stem cells obtained from patients who are homozygous for cardiovascular-risk alleles (R/R cardiomyocytes) or from healthy individuals who are homozygous for nonrisk alleles contracted synchronously, independently of genotype. After hydrogel stiffening to mimic fibrosis, only the R/R cardiomyocytes exhibited asynchronous contractions. These effects were associated with increased expression of the short ANRIL isoform in R/R cardiomyocytes, which induced a c-Jun N-terminal kinase (JNK) phosphorylation-based mechanism that impaired gap junctions (particularly, loss of connexin-43 expression) following stiffening. Deletion of the risk locus or treatment with a JNK antagonist was sufficient to maintain gap junctions and prevent asynchronous contraction of cardiomyocytes. Our findings suggest that mechanical changes in the microenvironment of cardiomyocytes can activate the regulation of their function by noncoding loci.

Long Noncoding RNA ANRIL: Lnc-ing Genetic Variation at the Chromosome 9p21 Locus to Molecular Mechanisms of Atherosclerosis

Ever since the first genome-wide association studies (GWAS) on coronary artery disease (CAD), the Chr9p21 risk locus has emerged as a top signal in GWAS of atherosclerotic cardiovascular disease, including stroke and peripheral artery disease. The CAD risk SNPs on Chr9p21 lie within a stretch of 58 kilobases of non-protein-coding DNA, containing the gene body of the long noncoding RNA (lncRNA) antisense non coding RNA in the INK4 locus (ANRIL). How risk is affected by the Chr9p21 locus in molecular detail is a matter of ongoing research. Here we will review recent advances in the understanding that ANRIL serves as a key risk effector molecule of atherogenesis at the locus. One focus of this review is the shift in understanding that genetic variation at Chr9p21 not only affects the abundance of ANRIL, and in some cases expression of the adjacent CDKN2A/B tumor suppressors, but also impacts ANRIL splicing, such that 3′-5′-linked circular noncoding ANRIL RNA species are produced. We describe how the balance of linear and circular ANRIL RNA, determined by the Chr9p21 genotype, regulates molecular pathways and cellular functions involved in atherogenesis. We end with an outlook on how manipulating circular ANRIL abundance may be exploited for therapeutic purposes.

Translating GWAS Findings to Novel Therapeutic Targets for Coronary Artery Disease

The success of genome-wide association studies (GWAS) has significantly advanced our understanding of the etiology of coronary artery disease (CAD) and opens new opportunities to reinvigorate the stalling CAD drug development. However, there exists remarkable disconnection between the CAD GWAS findings and commercialized drugs. While this could implicate major untapped translational and therapeutic potentials in CAD GWAS, it also brings forward extensive technical challenges. In this review we summarize the motivation to leverage GWAS for drug discovery, outline the critical bottlenecks in the field, and highlight several promising strategies such as functional genomics and network-based approaches to enhance the translational value of CAD GWAS findings in driving novel therapeutics

Genetics and Genomics of Coronary Artery Disease

Coronary artery disease (or coronary heart disease), is the leading cause of mortality in many of the developing as well as the developed countries of the world. Cholesterol-enriched plaques in the heart's blood vessels combined with inflammation lead to the lesion expansion, narrowing of blood vessels, reduced blood flow, and may subsequently cause lesion rupture and a heart attack. Even though several environmental risk factors have been established, such as high LDL-cholesterol, diabetes, and high blood pressure, the underlying genetic composition may substantially modify the disease risk; hence, genome composition and gene-environment interactions may be critical for disease progression. Ongoing scientific efforts have seen substantial advancements related to the fields of genetics and genomics, with the major breakthroughs yet to come. As genomics is the most rapidly advancing field in the life sciences, it is important to present a comprehensive overview of current efforts. Here, we present a summary of various genetic and genomics assays and approaches applied to coronary artery disease research.

Enhanced Megakaryopoiesis and Platelet Activity in Hypercholesterolemic, B6-Ldlr-/-, Cdkn2a-Deficient Mice

BACKGROUND

Genome-wide association studies for coronary artery disease/myocardial infarction revealed a 58 kb risk locus on 9p21.3. Refined genetic analyses revealed unique haplotype blocks conferring susceptibility to atherosclerosis per se versus risk for acute complications in the presence of underlying coronary artery disease. The cell proliferation inhibitor locus, CDKN2A, maps just upstream of the myocardial infarction risk block, is at least partly regulated by the noncoding RNA, ANRIL, overlapping the risk block, and has been associated with platelet counts in humans. Thus, we tested the hypothesis that CDKN2A deficiency predisposes to increased platelet production, leading to increased platelet activation in the setting of hypercholesterolemia.

METHODS AND RESULTS

Platelet production and activation were measured in B6-Ldlr(-/-)Cdkn2a(+/-) mice and a congenic strain carrying the region of homology with the human 9p21.3/CDKN2A locus. The strains exhibit decreased expression of CDKN2A (both p16(INK4a) and p19(ARF)) but not CDKN2B (p15(INK4b)). Compared with B6-Ldlr(-/-) controls, both Cdkn2a-deficient strains exhibited increased platelet counts and bone marrow megakaryopoiesis. The platelet overproduction phenotype was reversed by treatment with cyclin-dependent kinase 4/6 inhibitor, PD0332991/palbociclib, that mimics the endogenous effect of p16(INK4a). Western diet feeding resulted in increased platelet activation, increased thrombin/antithrombin complex, and decreased bleeding times in Cdkn2a-deficient mice compared with controls.

CONCLUSIONS

Together, the data suggest that one or more Cdkn2a transcripts modulate platelet production and activity in the setting of hypercholesterolemia, amenable to pharmaceutical intervention. Enhanced platelet production and activation may predispose to arterial thrombosis, suggesting an explanation, at least in part, for the association of 9p21.3 and myocardial infarction.

Experimental Biology for the Identification of Causal Pathways in Atherosclerosis

More than 60 genomic loci have been implicated by genome-wide association studies (GWAS) and exome-wide association studies as conferring an increased risk of myocardial infarction and coronary artery disease (CAD). However, the causal gene and variant is often unclear. Using the functional analysis of genetic variants in experimental animal models, we anticipate understanding which candidate gene at a specific locus is associated with atherosclerosis and revealing the underlying molecular and cellular mechanisms, ultimately leading to the identification of causal pathways in atherosclerosis and may provide novel therapeutic targets for the treatment of atherosclerotic cardiovascular disease.

Genetic Architecture of Atherosclerosis in Mice: A Systems Genetics Analysis of Common Inbred Strains

Common forms of atherosclerosis involve multiple genetic and environmental factors. While human genome-wide association studies have identified numerous loci contributing to coronary artery disease and its risk factors, these studies are unable to control environmental factors or examine detailed molecular traits in relevant tissues. We now report a study of natural variations contributing to atherosclerosis and related traits in over 100 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP). The mice were made hyperlipidemic by transgenic expression of human apolipoprotein E-Leiden (APOE-Leiden) and human cholesteryl ester transfer protein (CETP). The mice were examined for lesion size and morphology as well as plasma lipid, insulin and glucose levels, and blood cell profiles. A subset of mice was studied for plasma levels of metabolites and cytokines. We also measured global transcript levels in aorta and liver. Finally, the uptake of acetylated LDL by macrophages from HMDP mice was quantitatively examined. Loci contributing to the traits were mapped using association analysis, and relationships among traits were examined using correlation and statistical modeling. A number of conclusions emerged. First, relationships among atherosclerosis and the risk factors in mice resemble those found in humans. Second, a number of trait-loci were identified, including some overlapping with previous human and mouse studies. Third, gene expression data enabled enrichment analysis of pathways contributing to atherosclerosis and prioritization of candidate genes at associated loci in both mice and humans. Fourth, the data provided a number of mechanistic inferences; for example, we detected no association between macrophage uptake of acetylated LDL and atherosclerosis. Fifth, broad sense heritability for atherosclerosis was much larger than narrow sense heritability, indicating an important role for gene-by-gene interactions. Sixth, stepwise linear regression showed that the combined variations in plasma metabolites, including LDL/VLDL-cholesterol, trimethylamine N-oxide (TMAO), arginine, glucose and insulin, account for approximately 30 to 40% of the variation in atherosclerotic lesion area. Overall, our data provide a rich resource for studies of complex interactions underlying atherosclerosis.

While recent genetic association studies in human populations have succeeded in identifying genetic loci that contribute to coronary artery disease (CAD) and related phenotypes, these loci explain only a small fraction of the genetic variation in CAD and associated traits. Here, we present a complementary approach using association analysis of atherosclerotic traits among inbred strains of mice. A strength of this approach is that it enables in-depth phenotypic characterization including gene expression and metabolic profiling across a variety of tissues, and integration of these molecular phenotypes with coronary artery disease itself. A striking finding was the large fraction of atherosclerosis that was explained by genetic interactions. Association analysis allowed us to identify genetic loci for atherosclerotic lesion area as well as transcript, cytokine and metabolite levels, and relationships among the traits were examined by correlation and network modeling. The plasma metabolites associated with atherosclerosis in mice, namely, LDL/VLDL-cholesterol, TMAO, arginine, glucose and insulin, overlapped with those observed in humans and accounted for approximately 30 to 40% of the observed variation in atherosclerotic lesion area. In summary, our data provide a detailed overview of the genetic architecture of atherosclerosis in mice and a rich resource for studies of the complex genetic and metabolic interactions that underlie the disease.

Functional genomics of the CDKN2A/B locus in cardiovascular and metabolic disease: what have we learned from GWASs?

Genome-wide association studies (GWASs) provide an unprecedented opportunity to examine, on a large scale, the association of common genetic variants with complex diseases like type 2 diabetes (T2D) and cardiovascular disease (CVD), thus allowing the identification of new potential disease loci. Using this approach, numerous studies have associated SNPs on chromosome 9p21.3 situated near the cyclin-dependent kinase inhibitor 2A/B (CDKN2A/B) locus with the risk for coronary artery disease (CAD) and T2D. However, identifying the function of the nearby gene products (CDKN2A/B and ANRIL) in the pathophysiology of these conditions requires functional genomic studies. We review the current knowledge, from studies using human and mouse models, describing the function of CDKN2A/B gene products, which may mechanistically link the 9p21.3 risk locus with CVD and diabetes.

Targeting metabolic disturbance in the diabetic heart

Diabetic cardiomyopathy is defined as ventricular dysfunction initiated by alterations in cardiac energy substrates in the absence of coronary artery disease and hypertension. In addition to the demonstrated burden of cardiovascular events associated with diabetes, diabetic cardiomyopathy partly explains why diabetic patients are subject to a greater risk of heart failure and a worse outcome after myocardial ischemia. The raising prevalence and accumulating costs of cardiovascular disease in diabetic patients underscore the deficiencies of tertiary prevention and call for a shift in medical treatment. It is becoming increasingly clearer that the effective prevention and treatment of diabetic cardiomyopathy require measures to regulate the metabolic derangement occurring in the heart rather than merely restoring suitable systemic parameters. Recent research has provided deeper insight into the metabolic etiology of diabetic cardiomyopathy and numerous heart-specific targets that may substitute or reinforce current strategies. From both experimental and translational perspectives, in this review we first discuss the progress made with conventional therapies, and then focus on the need for prospective metabolic targets that may avert myocardial vulnerability and functional decline in next-generation diabetic care.

Functional genomics of the 9p21.3 locus for atherosclerosis: clarity or confusion?

The 9p21.3 locus was the first to yield to genome-wide association studies (GWAS) seeking common genetic variants predisposing to increased risk of coronary artery atherosclerotic disease (CAD). The 59 single nucleotide polymorphisms that show highest association with CAD are clustered in a region 100,000 to 150,000 base pairs 5' to the cyclin-dependent kinase inhibitors CDKN2B (coding for p15(ink4b)) and CDKN2A (coding for p16(ink4a) and p14(ARF)). This region also covers the 3' end of a long noncoding RNA transcribed antisense to CDKN2B (CDKN2BAS, aka ANRIL for antisense noncoding RNA at the ink4 locus) whose expression has been linked to chromatin remodeling at the locus. Despite intensive investigation over the past 7 years, the functional significance of the 9p21.3 locus remains elusive. Other variants at this locus have been associated with glaucoma, glioma, and type 2 diabetes mellitus, diseases that implicate tissue-resident macrophages. Here, we review the evidence that genetic variants at 9p21.3 disrupt tissue-specific enhancers and propose new insights to guide future studies.

The role of genetic risk factors in coronary artery disease

Genome-wide association studies for coronary artery disease utilizing the case control association study approach has identified 50 genetic risk variants associated with coronary artery disease or myocardial infarction. All of these genetic variants are of genome wide significance and replicated in an independent population. It is of note that 35 of these 50 genetic risk variants act through mechanisms as yet unknown. These findings have great implications for the pathogenesis of atherosclerosis, as well as new targets for the development of novel therapies for the prevention and treatment of CAD. The genetic variant PCSK9 has already led to the development of a monoclonal anti-body which is undergoing assessment in phases I, II, and III clinical trials. This therapy shows very promising results and since it increases removal of LDL-C, it is complementary to current statin therapy. Assessing the beneficial or deleterious effects of a lifelong exposure to a genetic risk variant (Mendelian randomization) will be an important adjunct to clinical trials.

Genetic variants at chromosome 9p21 and risk of first versus subsequent coronary heart disease events: a systematic review and meta-analysis

Objectives

The purpose of this analysis was to compare the association between variants at the chromosome 9p21 locus (Ch9p21) and risk of first versus subsequent coronary heart disease (CHD) events through systematic review and meta-analysis.

Background

Ch9p21 is a recognized risk factor for a first CHD event. However, its association with risk of subsequent events in patients with established CHD is less clear.

Methods

We searched PubMed and EMBASE for prospective studies reporting association of Ch9p21 with incident CHD events and extracted information on cohort type (individuals without prior CHD or individuals with established CHD) and effect estimates for risk of events.

Results

We identified 31 cohorts reporting on 193,372 individuals. Among the 16 cohorts of individuals without prior CHD (n = 168,209), there were 15,664 first CHD events. Ch9p21 was associated with a pooled hazard ratio (HR) of a first event of 1.19 (95% confidence interval: 1.17 to 1.22) per risk allele. In individuals with established CHD (n = 25,163), there were 4,436 subsequent events providing >99% and 91% power to detect a per-allele HR of 1.19 or 1.10, respectively. The pooled HR for subsequent events was 1.01 (95% confidence interval: 0.97 to 1.06) per risk allele. There was strong evidence of heterogeneity between the effect estimates for first and subsequent events (p value for heterogeneity = 5.6 × 10⁻¹¹). We found no evidence for biases to account for these findings.

Conclusions

Ch9p21 shows differential association with risk of first versus subsequent CHD events. This has implications for genetic risk prediction in patients with established CHD and for mechanistic understanding of how Ch9p21 influences risk of CHD.

Cyclin-dependent kinase inhibitor 2B regulates efferocytosis and atherosclerosis

Genetic variation at the chromosome 9p21 risk locus promotes cardiovascular disease; however, it is unclear how or which proteins encoded at this locus contribute to disease. We have previously demonstrated that loss of one candidate gene at this locus, cyclin-dependent kinase inhibitor 2B (Cdkn2b), in mice promotes vascular SMC apoptosis and aneurysm progression. Here, we investigated the role of Cdnk2b in atherogenesis and found that in a mouse model of atherosclerosis, deletion of Cdnk2b promoted advanced development of atherosclerotic plaques composed of large necrotic cores. Furthermore, human carriers of the 9p21 risk allele had reduced expression of CDKN2B in atherosclerotic plaques, which was associated with impaired expression of calreticulin, a ligand required for activation of engulfment receptors on phagocytic cells. As a result of decreased calreticulin, CDKN2B-deficient apoptotic bodies were resistant to efferocytosis and not efficiently cleared by neighboring macrophages. These uncleared SMCs elicited a series of proatherogenic juxtacrine responses associated with increased foam cell formation and inflammatory cytokine elaboration. The addition of exogenous calreticulin reversed defects associated with loss of Cdkn2b and normalized engulfment of Cdkn2b-deficient cells. Together, these data suggest that loss of CDKN2B promotes atherosclerosis by increasing the size and complexity of the lipid-laden necrotic core through impaired efferocytosis.

Cyclin-dependent kinase inhibitor 2B regulates efferocytosis and atherosclerosis

Genetic variation at the chromosome 9p21 risk locus promotes cardiovascular disease; however, it is unclear how or which proteins encoded at this locus contribute to disease. We have previously demonstrated that loss of one candidate gene at this locus, cyclin-dependent kinase inhibitor 2B (Cdkn2b), in mice promotes vascular SMC apoptosis and aneurysm progression. Here, we investigated the role of Cdnk2b in atherogenesis and found that in a mouse model of atherosclerosis, deletion of Cdnk2b promoted advanced development of atherosclerotic plaques composed of large necrotic cores. Furthermore, human carriers of the 9p21 risk allele had reduced expression of CDKN2B in atherosclerotic plaques, which was associated with impaired expression of calreticulin, a ligand required for activation of engulfment receptors on phagocytic cells. As a result of decreased calreticulin, CDKN2B-deficient apoptotic bodies were resistant to efferocytosis and not efficiently cleared by neighboring macrophages. These uncleared SMCs elicited a series of proatherogenic juxtacrine responses associated with increased foam cell formation and inflammatory cytokine elaboration. The addition of exogenous calreticulin reversed defects associated with loss of Cdkn2b and normalized engulfment of Cdkn2b-deficient cells. Together, these data suggest that loss of CDKN2B promotes atherosclerosis by increasing the size and complexity of the lipid-laden necrotic core through impaired efferocytosis.

Genetic analysis of the 9p21.3 CAD risk locus in Asian Indians

The 9p21.3 locus is the best replicated region to date for coronary artery disease (CAD). We investigated the association of 9p21.3 common variants with CAD, candidate gene expression including ANRIL, a non-coding RNA, followed by in vitro validation. Five variants, rs10757278, rs10757274, rs2383206, rs1333049 and rs4977574 were genotyped in 1,034 cases and 1,034 controls. Gene expression of C9orf5, MTAP1, MTAP 2, p16INK4a, p14ARF, p15INK4b and two ANRIL splice variants, NR_003529 and EU741058, were measured in 100 cases and 100 controls. Human aortic smooth muscle cells (HuAoSMCs) were transfected with siRNA targeting ANRIL exon 19 (siRNA-1) or exon 2 (siRNA-2) and consequent effect determined. rs2383206 showed the highest association with CAD (odds ratio [OR] 2.02, 95% confidence interval [CI] 1.56 -2.62) and an adjusted OR of 2.55, 1.33-2.88 along with rs10757278. Conventional risk factors (conventional RFs), rs2383206 and rs10757278 variants together yielded a higher c index (OR 0.790, 95% CI 0.770 -0.810) as compared to conventional RFs (OR 0.783, 95% CI 0.763-0.803) or genetic variants (OR 0.561, 95% CI 0.536-0.586) alone. GAAAA haplotype showed significant protective association with CAD compared to CGGGG risk haplotype (OR 0.45, 95% CI 0.27-0.77). Expression of p16INK4a, p14ARF and p15INK4b as well as plasma CDKN2A levels were lower in cases than controls. GG genotype was associated with higher EU741058 expression and lower p16INK4a expression. HuAoSMCs transfected with siRNA-1 showed lower NR_003529, p16INK4aand p14ARFexpression. Our study provides further evidence on the significance of 9p21.3 locus for CAD wherein the risk allele regulate the expression of ANRIL and adjacent tumour suppressor genes which in turn alter smooth muscle proliferation, a fundamental process in atherosclerosis.

Nutrigenetics and nutrigenomics of atherosclerosis

The latest genome-wide association studies (GWAS) have re-energized our effort to understand the genetic basis of atherosclerotic cardiovascular disease. Although the knowledge generated by GWAS has confirmed that mediators of inflammation and perturbed lipid metabolism are major players in cardiovascular disease (CVD) development, much of individual disease heritability remains unexplained by the variants identified through GWAS. Moreover, results from interventions that aim at the pharmaceutical modification of lipid parameters fall short of expectation. These elusive treatment goals based on heritability studies highlight a key supportive, and perhaps even primary, role of nutritional therapy to achieve better health outcomes. Nonetheless, effective and specific interventions for CVD prevention using principles of "personalized" nutrition require a better knowledge of gene-diet interactions, an area that remains poorly explored. Dietary fatty acids such as omega-3 polyunsaturated fatty acids (PUFAs) are an excellent example of a widely studied "environment" that interacts with the genetic makeup in relation to CVD. A thorough exploration of the nutrigenomics and nutrigenetics of omega-3 PUFAs is key to understanding the etiology, and developing effective preventive measures. In this review, we will summarize the current state of knowledge of genetic interactions with omega-3 PUFAs in modulating lipid metabolism and inflammation, and defining health outcomes. Nutrigenetics and nutrigenomics are still in their infancy with respect to CVD prediction and therapy. Integration of the progress in the omics, including metabolomics, lipidomics, transcriptomics, and proteomics, coupled with advances in nutrigenomic and nutrigenetic research will move us towards personalized medicine as the ultimate paradigm of responsible clinical practice.

Amplification and thrifty single-molecule sequencing of recurrent somatic structural variations

Deletion of tumor-suppressor genes as well as other genomic rearrangements pervade cancer genomes across numerous types of solid tumor and hematologic malignancies. However, even for a specific rearrangement, the breakpoints may vary between individuals, such as the recurrent CDKN2A deletion. Characterizing the exact breakpoints for structural variants (SVs) is useful for designating patient-specific tumor biomarkers. We propose AmBre (Amplification of Breakpoints), a method to target SV breakpoints occurring in samples composed of heterogeneous tumor and germline DNA. Additionally, AmBre validates SVs called by whole-exome/genome sequencing and hybridization arrays. AmBre involves a PCR-based approach to amplify the DNA segment containing an SV's breakpoint and then confirms breakpoints using sequencing by Pacific Biosciences RS. To amplify breakpoints with PCR, primers tiling specified target regions are carefully selected with a simulated annealing algorithm to minimize off-target amplification and maximize efficiency at capturing all possible breakpoints within the target regions. To confirm correct amplification and obtain breakpoints, PCR amplicons are combined without barcoding and simultaneously long-read sequenced using a single SMRT cell. Our algorithm efficiently separates reads based on breakpoints. Each read group supporting the same breakpoint corresponds with an amplicon and a consensus amplicon sequence is called. AmBre was used to discover CDKN2A deletion breakpoints in cancer cell lines: A549, CEM, Detroit562, MOLT4, MCF7, and T98G. Also, we successfully assayed RUNX1–RUNX1T1 reciprocal translocations by finding both breakpoints in the Kasumi-1 cell line. AmBre successfully targets SVs where DNA harboring the breakpoints are present in 1:1000 mixtures.

Molecular biology of atherosclerosis

At least 468 individual genes have been manipulated by molecular methods to study their effects on the initiation, promotion, and progression of atherosclerosis. Most clinicians and many investigators, even in related disciplines, find many of these genes and the related pathways entirely foreign. Medical schools generally do not attempt to incorporate the relevant molecular biology into their curriculum. A number of key signaling pathways are highly relevant to atherogenesis and are presented to provide a context for the gene manipulations summarized herein. The pathways include the following: the insulin receptor (and other receptor tyrosine kinases); Ras and MAPK activation; TNF-α and related family members leading to activation of NF-κB; effects of reactive oxygen species (ROS) on signaling; endothelial adaptations to flow including G protein-coupled receptor (GPCR) and integrin-related signaling; activation of endothelial and other cells by modified lipoproteins; purinergic signaling; control of leukocyte adhesion to endothelium, migration, and further activation; foam cell formation; and macrophage and vascular smooth muscle cell signaling related to proliferation, efferocytosis, and apoptosis. This review is intended primarily as an introduction to these key signaling pathways. They have become the focus of modern atherosclerosis research and will undoubtedly provide a rich resource for future innovation toward intervention and prevention of the number one cause of death in the modern world.

CDKN2B expression in adipose tissue of familial combined hyperlipidemia patients

The purpose of this study was to determine the core biological processes perturbed in the subcutaneous adipose tissue of familial combined hyperlipidemia (FCHL) patients. Annotation of FCHL and control microarray datasets revealed a distinctive FCHL transcriptome, characterized by gene expression changes regulating five overlapping systems: the cytoskeleton, cell adhesion and extracellular matrix; vesicular trafficking; lipid homeostasis; and cell cycle and apoptosis. Expression values for the cell-cycle inhibitor CDKN2B were increased, replicating data from an independent FCHL cohort. In 3T3-L1 cells, CDKN2B knockdown induced C/EBPα expression and lipid accumulation. The minor allele at SNP site rs1063192 (C) was predicted to create a perfect seed for the human miRNA-323b-5p. A miR-323b-5p mimic significantly reduced endogenous CDKN2B protein levels and the activity of a CDKN2B 3'UTR luciferase reporter carrying the rs1063192 C allele. Although the allele displayed suggestive evidence of association with reduced CDKN2B mRNA in the MuTHER adipose tissue dataset, family studies suggest the association between increased CDKN2B expression and FCHL-lipid abnormalities is driven by factors external to this gene locus. In conclusion, from a comparative annotation analysis of two separate FCHL adipose tissue transcriptomes and a subsequent focus on CDKN2B, we propose that dysfunctional adipogenesis forms an integral part of FCHL pathogenesis.

From genotype to phenotype in human atherosclerosis--recent findings

PURPOSE OF REVIEW

Since 2007, genome-wide association studies (GWAS) have led to the identification of numerous loci of atherosclerotic cardiovascular disease. The majority of these loci harbor genes previously not known to be involved in atherogenesis. In this review, we summarize the recent progress in understanding the pathophysiology of genetic variants in atherosclerosis.

RECENT FINDINGS

Fifty-eight loci with P < 10⁻⁷ have been identified in GWAS for coronary heart disease and myocardial infarction. Of these, 23 loci (40%) overlap with GWAS loci of classical risk factors such as lipids, blood pressure, and diabetes mellitus, suggesting a potential causal relation. The vast majority of the remaining 35 loci (60%) are at genomic regions where the mechanism in atherogenesis is unclear. Loci most frequently found in independent GWAS were at Chr9p21.3 (ANRIL/CDKN2B-AS1), Chr6p24.1 (PHACTR1), and Chr1p13.3 (CELSR2, PSRC1, MYBPHL, SORT1). Recent work suggests that Chr9p21.3 exerts its effects through epigenetic regulation of target genes, whereas mechanisms at Chr6p24.1 remain obscure, and Chr1p13.3 affects plasma LDL cholesterol.

SUMMARY

Novel GWAS loci indicate that our understanding of atherosclerosis is limited and implicate a role of hitherto unknown mechanisms, such as epigenetic gene regulation in atherogenesis.

Association between rs10118757(A/G) in methylthioadenosine phosphorylase gene and coronary artery disease in Chinese Hans

Studies focusing on the association of gene methylthioadenosine phosphorylase (MTAP) with the risk of coronary artery disease (CAD) and myocardial infarction (MI) are limited. In this study, we explored the effects of rs10118757 in MTAP gene on CAD and MI by performing association analysis in a Chinese Han population. rs10118757 was genotyped in 1007 CAD patients (including 338 MI patients) and 885 healthy controls. Allelic analysis showed that allele A of rs10118757 was associated with increased risk of CAD, with OR (95%CI)=1.193 (1.035-1.376), and P=0.015. After adjusted for age, BMI, gender, hypertension and smoking, rs10118757 was still significantly associated with CAD under additive and dominant models, with OR (95%CI)=1.252 (1.070-1.465), P=0.005, and OR (95%CI)=1.698 (1.168-2.467), P=0.006, respectively. Compared to additive model, dominant model may be the best-fitting model (P=6.63E-10 vs P=6.70E-10). As reported previously, rs10118757 was not associated with MI in the current study. Our study firstly reported that SNP rs10118757 was associated with CAD risk in a Chinese Han population, indicating that MTAP gene may play a potential role in the pathophysiological process of CAD.

Alu elements in ANRIL non-coding RNA at chromosome 9p21 modulate atherogenic cell functions through trans-regulation of gene networks

The chromosome 9p21 (Chr9p21) locus of coronary artery disease has been identified in the first surge of genome-wide association and is the strongest genetic factor of atherosclerosis known today. Chr9p21 encodes the long non-coding RNA (ncRNA) antisense non-coding RNA in the INK4 locus (ANRIL). ANRIL expression is associated with the Chr9p21 genotype and correlated with atherosclerosis severity. Here, we report on the molecular mechanisms through which ANRIL regulates target-genes in trans, leading to increased cell proliferation, increased cell adhesion and decreased apoptosis, which are all essential mechanisms of atherogenesis. Importantly, trans-regulation was dependent on Alu motifs, which marked the promoters of ANRIL target genes and were mirrored in ANRIL RNA transcripts. ANRIL bound Polycomb group proteins that were highly enriched in the proximity of Alu motifs across the genome and were recruited to promoters of target genes upon ANRIL over-expression. The functional relevance of Alu motifs in ANRIL was confirmed by deletion and mutagenesis, reversing trans-regulation and atherogenic cell functions. ANRIL-regulated networks were confirmed in 2280 individuals with and without coronary artery disease and functionally validated in primary cells from patients carrying the Chr9p21 risk allele. Our study provides a molecular mechanism for pro-atherogenic effects of ANRIL at Chr9p21 and suggests a novel role for Alu elements in epigenetic gene regulation by long ncRNAs.

Chromosome 9p21 is the strongest genetic factor for coronary artery disease and encodes the long non-coding RNA (ncRNA) ANRIL. Here, we show that increased ANRIL expression mediates atherosclerosis risk through trans-regulation of gene networks leading to pro-atherogenic cellular properties, such as increased proliferation and adhesion. ANRIL may act as a scaffold, guiding effector-proteins to chromatin. These functions depend on an Alu motif present in ANRIL RNA and mirrored several thousand-fold in the genome. Alu elements are a family of primate-specific short interspersed repeat elements (SINEs) and have been linked with genetic disease. Current models propose that either exonisation of Alu elements or changes of cis-regulation of adjacent genes are the underlying disease mechanisms. Our work extends the function of Alu transposons to regulatory components of long ncRNAs with a central role in epigenetic trans-regulation. Furthermore, it implies a pivotal role for Alu elements in genetically determined vascular disease and describes a plausible molecular mechanism for a pro-atherogenic function of ANRIL at chromosome 9p21.

Association between the chromosome 9p21 locus and angiographic coronary artery disease burden: a collaborative meta-analysis

OBJECTIVES

This study sought to ascertain the relationship of 9p21 locus with: 1) angiographic coronary artery disease (CAD) burden; and 2) myocardial infarction (MI) in individuals with underlying CAD.

BACKGROUND

Chromosome 9p21 variants have been robustly associated with coronary heart disease, but questions remain on the mechanism of risk, specifically whether the locus contributes to coronary atheroma burden or plaque instability.

METHODS

We established a collaboration of 21 studies consisting of 33,673 subjects with information on both CAD (clinical or angiographic) and MI status along with 9p21 genotype. Tabular data are provided for each cohort on the presence and burden of angiographic CAD, MI cases with underlying CAD, and the diabetic status of all subjects.

RESULTS

We first confirmed an association between 9p21 and CAD with angiographically defined cases and control subjects (pooled odds ratio [OR]: 1.31, 95% confidence interval [CI]: 1.20 to 1.43). Among subjects with angiographic CAD (n = 20,987), random-effects model identified an association with multivessel CAD, compared with those with single-vessel disease (OR: 1.10, 95% CI: 1.04 to 1.17)/copy of risk allele). Genotypic models showed an OR of 1.15, 95% CI: 1.04 to 1.26 for heterozygous carrier and OR: 1.23, 95% CI: 1.08 to 1.39 for homozygous carrier. Finally, there was no significant association between 9p21 and prevalent MI when both cases (n = 17,791) and control subjects (n = 15,882) had underlying CAD (OR: 0.99, 95% CI: 0.95 to 1.03)/risk allele.

CONCLUSIONS

The 9p21 locus shows convincing association with greater burden of CAD but not with MI in the presence of underlying CAD. This adds further weight to the hypothesis that 9p21 locus primarily mediates an atherosclerotic phenotype.

ANRIL: molecular mechanisms and implications in human health

ANRIL is a recently discovered long non-coding RNA encoded in the chromosome 9p21 region. This locus is a hotspot for disease-associated polymorphisms, and it has been consistently associated with cardiovascular disease, and more recently with several cancers, diabetes, glaucoma, endometriosis among other conditions. ANRIL has been shown to regulate its neighbor tumor suppressors CDKN2A/B by epigenetic mechanisms and thereby regulate cell proliferation and senescence. However, the clear role of ANRIL in the pathogenesis of these conditions is yet to be understood. Here, we review the recent findings on ANRIL molecular characterization and function, with a particular focus on its implications in human disease.

Loss of CDKN2B promotes p53-dependent smooth muscle cell apoptosis and aneurysm formation

OBJECTIVE

Genomewide association studies have implicated allelic variation at 9p21.3 in multiple forms of vascular disease, including atherosclerotic coronary heart disease and abdominal aortic aneurysm. As for other genes at 9p21.3, human expression quantitative trait locus studies have associated expression of the tumor suppressor gene CDKN2B with the risk haplotype, but its potential role in vascular pathobiology remains unclear.

METHODS AND RESULTS

Here we used vascular injury models and found that Cdkn2b knockout mice displayed the expected increase in proliferation after injury, but developed reduced neointimal lesions and larger aortic aneurysms. In situ and in vitro studies suggested that these effects were attributable to increased smooth muscle cell apoptosis. Adoptive bone marrow transplant studies confirmed that the observed effects of Cdkn2b were mediated through intrinsic vascular cells and were not dependent on bone marrow-derived inflammatory cells. Mechanistic studies suggested that the observed increase in apoptosis was attributable to a reduction in MDM2 and an increase in p53 signaling, possibly due in part to compensation by other genes at the 9p21.3 locus. Dual inhibition of both Cdkn2b and p53 led to a reversal of the vascular phenotype in each model.

CONCLUSIONS

These results suggest that reduced CDKN2B expression and increased smooth muscle cell apoptosis may be one mechanism underlying the 9p21.3 association with aneurysmal disease.