Convergence of coronary artery disease genes onto endothelial cell programs.

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

Protein interaction networks in the vasculature prioritize genes and pathways underlying coronary artery disease.

Population-based association studies have identified many genetic risk loci for coronary artery disease (CAD), but it is often unclear how genes within these loci are linked to CAD. Here, we perform interaction proteomics for 11 CAD-risk genes to map their protein-protein interactions (PPIs) in human vascular cells and elucidate their roles in CAD. The resulting PPI networks contain interactions that are outside of known biology in the vasculature and are enriched for genes involved in immunity-related and arterial-wall-specific mechanisms. Several PPI networks derived from smooth muscle cells are significantly enriched for genetic variants associated with CAD and related vascular phenotypes. Furthermore, the networks identify 61 genes that are found in genetic loci associated with risk of CAD, prioritizing them as the causal candidates within these loci. These findings indicate that the PPI networks we have generated are a rich resource for guiding future research into the molecular pathogenesis of CAD.

The Role of Endothelial Cells in Atherosclerosis: Insights from Genetic Association Studies.

Endothelial cells (ECs) mediate several biological functions that are relevant to atherosclerosis and coronary artery disease (CAD), regulating an array of vital processes including vascular tone, wound healing, reactive oxygen species, shear stress response, and inflammation. Although which of these functions is linked causally with CAD development and/or progression is not yet known, genome-wide association studies have implicated more than 400 loci associated with CAD risk, among which several have shown EC-relevant functions. Given the arduous process of mechanistically interrogating single loci to CAD, high-throughput variant characterization methods, including pooled Clustered Regularly Interspaced Short Palindromic Repeats screens, offer exciting potential to rapidly accelerate the discovery of bona fide EC-relevant genetic loci. These discoveries in turn will broaden the therapeutic avenues for CAD beyond lipid lowering and behavioral risk modification to include EC-centric modalities of risk prevention and treatment.

Leveraging polygenic enrichments of gene features to predict genes underlying complex traits and diseases.

Genome-wide association studies (GWASs) are a valuable tool for understanding the biology of complex human traits and diseases, but associated variants rarely point directly to causal genes. In the present study, we introduce a new method, polygenic priority score (PoPS), that learns trait-relevant gene features, such as cell-type-specific expression, to prioritize genes at GWAS loci. Using a large evaluation set of genes with fine-mapped coding variants, we show that PoPS and the closest gene individually outperform other gene prioritization methods, but observe the best overall performance by combining PoPS with orthogonal methods. Using this combined approach, we prioritize 10,642 unique gene-trait pairs across 113 complex traits and diseases with high precision, finding not only well-established gene-trait relationships but nominating new genes at unresolved loci, such as LGR4 for estimated glomerular filtration rate and CCR7 for deep vein thrombosis. Overall, we demonstrate that PoPS provides a powerful addition to the gene prioritization toolbox.

Polygenic prediction of preeclampsia and gestational hypertension.

Preeclampsia and gestational hypertension are common pregnancy complications associated with adverse maternal and child outcomes. Current tools for prediction, prevention and treatment are limited. Here we tested the association of maternal DNA sequence variants with preeclampsia in 20,064 cases and 703,117 control individuals and with gestational hypertension in 11,027 cases and 412,788 control individuals across discovery and follow-up cohorts using multi-ancestry meta-analysis. Altogether, we identified 18 independent loci associated with preeclampsia/eclampsia and/or gestational hypertension, 12 of which are new (for example, MTHFR-CLCN6, WNT3A, NPR3, PGR and RGL3), including two loci (PLCE1 and FURIN) identified in the multitrait analysis. Identified loci highlight the role of natriuretic peptide signaling, angiogenesis, renal glomerular function, trophoblast development and immune dysregulation. We derived genome-wide polygenic risk scores that predicted preeclampsia/eclampsia and gestational hypertension in external cohorts, independent of clinical risk factors, and reclassified eligibility for low-dose aspirin to prevent preeclampsia. Collectively, these findings provide mechanistic insights into the hypertensive disorders of pregnancy and have the potential to advance pregnancy risk stratification.

