Identification of CAD candidate genes in GWAS loci and their expression in vascular cells

Secreted protein profiling of human aortic smooth muscle cells identifies vascular disease associations

Background

Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular diseases (CVD), the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a "contractile" to a "synthetic" phenotype characterized by an increased proliferation, migration, production of extracellular matrix (ECM) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of CVD, including coronary artery disease (CAD), stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies (GWAS).

Methods

Using human aortic SMCs from 123 multi-ancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomic analysis of the conditioned media.

Results

We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping (pQTL) and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions.

Conclusions

Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Lipid phosphate phosphatase 3 maintains NO-mediated flow-mediated dilatation in human adipose resistance arterioles

Microvascular dysfunction predicts adverse cardiovascular events despite absence of large vessel disease. A shift in the mediator of flow-mediated dilatation (FMD) from nitric oxide (NO) to mitochondrial-derived hydrogen peroxide (H2 O2 ) occurs in arterioles from patients with coronary artery disease (CAD). The underlying mechanisms governing this shift are not completely defined. Lipid phosphate phosphatase 3 (LPP3) is a transmembrane protein that dephosphorylates lysophosphatidic acid, a bioactive lipid, causing a receptor-mediated increase in reactive oxygen species. A single nucleotide loss-of-function polymorphism in the gene coding for LPP3 (rs17114036) is associated with elevated risk for CAD, independent of traditional risk factors. LPP3 is suppressed by miR-92a, which is elevated in the circulation of patients with CAD. Repression of LPP3 increases vascular inflammation and atherosclerosis in animal models. We investigated the role of LPP3 and miR-92a as a mechanism for microvascular dysfunction in CAD. We hypothesized that modulation of LPP3 is critically involved in the disease-associated shift in mediator of FMD. LPP3 protein expression was reduced in left ventricle tissue from CAD relative to non-CAD patients (P = 0.004), with mRNA expression unchanged (P = 0.96). Reducing LPP3 expression (non-CAD) caused a shift from NO to H2 O2 (% maximal dilatation: Control 78.1 ± 11.4% vs. Peg-Cat 30.0 ± 11.2%; P < 0.0001). miR-92a is elevated in CAD arterioles (fold change: 1.9 ± 0.01 P = 0.04), while inhibition of miR-92a restored NO-mediated FMD (CAD), and enhancing miR-92a expression (non-CAD) elicited H2 O2 -mediated dilatation (P < 0.0001). Our data suggests LPP3 is crucial in the disease-associated switch in the mediator of FMD. KEY POINTS: Lipid phosphate phosphatase 3 (LPP3) expression is reduced in heart tissue patients with coronary artery disease (CAD). Loss of LPP3 in CAD is associated with an increase in the LPP3 inhibitor, miR-92a. Inhibition of LPP3 in the microvasculature of healthy patients mimics the CAD flow-mediated dilatation (FMD) phenotype. Inhibition of miR-92a restores nitric oxide-mediated FMD in the microvasculature of CAD patients.

A construction and comprehensive analysis of the immune-related core ceRNA network and infiltrating immune cells in peripheral arterial occlusive disease

Background: Peripheral arterial occlusive disease (PAOD) is a peripheral artery disorder that increases with age and often leads to an elevated risk of cardiovascular events. The purposes of this study were to explore the underlying competing endogenous RNA (ceRNA)-related mechanism of PAOD and identify the corresponding immune cell infiltration patterns.

Methods: An available gene expression profile (GSE57691 datasets) was downloaded from the GEO database. Differentially expressed (DE) mRNAs and lncRNAs were screened between 9 PAOD and 10 control samples. Then, the lncRNA-miRNA-mRNA ceRNA network was constructed on the basis of the interactions generated from the miRcode, TargetScan, miRDB, and miRTarBase databases. The functional enrichment and protein–protein interaction analyses of mRNAs in the ceRNA network were performed. Immune-related core mRNAs were screened out through the Venn method. The compositional patterns of the 22 types of immune cell fraction in PAOD were estimated through the CIBERSORT algorithm. The final ceRNA network and immune infiltration were validated using clinical tissue samples. Finally, the correlation between immune cells and mRNAs in the final ceRNA network was analyzed.

Results: Totally, 67 DE_lncRNAs and 1197 DE_mRNAs were identified, of which 130 DE_mRNAs (91 downregulated and 39 upregulated) were lncRNA-related. The gene ontology enrichment analysis showed that those down- and upregulated genes were involved in dephosphorylation and regulation of translation, respectively. The final immune-related core ceRNA network included one lncRNA (LINC00221), two miRNAs (miR-17-5p and miR-20b-5p), and one mRNA (CREB1). Meanwhile, we found that monocytes and M1 macrophages were the main immune cell subpopulations in PAOD. After verification, these predictions were consistent with experimental results. Moreover, CREB1 was positively correlated with naive B cells (R = 0.55, p = 0.035) and monocytes (R = 0.52, p = 0.049) and negatively correlated with M1 macrophages (R = −0.72, p = 0.004), resting mast cells (R = −0.66, p = 0.009), memory B cells (R = −0.55, p = 0.035), and plasma cells (R = −0.52, p = 0.047).

Conclusion: In general, we proposed that the immune-related core ceRNA network (LINC00221, miR-17-5p, miR-20b-5p, and CREB1) and infiltrating immune cells (monocytes and M1 macrophages) could help further explore the molecular mechanisms of PAOD.

Correlation between platelet-derived growth factor-B gene polymorphism and coronary heart disease

OBJECT

The aim of the present work was to investigate the correlation of plasma platelet-derived growth factor (PDGF)-BB level and single nucleotide polymorphism (SNP, rs1800817 and rs2285094) of PDGF-B gene with the onset and stability condition of coronary heart disease (CHD).

METHODS

Totally, 335 subjects were included in and divided into CHD (n = 247) and control group (n = 88) according to coronary angiography. Besides, the patients in the CHD group were divided into acute coronary syndrome (ACS) group (n = 165) and stable angina pectoria (SAP) group (n = 82), based on CHD stability condition. The plasma PDGF-BB level was measured by ELISA, and the genotype of PDGF-B was examined through qPCR assay.

RESULTS

The PDGF-BB level was positively correlated with hsCRP level (r = 0.149, p < 0.05). The genotype frequencies of SNP rs1800817 and rs2285094 match Hardy-Weinberg equilibrium. There was weak linkage disequilibrium between SNP rs1800817 and rs2285094: D' = 0.419, r²  = 0.04, which has no correlation with CHD. There was no statistical difference in plasma PDGF-BB level among different genotypes in rs1800817 and rs2285094. There were no differences in the plasma PDGF-BB level among patients with any genotype of SNP rs1800817 and rs2285094, no matter how it was grouped. Logistic regression results indicated that the plasma PDGF-BB level was the independent risk factor of CHD onset (OR = 1.003, 95% CI 1.001-1.006, p = 0.014).

CONCLUSIONS

High plasma PDGF-BB level is the risk factor of CHD and has correlation with instability of CHD. The plasma PDGF-BB level change may be related to inflammatory response. PDGF-B gene rs1800817 and rs2285094 polymorphisms are not correlated with CHD.

GALNT4 primes monocytes adhesion and transmigration by regulating O-Glycosylation of PSGL-1 in atherosclerosis

Atherosclerosis is a major underlying cause of cardiovascular disease. Genome wide association studies have predicted that GalNAc-T4 (GALNT4), which responsible for initiating step of mucin-type O-glycosylation, plays a causal role in the susceptibility to cardiovascular diseases, whereas the precise mechanism remains obscure. Thus, we sought to determine the role and mechanism of GALNT4 in atherosclerosis. Firstly, we found the expression of GALNT4 and protein O-glycosylation were both increased in plaque as atherosclerosis progressed in ApoE-/- mice by immunohistochemistry. And the expression of GALNT4 was also increased in human monocytes treated with ACS (acute coronary syndrome) sera and subjected to LPS and ox-LDL in vitro. Moreover, silencing expression of GALNT4 by shRNA lentivirus alleviated atherosclerotic plaque formation and monocyte/macrophage infiltration in ApoE-/- mice. Functional investigations demonstrate that GALNT4 knockdown inhibited P-selectin-induced activation of β2 integrin on the surface of monocytes, decreased monocytes adhesion under flow condition with P-selectin stimulation, as well as suppressed monocytes transmigration triggered by monocyte chemotactic protein- 1(MCP-1). In contrast, GALNT4 overexpression enhanced monocytes adhesion and transmigration. Furthermore, Vicia Villosa Lectin (VVL) pull down and PSGL-1 immunoprecipitation assays showed that GALNT4 overexpression increased O-Glycosylation of PSGL-1 and P-selectin induce phosphorylation of Akt/mTOR and IκBα/NFκB on monocytes. Conversely, knockdown of GALNT4 decreased VVL binding and attenuated the activation of Akt/mTOR and IκBα/NFκB. Additionally, mTOR inhibitor rapamycin blocked these effects of GALNT4 overexpression on monocytes. Collectively, GALNT4 catalyzed PSGL-1 O-glycosylation that involved in P-selectin induced monocytes adhesion and transmigration via Akt/mTOR and NFκB pathway. Thus, GALNT4 may be a potential therapeutic target for atherosclerosis.

