- 1Vascular and Endovascular Surgery, Huangyan Hospital Affiliated to Wenzhou Medical University, Taizhou First People’s Hospital, Taizhou, Zhejiang, China
- 2Department of Emergency Intensive Care Unit, Huangyan Hospital Affiliated to Wenzhou Medical University, Taizhou First People’s Hospital, Taizhou, Zhejiang, China
Introduction: Aortic dissection (AD) is often fatal, and its pathogenesis involves immune infiltration and pyroptosis, though the molecular pathways connecting these processes remain unclear. This study aimed to investigate the role of immune infiltration and pyroptosis in AD pathogenesis using bioinformatics analysis.
Methods: Two Gene Expression Omnibus datasets and a Gene Cards dataset of pyroptosis-related genes (PRGs) were utilized. Immunological infiltration was assessed using CIBERSORT, and AD diagnostic markers were identified through univariate logistic regression and least absolute shrinkage and selection operator regression. Interaction networks were constructed using STRING, and weighted gene correlation network analysis (WGCNA) was employed to identify important modules and essential genes. Single-sample gene set enrichment analysis determined immune infiltration, and Pearson correlation analysis assessed the association of key genes with infiltrating immune cells.
Results: Thirty-one PRGs associated with inflammatory response, vascular epidermal growth factor receptor, and Rap1 signaling pathways were identified. WGCNA revealed seven important genes within a critical module. CIBERSORT detected immune cell infiltration, indicating significant changes in immune cell infiltration and pyroptosis genes in AD and their connections.
Discussion: Our findings suggest that key PRGs may serve as indicators for AD or high-risk individuals. Understanding the role of pyroptosis and immune cell infiltration in AD pathogenesis may lead to the development of novel molecular-targeted therapies for AD.
Conclusion: This study provides insights into the molecular mechanisms underlying AD pathogenesis, highlighting the importance of immune infiltration and pyroptosis. Identification of diagnostic markers and potential therapeutic targets may improve the management of AD and reduce associated morbidity and mortality.
1 Introduction
Aortic dissection (AD) is a fatal condition for which there is currently no reliable drug treatment (Yang et al., 2020). Various risk factors for AD have been identified, including hypertension, atherosclerosis, and hypercholesterolemia (Hibino et al., 2022). In addition, smoking, drug abuse (especially cocaine), pregnancy, history of cardiac surgery, history of AD, aortic constriction, aortic aneurysm, and aortic mitral valve are other factors known to predispose an individual to AD (Parve et al., 2017). Typical AD is characterized by an endoluminal flap that divides the actual lumen from the fake lumen (Bossone et al., 2018). The Stanford approach distinguishes between two types of AD, types A and B, based on the location and degree of the trapping. Type A entrapment involves the ascending aorta (and may spread to the descending aorta), whereas type B entrapment exclusively affects the descending aorta. This classification is essential to determine the appropriate treatment strategy (Clouse et al., 2004). In the United States and Europe, the prevalence of AD ranges between 3.5 and 7.2 per 100,000 persons (Smedberg et al., 2020). The incidence of AD has shown a continuous increase in recent years, with a 50% and 30% increase in the annual incidence of AD reported for men and women, respectively (Olsson et al., 2006). Although the in-hospital mortality rate for AD has decreased significantly from 31% to 22% and the overall 5-year survival rate increased from 5% to 32%, the median survival of patients with AD is still only 3 days (Clouse et al., 2004; Bossone and Eagle, 2021). The high mortality rate of AD cannot be disregarded; thus, further research into the underlying factors affecting the condition will assist in our comprehension of its pathogenic processes, guiding clinical diagnosis and treatment and improving clinical prognosis.
Immune infiltration, as a key immune system response, is crucial in both the formation and remodeling phases of AD. Immune infiltration has been observed within the arterial wall of AD, and leukocyte cell adhesion, monocyte migration, and myeloid leukocyte migration are significantly increased in AD (Li et al., 2022). The monocyte-macrophage system plays a critical role in the immunological inflammatory response during the development of AD (Gao et al., 2021). Neutrophil infiltration is currently considered to predominate in the acute phase of the disease (Yoshida et al., 2019), whereas macrophages appear after 1 day and peak at 2–7 days (Xu and Burke, 2013). In an immune-high group of patients with AD, the percentages of infiltrating CD8+ T cells and M1-type macrophages were found to be considerably higher than those in the immune-low group (Li et al., 2022). Substantial lymphoplasmacytic infiltration has also been reported in AD (Uchida et al., 2018). Activated T cells and macrophages may enhance smooth muscle cell elimination (He et al., 2006), increased vascular inflammation (Lian et al., 2019), increased elastic fiber and extracellular matrix degradation (Gomez et al., 2016), and vascular smooth muscle cell death (Wu et al., 2022). These processes eventually lead to aortic dilatation and rupture. In contrast, immunological responses modulate AD arterial wall remodeling, which involves complicated interactions between cells and immune inflammatory factors. As a result, we suspect that immune infiltration contributes significantly to the onset of AD through controlling inflammation as well as the deterioration and remodeling of the artery walls.
Recent evidence has also pointed to a role of pyroptosis, a form of inflammatory cell death, in the development of AD (Zhang X. et al., 2020). AD occurs when inflammatory vesicles are activated in response to a stimulatory danger signal (Duan et al., 2020), leading to the release of pro-inflammatory factors to ultimately induce pyroptosis (Sun et al., 2021). Therefore, inhibition of the pyroptosis pathway offers a possible strategy to rescue AD (Duan et al., 2020). Pyroptosis is common in immune cells such macrophages, monocytes, and dendritic cells (Taabazuing et al., 2017). Pyroptotic cells are considered to be more inflammatory and immunogenic than apoptotic cells (Tsuchiya et al., 2019) because they release cellular contents such as inflammatory cytokines and damage-associated molecular patterns (Ruan et al., 2018). However, there is a lack of research on the potential correlations between pyroptosis and immune infiltration in AD.
To explore these associations, in this study, we used bioinformatics to analyze the Gene Expression Omnibus (GEO) database of genes associated with AD. We combined these datasets with a pyroptosis-related genes (PRGs) dataset to determine the probable functional processes and essential genes of pyroptosis. The link between PRG hub genes and invading immune cells was further investigated to elucidate the roles of pyroptosis and immunological processes during AD progression. A flow chart of the overall analysis process is provided in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g001.jpg)
FIGURE 1. Analysis flow chart. AD, aortic dissection; PRDEGs, pyroptosis-related differentially expressed genes; WGCNA, weighted gene correlation network analysis; PPI, protein-protein interaction network; LASSO, least absolute shrinkage and selection operator; SVM, support vector machine; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DO, Disease Ontology; ROC, receiver operating characteristic; ssGSEA, single-sample gene set enrichment analysis.
2 Materials and methods
2.1 Data acquisition and pre-processing
The GSE1534343 (AD = 10, control = 10) and GSE107844 (AD = 3, control = 3) gene expression profile datasets were obtained from the GEO database (Barrett et al., 2007) using the R program GEOquery (Davis and Meltzer, 2007). The species source of both datasets is human (Homo sapiens); the GSE153434 dataset was compiled from the HiSeq X Ten microarray sequencing platform (platform number GPL20795), and the tissue source is the ascending aorta, whereas the dataset GSE107844 was derived from the Illumina chip platform (platform number GPL20301), and the tissue source is the thoracic aorta. Details of the two datasets are provided in Table 1. For each dataset, platform annotation information was downloaded to convert probe names to gene names, and multiple expression results for specific genes were replaced with the mean of the expression values. The data were normalized using the normalizeBetweenArrays function (Further details about this function can be found at: https://web.mit.edu/∼r/current/arch/i386_linux26/lib/R/library/limma/html/05Normalization.html) in the limma package and merged with batch correction using the Combat function with the R package sva23 (More information about the Combat function can be found at: https://rdrr.io/bioc/sva/man/ComBat.html). After pre-processing and merging, the two datasets showed a similar distribution of gene expression values (Figures 2A, B), and the batch effects of the two datasets were corrected (Figures 2C, D). Finally, the GSE153434 and GSE107844 datasets were merged and used as the combined dataset for further analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g002.jpg)
FIGURE 2. Gene expression distribution before and after batch correction for the GSE153434 and GSE107844 datasets. (A) Box plots of gene expression distribution before normalization. (B) Box plots of gene expression distribution after normalization. (C) Principal components analysis (PCA) before batch correction. (D) PCA analysis after batch correction.
2.2 Differential gene expression analysis and acquisition of PRGs
The GeneCards database (http://www.genecards.org/) (Safran et al., 2010) incorporates gene-centric data from approximately 150 online sources, including genomic, transcriptional, protein biology, genetic, clinical, and functional data. By searching for “pyroptosis” in the GeneCards database, we obtained a total of 436 PRGs for subsequent analysis. Based on the grouping information in the combined dataset, we used the R package limma to obtain genes with differential expression between AD and controls according to log fold change (FC) > 1 (upregulated expression) or < −1 (downregulated expression) and adjusted p-value < 0.05.
2.3 Functional enrichment of PRGs
Gene Ontology (GO) enrichment (The gene ontology, 2019) is commonly used for large-scale functional enrichment studies of genes at various dimensions and levels, including biological process, molecular function, and cellular component (Ashburner et al., 2000). The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a popular database containing data related to genomes, biological processes, diseases, and medications (Kanehisa and Goto, 2000). The R package “clusterprofiler” (Yu et al., 2012; Wu et al., 2021) was used to perform GO functional annotation, pathway, and disease enrichment analyses, and the “pathview” package was used to display the KEGG pathways with significant alterations in AD (Luo and Brouwer, 2013). GO terms and pathways with p < 0.05 were deemed to be significantly enriched.
2.4 Immune cell analysis
CIBERSORT (Newman et al., 2019) is a web-based tool on the R platform for deconvolution of an expression matrix of human immune cell subtypes based on the principle of linear support vector regression. According to a gene expression signature collection comprising 22 known immune cell subtypes, CIBERSORT can identify the infiltration status of immune cells in sequenced samples. The CIBERSORT program was used in this study to evaluate the status of infiltrating immune cells in the combined data set, and Spearman correlation analysis was employed to calculate the interrelationships between different immune cells.
2.5 Gene set enrichment analysis (GSEA)
GSEA is a computational enrichment method developed at the Broad Institute (Subramanian et al., 2005) with corresponding analysis software and a gene set database (MSigDB), which is widely used to determine whether a set of predefined genes differs significantly between two biological states and to estimate changes in pathway and biological process activity in samples of expression data sets. The reference gene set “c2.all.v2022.1.Hs.symbols.gmt” was downloaded from the MSigDB database (Liberzon et al., 2015) based on the gene expression profile dataset, which was then enriched and visualized using the GSEA method included in the R package “clusterProfiler” to investigate the differences in biological processes between the AD and control samples. The following settings were used in this GSEA: seed number = 2022, calculations = 1000, at least 10 genes contained in each gene set, and a maximum number of genes featured of 500. The p-value was adjusted using Benjamini–Hochberg correction. The screening parameters for substantial enrichment were p < 0.05 and false discovery rate (q.value) < 0.05.
2.6 Single-factor logistic regression, least absolute shrinkage and selection operator (LASSO), and support vector machine (SVM) analyses
We used univariate logistic regression analysis, LASSO regression analysis, and the SVM-RFE technique for feature selection to search for AD diagnostic indicators. Initially, we used univariate logistic regression analysis to determine the relationship between the expression level of each differentially expressed gene in patients with AD, retaining genes with a p-value < 0.05. Subsequently, we utilized the SVM and LASSO (Friedman et al., 2010) algorithms of the glmnet package to downscale and select significant variables in the one-way logistic regression analysis. To assess the diagnostic performance of the model, we used the R package “pROC” (Robin et al., 2011) to plot the receiver operating characteristic (ROC) curve and computed the area under the curve (AUC) values.
2.7 Weighted gene correlation network analysis (WGCNA)
WGCNA (Langfelder and Horvath, 2008) seeks to find co-expressed gene modules, investigate the link between gene networks and phenotypes, and investigate the network’s core genes. The pickSoftTreshold function calculates the soft threshold, with 5 being the optimal soft threshold. A scale-free network is then constructed based on the soft threshold, followed by construction of a topology matrix and hierarchical clustering. The gene modules were identified by dynamic cutting using 50 as the minimum number of genes in the module, and eigengenes were computed. Inter-module correlations were built using the module eigengenes and hierarchical clustering, and Pearson correlation analysis was used to evaluate the relationships between modules and between modules and clinical characteristics. To identify important genes associated with pyroptosis in AD, genes from the most relevant modules for AD were intersected with the differentially expressed PRGs.
2.8 Construction of interaction networks
The STRING database (Szklarczyk et al., 2019), comprising 9.6 million proteins and 138 million protein–protein interactions (PPIs) from 2031 species, is used for searching known proteins and predicting PPIs. STRING includes experimental data, the results of text mining of PubMed abstracts, a synthesis of other database data results, and outcomes anticipated using bioinformatics approaches. The STRING database was used to build PPI networks for intersecting gene sets of modular and differentially expressed PRGs linked with AD, which were then visualized using Cytoscape (v3.9.1) software (Shannon et al., 2003). We also utilized Cytoscape’s Molecular Complex Detection (MCODE) tool to cluster and identify the essential modules in the PPI network. The microRNAs (miRNAs) linked with the identified hub genes were retrieved and intersected using starBase version 3.0 to investigate the association between hub genes and miRNAs (Huang et al., 2020). Lastly, Cytoscape software was used to visualize the mRNA–miRNA regulation network.
2.9 Molecular fractionation
Consensus clustering (Lock and Dunson, 2013) is a resampling-based algorithm that can be used to identify each member of a cluster and its subgroup number and for verifying the reasonableness of the generated clusters. Consistent clustering involves multiple iterations over subsamples of the dataset, indicating cluster stability and parameter decisions by using subsampling to induce sampling variability. We used the R package “ConsensusClusterPlus” (Wilkerson and Hayes, 2010) to perform consistent clustering of the dataset using the key genes associated with pyroptosis in AD to facilitate better differentiation between different subtypes of AD samples. In this process, 80% of the total samples were drawn in 1000 repetitions with the following settings: clusterAlg = “hc” and distance = “pearson.”
2.10 Immune infiltration correlation analysis
The quantity of certain infiltrating immune cells and the activity of specific immune responses can be estimated using single-sample GSEA (ssGSEA). According to published data of tumor immune infiltration, ssGSEA identified 28 gene sets for identifying distinct tumor-infiltrating immune cell types, containing numerous human immune cell subtypes (e.g., CD8+ T cells, dendritic cells, macrophages, regulatory T cells). Among ssGSEA methods, the R package GSVA protocol can yield enrichment scores for assessing the infiltration of each immune cell type within each sample. Therefore, we used the “GSVA” R package for ssGSEA on the immune cell infiltration gene set obtained from the literature (Zhang B. et al., 2020). The correlation of key genes associated with AD with infiltrating immune cells was determined by Pearson correlation analysis. Heat maps were created in R using the “ggplot2” visualization tool.
2.11 Statistical analysis
All data computations and statistical analyses were carried out using R programming (version 4.2.0, available at https://www.r-project.org/). The statistical significance of normally distributed variables was evaluated using an independent Student’s t-test, whereas differences between non-normally distributed variables were investigated using the Mann–Whitney U-test. Statistical significance was judged with a two-sided p-value < 0.05.
3 Results
3.1 Differential gene expression and enrichment analysis
We obtained 1398 differentially expressed genes in the AD group from analysis of the combined dataset (Figure 3A), 31 of which intersected with the PRGs, including 25 upregulated genes and 6 downregulated genes (Figure 3B). GO enrichment analysis showed that the differentially expressed PRGs in AD were mainly enriched in positive regulation of defense response, positive regulation of response to external stimulus, regulation of inflammatory process, and other biological processes; inflammasome complex, secretory granule membrane, secretory granule lumen as cellular components, and vascular endothelial growth factor receptor molecular function (Table 2). PRGs were linked to biological pathways such as the Rap1 signaling pathway (Figures 3D, F), rheumatoid arthritis (Figures 3D, G), and others, according to KEGG analysis. The PRGs were mostly associated with periodontal disease, bacterial infectious illness, and lung disease based on Disease Ontology enrichment (Figure 3E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g003.jpg)
FIGURE 3. Functional enrichment analysis of differentially expressed pyroptosis-related genes. (A) Volcano plot display of the differential gene expression analysis from the combined aortic dissection (AD) dataset. (B) Differentially expressed pyroptosis-related genes; orange circles indicate differentially expressed genes (DEGs), green circles indicate pyroptosis-related genes (PRGs), and the intersection indicates differentially expressed pyroptosis-related genes in AD. (C) GO enrichment pathway circle diagram. The diagram can be divided into two parts: the inner circle and outer circle. Each bar in the inner circle corresponds to an entry, and the height is inversely related to the adjusted p-value (p.adj). The color of the corresponding filled bar represents the Z-score value corresponding to the entry. (D) Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis; the horizontal coordinates are differentially expressed PRGs enriched to various entries and the vertical coordinate is the–log10 (p-value). Light blue indicates biological processes, red indicates cellular components, purple indicates molecular components, and yellow indicates enriched biological pathways. (E) Disease Ontology (DO) analysis; yellow nodes indicate the type of disease the differentially expressed PRGs are enriched in, and gray nodes indicate the genes enriched to different diseases. (F) Rap1 signaling pathway. (G) Rheumatoid arthritis pathway visualization results. Each node indicates a gene that plays an important role in the pathway, and the color of the node is determined by the log2 fold change (FC), with green indicating differentially downregulated genes and red indicating differentially upregulated genes.
3.2 Immune infiltration correlation analysis
The degree of infiltration of 22 immune cell types in different samples of the combined dataset is presented as a stacked bar graph in Figure 4A. Correlation analysis of infiltrating immune cells in AD revealed that M1 macrophages had a significant positive correlation with CD4+ memory activated T cells (r = 0.55), but had a significant positive correlation with naïve B cells (r = 0.64) and follicular helper T cells (r = −0.58) (Figure 4B). After removing immune cells with an immune abundance of zero (memory B cells and naïve CD4+ T cells), the Wilcoxon test showed that four immune cell types, including M0 macrophages, had significantly different infiltration levels between the control and AD groups from the combined dataset (p < 0.05). Neutrophils had considerably greater infiltration levels in the AD group, whereas B cells and M1 macrophages had significantly reduced infiltration levels in the AD group (Figure 4C), suggesting that these immune cells may play a key role in AD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g004.jpg)
FIGURE 4. Analysis of immune cell infiltration in the combined aortic dissection (AD) dataset of GSE153434 and GSE107844. (A) Bar graphs of 22 immune cell types stacked in different samples of the combined dataset, with different colored long bars representing different immune cells. (B) Correlation analysis between various immune cells in the dataset. The color represents the correlation strength, with darker orange indicating a stronger correlation. (C) Differences in the abundance of enrichment of 20 immune cells in the dataset; green indicates the control group, orange indicates the AD group, the horizontal axis indicates the 20 immune cells, and the vertical axis indicates the abundance of immune cell infiltration. (ns: p > 0.05, *: p < 0.05, **: p < 0.01, ***: p < 0.001, ****: p < 0.0001).
3.3 Gene set enrichment analysis (GSEA)
Gene Set Enrichment Analysis (GSEA) is conducted to elucidate the primary signaling pathways involved in Alzheimer’s Disease (AD). As depicted in Figure 5A, significant enrichment and upregulation (Normalized Enrichment Score = 1.944, nominal p-value = 0.031) of the KEGG pathway p53 signaling pathway are observed, indicating a strong overall association between apoptotic cell death and Thoracic Aortic Aneurysm and Dissection (TAAD). Furthermore, GSEA reveals significant enrichment of gene expression related to the cell cycle (Figure 5C), JAK-STAT signaling pathway (Figure 5E), and other crucial physiological pathways in AD (Table 3). Conversely, pathways such as Wnt signaling (Figure 5D) are downregulated in the AD group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g005.jpg)
FIGURE 5. Gene set enrichment analysis (GSEA). (A) Four main biological features of GSEA of genes among different subgroups in the aortic dissection (AD) group. The vertical coordinate is the gene set name. (B–E) Genes in the AD group were significantly enriched in the KEGG p53 signaling pathway (B), KEGG cell cycle (C), WikiPathways Wnt signaling (D), and KEGG JAK-STAT signaling pathway (E). Significant enrichment was assessed at p.adjust < 0.05 and FDR value (q.value) < 0.05.
3.4 Screening of key genes
Logistic regression of the 31 differentially expressed PRGs in AD is presented as a forest plot in Figure 6A. A total of 22 genes were screened as key genes (p < 0.05), including matrix metallopeptidase 1 (MMP1), basic helix-loop-helix family member E40 (BHLHE40), serpin family B member 1 (SERPINB1), vitamin D receptor (VDR), nuclear paraspeckle assembly transcript 1 (NEAT1), triggering receptor expressed on myeloid cells 1 (TREM1), formyl peptide receptor 2 (FPR2), and adenosine A3 receptor (ADORA3). The LASSO logistic regression algorithm was used to identify seven genes from the 22-gene set that could be used as diagnostic markers for disease control grouping (Figures 6B, C), and the SVM-RFE algorithm was used to classify these seven genes (Figures 6D, E). The gene markers generated by the two algorithms were then superimposed to yield five diagnostic marker genes (MMP1, BHLHE40, NEAT1, ADORA3, PPARG; Figure 6F) as candidate essential genes in AD for further investigation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g006.jpg)
FIGURE 6. Screening of key genes. (A) Single-factor logistic regression forest plot of pyroptosis-related genes from the merged dataset. (B) LASSO regression analysis. For each gene regression covariate, positive numbers are positively correlated, and negative numbers are negatively correlated. (C) Vertical coordinates are the evaluation index corresponding to each λ value, and the best covariate (λ) is selected. (D) Number of genes with the highest accuracy obtained by the SVM algorithm. © Number of genes with the lowest error rate obtained by the SVM algorithm. (E) The minimum number of genes with the lowest error rate obtained by the E.SVM algorithm. (F) Intersection of the genes screened by the SVM algorithm (green) and the genes screened by LASSO (orange). LASSO, least absolute shrinkage and selection operator; SVM, support vector machines.
3.5 Analysis of key genes
We then examined changes in the expression levels of the candidate key genes associated with AD (MMP1, BHLHE40, NEAT1, ADORA3, and PPARG) in the combined dataset, and discovered that MMP1, BHLHE40, NEAT1, and ADORA3 were all overexpressed in the AD group compared with their levels in the control group (Figures 7A–C, E), whereas PPARG showed reduced expression in the AD group as opposed to the control group (Figure 7D). In addition, ROC analysis showed that MMP1 (AUC = 0.953, Figure 7G), BHLHE40 (AUC = 0.959, Figure 7F), NEAT1 (AUC = 0.923, Figure 7H), and PPARG (AUC = 0.923, Figure 7J) all have good potential to differentiate AD patients from controls, and ADORA3 (AUC = 0.888, Figure 7I) also has moderate potential to distinguish AD, suggesting that these genes could be used as potential diagnostic markers for AD in the future.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g007.jpg)
FIGURE 7. Key genes analysis. (A–E) Differences in expression levels of BHLHE40 (A), MMP1 (B), NEAT1 (C), PPARG (D), and ADORA3© in the control (green) and aortic dissection (AD, orange). (F–J) Receiver operating characteristic (ROC) curves for BHLHE40 (F), MMP1 (G), NEAT1 (H), ADORA3 (I), and PPARG (J). The horizontal coordinate is 1—specificity and the vertical coordinate is sensitivity. AUC, area under the curve. (K) Functional similarity analysis of key genes. (ns: p > 0.05, *: p < 0.05, **: p < 0.01, ***: p < 0.001, ****: p < 0.0001); the closer the AUC of the ROC curve is to 1, the better the diagnosis prediction (AUC of 0.5–0.7 indicates low accuracy, AUC of 0.7–0.9 indicates moderate accuracy, and AUC > 0.9 indicates high accuracy).
For functional similarity of key genes (MMP1, BHLHE40, NEAT1, ADORA3, PPARG), we assessed the semantic similarity of GO terms, sets of GO terms, gene products, and gene clusters using the R package GOSemSim. Figure 7K shows a box-and-line graphic of these essential genes. Among the five critical genes, PPARG exhibited the highest functional similarity score with other key genes.
3.6 Profile of WGCNA analysis
To investigate the co-expression of genes, we constructed co-expression modules using the top 25% of the genes in the combined dataset mean expression profile using the R package WGCNA. The data were hierarchically grouped using the average technique in conjunction with the patients’ clinical information (Figure 8A). The soft threshold was set to 5 to create a scale-free network (Figure 8B), with R2 > 0.80 and good network connection. Ultimately, we identified 14 gene co-expression modules using the dynamic shear tree algorithm with each module represented by a distinct hue (Figure 8C). Correlation analysis of various modules with the AD group (Figure 8D) revealed that the turquoise module was most strongly and significantly positively correlated with the AD group (r = 0.94), which was therefore selected for further analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g008.jpg)
FIGURE 8. Weighted gene correlation network analysis (WGCNA) based on the expression spectrum dataset. (A) Sample clustering to detect abnormal samples; each branch indicates one sample and the red box in the sample information below indicates the category to which the sample belongs. (B) Scale-free topological model fit (left) and average connectivity (right) used to determine the best soft threshold. The horizontal coordinates of the left panel indicate the fitted soft threshold and the vertical coordinates indicate the level of fit (R2) in the scale-free topological model. The right panel horizontal coordinate indicates the fitted soft threshold and the vertical coordinate indicates the average connectivity between modules. (C) Dynamic shear clustering tree for different genes. The upper tree shows gene co-expression; each gene is represented by a leaf in the tree, and each module is represented by a trunk branch. The lower colored bars indicate the corresponding seven modules and are labeled with the indicated colors. (D) Heat map of correlations between different modules and aortic dissection (AD) patient groups. The x-axis is the sample classification, and the y-axis is the module color name. Red indicates a positive correlation, and dark blue color indicates a negative correlation. Color blocks indicate the Pearson correlation coefficient in the upper part and the p-value in the lower part.
3.7 PPI network
We intersected the genes in the turquoise module of WGCNA, which was closely related to the AD group, with the differentially expressed PRGs to obtain seven candidate hub genes for further analysis: BHLHE40, caspase 4 (CASP4), platelet and endothelial cell adhesion molecule 1 (PECAM1), pyruvate kinase M1/2 (PKM), serpin family B member 1 (SERPINB1), Toll-like receptor 2 (TLR2), and vascular endothelial growth factor A (VEGFA) (Figure 9A). Construction of the PPI network of these seven genes with the STRING database and Cytoscape (minimum required interaction score: medium confidence = 0.400; Figure 9B) revealed interactions among five proteins (CASP4, PECAM1, PKM, TLR2, VEGFA). The MCODE plugin identified a critical module of the PPI network comprising three proteins (PECAM1, TLR2, and VEGFA), as depicted in red in Figure 9B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g009.jpg)
FIGURE 9. Protein–protein interaction network analysis based on seven key genes. (A) Intersection of the turquoise module and differentially expressed pyroptosis-related genes; orange circles indicate genes in the turquoise module, and green circles indicate differentially expressed pyroptosis-related genes. (B) Key modules identified by MCODE plugin of Cytoscape software; red indicates the key modules found by MCODE. (C) mRNA–miRNA interaction network constructed based on starBase v3.0 database prediction and visualized by Cytoscape software. Green indicates key genes and orange indicates miRNAs. MCODE, Molecular Complex Detection; PRGs, pyroptosis-related genes; DEG, differentially expression gene.
We used mRNA–miRNA data from the starBase database to predict miRNAs interacting with the seven hub genes (BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, VEGFA) retaining only reliable mRNA–miRNA data pairs, and the interaction network was visualized with Cytoscape (Figure 9C). Our mRNA–miRNA interaction network included three hub genes (BHLHE40, PKM VEGFA), 30 miRNA molecules, and 38 pairs of mRNA–miRNA interaction connections (see Supplementary Table S2).
3.8 Immune infiltration correlation
We computed the correlations between major PRGs and infiltrating immune cells in AD to investigate the immunological microenvironment in AD. The connections between important (hub) PRGs (BHLHE40, CASP4, PKM, SERPINB1, TLR2, VEGFA) and numerous invading immune cells in AD were generally comparable; however, the correlations between PECAM1 expression and multiple immune infiltrating cells differed from those of the other genes (Figure 10A). There was a significant positive correlation between the expression of SERPINB1 (r = 0.79, Figure 10B) and PKM (r = 0.72, Figure 10E) and central memory CD4+ T cells, PECAM1 (r = 0.74, Figure 10C) had a significant positive correlation with type 1 T helper cells, and TLR2 (r = 0.74, Figure 10D) had a significant positive correlation with natural killer T cells. Furthermore, there was a significant negative correlation between PECAM1 and CD56dim natural killer cells (r = −0.73, Figure 10F) and central memory CD4+ T cells (r = −0.57, Figure 10G), between TLR2 and activated B cells (r = −0.55, Figure 10H), and between PKM and type 1 T helper cells (r = −0.54, Figure 10I). These results suggest that CD56dim natural killer cells, central memory CD4+ T cells, activated B cells, type 1 T helper cells, and natural killer T cells play a significant role in the emergence of AD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g010.jpg)
FIGURE 10. Immune infiltration correlations. (A) Heat map of the correlation analysis between key genes associated with pyroptosis in aortic dissection (AD) (vertical axis) and infiltrating immune cells (horizontal axis); a redder color indicates a stronger correlation between the gene and immune cells. *p < 0.05, **p < 0.01. (B–I) Scatter plot of the correlations between (B) SERPINB1 and central memory CD4 T cells, (E) PKM and central memory CD4 T cells, (C) PECAM1 and type 1 T helper cells, (D) TLR2 and natural killer T cells, (E) PKM and central memory CD4 T cells, (F) PCAM1 and CD56dim natural killer cells, (G) PCAM1 and central memory CD4 T cells, (H) TLR2 and activated B cells, and (I) PKM and type 1 T helper cells. The absolute values of correlation coefficients (r) in the correlation scatter plots above 0.8 indicate a strong correlation; r = 0.5–0.8 indicates a moderate correlation; r = 0.3–0.5 indicates a weak correlation; and r < 0.3 indicates a weak or no correlation.
3.9 Drug prediction and molecular Docking
We used expression datasets of hub genes related to pyroptosis in AD (BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, VEGFA) to categorize patients using into subtypes using unsupervised consistent clustering. In cluster 1, including two AD samples, and cluster 2, including 11 AD samples, we could distinguish between two distinct subtypes (Figure 11A). The cumulative distribution function (CDF) plot (Figure 11B) and the area under the CDF curve delta plot (Figure 11C) demonstrate that the best consistency clustering results for the AD dataset were attained when using k = 2 as the number of clusters for unsupervised clustering. Principal component analysis on the dataset expression matrix showed that the two subtype samples of the AD dataset could be clearly distinguished from each other (Figure 11D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g011.jpg)
FIGURE 11. Molecular typing. (A) Plots of consistent clustering (K = 2) for the aortic dissection (AD) dataset. (B,C) Plots of the cumulative distribution function (CDF) for consistent clustering (B) and the area under the CDF curve delta plot for different numbers of clusters in consistent clustering (C). (D) Principal component analysis for two AD subtypes (cluster 1 and cluster 2) in the AD dataset. (E) Comparison graph of the grouping of hub genes in different AD subtypes. (F) Heat map of the expression of hub genes in different AD subtypes. ns, not significant (p ≥ 0.05); *p < 0.05, **p < 0.01, ***p < 0.001.
We also analyzed the differences in expression of the seven hub genes (BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, VEGFA) between the two AD subtypes (cluster 1 and cluster 2) using the Wilcoxon rank sum test. Grouped comparison plots presented in Figure 11E show that the expression levels of BHLHE40, CASP4, and VEGFA were all significantly different between the two AD subtypes (cluster 1 and cluster 2), with higher expression in AD subtype 2 (cluster 2) than in AD subtype 1 (cluster 1). The heat maps in Figure 11F further demonstrate significant differences in expression of BHLHE40, CASP4, and VEGFA between the two AD subtypes.
3.10 Analysis of hub genes
Further analysis of the changes in expression levels of hub genes (BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, VEGFA) associated with pyroptosis in AD in the combined dataset revealed that BHLHE40, CASP4, PKM, SERPINB1, TLR2, and VEGFA were all more highly expressed in the AD group compared to the control group, whereas PECAM1 had a lower expression level in the AD group compared to the control group (Figure 12A). We calculated the semantic similarity between GO terms, sets of GO terms, gene products, and gene clusters of the hub genes using the GOSemSim R package for functional similarity analysis. The box line diagram of these genes is shown in Figure 12B, demonstrating that VEGFA had the highest functional similarity value with other key genes among the seven key genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1277818/fmolb-11-1277818-HTML/image_m/fmolb-11-1277818-g012.jpg)
FIGURE 12. Analysis of hub genes. (A) Differences in expression levels of BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, and VEGFA in the control (green) and aortic dissection (AD, orange) groups. (B) Functional similarity analysis of key genes. (C–I) Receiver operating characteristic (ROC) curves for CASP4 (C), PECAM1 (D), PKM (E), SERPINB1 (F), TLR2 (G), VEGFA (H), and BHLHE40 (I). The horizontal coordinate is 1—specificity and the vertical coordinate is sensitivity; AUC, area under the curve, with an AUC value closer to 1 indicating better diagnostic prediction ability. (ns: p > 0.05, *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).
ROC analysis showed that CASP4 (AUC = 1.0, Figure 12C), PECAM1 (AUC = 1.0, Figure 12D), PKM (AUC = 0.988, Figure 12E), SERPINB1 (AUC = 0.970, Figure 12F), TLR2 (AUC = 0.994, Figure 12G), VEGFA (AUC = 0.994, Figure 12H), and BHLHE40 (AUC = 0.959, Figure 12I) all have good potential to differentiate the AD group from normal controls, suggesting that these genes could serve as potential diagnostic markers for AD in the future.
4 Discussion
AD is a catastrophic event caused by a tear in the aortic intima or bleeding within the aortic wall (Nienaber et al., 2016). The incidence of AD has shown an increasing trend year by year (Howard et al., 2013). Approximately 20% of the patients suffering from AD die before reaching the hospital (Golledge and Eagle, 2008). In the absence of intervention, acute AD has a 90% mortality rate (YU et al., 2020). Surgery remains the gold-standard treatment approach (Erbel et al., 2014); however, early postoperative mortality is still high, at 9%–25% (Smedberg et al., 2020). Endovascular interventions offer new treatment options for descending AD (Erbel et al., 2014); however, this approach is also associated with several possible side effects such as acute lung damage, acute renal failure, ischemia of the spinal cord, and stroke (Fattori et al., 2008). Early diagnosis and treatment of AD can dramatically enhance a patient’s prognosis. Furthermore, research into the processes contributing to the development and progression of AD can facilitate earlier detection and treatment. Recent studies have found that the hallmarks of AD include infiltration of immune cells, hyperactivation of inflammation, and degradation of the extracellular matrix, leading to vascular remodeling and weakening of the aortic wall (Lian et al., 2019). Pyroptosis and immunological infiltration have also been found to play important roles in the pathophysiology of AD. However, there is a scarcity of studies investigating the relationship between immune cell infiltration and excessive inflammatory activation in AD.
In this study, we analyzed the GSE153434 and GSE107844 datasets and performed deep data mining of the AD-related literature and public databases using bioinformatics tools with the goal of providing insight into the roles of cellular pyroptosis and immune infiltration in the development of AD and their associations. Combining datasets identified 31 PRGs that are differentially expressed in AD. GO-KEGG enrichment analyses showed that these PRGs are mainly related to the immune response, inflammasome completion, and the Rap1 signaling pathway. Biologically relevant pathways such as the p53 signaling pathway, cell cycle, and JAK-STAT signaling pathway were upregulated in AD, whereas the Wnt signaling pathway was found to be downregulated. Among the PRGs, we found that MMP1, BHLHE40, NEAT1, and ADORA3 were highly expressed in the AD group, whereas the expression of PPARG was reduced compared with that of controls. These genes therefore have good potential to distinguish patients with AD from controls. Further analysis identified seven hub genes associated with AD (BHLHE40, CASP4, PECAM1, PKM, SERPINB1, TLR2, and VEGFA) that showed high diagnostic value. Among these hub PRGs, the expression of BHLHE40, CASP4, and VEGFA could also distinguish two AD subtypes characterized by dissection of the ascending aorta and thoracic aorta. Correlation analysis further showed enhanced infiltration of M0 macrophages and neutrophils, and decreased infiltration of naive B cells and M1 macrophages in AD samples compared with controls. The hub genes (BHLHE40, CASP4, PKM, SERPINB1, TLR2, VEGFA) were correlated with a variety of invading immune cells, suggesting the involvement of numerous immune cells, including CD56dim natural killer cells, central memory CD4+ T cells, activated B cells, type 1 T helper cells, and natural killer T cells, in the development of AD.
Our analysis identified multiple genes with good potential to differentiate AD, suggesting their potential as biomarkers for the diagnosis of AD. Among them, MMP1 is associated with signaling pathways of immune system cytokines that may promote persistent destruction of the aortic extracellular matrix and outer membrane degeneration (Zhang et al., 2014) to promote AD (Liao et al., 2018). BHLHE40 encodes a transcription factor with known regulatory effects on immune cells. BHLHE40 interferes with the cytokine production of CD4+ T cells and with the proliferation of macrophages and CD8+ T cells, promotes the production of interferon-gamma by natural killer T cells (Zafar et al., 2021), increases the recruitment of neutrophils, and promotes inflammatory and hypoxic responses (Wang L. et al., 2022). We discovered a positive link between BHLHE40 expression and a variety of immune cells, including CD4+ T cells, CD56dim natural killer cells, type 17 T helper cells, and memory B cells, and was inversely associated with activated B cells and T follicular helper cells. ADORA3 expression was found to be significantly increased in AD. ADORA3 has been reported to activate CD4+ T cells to release interleukin (IL)-10 (Durante et al., 2021), and to promote inflammatory response, apoptotic phagocytosis (Hodrea et al., 2012), and cell proliferation. PPARG, the only key gene that was downregulated in AD, is an anti-inflammatory gene that prevents the degeneration and death of the aortic stroma by inhibiting NF-κ-B-mediated pro-inflammatory responses (Wang M. et al., 2022) and regulates the role of the cardiovascular circadian rhythm (Aicher et al., 2021). CASP4 is an essential pyroptotic effector that facilitates the production of IL-1 and cell death, stimulates immune cell recruitment, and activates and causes mucosal inflammation (Wang et al., 2020). We found substantial positive correlations between CASP4 expression and the infiltration of CD4+ T cells, CD8+ T cells, macrophages, and natural killer T cells in AD. As a result, we may infer that CASP4 plays a significant role in immune cell activation in AD and is a crucial factor in pyroptosis. Platelets, monocytes, neutrophils, and certain subtypes of T cells express PECAM1 on their surface (Dasgupta et al., 2009), which has been suggested to play a role in leukocyte motility, angiogenesis, and integrin activation, thereby inhibiting macrophage-mediated phagocytosis (Brown et al., 2002). We found substantial negative correlations between PECAM1 expression and the infiltration of CD56dim natural killer cells and central memory CD4+ T cells in AD. TLR2 activation encourages the upregulation of inflammatory signaling pathways, promoting T cell/natural killer T cell activation and infiltration (Moonen et al., 2019; Imanishi et al., 2020), NF-κB activation, cytokine release, and an inflammatory response (Dutta et al., 2021) and further triggers apoptosis. In addition, we discovered a substantial negative correlation between activated B cells and TLR2 and a strong positive correlation between natural killer T cells and TLR2. VEGFA can induce endothelial cell proliferation, cell migration, apoptosis inhibition, and vascular permeabilization as a major contributor to both pathological and normal angiogenesis. VEGFA increases the amount of angiopoietin II and facilitates the recruitment of inflammatory cells, hence promoting inflammation (Glorioso et al., 2007). In the context of aortic dissection, there may be interactions and functional associations between VEGFA and PPARG. Regarding angiogenesis, VEGFA, as a crucial angiogenic factor, regulates angiogenesis and vascular permeability while inhibiting endothelial cell apoptosis (Melincovici et al., 2018). As a transcription factor, PPARG may modulate the expression of VEGFA (Blitek and Szymanska, 2019), thereby influencing angiogenesis and vascular permeability, further impacting the development and progression of aortic dissection. Additionally, in inflammation and immune responses, VEGFA (Chen et al., 2023) and PPAR (Toffoli et al., 2017) may mutually influence each other through the regulation of inflammation and immune responses, affecting the pathological processes of aortic dissection. Further research is warranted to deepen our understanding of their importance and mechanisms of action in this disease context. Our PPI network of key PRGs in AD showed an important module comprising interactions among PECAM1, TLR2, and VEGFA. We found a significant positive correlation between PKM expression and central memory CD4+ T cells and a significant negative correlation between PKM expression and type 1 T helper cells. PKM is engaged in several pathways of the innate immune system. Finally, SERPINB1 regulates the innate immune response, inflammation, and cellular homeostasis (Choi et al., 2019) and is significantly induced in the effector CD4 cell subpopulation (Hou et al., 2019). We identified a significant positive correlation between SERPINB1 expression and central memory CD4+ T cells in AD samples. Overall, our analysis demonstrated that these critical genes are involved in the control of inflammatory phenotype and immune cell activity, thereby having a significant impact on the pathogenesis and progression of AD. In particular, we found that BHLHE40, CASP4, and VEGFA can discriminate between two subtypes of AD (involvement of the thoracic or ascending aorta).
Analysis of the AD dataset demonstrated an assembly of inflammatory vesicles in AD, which stimulates the production of inflammatory cytokines as well as the death of inflammatory cells (i.e., pyroptosis) (Sharma and Kanneganti, 2021). Similarly, studies with animal models showed inflammatory vesicles in the aortic tract along with elevated serum levels of inflammatory cytokines, accompanied by pathological findings of immune infiltration, matrix degradation, and angiogenesis (Del Porto et al., 2014; Chen et al., 2022). Several signaling pathways have been implicated in the development of AD, in addition to increased inflammatory vesicles. The VEGFR signaling pathway plays an important role in aortic smooth muscle cell proliferation and migration (Wang et al., 2021). The use of VEGR tyrosine kinase inhibitors may be linked to the development of acute AD (Oshima et al., 2017; Ting and Lo, 2021). Similarly, the Rap1 signaling pathway is associated with cell migration, polarization, and proliferation (Cha et al., 2010). All of these cellular processes are implicated in the development of AD (Shah et al., 2019). Previous studies have also highlighted p53 signaling pathway-dependent proliferation inhibition, oxidative stress, and apoptosis among the potential mechanisms of AD (Wang et al., 2019; Wu et al., 2019). In addition, p53 plays a potential role in processes related to cell turnover during aortic intima-media degeneration (Ihling et al., 1999; Wang et al., 2019). The JAK-STAT pathway is extensively involved in aortic vascular injury, regulating cytokine expression and immune cell activation, and thus disease progression. This pathway stimulates the effector function of macrophages, promotes the differentiation of type 17 T helper lymphocytes, and enhances the expression of matrix metalloproteinases, ultimately leading to deterioration of the structural integrity of the vessel wall. Wnt signaling is essential for cell proliferation, differentiation, and migration. Wnt pathway inhibition enhances lipid and macrophage retention in the vessel wall as well as an increase in leukocyte-driven systemic inflammation (Borrell-pages et al., 2015). Molecules targeting the Wnt pathway attenuate intimal changes in the aorta caused by mechanical injury through the attenuation of Wnt signaling (Gay and Towler, 2017). Therefore, modulation of this common pathway might offer novel therapeutic interventions for inflammation-driven vasodilation in advanced aortic disease (Ijaz et al., 2016).
Among the multiple immune cells identified to be activated and infiltrated in AD, the infiltration of macrophages into the aortic wall is considered to be the main pathogenic mechanism of aortic injury (Guido et al., 2022). Amorphous macrophages (M0 macrophages) are driven to differentiate into a pro-inflammatory (M1) or anti-inflammatory (M2) phenotype (also known as macrophage polarization) (Xiao et al., 2020). M1 macrophages could increase and sustain the inflammatory response by secreting pro-inflammatory cytokines and promoting vascular injury (Sun et al., 2021). In addition, M0 macrophages play a key role in inflammation and immunity by recruiting more neutrophils and monocytes/macrophages and initiating further immune responses (Lian et al., 2019). The key PRGs identified in our study, including BHLHE40, MMP1, ADORA3, CASP4, and PECAM1, were closely associated with macrophage activation and infiltration. The levels of infiltration of M0 macrophages and neutrophil were considerably higher, whereas those of naïve B cells and M1 macrophages M1 infiltration levels were significantly lower in AD samples than in control samples, suggesting that neutrophil infiltration may dominate the acute phase of human AD. Consistent with this finding, previous studies showed that AD occurs when the expression of the outer membrane CXCL1/granulocyte colony-stimulating factor drives local neutrophil recruitment, activation, and infiltration (Yoshida et al., 2019), resulting in outer membrane inflammation, aortic dilatation, and rupture via IL-6 production (Anzai et al., 2015). We discovered a close correlation between the hub genes and immune cell infiltration, suggesting that CD56dim natural killer cells, central memory CD4+ T cells, activated B cells, type 1 T helper cells, and natural killer T cell immunity may play an important role in the development of AD.
However, this study has certain limitations to address. Despite performing a multifaceted bioinformatics analysis, further validation studies are needed to support the present findings, including in vivo and ex vivo investigations and potential clinical trials. In addition, we have not discovered the particular mechanism of action by which pyroptosis and immune cell infiltration contribute to AD; thus, more research is needed to understand the detailed process.
In conclusion, we discovered that pyroptosis and immune cell infiltration, along with their interactions, play crucial roles in the development of AD. Changes in the expression of genes associated with pyroptosis may encourage the infiltration of immune cells in various stages of AD development, whereas the pyroptosis of immune cells may result in the strong release of inflammatory factors, thereby exacerbating the damage in various regions of the aorta. MMP1, BHLHE40, NEAT1, PPARG, ADORA3, CASP4, and VEGFA may be useful diagnostic markers to distinguish AD from normal controls and facilitate early detection and intervention. Among these key genes, CASP4 and VEGFA could further distinguish between the two subtypes of AD. Overall, we expect that further research into the mechanisms of AD pyroptosis or immune cell infiltration will open up new avenues for the development of molecular targeted treatments for AD.
Impact statement
Aortic dissection (AD) is a life-threatening disorder that frequently leads in mortality. Recent research suggests that pyroptosis and immune cell infiltration play a role in the pathophysiology of AD; however, the precise relationship and molecular mechanisms mediating the links between immune cell infiltration and excessive inflammatory activation in AD remain unknown. In-depth studies into the mechanisms of AD pyroptosis and immune cell infiltration will pave the way for molecularly targeted AD treatment.
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 in the article/Supplementary Material.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
XG: Methodology, Software, Visualization, Writing–original draft. QC: Funding acquisition, Project administration, Writing–review and editing. YC: Conceptualization, Investigation, Writing–original draft. CM: Conceptualization, Formal Analysis, Investigation, Writing–original draft. JF: Data curation, Software, Writing–review and editing. FL: Resources, Supervision, Writing–review and editing.
Funding
The authors declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Medical Science and Technology Project of Zhejiang Province (grants 2020RC143 and 2023RC110) and the Science and Technology Plan Project of Taizhou (grant 23ywb44).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1277818/full#supplementary-material
References
Aicher, B. O., Zhang, J., Muratoglu, S. C., Galisteo, R., Arai, A. L., Gray, V. L., et al. (2021). Moderate aerobic exercise prevents matrix degradation and death in a mouse model of aortic dissection and aneurysm. Am. J. Physiol. Heart Circ. Physiol. 320 (5), H1786–H1801. doi:10.1152/ajpheart.00229.2020
Anzai, A., Shimoda, M., Endo, J., Kohno, T., Katsumata, Y., Matsuhashi, T., et al. (2015). Adventitial CXCL1/G-CSF expression in response to acute aortic dissection triggers local neutrophil recruitment and activation leading to aortic rupture. Circ. Res. 116 (4), 612–623. doi:10.1161/CIRCRESAHA.116.304918
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25 (1), 25–29. doi:10.1038/75556
Barrett, T., Troup, D. B., Wilhite, S. E., Ledoux, P., Rudnev, D., Evangelista, C., et al. (2007). NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 35, D760–D765. doi:10.1093/nar/gkl887
Blitek, A., and Szymanska, M. (2019). Regulation of expression and role of peroxisome proliferator-activated receptors (PPARs) in luminal epithelial and stromal cells of the porcine endometrium. Theriogenology 127, 88–101. doi:10.1016/j.theriogenology.2019.01.002
Borrell-Pages, M., Romero, J. C., and Badimon, L. (2015). LRP5 deficiency down-regulates Wnt signalling and promotes aortic lipid infiltration in hypercholesterolaemic mice. J. Cell Mol. Med. 19 (4), 770–777. doi:10.1111/jcmm.12396
Bossone, E., and Eagle, K. A. (2021). Epidemiology and management of aortic disease: aortic aneurysms and acute aortic syndromes. Nat. Rev. Cardiol. 18 (5), 331–348. doi:10.1038/s41569-020-00472-6
Bossone, E., Labounty, T. M., and Eagle, K. A. (2018). Acute aortic syndromes: diagnosis and management, an update. Eur. Heart J. 39 (9), 739–749d. doi:10.1093/eurheartj/ehx319
Brown, S., Heinisch, I., Ross, E., Shaw, K., Buckley, C. D., and Savill, J. (2002). Apoptosis disables CD31-mediated cell detachment from phagocytes promoting binding and engulfment. Nature 418 (6894), 200–203. doi:10.1038/nature00811
Cha, I., Lee, S. H., and Jeon, T. J. (2010). Chemoattractant-mediated Rap1 activation requires GPCR/G proteins. Mol. Cells 30 (6), 563–567. doi:10.1007/s10059-010-0153-5
Chen, Y., He, Y., Wei, X., and Jiang, D. S. (2022). Targeting regulated cell death in aortic aneurysm and dissection therapy. Pharmacol. Res. 176, 106048. doi:10.1016/j.phrs.2021.106048
Chen, Z., Wang, W., and Hua, Y. (2023). Identification and validation of BCL6 and VEGFA as biomarkers and ageing patterns correlating with immune infiltrates in OA progression. Sci. Rep. 13 (1), 2558. doi:10.1038/s41598-023-28000-9
Choi, Y. J., Kim, S., Choi, Y., Nielsen, T. B., Yan, J., Lu, A., et al. (2019). SERPINB1-mediated checkpoint of inflammatory caspase activation. Nat. Immunol. 20 (3), 276–287. doi:10.1038/s41590-018-0303-z
Clouse, W. D., Hallett, J. W., Schaff, H. V., Spittell, P. C., Rowland, C. M., Ilstrup, D. M., et al. (2004). Acute aortic dissection: population-based incidence compared with degenerative aortic aneurysm rupture. Mayo Clin. Proc. 79 (2), 176–180. doi:10.4065/79.2.176
Dasgupta, B., Dufour, E., Mamdouh, Z., and Muller, W. A. (2009). A novel and critical role for tyrosine 663 in platelet endothelial cell adhesion molecule-1 trafficking and transendothelial migration. J. Immunol. 182 (8), 5041–5051. doi:10.4049/jimmunol.0803192
Davis, S., and Meltzer, P. S. (2007). GEOquery: a bridge between the gene expression Omnibus (GEO) and BioConductor. Bioinforma. Oxf. Engl. 23 (14), 1846–1847. doi:10.1093/bioinformatics/btm254
Del Porto, F., Di Gioia, C., Tritapepe, L., Ferri, L., Leopizzi, M., Nofroni, I., et al. (2014). The multitasking role of macrophages in Stanford type A acute aortic dissection. Cardiology 127 (2), 123–129. doi:10.1159/000355253
Duan, H., Zhang, X., Song, R., Liu, T., Zhang, Y., and Yu, A. (2020). Upregulation of miR-133a by adiponectin inhibits pyroptosis pathway and rescues acute aortic dissection. Acta Biochim. Biophys. Sin. (Shanghai) 52 (9), 988–997. doi:10.1093/abbs/gmaa078
Durante, M., Squillace, S., Lauro, F., Giancotti, L. A., Coppi, E., Cherchi, F., et al. (2021). Adenosine A3 agonists reverse neuropathic pain via T cell-mediated production of IL-10. J. Clin. Invest. 131 (7), e139299. doi:10.1172/JCI139299
Dutta, D., Jana, M., Majumder, M., Mondal, S., Roy, A., and Pahan, K. (2021). Selective targeting of the TLR2/MyD88/NF-κB pathway reduces α-synuclein spreading in vitro and in vivo. Nat. Commun. 12 (1), 5382. doi:10.1038/s41467-021-25767-1
Erbel, R., Aboyans, V., Boileau, C., Bossone, E., Bartolomeo, R. D., Eggebrecht, H., et al. (2014). 2014 ESC Guidelines on the diagnosis and treatment of aortic diseases: document covering acute and chronic aortic diseases of the thoracic and abdominal aorta of the adult. The Task Force for the Diagnosis and Treatment of Aortic Diseases of the European Society of Cardiology (ESC). Eur. Heart J. 35 (41), 2873–2926. doi:10.1093/eurheartj/ehu281
Fattori, R., Tsai, T. T., Myrmel, T., Evangelista, A., Cooper, J. V., Trimarchi, S., et al. (2008). Complicated acute type B dissection: is surgery still the best option? a report from the International Registry of Acute Aortic Dissection. JACC Cardiovasc Interv. 1 (4), 395–402. doi:10.1016/j.jcin.2008.04.009
Friedman, J., Hastie, T., and Tibshirani, R. (2010). Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33 (1), 1–22. doi:10.18637/jss.v033.i01
Gao, H., Sun, X., Liu, Y., Liang, S., Zhang, B., Wang, L., et al. (2021). Analysis of hub genes and the mechanism of immune infiltration in stanford type a aortic dissection. Front. Cardiovasc Med. 8, 680065. doi:10.3389/fcvm.2021.680065
Gay, A., and Towler, D. A. (2017). Wnt signaling in cardiovascular disease: opportunities and challenges. Curr. Opin. Lipidol. 28 (5), 387–396. doi:10.1097/MOL.0000000000000445
Glorioso, N., Herrera, V. L., Bagamasbad, P., Filigheddu, F., Troffa, C., Argiolas, G., et al. (2007). Association of ATP1A1 and dear single-nucleotide polymorphism haplotypes with essential hypertension: sex-specific and haplotype-specific effects. Circ. Res. 100 (10), 1522–1529. doi:10.1161/01.RES.0000267716.96196.60
Golledge, J., and Eagle, K. A. (2008). Acute aortic dissection. Lancet 372 (9632), 55–66. doi:10.1016/S0140-6736(08)60994-0
Gomez, I., Ozen, G., Deschildre, C., Amgoud, Y., Boubaya, L., Gorenne, I., et al. (2016). Reverse regulatory pathway (H2S/PGE2/MMP) in human aortic aneurysm and saphenous vein varicosity. PLoS One 11 (6), e0158421. doi:10.1371/journal.pone.0158421
Guido, M. C., Lopes, N. M., Albuquerque, C. I., Tavares, E. R., Jensen, L., Carvalho, P. d. O., et al. (2022). Treatment with methotrexate associated with lipid core nanoparticles prevents aortic dilation in a murine model of marfan syndrome. Front. Cardiovasc Med. 9, 893774. doi:10.3389/fcvm.2022.893774
He, R., Guo, D. C., Estrera, A. L., Safi, H. J., Huynh, T. T., Yin, Z., et al. (2006). Characterization of the inflammatory and apoptotic cells in the aortas of patients with ascending thoracic aortic aneurysms and dissections. J. Thorac. Cardiovasc Surg. 131 (3), 671–678. doi:10.1016/j.jtcvs.2005.09.018
Hibino, M., Otaki, Y., Kobeissi, E., Pan, H., Hibino, H., Taddese, H., et al. (2022). Blood pressure, hypertension, and the risk of aortic dissection incidence and mortality: results from the J-SCH study, the UK biobank study, and a meta-analysis of cohort studies. Circulation 145 (9), 633–644. doi:10.1161/CIRCULATIONAHA.121.056546
Hodrea, J., Majai, G., Doro, Z., Zahuczky, G., Pap, A., Rajnavölgyi, É., et al. (2012). The glucocorticoid dexamethasone programs human dendritic cells for enhanced phagocytosis of apoptotic neutrophils and inflammatory response. J. Leukoc. Biol. 91 (1), 127–136. doi:10.1189/jlb.0511243
Hou, L., Rao, D. A., Yuki, K., Cooley, J., Henderson, L. A., Jonsson, A. H., et al. (2019). SerpinB1 controls encephalitogenic T helper cells in neuroinflammation. Proc. Natl. Acad. Sci. U. S. A. 116 (41), 20635–20643. doi:10.1073/pnas.1905762116
Howard, D. P., Banerjee, A., Fairhead, J. F., Perkins, J., Silver, L. E., Rothwell, P. M., et al. (2013). Population-based study of incidence and outcome of acute aortic dissection and premorbid risk factor control: 10-year results from the Oxford Vascular Study. Circulation 127 (20), 2031–2037. doi:10.1161/CIRCULATIONAHA.112.000483
Huang, H. Y., Lin, Y. C., Li, J., Huang, K. Y., Shrestha, S., Hong, H. C., et al. (2020). miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 48 (D1), D148–D54. doi:10.1093/nar/gkz896
Ihling, C., Szombathy, T., Nampoothiri, K., Haendeler, J., Beyersdorf, F., Uhl, M., et al. (1999). Cystic medial degeneration of the aorta is associated with p53 accumulation, Bax upregulation, apoptotic cell death, and cell proliferation. . Heart 82 (3), 286–293. doi:10.1136/hrt.82.3.286
Ijaz, T., Tilton, R. G., and Brasier, A. R. (2016). Cytokine amplification and macrophage effector functions in aortic inflammation and abdominal aortic aneurysm formation. J. Thorac. Dis. 8 (8), E746–E754. doi:10.21037/jtd.2016.06.37
Imanishi, T., Unno, M., Kobayashi, W., Yoneda, N., Akira, S., and Saito, T. (2020). mTORC1 signaling controls TLR2-mediated T-cell activation by inducing TIRAP expression. Cell Rep. 32 (3), 107911. doi:10.1016/j.celrep.2020.107911
Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids Res. 28 (1), 27–30. doi:10.1093/nar/28.1.27
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Lian, G., Li, X., Zhang, L., Zhang, Y., Sun, L., Zhang, L., 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
Liao, M., Zou, S., Bao, Y., Jin, J., Yang, J., Liu, Y., et al. (2018). Matrix metalloproteinases are regulated by MicroRNA 320 in macrophages and are associated with aortic dissection. Exp. Cell Res. 370 (1), 98–102. doi:10.1016/j.yexcr.2018.06.011
Liberzon, A., Birger, C., Thorvaldsdottir, H., Ghandi, M., Mesirov, J. P., and Tamayo, P. (2015). The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 1 (6), 417–425. doi:10.1016/j.cels.2015.12.004
Li, Z., Wang, J., Yu, Q., Shen, R., Qin, K., Zhang, Y., et al. (2022). Identification of immune-related gene signature in stanford type A aortic dissection. Front. Genet. 13, 911750. doi:10.3389/fgene.2022.911750
Lock, E. F., and Dunson, D. B. (2013). Bayesian consensus clustering. Bioinforma. Oxf. Engl. 29 (20), 2610–2616. doi:10.1093/bioinformatics/btt425
Luo, W., and Brouwer, C. (2013). Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinforma. Oxf. Engl. 29 (14), 1830–1831. doi:10.1093/bioinformatics/btt285
Melincovici, C. S., Bosca, A. B., Susman, S., Mărginean, M., Mihu, C., Istrate, M., et al. (2018). Vascular endothelial growth factor (VEGF) - key factor in normal and pathological angiogenesis. Rom. J. Morphol. Embryol. 59 (2), 455–467.
Moonen, C. G. J., Karlis, G. D., Schoenmaker, T., Forouzanfar, T., Loos, B. G., and de Vries, T. J. (2019). T cell proliferation is induced by chronically TLR2-stimulated gingival fibroblasts or monocytes. Int. J. Mol. Sci. 20 (24), 6134. doi:10.3390/ijms20246134
Newman, A. M., Steen, C. B., Liu, C. L., Gentles, A. J., Chaudhuri, A. A., Scherer, F., et al. (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 37 (7), 773–782. doi:10.1038/s41587-019-0114-2
Nienaber, C. A., Clough, R. E., Sakalihasan, N., Suzuki, T., Gibbs, R., Mussa, F., et al. (2016). Aortic dissection. Nat. Rev. Dis. Prim. 2, 16053. doi:10.1038/nrdp.2016.53
Olsson, C., Thelin, S., Stahle, E., Ekbom, A., and Granath, F. (2006). Thoracic aortic aneurysm and dissection: increasing prevalence and improved outcomes reported in a nationwide population-based study of more than 14,000 cases from 1987 to 2002. Circulation 114 (24), 2611–2618. doi:10.1161/CIRCULATIONAHA.106.630400
Oshima, Y., Tanimoto, T., Yuji, K., and Tojo, A. (2017). Association between aortic dissection and systemic exposure of vascular endothelial growth factor pathway inhibitors in the Japanese adverse drug event report database. Circulation 135 (8), 815–817. doi:10.1161/CIRCULATIONAHA.116.025144
Parve, S., Ziganshin, B. A., and Elefteriades, J. A. (2017). Overview of the current knowledge on etiology, natural history and treatment of aortic dissection. J. Cardiovasc Surg. (Torino) 58 (2), 238–251. doi:10.23736/S0021-9509.17.09883-4
Robin, X., Turck, N., Hainard, A., Tiberti, N., Lisacek, F., Sanchez, J. C., et al. (2011). pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinforma. 12, 77. doi:10.1186/1471-2105-12-77
Ruan, J., Xia, S., Liu, X., Lieberman, J., and Wu, H. (2018). Cryo-EM structure of the gasdermin A3 membrane pore. Nature 557 (7703), 62–67. doi:10.1038/s41586-018-0058-6
Safran, M., Dalah, I., Alexander, J., Rosen, N., Iny Stein, T., Shmoish, M., et al. (2010). GeneCards Version 3: the human gene integrator. Database (Oxford) 2010, baq020. doi:10.1093/database/baq020
Shah, S., Brock, E. J., Ji, K., and Mattingly, R. R. (2019). Ras and Rap1: a tale of two GTPases. Semin. Cancer Biol. 54, 29–39. doi:10.1016/j.semcancer.2018.03.005
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Sharma, B. R., and Kanneganti, T. D. (2021). NLRP3 inflammasome in cancer and metabolic diseases. Nat. Immunol. 22 (5), 550–559. doi:10.1038/s41590-021-00886-5
Smedberg, C., Steuer, J., Leander, K., and Hultgren, R. (2020). Sex differences and temporal trends in aortic dissection: a population-based study of incidence, treatment strategies, and outcome in Swedish patients during 15 years. Eur. Heart J. 41 (26), 2430–2438. doi:10.1093/eurheartj/ehaa446
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Sun, R., Zhou, Y., and Cui, Q. (2021). Comparative analysis of aneurysm subtypes associated genes based on protein-protein interaction network. BMC Bioinforma. 22 (1), 587. doi:10.1186/s12859-021-04513-w
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1), D607–D13. doi:10.1093/nar/gky1131
Taabazuing, C. Y., Okondo, M. C., and Bachovchin, D. A. (2017). Pyroptosis and apoptosis pathways engage in bidirectional crosstalk in monocytes and macrophages. Cell Chem. Biol. 24 (4), 507–514. doi:10.1016/j.chembiol.2017.03.009
The Gene Ontology, C. (2019). The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 47 (D1), D330–D338. doi:10.1093/nar/gky1055
Ting, J., and Lo, Z. J. (2021). Vascular endothelial growth factor tyrosine kinase inhibitor targeted therapy: a potential cause of an acute aortic dissection lesser known to the emergency physician. BMJ Case Rep. 14 (10), e245653. doi:10.1136/bcr-2021-245653
Toffoli, B., Gilardi, F., Winkler, C., Soderberg, M., Kowalczuk, L., Arsenijevic, Y., et al. (2017). Nephropathy in Pparg-null mice highlights PPARγ systemic activities in metabolism and in the immune system. PLoS One 12 (2), e0171474. doi:10.1371/journal.pone.0171474
Tsuchiya, K., Nakajima, S., Hosojima, S., Thi Nguyen, D., Hattori, T., Manh Le, T., et al. (2019). Caspase-1 initiates apoptosis in the absence of gasdermin D. Nat. Commun. 10 (1), 2091. doi:10.1038/s41467-019-09753-2
Uchida, T., Hamasaki, A., Kuroda, Y., Sadahiro, M., Tamazawa, N., Ohe, R., et al. (2018). Immunoglobulin G subclass 4-related lymphoplasmacytic thoracic aortitis in a patient with acute type A aortic dissection. Ann. Thorac. Cardiovasc Surg. 24 (4), 208–210. doi:10.5761/atcs.cr.17-00111
Wang, K., Sun, Q., Zhong, X., Zeng, M., Zeng, H., Shi, X., et al. (2020). Structural mechanism for GSDMD targeting by autoprocessed caspases in pyroptosis. Cell 180 (5), 941–955. doi:10.1016/j.cell.2020.02.002
Wang, L., Liu, Y., Dai, Y., Tang, X., Yin, T., Wang, C., et al. (2022a). Single-cell RNA-seq analysis reveals BHLHE40-driven pro-tumour neutrophils with hyperactivated glycolysis in pancreatic tumour microenvironment. Gut 72, 958–971. doi:10.1136/gutjnl-2021-326070
Wang, M., Cui, B., Gong, M., Liu, Q., Zhuo, X., Lv, J., et al. (2022b). Arctium lappa leaves based on network pharmacology and experimental validation attenuate atherosclerosis by targeting the AMPK-mediated PPARG/LXRα pathway. Biomed. Pharmacother. 153, 113503. doi:10.1016/j.biopha.2022.113503
Wang, W., Liu, Q., Wang, Y., Piao, H., Li, B., Zhu, Z., et al. (2019). Verification of hub genes in the expression profile of aortic dissection. PLoS One 14 (11), e0224922. doi:10.1371/journal.pone.0224922
Wang, Z., Zhuang, X., Chen, B., and Wei, M. (2021). Osteoglycin knockdown promotes vascular smooth muscle cell proliferation and migration in aortic dissection via the VEGF/VEGFR2 axis. Mol. Med. Rep. 23 (1), 65. doi:10.3892/mmr.2020.11703
Wilkerson, M. D., and Hayes, D. N. (2010). ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinforma. Oxf. Engl. 26 (12), 1572–1573. doi:10.1093/bioinformatics/btq170
Wu, Q., Hong, J., Wang, Z., Hu, J., Chen, R., Hu, Z., et al. (2019). Abnormal ribosome biogenesis partly induced p53-dependent aortic medial smooth muscle cell apoptosis and oxidative stress. Oxid. Med. Cell Longev. 2019, 7064319. doi:10.1155/2019/7064319
Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. . Innov. (Camb) 2 (3), 100141. doi:10.1016/j.xinn.2021.100141
Wu, Y., Jiang, D., Zhang, H., Yin, F., Guo, P., Zhang, X., et al. (2022). N1-Methyladenosine (m1A) regulation associated with the pathogenesis of abdominal aortic aneurysm through YTHDF3 modulating macrophage polarization. Front. Cardiovasc Med. 9, 883155. doi:10.3389/fcvm.2022.883155
Xiao, J., Wei, Z., Chen, X., Chen, W., Zhang, H., Yang, C., et al. (2020). Experimental abdominal aortic aneurysm growth is inhibited by blocking the JAK2/STAT3 pathway. Int. J. Cardiol. 312, 100–106. doi:10.1016/j.ijcard.2020.03.072
Xu, L., and Burke, A. (2013). Acute medial dissection of the ascending aorta: evolution of reactive histologic changes. Am. J. Surg. Pathol. 37 (8), 1275–1282. doi:10.1097/PAS.0b013e318294adc3
Yang, K., Ren, J., Li, X., Wang, Z., Xue, L., Cui, S., et al. (2020). Prevention of aortic dissection and aneurysm via an ALDH2-mediated switch in vascular smooth muscle cell phenotype. Eur. Heart J. 41 (26), 2442–2453. doi:10.1093/eurheartj/ehaa352
Yoshida, S., Yamamoto, M., Aoki, H., Fukuda, H., Akasu, K., Takagi, K., et al. (2019). STAT3 activation correlates with adventitial neutrophil infiltration in human aortic dissection. Ann. Vasc. Dis. 12 (2), 187–193. doi:10.3400/avd.oa.19-00007
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yu, X., Suki, B., and Zhang, Y. (2020). Avalanches and power law behavior in aortic dissection propagation. Sci. Adv. 6 (21), eaaz1173. doi:10.1126/sciadv.aaz1173
Zafar, A., Ng, H. P., Kim, G. D., Chan, E. R., and Mahabeleshwar, G. H. (2021). BHLHE40 promotes macrophage pro-inflammatory gene expression and functions. FASEB J. 35 (10), e21940. doi:10.1096/fj.202100944R
Zhang, B., Wu, Q., Li, B., Wang, D., Wang, L., and Zhou, Y. L. (2020b). m(6)A regulator-mediated methylation modification patterns and tumor microenvironment infiltration characterization in gastric cancer. Mol. Cancer 19 (1), 53. doi:10.1186/s12943-020-01170-0
Zhang, X., Li, F., Wang, W., Ji, L., Sun, B., Xiao, X., et al. (2020a). Macrophage pyroptosis is mediated by immunoproteasome subunit β5i (LMP7) in abdominal aortic aneurysm. Biochem. Biophys. Res. Commun. 533 (4), 1012–1020. doi:10.1016/j.bbrc.2020.09.082
Keywords: aortic dissection, bioinformatics, pyroptosis, immune infiltration, database
Citation: Ge X, Cai Q, Cai Y, Mou C, Fu J and Lin F (2024) Roles of pyroptosis and immune infiltration in aortic dissection. Front. Mol. Biosci. 11:1277818. doi: 10.3389/fmolb.2024.1277818
Received: 15 August 2023; Accepted: 21 February 2024;
Published: 19 March 2024.
Edited by:
Debipreeta Bhowmik, Emory University, United StatesReviewed by:
Ruheena Javed, University of New Mexico, United StatesXiangbo An, Dalian Medical University, China
Copyright © 2024 Ge, Cai, Cai, Mou, Fu and Lin. 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: Qiqi Cai, caiqiqi302@126.com