Multiomic Analysis and CRISPR Perturbation Screens Identify Endothelial Cell Programs and Novel Therapeutic Targets for Coronary Artery Disease.

Endothelial cells (EC) are an important mediator of atherosclerosis and vascular disease. Their exposure to atherogenic risk factors such as hypertension and serum cholesterol leads to endothelial dysfunction and many disease-associated processes. Identifying which of these multiple EC functions is causally related to disease risk has been challenging. There is evidence from in vivo models and human sequencing studies that dysregulation of nitric oxide production directly affects risk of coronary artery disease. Human genetics can help prioritize the other EC functions with causal relationships because germline mutations are acquired at birth and serve as a randomized test of which pathways affect disease risk. Though several coronary artery disease risk variants have been linked to EC function, this process has been slow and laborious. Unbiased analyses of EC dysfunction using multiomic approaches promise to identify the causal genetic mechanisms responsible for vascular disease. Here, we review the data from genomic, epigenomic, and transcriptomic studies that prioritize EC-specific causal pathways. New methods that CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technology with genomic, epigenomic, and transcriptomic analysis promise to speed up the characterization of disease-associated genetic variation. We summarize several recent studies in ECs which use high-throughput genetic perturbation to identify disease-relevant pathways and novel mechanisms of disease. These genetically validated pathways can accelerate the identification of drug targets for the prevention and treatment of atherosclerosis.

Endothelial cell-specific deletion of a microRNA accelerates atherosclerosis.

BACKGROUND AND AIMS: Chronic vascular endothelial inflammation predisposes to atherosclerosis; however, the cell-autonomous roles for endothelial-expressing microRNAs (miRNAs) are poorly understood in this process. MiR-181b is expressed in several cellular constituents relevant to lesion formation. The aim of this study is to examine the role of genetic deficiency of the miR-181b locus in endothelial cells during atherogenesis. METHODS AND RESULTS: Using a proprotein convertase subtilisin/kexin type 9 (PCSK9)-induced atherosclerosis mouse model, we demonstrated that endothelial cell (EC)-specific deletion of miR-181a2b2 significantly promoted atherosclerotic lesion formation, cell adhesion molecule expression, and the influx of lesional macrophages in the vessel wall. Yet, endothelium deletion of miR-181a2b2 did not affect body weight, lipid metabolism, anti-inflammatory Ly6Clow or the pro-inflammatory Ly6Cinterm and Ly6Chigh fractions in circulating peripheral blood mononuclear cells (PBMCs), and pro-inflammatory or anti-inflammatory mediators in both bone marrow (BM) and PBMCs. Mechanistically, bulk RNA-seq and gene set enrichment analysis of ECs enriched from the aortic arch intima, as well as single cell RNA-seq from atherosclerotic lesions, revealed that endothelial miR-181a2b2 serves as a critical regulatory hub in controlling endothelial inflammation, cell adhesion, cell cycle, and immune response during atherosclerosis. CONCLUSIONS: Our study establishes that deficiency of a miRNA specifically in the vascular endothelium is sufficient to profoundly impact atherogenesis. Endothelial miR-181a2b2 deficiency regulates multiple key pathways related to endothelial inflammation, cell adhesion, cell cycle, and immune response involved in the development of atherosclerosis.

Deep learning enables genetic analysis of the human thoracic aorta.

Enlargement or aneurysm of the aorta predisposes to dissection, an important cause of sudden death. We trained a deep learning model to evaluate the dimensions of the ascending and descending thoracic aorta in 4.6 million cardiac magnetic resonance images from the UK Biobank. We then conducted genome-wide association studies in 39,688 individuals, identifying 82 loci associated with ascending and 47 with descending thoracic aortic diameter, of which 14 loci overlapped. Transcriptome-wide analyses, rare-variant burden tests and human aortic single nucleus RNA sequencing prioritized genes including SVIL, which was strongly associated with descending aortic diameter. A polygenic score for ascending aortic diameter was associated with thoracic aortic aneurysm in 385,621 UK Biobank participants (hazard ratio = 1.43 per s.d., confidence interval 1.32-1.54, P = 3.3 × 10-20). Our results illustrate the potential for rapidly defining quantitative traits with deep learning, an approach that can be broadly applied to biomedical images.

