Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease.

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1-CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease.

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1-CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

Genetic influence on vascular smooth muscle cell apoptosis.

Vascular smooth muscle cell (VSMC) proliferation, migration, and apoptosis play important roles in many physiological processes and pathological conditions. To identify genetic influences on VSMC behavior, we measured these traits and undertook genome-wide association studies in primary umbilical artery-derived VSMCs from >2000 individuals. Although there were no genome-wide significant associations for VSMC proliferation or migration, genetic variants at two genomic loci (7p15.3 and 7q32.3) showed highly significant associations with VSMC apoptosis (P = 1.95 × 10-13 and P = 7.47 × 10-9, respectively). The lead variant at the 7p51.3 locus was associated with increased expression of the GSDME and PALS2 genes in VSMCs. Knockdown of GSDME or PALS2 in VSMCs attenuated apoptotic cell death. A protein co-immunoprecipitation assay indicated that GSDME complexed with PALS2. PALS2 knockdown attenuated activated caspase-3 and GSDME fragmentation, whilst GSDME knockdown also reduced activated caspase-3. These findings provide new insights into the genetic regulation of VSMC apoptosis, with potential utility for therapeutic development.

Genome-wide association meta-analysis of spontaneous coronary artery dissection identifies risk variants and genes related to artery integrity and tissue-mediated coagulation.

Spontaneous coronary artery dissection (SCAD) is an understudied cause of myocardial infarction primarily affecting women. It is not known to what extent SCAD is genetically distinct from other cardiovascular diseases, including atherosclerotic coronary artery disease (CAD). Here we present a genome-wide association meta-analysis (1,917 cases and 9,292 controls) identifying 16 risk loci for SCAD. Integrative functional annotations prioritized genes that are likely to be regulated in vascular smooth muscle cells and artery fibroblasts and implicated in extracellular matrix biology. One locus containing the tissue factor gene F3, which is involved in blood coagulation cascade initiation, appears to be specific for SCAD risk. Several associated variants have diametrically opposite associations with CAD, suggesting that shared biological processes contribute to both diseases, but through different mechanisms. We also infer a causal role for high blood pressure in SCAD. Our findings provide novel pathophysiological insights involving arterial integrity and tissue-mediated coagulation in SCAD and set the stage for future specific therapeutics and preventions.

The FES Gene at the 15q26 Coronary-Artery-Disease Locus Inhibits Atherosclerosis.

BACKGROUND: Genome-wide association studies have discovered a link between genetic variants on human chromosome 15q26.1 and increased coronary artery disease (CAD) susceptibility; however, the underlying pathobiological mechanism is unclear. This genetic locus contains the FES (FES proto-oncogene, tyrosine kinase) gene encoding a cytoplasmic protein-tyrosine kinase involved in the regulation of cell behavior. We investigated the effect of the 15q26.1 variants on FES expression and whether FES plays a role in atherosclerosis. METHODS AND RESULTS: Analyses of isogenic monocytic cell lines generated by CRISPR (clustered regularly interspaced short palindromic repeats)-mediated genome editing showed that monocytes with an engineered 15q26.1 CAD risk genotype had reduced FES expression. Small-interfering-RNA-mediated knockdown of FES promoted migration of monocytes and vascular smooth muscle cells. A phosphoproteomics analysis showed that FES knockdown altered phosphorylation of a number of proteins known to regulate cell migration. Single-cell RNA-sequencing revealed that in human atherosclerotic plaques, cells that expressed FES were predominately monocytes/macrophages, although several other cell types including smooth muscle cells also expressed FES. There was an association between the 15q26.1 CAD risk genotype and greater numbers of monocytes/macrophage in human atherosclerotic plaques. An animal model study demonstrated that Fes knockout increased atherosclerotic plaque size and within-plaque content of monocytes/macrophages and smooth muscle cells, in apolipoprotein E-deficient mice fed a high fat diet. CONCLUSIONS: We provide substantial evidence that the CAD risk variants at the 15q26.1 locus reduce FES expression in monocytes and that FES depletion results in larger atherosclerotic plaques with more monocytes/macrophages and smooth muscle cells. This study is the first demonstration that FES plays a protective role against atherosclerosis and suggests that enhancing FES activity could be a potentially novel therapeutic approach for CAD intervention.

