Integrative analysis of single-cell transcriptomics and genetic associations identify cell states associated with vascular disease

Abstract

Background

Vascular diseases are accompanied by alterations in cellular phenotypes which underlie disease pathogenesis, with single-cell technologies aiding in the discovery of cellular heterogeneity among endothelial cell (EC) and vascular smooth muscle cell (VSMC) populations. In atherosclerotic disease, VSMCs are hypothesized to transition between contractile and synthetic states; however, the specific vascular subpopulations and intermediate cell states responsible for early vascular dysfunction remain unclear.

Methods

We integrated newly generated and published single-nuclear RNA-sequencing (snRNA-seq) datasets to analyze normal (n = 7), aneurysmal (n = 9), and atherosclerotic (n = 2) flash-frozen human ascending thoracic aortas. Cell types and subtypes were defined using both top marker genes and canonical gene markers. Disease enrichment and relevant cell types were identified using newly developed computational tools to integrate GWAS data from multiple vascular disease-relevant studies with the single nuclei aortic expression profiles.

Results

Nuclear dissociation and snRNA-seq identified ten distinct transcriptomic clusters from the integrated analysis representing all major vascular cell populations. Three distinct VSMC populations emerged that exhibited differential expression of extracellular matrix, contractile and pro-proliferative genes. Aneurysmal specimens were enriched for one fibroblast and one VSMC subpopulation compared to healthy tissue. RNA-trajectory analysis inferred a phenotypic continuum of gene expression between VSMC A and VSMC B or C and between two of the identified fibroblast types. VSMCs and Fibroblast C exhibited the greatest cell type-specific enrichment of genes mapped to GWAS loci for coronary artery disease (CAD), blood pressure, and migraine. Cell type-specific enrichment scores were more robust among the transcriptional profiles from non-diseased vascular tissue.

Conclusions

Our use of single-cell isolation and new computational methods prioritizes the cell types that most contribute to vascular disease pathogenesis. Specifically, tissue dissociation and single-nuclear transcriptomics better represent all vascular cell types, from which we demonstrate enrichment of pro-proliferative VSMCs in TAA and further implicate phenotypic switching as a likely pathologic mechanism. Integrated analysis of cell-specific gene expression and vascular disease GWAS data implicate genes and pathways associated with fibroblast and VSMC cell-state transitions.