Bioactive lipids and metabolic syndrome-a symposium report

Recent research has shed light on the cellular and molecular functions of bioactive lipids that go far beyond what was known about their role as dietary lipids. Bioactive lipids regulate inflammation and its resolution as signaling molecules. Genetic studies have identified key factors that can increase the risk of cardiovascular diseases and metabolic syndrome through their effects on lipogenesis. Lipid scientists have explored how these signaling pathways affect lipid metabolism in the liver, adipose tissue, and macrophages by utilizing a variety of techniques in both humans and animal models, including novel lipidomics approaches and molecular dynamics models. Dissecting out these lipid pathways can help identify mechanisms that can be targeted to prevent or treat cardiometabolic conditions. Continued investigation of the multitude of functions mediated by bioactive lipids may reveal additional components of these pathways that can provide a greater understanding of metabolic homeostasis.

ZEB2 Shapes the Epigenetic Landscape of Atherosclerosis

BACKGROUND

Smooth muscle cells (SMCs) transition into a number of different phenotypes during atherosclerosis, including those that resemble fibroblasts and chondrocytes, and make up the majority of cells in the atherosclerotic plaque. To better understand the epigenetic and transcriptional mechanisms that mediate these cell state changes, and how they relate to risk for coronary artery disease (CAD), we have investigated the causality and function of transcription factors at genome-wide associated loci.

METHODS

We used CRISPR-Cas 9 genome and epigenome editing to identify the causal gene and cells for a complex CAD genome-wide association study signal at 2q22.3. Single-cell epigenetic and transcriptomic profiling in murine models and human coronary artery smooth muscle cells were used to understand the cellular and molecular mechanism by which this CAD risk gene exerts its function.

RESULTS

CRISPR-Cas 9 genome and epigenome editing showed that the complex CAD genetic signals within a genomic region at 2q22.3 lie within smooth muscle long-distance enhancers for ZEB2, a transcription factor extensively studied in the context of epithelial mesenchymal transition in development of cancer. Zeb2 regulates SMC phenotypic transition through chromatin remodeling that obviates accessibility and disrupts both Notch and transforming growth factor β signaling, thus altering the epigenetic trajectory of SMC transitions. SMC-specific loss of Zeb2 resulted in an inability of transitioning SMCs to turn off contractile programing and take on a fibroblast-like phenotype, but accelerated the formation of chondromyocytes, mirroring features of high-risk atherosclerotic plaques in human coronary arteries.

CONCLUSIONS

These studies identify ZEB2 as a new CAD genome-wide association study gene that affects features of plaque vulnerability through direct effects on the epigenome, providing a new therapeutic approach to target vascular disease.

Large-Scale Whole Genome Sequencing Study Reveals Genetic Architecture and Key Variants for Breast Muscle Weight in Native Chickens

Breast muscle weight (BrW) is one of the most important economic traits in chicken, and directional breeding for that results in both phenotypic and genetic changes. The Jingxing yellow chicken, including an original (without human-driven selection) line and a selected line (based on selection for increased intramuscular fat content), were used to dissect the genetic architecture and key variants associated with BrW. We detected 1069 high-impact single nucleotide polymorphisms (SNPs) with high conserved score and significant frequency difference between two lines. Based on the annotation result, the ECM-receptor interaction and fatty acid biosynthesis were enriched, and muscle-related genes, including MYOD1, were detected. By performing genome-wide association study for the BrW trait, we defined a major haplotype and two conserved SNPs that affected BrW. By integrated genomic and transcriptomic analysis, IGF2BP1 was identified as the crucial gene associated with BrW. In conclusion, these results offer a new insight into chicken directional selection and provide target genetic markers by which to improve chicken BrW.

HMOX1 upregulation promotes ferroptosis in diabetic atherosclerosis

OBJECTIVE

Atherosclerotic vascular disease remains the principal cause of death and disability among patients with type 2 diabetes. Unfortunately, the problem is not adequately resolved by therapeutic strategies with currently available drugs or approaches that solely focus on optimal glycemic control. To identify the key contributors and better understand the mechanism of diabetic atherosclerotic vascular disease, we aimed to elucidate the key genetic characteristics and pathological pathways in atherosclerotic vascular disease through nonbiased bioinformatics analysis and subsequent experimental demonstration and exploration in diabetic atherosclerotic vascular disease.

METHODS AND RESULTS

Sixty-eight upregulated and 23 downregulated genes were identified from the analysis of gene expression profiles (GSE30169 and GSE6584). A comprehensive bioinformatic assay further identified that ferroptosis, a new type of programmed cell death and HMOX1 (a gene that encodes heme oxygenase), were vital factors in atherosclerotic vascular disease. We further demonstrated that diabetes significantly increased ferroptosis and HMOX1 levels compared to normal controls. Importantly, the ferroptosis inhibitor ferrostatin-1 (Fer-1) effectively attenuated diabetic atherosclerosis, suggesting the causative role of ferroptosis in diabetic atherosclerosis development. At the cellular level, Fer-1 ameliorated high glucose high lipid-induced lipid peroxidation and downregulated ROS production. More importantly, HMOX1 knockdown attenuated Fe2+ overload, reduced iron content and ROS, and alleviated lipid peroxidation, which led to a reduction in ferroptosis in diabetic human endothelial cells.

CONCLUSIONS

We demonstrated that HMOX1 upregulation is responsible for the increased ferroptosis in diabetic atherosclerosis development, suggesting that HMOX1 may serve as a potential therapeutic or drug development target for diabetic atherosclerosis.

Multi-omics approaches for revealing the complexity of cardiovascular disease

The development and progression of cardiovascular disease (CVD) can mainly be attributed to the narrowing of blood vessels caused by atherosclerosis and thrombosis, which induces organ damage that will result in end-organ dysfunction characterized by events such as myocardial infarction or stroke. It is also essential to consider other contributory factors to CVD, including cardiac remodelling caused by cardiomyopathies and co-morbidities with other diseases such as chronic kidney disease. Besides, there is a growing amount of evidence linking the gut microbiota to CVD through several metabolic pathways. Hence, it is of utmost importance to decipher the underlying molecular mechanisms associated with these disease states to elucidate the development and progression of CVD. A wide array of systems biology approaches incorporating multi-omics data have emerged as an invaluable tool in establishing alterations in specific cell types and identifying modifications in signalling events that promote disease development. Here, we review recent studies that apply multi-omics approaches to further understand the underlying causes of CVD and provide possible treatment strategies by identifying novel drug targets and biomarkers. We also discuss very recent advances in gut microbiota research with an emphasis on how diet and microbial composition can impact the development of CVD. Finally, we present various biological network analyses and other independent studies that have been employed for providing mechanistic explanation and developing treatment strategies for end-stage CVD, namely myocardial infarction and stroke.

An artificial neural network approach integrating plasma proteomics and genetic data identifies PLXNA4 as a new susceptibility locus for pulmonary embolism

Venous thromboembolism is the third common cardiovascular disease and is composed of two entities, deep vein thrombosis (DVT) and its potential fatal form, pulmonary embolism (PE). While PE is observed in ~ 40% of patients with documented DVT, there is limited biomarkers that can help identifying patients at high PE risk. To fill this need, we implemented a two hidden-layers artificial neural networks (ANN) on 376 antibodies and 19 biological traits measured in the plasma of 1388 DVT patients, with or without PE, of the MARTHA study. We used the LIME algorithm to obtain a linear approximate of the resulting ANN prediction model. As MARTHA patients were typed for genotyping DNA arrays, a genome wide association study (GWAS) was conducted on the LIME estimate. Detected single nucleotide polymorphisms (SNPs) were tested for association with PE risk in MARTHA. Main findings were replicated in the EOVT study composed of 143 PE patients and 196 DVT only patients. The derived ANN model for PE achieved an accuracy of 0.89 and 0.79 in our training and testing sets, respectively. A GWAS on the LIME approximate identified a strong statistical association peak (rs1424597: p = 5.3 × 10^-7) at the PLXNA4 locus. Homozygote carriers for the rs1424597-A allele were then more frequently observed in PE than in DVT patients from the MARTHA (2% vs. 0.4%, p = 0.005) and the EOVT (3% vs. 0%, p = 0.013) studies. In a sample of 112 COVID-19 patients known to have endotheliopathy leading to acute lung injury and an increased risk of PE, decreased PLXNA4 levels were associated (p = 0.025) with worsened respiratory function. Using an original integrated proteomics and genetics strategy, we identified PLXNA4 as a new susceptibility gene for PE whose exact role now needs to be further elucidated.