CRISPR-Cas9 Genome Editing of Primary Human Vascular Cells In Vitro.

Genome editing of primary human cells with CRISPR-Cas9 is a powerful tool to study gene function. For many cell types, there are efficient protocols for editing with optimized plasmids for Cas9 and sgRNA expression. Vascular cells, however, remain refractory to plasmid-based delivery of CRISPR machinery for in vitro genome editing due to low transfection efficiency, poor expression of the Cas9 machinery, and toxic effects of the selection antibiotics. Here, we describe a method for high-efficiency editing of primary human vascular cells in vitro using nucleofection for direct delivery of sgRNA:Cas9-NLS ribonucleoprotein complexes. This method is more rapid and its high editing efficiency eliminates the need for additional selection steps. The edited cells can be employed in diverse applications, such as gene expression measurement or functional assays to assess various genetic perturbation effects in vitro. This method proves effective in vascular cells that are refractory to standard genome manipulation techniques using viral plasmid delivery. We anticipate that this technique will be applied to other non-vascular cell types that face similar barriers to efficient genome editing. © 2021 Wiley Periodicals LLC. Basic Protocol: CRISPR-Cas9 genome editing of primary human vascular cells in vitro.

Rare, Damaging DNA Variants in CORIN and Risk of Coronary Artery Disease: Insights From Functional Genomics and Large-Scale Sequencing Analyses.

BACKGROUND: Corin is a protease expressed in cardiomyocytes that plays a key role in salt handling and intravascular volume homeostasis via activation of natriuretic peptides. It is unknown if Corin loss-of-function (LOF) is causally associated with risk of coronary artery disease (CAD). METHODS: We analyzed all coding CORIN variants in an Italian case-control study of CAD. We functionally tested all 64 rare missense mutations in Western Blot and Mass Spectroscopy assays for proatrial natriuretic peptide cleavage. An expanded rare variant association analysis for Corin LOF mutations was conducted in whole exome sequencing data from 37 799 CAD cases and 212 184 controls. RESULTS: We observed LOF variants in CORIN in 8 of 1803 (0.4%) CAD cases versus 0 of 1725 controls (P, 0.007). Of 64 rare missense variants profiled, 21 (33%) demonstrated <30% of wild-type activity and were deemed damaging in the 2 functional assays for Corin activity. In a rare variant association study that aggregated rare LOF and functionally validated damaging missense variants from the Italian study, we observed no association with CAD-21 of 1803 CAD cases versus 12 of 1725 controls with adjusted odds ratio of 1.61 ([95% CI, 0.79-3.29]; P=0.17). In the expanded sequencing dataset, there was no relationship between rare LOF variants with CAD was also observed (odds ratio, 1.15 [95% CI, 0.89-1.49]; P=0.30). Consistent with the genetic analysis, we observed no relationship between circulating Corin concentrations with incident CAD events among 4744 participants of a prospective cohort study-sex-stratified hazard ratio per SD increment of 0.96 ([95% CI, 0.87-1.07], P=0.48). CONCLUSIONS: Functional testing of missense mutations improved the accuracy of rare variant association analysis. Despite compelling pathophysiology and a preliminary observation suggesting association, we observed no relationship between rare damaging variants in CORIN or circulating Corin concentrations with risk of CAD.

Single-Cell Analysis of the Normal Mouse Aorta Reveals Functionally Distinct Endothelial Cell Populations.

