Single-Cell Epigenomics and Functional Fine-Mapping of Atherosclerosis GWAS Loci

Rationale: Genome-wide association studies have identified hundreds of loci associated with coronary artery disease (CAD). Many of these loci are enriched in cisregulatory elements but not linked to cardiometabolic risk factors nor to candidate causal genes, complicating their functional interpretation. Objective: Single-nucleus chromatin accessibility profiling of the human atherosclerotic lesions was used to investigate cell type-specific patterns of cisregulatory elements, to understand transcription factors establishing cell identity, and to interpret CAD-relevant, noncoding genetic variation. Methods and Results: We used single-nucleus ATAC-seq (assay for transposase-accessible chromatin with sequencing) to generate DNA accessibility maps in >7000 cells derived from human atherosclerotic lesions. We identified 5 major lesional cell types including endothelial cells, smooth muscle cells, monocyte/macrophages, natural killer/T cells, and B cells and further investigated subtype characteristics of macrophages and smooth muscle cells transitioning into fibromyocytes. We demonstrated that CAD-associated genetic variants are particularly enriched in endothelial and smooth muscle cell-specific open chromatin. Using single-cell coaccessibility and cis-expression quantitative trait loci information, we prioritized putative target genes and candidate regulatory elements for approximate to 30% of all known CAD loci. Finally, we performed genome-wide experimental fine-mapping of the CAD variants identified in genome-wide association studies using epigenetic quantitative trait loci analysis in primary human aortic endothelial cells and self-transcribing active regulatory region sequencing (STARR-Seq) massively parallel reporter assay in smooth muscle cells. This analysis identified potential causal single-nucleotide polymorphisms (SNPs) and the associated target gene for over 30 CAD loci. We present several examples where the chromatin accessibility and gene expression could be assigned to one cell type predicting the cell type of action for CAD loci. Conclusions: These findings highlight the potential of applying single-nucleus ATAC-seq to human tissues in revealing relative contributions of distinct cell types to diseases and in identifying genes likely to be influenced by noncoding genome-wide association study variants.</p

Female gene networks are expressed in myofibroblast-like smooth muscle cells in vulnerable atherosclerotic plaques

Women presenting with coronary artery disease (CAD) more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. Here, we show sex-stratified gene regulatory networks (GRNs) from human carotid atherosclerotic tissue. Prioritization of these networks identified two main SMC GRNs in late-stage atherosclerosis. Single-cell RNA-sequencing mapped these GRNs to two SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like GRN was mostly expressed in plaques that were vulnerable in females. Finally, mice orthologs of the female myofibroblast-like genes showed retained expression in advanced plaques from female mice but were downregulated in male mice during atherosclerosis progression. Female atherosclerosis is driven by GRNs that promote a fibrous vulnerable plaque rich in myofibroblast-like SMCs

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques

BACKGROUND: Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. METHODS: Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe -/- mice fed a Western diet for 18 and 30 weeks. RESULTS: Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression. CONCLUSIONS: Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs

Single-Cell Gene-Regulatory Networks of Advanced Symptomatic Atherosclerosis