Shared genetic etiology and causality between body fat percentage and cardiovascular diseases: a large-scale genome-wide cross-trait analysis

BACKGROUND

Accumulating evidences have suggested that high body fat percentage (BF%) often occurs in parallel with cardiovascular diseases (CVDs), implying a common etiology between them. However, the shared genetic etiology underlying BF% and CVDs remains unclear.

METHODS

Using large-scale genome-wide association study (GWAS) data, we investigated shared genetics between BF% (N = 100,716) and 10 CVD-related traits (n = 6968-977,323) with linkage disequilibrium score regression, multi-trait analysis of GWAS, and transcriptome-wide association analysis, and evaluated causal associations using Mendelian randomization.

RESULTS

We found strong positive genetic correlations between BF% and heart failure (HF) (Rg = 0.47, P = 1.27 × 10^- 22) and coronary artery disease (CAD) (Rg = 0.22, P = 3.26 × 10^- 07). We identified 5 loci and 32 gene-tissue pairs shared between BF% and HF, as well as 16 loci and 28 gene-tissue pairs shared between BF% and CAD. The loci were enriched in blood vessels and brain tissues, while the gene-tissue pairs were enriched in the nervous, cardiovascular, and exo-/endocrine system. In addition, we observed that BF% was causally related with a higher risk of HF (odds ratio 1.63 per 1-SD increase in BF%, P = 4.16 × 10-04) using a MR approach.

CONCLUSIONS

Our findings suggest that BF% and CVDs have shared genetic etiology and targeted reduction of BF% may improve cardiovascular outcomes. This work advances our understanding of the genetic basis underlying co-morbid obesity and CVDs and opens up a new way for early prevention of CVDs.

Sex-Stratified Gene Regulatory Networks Reveal Female Key Driver Genes of Atherosclerosis Involved in Smooth Muscle Cell Phenotype Switching

BACKGROUND

Although sex differences in coronary artery disease are widely accepted with women developing more stable atherosclerosis than men, the underlying pathobiology of such differences remains largely unknown. In coronary artery disease, recent integrative systems biological studies have inferred gene regulatory networks (GRNs). Within these GRNs, key driver genes have shown great promise but have thus far been unidentified in women.

METHODS

We generated sex-specific GRNs of the atherosclerotic arterial wall in 160 women and age-matched men in the STARNET study (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task). We integrated the female GRNs with single-cell RNA-sequencing data of the human atherosclerotic plaque and single-cell RNA sequencing of advanced atherosclerotic lesions from wild type and Klf4 knockout atherosclerotic smooth muscle cell (SMC) lineage-tracing mice.

RESULTS