Effects of Coronary Artery Disease-Associated Variants on Vascular Smooth Muscle Cells.

BACKGROUND: Genome-wide association studies have identified many genetic loci that are robustly associated with coronary artery disease (CAD). However, the underlying biological mechanisms are still unknown for most of these loci, hindering the progress to medical translation. Evidence suggests that the genetic influence on CAD susceptibility may act partly through vascular smooth muscle cells (VSMCs). METHODS: We undertook genotyping, RNA sequencing, and cell behavior assays on a large bank of VSMCs (n>1499). Expression quantitative trait locus and splicing quantitative trait locus analyses were performed to identify genes with an expression that was influenced by CAD-associated variants. To identify candidate causal genes for CAD, we ascertained colocalizations of VSMC expression quantitative trait locus signals with CAD association signals by performing causal variants identification in associated regions analysis and the summary data-based mendelian randomization test. Druggability analysis was then performed on the candidate causal genes. CAD risk variants were tested for associations with VSMC proliferation, migration, and apoptosis. Collective effects of multiple CAD-associated variants on VSMC behavior were estimated by polygenic scores. RESULTS: Approximately 60% of the known CAD-associated variants showed statistically significant expression quantitative trait locus or splicing quantitative trait locus effects in VSMCs. Colocalization analyses identified 84 genes with expression quantitative trait locus signals that significantly colocalized with CAD association signals, identifying them as candidate causal genes. Druggability analysis indicated that 38 of the candidate causal genes were druggable, and 13 had evidence of drug-gene interactions. Of the CAD-associated variants tested, 139 showed suggestive associations with VSMC proliferation, migration, or apoptosis. A polygenic score model explained up to 5.94% of variation in several VSMC behavior parameters, consistent with polygenic influences on VSMC behavior. CONCLUSIONS: This comprehensive analysis shows that a large percentage of CAD loci can modulate gene expression in VSMCs and influence VSMC behavior. Several candidate causal genes identified are likely to be druggable and thus represent potential therapeutic targets.

Exploring the Genetic Architecture of Spontaneous Coronary Artery Dissection Using Whole-Genome Sequencing.

BACKGROUND: Spontaneous coronary artery dissection (SCAD) is a cause of acute coronary syndrome that predominantly affects women. Its pathophysiology remains unclear but connective tissue disorders (CTD) and other vasculopathies have been observed in many SCAD patients. A genetic component for SCAD is increasingly appreciated, although few genes have been robustly implicated. We sought to clarify the genetic cause of SCAD using targeted and genome-wide methods in a cohort of sporadic cases to identify both common and rare disease-associated variants. METHODS: A cohort of 91 unrelated sporadic SCAD cases was investigated for rare, deleterious variants in genes associated with either SCAD or CTD, while new candidate genes were sought using rare variant collapsing analysis and identification of novel loss-of-function variants in genes intolerant to such variation. Finally, 2 SCAD polygenic risk scores were applied to assess the contribution of common variants. RESULTS: We identified 10 cases with at least one rare, likely disease-causing variant in CTD-associated genes, although only one had a CTD phenotype. No genes were significantly associated with SCAD from genome-wide collapsing analysis, however, enrichment for TGF (transforming growth factor)-ß signaling pathway genes was found with analysis of 24 genes harboring novel loss-of-function variants. Both polygenic risk scores demonstrated that sporadic SCAD cases have a significantly elevated genetic SCAD risk compared with controls. CONCLUSIONS: SCAD shares some genetic overlap with CTD, even in the absence of any major CTD phenotype. Consistent with a complex genetic architecture, SCAD patients also have a higher burden of common variants than controls.