BACKGROUND: While our understanding of the single-cell gene expression patterns underlying the transformation of vascular cell types during the progression of atherosclerosis is rapidly improving, the clinical and pathophysiological relevance of these changes remains poorly understood. METHODS: Single-cell RNA sequencing data generated with SmartSeq2 (≈8000 genes/cell) in 16 588 single cells isolated during atherosclerosis progression in Ldlr-/-Apob100/100 mice with human-like plasma lipoproteins and from humans with asymptomatic and symptomatic carotid plaques was clustered into multiple subtypes. For clinical and pathophysiological context, the advanced-stage and symptomatic subtype clusters were integrated with 135 tissue-specific (atherosclerotic aortic wall, mammary artery, liver, skeletal muscle, and visceral and subcutaneous, fat) gene-regulatory networks (GRNs) inferred from 600 coronary artery disease patients in the STARNET (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task) study. RESULTS: Advanced stages of atherosclerosis progression and symptomatic carotid plaques were largely characterized by 3 smooth muscle cells (SMCs), and 3 macrophage subtype clusters with extracellular matrix organization/osteogenic (SMC), and M1-type proinflammatory/Trem2-high lipid-associated (macrophage) phenotypes. Integrative analysis of these 6 clusters with STARNET revealed significant enrichments of 3 arterial wall GRNs: GRN33 (macrophage), GRN39 (SMC), and GRN122 (macrophage) with major contributions to coronary artery disease heritability and strong associations with clinical scores of coronary atherosclerosis severity. The presence and pathophysiological relevance of GRN39 were verified in 5 independent RNAseq data sets obtained from the human coronary and aortic artery, and primary SMCs and by targeting its top-key drivers, FRZB and ALCAM in cultured human coronary artery SMCs. CONCLUSIONS: By identifying and integrating the most gene-rich single-cell subclusters of atherosclerosis to date with a coronary artery disease framework of GRNs, GRN39 was identified and independently validated as being critical for the transformation of contractile SMCs into an osteogenic phenotype promoting advanced, symptomatic atherosclerosis

Genetic regulation of circular RNA expression in human aortic smooth muscle cells and vascular traits

Summary: Circular RNAs (circRNAs) are a class of non-coding RNAs that have cell-type-specific expression and are relevant in cardiovascular disease. Aortic smooth muscle cells (SMCs) play a crucial role in cardiovascular disease. In this study, we employed a systems genetics approach to identify SMC circRNA transcripts and their relevance in cardiovascular traits across the genome. We quantified circRNA expression across 151 quiescent and proliferative human aortic SMCs from donors of various genetic ancestries. We identified 1,589 expressed circRNAs. Between quiescent and proliferative SMCs, we identified 173 differentially expressed circRNAs. To characterize the genetic regulation of circRNA expression, we associated the genotypes of 6.3 million single nucleotide polymorphisms (SNPs) with circRNA abundance and found 96 circRNAs that were associated with genetic loci. Three SNPs were associated with circRNA expression in proliferative SMCs but not quiescent SMCs. We identified six SNPs that had distinct association directions with circRNA isoforms from the same gene. Lastly, to identify the relevance of circRNAs in cardiovascular disease, we overlapped genetic loci associated with circRNA expression with vascular disease-related genome-wide association studies loci. We identified 14 blood pressure, one myocardial infarction, and three coronary artery disease loci, which were associated with a circRNA transcript but not an mRNA transcript. Overall, our results provide insight into the genetic basis of vascular disease traits mediated by circRNA expression

Effect of Differences in the Microbiome of Cyp17a1-Deficient Mice on Atherosclerotic Background

CYP17A1 is a cytochrome P450 enzyme that has 17-alpha-hydroxylase and C17,20-lyase activities. Cyp17a11 deficiency is associated with high body mass and visceral fat deposition in atherosclerotic female ApoE knockout (KO, d/d or −/−) mice. In the present study, we aimed to investigate the effects of diet and Cyp17a1 genotype on the gut microbiome. Female Cyp17a1 (d/d) × ApoE (d/d) (DKO) and ApoE (d/d) (controls) were fed either standard chow or a Western-type diet (WTD), and we demonstrated the effects of genetics and diet on the body mass of the mice and composition of their gut microbiome. We found a significantly lower alpha diversity after accounting for the ecological network structure in DKO mice and WTD-fed mice compared with chow-fed ApoE(d/d). Furthermore, we found a strong significant positive association of the Firmicutes vs. Bacteroidota ratio with body mass and the circulating total cholesterol and triglyceride concentrations of the mice when feeding the WTD, independent of the Cyp17a1 genotype. Further pathway enrichment and network analyses revealed a substantial effect of Cyp17a1 genotype on associated cardiovascular and obesity-related pathways involving aspartate and L-arginine. Future studies are required to validate these findings and further investigate the role of aspartate/L-arginine pathways in the obesity and body fat distribution in our mouse model

Media 1: Ultrahigh-resolution, high-speed spectral domain optical coherence phase microscopy

Originally published in Optics Letters on 01 January 2014 (ol-39-1-45