By comparing sex-specific GRNs, we observed clear sex differences in network activity within the atherosclerotic tissues. Genes more active in women were associated with mesenchymal cells and endothelial cells, whereas genes more active in men were associated with the immune system. We determined that key drivers of GRNs active in female coronary artery disease were predominantly found in (SMCs by single-cell sequencing of the human atherosclerotic plaques, and higher expressed in female plaque SMCs, as well. To study the functions of these female SMC key drivers in atherosclerosis, we examined single-cell RNA sequencing of advanced atherosclerotic lesions from wild type and Klf4 knockout atherosclerotic SMC lineage-tracing mice. The female key drivers were found to be expressed by phenotypically modulated SMCs and affected by Klf4, suggesting that sex differences in atherosclerosis involve phenotypic switching of plaque SMCs.

CONCLUSIONS

Our systems approach provides novel insights into molecular mechanisms that underlie sex differences in atherosclerosis. To discover sex-specific therapeutic targets for atherosclerosis, an increased emphasis on sex-stratified approaches in the analysis of multi-omics data sets is warranted.

Association of PHACTR1 intronic variants with the first myocardial infarction and their effect on PHACTR1 mRNA expression in PBMCs

BACKGROUND

Myocardial infarction (MI) and underlining atherosclerosis are the main causes of death worldwide. Phosphatase and actin regulator 1 (PHACTR1) variants have been associated with early onset MI, coronary artery disease and carotid dissection. PHACTR1 mRNA expression has been detected in tissues and cells related to atherosclerosis. Nonetheless, the true effect of PHACTR1 on vascular diseases is still unknown. Our aim was to examine the association of PHACTR1 intronic variants, rs9349379, rs2026458 and rs2876300, with MI and multi-vessel disease (MVD) and to assess their effect on PHACTR1 and EDN1 mRNA expression in PBMCs of patients six months after MI.

METHODS

The study enrolled 537 patients with the first MI and 310 controls. Gene expression was assessed in 74 patients six months after MI and 37 healthy controls. Rs9349379, rs2026458, rs2876300 and relative mRNA expressions were detected by TaqMan® technology.

RESULTS

The significant association between PHACTR1 variants and MI was not found, either individually or in haplotype. A higher frequency of rs2876300G-allele in MVD was rendered not significant after Bonferroni correction. PHACTR1 mRNA was significantly increased in PBMCs of patients six months after MI compared to controls (p = 0.02). Patients that carry ACG haplotype have increased PHACTR1 mRNA expression in PBMCs (p = 0.04). There was no effect of PHACTR1 variants on EDN1 mRNA expression.

CONCLUSION

Our findings suggest that PHACTR1 intronic variants may have a role in severity and progression of coronary atherosclerosis. Future research is needed to clarify the mechanism underlying the role of PHACTR1 in coronary atherosclerosis and MI.

Genetic Regulation of Atherosclerosis-Relevant Phenotypes in Human Vascular Smooth Muscle Cells

RATIONALE

Coronary artery disease (CAD) is a major cause of morbidity and mortality worldwide. Recent genome-wide association studies revealed 163 loci associated with CAD. However, the precise molecular mechanisms by which the majority of these loci increase CAD risk are not known. Vascular smooth muscle cells (VSMCs) are critical in the development of CAD. They can play either beneficial or detrimental roles in lesion pathogenesis, depending on the nature of their phenotypic changes.

OBJECTIVE

To identify genetic variants associated with atherosclerosis-relevant phenotypes in VSMCs.

METHODS AND RESULTS

We quantified 12 atherosclerosis-relevant phenotypes related to calcification, proliferation, and migration in VSMCs isolated from 151 multiethnic heart transplant donors. After genotyping and imputation, we performed association mapping using 6.3 million genetic variants. We demonstrated significant variations in calcification, proliferation, and migration. These phenotypes were not correlated with each other. We performed genome-wide association studies for 12 atherosclerosis-relevant phenotypes and identified 4 genome-wide significant loci associated with at least one VSMC phenotype. We overlapped the previously identified CAD loci with our data set and found nominally significant associations at 79 loci. One of them was the chromosome 1q41 locus, which harbors MIA3. The G allele of the lead risk single nucleotide polymorphism (SNP) rs67180937 was associated with lower VSMC MIA3 expression and lower proliferation. Lentivirus-mediated silencing of MIA3 (melanoma inhibitory activity protein 3) in VSMCs resulted in lower proliferation, consistent with human genetics findings. Furthermore, we observed a significant reduction of MIA3 protein in VSMCs in thin fibrous caps of late-stage atherosclerotic plaques compared to early fibroatheroma with thick and protective fibrous caps in mice and humans.

CONCLUSIONS

Our data demonstrate that genetic variants have significant influences on VSMC function relevant to the development of atherosclerosis. Furthermore, high MIA3 expression may promote atheroprotective VSMC phenotypic transitions, including increased proliferation, which is essential in the formation or maintenance of a protective fibrous cap.

A Platelet Function Modulator of Thrombin Activation Is Causally Linked to Cardiovascular Disease and Affects PAR4 Receptor Signaling

Dual antiplatelet therapy reduces ischemic events in cardiovascular disease, but it increases bleeding risk. Thrombin receptors PAR1 and PAR4 are drug targets, but the role of thrombin in platelet aggregation remains largely unexplored in large populations. We performed a genome-wide association study (GWAS) of platelet aggregation in response to full-length thrombin, followed by clinical association analyses, Mendelian randomization, and functional characterization including iPSC-derived megakaryocyte and platelet experiments. We identified a single sentinel variant in the GRK5 locus (rs10886430-G, p = 3.0 × 10^-42) associated with increased thrombin-induced platelet aggregation (β = 0.70, SE = 0.05). We show that disruption of platelet GRK5 expression by rs10886430-G is associated with enhanced platelet reactivity. The proposed mechanism of a GATA1-driven megakaryocyte enhancer is confirmed in allele-specific experiments. Utilizing further data, we demonstrate that the allelic effect is highly platelet- and thrombin-specific and not likely due to effects on thrombin levels. The variant is associated with increased risk of cardiovascular disease outcomes in UK BioBank, most strongly with pulmonary embolism. The variant associates with increased risk of stroke in the MEGASTROKE, UK BioBank, and FinnGen studies. Mendelian randomization analyses in independent samples support a causal role for rs10886430-G in increasing risk for stroke, pulmonary embolism, and venous thromboembolism through its effect on thrombin-induced platelet reactivity. We demonstrate that G protein-coupled receptor kinase 5 (GRK5) promotes platelet activation specifically via PAR4 receptor signaling. GRK5 inhibitors in development for the treatment of heart failure and cancer could have platelet off-target deleterious effects. Common variants in GRK5 may modify clinical outcomes with PAR4 inhibitors, and upregulation of GRK5 activity or signaling in platelets may have therapeutic benefits.

Intrinsic transcriptomic sex differences in human endothelial cells at birth and in adults are associated with coronary artery disease targets

Sex differences in endothelial cell (EC) biology may reflect intrinsic differences driven by chromosomes or sex steroid exposure and gender differences accumulated over life. We analysed EC gene expression data from boy-girl twins at birth and in non-twin adults to detect sex differences at different stages of life, and show that 14-25% of the EC transcriptome is sex-biased. By combining data from both stages of life, we identified sex differences that are present at birth and maintained throughout life, and those that are acquired over life. Promisingly, we found that genes that present with an acquired sex difference in ECs are more likely to be targets of sex steroids. Annotating both gene sets with data from multiple genome-wide association studies (GWAS) revealed that genes with an intrinsic sex difference in ECs are enriched for coronary artery disease GWAS hits. This study underscores the need for treating sex as a biological variable.

Detecting rare haplotypes associated with complex diseases using both population and family data: Combined logistic Bayesian Lasso

Haplotype-based association methods have been developed to understand the genetic architecture of complex diseases. Compared to single-variant-based methods, haplotype methods are thought to be more biologically relevant, since there are typically multiple non-independent genetic variants involved in complex diseases, and the use of haplotypes implicitly accounts for non-independence caused by linkage disequilibrium. In recent years, with the focus moving from common to rare variants, haplotype-based methods have also evolved accordingly to uncover the roles of rare haplotypes. One particular approach is regularization-based, with the use of Bayesian least absolute shrinkage and selection operator (Lasso) as an example. This type of methods has been developed for either case-control population data (the logistic Bayesian Lasso (LBL)) or family data (family-triad-based logistic Bayesian Lasso (famLBL)). In some situations, both family data and case-control data are available; therefore, it would be a waste of resources if only one of them could be analyzed. To make full usage of available data to increase power, we propose a unified approach that can combine both case-control and family data (combined logistic Bayesian Lasso (cLBL)). Through simulations, we characterized the performance of cLBL and showed the advantage of cLBL over existing methods. We further applied cLBL to the Framingham Heart Study data to demonstrate its utility in real data applications.

Genetic Insights Into Smooth Muscle Cell Contributions to Coronary Artery Disease

Coronary artery disease is a complex cardiovascular disease involving an interplay of genetic and environmental influences over a lifetime. Although considerable progress has been made in understanding lifestyle risk factors, genetic factors identified from genome-wide association studies may capture additional hidden risk undetected by traditional clinical tests. These genetic discoveries have highlighted many candidate genes and pathways dysregulated in the vessel wall, including those involving smooth muscle cell phenotypic modulation and injury responses. Here, we summarize experimental evidence for a few genome-wide significant loci supporting their roles in smooth muscle cell biology and disease. We also discuss molecular quantitative trait locus mapping as a powerful discovery and fine-mapping approach applied to smooth muscle cell and coronary artery disease-relevant tissues. We emphasize the critical need for alternative genetic strategies, including cis/trans-regulatory network analysis, genome editing, and perturbations, as well as single-cell sequencing in smooth muscle cell tissues and model organisms, under both normal and disease states. By integrating multiple experimental and analytical modalities, these multidimensional datasets should improve the interpretation of coronary artery disease genome-wide association studies and molecular quantitative trait locus signals and inform candidate targets for therapeutic intervention or risk prediction.

A Genome Wide Association Study on plasma FV levels identified PLXDC2 as a new modifier of the coagulation process

BACKGROUND

Factor V (FV) is a circulating protein primarily synthesized in the liver, and mainly present in plasma. It is a major component of the coagulation process.

OBJECTIVE

To detect novel genetic loci participating to the regulation of FV plasma levels.

METHODS

We conducted the first Genome Wide Association Study on FV plasma levels in a sample of 510 individuals and replicated the main findings in an independent sample of 1156 individuals.

RESULTS

In addition to genetic variations at the F5 locus, we identified novel associations at the PLXDC2 locus, with the lead PLXDC2 rs927826 polymorphism explaining ~3.7% (P = 7.5 × 10-15 in the combined discovery and replication samples) of the variability of FV plasma levels. In silico transcriptomic analyses in various cell types confirmed that PLXDC2 expression is positively correlated to F5 expression. SiRNA experiments in human hepatocellular carcinoma cell line confirmed the role of PLXDC2 in modulating factor F5 gene expression, and revealed further influences on F2 and F10 expressions.

CONCLUSION

Our study identified PLXDC2 as a new molecular player of the coagulation process.

Coronary Artery Disease Risk-Associated Plpp3 Gene and Its Product Lipid Phosphate Phosphatase 3 Regulate Experimental Atherosclerosis

OBJECTIVE

Genome-wide association studies identified novel loci in PLPP3(phospholipid phosphatase 3) that associate with coronary artery disease risk independently of traditional risk factors. PLPP3 encodes LPP3 (lipid phosphate phosphatase 3), a cell-surface enzyme that can regulate the availability of bioactive lysophopsholipids including lysophosphatidic acid (LPA). The protective allele of PLPP3 increases LPP3 expression during cell exposure to oxidized lipids, however, the role of LPP3 in atherosclerosis remains unclear. Approach and Results: In this study, we sought to validate LPP3 as a determinate of the development of atherosclerosis. In experimental models of atherosclerosis, LPP3 is upregulated and co-localizes with endothelial, smooth muscle cell, and CD68-positive cell markers. Global post-natal reductions in Plpp3 expression in mice substantially increase atherosclerosis, plaque-associated LPA, and inflammation. Although LPP3 expression increases during ox-LDL (oxidized low-density lipoprotein)-induced phenotypic modulation of bone marrow-derived macrophages, myeloid Plpp3 does not appear to regulate lesion formation. Rather, smooth muscle cell LPP3 expression is a critical regulator of atherosclerosis and LPA content in lesions. Moreover, mice with inherited deficiency in LPA receptor signaling are protected from experimental atherosclerosis.

CONCLUSIONS

Our results identify a novel lipid signaling pathway that regulates inflammation in the context of atherosclerosis and is not related to traditional risk factors. Pharmacological targeting of bioactive LPP3 substrates, including LPA, may offer an orthogonal approach to lipid-lowering drugs for mitigation of coronary artery disease risk.

Effects of diet and hyperlipidemia on levels and distribution of circulating lysophosphatidic acid

Lysophosphatidic acids (LPAs) are bioactive radyl hydrocarbon-substituted derivatives of glycerol 3-phosphate. LPA metabolism and signaling are implicated in heritable risk of coronary artery disease. Genetic and pharmacological inhibition of these processes attenuate experimental atherosclerosis. LPA accumulates in atheromas, which may be a consequence of association with LDLs. The source, regulation, and biological activity of LDL-associated LPA are unknown. We examined the effects of experimental hyperlipidemia on the levels and distribution of circulating LPA in mice. The majority of plasma LPA was associated with albumin in plasma from wild-type mice fed normal chow. LDL-associated LPA was increased in plasma from high-fat Western diet-fed mice that are genetically prone to hyperlipidemia (LDL receptor knockout or activated proprotein convertase subtilisin/kexin type 9-overexpressing C57Bl6). Adipose-specific deficiency of the ENPP2 gene encoding the LPA-generating secreted lysophospholipase D, autotaxin (ATX), attenuated these Western diet-dependent increases in LPA. ATX-dependent increases in LDL-associated LPA were observed in isolated incubated plasma. ATX acted directly on LDL-associated lysophospholipid substrates in vitro. LDL from all human subjects examined contained LPA and was decreased by lipid-lowering drug therapies. Human and mouse plasma therefore contains a diet-sensitive LDL-associated LPA pool that might contribute to the cardiovascular disease-promoting effects of LPA.

Regulation of PLPP3 gene expression by NF-κB family transcription factors

Lipid phosphate phosphatase 3 (LPP3), encoded by the PLPP3 gene, is an integral membrane enzyme that dephosphorylates phosphate esters of glycero- and sphingophospholipids. Cell surface LPP3 can terminate the signaling actions of bioactive lysophosphatidic acid (LPA) and sphingosine 1 phosphate, which likely explains its role in developmental angiogenesis, vascular injury responses, and cell migration. Heritable variants in the final intron PLPP3 associate with interindividual variability in coronary artery disease risk that may result from disruption of enhancer sequences that normally act in cis to increase expression of the gene. However, the mechanisms regulating PLPP3 expression are not well understood. We show that the human PLPP3 promoter contains three functional NF-κB response elements. All of these are required for maximal induction of PLPP3 promoter activity in reporter assays. The identified sequences recruit RelA and RelB components of the NF-κB transcription complex to chromatin, and these transcription factors bind to the identified target sequences in two different cell types. LPA promotes binding of Rel family transcription factors to the PLPP3 promoter and increases PLPP3 gene expression through mechanisms that are attenuated by an NF-κB inhibitor, LPA receptor antagonists, and inhibitors of phosphoinositide 3 kinase. These findings indicate that up-regulation of PLPP3 during inflammation and atherosclerosis results from canonical activation of the NF-κB signaling cascade to increase PLPP3 expression through nuclear import and binding of RelA and RelB transcription factors to the PLPP3 promoter and suggest a mechanism by which the LPP3 substrate, LPA, can regulate PLPP3 expression.

Challenges and opportunities in stroke genetics

Stroke, ischaemic stroke and subtypes of ischaemic stroke display substantial heritability. When compared with related vascular conditions, the number of established risk loci reaching genome-wide significance for association with stroke is still in the lower range, particularly for aetiological stroke subtypes such as large artery atherosclerotic stroke or small vessel stroke. Nevertheless, for individual loci substantial progress has been made in determining the specific mechanisms mediating stroke risk. In this review, we present a roadmap for functional follow-up of common risk variants associated with stroke. First, we discuss in silico strategies for characterizing signals in non-coding regions and highlight databases providing information on quantitative trait loci for mRNA and protein expression, as well as methylation, focussing on those with presumed relevance for stroke. Next, we discuss experimental strategies for following up on non-coding risk variants and regions such as massively parallel reporter assays, proteome-wide association studies, and chromatin conformation capture (3C) assays. These and other approaches are relevant for gaining insight into the specific variants and mechanisms mediating genetic stroke risk. Finally, we discuss how genetic findings could influence clinical practice by adding to diagnostic algorithms and eventually improve treatment options for stroke.

MiR-4262 promotes cell apoptosis and inhibits proliferation of colon cancer cells: involvement of GALNT4

A role of microRNA-4262 (miR-4262) in the carcinogenesis of colon cancer remains undetermined. In this study, we studied the effects and mechanisms of miR-4262 to the colon cancer cell proliferation and apoptosis. We found that the levels of miR-4262 significantly down-regulated in colon cancer tissue, compared to the paired adjacent non-tumor colon tissue. The miR-4262 levels in colon cancer cell lines were significantly lower than those in control normal colon tissues. Transfection with the miR-4262 mimic decreased the cell proliferation and increased cell apoptosis in colon cancer cells, while transfection with the antisense of miR-4262 (as-miR-4262) increased cell proliferation and suppressed cell apoptosis in colon cancer cells. Bioinformatics analyses showed that GALNT4 was a potential target gene of miR-4262. The luciferase activities assay and Western blot verified that miR-4262 targeted GALNT4 mRNA to modulate its protein levels. When we treated cells with miR-4262 and GALN4 siRNA, the cell viability was significantly decreased. Together, our study suggests that aberrantly expressed miR-4262 may affect cell apoptosis and proliferation of human colon cancer cells via GALNT4, which appears to be a promising therapeutic target for colon cancer.

Lysophosphatidic acid acts on LPA1 receptor to increase H2 O2 during flow-induced dilation in human adipose arterioles

BACKGROUND AND PURPOSE

NO produces arteriolar flow-induced dilation (FID) in healthy subjects but is replaced by mitochondria-derived hydrogen peroxide (mtH2 O2 ) in patients with coronary artery disease (CAD). Lysophosphatidic acid (LPA) is elevated in patients with risk factors for CAD, but its functional effect in arterioles is unknown. We tested whether elevated LPA changes the mediator of FID from NO to mtH2 O2 in human visceral and subcutaneous adipose arterioles.

EXPERIMENTAL APPROACH

Arterioles were cannulated on glass micropipettes and pressurized to 60 mmHg. We recorded lumen diameter after graded increases in flow in the presence of either NOS inhibition (L-NAME) or H2 O2 scavenging (Peg-Cat) ± LPA (10 μM, 30 min), ±LPA1 /LPA3 receptor antagonist (Ki16425) or LPA2 receptor antagonist (H2L5186303). We analysed LPA receptor RNA and protein levels in human arterioles and human cultured endothelial cells.

KEY RESULTS

FID was inhibited by L-NAME but not Peg-Cat in untreated vessels. In vessels treated with LPA, FID was of similar magnitude but inhibited by Peg-Cat while L-NAME had no effect. Rotenone attenuated FID in vessels treated with LPA indicating mitochondria as a source of ROS. RNA transcripts from LPA1 and LPA2 but not LPA3 receptors were detected in arterioles. LPA1 but not LPA3 receptor protein was detected by Western blot. Pretreatment of vessels with an LPA1 /LPA3 , but not LPA2 , receptor antagonist prior to LPA preserved NO-mediated dilation.

CONCLUSIONS AND IMPLICATIONS

These findings suggest an LPA1 receptor-dependent pathway by which LPA increases arteriolar release of mtH2 O2 as a mediator of FMD.

A systems-approach reveals human nestin is an endothelial-enriched, angiogenesis-independent intermediate filament protein

The intermediate filament protein nestin is expressed during embryonic development, but considered largely restricted to areas of regeneration in the adult. Here, we perform a body-wide transcriptome and protein-profiling analysis to reveal that nestin is constitutively, and highly-selectively, expressed in adult human endothelial cells (EC), independent of proliferative status. Correspondingly, we demonstrate that it is not a marker for tumour EC in multiple malignancy types. Imaging of EC from different vascular beds reveals nestin subcellular distribution is shear-modulated. siRNA inhibition of nestin increases EC proliferation, and nestin expression is reduced in atherosclerotic plaque neovessels. eQTL analysis reveals an association between SNPs linked to cardiovascular disease and reduced aortic EC nestin mRNA expression. Our study challenges the dogma that nestin is a marker of proliferation, and provides insight into its regulation and function in EC. Furthermore, our systems-based approach can be applied to investigate body-wide expression profiles of any candidate protein.

Transcription Factor Zhx2 Deficiency Reduces Atherosclerosis and Promotes Macrophage Apoptosis in Mice

Objective- The objective of this study was to determine the basis of resistance to atherosclerosis of inbred mouse strain BALB/cJ. Approach and Results- BALB/cJ mice carry a naturally occurring null mutation of the gene encoding the transcription factor Zhx2, and genetic analyses suggested that this may confer resistance to atherosclerosis. On a hyperlipidemic low-density lipoprotein receptor null background, BALB/cJ mice carrying the mutant allele for Zhx2 exhibited up to a 10-fold reduction in lesion size as compared with an isogenic strain carrying the wild-type allele. Several lines of evidence, including bone marrow transplantation studies, indicate that this effect of Zhx2 is mediated, in part, by monocytes/macrophages although nonbone marrow-derived pathways are clearly involved as well. Both in culture and in atherosclerotic lesions, macrophages from Zhx2 null mice exhibited substantially increased apoptosis. Zhx2 null macrophages were also enriched for M2 markers. Effects of Zhx2 on proliferation and other bone marrow-derived cells, such as lymphocytes, were at most modest. Expression microarray analyses identified >1000 differentially expressed transcripts between Zhx2 wild-type and null macrophages. To identify the global targets of Zhx2, we performed ChIP-seq (chromatin immunoprecipitation sequencing) studies with the macrophage cell line RAW264.7. The ChIP-seq peaks overlapped significantly with gene expression and together suggested roles for transcriptional repression and apoptosis. Conclusions- A mutation of Zhx2 carried in BALB/cJ mice is responsible in large part for its relative resistance to atherosclerosis. Our results indicate that Zhx2 promotes macrophage survival and proinflammatory functions in atherosclerotic lesions, thereby contributing to lesion growth.

Integrated Human Evaluation of the Lysophosphatidic Acid Pathway as a Novel Therapeutic Target in Atherosclerosis

Variants in the PLPP3 gene encoding for lipid phosphate phosphohydrolase 3 have been associated with susceptibility to atherosclerosis independently of classical risk factors. PLPP3 inactivates lysophosphatidic acid (LPA), a pro-inflammatory, pro-thrombotic product of phospholipase activity. Here we performed the first exploratory analysis of PLPP3, LPA, and LPA receptors (LPARs 1–6) in human atherosclerosis. PLPP3 transcript and protein were repressed when comparing plaques versus normal arteries and plaques from symptomatic versus asymptomatic patients, and they were negatively associated with risk of adverse cardiovascular events. PLPP3 localized to macrophages, smooth muscle, and endothelial cells (ECs) in plaques. LPAR 2, 5, and especially 6 showed increased expression in plaques, with LPAR6 localized in ECs and positively correlated to PLPP3. Utilizing in situ mass spectrometry imaging, LPA and its precursors were found in the plaque fibrous cap, co-localizing with PLPP3 and LPAR6. In vitro, PLPP3 silencing in ECs under LPA stimulation resulted in increased expression of adhesion molecules and cytokines. LPAR6 silencing inhibited LPA-induced cell activation, but not when PLPP3 was silenced simultaneously. Our results show that repression of PLPP3 plays a key role in atherosclerosis by promoting EC activation. Altogether, the PLPP3 pathway represents a suitable target for investigations into novel therapeutic approaches to ameliorate atherosclerosis.

Oxidized phospholipids regulate amino acid metabolism through MTHFD2 to facilitate nucleotide release in endothelial cells

Oxidized phospholipids (oxPAPC) induce endothelial dysfunction and atherosclerosis. Here we show that oxPAPC induce a gene network regulating serine-glycine metabolism with the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) as a causal regulator using integrative network modeling and Bayesian network analysis in human aortic endothelial cells. The cluster is activated in human plaque material and by atherogenic lipoproteins isolated from plasma of patients with coronary artery disease (CAD). Single nucleotide polymorphisms (SNPs) within the MTHFD2-controlled cluster associate with CAD. The MTHFD2-controlled cluster redirects metabolism to glycine synthesis to replenish purine nucleotides. Since endothelial cells secrete purines in response to oxPAPC, the MTHFD2-controlled response maintains endothelial ATP. Accordingly, MTHFD2-dependent glycine synthesis is a prerequisite for angiogenesis. Thus, we propose that endothelial cells undergo MTHFD2-mediated reprogramming toward serine-glycine and mitochondrial one-carbon metabolism to compensate for the loss of ATP in response to oxPAPC during atherosclerosis.