Contribution of NOTCH1 genetic variants to bicuspid aortic valve and other congenital lesions.

INTRODUCTION: Bicuspid aortic valve (BAV) affects 1% of the general population. NOTCH1 was the first gene associated with BAV. The proportion of familial and sporadic BAV disease attributed to NOTCH1 mutations has not been estimated. AIM: The aim of our study was to provide an estimate of familial and sporadic BAV disease attributable to NOTCH1 mutations. METHODS: The population of our study consisted of participants of the University of Leicester Bicuspid aoRtic vAlVe gEnetic research-8 pedigrees with multiple affected family members and 381 sporadic patients. All subjects underwent NOTCH1 sequencing. A systematic literature search was performed in the NCBI PubMed database to identify publications reporting NOTCH1 sequencing in context of congenital heart disease. RESULTS: NOTCH1 sequencing in 36 subjects from 8 pedigrees identified one variant c.873C>G/p.Tyr291* meeting the American College of Medical Genetics and Genomics criteria for pathogenicity. No pathogenic or likely pathogenic NOTCH1 variants were identified in 381 sporadic patients. Literature review identified 64 relevant publication reporting NOTCH1 sequencing in 528 pedigrees and 9449 sporadic subjects. After excluding families with syndromic disease pathogenic and likely pathogenic NOTCH1 variants were detected in 9/435 (2.1%; 95% CI: 0.7% to 3.4%) of pedigrees and between 0.05% (95% CI: 0.005% to 0.10%) and 0.08% (95% CI: 0.02% to 0.13%) of sporadic patients. Incomplete penetrance of definitely pathogenic NOTCH1 mutations was observed in almost half of reported pedigrees. CONCLUSIONS: Pathogenic and likely pathogenic NOTCH1 genetic variants explain 2% of familial and <0.1% of sporadic BAV disease and are more likely to associate with tetralogy of Fallot and hypoplastic left heart.

Novel LOX Variants in Five Families with Aortic/Arterial Aneurysm and Dissection with Variable Connective Tissue Findings.

Thoracic aortic aneurysm and dissection (TAAD) is a major cause of cardiovascular morbidity and mortality. Loss-of-function variants in LOX, encoding the extracellular matrix crosslinking enzyme lysyl oxidase, have been reported to cause familial TAAD. Using a next-generation TAAD gene panel, we identified five additional probands carrying LOX variants, including two missense variants affecting highly conserved amino acids in the LOX catalytic domain and three truncating variants. Connective tissue manifestations are apparent in a substantial fraction of the variant carriers. Some LOX variant carriers presented with TAAD early in life, while others had normal aortic diameters at an advanced age. Finally, we identified the first patient with spontaneous coronary artery dissection carrying a LOX variant. In conclusion, our data demonstrate that loss-of-function LOX variants cause a spectrum of aortic and arterial aneurysmal disease, often combined with connective tissue findings.

Rare loss-of-function mutations of PTGIR are enriched in fibromuscular dysplasia.

AIMS: Fibromuscular dysplasia (FMD) and spontaneous coronary artery dissection (SCAD) are related, non-atherosclerotic arterial diseases mainly affecting middle-aged women. Little is known about their physiopathological mechanisms. We aimed to identify rare genetic causes to elucidate molecular mechanisms implicated in FMD and SCAD. METHODS AND RESULTS: We analysed 29 exomes that included familial and sporadic FMD. We identified one rare loss-of-function variant (LoF) (frequencygnomAD = 0.000075) shared by two FMD sisters in the prostaglandin I2 receptor gene (PTGIR), a key player in vascular remodelling. Follow-up was conducted by targeted or Sanger sequencing (1071 FMD and 363 SCAD patients) or lookups in exome (264 FMD) or genome sequences (480 SCAD), all independent and unrelated. It revealed four additional LoF allele carriers, in addition to several rare missense variants, among FMD patients, and two LoF allele carriers among SCAD patients, including one carrying a rare splicing mutation (c.768 + 1C>G). We used burden test to test for enrichment in patients compared to gnomAD controls, which detected a putative enrichment in FMD (PTRAPD = 8 × 10-4), but not a significant enrichment (PTRAPD = 0.12) in SCAD. The biological effects of variants on human prostaclycin receptor (hIP) signalling and protein expression were characterized using transient overexpression in human cells. We confirmed the LoFs (Q163X and P17RfsX6) and one missense (L67P), identified in one FMD and one SCAD patient, to severely impair hIP function in vitro. CONCLUSIONS: Our study shows that rare genetic mutations in PTGIR are enriched among FMD patients and found in SCAD patients, suggesting a role for prostacyclin signalling in non-atherosclerotic stenosis and dissection.

