- Research article
- Open access
- Published:
Autophagy-related genes analysis reveals potential biomarkers for prediction of the impaired walking capacity of peripheral arterial disease
BMC Medicine volume 21, Article number: 186 (2023)
Abstract
Background
The role of autophagy and autophagy-related genes in peripheral arterial disease (PAD) remains unknown and may be of diagnostic and prognostic value. The aim of this study is to investigate the relationship between autophagy and PAD, and identify potential diagnostic or prognostic biomarkers for medical practice.
Methods
Differentially expressed autophagy-related genes in PAD were explored from GSE57691 and validated in our WalkByLab registry participants by quantitative real-time polymerase chain reaction (qRT-PCR). The level of autophagy in peripheral blood mononuclear cells (PBMCs) of WalkByLab participants was assessed by analyzing autophagic marker proteins (beclin-1, P62, LC3B). Single sample gene set enrichment analysis (ssGSEA) was used to evaluate the immune microenvironment within the artery wall of PAD patients and healthy persons. Chemokine antibody array and enzyme-linked immunosorbent assay were used to assess the chemokines in participants’ plasma. Treadmill testing with Gardner protocol was used to evaluate participants’ walking capacity. Pain-free walking distance, maximum walking distance, and walking time were recorded. Finally, a nomogram model based on logistic regression was built to predict impaired walking performance.
Results
A total of 20 relevant autophagy-related genes were identified, and these genes were confirmed to be expressed at low levels in our PAD participants. Western blotting demonstrated that the expression of autophagic marker proteins beclin-1 and LC3BII were significantly reduced in PAD patients’ PBMCs. ssGSEA revealed that most of the autophagy-related genes were strongly correlated with immune function, with the largest number of associated genes showing interaction between cytokine-and-cytokine receptors (CCR). In this context, the chemokines growth-related oncogene (GRO) and neutrophil activating protein2 (NAP2) are highly expressed in the plasma of WalkByLab PAD patients and were significantly negatively correlated with the walking distance assessed by Gardner treadmill testing. Finally, the plasma NAP2 level (AUC: 0.743) and derived nomogram model (AUC: 0.860) has a strong predictive potential to identify a poor walking capacity.
Conclusions
Overall, these data highlight both the important role of autophagy and autophagy-related genes in PAD and link them to vascular inflammation (expression of chemokines). In particular, chemokine NAP2 emerged as a novel biomarker that can be used to predict the impaired walking capacity in PAD patients.
Background
Atherosclerosis of the arteries supplying the legs is commonly labeled as lower extremity peripheral arterial disease (PAD). PAD is a common disorder and the third leading cause of atherosclerotic cardiovascular morbidity, following coronary artery disease (CAD) and stroke [1]. It can manifest as intermittent claudication when walking or, in the most severe form, critical limb ischemia. Limb ischemia defined as chronic ischemic rest pain, ulcers, or gangrene of the lower extremity, which is a major cause of limb amputation, and a harbinger of cardiovascular mortality in PAD [1, 2]. Age, tobacco use, diabetes, hypertension, hypercholesterolemia, and a sedentary lifestyle are the major risk factors for PAD [3].
Autophagy is a cellular recycling process in response to various stressors (e.g., oxidative stress, hypoxia, and starvation). It has been implicated in several fundamental biological processes, including aging, immunity, development, tumorigenesis, vascular disease, cell death and differentiation [4, 5]. In the cardiovascular system, autophagy is a key regulator of homeostasis, responding to physiological and pathophysiological stimuli. Although recycling of cellular organelles is generally viewed as a beneficial process, insufficient and excessive levels of autophagy can lead to premature cell death (apoptosis) [6]. To date, there is increasing evidence that basal autophagy plays an essential role in protecting endothelial cells and smooth muscle cells from cell death and the development of vascular disease, particularly heart failure and atherosclerosis [5, 7]. Current research suggest that autophagy is also a fundamental process for cardiovascular homeostasis, health, and aging [8]. To the best of our knowledge, despite the protective role of autophagy in various diseases, its role in the vascular system is poorly understood and is entirely unknown in the context of PAD [9]. Exploring and uncovering potential autophagy-related genes and the autophagy level in PAD as well as a functional linkage with relevant molecular processes of vascular pathophysiology may provide potential biomarkers for assessing and monitoring disease risk and prognosis.
Recently, Biros et al. [10, 11] uploaded a gene expression omnibus (GEO) dataset GSE57691, which revealed the differentially expressed genes in the artery wall tissue from 49 patients with abdominal aortic aneurysm, 9 patients with chronic lower limb ischemia and 10 healthy controls. In the present study, we reevaluated this dataset for potential implications for lower limb ischemia, which served as the basis for an extensive new analysis. First, we explored the differentially expressed autophagy-related genes in PAD patients from the GSE57691 data set and validated target genes expression in peripheral blood samples from our Lauflabor (WalkByLab) registry participants. Then protein–protein interaction (PPI) analyses, Gene Ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis were performed to investigate the biological functions of these autophagy-related genes. Subsequently, western blot was performed to assess the level of autophagy in peripheral blood mononuclear cells (PBMCs) of WalkByLab participants by analyzing the degradation of autophagic marker proteins (beclin-1, P62, and LC3B). Here, a `pro-inflammatory´ phenotype and immune activation is an important feature of PAD [12]. Therefore, we next performed single sample gene set enrichment analysis (ssGSEA) to investigate the immune microenvironment of the artery wall of PAD patients and link molecular inflammatory processes with autophagy-related genes. Finally, we analyzed and validated potential biomarkers for predicting impaired walking capacity as assessed by Gardner treadmill testing used in clinical practice. In particular, the Gardner treadmill testing is regarded as the gold standard for physical performance in the diagnoses of PAD status, but it requires appropriate equipment, time, and a professionally trained team, which is why it is rarely used in clinical practice. In summary, the aim of this study was to (1) gain insight into the relationship between autophagy effector molecules and PAD, (2) link autophagy to inflammatory processes in the development of PAD, and (3) identify novel biomarkers that can reliably predict walking ability and PAD status in clinical practice.
Methods
Autophagy-related genes datasets and microarray data
A total of 232 genes were obtained from The Human Autophagy Database (http://www.autophagy.lu/index.html). The mRNA expression profile dataset of GSE57691 was downloaded from GEO (http://www.ncbi.nlm.nih.gov/geo/). GSE57691 is in GPL10558 platform (Illumina HumanHT-12 V4.0 expression beadchip), which included a gene expression analysis of artery wall tissue from 49 patients with abdominal aortic aneurysm, 9 patients with chronic lower limb ischemia and 10 healthy controls.
Differentially expressed analysis of autophagy-related genes
The normalized expression matrix of microarray data was downloaded from the GSE57691 dataset. Then the probes were annotated with the annotation files from the dataset. The repeatability of data in GSE57691 was verified by principal component analysis (PCA). The “limma” package of R software was used to identify the differentially expressed autophagy-related genes. Genes with an adjusted P-value < 0.05 and |log2fold change (FC)|≥ 1 were considered as differentially expressed genes. The heatmap and volcano plot were conducted using “heatmap” and “ggplot2” packages of R software.
PPI analysis and correlation analysis of the differentially expressed autophagy-related genes
PPI analysis of differentially expressed autophagy-related genes was analyzed using STRING database (https://string-db.org/) and Cytoscape software (version 3.8.2). The correlation analysis of the differentially expressed autophagy-related genes was identified using Spearman correlation in the “corrplot” package of R software.
GO and KEGG pathway enrichment analysis of autophagy-related genes
GO and KEGG pathway enrichment analysis were conducted in R software using the package “GO plot”. The GO analysis consisted of cellular component (CC), biological process (BP) and molecular function (MF).
Exploration of immune microenvironment
The single sample gene set enrichment analysis (ssGSEA) was performed to quantify the immune function levels. The normalized enrichment score that calculated from the ssGSEA was used in the ‘GSVA’ and ‘GSEABase’ R package [13, 14]. The annotated gene set file was obtained from the study of Jie-Ying Liang et al. [15]. We finally quantified the enrichment levels of the 16 immune cells and 13 immune-related pathways in each sample to reveal the immune function of PAD, and the results were expressed as immune scores.
The heat map of the involved immune cells and immune-related pathways was made using the "pheatmap" R package. The boxplot showed the level of immune infiltration in PAD group and healthy control group was made using the "ggpubr" and "reshape2" R package. A Spearman’s rank correlation analysis heatmap, which revealed the correlation of different immune cells infiltration levels and immune-related pathways, was performed using the "corrplot" R package. The correlation of autophagy-related genes expression and immune function was analyzed using the "ggcorrplot" R package.
Participants and their clinical characteristics
A total of 179 participants were selected from the WalkByLab-Registry database. These are the patients who presented to the WalkByLab center Brandenburg/Havel (Brandenburg Clinic, Brandenburg Medical School) Germany between October 2020 and August 2022. The WalkyByLab aims to interdisciplinarily screen, diagnose and follow-up patients with cardiovascular disease (www.lauflab.de). The WalkByLab register trial protocol was reviewed and approved by the ethical committee of the Cottbus Medical Association (Landesärztekammer Cottbus, study number of the ethics committee: AS 74(bB)/2018). The screening trial is performed in accordance with the principles of the declaration of Helsinki. The demographic and clinical characteristics of participates in this study are shown in Table 1. Propensity score matching analysis was performed to adjust for differences in baseline characteristics with the covariate of age using the "MatchIt" R package.
Isolation of peripheral blood mononuclear cells (PBMCs), RNA Extraction and quantitative real-time polymerase chain reaction (qRT-PCR)
The blood samples from each participant were processed to isolate PBMCs using Ficoll solution as described previously [16]. Briefly, 15 ml peripheral venous blood samples were collected in three vacutainer EDTA tubes. Next, 10 ml samples were diluted 1:1 with PBS, and carefully layered onto the Ficoll-Paque density gradient media (GE Healthcare) at a ratio of 4:3 and centrifuged at 400 g for 25 min. The middle layer containing PBMCs was collected and washed twice with PBS, and stored at -80 °C for RNA isolation. 5 ml samples were centrifuged at 1000 g for 10 min, the plasma was collected at -80 °C until use.
Total RNA was extracted from PBMCs by using the Trizol reagent (Thermo Fisher Scientific) according to the manufacturer’s instructions. The quality and quantity of total RNA were measured using NanoDrop-2000 spectrophotometer (NanoDrop Technologies). RNA samples with a 260/280 ratio between 1.8—2.0 were considered qualified. 1 μg of total RNA was reverse transcribed into cDNA using QuantiTect Reverse Transcription Kit (QIAGEN). The mRNA levels of the target gene were analyzed by real-time polymerase chain reaction using the LightCycler® 96 Real-Time PCR System (Roche). The relative expression of Mrna was calculated by 2−ΔΔCt method. All the primers were synthesized and purchased from Eurofins Genomics Germany GmbH (available in Additional file 1: Table S1).
Western blotting
PBMCs were lysed with RIPA lysis buffer and protein concentrations determined with the BCA Protein Assay kit according to the manufacturer’s protocol. Subsequently, 20 µg of protein/well were separated using 12% or 15% SDS‑PAGE. The proteins were subsequently transferred to nitrocellulose membranes at 250 mA for 15 to 60 min. Then membranes were blocked with Tris‑buffered saline containing 0.1% Tween 20 (TBST) and 5% non‑fat milk powder for 2 h at room temperature. After that, the membranes were incubated with primary antibodies overnight at 4˚C on a shaker set at a slow speed, then washed thrice with TBST and incubated with secondary antibody for 2 h at room temperature. After washing thrice, ECL substrate was added to the nitrocellulose membrane. The signal was detected using VWR Imaging Capture System. Band density was quantified with ImageJ software. All antibodies (P62, beclin1, LC3B, GAPDH, and HRP) were purchased from Abcam.
Human chemokine antibody array
The Human Chemokine Antibody Array-Membrane from Abcam was used to detect 38 chemokines in plasma according to the manufacturer’s instructions. In brief, chemokines receptor-coated membranes were incubated with diluted plasma overnight at 4 °C. After appropriate washing and antibody incubating using reagents provided by the manufacturers, membranes were imaged on the VWR Imaging Capture System. The densitometric analysis of chemokine spots was performed using the MicroArray Profile plugin for ImageJ.
Enzyme-linked immunosorbent assay (ELISA)
Plasma GRO and NAP2 levels were measured using the human GRO ELISA kit (RayBio) and human NAP2 ELISA kit (Abcam) according to the manufacturer’s instructions.
Flow mediated dilation (FMD) measurement
FMD measurement is regarded as the gold standard method for evaluation of endothelial function. Here, we used the AngioDefender (Everist Health) medical device to measure participants’ FMD, which allows for a precise, standardized, and automated value by oscillation technique and has been proven to show less measurement errors in compared with classical ultrasound techniques measurement.
Gardner treadmill testing
A treadmill and the Gardner treadmill protocol were used to assess patients’ walking capacity as previously reported [16]. Briefly, treadmill speed was maintained at 3.2 km/h with gradient increased 2 degrees every 2 min. This testing was performed under the supervision of a physician. The patient was secured with a drop stop device (safety bar with chest harness). Pain-free walking distance, maximum walking distance and walking time were recorded.
The construction of nomogram
Reduced walking capacity is one of the main clinical features of PAD. Here, we construct a nomogram to predict the impaired walking capacity. Logistic regression model with least absolute shrinkage and selection operator (LASSO) selection method was used to identify relevant factors, and then, the nomogram was constructed using the “rms” R packages. The calibration plot was used to evaluate the accuracy of the nomogram, and receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) value were used to evaluate the predictive efficacy.
Statistical analysis
The statistical analyses and plots were performed using R software (version 4.1.0) and R Studio. For independent-samples, Student’s t test was used. Continuous variables were compared using the t test if data distribution met the criteria for normality; otherwise, the Wilcoxon rank sum test was used. P < 0.05 was considered statistically significant.
Results
Differentially expressed autophagy-related genes in PAD
First, we evaluated the quality of the data, and performed a PCA. The data in GSE57691 exhibit satisfactory repeatability (Fig. 1A). In order to screen for differentially expressed autophagy-related genes, we next analyzed the expression of 232 autophagy-related genes in GSE57691, and 20 autophagy-related genes were identified using the criteria of adjusted P-value < 0.05 and |log2FC|≥ 1, including 1 up-regulated genes and 19 down-regulated genes, which were presented in heatmap and volcano plot (Fig. 1B and C).
PPI network and correlation analysis of the differentially expressed autophagy-related genes
To determine the interactions among differentially expressed autophagy-related genes, we performed PPI analysis. The proteins corresponding to the genes interact biochemically with each other in at least one and up to 10 different ways (Fig. 2A). To explore the expression correlation of these autophagy-related genes, Spearman correlation analysis was performed, yielding a potential interaction of 20 differentially expressed autophagy-related genes in the GSE57691 dataset (Fig. 2B).
Validation of differentially expressed autophagy-related genes in WalkByLab registry PAD patients
The expression levels of these differentially expressed autophagy-related genes were further validated by qRT-PCR analyses of PBMCs from our cohort of WalkByLab Brandenburg (PAD, n = 18; non-PAD controls, n = 18. Clinical characteristics available in Additional file 1: Table S2). Expression levels of BNIP3, EEF2, FOXO1, HSPA5, KLHL24, MAP1LC3B, PEA15, PTEN, RAB11A, RB1CC1, RPS6KB1, SAR1A, ST13, and VAMP3 genes were significantly lower in PAD patients as compared to non-PAD controls (Fig. 2C). The BAG3, CCL2, CSTB, and EEF2K genes tended to be less expressed in the PAD group but did not reach statistical significance. However, Expression levels of CX3CL1 and HSPB8 genes were very low, so we did not include them in the statistical analysis.
GO and KEGG enrichment analysis of the differentially expressed autophagy-related genes
To analyze the potential biological functions of these differentially expressed autophagy-related genes, we conducted GO and KEGG enrichment analysis. The results revealed that the most significant GO enriched terms are linked to the biological functions as follows: 1) response to nutrient levels, process utilizing autophagic mechanism, autophagy (biological process); 2) transport vesicle, aggresome, chaperone complex (cellular component); and 3) cell adhesion molecule binding, GTPase activity, ubiquitin-like protein ligase binding (molecular function) (Fig. 3A; Additional file 1: Table S3). In KEGG enrichment analysis, the differentially expressed autophagy-related genes are mainly involved in the process of autophagy and AMPK signaling pathway (Fig. 3B and C; Additional file 1: Table S4).
Autophagy levels of PBMC in WalkByLab registry PAD patients
So far, our analysis revealed that autophagy-related genes’ mRNAs are expressed at low levels in PBMCs of PAD patients and their potential biological functions are mainly associated with autophagy according to GO and KEGG analysis. Next, we assessed the expression of important autophagy markers (p62, beclin1, and LC3B) to explore the autophagy level between PAD and health.
First, the mRNA levels of p62, beclin1, and LC3B were examined in GSE57691 No significant differences in the mRNA levels of p62 and beclin1 were detected between PAD patients and healthy controls. However, the gene expression level of LC3B was found to be significantly decreased in PAD patients (Fig. 4A). Next, we validated a possible deregulation of LC3B mRNA expression in PBMCs from WalkByLab patients using qRT-PCR. LC3B in WalkByLab participates was less expressed in the PBMCs of patients with PAD (Fig. 4B). Finally, we assessed the protein levels of p62, beclin1 and LC3B in the PBMCs of WalkByLab participates. It turns out that there was a significant decrease in beclin1 and LC3B-II protein levels in PAD patients as compared to healthy controls (Fig. 4C). There was a trend of increased p62 level, but it was not statistically significant. In summary, results showed that protein levels of autophagy markers were decreased in PBMCs from patients with PAD.
Quantitative analysis of immune microenvironment
PAD is a systemic disease of the arteries, which manifests in the periphery (e.g. legs, feet) but usually affects the majority of all body arteries. PAD is also an inflammatory disease, and the degree of systemic inflammation is even greater than in CAD due to the large number of arteries affected [17]. Thus, we performed a ssGSEA here to quantify the immune microenvironment of the artery wall and compared tissue from PAD patients with healthy controls. ssGSEA revealed 16 types of immune cells and 13 types of immune-related pathways in each sample, which is shown in the heat map (Fig. 5A). A higher level of immune cell infiltration (including aDCs, iDCs, Mast cells, neutrophils, NK cells, pDCs, T helper cells, Tfh, Th1, and Th2 cells) was identified in PAD as compared with healthy controls (Fig. 5B). Moreover, more active immune-related pathways (including APC co-stimulation, cytokine-cytokine receptor (CCR) interaction, check point, MHC class I, and T cell co-stimulation activity) were observed in PAD as compared with healthy controls (Fig. 5C). Further Spearman’s correlation analysis showed that the identified subset of infiltrating immune cell subgroups and active immune-related pathways were moderately to strongly correlated, respectively (Fig. 5D and E).
Next, we analyzed whether these 20 differentially expressed autophagy-related genes were correlated with immune function in PAD. We found that most autophagy-related genes (except for KLHL24, HSPB8, CX3CL1 and CCL2) correlated with immune cells infiltration and immune-related pathways (Fig. 5F). Among them, RB1CC1 and FOXO1 genes showed the highest reported number of relationships with immune microenvironment. Most infiltrating immune cells linked to autophagy-related genes were PBMCs such as iDCs (involved genes n = 6), pDCs (n = 8) and Th2 cells (n = 9). The top immune-related pathways associated with autophagy-related genes were found to be CCR interactions (n = 10), as demonstrated in Fig. 5F. Specifically, the following genes were identified as being involved in these pathways: VAMP3, ST13, RPS6KB1, RB1CC1, PEA15, MAP1LC3B, HSPA5, FOXO1, EEF2K, and BNIP3.
Expression level of chemokines in plasma
Since CCR interaction was the top active immune-related pathway that linked to autophagy-related genes in the ssGSEA analysis, and combined with the results that autophagy-related genes and autophagy level were decreased in PAD patients’ PBMC, we suspect that CCR interaction might play an important role in the development and progression of PAD. Among CCR interactions, chemokines are best known for their ability to stimulate the migration of leukocytes. Since leukocytes are central to the progression of PAD and to our analysis of autophagy in PBMCs, chemokines play a particularly central role in our investigation.
Hence, we used the human chemokine antibody array to explore the plasma levels of chemokines (3 PAD and 3 health controls) and found that growth-related oncogene (GRO) and neutrophil activating protein2 (NAP2) chemokines had highest fold change expression between the two groups (Fig. 6A and B). The mRNA levels of GRO family (GROα, GROβ and GROγ) and NAP2 were also highly expressed in PAD patients in GSE57691 (Fig. 6C). Thus, we next utilized ELISA kits to measure the plasma expression levels of GRO and NAP2 in WalkByLab registry cohort of 179 participants. GRO and NAP2 levels were significantly higher in PAD patients’ plasma, and NAP2 (AUC: 0.752) demonstrated a better diagnostic performance for PAD patients as compared with GRO (AUC: 0.670) (Fig. 6D). Similar results were also shown in the propensity score matching analysis cohort after excluding the potential influence of age (Fig. 6E).
The relationship between autophagy-related genes, chemokines and endothelial function
Spearman correlation analysis, which investigate the potential pairwise correlation between these three factors, showed that autophagy-related genes KLHL24, VAMP3, HSPA5, and ST13 were moderately negative correlated with chemokine GRO (|Cor|> 0.4) (Fig. 6F), and showed mild to moderate negative correlation with chemokine NAP2 (|Cor|> 0.3) (Fig. 6G). Fourteen autophagy-related genes were moderately positive correlated with FMD (Cor > 0.4) (Fig. 6H). However, the chemokines GRO and NAP2 were not correlated with FMD (P > 0.05) (Fig. 6I).
The correlation between plasma chemokines and walking capacity
Gardner treadmill testing is a gold standard exercise stress test, providing clinically important diagnostic and prognostic information for patients with symptomatic PAD [18]. The results of the Spearman correlation analysis revealed that there were significant correlations between the levels of GRO and NAP2 and the walking capacity. Specifically, the GRO level showed a mild negative correlation with pain-free walking distance (Correlation coefficient, Cor = -0.331) and maximum walking distance (Cor = -0.340). The NAP2 level demonstrated a mild to moderate negative correlation with pain-free walking distance (Cor = -0.387) and maximum walking distance (Cor = -0.403). Besides, there was a moderate positive correlation between the levels of GRO and NAP2 (Cor = 0.434), and a strong positive correlation between pain-free walking distance and maximum walking distance (Cor = 0.913). The correlation network plots for pairwise correlations among the studied variables are showed in Fig. 7A (all P < 0.05).
Furthermore, participates which showed during the walking test a walking time less than 6 min had higher levels of GRO and NAP2 protein (Fig. 7B). ROC analysis to determine the diagnostic performance in walking time less than 6 min showed that NAP2 (AUC: 0.743) had a higher efficacy than GRO (AUC:0.701) (Fig. 7C).
The correlation between autophagy-related genes and walking capacity
Spearman correlation analysis showed that, except for CCL2 and CSTB genes, all other autophagy-related genes exhibited a significant positive correlation with walking distance, with correlation coefficients ranging from moderate to strong (0.40 < Cor < 0.73) (Fig. 7D). Specifically, MAP1LC3B gene demonstrated strong positive correlations with pain-free (Cor = 0.70) and maximum walking distance (Cor = 0.73).
The correlation between plasma chemokines and major adverse cardiovascular events (MACE)
We reviewed the history of MACE occurrence, including stroke, transient ischemic attack, and myocardial infarction, among the enrolled participants. Our results showed that there was no statistical difference in NAP2, GRO, and ankle-brachial index (ABI) levels between the groups with and without the history of MACE occurrence (P > 0.5) (Additional file 2: Figure S1). The AUC of GRO was 0.529 and that of NAP2 was 0.516, which is marginally higher than that of ABI (AUC: 0.508). Overall, these biomarkers almost had little predictive power for MACE.
Construction of a nomogram to predict impaired exercise tolerance
Nomogram is a powerful tool used for predicting the incidence of diseases or patients’ prognosis. Here, we constructed a nomogram model to predict impaired walking capacity. A walking time < 6 min in Gardner treadmill testing was defined as the endpoint. LASSO regression was performed firstly to select the most effective variables from all clinical features and ten clinical parameters were obtained (Fig. 8A). The respective AUC values of these parameters were shown in Fig. 8B.
Given the significant correlation between the levels of GRO and NAP2 and walking capacity, we selected these two factors as candidate variables for multivariate logistic regression analysis and nomogram construction (Fig. 8C and Additional file 2: Figure S2). Due to the limited sample size available for the autophagy-related genes in our study, we opted not to include them as candidate variables for analysis. ROC curves analysis showed that the nomogram integrated with NAP2 (AUC: 0.860) had a better predicted efficacy than nomogram with GRO (AUC: 0.854) (Fig. 8D). Finally, calibration plots were used to visualize the performances of the nomograms. The calibration plot confirmed the performance of our predictive model (Fig. 8E).
Discussion
Autophagy is a cellular self-recycling process, which is involved in both vascular homeostasis as well as physiological and pathological arterial remodeling. Physiological autophagy serves as a protective mechanism to maintain normal cardiovascular function, whereas disordered autophagy contributes to the development and progression of diseases [19]. It seems to play a protective role in the early phase of atherosclerosis, as deletion of autophagy genes (Atg5) accelerates the progression of atherosclerosis in murine models [20]. In clinical studies, protein levels of autophagy‑associated protein beclin1 and light chain 3 (LC3) were decreased in PBMCs from patients with acute myocardial infarction compared with healthy controls or patients suffering from stable angina pectoris [21]. Based on these findings, it is speculated that acute cardiovascular events and increased severity of disease might lead to a reduction in autophagy and vice versa [6, 21]. However, the role of autophagy in cardiovascular disease development is not entirely clear, in particular, there are little study related to PAD.
In this study, we identified 20 potential autophagy-related genes in PAD through bioinformatic analysis, and validated changes in the expression of these genes in peripheral blood samples from the WalkByLab registry trial. All identified genes were expressed at low levels or had the tendency of decreasing in PAD patients’ PBMCs, which indicated a role of autophagy-related genes in PAD. The GO and KEGG enrichment analysis further revealed that these autophagy-related genes might be involved in autophagy activity and control of autophagy activation. Indeed, western blot analysis demonstrated that autophagic marker proteins beclin1 and LC3B levels were significantly decreased in the PBMCs of PAD patients. Beclin1 and LC3B are common biomarkers of autophagy, which are widely used to evaluate the level of autophagy. Our study now provided evidence that PAD is linked to a reduced expression level of these autophagy-related genes and a decreased autophagy marker protein level.
Inflammation is one of the main features of PAD. Autophagy plays an important role in inflammation by regulating the development, homeostasis, and survival of inflammatory cells such as macrophages, neutrophils, and lymphocytes. Autophagy also influences the transcription, processing, and release of a variety of cytokines and is regulated by cytokines [22]. Autophagy can trigger a dynamic response to inflammation. Active autophagy suppresses excessive inflammatory response by inhibiting inflammasome activation, facilitating damage-associated molecular patterns and damaged mitochondria clearance, and destroying inflammatory mediators [23]. Suppressed autophagy resulted in an enhanced inflammatory response [24].
Thus, we explored the overall characteristics of the immune microenvironment of artery tissue from PAD patients and healthy populations and the association with autophagy-related genes. We found higher level of immune cell infiltration and more active immune-related pathways in PAD and that the CCR interaction was the most active immune-related pathway negatively associated to autophagy-related genes. This seemed to imply that the decreased autophagy will increase inflammation in PAD. Here, the cytokines include the chemokines, the PDGF family, and the TGF-β family [25]. Among CCR interactions, chemokines are best known for their ability to stimulate the migration of leukocytes (e.g. PBMCs) and chemokines in particular are known to play an important role in the pathophysiology of cardiovascular disease [26]. Chemokines bind to specific G protein-coupled receptors and exert distinct functions, such as regulating the activation of leukocytes and coordinating their trafficking to the sites of inflammation [27]. For example, chemokine RANTES level is a useful marker of CAD severity because elevated chemokine RANTES level in patients with stable angina may predict the high risk of plaque formation in early stages of atherosclerosis [28]. Therefore, we measured the expression profiles of human chemokines in plasma from PAD patients and healthy individuals to identify a potential role of chemokines in PAD pathophysiology. ELISA analysis of the WalkByLab participants showed that human chemokine GRO and NAP2 was highly expressed in PAD which was also demonstrated by a strict propensity score-matched analysis after adjustment for covariate of age. In addition, our analysis indicated that GRO and NAP2 may be promising diagnostic biomarkers for PAD, as reflected by their respective AUC values of 0.670 and 0.752. Moreover, our analysis showed that the levels of GRO and NAP2 were not associated with MACE (acute events), therefore indicating their specificity for a systemic PAD-related pathophysiology. Interestingly, our study also revealed a negative correlation between GRO and NAP2 levels and the expression of autophagy-related genes KLHL24, VAMP3, HSPA5, and ST13, suggesting a potential link between the dysregulation of autophagy and chemokine-mediated inflammation in PAD. Further investigation is needed to elucidate the mechanistic relationship between autophagy and GRO/NAP2 signaling in the context of PAD pathophysiology.
We next examined other crucial characteristics of PAD pathophysiology. Important cofactors and clinical parameters for the evaluation of PAD are endothelial function measured by FMD and walking distance measured by a Gardner treadmill test. Our study demonstrated that plasma GRO and NAP2 levels were not correlated with participants’ FMD, but were negatively correlated with the pain-free walking and maximum walking distance in Gardner treadmill testing. Besides, GRO and NAP2 were able to predict poor walking capacity (walking time < 6 min in Gardner treadmill testing) with respective AUC values of 0.701 and 0.743. Although FMD is a very valuable clinical parameter, it is also very sensitive to a variety of external influences, such as biological age, diet, or exercise level of the patient [29]. In clinical practice FMD is therefore not used for the primary diagnosis of PAD. However, the Gardner treadmill walk test is considered the gold standard to accurately determine the stage of PAD according to the Fontaine stage. To our knowledge, a role of GRO and NAP2 for PAD in the context of walking distance and walking capacity has never been demonstrated before. Previous studies have mainly focused on the role of GRO and NAP2 in cancer development and metastasis of various malignancies [30], whereas few reports have focused on cardiovascular disease, and even fewer in the context of PAD. In the study of Ku et al. [31], upregulated NAP2 was significantly associated with severity of coronary artery stenosis in patients with diabetes. Recently Wang et al. showed that elevated plasma NAP2 level was independently associated with critical limb ischemia, not with diabetes [32], however a connection with the walking capacity was not considered. Until now the deregulation of GRO was never shown before in the context of PAD. Chemokines GROα, GROβ, and GROγ constitute the three members of the GRO family [33], which together with NAP2, share the chemokine receptor CXCR2 as a common binding receptor [34]. GROγ was reported to be upregulated in acute coronary syndromes and related to the severity of the inflammatory response [35]. Moreover, a recent study demonstrated that CXCL3 (the coding gene of GROγ protein) gene was highly expressed in acute ischemic stroke and it can be recognized as a predictor for more infarct volume [36]. In summary, chemokines GRO and NAP2 appear to be useful biomarkers in PAD, particularly in predicting impaired walking ability.
Finally, we further analyzed and processed plasma GRO or NAP2 levels with the relevant clinical parameters obtained by LASSO analysis to build nomogram models for the prediction of poor walking capacity (walking time < 6 min in Gardner treadmill testing). The NAP2-derived nomogram showed a well predictive performance, with an AUC of 0.860, while the GRO-derived nomogram had a slightly lower AUC of 0.854. These evidence demonstrated that plasma NAP2 indeed could be used as a predictive biomarker and the derived nomogram model simply and accurately predicted impaired walking capacity in patients with PAD.
The here presented study has certain limitations that need to be acknowledged. Firstly, the sample size was relatively small, and therefore, it is important to exercise caution when generalizing the results. Secondly, this study was retrospective, and as such, there is a possibility of selection bias. Nevertheless, the retrospective nature of the study was necessary to pave the way for future prospective studies. Thirdly, a correlation between autophagy-related genes and chemokines with walking ability was discovered, suggesting that autophagy may play a role in the progression of PAD by regulating chemokines and inflammation. However, this potential mechanism requires further validation in animal models, and thus, future studies are necessary to shed more light on this issue.
Conclusions
In conclusion, our study shed light on the interrelationships of autophagy and autophagy-related genes in the development of PAD, whcih were previously unknown. Our comprehensive analysis identified a set of 20 autophagy-related genes that were expressed at lower levels in PAD and demonstrated that the autophagy levels of PBMCs were significantly decreased in patients with PAD. This finding has significant implications for future clinical strategies and for deepening the understanding of disease onset and progression, as well as providing novel biomarkers. These results will serve to guide the clinic in developing new prognostic strategies for the evaluation of PAD. Our bioinformatics analysis further revealed a link between vascular inflammation and the effect of cytokines, particularly chemokines, with autophagy-related genes. This was also further validated in subsequent experiments. We found that GRO and NAP2 levels were elevated in the plasma of PAD patients, with a mild-to-moderate negative correlation between them and the autophagy-related genes KLHL24, VAMP3, HSPA5, and ST13. The effects of GRO and NAP2 on walking capacity in PAD have never been proven before. Our study reveals a mild-to-moderate negative correlation between chemokines and pain-free walking distance and maximal walking distance in the Gardner treadmill test. This suggests that reduced autophagy may influence the onset and progression of PAD (especially, the walking capacity) by affecting chemokines (inflammation). Finally, we constructed nomogram models that can predict poor walking capacity. The results showed that NAP2 is a promising biomarker for PAD patients, and its derived nomogram had stronger predictive potential than the GRO derived nomogram. This nomogram provides a valuable tool for clinicians to improve their understanding of a patient's disease risk and prognosis, thereby facilitating informed clinical decision-making. In all, our study provides new insights into the molecular mechanisms underlying PAD and may have implications for the development of new diagnostic and therapeutic approaches for PAD.
Availability of data and materials
Publicly available datasets were analyzed in this study. This data can be found here: GEO database, accession number: GSE57691.
Abbreviations
- ABI:
-
Ankle brachial index
- AUC:
-
Area under the ROC curve
- BP:
-
Biological process
- CAD:
-
Coronary artery disease
- CC:
-
Cellular component
- CCR:
-
Cytokine and cytokine receptors
- ELISA:
-
Enzyme linked immunosorbent assay
- FC:
-
Fold change
- FMD:
-
Flow-mediated dilation
- GEO:
-
Gene expression omnibus
- GO:
-
Gene Ontology
- GRO:
-
Growth-related oncogene
- KEGG:
-
Kyoto encyclopedia of genes and genomes
- LASSO:
-
Least absolute shrinkage and selection operator
- LC3:
-
Light chain
- MACE:
-
Major adverse cardiovascular events
- MF:
-
Molecular function
- NAP2:
-
Neutrophil activating protein2
- PAD:
-
Peripheral arterial disease
- PBMCs:
-
Peripheral blood mononuclear cells
- PCA:
-
Principal component analysis
- PPI:
-
Protein–protein interaction
- qRT-PCR:
-
Quantitative real-time polymerase chain reaction
- ROC:
-
Receiver operating characteristic
- ssGSEA:
-
Single sample gene set enrichment analysis
- TBST:
-
Tris‑buffered saline containing 0.1% Tween 20
References
Cooke JP, Meng S. Vascular regeneration in peripheral artery disease. Arterioscler Thromb Vasc Biol. 2020;40(7):1627–34.
Fowkes FG, Aboyans V, Fowkes FJ, McDermott MM, Sampson UK, Criqui MH. Peripheral artery disease: epidemiology and global perspectives. Nat Rev Cardiol. 2017;14(3):156–70.
Fowkes FG, Rudan D, Rudan I, Aboyans V, Denenberg JO, McDermott MM, Norman PE, Sampson UK, Williams LJ, Mensah GA, et al. Comparison of global estimates of prevalence and risk factors for peripheral artery disease in 2000 and 2010: a systematic review and analysis. Lancet (London, England). 2013;382(9901):1329–40.
Ryter SW, Lee SJ, Smith A, Choi AM. Autophagy in vascular disease. Proc Am Thorac Soc. 2010;7(1):40–7.
De Meyer GR, Grootaert MO, Michiels CF, Kurdi A, Schrijvers DM, Martinet W. Autophagy in vascular disease. Circ Res. 2015;116(3):468–79.
Hughes WE, Beyer AM, Gutterman DD. Vascular autophagy in health and disease. Basic Res Cardiol. 2020;115(4):41.
Yan Y, Finkel T. Autophagy as a regulator of cardiovascular redox homeostasis. Free Radical Biol Med. 2017;109:108–13.
Abdellatif M, Sedej S, Carmona-Gutierrez D, Madeo F, Kroemer G. Autophagy in cardiovascular aging. Circ Res. 2018;123(7):803–24.
Jeong IH, Bae WY, Choi JS, Jeong JW. Ischemia induces autophagy of endothelial cells and stimulates angiogenic effects in a hindlimb ischemia mouse model. Cell Death Dis. 2020;11(8):624.
Biros E, Gäbel G, Moran CS, Schreurs C, Lindeman JH, Walker PJ, Nataatmadja M, West M, Holdt LM, Hinterseher I, et al. Differential gene expression in human abdominal aortic aneurysm and aortic occlusive disease. Oncotarget. 2015;6(15):12984–96.
Gäbel G LJ, Walker PJ, Golledge J, Biros E: Differential gene expression in human abdominal aortic aneurysm and atherosclerosis. GEO https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE57691, 2015.
Chu LH, Annex BH, Popel AS. Computational drug repositioning for peripheral arterial disease: prediction of anti-inflammatory and pro-angiogenic therapeutics. Front Pharmacol. 2015;6:179.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1–2):48–61.
Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, Li YH. A novel ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma. Int J Biol Sci. 2020;16(13):2430–41.
Li K, Zemmrich C, Bramlage P, Persson AB, Sacirovic M, Ritter O, et al. Effect of ACEI and ARB treatment on nitric oxide-dependent endothelial function. VASA Zeitschrift fur Gefasskrankheiten. 2021;50(6):413-22.
Mizobuchi K, Jujo K, Minami Y, Ishida I, Nakao M, Hagiwara N. The Baseline Nutritional Status Predicts Long-Term Mortality in Patients Undergoing Endovascular Therapy. Nutrients. 2019;11(8):1745.
Gardner AW, Montgomery PS, Wang M, Xu C. Predictors of health-related quality of life in patients with symptomatic peripheral artery disease. J Vasc Surg. 2018;68(4):1126–34.
Mei Y, Thompson MD, Cohen RA, Tong X. Autophagy and oxidative stress in cardiovascular diseases. Biochem Biophys Acta. 2015;1852(2):243–51.
Razani B, Feng C, Coleman T, Emanuel R, Wen H, Hwang S, Ting JP, Virgin HW, Kastan MB, Semenkovich CF. Autophagy links inflammasomes to atherosclerotic progression. Cell Metab. 2012;15(4):534–44.
Zhao K, Xu XS, Meng X, Li YL, Li JF, Chen WQ. Autophagy of monocytes attenuates the vulnerability of coronary atherosclerotic plaques. Coron Artery Dis. 2013;24(8):651–6.
Qian M, Fang X, Wang X. Autophagy and inflammation. Clin Transl Med. 2017;6(1):24.
Gong L, Pan Q, Yang N. Autophagy and inflammation regulation in acute kidney injury. Front Physiol. 2020;11: 576463.
Nakajima T, Wang B, Ono T, Uehara Y, Nakamura S, Ichinohe K, Braga-Tanaka I 3rd, Tanaka S, Tanaka K, Nenoi M. Differences in sustained alterations in protein expression between livers of mice exposed to high-dose-rate and low-dose-rate radiation. J Radiat Res. 2017;58(4):421–9.
Zhou Z, Zhang J, Li X, Xia C, Han Y, Chen H. Protein microarray analysis identifies key cytokines associated with malignant middle cerebral artery infarction. Brain Behav. 2017;7(8):e00746–e00746.
Koper-Lenkiewicz OM, Kamińska J, Lisowska A, Milewska A, Hirnle T, Dymicka-Piekarska V. Factors Associated with RANTES Concentration in Cardiovascular Disease Patients. Biomed Res Int. 2019;2019:3026453–3026453.
Vajen T, Koenen RR, Werner I, Staudt M, Projahn D, Curaj A, Sönmez TT, Simsekyilmaz S, Schumacher D, Möllmann J, et al. Blocking CCL5-CXCL4 heteromerization preserves heart function after myocardial infarction by attenuating leukocyte recruitment and NETosis. Sci Rep. 2018;8(1):10647–10647.
Podolec J, Kopec G, Niewiara L, Komar M, Guzik B, Bartus K, Tomkiewicz-Pajak L, Guzik TJ, Plazak W, Zmudka K. Chemokine RANTES is increased at early stages of coronary artery disease. J Physiol Pharmacol. 2016;67(2):321–8.
Mućka S, Miodońska M, Jakubiak GK, Starzak M, Cieślar G, Stanek A. Endothelial Function Assessment by Flow-Mediated Dilation Method: A Valuable Tool in the Evaluation of the Cardiovascular System. Int J Environ Res Public Health. 2022;19(18):11242.
Wang H, Wang T, Dai L, Cao W, Ye L, Gao L, Zhou B, Zhou R. Effects of CXCL3 on migration, invasion, proliferation and tube formation of trophoblast cells. Placenta. 2018;66:47–56.
Ku EJ, Cho KC, Lim C, Kang JW, Oh JW, Choi YR, Park JM, Han NY, Oh JJ, Oh TJ, et al. Discovery of plasma biomarkers for predicting the severity of coronary artery atherosclerosis by quantitative proteomics. BMJ Open Diabetes Res Care. 2020;8(1):e001152.
Wang X, Li J, Gan L, Liu Q. Plasma NAP-2 levels are associated with critical limb ischemia in peripheral arterial disease patients. Exp Biol Med (Maywood). 2019;244(1):22–7.
Kent LM, Smyth LJC, Plumb J, Clayton CL, Fox SM, Ray DW, Farrow SN, Singh D. Inhibition of lipopolysaccharide-stimulated chronic obstructive pulmonary disease macrophage inflammatory gene expression by dexamethasone and the p38 mitogen-activated protein kinase inhibitor N-cyano-N’-(2-{[8-(2,6-difluorophenyl)-4-(4-fluoro-2-methylphenyl)-7-oxo-7,8-dihydropyrido[2,3-d] pyrimidin-2-yl]amino}ethyl)guanidine (SB706504). J Pharmacol Exp Ther. 2009;328(2):458–68.
Tecimer T, Dlott J, Chuntharapai A, Martin AW, Peiper SC. Expression of the chemokine receptor CXCR2 in normal and neoplastic neuroendocrine cells. Arch Pathol Lab Med. 2000;124(4):520–5.
Zhang L, Li J, Liang A, Liu Y, Deng B, Wang H. Immune-related chemotactic factors were found in acute coronary syndromes by bioinformatics. Mol Biol Rep. 2014;41(7):4389–95.
Martha SR, Cheng Q, Fraser JF, Gong L, Collier LA, Davis SM, Lukins D, Alhajeri A, Grupke S, Pennypacker KR. Expression of cytokines and chemokines as predictors of stroke outcomes in acute ischemic stroke. Front Neurol. 2020;10:1391–1391.
Funding
Open Access funding enabled and organized by Projekt DEAL. The authors gratefully acknowledge financial support from the State of Brandenburg (MWFK), the Federal Ministry of Education and Research (Bundesministerium für Bildung und Forschung), and the Chinese Scholarship Council (CSC).
Author information
Authors and Affiliations
Contributions
All authors made a significant contribution to the work reported, that is in the conception and study design (AP, IB, and PH); execution (MD, KL, MS, and CZ); acquisition of data, analysis, and interpretation (EB, OR, PB, and PH); took part in drafting (MD, KL and MS), revising or critically reviewing the article (AP, IB, and PH). All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The studies involving human participants were reviewed and approved by the ethical committee of the Cottbus Medical Association (Landesärztekammer Cottbus, study number of the ethics committee: AS 74(bB)/2018). The screening trial is performed in accordance with the principles of the declaration of Helsinki. The patients/participants provided their written informed consent to participate in this study.
Consent for publication
All authors gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Table S1. Primer sequences. Table S2. Clinical characteristics of PAD patients and controls in the qRT-PCR analysis. Table S3. GO enrichment analysis of 20 differentially expressed autophagy-related genes. Table S4. KEGG pathway enrichment analysis of 20 differentially expressed autophagy-related genes.
Additional file 2:
Figure S1. The correlation between plasma GRO/NAP2 chemokines and major adverse cardiovascular events. Figure S2. Nomogram integrated with GRO predicting walking time<6min in Gardner treadmill testing.
Additional file 3:
Figure S1. The original blot image of P62. Figure S2. The original blot image of Beclin1. Figure S3. The original blot image of LC3B. Figure S4. The original blot image of GAPDH.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Dai, M., Li, K., Sacirovic, M. et al. Autophagy-related genes analysis reveals potential biomarkers for prediction of the impaired walking capacity of peripheral arterial disease. BMC Med 21, 186 (2023). https://doi.org/10.1186/s12916-023-02889-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12916-023-02889-5