Blood lipid-related low-frequency variants in LDLR and PCSK9 are associated with onset age and risk of myocardial infarction in Japanese

Recent studies have revealed the importance of rare variants in myocardial infarction (MI) susceptibility in European populations. Because genetic architectures vary in different populations, we investigated how they contribute to MI susceptibility in Japanese subjects. We performed targeted sequencing of 36 coronary artery disease risk genes, identified by genome-wide association studies, in 9,956 cases and 8,373 controls. Gene-based association tests identified significant enrichment of rare variants in LDLR and PCSK9 in MI cases. We identified 52 (novel 22) LDLR variants predicted to be damaging. Carriers of these variants showed a higher risk of MI (carriers/non-carriers 89/9867 in cases, 17/8356 controls, OR = 4.4, P = 7.2 × 10⁻¹⁰), higher LDL-cholesterol levels and younger age of onset for MI. With respect to PCSK9, E32K carriers showed higher LDL-cholesterol levels and younger age of onset for MI, whereas R93C carriers had lower LDL-cholesterol levels. A significant correlation between LDL-cholesterol levels and onset age of MI was observed in these variant carriers. In good agreement with previous studies in patients with familial hypercholesterolaemia, our study in the Japanese general population showed that rare variants in LDLR and PCSK9 were associated with the onset age of MI by altering LDL-cholesterol levels.

Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