Functional investigation of the coronary artery disease gene SVEP1.

A missense variant of the sushi, von Willebrand factor type A, EGF and pentraxin domain containing protein 1 (SVEP1) is genome-wide significantly associated with coronary artery disease. The mechanisms how SVEP1 impacts atherosclerosis are not known. We found endothelial cells (EC) and vascular smooth muscle cells to represent the major cellular source of SVEP1 in plaques. Plaques were larger in atherosclerosis-prone Svep1 haploinsufficient (ApoE-/-Svep1+/-) compared to Svep1 wild-type mice (ApoE-/-Svep1+/+) and ApoE-/-Svep1+/- mice displayed elevated plaque neutrophil, Ly6Chigh monocyte, and macrophage numbers. We assessed how leukocytes accumulated more inside plaques in ApoE-/-Svep1+/- mice and found enhanced leukocyte recruitment from blood into plaques. In vitro, we examined how SVEP1 deficiency promotes leukocyte recruitment and found elevated expression of the leukocyte attractant chemokine (C-X-C motif) ligand 1 (CXCL1) in EC after incubation with missense compared to wild-type SVEP1. Increasing wild-type SVEP1 levels silenced endothelial CXCL1 release. In line, plasma Cxcl1 levels were elevated in ApoE-/-Svep1+/- mice. Our studies reveal an atheroprotective role of SVEP1. Deficiency of wild-type Svep1 increased endothelial CXCL1 expression leading to enhanced recruitment of proinflammatory leukocytes from blood to plaque. Consequently, elevated vascular inflammation resulted in enhanced plaque progression in Svep1 deficiency.

Spontaneous Coronary Artery Dissection: Insights on Rare Genetic Variation From Genome Sequencing.

BACKGROUND: Spontaneous coronary artery dissection (SCAD) occurs when an epicardial coronary artery is narrowed or occluded by an intramural hematoma. SCAD mainly affects women and is associated with pregnancy and systemic arteriopathies, particularly fibromuscular dysplasia. Variants in several genes, such as those causing connective tissue disorders, have been implicated; however, the genetic architecture is poorly understood. Here, we aim to better understand the diagnostic yield of rare variant genetic testing among a cohort of SCAD survivors and to identify genes or gene sets that have a significant enrichment of rare variants. METHODS: We sequenced a cohort of 384 SCAD survivors from the United Kingdom, alongside 13 722 UK Biobank controls and a validation cohort of 92 SCAD survivors. We performed a research diagnostic screen for pathogenic variants and exome-wide and gene-set rare variant collapsing analyses. RESULTS: The majority of patients within both cohorts are female, 29% of the study cohort and 14% validation cohort have a remote arteriopathy. Four cases across the 2 cohorts had a diagnosed connective tissue disorder. We identified pathogenic or likely pathogenic variants in 7 genes (PKD1, COL3A1, SMAD3, TGFB2, LOX, MYLK, and YY1AP1) in 14/384 cases in the study cohort and in 1/92 cases in the validation cohort. In our rare variant collapsing analysis, PKD1 was the highest-ranked gene, and several functionally plausible genes were enriched for rare variants, although no gene achieved study-wide statistical significance. Gene-set enrichment analysis suggested a role for additional genes involved in renal function. CONCLUSIONS: By studying the largest sequenced cohort of SCAD survivors, we demonstrate that, based on current knowledge, only a small proportion have a pathogenic variant that could explain their disease. Our findings strengthen the overlap between SCAD and renal and connective tissue disorders, and we highlight several new genes for future validation.

