Maria Elishaev1
Chani J. Hodonsky2
Saikat Kumar B. Ghosh3
Aloke V. Finn3Moritz von Scheidt4,5
Ying Wang1*- 1Department of Pathology and Laboratory Medicine, Center for Heart Lung Innovation, University of British Columbia, Vancouver, BC, Canada
- 2Center for Public Health Genomics, Department of Public Health Sciences, University of Virginia, Charlottesville, VA, United States
- 3Cardiovascular Pathology Institute, Gaithersburg, MD, United States
- 4Department of Cardiology, Deutsches Herzzentrum München, Technische Universität München, Munich, Germany
- 5Deutsches Zentrum für Herz-Kreislauf-Forschung, Partner Site Munich Heart Alliance, Munich, Germany
Over the last few years, new high-throughput biotechnologies and bioinformatic methods are revolutionizing our way of deep profiling tissue specimens at the molecular levels. These recent innovations provide opportunities to advance our understanding of atherosclerosis using human lesions aborted during autopsies and cardiac surgeries. Studies on human lesions have been focusing on understanding the relationship between molecules in the lesions with tissue morphology, genetic risk of atherosclerosis, and future adverse cardiovascular events. This review will highlight ways to utilize human atherosclerotic lesions in translational research by work from large cardiovascular biobanks to tissue registries. We will also discuss the opportunities and challenges of working with human atherosclerotic lesions in the era of next-generation sequencing.
Introduction
Ischemic cardiovascular events, including heart attack and strokes, is the leading cause of mortality and morbidity worldwide (1). Atherosclerosis, the build-up of lesion cells on the wall of the blood vessels, is a chronic process underlying most ischemic cardiovascular events. Studies of human biospecimens have played an indispensable part in understanding the pathophysiology of atherosclerosis because neither cell culture nor animal models can recapitulate the complex components and structure of advanced human atherosclerotic lesions (2, 3). Before next-generation sequencing was invented, human biospecimen studies focused on the morphology and a few lesion components such as collagen, foam cells, and smooth muscle cells. Even with limited dimension, the classic histology methods have set the widely-accepted lesion classification standards and found the association between lesion morphology and ischemic events (4). The era of next-generation sequencing brought opportunities to explore the molecular features of atherosclerotic lesions in depth and in situ. Hence, researchers have multiple options to apply the classic and new technologies to human biospecimens. The question is how to match the research scope and approach with the human biospecimens available in large biobanks or small-scale tissue registries. This review will provide a few examples to showcase current advances in utilizing human atherosclerotic lesions, share our views of challenges in knowledge translation and vision of future needs to advance research on human biospecimens in the area of atherosclerosis.
Current Advances in Human Biospecimen Studies
Morphological Insights
Histology assessment is a routine in most tissue biobanks. Traditional hematoxylin-eosin and Movat's Pentachrome stainings visualize the structure and basic components of the lesions. These widely adopted methods revealed not only the association between morphology and cardiovascular events, but also risk factors that contribute to disease progression. Biopsy studies found that the morphological features of lesions underneath cardiac thrombosis fall into three categories: rupture, erosion, and calcified nodule (5), indicating the connection between morphology and risk of ischemic cardiovascular events. Following this notion, Virmani et al. developed a comprehensive assessment criteria using morphology to define the trajectory of lesion progression (5). Correlation studies can further infer the pathogenic drive of lesion progression by connecting the traditional risk factors of atherosclerosis with morphology. Burke et al. found that patients having plaque rupture also had higher cholesterol levels as recorded in their postmortem toxicological tests (6). It suggests that the traditional risk factor, cholesterol levels, drives lesion progression toward a rupture-prone structure and lipid-lowering drugs will benefit patients having this type of lesions. In addition to studying lesions at late stages, Nakashima et al. observed that vascular beds that develop diffuse intimal thickening in early childhood are also the “hot spots” for atherogenesis in later life (7). Large biobanks have a statistically powerful amount of samples to explore how lesion morphology is affected by genetic risk factors in patients of different ethnicities. CVPath Institute has the world's largest and most comprehensive collection of heart samples from more than 8,500 sudden coronary death cases. Their studies have shown a higher risk of sudden coronary death in the African American population compared to the Caucasian population (8). Guo et al. explored the genetic reasons behind this phenomenon and found that more African Americans carry a single nucleotide polymorphism rs7136716 than the Caucasian population and this genetic variant is correlated with increased expression of CD163 macrophages in ruptured lesions (9). Then they determined that CD163+ macrophages contribute to intraplaque microvessels and inflammation, remodeling lesion structure to be rupture-prone. This data suggests that rs7136716 is a genetic risk variant particularly enriched in patients of an African ancestry. Hence, the integration of pathology and genetics from large biobank studies can provide new insights into the discovery of biomarkers and personalized medicine.
Molecular Phenotyping
There has been an explosive increase in atherosclerosis research using single-cell RNA sequencing (scRNA-seq). These studies have revealed the vast heterogeneity in the cell components of atherosclerotic lesions, including more than six phenotypes of SMCs (10–12) and at least three macrophage subsets (13, 14) among eleven distinct leukocyte populations (15). Most previous studies applied next-generation sequencing to animal models, in which disease stages can be controlled, and cell lineage tracing is available. These advantages made them better tools for understanding the trajectory of phenotypic changes and trans-differentiation of lesion cells in the early disease stage (10, 16). Applying next-generation sequencing to human lesions has its own technical challenges and limitations, discussed in several recent reviews (17–19). Nevertheless, single-cell sequencing of human lesions will distinguish cell phenotypes and signaling pathways in categorized patient cohorts, the differences in which can be used to predict therapeutic targets and treatment outcomes. Slenders et al. projected GWAS loci into scRNA-seq data of human carotid lesions and defined cell-specific risk genes that can be translated to therapeutic targets (20). Fernandez et al. applied several single-cell sequencing technologies to human carotid lesions. They found that differences in the interleukin-1β (IL-1β) signaling pathway between symptomatic and asymptomatic patients may lead to diverse treatment outcomes in the CANTOS trial, which aims to block IL-1β to reduce the risks of cardiovascular events (21). One limitation of using scRNA-seq to predict cellular function is that it ignores all the post-transcriptional regulation of gene expression. Quantitative proteomics uses technologies such as mass spectrometry to characterize proteins, molecules that eventually execute gene function. Since traditional mass spectrometry usually requires a few milligrams of input material, most previous studies focus on cultured cells instead of digested lesions to obtain the protein atlas of a single cell type. Okui et al. applied mass spectrometry to cultured human coronary artery SMCs and found that carnitine O-octanoyltransferase increases during osteogenic transition of SMCs, which will lead to calcification during atherogenesis (22). The Athero-Express study used proteomics to analyze the protein features of human carotid lesions and found that osteopontin was strongly associated with cardiovascular events during the follow-up post endarterectomy surgery. suggesting that pathologists can use osteopontin as a biomarker to predict patients' future risk of cardiovascular events (23). The Athero-Express study collects carotid lesions from endarterectomy surgeries and follows up with the clinical records of donors to determine the relationship between lesion composition and future adverse cardiovascular events (23). It embraces the idea that atherosclerosis is a systemic disease and biopsy tests of surgically aborted lesions will inform the progression of lesions in other vascular beds (24). The Athero-Express study has recently applied state-of-the-art scRNA-seq technology to their biobank to advance our understanding of the transcriptome landscape of human lesions (25), and more specific biomarkers are expected with accumulated data of clinical follow-up. Hence, when the molecular features of human lesions are linked to clinical follow-up of the specimen donors, it will support the discovery of biomarkers.
Proof of Concept
Compared to large biobanks, sample availability in small-scale tissue registries is often too low to support genetic studies such as investigating the polygenic risk factors. Without follow-up of clinical outcomes, the archived tissue samples are not ideal for biomarker studies. However, these tissues can be readily used to prove the concepts of basic sciences: “reality check” that observations in animal studies or cultured cells will apply to human atherogenesis. It is now well-known that vascular smooth muscle cells (SMCs) undergo phenotypic changes during atherogenesis and most of them do not express lineage markers such as myosin heavy chain 11 and smooth muscle alpha-actin, but it was not till the invention of SMC-lineage tracing mice when researchers realized about this (26). To prove the presence of de-differentiated SMCs in human atherogenesis, Gomez et al. developed a staining method to visualize cells of SMC origin using SMC-lineage tracing mice and then applied the technique to human carotid lesions (n = 5) (27). Clinical information can be completely detached from the biospecimen studies in such proof of concept research and small sample size is acceptable. Working with archived samples shortens the time of ethical approval for researchers. Proof of concept studies allow archived specimens to be repurposed and maximized in the tissue registries. For example, the Cardiovascular Tissue Registry at the Centre for Heart Lung Innovation at the University of British Columbia collects hearts donated by heart transplant patients, and the myocardium has been utilized to support studies of cardiac allograft vasculopathy (28). In parallel, Allahverdian et al. characterized the atherosclerotic coronary arteries from these hearts and discovered that an underestimated amount of foam cells are derived from SMCs instead of macrophages (29). This study, along with growing evidence from GWAS of coronary artery disease (30–32), has changed the traditional dogma that atherosclerosis is mainly a macrophage-driven disease.
Current Gaps in Knowledge Translation
Despite the current advances in new technologies, guidelines of research design (19) and customized protocols (33) to maximize the utilization of human biospecimen, translation of biobank research to biomarkers and therapeutic targets is still limited. Current gaps in knowledge translation include the lack of tailored bioinformatics tools to interpret “Omics” data, “Omics” in the tissue context, and connections of traits in the lesions with blood-borne biomarkers, which is more applicable for clinical testing.
Bioinformatic Tools
As more human biospecimen studies have embraced scRNA-seq technology, bioinformatics tools that are primarily R- and python-based have been developed, including Seurat (34), Signac (35), archR (36), singleR (37), scCATH (38), and Garnett (39). Additionally, the less robust single-cell sequencing data can intersect with large bulk sequencing datasets, such as the Human Cell Atlas (40) through Cibersort (41) and MuSiC (42), to estimate the portion of specific cell types in the bulk-seq data. However, only a few tools incorporate established pipelines into webapps or standalone applications to allow scientists without programming experience to analyse their data. Most bioinformatics tools are also incompatible with integrating datasets published in various formats. This problem is accompanied by bioethical considerations of what clinical information affiliated with the human biospecimen can be shared and how to oversee the appropriate use of shared datasets (43–46). Furthermore, myriad analysis options are available: over 1,000 tools are developed for scRNA-seq alone (47), with limited benchmarking for atherosclerotic lesions. The nature of cells in the atherosclerotic lesion: cell plasticity and dynamic phenotypic changes during atherogenesis (10, 48), is especially challenging for bioinformatics tools that rely on: (1) reference datasets to define a cell type, and (2) the persistent presence of a phenotype in transition to trace the trajectory of atherogenic changes (49–52).
Lack of Spatial “Omics” Data
Gene ontology analyses of scRNA-seq data predict the function of each cell subset based on signaling pathways and biological processes involving genes specific to that subset. However, cells are unevenly distributed among the atherosclerotic lesions. In the traditional diagram of a fibroatheroma, endothelial cells are located at the luminal side with SMCs underneath and macrophages form foam cells in the intima and shoulder region. After the discovery of cell subsets and phenotypic transition of lesion cells, now we know that the fibrous cap is derived from SMCs and endothelial cells that have undergone endothelial-to-mesenchymal transition (16). CD68+ “macrophages” in the deep intima are made of leukocytes and SMCs (29). The stability of lesion structure is affected by the distribution of these cell subsets and their interactions with the others. Within the several types of SMCs in atherosclerotic lesions, some are fibroblast-like, potentially forming a protective fibrous cap (53), whereas pro-inflammatory macrophages and SMCs may enlarge the necrotic core by altering the function of each other (54). These assumptions can be validated only when we see fibroblast-like SMCs on the fibrous cap of stable lesions and more impaired macrophages are close to the pro-inflammatory SMCs near the necrotic core. In the absence of lesion context, it can be misleading to predict lesion progression purely based on signaling pathways and biological processes represented by gene expression. The same mitogenic signaling pathway turned on by IL-1β in SMCs may cause pathogenic expansion of the atherosclerotic lesions (55) when these SMCs are located in the intima, but it may be athero-protective when these cells are in the fibrous cap. Blocking IL-1β in animals with established fibroatheroma resulted in lesion destabilization, suggesting that at least a subset of SMCs executed IL-1β's mitogenic effects in an athero-protective manner by investing into the fibrous cap (56). Hence, it is critical to consider spatial “Omics” data when translating single-cell “Omics” data to therapeutic targets and biomarkers.
Connection to Traits in the Blood
Tissue biopsy examinations are well-accepted for phenotyping tumors but are not practical in the clinical practice of atherosclerotic disease where not all patients require surgical treatment. Connecting the traits of atherosclerotic lesions to those in the blood will help the discovery of blood-borne biomarkers that can be easily implemented in clinical laboratories. Phenotyping blood cells has the potential to reveal cellular biomarkers in the circulation (57). Hamers et al. explored protein expression of circulating monocytes using a 39-plex CyTOF panel and found a positive correlation between Slan+CXCR6+ monocyte population in the blood and the severity of atherosclerotic disease estimated by medical imaging, suggesting that this monocyte population could be a biomarker for disease progression (57). More importantly, studying biopsy and blood samples from the same donor will determine whether molecules in the blood mirror biological activities in the lesion. The interleukin-6 (IL-6) signaling pathway is a crucial mediator of inflammation in the lesions (58), and elevated plasma IL-6 in patients with acute coronary syndrome is currently the most powerful predictor of increased mortality in both short and long terms (59, 60). To explore its utilization as a biomarker for risk of cardiovascular events, the Biobank of Karolinska Endarterectomies study found that plasma IL-6 is positively correlated with those in the carotid lesions (61), but no significant differences were observed between symptomatic and asymptomatic patients (62). Hence, plasma IL-6 does not precisely reflect lesion progression and risk of adverse cardiovascular events. Until now, we have not found clinically applicable biomarkers to mirror biological activities in atherosclerotic lesions. Clonal hematopoiesis of indeterminate potential (CHIP), a phenomenon by which blood precursor cells in bone marrow obtain mutations during aging, is associated with a higher risk of atherosclerosis (63). In this original study, the authors found more than one hundred CHIP mutations in peripheral blood cells but did not have access to any molecular information about the atherosclerotic lesions. This study left the knowledge gap on how specific CHIP mutations will reflect lesion progression and can be promising biomarkers for risk prediction.
Opportunities and Challenges in Future Studies
Collaborations are the key to advancing knowledge translation of human biospecimen studies. Publicly available “Omics” data from multiple studies can be pooled together to enhance the diversity of the studied population given that most data are based on European descendants. For researchers who already have the datasets, meta-analysis will combine different studies, as applied in the Coronary Artery Disease Genome-wide Replication and Meta-analysis plus The Coronary Artery Disease Genetics consortium. It is critical for researchers working on human atherosclerotic lesions to assess the potential pitfalls of applying a bioinformatic pipeline tested on other tissues to work on atherosclerotic lesions, which requires communication between the dry and wet labs to reflect on the principles of the pipelines and biological features of the tissues. For example, a spectrum of early to late staged atherosclerotic lesions need to be incorporated into the research design to answer questions about cell trajectories. For bioinformaticians, a consortium that integrates RNA-seq datasets has not been established yet, but initiatives have been made for sharing datasets in a user-friendly way. PlaqView is the first platform that empowers researchers who do not have access to human atherosclerotic lesions or who do not have bioinformatics expertise to re-analyze published scRNA-seq data (64). More biobanks have started to investigate blood and lesion samples from the same donor side-by-side to develop blood-borne biomarkers. The recently established Munich Cardiovascular Studies Biobank at the German Heart Center Munich has collected more than 800 pairs of blood and lesion samples to connect molecular traits in atherosclerotic lesions to those in blood. Small-scale tissue registries have the opportunity to collaborate with biobanks that have access to blood samples and patients' clinical follow-up records to validate molecular targets or biomarkers found in tissue studies. Spatial gene expression (Table 1) and multiplex imaging technologies (Table 2) have been commercialized to optimize the use of archived tissue blocks by measuring gene and protein expression in a high-throughput fashion. Roadblocks to applying them to human atherosclerotic lesions are mainly technical. Regarding histology integrity, lesion sections can easily fold and tear during sample preparation due to the natural curvature of the lumen and the presence of calcification and large necrotic cores. Regarding quality control of the RNA, atherosclerotic lesions, especially the highly necrotic ones, have low cellularity compared to the more frequently reported tumor, brain, and myocardium tissues. Based on our experience and previous research (65), a large part of one lesion (at least 100 μm long) is required to extract an adequate amount of RNA for quality assessment: 5 ng for the 2100 Agilent bioanalyzer (RNA 6000 Nano kit) or RT-PCR to amplify housekeeping genes. Moreover, Formalin-Fixed Paraffin-Embedded is the standard archiving format in biobanks to preserve lesion morphology. RNA stability remains unchanged for up to 10 years in this format (65). Previous research found that RNA is unevenly degraded within the same lesion section (66). Hence, we think that in situ assessment of RNA quality (67) is more suitable than extracting RNA from the whole section to select samples for spatial gene expression. Such a method can advise which lesion region has high quality input material to generate sequencing reads of high fidelity. Regarding the interpretation of data, cell segmentation has been the most problematic step for all the spatial biology tools. So far, most software and analysis pipelines assign the signals to a cell nucleus nearby based on algorithms trained in cancer tissues (68–70), assuming that cells in other tissues have similar size and shape. This may lead to inaccurate cell segmentation in lesion sections when the spindle-shaped 200 μm long SMCs, spherical lymphocytes of 7–10 μm in diameter, and foam cells with various sizes are all in close proximity. For spatial gene expression, transcriptomes from multiple cells can be mixed in one captured region and devolution of cell types requires reference scRNA-seq data (71, 72). Unlike working with dissociated cells, multiple sections are required to assess the entire lesion using spatial biology technologies, significantly increasing the cost. Removing these roadblocks requires efforts from both the academic and the technology industries.
Table 1. Spatial gene expression technologies.
Table 2. Spatial protein expression technologies.
Conclusion
Human atherosclerotic lesions carry valuable morphological and molecular information to decipher the mechanism of atherogenesis and reveal therapeutic targets and biomarkers for patients. While large biobanks with adequate sample numbers and clinical data can perform genetic and biomarker studies, a smaller sample pool in the tissue registry also plays a vital role in translating basic sciences to a human disease scenario. New biotechnologies, along with bioinformatic tools to process the data, are modernizing biobank-based research in atherosclerosis, providing opportunities to maximize the utilization of human biospecimen. However, these new tools validated in other tissues are not one-size-fits-all, given the complex cell components in human atherosclerotic lesions. Collaborations from researchers in cardiovascular disease, bioinformaticians, and technology developers are essential for benchmarking and customization of available tools to address the unique challenges in future studies of human atherosclerotic lesions.
Author Contributions
YW designed the topic of this review article and wrote 70% of the manuscript with ME. ME performed literature review and designed the tables. SG, CH, AF, and MS wrote the rest of the manuscript and contributed to the design of review topic. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Canadian Institutes of Health Research (GR022678 to YW) and the Fondation Leducq (PlaqOmics 18CVD02 to AF and its Junior Investigator Award to CH, SG, MS, and YW).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors wish to acknowledge Gurpreet Singhera and the Cardiovascular Tissue Registry at Centre for Heart Lung Innovation for critical input on the manuscript.
References
1. World Health Organization. The Top 10 Causes of Death. (2020). Available online at: https://www.who.int/news-room/fact-sheets/detail/the-top-10-causes-of-death (accessed April 28, 2022].
2. Silvestre-Roig C, de Winther MP. Weber C., Daemen MJ, Lutgens E, Soehnlein O. Atherosclerotic plaque destabilization: mechanisms, models, therapeutic strategies. Circ Res. (2014) 114:214–26. doi: 10.1161/CIRCRESAHA.114.302355
3. Hartwig H, Silvestre-Roig C, Hendrikse J, Beckers L, Paulin N, Van der Heiden K, et al. Atherosclerotic plaque destabilization in mice: a comparative study. PLoS ONE. (2015) 10:e0141019. doi: 10.1371/journal.pone.0141019
4. Sakakura K, Nakano M, Otsuka F, Ladich E, Kolodgie FD, Virmani R. Pathophysiology of atherosclerosis plaque progression. Heart Lung Circ. (2013) 22:399–411. doi: 10.1016/j.hlc.2013.03.001
5. Virmani R, Kolodgie FD, Burke AP, Farb A, Schwartz SM. Lessons from sudden coronary death: a comprehensive morphological classification scheme for atherosclerotic lesions. Arterioscler Thromb Vasc Biol. (2000) 20:1262–75. doi: 10.1161/01.ATV.20.5.1262
6. Burke AP, Farb A, Malcom GT, Liang YH, Smialek J, Virmani R. Coronary risk factors and plaque morphology in men with coronary disease who died suddenly. N Engl J Med. (1997) 336:1276–82. doi: 10.1056/NEJM199705013361802
7. Nakashima Y, Chen YX, Kinukawa N, Sueishi K. Distributions of diffuse intimal thickening in human arteries: preferential expression in atherosclerosis-prone arteries from an early age. Virchows Arch. (2002) 441:279–88. doi: 10.1007/s00428-002-0605-1
8. Burke AP, Farb A, Pestaner J, Malcom GT, Zieske A. Kutys R, et al. Traditional risk factors and the incidence of sudden coronary death with and without coronary thrombosis in blacks. Circulation. (2002) 105:419–24. doi: 10.1161/hc0402.102952
9. Guo L, Akahori H, Harari E, Smith SL, Polavarapu R. Karmali V, et al. CD163+ macrophages promote angiogenesis and vascular permeability accompanied by inflammation in atherosclerosis. J Clin Invest. (2018) 128:1106–24. doi: 10.1172/JCI93025
10. Pan H, Xue C, Auerbach BJ, Fan J, Bashore AC, Cui J, et al. Single-cell genomics reveals a novel cell state during smooth muscle cell phenotypic switching and potential therapeutic targets for atherosclerosis in mouse and human. Circulation. (2020) 142:2060–75. doi: 10.1161/CIRCULATIONAHA.120.048378
11. Yap C, Mieremet A, de Vries CJM, Micha D, de Waard V. Six shades of vascular smooth muscle cells illuminated by KLF4 (Kruppel-Like Factor 4). Arterioscler Thromb Vasc Biol. (2021) 41:2693–707. doi: 10.1161/ATVBAHA.121.316600
12. Alencar GF, Owsiany KM, Karnewar S, Sukhavasi K, Mocci G, Nguyen AT, et al. Stem cell pluripotency genes Klf4 and Oct4 regulate complex SMC phenotypic changes critical in late-stage atherosclerotic lesion pathogenesis. Circulation. (2020) 142:2045–59. doi: 10.1161/CIRCULATIONAHA.120.046672
13. Cochain C, Vafadarnejad E, Arampatzi P, Pelisek J, Winkels H, Ley K, et al. Single-cell RNA-Seq reveals the transcriptional landscape and heterogeneity of aortic macrophages in murine. Atheroscler Circ Res. (2018) 122:1661–74. doi: 10.1161/CIRCRESAHA.117.312509
14. Zernecke A, Winkels H, Cochain C, Williams JW, Wolf D, Soehnlein O, et al. Meta-analysis of leukocyte diversity in atherosclerotic mouse. Aortas Circ Res. (2020) 127:402–26. doi: 10.1161/CIRCRESAHA.120.316903
15. Winkels H, Ehinger E, Vassallo M, Buscher K, Dinh HQ, Kobiyama K, et al. Atlas of the immune cell repertoire in mouse atherosclerosis defined by single-cell RNA-sequencing and mass cytometry. Circ Res. (2018) 122:1675–88. doi: 10.1161/CIRCRESAHA.117.312513
16. Newman AAC. Serbulea V, Baylis RA, Shankman LS, Bradley X, Alencar GF, et al. Multiple cell types contribute to the atherosclerotic lesion fibrous cap by PDGFRbeta and bioenergetic mechanisms. Nat Metab. (2021) 3:166–81. doi: 10.1038/s42255-020-00338-8
17. Iqbal F, Lupieri A, Aikawa M, Aikawa E. Harnessing single-cell RNA sequencing to better understand how diseased cells behave the way they do in cardiovascular disease. Arterioscler Thromb Vasc Biol. (2021) 41:585–600. doi: 10.1161/ATVBAHA.120.314776
18. Slenders L, Tessels DE, van der Laan SW, Pasterkamp G, Mokry M. The applications of single-cell RNA sequencing in atherosclerotic disease. Front Cardiovasc Med. (2022) 9:826103. doi: 10.3389/fcvm.2022.826103
19. Williams JW, Winkels H, Durant CP, Zaitsev K, Ghosheh Y, Ley K. Single cell RNA sequencing in atherosclerosis research. Circ Res. (2020) 126:1112–26. doi: 10.1161/CIRCRESAHA.119.315940
20. Slenders L, Landsmeer LPL, Cui K, Depuydt MAC, Verwer M, Mekke J, et al. Intersecting single-cell transcriptomics and genome-wide association studies identifies crucial cell populations and candidate genes for atherosclerosis. Eur Heart J Open. (2022) 2:oeab043. doi: 10.1093/ehjopen/oeab043
21. Fernandez DM, Rahman AH, Fernandez NF, Chudnovskiy A, Amir ED. Amadori L, et al. Single-cell immune landscape of human atherosclerotic plaques. Nat Med. (2019) 25:1576–88. doi: 10.1038/s41591-019-0590-4
22. Okui T, Iwashita M, Rogers MA, Halu A, Atkins SK, Kuraoka S, et al. CROT (Carnitine O-Octanoyltransferase) is a novel contributing factor in vascular calcification via promoting fatty acid metabolism and mitochondrial dysfunction. Arterioscler Thromb Vasc Biol. (2021) 41:755–68. doi: 10.1161/ATVBAHA.120.315007
23. de Kleijn DP, Moll FL, Hellings WE, Ozsarlak-Sozer G, de Bruin P, Doevendans PA, et al. Local atherosclerotic plaques are a source of prognostic biomarkers for adverse cardiovascular events. Arterioscler Thromb Vasc Biol. (2010) 30:612–9. doi: 10.1161/ATVBAHA.109.194944
24. Hellings WE, Moll FL, de Kleijn DP, Pasterkamp G. 10-years experience with the Athero-Express study. Cardiovasc Diagn Ther. (2012) 2:63–73. doi: 10.3978/j.issn.2223-3652.2012.02.01
25. Depuydt MAC, Prange KHM, Slenders L, Ord T, Elbersen D, Boltjes A, et al. Microanatomy of the human atherosclerotic plaque by single-cell transcriptomics. Circ Res. (2020) 127:1437–55. doi: 10.1161/CIRCRESAHA.120.316770
26. Feil S, Fehrenbacher B, Lukowski R, Essmann F, Schulze-Osthoff K, Schaller M, et al. Transdifferentiation of vascular smooth muscle cells to macrophage-like cells during atherogenesis. Circ Res. (2014) 115:662–7. doi: 10.1161/CIRCRESAHA.115.304634
27. Gomez D, Shankman LS, Nguyen AT, Owens GK. Detection of histone modifications at specific gene loci in single cells in histological sections. Nat Methods. (2013) 10:171–7. doi: 10.1038/nmeth.2332
28. Wong BW, Rahmani M, Luo Z, Yanagawa B, Wong D, Luo H, et al. Vascular endothelial growth factor increases human cardiac microvascular endothelial cell permeability to low-density lipoproteins. J Heart Lung Transplant. (2009) 28:950–7. doi: 10.1016/j.healun.2009.05.005
29. Allahverdian S, Chehroudi AC, McManus BM, Abraham T, Francis GA. Contribution of intimal smooth muscle cells to cholesterol accumulation and macrophage-like cells in human atherosclerosis. Circulation. (2014) 129:1551–9. doi: 10.1161/CIRCULATIONAHA.113.005015
30. Wong D, Turner AW, Miller CL. Genetic insights into smooth muscle cell contributions to coronary artery disease. Arterioscler Thromb Vasc Biol. (2019) 39:1006–17. doi: 10.1161/ATVBAHA.119.312141
31. Erdmann J, Kessler T, Munoz Venegas L, Schunkert H. A decade of genome-wide association studies for coronary artery disease: the challenges ahead. Cardiovasc Res. (2018) 114:1241–57. doi: 10.1093/cvr/cvy084
32. Hao K, Ermel R, Sukhavasi K, Cheng H, Ma L, Li L, et al. Integrative prioritization of causal genes for coronary artery disease. Circ Genom Precis Med. (2022) 15:e003365. doi: 10.1161/CIRCGEN.121.003365
33. Corces MR, Trevino AE, Hamilton EG, Greenside PG, Sinnott-Armstrong NA, Vesuna S, et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat Methods. (2017) 14:959–62. doi: 10.1038/nmeth.4396
34. Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. (2015) 33:495–502. doi: 10.1038/nbt.3192
35. Stuart T, Srivastava A, Madad S, Lareau CA, Satija R. Single-cell chromatin state analysis with Signac. Nat Methods. (2021) 18:1333–41. doi: 10.1038/s41592-021-01282-5
36. Granja JM, Corces MR, Pierce SE, Bagdatli ST, Choudhry H, Chang HY, et al. ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis. Nat Genet. (2021) 53:403–11. doi: 10.1038/s41588-021-00790-6
37. Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol. (2019) 20:163–72. doi: 10.1038/s41590-018-0276-y
38. Shao X, Liao J, Lu X, Xue R, Ai N, Fan X. scCATCH: automatic annotation on cell types of clusters from single-cell RNA sequencing data. iScience. (2020) 23:100882. doi: 10.1016/j.isci.2020.100882
39. Pliner HA, Shendure J, Trapnell C. Supervised classification enables rapid annotation of cell atlases. Nat Methods. (2019) 16:983–6. doi: 10.1038/s41592-019-0535-3
40. Regev A, Teichmann SA, Lander ES, Amit I, Benoist C, Birney E, et al. The human cell atlas. Elife. (2017) 6. doi: 10.7554/eLife.27041
41. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W. Xu, Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
42. Wang X, Park J, Susztak K, Zhang NR Li M. Bulk tissue cell type deconvolution with multi-subject single-cell expression reference. Nat Commun. (2019) 10:380. doi: 10.1038/s41467-018-08023-x
43. Powell K. The broken promise that undermines human genome research. Nature. (2021) 590:198–201. doi: 10.1038/d41586-021-00331-5
44. Tsosie KS, Fox K, Yracheta JM. Genomics data: the broken promise is to Indigenous people. Nature. (2021) 591:529. doi: 10.1038/d41586-021-00758-w
45. Fox K. The illusion of inclusion - the “All of Us” Research Program and Indigenous Peoples' DNA. N Engl J Med. (2020) 383:411–3. doi: 10.1056/NEJMp1915987
46. Phillips M. International data-sharing norms: from the OECD to the General Data Protection Regulation (GDPR). Hum Genet. (2018) 137:575–82. doi: 10.1007/s00439-018-1919-7
47. Zappia L, Theis FJ. Over 1000 tools reveal trends in the single-cell RNA-seq analysis landscape. Genome Biol. (2021) 22:301. doi: 10.1186/s13059-021-02519-4
48. Wang Y, Gao H, Wang F, Ye Z, Mokry M, Turner AW, et al. Dynamic changes in chromatin accessibility are associated with the atherogenic transitioning of vascular smooth muscle cells. Cardiovasc Res. (2021). cvab347. doi: 10.1093/cvr/cvab347
49. Wagner DE, Klein AM. Lineage tracing meets single-cell omics: opportunities and challenges. Nat Rev Genet. (2020) 21:410–27. doi: 10.1038/s41576-020-0223-2
50. Zaccaria S, Raphael BJ. Characterizing allele- and haplotype-specific copy numbers in single cells with CHISEL. Nat Biotechnol. (2021) 39:207–14. doi: 10.1038/s41587-020-0661-6
51. Cao J, Spielmann M, Qiu X, Huang X, Ibrahim DM, Hill AJ, et al. The single-cell transcriptional landscape of mammalian organogenesis. Nature. (2019) 566:496–502. doi: 10.1038/s41586-019-0969-x
52. Bergen V, Lange M, Peidli S, Wolf FA, Theis FJ. Generalizing RNA velocity to transient cell states through dynamical modeling. Nat Biotechnol. (2020) 38:1408–14. doi: 10.1038/s41587-020-0591-3
53. Wirka RC, Wagh D, Paik DT, Pjanic M, Nguyen T, Miller CL, et al. Atheroprotective roles of smooth muscle cell phenotypic modulation and the TCF21 disease gene as revealed by single-cell analysis. Nat Med. (2019) 25:1280–9. doi: 10.1038/s41591-019-0512-5
54. Wang Y, Nanda V, Direnzo D, Ye J, Xiao S, Kojima Y, et al. Clonally expanding smooth muscle cells promote atherosclerosis by escaping efferocytosis and activating the complement cascade. Proc Natl Acad Sci USA. (2020) 117:15818–26. doi: 10.1073/pnas.2006348117
55. Libby P, Warner SJ, Friedman GB. Interleukin 1: a mitogen for human vascular smooth muscle cells that induces the release of growth-inhibitory prostanoids. J Clin Invest. (1988) 81:487–98. doi: 10.1172/JCI113346
56. Gomez D, Baylis RA, Durgin BG, Newman AAC, Alencar GF, Mahan S, et al. Interleukin-1beta has atheroprotective effects in advanced atherosclerotic lesions of mice. Nat Med. (2018) 24:1418–29. doi: 10.1038/s41591-018-0124-5
57. Hamers AAJ, Dinh HQ, Thomas GD, Marcovecchio P, Blatchley A, Nakao CS, et al. Human monocyte heterogeneity as revealed by high-dimensional mass cytometry. Arterioscler Thromb Vasc Biol. (2019) 39:25–36. doi: 10.1161/ATVBAHA.118.311022
58. Tousoulis D, Oikonomou E, Economou EK, Crea F, Kaski JC. Inflammatory cytokines in atherosclerosis: current therapeutic approaches. Eur Heart J. (2016) 37:1723–32. doi: 10.1093/eurheartj/ehv759
59. Andrie RP, Becher UM, Frommold R, Tiyerili V, Schrickel JW, Nickenig G, et al. Interleukin-6 is the strongest predictor of 30-day mortality in patients with cardiogenic shock due to myocardial infarction. Crit Care. (2012) 16:R152. doi: 10.1186/cc11467
60. Gager GM, Biesinger B, Hofer F, Winter MP, Hengstenberg C, Jilma B, et al. Interleukin-6 level is a powerful predictor of long-term cardiovascular mortality in patients with acute coronary syndrome. Vascul Pharmacol. (2020) 135:106806. doi: 10.1016/j.vph.2020.106806
61. Ziegler L, Lundqvist J, Dreij K, Wallen H, de Faire U, Paulsson-Berne G, et al. Expression of Interleukin 6 signaling receptors in carotid atherosclerosis. Vasc Med. (2021) 26:3–10. doi: 10.1177/1358863X20977662
62. Debing E, Peeters E, Demanet C, De Waele M, Van den Brande P. Markers of inflammation in patients with symptomatic and asymptomatic carotid artery stenosis: a case-control study. Vasc Endovascular Surg. (2008) 42:122–7. doi: 10.1177/1538574407307406
63. Jaiswal S, Natarajan P, Silver AJ, Gibson CJ, Bick AG, Shvartz E, et al. Clonal hematopoiesis and risk of atherosclerotic cardiovascular disease. N Engl J Med. (2017) 377:111–21. doi: 10.1056/NEJMoa1701719
64. Ma WF, Hodonsky CJ, Turner AW, Wong D, Song Y, Mosquera JV, et al. Enhanced single-cell RNA-seq workflow reveals coronary artery disease cellular cross-talk and candidate drug targets. Atherosclerosis. (2022) 340:12–22. doi: 10.1016/j.atherosclerosis.2021.11.025
65. Pelisek J, Hegenloh R, Bauer S, Metschl S, Pauli J, Glukha N, et al. Biobanking: objectives, requirements, future challenges-experiences from the munich vascular biobank. J Clin Med. (2019) 8:251. doi: 10.3390/jcm8020251
66. Martinet W, De Meyer GR, Herman AG, Kockx MM. RNA damage in human atherosclerosis: pathophysiological significance and implications for gene expression studies. RNA Biol. (2005) 2:4–7. doi: 10.4161/rna.2.1.1430
67. Kvastad L, Carlberg K, Larsson L, Villacampa EG, Stuckey A, Stenbeck L, et al. The spatial RNA integrity number assay for in situ evaluation of transcriptome quality. Commun Biol. (2021) 4:57. doi: 10.1038/s42003-020-01573-1
68. Lee MY, Bedia JS, Bhate SS, Barlow GL, Phillips D, Fantl WJ, et al. CellSeg: a robust, pre-trained nucleus segmentation and pixel quantification software for highly multiplexed fluorescence images. BMC Bioinformatics. (2022) 23:46. doi: 10.1186/s12859-022-04570-9
69. Baars MJD. Sinha N, Amini M, Pieterman-Bos A, van Dam S, Ganpat MMP, et al. MATISSE: a method for improved single cell segmentation in imaging mass cytometry. BMC Biol. (2021) 19:99. doi: 10.1186/s12915-021-01043-y
70. Keren L, Bosse M, Marquez D, Angoshtari R, Jain S, Varma S, et al. A structured tumor-immune microenvironment in triple negative breast cancer revealed by multiplexed ion beam imaging. Cell. (2018) 174:1373–87 e1319. doi: 10.1016/j.cell.2018.08.039
71. Cable DM, Murray E, Zou LS, Goeva A, Macosko EZ, Chen F, et al. Robust decomposition of cell type mixtures in spatial transcriptomics. Nat Biotechnol. (2022) 40:517–26. doi: 10.1038/s41587-021-00830-w
72. Moncada R, Barkley D, Wagner F, Chiodin M, Devlin JC, Baron M, et al. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol. (2020) 38:333–42. doi: 10.1038/s41587-019-0392-8
73. Sadeghi Rad H, Bazaz SR, Monkman J, Ebrahimi Warkiani M, Rezaei N, O'Byrne K, et al. The evolving landscape of predictive biomarkers in immuno-oncology with a focus on spatial technologies. Clin Transl Immunol. (2020) 9:e1215. doi: 10.1002/cti2.1215
74. Wu T, Liang Z, Zhang Z, Liu C, Zhang L, Gu Y, et al. PRDM16 is a compact myocardium-enriched transcription factor required to maintain compact myocardial cardiomyocyte identity in left ventricle. Circulation. (2022) 145:586–602. doi: 10.1161/CIRCULATIONAHA.121.056666
75. Porritt RA, Zemmour D, Abe M, Lee Y, Narayanan M, Carvalho TT, et al. NLRP3 inflammasome mediates immune-stromal interactions in vasculitis. Circ Res. (2021) 129:e183–200. doi: 10.1161/CIRCRESAHA.121.319153
76. Asp M, Salmén F, Ståhl PL, Vickovic S, Felldin U, Löfling M, et al. Spatial detection of fetal marker genes expressed at low level in adult human heart tissue. Sci Rep. (2017) 7:12941. doi: 10.1038/s41598-017-13462-5
77. Guilliams M, Bonnardel J, Haest B, Vanderborght B, Wagner C, Remmerie A, et al. Spatial proteogenomics reveals distinct and evolutionarily conserved hepatic macrophage niches. Cell. (2022) 185:379–96.e38. doi: 10.1016/j.cell.2021.12.018
78. Tavares-Ferreira D, Shiers S, Ray PR, Wangzhou A, Jeevakumar V, Sankaranarayanan I, et al. Spatial transcriptomics of dorsal root ganglia identifies molecular signatures of human nociceptors. Sci Transl Med. (2022) 14:eabj8186. doi: 10.1126/scitranslmed.abj8186
79. Zhu Y, Wu Z, Yan W, Shao F, Ke B, Jiang X, et al. Allosteric inhibition of SHP2 uncovers aberrant TLR7 trafficking in aggravating psoriasis. EMBO Mol Med. (2022) 14:e14455. doi: 10.15252/emmm.202114455
80. Andersson A, Larsson L, Stenbeck L, Salmén F, Ehinger A, Wu SZ, et al. Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions. Nat Commun. (2021) 12:6012. doi: 10.1038/s41467-021-26271-2
81. Gracia Villacampa E, Larsson L, Mirzazadeh R, Kvastad L, Andersson A, Mollbrink A, et al. Genome-wide spatial expression profiling in formalin-fixed tissues. Cell Genomics. (2021) 1:100065. doi: 10.1016/j.xgen.2021.100065
82. Salmén F, Ståhl PL, Mollbrink A, Navarro JF, Vickovic S, Frisén J, et al. Barcoded solid-phase RNA capture for spatial transcriptomics profiling in mammalian tissue sections. Nat Protoc. (2018) 13:2501–34. doi: 10.1038/s41596-018-0045-2
83. Merritt CR, Ong GT, Church SE, Barker K, Danaher P, Geiss G, et al. Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Nat Biotechnol. (2020) 38:586–99. doi: 10.1038/s41587-020-0472-9
84. Bergholtz H, Carter JM, Cesano A, Cheang MCU, Church SE, Divakar P, et al. Best practices for spatial profiling for breast cancer research with the GeoMx® digital spatial profiler. Cancers. (2021) 13:4456. doi: 10.3390/cancers13174456
85. Smith KD, Prince DK, Henriksen KJ, Nicosia RF, Alpers CE, Akilesh S. Digital spatial profiling of collapsing glomerulopathy. Kidney Int. (2022) 101:1017–26. doi: 10.1016/j.kint.2022.01.033
86. Dottermusch M, Schumann Y, Kordes U, Hasselblatt M, Neumann JE. Spatial molecular profiling of a central nervous system low-grade diffusely infiltrative tumour with INI1 deficiency featuring a high-grade atypical teratoid/rhabdoid tumour component. Neuropathol Appl Neurobiol. (2022) 48:e12777. doi: 10.1111/nan.12777
87. Delorey TM, Ziegler CGK, Heimberg G, Normand R, Yang Y, Segerstolpe Å, et al. COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets. Nature. (2021) 595:107–13. doi: 10.1038/s41586-021-03570-8
88. Schlam I, Church SE, Hether TD, Chaldekas K, Hudson BM, White AM, et al. The tumor immune microenvironment of primary and metastatic HER2- positive breast cancers utilizing gene expression and spatial proteomic profiling. J Transl Med. (2021) 19:480. doi: 10.1186/s12967-021-03113-9
89. Bodenmiller B. Multiplexed epitope-based tissue imaging for discovery and healthcare applications. Cell Syst. (2016) 2:225–38. doi: 10.1016/j.cels.2016.03.008
90. Baharlou H, Canete NP, Cunningham AL, Harman AN, Patrick E. Mass cytometry imaging for the study of human diseases-applications and data analysis strategies. Front Immunol. (2019) 10:2657. doi: 10.3389/fimmu.2019.02657
91. Schürch CM, Bhate SS, Barlow GL, Phillips DJ, Noti L, Zlobec I, et al. Coordinated cellular neighborhoods orchestrate antitumoral immunity at the colorectal cancer invasive. Front Cell. (2020) 182:1341–59.e19. doi: 10.1016/j.cell.2020.07.005
92. Gouin KH, Ing N, Plummer JT, Rosser CJ. ben Cheikh B, Oh C, et al. An N-Cadherin 2 expressing epithelial cell subpopulation predicts response to surgery, chemotherapy and immunotherapy in bladder cancer. Nat Commun. (2021) 12:4906. doi: 10.1038/s41467-021-25103-7
93. Palla AR, Ravichandran M, Wang YX, Alexandrova L, Yang A v, Kraft P, et al. Inhibition of prostaglandin-degrading enzyme 15-PGDH rejuvenates aged muscle mass and strength. Science. (2021) 371. doi: 10.1126/science.abc8059
94. Goltsev Y, Samusik N, Kennedy-Darling J, Bhate S, Hale M, Vazquez G, et al. Deep profiling of mouse splenic architecture with CODEX multiplexed imaging. Cell. (2018) 174:968–81.e15. doi: 10.1016/j.cell.2018.07.010
95. Fluidigm Web Page (2022). Available online at: https://www.fluidigm.com (accessed April 28, 2022).
96. Abdulrahman Z, Santegoets SJ, Sturm G, Charoentong P, Ijsselsteijn ME, Somarakis A, et al. Tumor-specific T cells support chemokine-driven spatial organization of intratumoral immune microaggregates needed for long survival. J Immunother Cancer. (2022) 10:e004346. doi: 10.1136/jitc-2021-004346
97. Zheng Y, Han L, Chen Z, Li Y, Zhou B, Hu R, et al. PD-L1+CD8+ T cells enrichment in lung cancer exerted regulatory function and tumor-promoting tolerance. IScience. (2022) 25:103785. doi: 10.1016/j.isci.2022.103785
98. Feng C, Wang X, Tao Y, Xie Y, Lai Z, Li Z, et al. Single-cell proteomic analysis dissects the complexity of tumor microenvironment in muscle invasive bladder cancer. Cancers. (2021) 13:5440. doi: 10.3390/cancers13215440
99. Chen C, McDonald D, Blain A, Sachdeva A, Bone L, Smith ALM, et al. Imaging mass cytometry reveals generalised deficiency in OXPHOS complexes in Parkinson's disease. NPJ Parkinsons Dis. (2021) 7:39. doi: 10.1038/s41531-021-00182-x
100. Ionpath Web Page (2022). Available online at: https://www.ionpath.com (accessed April 28, 2022).
101. Rovira-Clavé X, Jiang S, Bai Y, Zhu B, Barlow G, Bhate S, et al. Subcellular localization of biomolecules and drug distribution by high-definition ion beam imaging. Nat Commun. (2021) 12:4628. doi: 10.1038/s41467-021-24822-1
102. Risom T, Glass DR, Averbukh I, Liu CC, Baranski A, Kagel A, et al. Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma. Cell. (2022) 185:299–310.e18. doi: 10.1016/j.cell.2021.12.023
103. McCaffrey EF, Donato M, Keren L, Chen Z, Delmastro A, Fitzpatrick MB, et al. The immunoregulatory landscape of human tuberculosis granulomas. Nat Immunol. (2022) 23:318–29. doi: 10.1038/s41590-021-01121-x
104. Liu CC, Bosse M, Kong A, Kagel A, Kinders R, Hewitt SM, et al. Reproducible, high-dimensional imaging in archival human tissue by multiplexed ion beam imaging by time-of-flight (MIBI-TOF). Lab Invest. (2022). doi: 10.1038/s41374-022-00778-8
Keywords: atherosclerosis, biobanked human biospecimens, next-generation sequencing, bioinformatic analyses, spatial biology
Citation: Elishaev M, Hodonsky CJ, Ghosh SKB, Finn AV, von Scheidt M and Wang Y (2022) Opportunities and Challenges in Understanding Atherosclerosis by Human Biospecimen Studies. Front. Cardiovasc. Med. 9:948492. doi: 10.3389/fcvm.2022.948492
Received: 19 May 2022; Accepted: 07 June 2022;
Published: 07 July 2022.
Edited by:
Vicente Andres, Spanish National Centre for Cardiovascular Research, SpainReviewed by:
Jose Luis Martin-Ventura, Universidade Autónoma de Lisboa, PortugalCopyright © 2022 Elishaev, Hodonsky, Ghosh, Finn, von Scheidt and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ying Wang, Wang85@mail.ubc.ca
orcid.org/0000-0002-1444-5778