Atherosclerosis Is an Inflammatory Disease which Lacks a Common Anti-inflammatory Therapy: How Human Genetics Can Help to This Issue. A Narrative Review

Atherosclerosis is a multifactorial disease triggered and sustained by different risk factors such as dyslipidemia, arterial hypertension, diabetes mellitus, smoke, etc. Since a couple of decades, a pivotal role for inflammation in its pathogenesis has been recognized and proved at molecular levels, and already described in many animal models. Despite all this knowledge, due to the complexity of the specific inflammatory process subtending atherosclerosis and to the fact that inflammation is also a protective response against microorganisms, no anti-inflammatory therapy has been rendered available in the therapeutic armamentarium against atherosclerosis and vascular events till 2017 when canakinumab in the first ad-hoc randomized clinical trial (RCT) proved for the first time that targeting specifically inflammation lowers cardiovascular (CV) events. From the genetic side, in the 90's and early 2000, several genetic markers in inflammatory pathway have been explored searching for an association with athero-thrombosis which gave seldom consistent results. Then, in the genomic era, plenty of genetic markers covering most of the genome have been analyzed at once without a priori information. The results coming from genome wide association studies (GWAS) have pinpointed some loci closed to inflammatory molecules consistently associated with atherosclerosis and CV consequences revamping the strict link between inflammation and atherosclerosis and suggesting some tailored target therapy. Whole-exome and whole-genome sequencing will come soon showing new and old loci associated with atherosclerosis suggesting new molecular targets or underlying which inflammatory pathway could be most attractive to target for blocking atherosclerosis even in its early stages.

Monocyte-Macrophages and T Cells in Atherosclerosis

Atherosclerosis is an arterial disease process characterized by the focal subendothelial accumulation of apolipoprotein-B-containing lipoproteins, immune and vascular wall cells, and extracellular matrix. The lipoproteins acquire features of damage-associated molecular patterns and trigger first an innate immune response, dominated by monocyte-macrophages, and then an adaptive immune response. These inflammatory responses often become chronic and non-resolving and can lead to arterial damage and thrombosis-induced organ infarction. The innate immune response is regulated at various stages, from hematopoiesis to monocyte changes and macrophage activation. The adaptive immune response is regulated primarily by mechanisms that affect the balance between regulatory and effector T cells. Mechanisms related to cellular cholesterol, phenotypic plasticity, metabolism, and aging play key roles in affecting these responses. Herein, we review select topics that shed light on these processes and suggest new treatment strategies.

Consensus Genome-Wide Expression Quantitative Trait Loci and Their Relationship with Human Complex Trait Disease

Most of the risk loci identified from genome-wide association (GWA) studies do not provide direct information on the biological basis of a disease or on the underlying mechanisms. Recent expression quantitative trait locus (eQTL) association studies have provided information on genetic factors associated with gene expression variation. These eQTLs might contribute to phenotype diversity and disease susceptibility, but interpretation is handicapped by low reproducibility of the expression results. To address this issue, we have generated a set of consensus eQTLs by integrating publicly available data for specific human populations and cell types. Overall, we find over 4000 genes that are involved in high-confidence eQTL relationships. To elucidate the role that eQTLs play in human common diseases, we matched the high-confidence eQTLs to a set of 335 disease risk loci identified from the Wellcome Trust Case Control Consortium GWA study and follow-up studies for 7 human complex trait diseases-bipolar disorder (BD), coronary artery disease (CAD), Crohn's disease (CD), hypertension (HT), rheumatoid arthritis (RA), type 1 diabetes (T1D), and type 2 diabetes (T2D). The results show that the data are consistent with ∼50% of these disease loci arising from an underlying expression change mechanism.

Sex-specific association of SH2B3 and SMARCA4 polymorphisms with coronary artery disease susceptibility

To determine whether sex differences affect the association between genetic polymorphisms and coronary artery disease (CAD) in the Chinese Han population, we conducted a study comparing the frequency of SH2B3 and SMARCA4 variants in 456 CAD patients (291 men, 165 women) and 685 age-matched controls (385 men, 300 women). Ten single nucleotide polymorphisms (SNPs) in SH2B3 and SMARCA4 were genotyped using MassARRAY technology. Allelic and genotypic models and haplotype frequencies were compared between groups. Logistic regression was used to estimate the CAD risk associated with the genotypes. We found that the “A” alleles in both rs11879293 and rs12232780 of SMARCA4 were associated with CAD risk in men (p = 0.036 and p = 0.001, respectively). The genetic model showed that SH2B3 was associated with CAD susceptibility in both women and men, while SMARCA4 was associated with reduced odds of CAD in men. SH2B3 haplotypes were associated with decreased CAD risk in women (p = 0.007) and increased CAD risk in men (p = 0.047). By providing evidence for the sex-related association between SH2B3 and SMARCA4 gene variants and CAD susceptibility in the Chinese Han population, this study may help define useful diagnostic and preventive markers for these patients.

Loss of Cardioprotective Effects at the ADAMTS7 Locus as a Result of Gene-Smoking Interactions

BACKGROUND