Novel loss of function mutation in NOTCH1 in a family with bicuspid aortic valve, ventricular septal defect, thoracic aortic aneurysm, and aortic valve stenosis.

BACKGROUND: Bicuspid aortic valve is the most common congenital valvular heart defect in the general population. BAV is associated with significant morbidity due to valve failure, formation of thoracic aortic aneurysm, and increased risk of infective endocarditis and aortic dissection. Loss of function mutations in NOTCH1 (OMIM 190198) has previously been associated with congenital heart disease involving the aortic valve, left ventricle outflow tract, and mitral valve that segregates in affected pedigrees as an autosomal dominant trait with variable expressivity. METHODS: We performed whole-exome sequencing in four members of a three-generational family (three affected and one unaffected subject) with clinical phenotypes including aortic valve stenosis, thoracic aortic aneurysm, and ventricular septal defect. RESULTS: We identified 16 potentially damaging genetic variants (one stop variant, one splice variant, and 14 missense variants) cosegregating with the phenotype. Of these variants, the nonsense mutation (p.Tyr291*) in NOTCH1 was the most deleterious variant identified and the most likely variant causing the disease. CONCLUSION: Inactivating NOTCH1 mutations are a rare cause of familial heart disease involving predominantly left ventricular outflow tract lesions and characterized by the heterogeneity of clinical phenotype.

Effect of a coronary-heart-disease-associated variant of ADAMTS7 on endothelial cell angiogenesis.

BACKGROUND AND AIMS: Recent studies have unveiled an association between ADAMTS7 gene variation and coronary artery disease (CAD) caused by atherosclerosis. We investigated if the ADAMTS7 Serine214-to-Proline substitution arising from a CAD-associated variant affected angiogenesis, since neovascularization plays an important role in atherosclerosis. METHODS AND RESULTS: ADAMTS7 knockdown in vascular endothelial cells (ECs) attenuated their angiogenesis potential, whereas augmented ADAMTS7-Ser214 expression had the opposite effect, leading to increased ECs migratory and tube formation ability. Proteomics analysis showed an increase in thrombospondin-1, a reported angiogenesis inhibitor, in culture media conditioned by ECs with ADAMTS7 knockdown and a decrease of thrombospondin-1 in media conditioned by ECs with ADAMTS7-Ser214 overexpression. Cleavage assay indicated that ADAMTS7 possessed thrombospondin-1 degrading activity, which was reduced by the Ser214-to-Pro substitution. The pro-angiogenic effect of ADAMTS7-Ser214 diminished in the presence of a thrombospondin-1 blocking antibody. CONCLUSIONS: The ADAMTS7 Ser217-to-Pro substitution as a result of ADAMTS7 polymorphism affects thrombospondin-1 degradation, thereby promoting atherogenesis through increased EC migration and tube formation.

HHIPL1, a Gene at the 14q32 Coronary Artery Disease Locus, Positively Regulates Hedgehog Signaling and Promotes Atherosclerosis.