BACKGROUND: The cells that form the arterial wall contribute to multiple vascular diseases. The extent of cellular heterogeneity within these populations has not been fully characterized. Recent advances in single-cell RNA-sequencing make it possible to identify and characterize cellular subpopulations. METHODS: We validate a method for generating a droplet-based single-cell atlas of gene expression in a normal blood vessel. Enzymatic dissociation of 4 whole mouse aortas was followed by single-cell sequencing of >10 000 cells. RESULTS: Clustering analysis of gene expression from aortic cells identified 10 populations of cells representing each of the main arterial cell types: fibroblasts, vascular smooth muscle cells, endothelial cells (ECs), and immune cells, including monocytes, macrophages, and lymphocytes. The most significant cellular heterogeneity was seen in the 3 distinct EC populations. Gene set enrichment analysis of these EC subpopulations identified a lymphatic EC cluster and 2 other populations more specialized in lipoprotein handling, angiogenesis, and extracellular matrix production. These subpopulations persist and exhibit similar changes in gene expression in response to a Western diet. Immunofluorescence for Vcam1 and Cd36 demonstrates regional heterogeneity in EC populations throughout the aorta. CONCLUSIONS: We present a comprehensive single-cell atlas of all cells in the aorta. By integrating expression from >1900 genes per cell, we are better able to characterize cellular heterogeneity compared with conventional approaches. Gene expression signatures identify cell subpopulations with vascular disease-relevant functions.

Induced Pluripotent Stem Cells for Cardiovascular Disease Modeling and Precision Medicine: A Scientific Statement From the American Heart Association.

Induced pluripotent stem cells (iPSCs) offer an unprece-dented opportunity to study human physiology and disease at the cellular level. They also have the potential to be leveraged in the practice of precision medicine, for example, personalized drug testing. This statement comprehensively describes the provenance of iPSC lines, their use for cardiovascular disease modeling, their use for precision medicine, and strategies through which to promote their wider use for biomedical applications. Human iPSCs exhibit properties that render them uniquely qualified as model systems for studying human diseases: they are of human origin, which means they carry human genomes; they are pluripotent, which means that in principle, they can be differentiated into any of the human body's somatic cell types; and they are stem cells, which means they can be expanded from a single cell into millions or even billions of cell progeny. iPSCs offer the opportunity to study cells that are genetically matched to individual patients, and genome-editing tools allow introduction or correction of genetic variants. Initial progress has been made in using iPSCs to better understand cardiomyopathies, rhythm disorders, valvular and vascular disorders, and metabolic risk factors for ischemic heart disease. This promising work is still in its infancy. Similarly, iPSCs are only just starting to be used to identify the optimal medications to be used in patients from whom the cells were derived. This statement is intended to (1) summarize the state of the science with respect to the use of iPSCs for modeling of cardiovascular traits and disorders and for therapeutic screening; (2) identify opportunities and challenges in the use of iPSCs for disease modeling and precision medicine; and (3) outline strategies that will facilitate the use of iPSCs for biomedical applications. This statement is not intended to address the use of stem cells as regenerative therapy, such as transplantation into the body to treat ischemic heart disease or heart failure.

A Genetic Variant Associated with Five Vascular Diseases Is a Distal Regulator of Endothelin-1 Gene Expression.

Genome-wide association studies (GWASs) implicate the PHACTR1 locus (6p24) in risk for five vascular diseases, including coronary artery disease, migraine headache, cervical artery dissection, fibromuscular dysplasia, and hypertension. Through genetic fine mapping, we prioritized rs9349379, a common SNP in the third intron of the PHACTR1 gene, as the putative causal variant. Epigenomic data from human tissue revealed an enhancer signature at rs9349379 exclusively in aorta, suggesting a regulatory function for this SNP in the vasculature. CRISPR-edited stem cell-derived endothelial cells demonstrate rs9349379 regulates expression of endothelin 1 (EDN1), a gene located 600 kb upstream of PHACTR1. The known physiologic effects of EDN1 on the vasculature may explain the pattern of risk for the five associated diseases. Overall, these data illustrate the integration of genetic, phenotypic, and epigenetic analysis to identify the biologic mechanism by which a common, non-coding variant can distally regulate a gene and contribute to the pathogenesis of multiple vascular diseases.