Common diseases such as coronary heart disease (CHD) are complex in etiology. The interaction of genetic susceptibility with lifestyle factors may play a prominent role. However, gene-lifestyle interactions for CHD have been difficult to identify. Here, we investigate interaction of smoking behavior, a potent lifestyle factor, with genotypes that have been shown to associate with CHD risk.

METHODS

We analyzed data on 60 919 CHD cases and 80 243 controls from 29 studies for gene-smoking interactions for genetic variants at 45 loci previously reported to be associated with CHD risk. We also studied 5 loci associated with smoking behavior. Study-specific gene-smoking interaction effects were calculated and pooled using fixed-effects meta-analyses. Interaction analyses were declared to be significant at a P value of <1.0×10-3 (Bonferroni correction for 50 tests).

RESULTS

We identified novel gene-smoking interaction for a variant upstream of the ADAMTS7 gene. Every T allele of rs7178051 was associated with lower CHD risk by 12% in never-smokers (P=1.3×10-16) in comparison with 5% in ever-smokers (P=2.5×10-4), translating to a 60% loss of CHD protection conferred by this allelic variation in people who smoked tobacco (interaction P value=8.7×10-5). The protective T allele at rs7178051 was also associated with reduced ADAMTS7 expression in human aortic endothelial cells and lymphoblastoid cell lines. Exposure of human coronary artery smooth muscle cells to cigarette smoke extract led to induction of ADAMTS7. CONCLUSIONS: Allelic variation at rs7178051 that associates with reduced ADAMTS7 expression confers stronger CHD protection in never-smokers than in ever-smokers. Increased vascular ADAMTS7 expression may contribute to the loss of CHD protection in smokers.

Association between Helicobacter pylori infection and angiographically demonstrated coronary artery disease: A meta-analysis

Coronary artery disease (CAD) is a leading cause of mortality globally. However, the etiology and pathogenesis of CAD are not fully understood. The aim of the present meta-analysis was to estimate the association between the risk of CAD and Helicobacter pylori (H. pylori) infection. A literature search was performed to identify eligible studies published prior to August 14, 2014. Fixed or random effect meta-analytical methods were used to pool the data and perform the subgroup analyses. The effect measures estimated were the odds ratios (OR) for dichotomous data reported with 95% confidence intervals (95% CI). Of the 109 studies identified using the search parameters, 26 cross-sectional studies were eligible involving 3,901 CAD patients and 2,751 controls. H. pylori infection was associated with an increased risk of CAD (OR: 1.96, 95% CI: 1.47-2.63, P<0.00001). When the adjusted ORs were used to conduct another meta-analysis, the OR value decreased, but the association remained significant (OR: 1.42, 95% CI: 1.09-1.86, P=0.008). The association between H. pylori infection and CAD risk was stronger in younger individuals than in older individuals (OR: 2.36, 95% CI 1.50-3.73 vs. OR: 1.59, 95% CI: 1.19-2.11). A significant association was observed in studies from Europe (OR: 2.11, 95% CI: 1.54-2.88, P=0.01) and the USA (OR: 1.43, 95% CI: 1.08-1.91, P=0.36). There is a potential association between H. pylori infection and the risk of CAD. The association may be influenced by age and ethnicity.

Functional Validation of a Common Nonsynonymous Coding Variant in ZC3HC1 Associated With Protection From Coronary Artery Disease

Background—

Although virtually all coronary artery disease associated single-nucleotide polymorphisms identified by genome-wide association studies (GWAS) are in noncoding regions of the genome, a common polymorphism in ZC3HC1 (rs11556924), resulting in an arginine (Arg) to histidine (His) substitution in its encoded protein, NIPA (Nuclear Interacting Partner of Anaplastic Lyphoma Kinase) is linked to a protection from coronary artery disease. NIPA plays a role in cell cycle progression, but the functional consequences of this polymorphism have not been established.

Methods and Results—

Here we demonstrate that total ZC3HC1 expression in whole blood is similar across genotypes, despite expression being slightly biased toward the risk allele in heterozygotes. At the protein level, the protective His363 NIPA variant exhibits increased phosphorylation of a critical serine residue (Ser354) and higher protein expression as compared with the Arg363 variant. Binding experiments indicate that neither SKP1 (S-phase kinase-associated protein 1) nor CCNB1 binding were affected by the polymorphism. Despite similar nuclear distribution, NIPA His363 exhibits greater nuclear mobility. NIPA suppression results in a modest reduction of proliferation in vascular smooth muscle cells, but given low proliferative capacity, a significant effect of the variant was not noted. By contrast, we demonstrate that the protective variant reduces cell proliferation in HeLa cells.

Conclusions—

These findings extend the genetic association between rs11556924 and coronary artery disease risk by characterizing its effects on the encoded protein, NIPA. The resulting amino acid change Arg363His is associated with increased expression and nuclear mobility, as well as lower rates of cell growth in HeLa cells, further supporting a role for cell proliferation in atherosclerosis and its clinical consequences.

Different origins of lysophospholipid mediators between coronary and peripheral arteries in acute coronary syndrome

Lysophosphatidic acids (LysoPAs) and lysophosphatidylserine (LysoPS) are emerging lipid mediators proposed to be involved in the pathogenesis of acute coronary syndrome (ACS). In this study, we attempted to elucidate how LysoPA and LysoPS become elevated in ACS using human blood samples collected simultaneously from culprit coronary arteries and peripheral arteries in ACS subjects. We found that: 1) the plasma LysoPA, LysoPS, and lysophosphatidylglycerol levels were not different, while the lysophosphatidylcholine (LysoPC), lysophosphatidylinositol, and lysophosphatidylethanolamine (LysoPE) levels were significantly lower in the culprit coronary arteries; 2) the serum autotaxin (ATX) level was lower and the serum phosphatidylserine-specific phospholipase A₁ (PS-PLA₁) level was higher in the culprit coronary arteries; 3) the LysoPE and ATX levels were significant explanatory factors for the mainly elevated species of LysoPA, except for 22:6 LysoPA, in the peripheral arteries, while the LysoPC and LysoPE levels, but not the ATX level, were explanatory factors in the culprit coronary arteries; and 4) 18:0 and 18:1 LysoPS were significantly correlated with PS-PLA₁ only in the culprit coronary arteries. In conclusion, the origins of LysoPA and LysoPS might differ between culprit coronary arteries and peripheral arteries, and substrates for ATX, such as LysoPC and LysoPE, might be important for the generation of LysoPA in ACS.

Identification of Risk Pathways and Functional Modules for Coronary Artery Disease Based on Genome-wide SNP Data

Coronary artery disease (CAD) is a complex human disease, involving multiple genes and their nonlinear interactions, which often act in a modular fashion. Genome-wide single nucleotide polymorphism (SNP) profiling provides an effective technique to unravel these underlying genetic interplays or their functional involvements for CAD. This study aimed to identify the susceptible pathways and modules for CAD based on SNP omics. First, the Wellcome Trust Case Control Consortium (WTCCC) SNP datasets of CAD and control samples were used to assess the joint effect of multiple genetic variants at the pathway level, using logistic kernel machine regression model. Then, an expanded genetic network was constructed by integrating statistical gene–gene interactions involved in these susceptible pathways with their protein–protein interaction (PPI) knowledge. Finally, risk functional modules were identified by decomposition of the network. Of 276 KEGG pathways analyzed, 6 pathways were found to have a significant effect on CAD. Other than glycerolipid metabolism, glycosaminoglycan biosynthesis, and cardiac muscle contraction pathways, three pathways related to other diseases were also revealed, including Alzheimer’s disease, non-alcoholic fatty liver disease, and Huntington’s disease. A genetic epistatic network of 95 genes was further constructed using the abovementioned integrative approach. Of 10 functional modules derived from the network, 6 have been annotated to phospholipase C activity and cell adhesion molecule binding, which also have known functional involvement in Alzheimer’s disease. These findings indicate an overlap of the underlying molecular mechanisms between CAD and Alzheimer’s disease, thus providing new insights into the molecular basis for CAD and its molecular relationships with other diseases.

Single nucleotide polymorphisms in an intergenic chromosome 2q region associated with tissue factor pathway inhibitor plasma levels and venous thromboembolism