BACKGROUND: Genome-wide association studies have identified chromosome 14q32 as a locus for coronary artery disease. The disease-associated variants fall in a hitherto uncharacterized gene called HHIPL1 (hedgehog interacting protein-like 1), which encodes a sequence homolog of an antagonist of hedgehog signaling. The function of HHIPL1 and its role in atherosclerosis are unknown. METHODS: HHIPL1 cellular localization, interaction with sonic hedgehog (SHH), and influence on hedgehog signaling were tested. HHIPL1 expression was measured in coronary artery disease-relevant human cells, and protein localization was assessed in wild-type and Apoe-/- (apolipoprotein E deficient) mice. Human aortic smooth muscle cell phenotypes and hedgehog signaling were investigated after gene knockdown. Hhipl1-/- mice were generated and aortic smooth muscle cells collected for phenotypic analysis and assessment of hedgehog signaling activity. Hhipl1-/- mice were bred onto both the Apoe-/- and Ldlr-/- (low-density lipoprotein receptor deficient) knockout strains, and the extent of atherosclerosis was quantified after 12 weeks of high-fat diet. Cellular composition and collagen content of aortic plaques were assessed by immunohistochemistry. RESULTS: In vitro analyses revealed that HHIPL1 is a secreted protein that interacts with SHH and increases hedgehog signaling activity. HHIPL1 expression was detected in human smooth muscle cells and in smooth muscle within atherosclerotic plaques of Apoe-/- mice. The expression of Hhipl1 increased with disease progression in aortic roots of Apoe-/- mice. Proliferation and migration were reduced in Hhipl1 knockout mouse and HHIPL1 knockdown aortic smooth muscle cells, and hedgehog signaling was decreased in HHIPL1-deficient cells. Hhipl1 knockout caused a reduction of >50% in atherosclerosis burden on both Apoe-/- and Ldlr-/- knockout backgrounds, and lesions were characterized by reduced smooth muscle cell content. CONCLUSIONS: HHIPL1 is a secreted proatherogenic protein that enhances hedgehog signaling and regulates smooth muscle cell proliferation and migration. Inhibition of HHIPL1 protein function might offer a novel therapeutic strategy for coronary artery disease.

Long noncoding RNA NEXN-AS1 mitigates atherosclerosis by regulating the actin-binding protein NEXN.

Noncoding RNAs are emerging as important players in gene regulation and disease pathogeneses. Here, we show that a previously uncharacterized long noncoding RNA, nexilin F-actin binding protein antisense RNA 1 (NEXN-AS1), modulates the expression of the actin-binding protein NEXN and that NEXN exerts a protective role against atherosclerosis. An expression microarray analysis showed that the expression of both NEXN-AS1 and NEXN was reduced in human atherosclerotic plaques. In vitro experiments revealed that NEXN-AS1 interacted with the chromatin remodeler BAZ1A and the 5' flanking region of the NEXN gene and that it also upregulated NEXN expression. Augmentation of NEXN-AS1 expression inhibited TLR4 oligomerization and NF-κB activity, downregulated the expression of adhesion molecules and inflammatory cytokines by endothelial cells, and suppressed monocyte adhesion to endothelial cells. These inhibitory effects of NEXN-AS1 were abolished by knockdown of NEXN. In vivo experiments using ApoE-knockout mice fed a Western high-fat diet demonstrated that NEXN deficiency promoted atherosclerosis and increased macrophage abundance in atherosclerotic lesions, with heightened expression of adhesion molecules and inflammatory cytokines, whereas augmented NEXN expression deterred atherosclerosis. Patients with coronary artery disease were found to have lower blood NEXN levels than healthy individuals. These results indicate that NEXN-AS1 and NEXN represent potential therapeutic targets in atherosclerosis-related diseases.

Influence of a Coronary Artery Disease-Associated Genetic Variant on FURIN Expression and Effect of Furin on Macrophage Behavior.

