- 1Department of Medical and Surgical Sciences (DIMEC), Alma Mater Studiorum-University of Bologna, Bologna, Italy
- 2Interdepartmental Centre - Alma Mater Research Institute on Global Challenges and Climate Change, University of Bologna, Bologna, Italy
- 3Vascular Surgery Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
- 4Department of Clinical and Molecular Sciences, Università Politecnica delle Marche, Ancona, Italy
- 5Clinic of Laboratory and Precision Medicine, IRCCS INRCA, Ancona, Italy
- 6Nephrology, Dialysis and Renal Transplant Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
- 7IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
- 8Pathology Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy
Background: Carotid endarterectomy (CEA) for the prevention of upcoming vascular and cerebral events is necessary in patients with high-grade stenosis (≥70%). In the framework of the Italian National project Age.It, a pilot study was proposed aiming at the discovery of a molecular signature with predictive potential of carotid stenosis comparing 65+ asymptomatic and symptomatic inpatients.
Methods: A total of 42 inpatients have been enrolled, including 26 men and 16 women, with a mean age of 74 ± 6 years. Sixteen symptomatic and 26 asymptomatic inpatients with ≥70% carotid stenosis underwent CEA, according to the recommendations of the European Society for Vascular Surgery and the Society for Vascular Surgeons. Plaque biopsies and peripheral blood samples from the same individuals were obtained. Hematobiochemical analyses were conducted on all inpatients, and plasma cytokines/molecules, such as microRNAs (miRs), IL-6, sIL-6Ralpha, sgp130, myostatin (GDF8), follistatin, activin A, CXCL9, FGF21, and fibronectin, were measured using the ELISA standard technique. MiR profiles were obtained in the discovery phase including four symptomatic and four asymptomatic inpatients (both plasma and plaque samples), testing 734 miRs. MiRs emerging from the profiling comparison were validated through RT-qPCR analysis in the total cohort.
Results and conclusion: The two groups of inpatients differ in the expression levels of blood c-miRs-126–5p and -1271–5p (but not in their plaques), which are more expressed in symptomatic subjects. Three cytokines were significant between the two groups: IL-6, GDF8, and CXCL9. Using receiver operating characteristic (ROC) analysis with a machine learning-based approach, the most significant blood molecular signature encompasses albumin, C-reactive protein (CRP), the percentage of monocytes, and CXCL9, allowing for the distinction of the two groups (AUC = 0.83, 95% c.i. [0.85, 0.81], p = 0.0028). The potential of the molecular signature will be tested in a second cohort of monitored patients, allowing the application of a predictive model and the final evaluation of cost/benefit for an assessable screening test.
1 Introduction
Epidemiological studies reveal that 21% of individuals aged 30–79 years have carotid plaque, defined as a focal carotid intima-media thickness of 0.5–1.5 mm or more invading into the lumen (≤50 of stenosis), and 1.5% of individuals have carotid stenosis (≥50%), equivalent to approximately 816 million people with carotid plaque and 58 million with carotid stenosis (1). As is well recognized, people with carotid plaque or carotid stenosis are at an increased risk of developing systemic cardiovascular diseases (CVDs) (2–4), and carotid endarterectomy (CEA) is deemed necessary in patients with high-grade stenosis (≥70%) to prevent forthcoming vascular and cerebral events.
In this perspective, the possibility to monitor the development and the stability of carotid plaque is clinically relevant especially for asymptomatic and elderly patients. On the other hand, symptomatic patients are those subjects who have manifested signs of prior cerebrovascular events, such as transient ischemic attack, ischemic stroke, and transient monocular vision loss, and urgently require surgical intervention of carotid artery (stenosis ≥70%). Concurrently, the identification of asymptomatic patients with carotid plaque has been recommended in assessing cardiovascular risk (5), and various circulating blood biomarkers may aid in the faster identification of carotid atherosclerosis. This research area has been growing during last decade, focusing on immune system cell sub-populations, pro-inflammatory molecules, and tissue/blood circulating microRNAs (c-miRs) (6).
Our team has recently investigated the differences in miR profiling expression within the human artery wall, revealing miR-related differences in the carotid, femoral, and abdominal aortic arteries in both health and pathological conditions (atheroma) (7, 8), and other researchers have deeply investigated plaque stability/instability and associated miR expression (9–11). In the light of these previous results, the primary object of the current work has been to identify a miR-based molecular signature to distinguish asymptomatic and symptomatic older (≥65 years) inpatients in blood and plaques, plausibly.
The use of c-miRs and cytokines/chemokines as potential biomarkers has been well recognized for both the aging process and different pathologies/age-related diseases (6, 12–17). C-miRs are protected from plasma/serum RNase activity through their transport on RNA–protein complexes with lipoproteins or Ago2 proteins (18, 19), as well as through their inclusion inside extracellular nano-microvesicles. In the current work, c-miRs have been obtained by plasma RNA extraction including all the components and assessed by standard RT-qPCR. Furthermore, the selection of specific cytokines/chemokines, in accordance with the ongoing Italian National project (PNRR-PE8: Age.It at the website: https://ageit.eu/wp/), has been proposed for the implementation of the molecular signature. Cytokines/chemokines, including IL-6, CCL2, CXCL9, and CXCL12, are recognized to be important drivers of atheroma development and the inflammatory microenvironment (20–24).
The first hypothesis testing aimed to identify a common miR-based signature between blood c-miRs and carotid atheroma, likely able to distinguish between asymptomatic and symptomatic patients. The second hypothesis testing focused on specific cytokines/molecules, such as IL-6, sIL-6Ralpha, sgp130, GDF8, follistatin, activin A, CXCL9, FGF21, and fibronectin. The objective was to enhance the potential miR-based signature, with the aim of better distinguishing between asymptomatic and symptomatic patients.
As part of the Italian National project (Age.It), the envisioned final achievement will be the validation of the identified molecular signature in patients with different carotid stenosis degree during the monitoring period and the possibility of applying a predictive model.
2 Materials and methods
2.1 Pilot study and patients’ recruitment
The current pilot study focused on 42 individuals, comprising 26 asymptomatic (11 women and 15 men) and 16 symptomatic (5 women and 11 men) inpatients diagnosed with carotid artery disease. The two groups had a mean age of 72 ± 5 years and 77 ± 7 years, respectively. Recruitment was carried out after obtaining approval from the regional ethical committee (88/2019/Sper/AOUBo) and following an amendment (n. EM511–2023) for the project inclusion in Age.it (PNRR-PE8).
All recruited inpatients showed carotid artery stenosis ≥70% and underwent CEA according to the recommendations of the European Society for Vascular Surgery and the Society for Vascular Surgeons, since 70% of stenosis is a critical threshold for urgent surgical intervention. Symptomatic carotid stenosis was defined as the occurrence of ipsilateral cerebral ischemic events (major or minor stroke, TIA, or amaurosis fugax) within the last 6 months. CEA-derived biopsies (plaques) and blood samples (before CEA) were collected from each inpatient. Neurological symptoms (amaurosis fugax, TIA, and minor and major stroke), vascular risk factors, morbidities or comorbidities (hypertension, coronary artery disease, chronic obstructive pulmonary disease, dyslipidemia, type II diabetes, current smoking, and chronic renal failure), and current therapies were recorded. Exclusion criteria were current acute illnesses; hepatic, severe renal, or cardiac insufficiency; obesity; and type II diabetes with severe complications.
All inpatients were also tested for hemato-biochemical parameters, i.e., white blood cell count (WBC), erythrocytes (RBC), hemoglobin (HGB), hematocrit (HCT), mean corpuscular volume (MCV), the average hemoglobin content in erythrocytes (MCH), mean hemoglobin concentration in erythrocytes (MCHC), red blood cell distribution (% RDW-CV, RDW-SD), neutrophils, lymphocytes, monocytes, eosinophils, basophils, platelets, mean platelet volume (MPV), urea, creatinine, estimated glomerular filtration rate (eGFR), uric acid, sodium, potassium, total calcium, total bilirubin, direct bilirubin, indirect bilirubin, aspartate aminotransferase (AST), alanine aminotransferase (ALT), γ-glutamyl transferase (GGT), total amylase, creatine kinase (CK), lactate dehydrogenase (LDH), C-reactive protein, albumin, total proteins, glucose, total cholesterol, high-density lipoprotein (HDL) cholesterol, LDL (low-density lipoprotein) cholesterol, and triglycerides (TGC).
2.2 Sample processing
Whole blood was collected from 42 inpatients before surgery in vacuum tubes following standard procedures and processed within 1 h. The tubes were centrifuged at 2,500 × g for 20 min at 4°C for plasma separation and aliquoted in cryotubes for long-term storage at −80°C. Carotid plaques were completely removed during surgery to preserve the plaque structure. After 24 h of decalcification, samples were cut into serial sections and the area with the highest percentage of stenosis was identified and defined for the analysis. Biopsies from 41 inpatients were formalin-fixed paraffin-embedded (FFPE) and routinely processed. Carotid atherosclerotic lesions were defined according to the American Heart Association classification and grouped as hemorrhagic or non-hemorrhagic plaque.
2.3 RNA extraction from plasma
Total RNA was extracted from plasma samples using the Total RNA Purification Kit (Norgen Biotek Corporation, Thorold, Ontario, Canada) and strictly adhering to the manufacturer’s instructions. Additionally, cel-miR-39 (Norgen) was spiked into each sample as an internal control for RNA extraction in plasma samples. On average, 100 μL of plasma produced 18 ± 5 ng/μL total RNA. The extracted RNA was quantified using NanodropTM One/OneC Microvolume UV-Vis 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).
2.4 RNA extraction from formalin-fixed, paraffin-embedded tissue
Total RNA was extracted from formalin-fixed paraffin-embedded (FFPE) sections. Each slice of tissue (20 µm thick) produced 19 ± 7 ng/µL total RNA on average. Briefly, four slices from each biopsy block were deparaffinized (∼20 min) and digested with protease (overnight at 50°C and 15 min at 80°C). RNA extraction was obtained using a commercial kit (RecoverAll™ Total Nucleic Acid Isolation Kit for FFPE, Thermo Fisher Scientific, Waltham, MA, USA), which allows isolation of RNA including miRs, following the manufacturer’s instructions. Subsequently, the extracted RNA was quantified using a NanodropTM One/OneC Microvolume UV-Vis spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).
2.5 Discovery phase with miR profiling
MiR profiles were obtained from eight selected inpatients, i.e., four symptomatic (two M, two F; age 75 ± 9 years; max stenosis of 80%) and four asymptomatic (two M, two F; age 76 ± 4 years; max stenosis of 95%), using human miR microfluidic cards for the assessment of 754 miRs (TaqMan Array Human MicroRNA A+B Cards/arrays; Applied Biosystems by Thermo Fisher Scientific, Waltham, MA, USA).
The inpatients were chosen for advanced age and maximum degree of occlusion in both carotid arteries. Carotid plaque specimens and plasma samples from the same patient were analyzed.
Six microliters of total RNA at the concentration of 4.2 ng/µL were used for the card/array analysis of the biopsies, and similarly, 6 µL of total RNA obtained from plasma samples was used for the card/array analysis. RNA was converted to cDNA by priming with a mixture of looped primers and then pre-amplified using the MegaPlex™ primer pools (Applied Biosystems by Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions, as previously published (25). The assay was performed using Applied Biosystems 7900 HT real-time PCR instrument. All the arrays/cards were normalized using the median of the overall miR expression on each array, thus obtaining ΔCt values. Only miRs expressed in all samples were selected for analyses and Ct values ≤ 30 were set as cutoff. To compare symptomatic and asymptomatic groups, the fold change (FC) were calculated. It was calculated on the estimated mean difference of ΔCt values between the two groups as 2−ΔΔCt, as applied in previous works (25). FC ≥ 2 and ≤ −2 were selected, thus obtaining miR profiling for further validation.
2.6 Validation phase with RT-qPCR in both plasma and plaques
RT-qPCR assays were performed on selected miRs that showed variation in FC values (FC ≥ 2 or ≤ −2) comparing symptomatic and asymptomatic inpatients in both plaque and plasma samples, considering a plausible cross-talk between blood circulating molecules and carotid artery wall tissue.
The analysis was extended to all subjects using TaqMan technologies (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s protocol. Five microliters of extracted RNA (both from plaque and plasma) were transcribed to cDNA using the TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, by Thermo Fisher Scientific, Waltham, MA, USA), followed by RT-qPCR with TaqMan MicroRNA Assays (Thermo Fisher Scientific cod. 4427975). The following assays were used: hsa-miR-134–5p (ID cod. 000459), hsa-miR-145–5p (ID cod. 002278), hsa-miR-151a-5p (ID cod. 002642), hsa-miR-34b-3p (ID cod. 002102), hsa-miR-451 (ID cod. 001105), hsa-miR-720 (ID cod. 002895), hsa-miR-126–5p (ID cod. 000451), and hsa-miR-1271–5p (ID cod. 002779). Data normalization was performed using 20 fmol of cel-mir-39 (Qiagen, Hilden, Germany) for plasma samples while hsa-miR-16 (ID cod. 000391) was adopted for biopsies among other housekeepers tested. The relative expression of all miRs was determined using the standard formula: 2−ΔCt.
2.7 Plasma cytokine/molecule assessment
All the cytokines/molecules, i.e., IL-6, sIL-6Ralpha, sgp130, GDF8, follistatin, activin A, CXCL9, FGF21, and fibronectin, were assessed with the Human Quantikine ELISA Kit in plasma samples using commercial kits (R&D, Minneapolis, USA) with catalog numbers D6050, DR600, DGP00, DGDF80, DFN00, DAC00B, DCX900, DF2100, and DFBN10, respectively. All samples were analyzed in duplicate according to the manufacturer’s instructions.
2.8 Statistical analysis
Statistical analysis was performed using both the software SPSS v.26 and Python v 3.11.7/scikit-learn v 1.3.2. Statistical significance was defined as p-value ≤ 0.05; Student’s t-test for unpaired and parametric data was applied to compare variables between symptomatic and asymptomatic groups.
The model for patient classification and molecular signature identification was developed from the full dataset (hemato-biochemical analysis, cytokines, and miRs). To ensure the model robustness, all analyses were performed using a leave-one-out cross-validation approach, where the model was fitted repeatedly on a subset of all the patients except for one, and the prediction was assessed on the left out one, repeating this procedure to obtain a prediction for every subject. Missing data were estimated using an Iterative Imputation algorithm (26) only on the training subset, then each feature was standardized to mean 0 and standard deviation 1. The most important features were selected using a Sequential Feature Selection method (27), based on the usefulness of the features for a Random Forest Classifier, starting with no features and adding them progressively based on the obtained improvement in a secondary leave-one-out cross-validation. The actual prediction was performed training a Random Forest Classifier with 100 decision trees. This prediction method allows obtaining both a binary prediction (symptomatic or asymptomatic) and a continuous one, with a degree of similarity to each class (going from 0, totally different, to 1, completely similar). The quality of the prediction for each set of features was estimated using the Area Under the Receiver Operating Curve (AUC-ROC) with 10 replications for each model, thus obtaining an average AUC value ± SD with the related confidence interval for each model, to account for the randomized initialization of the components of the whole prediction pipeline (imputation, sequential feature selection, and Random Forest Classifier).
3 Results
3.1 Pilot study
The two cohorts of the 26 asymptomatic and 16 symptomatic inpatients were assessed for hemato-biochemical parameters (Table 1). By comparing the two groups, the percentage of monocytes and albumin were significantly different (p = 0.036 and p = 0.005, respectively).
Table 1 Asymptomatic and symptomatic recruited inpatients and hemato-biochemical analysis (Student’s t-test).
3.2 Carotid plaques between asymptomatic and symptomatic
Histological analysis showed the presence of an atheromatous plaque in 41 obtained biopsies with minimal differences in the histological features of the atherosclerotic specimens. For this reason, attention was focused on the presence or absence of intraplaque hemorrhage as evidence of plaque vulnerability and susceptibility to rupture, potentially associated with the clinical setting of symptomatic/asymptomatic nature.
Carotid plaque specimens were therefore grouped as hemorrhagic plaque (N = 20) or non-hemorrhagic plaque (N = 21). The hemorrhagic plaques included 10 specimens from asymptomatic inpatients and 10 specimens from symptomatic inpatients. Conversely, the group of non-hemorrhagic plaques consisted of 15 specimens from asymptomatic and 6 specimens from symptomatic inpatients. Comparing the hemorrhagic with the non-hemorrhagic group, no significant differences in the selected, plaque-derived miRs and c-miRs expression levels were observed (data not shown).
3.3 MiRs profiling in carotid plaque and plasma from the discovery phase
MiRs profiling was performed to identify miR changes in plasma and carotid plaque from the same individual, comparing extreme phenotypes of symptomatic and asymptomatic inpatients. The average of detected miRs with card A was 205 in plasma samples and 279 in biopsies (miRs with Ct < 32 were 159 and 250, respectively), while the average detected with card B was 124 in plasma samples and 188 in biopsies (miRs with Ct < 32 were 71 and 147, respectively). Putative differentially expressed miRs were identified according to FC values (FC ≥ 2 and FC ≤ −2). A total of 45 and 63 miRs were found to satisfy such prerequisites in biopsy and plasma, respectively. All the data related to profiling and comparison of asymptomatic versus symptomatic inpatients are included as Supplementary Material (Supplementary Tables S1, S2) for both atheroma biopsies and plasma. When putative differently expressed miRs were compared between carotid plaque and plasma, eight shared miRs were identified: miR-126–5p, miR-134–5p, miR-145–5p, miR-151a-5p, miR-34b, miR-451a, miR-720, and miR-1271–5p (Figure 1). The figure also shows the FC between symptomatic and asymptomatic groups, while the relative expression average ( ± SD) is reported in Supplementary Table S3.
Figure 1 Venn diagram of common miRs between carotid plaque and plasma in the selected asymptomatic and symptomatic inpatients. Results of miR profiling in terms of shared miRs and FC (comparing asymptomatic and symptomatic inpatients, FC ≥ 2 and FC ≤ −2) in both carotid plaque and plasma. Data were validated in single RT-qPCR; only miR-126–5p and miR-1271–5p have been confirmed in plasma samples (see Supplementary Table S3; Table 2).
3.4 MiRs validation in carotid plaque and plasma
Validation through RT-qPCR analysis was performed for the selected miRs in the total cohort of subjects, i.e., 26 asymptomatic and 16 symptomatic inpatients. Carotid plaques of all inpatients were tested by RT-qPCR for single miR analysis of miR-126–5p, miR-1271–5p, miR-145–5p, miR-151a-5p, miR-451a, and miR-720, but results did not confirm data obtained by the miR profiling, and no significant differences were observed (Table 2). Plasma samples were also tested by RT-qPCR for single miR analysis of the above mentioned miRs comparing asymptomatic and symptomatic groups, and borderline and/or significant difference were reported for miR-126–5p (p = 0.0574) and miR-1271–5p (p = 0.0396), being more expressed in symptomatic inpatients (Table 2), in accordance with the trend observed through card/array technology. Conversely, miR-145–5p, miR-151a-5p, miR-451a, and miR-720 showed no significant difference in expression levels between asymptomatic and symptomatic inpatients (Table 2).
Table 2 Relative expression of the selected miRs obtained through RT-qPCR both in plaque biopsies and in plasma from the same inpatients.
3.5 Plasma cytokines and the molecular signature
Among all the assessed cytokines, IL-6, GDF8, and CXCL9 were statistically significant (p = 0.027, p = 0.015, and p = 0.0018, respectively) comparing the two groups, as reported in Figure 2.
Figure 2 CXCL9, GDF8, and IL-6 differ between asymptomatic and symptomatic inpatients. Cytokine concentration was evaluated in plasma samples from asymptomatic inpatients (n = 26) and compared to symptomatic inpatients (n = 16). Data are reported as mean ± SD. Statistical analysis was performed with Student’s t-test. *p ≤ 0.05; **p ≤ 0.01.
MiR-126–5p and miR-1271–5p ΔCt (or expression levels) combined with albumin concentration, percentage of monocytes, CRP, and the three significant cytokines were tested with the leave-one-out cross-validation approach for the best model by means of AUC-ROC values (with standard deviation and p-values), as reported in Table 3. The optimal model was determined to include CXCL9, percentage of monocyte, albumin, and CRP. Together, these variables were capable of distinguishing between symptomatic and asymptomatic inpatients (AUC = 0.83, 95% c.i. [0.85, 0.81], p = 0.028), as depicted in Figure 3, where the two AUC-ROC curves are displayed with and without the two selected miRs.
Table 3 AUC-ROC values with models including or excluding the selected c-miRs (miR-126–5p and miR-1271–5p).
Figure 3 ROC analyses with the model: CXCL9, % of monocytes, albumin, CRP, with and without miR-126–5p and miR-1271–5p. ROC analyses were obtained with a machine learning-based approach (see Methods). A 10× repetition of the analyses resulted in the SD calculation and p-values. The model with CXCL9 and without the identified miRs resulted in the most significant one (AUC = 0.83, 95% c.i. [0.85, 0.81], p = 0.0028).
4 Discussion
Atherosclerosis is a chronic arterial disease that typically manifests more frequently as individuals age, although early stages can be observed in children and teenagers, highlighting the pervasive nature of artery disease development across different age groups (28, 29). This phenomenon poses a significant global health burden.
As highlighted in a recent review by our team (6), carotid plaque lesions can be viewed as tissue-specific manifestations of accelerated aging, akin to chronic pro-inflammatory macro-niches. The molecular and cellular events occurring within these niches encompass well-recognized mechanisms of the aging process, including cellular senescence characterized by a senescence-associated secretory phenotype.
Given these molecular alterations, identifying blood-circulating biomarkers predictive of carotid stenosis is a reasonable challenge. This endeavor is in line with the goals of the current Age.It project, which embraces significant translational power for preventive medicine within the CVD field. The successful identification of such biomarkers could have a substantial impact on healthcare systems, facilitating early detection and intervention strategies to mitigate the burden of CVDs. In fact, our work wants to highlight, as a pilot study, new assessable biomarkers, such as c-miRs and cytokines/chemokines, for carotid artery disease.
As a first approach, the difference between the two groups of asymptomatic and symptomatic inpatients (with carotid occlusion ≥70%) was considered to identify blood molecules shared with different concentrations, and symptomatic inpatients were regarded as those at the highest cardiovascular risk. These molecules will be further assessed as a predictive model in patients with different percentages of carotid stenoses in follow-up monitoring, as proposed in the ongoing National Italian project Age.It. Thus, the molecular signature could likely be applicable for the identification of subjects with a high probability of developing cerebrovascular events in the near future. The simultaneous detection of miRs in the two districts (plaque and blood) could give informative results on circulating miR-based signature in association with a specific plaque phenotype. However, results obtained in the discovery phase were not confirmed in the validation phase, likely due to the increase of variability among inpatients, including the increase of Ct (>32) in some samples (no detection). Similarly, our previous studies, focused on tissue versus blood miR-levels, led to comparable conclusions (25); i.e., card/arrays technology results are often not reproducible in RT-qPCR. This is an additional reason to promote the validation phase. However, the two c-miRs, miR-126–5p and miR-1271–5p, showed similar trends when comparing the two different technologies.
Thus, we identified two c-miRs, i.e., miR-126–5p and miR-1271-5p, differently expressed in plasma from asymptomatic vs. symptomatic inpatients, with both levels being higher in symptomatic individuals. Actually, the two miRs did not correlate with the plaque state (hemorrhagic vs. non-hemorrhagic), but further specific investigations on the various features and parameters of plaques in asymptomatic and symptomatic are currently underway (manuscript in preparation). These investigations are relevant due to the well-recognized role of plaque characteristics in determining cerebrovascular risk (30–33).
As far as miR-126–5p is concerned, it has received less attention than miR-126–3p even if both derive from the same pre-miRNA-126, which maps on human chromosome 9, and it is the principal miR expressed in endothelial cells. MiR-126 has a well-recognized role in the function, integrity, and proliferation of endothelial cells (34–36). In particular, recent data suggest the anti-atherosclerotic role of miR-126–5p through non-canonical, nuclear inhibition of caspase-3, thus conferring endothelial protection from stress through autophagy (37). The current work does not identify differences of miR-126–5p levels between plaques from symptomatic and asymptomatic patients, thus suggesting that higher levels of this miR could be necessary for its function. In this respect, miR-126–5p has been found to be downregulated in human atherosclerotic lesions (34), thus supporting the general hypothesis of the crucial role of the level of miR-126–5p expression. In addition, the-3p has been defined as inflamma-miR (12) and is involved in the regulation of the inflammatory pathway. Interestingly, both c-miR-126–3p and -5p have been proposed as biomarkers of cardiovascular risk in association with stable or vulnerable coronary artery plaque (38). In this perspective, the results here reported are aligned with those previously published, with c-miR-126–5p blood level being higher in symptomatic patients and assuming their highest cardiovascular risk.
As far as c-miR-1271–5p is concerned, a recent article suggests that blood c-miR-1271–5p is involved with other regulatory circular RNA in the progression of thoracic aortic dissection (39), but its role in atherosclerosis and carotid artery disease has yet to be explored.
However, the current results suggest that the two identified c-miRs are not able to fully distinguish (with sensitivity and specificity) the two groups of inpatients, and only by adding other biomarkers is it possible to obtain a molecular signature able to improve the distinction of asymptomatic and symptomatic inpatients.
The identified signature comprises three hemato-biochemical parameters i.e., the percentage of monocyte, albumin, and CRP. Indeed, the comparison of asymptomatic and symptomatic inpatients through hemato-biochemical analysis revealed striking similarities between the two groups. In particular, CRP was not statistically significant comparing directly the two groups, but the machine-learning-based approach was able to select this parameter as informative for the best model, based on the whole dataset.
Among the plethora of various cytokines/molecules assessed within the framework of the ongoing Italian National Project (PNRR-PE8: Age.It), only three of them were deemed significant and able to improve the molecular signature, thereby enhancing the distinction between asymptomatic and symptomatic inpatients. These pivotal molecules include IL-6, GDF8, and CXCL9.
These three cytokines significantly enhance sensitivity and specificity in distinguishing between the two groups. Among them, CXCL9 currently stands out as the most promising chemokine to be included in the signature, alongside the percentage of monocyte, albumin, and CRP (AUC = 0.83, p = 0.0028). This finding suggests the potential exclusion of the two identified miRs.
As far as the three circulating cytokines are concerned, IL-6 is a well-known pro-inflammatory cytokine at the systemic level (40); increased serum level of GDF8 has been recently associated with brachial diastolic pressure and with carotid-femoral pulse wave velocity, a measure of aortic stiffness, in healthy young male adolescents (41). Actually, a recent literature suggests a role of GDF8 in the plaque development (42), and a relationship of plasma GDF8 concentration with chronic kidney disease has been found (43–45). Our data revealed a lower plasma concentration of GDF8 in symptomatic than in asymptomatic inpatients, but activin A/follistatin, recognized to be involved together with GDF8 in the muscle metabolism, were not statistically different between the two groups. Thus, the biological meaning of these results remains to be further explored. As far as CXCL9 is concerned, it is a well-known pro-inflammatory chemokine, also termed monokine induced by gamma interferon (MIG). Interestingly, this chemokine has been associated with old phenotype/poor vascular function and multimorbidity in the so-called inflammatory aging clock (iAge), as the strongest contributor to iAge (46). In addition, the same chemokine has been proposed to be involved in plaque development and its inflammatory environment (24). In this respect, recent research supports the role of the receptor CXCR3 in plaque development. CXCR3, which binds to CXCL9, CXCL10, and CXCL11, has been shown to play a role in human smooth muscle cells in an in vitro model (47), highlighting the potential duality of CXCL9 as both an inflammatory plaque driver and a circulating biomarker. Various chemokines, such as CXCL10 and CCL2, have been extensively investigated for their roles in plaque development over many years, with recent studies also noting the protective role of CXCR4 (48–50). However, the current work focuses on identifying circulating biomarkers, and CXCL9 appears to be a very promising candidate for further testing.
In fact, the next challenge will be to test the identified molecular signature in patients with different carotid stenosis aiming to estimate the prediction of artery occlusion. In fact, the next stage of the current pilot study will be the possibility to assess the identified molecules in an enlarged cohort with monitored outpatients having a different level of carotid stenosis and different comorbidities in the framework of the National Plan of Resilience and Research (PNRR) i.e., Age.it.
Nevertheless, the current article exhibits some weaknesses due to the relatively small number of inpatients, as it is a pilot study and it requires replication in an independent cohort. Additionally, card/array results were not confirmed in RT-qPCR, except the trends of the two c-miRs (miR-126–5p and miR-1271–5p), which were finally validated. On the other hand, the two identified miRs do not appear to be strictly necessary for the molecular signature or for potential application as a screening test. In this regard, the use of c-miRs as biomarkers has yet to gain consensus within the scientific community and particularly concerning CVDs (51). However, other techniques are rapidly advancing to provide more reproducible measurements.
Overall, the challenge of biomarkers for the identification of biological aging (52) and predictive models of age-associated pathology development is still a growing field with expected results not only for each citizen in terms of healthy life span increase, but also for the economic impact on the healthcare system.
Data availability statement
The data supporting the conclusions will be made available upon private request, without undue reservation.
Ethics statement
The studies involving humans were approved by Comitato Etico Area Vasta Emilia Centro (AVEC). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
MCa: Conceptualization, Funding acquisition, Supervision, Writing – original draft. SF: Data curation, Project administration, Writing – review & editing. SCo: Data curation, Investigation, Methodology, Writing – review & editing. EG: Data curation, Visualization, Writing – review & editing. SCa: Data curation, Investigation, Visualization, Writing – review & editing. FF: Data curation, Writing – review & editing. EC: Data curation, Supervision, Writing – review & editing. MCo: Visualization, Writing – review & editing. FO: Supervision, Writing – review & editing. IU: Data curation, Writing – review & editing. RP: Data curation, Writing – review & editing. AV: Data curation, Writing – review & editing. AA: Data curation, Writing – review & editing. FV: Data curation, Writing – review & editing. GL: Data curation, Supervision, Writing – review & editing. GP: Methodology, Supervision, Writing – original draft. MG: Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors acknowledge co-funding from Next Generation EU, in the context of the National Recovery and Resilience Plan, Investment PE8 – Project Age-It: “Ageing Well in an Ageing Society”. This resource was co-financed by the Next Generation EU (DM 1557 11.10.2022). The authors also acknowledge co-funding from FONDAZIONE CARISBO for CAROMIRNA-19 (n 2019.0538) to MC.
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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.
Author disclaimer
The views and opinions expressed are only those of the authors and do not necessarily reflect those of the European Union or the European Commission. Neither the European Union nor the European Commission can be held responsible for them.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1407396/full#supplementary-material
References
1. Song P, Fang Z, Wang H, Cai Y, Rahimi K, Zhu Y, et al. Global and regional prevalence, burden, and risk factors for carotid atherosclerosis: a systematic review, meta-analysis, and modelling study. Lancet Global Health. (2020) 8:e721–9. doi: 10.1016/S2214–109X(20)30117–0
2. Johnsen SH, Mathiesen EB, Joakimsen O, Stensland E, Wilsgaard T, Løchen M-L, et al. Carotid atherosclerosis is a stronger predictor of myocardial infarction in women than in men: A 6-year follow-up study of 6226 persons: the Tromsø Study. Stroke. (2007) 38:2873–80. doi: 10.1161/STROKEAHA.107.487264
3. Inaba Y, Chen JA, Bergmann SR. Carotid plaque, compared with carotid intima-media thickness, more accurately predicts coronary artery disease events: A meta-analysis. Atherosclerosis. (2012) 220:128–33. doi: 10.1016/j.atherosclerosis.2011.06.044
4. Touboul P-J, Hennerici MG, Meairs S, Adams H, Amarenco P, Bornstein N, et al. Mannheim carotid intima-media thickness and plaque consensus (2004–2006–2011). Cerebrovasc Dis. (2012) 34:290–6. doi: 10.1159/000343145
5. Piepoli MF, Hoes AW, Agewall S, Albus C, Brotons C, Catapano AL, et al. 2016 European Guidelines on cardiovascular disease prevention in clinical practice: The Sixth Joint Task Force of the European Society of Cardiology and Other Societies on Cardiovascular Disease Prevention in Clinical Practice (constituted by representatives of 10 societies and by invited experts)Developed with the special contribution of the European Association for Cardiovascular Prevention & Rehabilitation (EACPR). Eur Heart J. (2016) 37:2315–81. doi: 10.1093/eurheartj/ehw106
6. Collura S, Morsiani C, Vacirca A, Fronterrè S, Ciavarella C, Vasuri F, et al. The carotid plaque as paradigmatic case of site-specific acceleration of aging process: The microRNAs and the inflammaging contribution. Ageing Res Rev. (2020) 61:101090. doi: 10.1016/j.arr.2020.101090
7. Collura S, Ciavarella C, Morsiani C, Motta I, Valente S, Gallitto E, et al. MicroRNA profiles of human peripheral arteries and abdominal aorta in normal conditions: MicroRNAs-27a-5p, -139–5p and -155–5p emerge and in atheroma too. Mech Ageing Dev. (2021) 198:111547. doi: 10.1016/j.mad.2021.111547
8. Ciavarella C, Motta I, Capri M, Gargiulo M, Pasquinelli G. Heterogeneity and differentiation of the human arterial tree: focus on microRNA expression in vascular disease. Biomolecules. (2024) 14:343. doi: 10.3390/biom14030343
9. Eken SM, Jin H, Chernogubova E, Li Y, Simon N, Sun C, et al. MicroRNA-210 enhances fibrous cap stability in advanced atherosclerotic lesions. Circ Res. (2017) 120:633–44. doi: 10.1161/CIRCRESAHA.116.309318
10. Maitrias P, Metzinger-Le Meuth V, Massy ZA, M’Baya-Moutoula E, Reix T, Caus T, et al. MicroRNA deregulation in symptomatic carotid plaque. J Vasc Surg. (2015) 62:1245–1250.e1. doi: 10.1016/j.jvs.2015.06.136
11. Cipollone F, Felicioni L, Sarzani R, Ucchino S, Spigonardo F, Mandolini C, et al. A unique microRNA signature associated with plaque instability in humans. Stroke. (2011) 42:2556–63. doi: 10.1161/STROKEAHA.110.597575
12. Olivieri F, Rippo MR, Monsurrò V, Salvioli S, Capri M, Procopio AD, et al. MicroRNAs linking inflamm-aging, cellular senescence and cancer. Ageing Res Rev. (2013) 12:1056–68. doi: 10.1016/j.arr.2013.05.001
13. Olivieri F, Rippo MR, Procopio AD, Fazioli F. Circulating inflamma-miRs in aging and age-related diseases. Front Genet. (2013) 4:121. doi: 10.3389/fgene.2013.00121
14. Olivieri F, Capri M, Bonafè M, Morsiani C, Jung HJ, Spazzafumo L, et al. Circulating miRNAs and miRNA shuttles as biomarkers: Perspective trajectories of healthy and unhealthy aging. Mech Ageing Dev. (2017) 165:162–70. doi: 10.1016/j.mad.2016.12.004
15. Morsiani C, Terlecki-Zaniewicz L, Skalicky S, Bacalini MG, Collura S, Conte M, et al. Circulating miR-19a-3p and miR-19b-3p characterize the human aging process and their isomiRs associate with healthy status at extreme ages. Aging Cell. (2021) 20:e13409. doi: 10.1111/acel.13409
16. Morsiani C, Bacalini MG, Collura S, Moreno-Villanueva M, Breusing N, Bürkle A, et al. Blood circulating miR-28–5p and let-7d-5p associate with premature ageing in Down syndrome. Mech Ageing Dev. (2022) 206:111691. doi: 10.1016/j.mad.2022.111691
17. Morsiani C, Collura S, Sevini F, Ciurca E, Bertuzzo VR, Franceschi C, et al. Circulating miR-122–5p, miR-92a-3p, and miR-18a-5p as Potential Biomarkers in Human Liver Transplantation Follow-Up. IJMS. (2023) 24:3457. doi: 10.3390/ijms24043457
18. Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci U S A. (2011) 108:5003–8. doi: 10.1073/pnas.1019055108
19. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. (2011) 13:423–33. doi: 10.1038/ncb2210
20. Prapiadou S, Živković L, Thorand B, George MJ, van der Laan SW, Malik R, et al. Proteogenomic data integration reveals CXCL10 as a potentially downstream causal mediator for IL-6 signaling on atherosclerosis. Circulation. (2024) 149:669–83. doi: 10.1161/CIRCULATIONAHA.123.064974
21. Georgakis MK, Bernhagen J, Heitman LH, Weber C, Dichgans M. Targeting the CCL2–CCR2 axis for atheroprotection. Eur Heart J. (2022) 43:1799–808. doi: 10.1093/eurheartj/ehac094
22. Döring Y, van der Vorst EPC, Duchene J, Jansen Y, Gencer S, Bidzhekov K, et al. CXCL12 derived from endothelial cells promotes atherosclerosis to drive coronary artery disease. Circulation. (2019) 139:1338–40. doi: 10.1161/CIRCULATIONAHA.118.037953
23. Döring Y, Noels H, van der Vorst EPC, Neideck C, Egea V, Drechsler M, et al. Vascular CXCR4 limits atherosclerosis by maintaining arterial integrity: evidence from mouse and human studies. Circulation. (2017) 136:388–403. doi: 10.1161/CIRCULATIONAHA.117.027646
24. Zernecke A, Shagdarsuren E, Weber C. Chemokines in atherosclerosis: an update. ATVB. (2008) 28:1897–908. doi: 10.1161/ATVBAHA.107.161174
25. Kangas R, Morsiani C, Pizza G, Lanzarini C, Aukee P, Kaprio J, et al. Menopause and adipose tissue: miR-19a-3p is sensitive to hormonal replacement. Oncotarget. (2017) 9:2279–94. doi: 10.18632/oncotarget.23406
26. Buuren SV, Groothuis-Oudshoorn K. Multivariate imputation by chained equations in. R J Stat Soft. (2011) 45:1–67. doi: 10.18637/jss.v045.i03
27. Ferri FJ, Pudil P, Hatef M, Kittler J. Comparative study of techniques for large-scale feature selection* *This work was suported by a SERC grant GR/E 97549. The first author was also supported by a FPI grant from the Spanish MEC, PF92 73546684. In: Machine Intelligence and Pattern Recognition. North-Holland: Elsevier (1994). p. 403–13. doi: 10.1016/B978–0-444–81892-8.50040–7
28. Hong YM. Atherosclerotic cardiovascular disease beginning in childhood. Korean Circ J. (2010) 40:1. doi: 10.4070/kcj.2010.40.1.1
29. McGill HC, McMahan CA, Herderick EE, Tracy RE, Malcom GT, Zieske AW, et al. Effects of coronary heart disease risk factors on atherosclerosis of selected regions of the aorta and right coronary artery. ATVB. (2000) 20:836–45. doi: 10.1161/01.ATV.20.3.836
30. Pini R, Faggioli G, Fittipaldi S, Vasuri F, Longhi M, Gallitto E, et al. Relationship between calcification and vulnerability of the carotid plaques. Ann Vasc Surg. (2017) 44:336–42. doi: 10.1016/j.avsg.2017.04.017
31. Zazzeroni L, Faggioli G, Pasquinelli G. Mechanisms of arterial calcification: the role of matrix vesicles. Eur J Vasc Endovascular Surg. (2018) 55:425–32. doi: 10.1016/j.ejvs.2017.12.009
32. Vasuri F, Ciavarella C, Fittipaldi S, Pini R, Vacirca A, Gargiulo M, et al. Different histological types of active intraplaque calcification underlie alternative miRNA-mRNA axes in carotid atherosclerotic disease. Virchows Arch. (2020) 476:307–16. doi: 10.1007/s00428-019-02659-w
33. Vasuri F, De Biase D, Vacirca A, Acquaviva G, Sanza V, Gargiulo M, et al. Gene polymorphism in tissue epidermal growth factor receptor (EGFR) influences clinical and histological vulnerability of carotid plaques. Pathol - Res Pract. (2022) 229:153721. doi: 10.1016/j.prp.2021.153721
34. Schober A, Nazari-Jahantigh M, Wei Y, Bidzhekov K, Gremse F, Grommes J, et al. MicroRNA-126–5p promotes endothelial proliferation and limits atherosclerosis by suppressing Dlk1. Nat Med. (2014) 20:368–76. doi: 10.1038/nm.3487
35. Harris TA, Yamakuchi M, Ferlito M, Mendell JT, Lowenstein CJ. MicroRNA-126 regulates endothelial expression of vascular cell adhesion molecule 1. Proc Natl Acad Sci U S A. (2008) 105:1516–21. doi: 10.1073/pnas.0707493105
36. Fish JE, Santoro MM, Morton SU, Yu S, Yeh R-F, Wythe JD, et al. miR-126 regulates angiogenic signaling and vascular integrity. Dev Cell. (2008) 15:272–84. doi: 10.1016/j.devcel.2008.07.008
37. Santovito D, Egea V, Bidzhekov K, Natarelli L, Mourão A, Blanchet X, et al. Noncanonical inhibition of caspase-3 by a nuclear microRNA confers endothelial protection by autophagy in atherosclerosis. Sci Transl Med. (2020) 12:eaaz2294. doi: 10.1126/scitranslmed.aaz2294
38. Rozhkov AN, Shchekochikhin D, Ashikhmin YI, Mitina YO, Evgrafova VV, Zhelankin AV, et al. The profile of circulating blood microRNAs in outpatients with vulnerable and stable atherosclerotic plaques: association with cardiovascular risks. ncRNA. (2022) 8:47. doi: 10.3390/ncrna8040047
39. Xie X, Hong X, Hong S, Huang Y, Chen G, Chen Y, et al. Progression of Thoracic Aortic Dissection Is Aggravated by the hsa_circ_0007386/miR-1271–5P/IGF1R/AKT Axis via Induction of Arterial Smooth Muscle Cell Apoptosis. Biomedicines. (2023) 11:571. doi: 10.3390/biomedicines11020571
40. Forcina L, Franceschi C, Musarò A. The hormetic and hermetic role of IL-6. Ageing Res Rev. (2022) 80:101697. doi: 10.1016/j.arr.2022.101697
41. Pucci G, Ministrini S, Nulli Migliola E, Nunziangeli L, Battista F, D’Abbondanza M, et al. Relationship between serum myostatin levels and carotid-femoral pulse wave velocity in healthy young male adolescents: the MACISTE study. J Appl Physiol. (2021) 130:987–92. doi: 10.1152/japplphysiol.00782.2020
42. Verzola D, Milanesi S, Bertolotto M, Garibaldi S, Villaggio B, Brunelli C, et al. Myostatin mediates abdominal aortic atherosclerosis progression by inducing vascular smooth muscle cell dysfunction and monocyte recruitment. Sci Rep. (2017) 7:46362. doi: 10.1038/srep46362
43. Kuo C-S, Lu Y-W, Hsu C-Y, Chang C-C, Chou R-H, Liu L-K, et al. Increased activin A levels in prediabetes and association with carotid intima-media thickness: a cross-sectional analysis from I-Lan Longitudinal Aging Study. Sci Rep. (2018) 8:9957. doi: 10.1038/s41598–018-27795–2
44. Yano S, Nagai A, Isomura M, Yamasaki M, Kijima T, Takeda M, et al. Relationship between blood myostatin levels and kidney function : Shimane CoHRE study. PloS One. (2015) 10:e0141035. doi: 10.1371/journal.pone.0141035
45. Esposito P, Verzola D, Picciotto D, Cipriani L, Viazzi F, Garibotto G. Myostatin/activin-A signaling in the vessel wall and vascular calcification. Cells. (2021) 10:2070. doi: 10.3390/cells10082070
46. Sayed N, Huang Y, Nguyen K, Krejciova-Rajaniemi Z, Grawe AP, Gao T, et al. An inflammatory aging clock (iAge) based on deep learning tracks multimorbidity, immunosenescence, frailty and cardiovascular aging. Nat Aging. (2021) 1:598–615. doi: 10.1038/s43587-021-00082-y
47. Hayderi A, Kumawat AK, Shavva VS, Dreifaldt M, Sigvant B, Petri MH, et al. RSAD2 is abundant in atherosclerotic plaques and promotes interferon-induced CXCR3-chemokines in human smooth muscle cells. Sci Rep. (2024) 14:8196. doi: 10.1038/s41598–024-58592–9
48. Sjaarda J, Gerstein H, Chong M, Yusuf S, Meyre D, Anand SS, et al. Blood CSF1 and CXCL12 as causal mediators of coronary artery disease. J Am Coll Cardiol. (2018) 72:300–10. doi: 10.1016/j.jacc.2018.04.067
49. Winter C, Silvestre-Roig C, Ortega-Gomez A, Lemnitzer P, Poelman H, Schumski A, et al. Chrono-pharmacological targeting of the CCL2-CCR2 axis ameliorates atherosclerosis. Cell Metab. (2018) 28:175–82.e5. doi: 10.1016/j.cmet.2018.05.002
50. Bot I, Daissormont ITMN, Zernecke A, Van Puijvelde GHM, Kramp B, De Jager SCA, et al. CXCR4 blockade induces atherosclerosis by affecting neutrophil function. J Mol Cell Cardiol. (2014) 74:44–52. doi: 10.1016/j.yjmcc.2014.04.021
51. Wronska A, Kurkowska-Jastrzebska I, Santulli G. Application of micro RNA s in diagnosis and treatment of cardiovascular disease. Acta Physiologica. (2015) 213:60–83. doi: 10.1111/apha.12416
Keywords: microRNAs, carotid artery disease, biomarkers, asymptomatic patients, symptomatic patients, elders, CXCL9, IL-6
Citation: Capri M, Fronterrè S, Collura S, Giampieri E, Carrino S, Feroldi FM, Ciurca E, Conte M, Olivieri F, Ullo I, Pini R, Vacirca A, Astolfi A, Vasuri F, La Manna G, Pasquinelli G and Gargiulo M (2024) Circulating CXCL9, monocyte percentage, albumin, and C-reactive protein as a potential, non-invasive, molecular signature of carotid artery disease in 65+ patients with multimorbidity: a pilot study in Age.It. Front. Endocrinol. 15:1407396. doi: 10.3389/fendo.2024.1407396
Received: 26 March 2024; Accepted: 31 May 2024;
Published: 23 July 2024.
Edited by:
Francesco Prattichizzo, MultiMedica Holding SpA (IRCCS), ItalyReviewed by:
Elettra Mancuso, University Magna Graecia of Catanzaro, ItalyDonato Santovito, Ludwig Maximilian University of Munich, Germany
Copyright © 2024 Capri, Fronterrè, Collura, Giampieri, Carrino, Feroldi, Ciurca, Conte, Olivieri, Ullo, Pini, Vacirca, Astolfi, Vasuri, La Manna, Pasquinelli and Gargiulo. 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: Miriam Capri, miriam.capri@unibo.it
†These authors share senior authorship