Essentials Tissue factor pathway inhibitor (TFPI) regulates the blood coagulation cascade. We replicated previously reported linkage of TFPI plasma levels to the chromosome 2q region. The putative causal locus, rs62187992, was associated with TFPI plasma levels and thrombosis. rs62187992 was marginally associated with TFPI expression in human aortic endothelial cells. Click to hear Ann Gil's presentation on new insights into thrombin activatable fibrinolysis inhibitor SUMMARY: Background Tissue factor pathway inhibitor (TFPI) regulates fibrin clot formation, and low TFPI plasma levels increase the risk of arterial thromboembolism and venous thromboembolism (VTE). TFPI plasma levels are also heritable, and a previous linkage scan implicated the chromosome 2q region, but no specific genes. Objectives To replicate the finding of the linkage region in an independent sample, and to identify the causal locus. Methods We first performed a linkage analysis of microsatellite markers and TFPI plasma levels in 251 individuals from the F5L Family Study, and replicated the finding of the linkage peak on chromosome 2q (LOD = 3.06). We next defined a follow-up region that included 112 603 single nucleotide polymorphisms (SNPs) under the linkage peak, and meta-analyzed associations between these SNPs and TFPI plasma levels across the F5L Family Study and the Marseille Thrombosis Association (MARTHA) Study, a study of 1033 unrelated VTE patients. SNPs with false discovery rate q-values of < 0.10 were tested for association with TFPI plasma levels in 892 patients with coronary artery disease in the AtheroGene Study. Results and Conclusions One SNP, rs62187992, was associated with TFPI plasma levels in all three samples (β = + 0.14 and P = 4.23 × 10-6 combined; β = + 0.16 and P = 0.02 in the F5L Family Study; β = + 0.13 and P = 6.3 × 10-4 in the MARTHA Study; β = + 0.17 and P = 0.03 in the AtheroGene Study), and contributed to the linkage peak in the F5L Family Study. rs62187992 was also associated with clinical VTE (odds ratio 0.90, P = 0.03) in the INVENT Consortium of > 7000 cases and their controls, and was marginally associated with TFPI expression (β = + 0.19, P = 0.08) in human aortic endothelial cells, a primary site of TFPI synthesis. The biological mechanisms underlying these associations remain to be elucidated.

Role of lipid phosphate phosphatase 3 in human aortic endothelial cell function

AIMS

Lipid phosphate phosphatase 3; type 2 phosphatidic acid phosphatase β (LPP3; PPAP2B) is a transmembrane protein dephosphorylating and thereby terminating signalling of lipid substrates including lysophosphatidic acid (LPA) and sphingosine-1-phosphate (S1P). Human LPP3 possesses a cell adhesion motif that allows interaction with integrins. A polymorphism (rs17114036) in PPAP2B is associated with coronary artery disease, which prompted us to investigate the possible role of LPP3 in human endothelial dysfunction, a condition promoting atherosclerosis.

METHODS AND RESULTS

To study the role of LPP3 in endothelial cells we used human primary aortic endothelial cells (HAECs) in which LPP3 was silenced or overexpressed using either wild type or mutated cDNA constructs. LPP3 silencing in HAECs enhanced secretion of inflammatory cytokines, leucocyte adhesion, cell survival, and migration and impaired angiogenesis, whereas wild-type LPP3 overexpression reversed these effects and induced apoptosis. We also demonstrated that LPP3 expression was negatively correlated with vascular endothelial growth factor expression. Mutations in either the catalytic or the arginine-glycine-aspartate (RGD) domains impaired endothelial cell function and pharmacological inhibition of S1P or LPA restored it. LPA was not secreted in HAECs under silencing or overexpressing LPP3. However, the intra- and extra-cellular levels of S1P tended to be correlated with LPP3 expression, indicating that S1P is probably degraded by LPP3.

CONCLUSIONS

We demonstrated that LPP3 is a negative regulator of inflammatory cytokines, leucocyte adhesion, cell survival, and migration in HAECs, suggesting a protective role of LPP3 against endothelial dysfunction in humans. Both the catalytic and the RGD functional domains were involved and S1P, but not LPA, might be the endogenous substrate of LPP3.

The Coronary Artery Disease-associated Coding Variant in Zinc Finger C3HC-type Containing 1 (ZC3HC1) Affects Cell Cycle Regulation

Genome-wide association studies have to date identified multiple coronary artery disease (CAD)-associated loci; however, for most of these loci the mechanism by which they affect CAD risk is unclear. The CAD-associated locus 7q32.2 is unusual in that the lead variant, rs11556924, is not in strong linkage disequilibrium with any other variant and introduces a coding change in ZC3HC1, which encodes NIPA. In this study, we show that rs11556924 polymorphism is associated with lower regulatory phosphorylation of NIPA in the risk variant, resulting in NIPA with higher activity. Using a genome-editing approach we show that this causes an effective decrease in cyclin-B1 stability in the nucleus, thereby slowing its nuclear accumulation. By perturbing the rate of nuclear cyclin-B1 accumulation, rs11556924 alters the regulation of mitotic progression resulting in an extended mitosis. This study shows that the CAD-associated coding polymorphism in ZC3HC1 alters the dynamics of cell-cycle regulation by NIPA.

Genetic and environmental risk factors for atherosclerosis regulate transcription of phosphatase and actin regulating gene PHACTR1

Background and aims

Coronary artery disease (CAD) risk is associated with non-coding genetic variants at the phosphatase and actin regulating protein 1(PHACTR1) gene locus. The PHACTR1 gene encodes an actin-binding protein with phosphatase regulating activity. The mechanism whereby PHACTR1 influences CAD risk is unknown. We hypothesized that PHACTR1 would be expressed in human cell types relevant to CAD and regulated by atherogenic or genetic factors.

Methods and results

Using immunohistochemistry, we demonstrate that PHACTR1 protein is expressed strongly in human atherosclerotic plaque macrophages, lipid-laden foam cells, adventitial lymphocytes and endothelial cells. Using a combination of genomic analysis and molecular techniques, we demonstrate that PHACTR1 is expressed as multiple previously uncharacterized transcripts in macrophages, foam cells, lymphocytes and endothelial cells. Immunoblotting confirmed a total absence of PHACTR1 in vascular smooth muscle cells. Real-time quantitative PCR showed that PHACTR1 is regulated by atherogenic and inflammatory stimuli. In aortic endothelial cells, oxLDL and TNF-alpha both upregulated an intermediate length transcript. A short transcript expressed only in immune cells was upregulated in macrophages by oxidized low-density lipoprotein, and oxidized phospholipids but suppressed by lipopolysaccharide or TNF-alpha. In primary human macrophages, we identified a novel expression quantitative trait locus (eQTL) specific for this short transcript, whereby the risk allele at CAD risk SNP rs9349379 is associated with reduced PHACTR1 expression, similar to the effect of an inflammatory stimulus.

Conclusions

Our data demonstrate that PHACTR1 is a key atherosclerosis candidate gene since it is regulated by atherogenic stimuli in macrophages and endothelial cells and we identify an effect of the genetic risk variant on PHACTR1 expression in macrophages that is similar to that of an inflammatory stimulus.

•

PHACTR1 is expressed as two transcripts in both immune and endothelial cells in human atherosclerotic plaque.

•

Oxidized-LDL upregulates a short PHACTR1 transcript, but suppresses an intermediate length transcript in macrophages.

•

Lipopolysaccharide and TNF-alpha cause the opposite effect with strong suppression of the short transcript in macrophages.

•

The coronary artery disease risk SNP, rs9349379, is associated with expression of the short transcript in macrophages.

•

The effect of the coronary artery disease risk allele on PHACTR1 mirrors that of inflammatory stimuli.

Network-Based Identification and Prioritization of Key Regulators of Coronary Artery Disease Loci

OBJECTIVE

Recent genome-wide association studies of coronary artery disease (CAD) have revealed 58 genome-wide significant and 148 suggestive genetic loci. However, the molecular mechanisms through which they contribute to CAD and the clinical implications of these findings remain largely unknown. We aim to retrieve gene subnetworks of the 206 CAD loci and identify and prioritize candidate regulators to better understand the biological mechanisms underlying the genetic associations.

APPROACH AND RESULTS

We devised a new integrative genomics approach that incorporated (1) candidate genes from the top CAD loci, (2) the complete genetic association results from the 1000 genomes-based CAD genome-wide association studies from the Coronary Artery Disease Genome Wide Replication and Meta-Analysis Plus the Coronary Artery Disease consortium, (3) tissue-specific gene regulatory networks that depict the potential relationship and interactions between genes, and (4) tissue-specific gene expression patterns between CAD patients and controls. The networks and top-ranked regulators according to these data-driven criteria were further queried against literature, experimental evidence, and drug information to evaluate their disease relevance and potential as drug targets. Our analysis uncovered several potential novel regulators of CAD such as LUM and STAT3, which possess properties suitable as drug targets. We also revealed molecular relations and potential mechanisms through which the top CAD loci operate. Furthermore, we found that multiple CAD-relevant biological processes such as extracellular matrix, inflammatory and immune pathways, complement and coagulation cascades, and lipid metabolism interact in the CAD networks.

CONCLUSIONS

Our data-driven integrative genomics framework unraveled tissue-specific relations among the candidate genes of the CAD genome-wide association studies loci and prioritized novel network regulatory genes orchestrating biological processes relevant to CAD.