Objective- Genome-wide association studies have revealed a robust association between genetic variation on chromosome 15q26.1 and coronary artery disease (CAD) susceptibility; however, the underlying biological mechanism is still unknown. The lead CAD-associated genetic variant (rs17514846) at this locus resides in the FURIN gene. In advanced atherosclerotic plaques, furin is expressed primarily in macrophages. We investigated whether this CAD-associated variant alters FURIN expression and whether furin affects monocyte/macrophage behavior. Approach and Results- A quantitative reverse transcription polymerase chain reaction analysis showed that leukocytes from individuals carrying the CAD risk allele (A) of rs17514846 had increased FURIN expression. A chromatin immunoprecipitation assay revealed higher RNA polymerase II occupancy in the FURIN gene in mononuclear cells of individuals carrying this allele. A reporter gene assay in transiently transfected monocytes/macrophages indicated that the CAD risk allele had higher transcriptional activity than the nonrisk allele (C). An analysis of isogenic monocyte cell lines created by CRISPR (clustered regularly interspaced short palindromic repeats)-mediated genome editing showed that isogenic cells with the A/A genotype for rs17514846 had higher FURIN expression levels than the isogenic cells with the C/C genotype. An electrophoretic mobility shift assay exhibited preferential binding of a nuclear protein to the risk allele. Studies of monocytes/macrophages with lentivirus-mediated furin overexpression or shRNA (short hairpin RNA)-induced furin knockdown showed that furin overexpression promoted monocyte/macrophage migration, increased proliferation, and reduced apoptosis whereas furin knockdown had the opposite effects. Conclusions- Our study shows that the CAD-associated genetic variant increases FURIN expression and that furin promotes monocyte/macrophage migration and proliferation while inhibiting apoptosis, providing a biological mechanism for the association between variation at the chromosome 15q26.1 locus and CAD risk.

JCAD, a Gene at the 10p11 Coronary Artery Disease Locus, Regulates Hippo Signaling in Endothelial Cells.

Objective- A large number of genetic loci have been associated with risk of coronary artery disease (CAD) through genome-wide association studies, however, for most loci the underlying biological mechanism is unknown. Determining the molecular pathways and cellular processes affected by these loci will provide new insights into CAD pathophysiology and may lead to new therapies. The CAD-associated variants at 10p11.23 fall in JCAD, which encodes an endothelial junction protein, however, its molecular function in endothelial cells is not known. In this study, we characterize the molecular role of JCAD (junctional cadherin 5 associated) in endothelial cells. Approach and Results- We show that JCAD knockdown in endothelial cells affects key phenotypes related to atherosclerosis including proliferation, migration, apoptosis, tube formation, and monocyte binding. We demonstrate that JCAD interacts with LATS2 (large tumor suppressor kinase 2) and negatively regulates Hippo signaling leading to increased activity of YAP (yes-associated protein), the transcriptional effector of the pathway. We also show by double siRNA knockdown that the phenotypes caused by JCAD knockdown require LATS2 and that JCAD is involved in transmission of RhoA-mediated signals into the Hippo pathway. In human tissues, we find that the CAD-associated lead variant, rs2487928, is associated with expression of JCAD in arteries, including atherosclerotic arteries. Gene co-expression analyses across disease-relevant tissues corroborate our phenotypic findings and support the link between JCAD and Hippo signaling. Conclusions- Our results show that JCAD negatively regulates Hippo signaling in endothelial cells and we suggest that JCAD contributes to atherosclerosis by mediating YAP activity and contributing to endothelial dysfunction.

Network analysis of coronary artery disease risk genes elucidates disease mechanisms and druggable targets.

Genome-wide association studies (GWAS) have identified over two hundred chromosomal loci that modulate risk of coronary artery disease (CAD). The genes affected by variants at these loci are largely unknown and an untapped resource to improve our understanding of CAD pathophysiology and identify potential therapeutic targets. Here, we prioritized 68 genes as the most likely causal genes at genome-wide significant loci identified by GWAS of CAD and examined their regulatory roles in 286 metabolic and vascular tissue gene-protein sub-networks ("modules"). The modules and genes within were scored for CAD druggability potential. The scoring enriched for targets of cardiometabolic drugs currently in clinical use and in-depth analysis of the top-scoring modules validated established and revealed novel target tissues, biological processes, and druggable targets. This study provides an unprecedented resource of tissue-defined gene-protein interactions directly affected by genetic variance in CAD risk loci.