- 1Department of Cardiology, Taizhou Hospital Affiliated to Wenzhou Medical University, Taizhou, China
- 2Department of Cardiology, Zhejiang Hospital, Hangzhou, China
- 3Department of Medicine, Crozer-Chester Medical Center, Upland, PA, United States
- 4Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 5Department of Surgery, Stanford University School of Medicine, Stanford, CA, United States
- 6Department of Radiation Oncology, Indiana University School of Medicine, Indianapolis, IN, United States
Purpose: This study aims to illustrate the cellular landscape in the aorta of experimental aortic dissection (AD) and elaborate on the smooth muscle cells (SMCs) heterogeneity and functions among various cell types.
Methods: Male Apolipoprotein deficient (ApoE−/−) mice at 28 weeks of age were infused with Ang II (2,500 ng/kg/min) to induce AD. Aortas from euthanized mice were harvested after 7 days for 10×Genomics single-cell RNA sequencing (scRNA-seq), followed by the identification of cell types and differentially expressed genes (DEGs). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was conducted.
Results: AD was successfully induced in ApoE−/− mice. scRNA-seq identified 15 cell clusters and nine cell types, including non-immune cells (endothelials, fibroblasts, and SMCs) and immune cells (B cells, natural killer T cell, macrophages, dendritic cells, neutrophils, and mast cells). The relative numbers of SMCs were remarkably changed, and seven core DEGs (ACTA2,IL6,CTGF,BGN,ITGA8,THBS1, and CDH5) were identified in SMCs. Moreover, we found SMCs can differentiate into 8 different subtypes through single-cell trajectory analysis.
Conclusion: scRNA-seq technology can successfully identify unique cell composition in experimental AD. To our knowledge, this is the first study that provided the complete cellular landscape in AD tissues from mice, seven core DEGs and eight subtypes of SMCs were identified, and the SMCs have evolution from matrix type to inflammatory type.
Introduction
Aortic dissection (AD) is a life-threatening aortic dilatation disease with a significant fatality rate as high as 1–2% per hour after symptom onset in untreated patients (Gawinecka et al., 2017). Epidemiological data revealed that 93.5% of AD patients in China are under the age of 70, with an average age of onset of 52 years-old (Wang et al., 2014). The pathophysiology behind remains unclear and the main clinical treatment modalities are surgery and endovascular repair as effective drug treatment are still unavailable to stop the progression (Wang et al., 2018).It necessary to have a deeper understanding of the pathogenesis of AD to improve clinical outcomes. The occurrence and development of AD are complex, comprising of gene mutation, hypertension, blood lipid level, angiotensin II (Ang II) and inflammation (Fan et al., 2014; Qi et al., 2018; Renard et al., 2018; Tanaka et al., 2018; Wolford et al., 2019). Common pathological changes of AD patients include extensive destruction of vascular elastic fiber dimension, loss of the adventitia elastic fiber and neutrophil infiltration (Yuan and Wu, 2018).However, the cell types contributing to AD progression have not been fully elucidated. Single-cell RNA sequencing (scRNA-seq) technology is to sequence the RNA of a single cell, which can generate the unique transcription of an individual cell at single-cell resolution and establish a gene regulation network (Lambrechts et al., 2018; Han et al., 2020). The cellular heterogeneity of aneurysmal abdominal aorta has been described by scRNA-seq (Zhao et al., 2021). A mouse endothelial cell map and the heterogeneity of mouse endothelial cells has also been reported recently (Kalucka et al., 2020). In our current study, we successfully constructed a mouse model of AD by using Ang II followed by scRNA-seq technology to clarify the diversity of the cell profiles of vascular tissue and identify marker genes for different cell clusters. We further analyzed the smooth muscle cells (SMCs) and clarified the core genes of AD formation and development trajectory, which provided a new perspective on the complex biological systems and key functions associated with disease development.
Methods
Mouse Model
28-weeks-old male ApoE−/− mice were purchased from Zhejiang Academy of Medical Sciences and raised in the animal room of Zhejiang Academy of Medical Sciences under specific pathogen-free (SPF) environment. Mice were housed in a temperature-controlled room at 25 ± 2°C and humidity of 55 ± 3%, under 12 h light and dark cycle. ApoE-/- mice were randomly divided into two groups, with 15 mice in each control and study group, infused for 7 days with Ang II (2500 ng/kg/min) and saline respectively.
Ultrasonography Examination
After 7 days of aorta modeling, ultrasound was used to evaluate the formation of abdominal AD. A MS550 S probe on the small animal ultrasonic imaging system was adjusted at center frequency of 40MHz, focal length of 6 mm, fixed focus depth of 10 mm, and width of 12.00 mm. After the mice were anesthetized with 1.5–2% isoflurane gas, the hair from the chest to the abdomen of the mice was removed with depilation cream followed by application of a warm paper towel to wipe dry the depilated skin. After fixating the mice on the examination table, ultrasonic gel was applied to the area of interest prior to the ultrasonography examination.
Specimen Collection and Pathological Examination
After the ultrasonography examination was performed, the aorta tissues were removed and preserved in 4% paraformaldehyde solution for 24 h. The tissues were then embedded in paraffin solution and subsequently sectioned evenly to the pieces of 5 μm thickness. The specimens were then stained with hematoxylin-eosin (H&E) and van Gieson (EVG)stains and examined using a light microscope.
Aortic Vascular Tissue Procurement and Isolation
The harvested aortic tissue was washed three times with 0.9% normal saline on ice. The samples were washed with a medium and appropriate amount of enzyme (1X, Gibco, 12605010) were added followed by incubation for 1 h. After enzyme digestion, the samples were placed on a stainless steel cell filter and washed with culture medium for two times followed by centrifugation at 4°C and 300 g for 5 min. Red blood cells were subsequently removed using 1X RBC lysis buffer (Cat: R1010, Solarbio, China) for 5 min and cells were centrifuged at 4°C and 300 g for 5 min. The final samples were washed with medium for two times.
Single-Cell RNA Sequencing
Sc-RNA seq was performed using the 10x Genomics’ single-cell solution. Single-cell suspensions from 4 aortas in control group were pooled together as one sample and three aortas in experimental group were pooled together as another sample respectively. The emulsion and the library were prepared according to the user guide of 10x Genomics Chromium Single Cell 3’v3 Reagent Kit. GEMs (Gel Bead-in-Emulsions) are generated by combining specimens, a Master Mix, gel beads and partitioning oil, onto the chip in the chromium single cell controller kit. After the cell is lysed and RNA reverse transcription occurs, the cDNA was purified from the mixture. The cDNA was amplified, fragmented, and spliced to generate sufficient mass for library construction. The cDNA was subsequently sequenced by Illumina sequencer.
Single-Cell RNA Sequencing Data Analysis
A criterion to filter out cells with UMI/gene numbers out of the limit of mean value +/- 2 fold of standard deviations was applied to remove cells of low quality and likely multiple captures, which is a major concern in microdroplet-based experiments. Low-quality cells where >10% of the counts belonged to mitochondrial genes were further discarded following visual inspection of the distribution of cells by the fraction of mitochondrial genes expressed.
Identification of Cluster Marker Genes and Cell Types
Principal component analysis (PCA) was performed to reduce the dimensionality with RunPCA function in Seurat (Butler et al., 2018).A 2-dimensional t-distributed stochastic neighbor embedding (t-SNE) algorithm with the RunTSNE function in Seurat was used to identify cell clusters (Butler et al., 2018). The FindAllMarkers function (test.use = bimod) in Seurat helps to identify marker genes of each cluster. We used the R package SingleR (Aran et al., 2019), a novel computational method for unbiased cell type recognition of scRNA-seq with Mapping the Mouse Cell Atlas by Microwell-Seq as a reference to infer the cell origin of each of the single cells independently and identify cell types (Han et al., 2018).
Identification of DEGs
The FindMarkers function (test.use = MAST) in Seurat (Butler et al., 2018) was utilized to identify differentially expressed genes (DEGs) between aortic dissection group and normal control group samples in RStudio.The criteria were as follows: (1) the p < 0:01; (2) log2FoldChange≥1.5 or ≤ -1.5.
Protein-Protein Interaction Analysis
In order to determine the core genes in the DEGs of SMCs, we used online tool STRING (version 11.5) (http://string-db.org/) to construct the PPI network, and the parameters were set to the default value.
Hub Genes Selection and Analyses in SMCs
The cytohubba (Chin et al., 2014) plug-in in Cytoscape was used to screen the core genes in PPI network, We randomly selected 5 of the 12 algorithms in the cytohubba plugin, and then used the intersection of the results of the 4 algorithms to determine the hub gene. After screening the core gene, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to determine the function of the core gene.
GO Enrichment Analysis and KEGG Pathway Analysis
GO and KEGG pathway enrichment analyses for DEGs in SMCs was performed using the R package clusterProfile (Yu et al., 2012), which provided a comprehensive set of functional annotation tools allowing investigators to understand the biological meaning behind a large list of genes.
Reconstructing SMCs Differentiation Trajectories by Monocle2
SMCs pseudotime trajectories were reconstructed by the Monocle2 R package (version 4.1.2). Briefly, the work flow chart is as follows: (1) the types of SMCs were selected by Seurat; (2) Reduces the high-dimensional expression profiles to a low-dimensional space; (3) Single cells are projected onto this space and ordered into a trajectory with branch points.
Results
Ultrasonographic Results
After Ang II infusion, there was dilatation of the aorta, progressing to separation of the blood vessel walls, and subsequently forming a false lumen. This was consistent with the ultrasound findings of aortic dissection (Figures 1A,C), while the vascular wall of the abdominal aorta was normal in the control group (Figures 1B,D).
FIGURE 1. The aortic dissection induced by angiotensin II infusion. (A) Gross morphology of aortic dissection. (B) Gross morphology of normal blood vessels. (C)Ultrasound image of aortic dissection in the experimental group. (D) Ultrasound image of normal vessels in the control group.
Pathological Examination
Aortic sections were stained with H&E and EVG. In the control group, vascular wall structure had no changes with a well-organized smooth muscle layer and no change in wall thickness or formation of pseudo-lumen (Figures 2A,C). In the study group, the vascular intima was torn with a decreased wall thickness. Medial elastic fibers were disordered in conjunction with adventitial inflammatory cell infiltration, intramural thrombus and pseudo-lumen (Figures 2B,D).
FIGURE 2. EVG and H&E staining of aortic vessels in mice from each experimental group. (A) H&E staining of normal blood vessels. (B) H&E staining of aortic dissection blood vessels. (C) EVG staining of normal blood vessels. (D) EVG staining of aortic dissection blood vessels.
Single-Cell RNA Sequencing Clustering Analysis
SCRNA-seq analysis was performed on vascular tissue from 2 groups as illustrated in Figure 3A. 18,835 cells were used for integrated single-cell RNA-Seq analysis. Based on the similarity of gene expression, 15 clusters of cells and nine cell types were identified (Figures 3B,C). The key marker gene of each cell type was shown in Table 1. The major cell types comprised of: (i) Fibroblasts (Clusters1,2); (ii) Endothelial cells (Clusters 9,13,14); (iii) Smooth muscle cells (Cluster 6,12); (iv) B cells (Clusters 3); (v) Macrophage & Monocyte (Cluster 4,8); (vi) DC (Cluster 11); (vii) NK-T cells (Cluster 5,15) (viii) Mast cells (Cluster 8); (ix) Neutrophils (Cluster 7), the key marker genes for each cell cluster are presented in Figure 3D.
FIGURE 3. Identification of cell clusters and types by single-cell RNA sequencing (scRNA-seq). (A) Schematic of the experimental design for single-cell RNA sequencing. (B) Distribution of cells in AD and control groups.(C) t-SNE plot demonstrated nine cell types. (D) the key marker genes in 15 clusters.(E) Pie chart of cell types and clusters.
Differences in Cell Type Among Control and AD Groups
The 6,323 cells from the control group were assigned to nine different cell types (Figure 3E) with the following proportions: 1.20% of B cells, 1.5% of DC, 7.3% of endothelial cells, 51.2% of fibroblast, 17.1% of macrophages and monocyte, 0.16% of mast cells, 0.35% of neutrophils, 2.78% of NK-T cell, and 18.31% of SMCs. In AD group, 10756 cells were clustered depending on the similarity of gene expression, and 9 cell types were identified. The main cell types with their corresponding proportions were as follows: 22.1% of B cells, 2.5% of DC, 4.7% of endothelial cells, 34.2% of fibroblasts, 13.8% of macrophages and monocytes, 0.30% of mast cells, 6.56% of neutrophils, 13.85% of NK-T cells, and 2.06% of SMCs. Besides, we noticed significant changes in SMCs, which decreased from 18.3% of normal tissue to 2.1% of dissection tissue, In our previous study, the expression of Extracellular matrix metalloproteinase inducer (EMMPRIN) in SMCs will significantly affect the progress of AD (Chen et al., 2009), this suggests that SMCs may play an important role in AD and deserves priority attention.
Identification of DEGs and PPI Network Construction in SMCs
In SMCs, we identified 59 up-regulated genes and 51 down-regulated genes (Figure 4A), After importing DEGs into STRING database, an interaction network with 86 nodes and 440 sides was obtained (Figure 4B).
FIGURE 4. Functional analysis of smooth muscle cells. (A) Volcanic map of DEGs in smooth muscle cells. (B) PPI network of DEGs. (C) seven hub genes were identified by overlapping the first 15 genes in the four classification methods of cytoHubba. (D) GO enrichment analysis of seven hub genes. (E) KEGG enrichment analysis of seven hub genes.
Hub Gene Selection and Analysis
In this study, we used cytohubba to identify core genes, according to the four algorithms in cytohubba (Bottleneck, Degree, EPC, Closeness), the ranking of the top 15 core genes selected by each algorithm is shown in Table 2. Finally, seven central genes (ACTA2, IL6, CTGF, BGN,ITGA8,THBS1, and CDH5) were identified by overlapping the15 genes (Figure 4C). GO and KEGG enrichment analyses were conducted based on the seven central genes of SMCs. The significant GO terms of core genes were “extracellular matrix organization”, “extracellular structure organization”, and “regulation of epithelial cell apoptotic process” (Figure 4D), and indicates that these seven core genes play an significant role in changing the function of extracellular matrix. In KEGG enrichment analysis, only three core genes (ITGA8, THBS1, and IL6) were enriched in four signal pathways: “malaria”, “PI3K−Akt signaling pathway”, “ECM−receiver interaction”, and “hypertrophic cardiopathway” (Figure 4E), which imply that these genes play a role in promoting metabolism, growth and angiogenesis in SMCs.
The Trajectory Analysis of SMCs in AD
The changes of cell phenotype and function occur gradually. Due to the great changes in the number of SMCs, we applied Monocle2 to perform pseudotime trajectories of all SMCs and showed the differentiation direction between them. Through the expression of characteristic genes, all SMCs could be divided into 8 types (Figures 5A,B), and the marker genes of each subtype are shown in Supplementary Material. SMCs mainly differentiate from SMC2 to SMC5 and undergo a variety of phenotypic changes (Figures 5E,F). Besides, we divided 8 kinds of SMCs into three modules according to the changes of gene expression. Module1:The expression of genes first increased and then decreased with the development of AD; Module2: The expression of genes decreased with the development of AD; Module3: The expression of genes increased with the development of AD (Figure 5C). In addition, it can be seen from Figure 5F that SMCs differentiate from SMC2 and are finally divided into SMC5, SMC6 and SMC7. In order to understand the role of each SMC subtype, we further performed GO function analysis.GO enrichment analysis showed that SMC2 was mainly related to “regulation of actin cytoskeleton organization”, “cell junction”, and “cell adhesion molecular binding”, (Figure 6A), suggesting that it was a stromal SMC; SMC6 is mainly associated with “regulation of calcium ion transmembrane transport”, “calcium signaling pathway”, and “positive regulation of calcium ion-dependent exocytosis” (Figure 6B), suggesting that it is a signal transduction type smooth muscle cell; SMC7 is mainly related to “endoderm formation”, “cytoplastmm” and “analytical structure morphogenesis” (Figure 6C). We believe that SMC7 is a transitional SMC after SMC6 signal transduction. SMC5 is related to “cellular response to tumor necrosis factor”, “immune response”, “cellular response to interleukin −1″, suggesting that it is an inflammatory SMC (Figure 6D).
FIGURE 5. Analysis of smooth muscle cell subgroup .(A)Heatmap of the top 10 enriched genes in each SMCs population. (B) Distribution of SMCs Subgroup between the two groups. (C)The heatmap of top 50 genes which had most critical influences on cell transformation. (D) Time progression of SMCs differentiation. (E) and (F) SMCs differentiation trajectory.
FIGURE 6. GOanalysis of smooth muscle cell subtypes. (A) GO enrichment analysis of SMC2. (B) GO enrichment analysis of SMC6. (C) GO enrichment analysis of SMC7. (D) GO enrichment analysis of SMC5.
Discussion
AD is a medical emergency. However, the underlying mechanism of the disease has remained unclear. Currently, there is no identifiable effective pharmacologic management to treat AD. In this study, we successfully constructed an AD model using Ang II with an input of 2,500 ng/kgmin within 7 days. The histopathological results of the model constructed demonstrated similar pathological manifestations of AD in humans, including tearing of tunica intima, rupturing of elastic intima, and infiltration of inflammatory cells on the adventitia (Del Porto et al., 2010; Wang et al., 2020).
The present study aims to explore the specificity of cell types in the aorta of experimental AD using scRNA-seq. We successfully identified 15 clusters of cells which were further divided into nine types of cells based on their marker genes, including Fibroblast, Endothelial, SMC, B cells, Macrophage, DC,T-NK cells, Mast cells, and Neutrophil. We noticed that there were significant changes in SMCs in the formation of AD. Besides, our previous study found that the expression of EMMPRIN in SMCs can promote the formation and progression of AD (Chen et al., 2009). Therefore, we focus on SMCs for further analysis.
Seven core genes were identified in SMCs. ACTA2, is a smooth muscle actin, ranks first among the four algorithms. The mutation of ACTA2 is closely related to aortic dissection (Cooper and Brown, 2017; van de Laar et al., 2019).The second key gene IL6, encodes a cytokine that functions in inflammation and the maturation of B cells, the inflammatory mechanism plays a key role in the development and progression of AD (Lu et al., 2021).CTGF, is a Protein Coding gene, and has important pro-inflammatory and pro-fibrotic characteristic (Morales et al., 2018), which play an important role in the progression of AD (Bode-Jänisch et al., 2012; Bersi et al., 2017; Zhao and Bie, 2020). BGN, encodes a member of the small leucine-rich proteoglycan (SLRP) family of proteins, which plays a role in bone growth, muscle development and regeneration, and collagen fibril assembly in multiple tissues (Fisher et al., 1991; Didangelos et al., 2010). ITGA8 mediate numerous cellular processes including cell adhesion, cytoskeletal rearrangement, and activation of cell signaling pathways (Gaudet et al., 2011), The study of organ injury and fibrosis in rat and mouse models showed that ITGA8 was expressed in vascular SMCs and stromal cells after injury (Nishimura et al., 1991), this suggests that it may play a similar role in the formation of AD. THBS1 is an adhesive glycoprotein that mediates cell-to-cell and cell-to-matrix interactions (Mumby et al., 1984; Staniszewska et al., 2007), In addition, THBS1 can activate TGFβ and matrix metalloproteinase (MMP), which mediates the interaction between cells and matrix and participates in angiogenesis, proliferation and platelet aggregation (Savarimuthu Francis et al., 2011).The finally gene CDH5, encodes a classical cadherin of the cadherin superfamily, it is essential for vascular integrity and endothelial function (Mao et al., 2013). The GO and KEGG enrichment analysis also displayed that the above genes play an crucial role in regulating extracellular matrix and mediating inflammatory signaling pathways.
In single-cell trajectory analysis, we identified 8 subtypes of SMCs. SMC5, SMC6, and SMC7 are the main component of vascular smooth muscle in AD, After GO enrichment analysis of vascular SMCs, we consider that SMC2 is a stromal vascular smooth muscle cell, which is mainly involved in the change of extracellular matrix function. Extracellular matrix remodeling occurs with the increase of aortic stiffness and the gradual weakening of aortic wall, resulting in AD (Bhushan et al., 2019),Therefore, SMC2 is an early cell in the formation of AD.SMC6 is related to the signal transduction of calcium ions, so we consider it as a signal transduction type SMC. In a clinical study, calcium levels in elderly patients with acute AD may help to change normal endothelial physiological function (Vianello et al., 2017), In addition, studies have confirmed human vascular smooth muscle cells (HASMCs) can regulate the functional activity of SMCs through intracellular calciu (Sun et al., 2017). Therefore, SMC6 is considered to play a role in phenotypic transformation in the middle stage of AD formation.SMC7 is related to endoderm formation, cytolastmm and structure morphogenesis, so it may be a transitional SMC before final morphology. SMC5 is associated to the secretion of inflammatory factors, and macrophage-induced inflammation are critical for aortic dissection (Li et al., 2017), exacerbated inflammation contributes to the formation of AD and eventually leads to the tearing of blood vessels (Lian et al., 2019; Wu et al., 2021),Therefore, SMC5 may play a major role in AD rupture.
In summary, our study successfully constructed a mouse model of AD using Ang II, and to our knowledge, this is the first study that using scRNA-seq to reveal the cellular components in mouse AD tissues at a single-cell level. Our research confirmed that SMCs contribute a significant role in the formation of AD, and successfully identified seven key genes that associated with extracellular matrix function in SMCs, namely ACTA2, IL6, CTGF, BGN, ITGA8, THBS1, and CDH5, providing us new therapeutic targets for treatment of AD. We also found the changes of four main cell states of SMCs from matrix type to inflammatory type during the formation of AD, which improved our understanding of the pathogenesis and underlying molecular mechanism in the formation of AD.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Gene Expression Omnibus, accession number GSE203594.
Ethics Statement
The animal experiments were approved by the Laboratory Animal Welfare Ethics Committee of Zhejiang Academy of Medical Sciences (No. 20200271). Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
XML interpretation of data for the work, drafting the work, revising it for important intellectual content; and final approval of the version to be published. JR and BX visualization, supervision the work, agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy and final approval of the version to be published.
Funding
This study was supported by the National Natural Science Foundation of China (Grant No. 81770475).
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 thank Dr. Shaomeng Yu and Wenya Lu (OE Biotech Co., Ltd., Shanghai, China) for assistance with bioinformatics analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.836593/full#supplementary-material
Supplementary Figure S1 | Smooth muscle cell subtype marker.
References
Aran, D., Looney, A. P., Liu, L., Wu, E., Fong, V., Hsu, A., et al. (2019). Reference-based Analysis of Lung Single-Cell Sequencing Reveals a Transitional Profibrotic Macrophage. Nat. Immunol. 20 (2), 163–172. doi:10.1038/s41590-018-0276-y
Bersi, M. R., Khosravi, R., Wujciak, A. J., Harrison, D. G., and Humphrey, J. D. (2017). Differential Cell-Matrix Mechanoadaptations and Inflammation Drive Regional Propensities to Aortic Fibrosis, Aneurysm or Dissection in Hypertension. J. R. Soc. Interface 14 (136), 20170327. doi:10.1098/rsif.2017.0327
Bhushan, R., Altinbas, L., Jäger, M., Zaradzki, M., Lehmann, D., Timmermann, B., et al. (2019). An Integrative Systems Approach Identifies Novel Candidates in Marfan Syndrome‐related Pathophysiology. J. Cell. Mol. Medi 23 (4), 2526–2535. doi:10.1111/jcmm.14137
Bode-Jänisch, S., Schmidt, A., Günther, D., Stuhrmann, M., and Fieguth, A. (2012). Aortic Dissecting Aneurysms-Hhistopathological Findings. Forensic Sci. Int. 214 (1-3), 13–17. doi:10.1016/j.forsciint.2011.07.006
Butler, A., Hoffman, P., Smibert, P., Papalexi, E., and Satija, R. (2018). Integrating Single-Cell Transcriptomic Data across Different Conditions, Technologies, and Species. Nat. Biotechnol. 36 (5), 411–420. doi:10.1038/nbt.4096
Chen, X.-f., Wang, J.-a., Hou, J., Gui, C., Tang, L.-j., Chen, X.-q., et al. (2009). Extracellular Matrix Metalloproteinase Inducer (EMMPRIN) Is Present in Smooth Muscle Cells of Human Aneurysmal Aorta and Is Induced by Angiotensin II In Vitro. Clin. Sci. (Lond) 116 (11), 819–826. doi:10.1042/cs20080235
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol. 8 Suppl 4 (Suppl. 4), S11. doi:10.1186/1752-0509-8-S4-S11
Cooper, K., and Brown, S. (2017). ACTA2 Mutation and Postpartum Hemorrhage: a Case Report. BMC Med. Genet. 18 (1), 143. doi:10.1186/s12881-017-0505-5
Del Porto, F., Proietta, M., Tritapepe, L., Miraldi, F., Koverech, A., Cardelli, P., et al. (2010). Inflammation and Immune Response in Acute Aortic Dissection. Ann. Med. 42 (8), 622–629. doi:10.3109/07853890.2010.518156
Didangelos, A., Yin, X., Mandal, K., Baumert, M., Jahangiri, M., and Mayr, M. (2010). Proteomics Characterization of Extracellular Space Components in the Human Aorta. Mol. Cell. Proteomics 9 (9), 2048–2062. doi:10.1074/mcp.m110.001693
Fan, L. M., Douglas, G., Bendall, J. K., McNeill, E., Crabtree, M. J., Hale, A. B., et al. (2014). Endothelial Cell-specific Reactive Oxygen Species Production Increases Susceptibility to Aortic Dissection. Circulation 129 (25), 2661–2672. doi:10.1161/circulationaha.113.005062
Fisher, L. W., Heegaard, A. M., Vetter, U., Vogel, W., Just, W., Termine, J. D., et al. (1991). Human Biglycan Gene. Putative Promoter, Intron-Exon Junctions, and Chromosomal Localization. J. Biol. Chem. 266 (22), 14371–14377. doi:10.1016/s0021-9258(18)98694-1
Gaudet, P., Livstone, M. S., Lewis, S. E., and Thomas, P. D. (2011). Phylogenetic-based Propagation of Functional Annotations within the Gene Ontology Consortium. Briefings Bioinforma. 12 (5), 449–462. doi:10.1093/bib/bbr042
Gawinecka, J., Schönrath, F., and von Eckardstein, A. (2017). Acute Aortic Dissection: Pathogenesis, Risk Factors and Diagnosis. Swiss Med. Wkly. 147, w14489. doi:10.4414/smw.2017.14489
Han, X., Wang, R., Zhou, Y., Fei, L., Sun, H., Lai, S., et al. (2018). Mapping the Mouse Cell Atlas by Microwell-Seq. Cell 172 (5), 1091–1107. doi:10.1016/j.cell.2018.02.001
Han, X., Zhou, Z., Fei, L., Sun, H., Wang, R., Chen, Y., et al. (2020). Construction of a Human Cell Landscape at Single-Cell Level. Nature 581 (7808), 303–309. doi:10.1038/s41586-020-2157-4
Kalucka, J., de Rooij, L. P. M. H., Goveia, J., Rohlenova, K., Dumas, S. J., Meta, E., et al. (2020). Single-Cell Transcriptome Atlas of Murine Endothelial Cells. Cell 180 (4), 764–779. doi:10.1016/j.cell.2020.01.015
Lambrechts, D., Wauters, E., Boeckx, B., Aibar, S., Nittner, D., Burton, O., et al. (2018). Phenotype Molding of Stromal Cells in the Lung Tumor Microenvironment. Nat. Med. 24 (8), 1277–1289. doi:10.1038/s41591-018-0096-5
Li, X., Fang, Q., Tian, X., Wang, X., Ao, Q., Hou, W., et al. (2017). Curcumin Attenuates the Development of Thoracic Aortic Aneurysm by Inhibiting VEGF Expression and Inflammation. Mol. Med. Rep. 16 (4), 4455–4462. doi:10.3892/mmr.2017.7169
Lian, G., Li, X., Zhang, L., Zhang, Y., Sun, L., Zhang, X., et al. (2019). Macrophage Metabolic Reprogramming Aggravates Aortic Dissection through the HIF1α-ADAM17 Pathway✰. EBioMedicine 49, 291–304. doi:10.1016/j.ebiom.2019.09.041
Lu, H. Y., Hsu, H. L., Li, C. H., Li, S. J., Lin, S. J., Shih, C. M., et al. (2021). Hydrogen Sulfide Attenuates Aortic Remodeling in Aortic Dissection Associating with Moderated Inflammation and Oxidative Stress through a NO-dependent Pathway. Antioxidants (Basel) 10 (5), 682. doi:10.3390/antiox10050682
Mao, X.-g., Xue, X.-y., Wang, L., Zhang, X., Yan, M., Tu, Y.-y., et al. (2013). CDH5 Is Specifically Activated in Glioblastoma Stemlike Cells and Contributes to Vasculogenic Mimicry Induced by Hypoxia. Neuro-Oncology 15 (7), 865–879. doi:10.1093/neuonc/not029
Morales, M. G., Acuña, M. J., Cabrera, D., Goldschmeding, R., and Brandan, E. (2018). The Pro-fibrotic Connective Tissue Growth Factor (CTGF/CCN2) Correlates with the Number of Necrotic-Regenerative Foci in Dystrophic Muscle. J. Cell Commun. Signal. 12 (1), 413–421. doi:10.1007/s12079-017-0409-3
Mumby, S. M., Raugi, G. J., and Bornstein, P. (1984). Interactions of Thrombospondin with Extracellular Matrix Proteins: Selective Binding to Type V Collagen. J. Cell Biol. 98 (2), 646–652. doi:10.1083/jcb.98.2.646
Nishimura, R. N., Dwyer, B. E., Clegg, K., Cole, R., and de Vellis, J. (1991). Comparison of the Heat Shock Response in Cultured Cortical Neurons and Astrocytes. Brain Res. Mol. Brain Res. 9 (1-2), 39–45. doi:10.1016/0169-328x(91)90128-k
Qi, Y. F., Shu, C., Xiao, Z. X., Luo, M. Y., Fang, K., Guo, Y. Y., et al. (2018). Post-transcriptional Control of Tropoelastin in Aortic Smooth Muscle Cells Affects Aortic Dissection Onset. Mol. Cells 41 (3), 198–206. doi:10.14348/molcells.2018.2193
Renard, M., Francis, C., Ghosh, R., Scott, A. F., Witmer, P. D., Adès, L. C., et al. (2018). Clinical Validity of Genes for Heritable Thoracic Aortic Aneurysm and Dissection. J. Am. Coll. Cardiol. 72 (6), 605–615. doi:10.1016/j.jacc.2018.04.089
Savarimuthu Francis, S. M., Larsen, J. E., Pavey, S. J., Duhig, E. E., Clarke, B. E., Bowman, R. V., et al. (2011). Genes and Gene Ontologies Common to Airflow Obstruction and Emphysema in the Lungs of Patients with COPD. PLoS One 6 (3), e17442. doi:10.1371/journal.pone.0017442
Staniszewska, I., Zaveri, S., Valle, L. D., Oliva, I., Rothman, V. L., Croul, S. E., et al. (2007). Interaction of α9β1 Integrin with Thrombospondin-1 Promotes Angiogenesis. Circulation Res. 100 (9), 1308–1316. doi:10.1161/01.res.0000266662.98355.66
Sun, Y., Zhao, Z., Hou, L., Xiao, Y., Qin, F., Yan, J., et al. (2017). The Regulatory Role of Smooth Muscle 22 on the Proliferation of Aortic Smooth Muscle Cells Participates in the Development of Aortic Dissection. J. Vasc. Surg. 66 (3), 875–882. doi:10.1016/j.jvs.2016.02.063
Tanaka, H., Iida, Y., Iwaki, T., Suzuki, Y., Sano, H., Miyajima, C., et al. (2018). Elevated Plasma Levels of LDL Cholesterol Promote Dissecting Thoracic Aortic Aneurysms in Angiotensin II-Induced Mice. Ann. Vasc. Surg. 48, 204–213. doi:10.1016/j.avsg.2017.10.006
van de Laar, I. M. B. H., Arbustini, E., Loeys, B., Björck, E., Murphy, L., Groenink, M., et al. (2019). European Reference Network for Rare Vascular Diseases (VASCERN) Consensus Statement for the Screening and Management of Patients with Pathogenic ACTA2 Variants. Orphanet J. Rare Dis. 14 (1), 264. doi:10.1186/s13023-019-1186-2
Vianello, E., Dozio, E., Barassi, A., Sammarco, G., Tacchini, L., Marrocco-Trischitta, M. M., et al. (2017). A Pilot Observational Study on Magnesium and Calcium Imbalance in Elderly Patients with Acute Aortic Dissection. Immun. Ageing 14, 1. doi:10.1186/s12979-016-0083-y
Wang, W., Duan, W., Xue, Y., Wang, L., Liu, J., Yu, S., et al. (2014). Clinical Features of Acute Aortic Dissection from the Registry of Aortic Dissection in China. J. Thorac. Cardiovasc. Surg. 148 (6), 2995–3000. doi:10.1016/j.jtcvs.2014.07.068
Wang, X., Zhang, H., Cao, L., He, Y., Ma, A., and Guo, W. (2020). The Role of Macrophages in Aortic Dissection. Front. Physiol. 11, 54. doi:10.3389/fphys.2020.00054
Wang, Y.-d., Liu, Z.-j., Ren, J., and Xiang, M.-X. (2018). Pharmacological Therapy of Abdominal Aortic Aneurysm: An Update. Cvp 16 (2), 114–124. doi:10.2174/1570161115666170413145705
Wolford, B. N., Hornsby, W. E., Guo, D., Zhou, W., Lin, M., Farhat, L., et al. (2019). Clinical Implications of Identifying Pathogenic Variants in Individuals with Thoracic Aortic Dissection. Circ. Genom Precis. Med. 12 (6), e002476. doi:10.1161/CIRCGEN.118.002476
Wu, X.-w., Li, G., Cheng, X.-b., Wang, M., Wang, L.-l., Wang, H.-h., et al. (2021). Association of Angiotensin II Type 1 Receptor Agonistic Autoantibodies with Outcomes in Patients with Acute Aortic Dissection. JAMA Netw. Open 4 (10), e2127587. doi:10.1001/jamanetworkopen.2021.27587
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yuan, S. M., and Wu, N. (2018). Aortic α-smooth Muscle Actin Expressions in Aortic Disorders and Coronary Artery Disease: An Immunohistochemical Study. Anatol. J. Cardiol. 19 (1), 11–16. doi:10.14744/AnatolJCardiol.2017.7839
Zhao, G., Lu, H., Chang, Z., Zhao, Y., Zhu, T., Chang, L., et al. (2021). Single-cell RNA Sequencing Reveals the Cellular Heterogeneity of Aneurysmal Infrarenal Abdominal Aorta. Cardiovasc Res. 117 (5), 1402–1416. doi:10.1093/cvr/cvaa214
Keywords: aortic dissection, single-cell RNA sequencing, gene ontology enrichment analysis, smooth muscle cells, KEGG enrichment analysis
Citation: Xu C, Liu X, Fang X, Yu L, Lau HC, Li D, Liu X, Li H, Ren J, Xu B, Jiang J, Tang L and Chen X (2022) Single-Cell RNA Sequencing Reveals Smooth Muscle Cells Heterogeneity in Experimental Aortic Dissection. Front. Genet. 13:836593. doi: 10.3389/fgene.2022.836593
Received: 15 December 2021; Accepted: 16 June 2022;
Published: 11 August 2022.
Edited by:
Wei Xu, Texas A&M University Corpus Christi, United StatesReviewed by:
Frauke Seemann, Texas A&M University Corpus Christi, United StatesChunhou Zheng, Anhui University, China
Copyright © 2022 Xu, Liu, Fang, Yu, Lau, Li, Liu, Li, Ren, Xu, Jiang, Tang and Chen. 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: Lijiang Tang, emp5eXRhbmdAMTI2LmNvbQ==; Xiaofeng Chen, eGlhb2YuY2hlbkB5YWhvby5jb20=
†These authors have contributed equally to this work