Next Article in Journal
The Prevalence of Sjögren’s Disease in Dental Clinics in the Netherlands Compared with the Prevalence in a Systematic Literature Review of Studies in Other Countries
Previous Article in Journal
Extraperitoneal Robot-Assisted Radical Prostatectomy with the Hugo™ RAS System: Initial Experience at a High-Volume Robotic Centre
Previous Article in Special Issue
Self-Management among Stroke Survivors in the United States, 2016 to 2021
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 13
Issue 19
10.3390/jcm13195917
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving Stroke Outcome Prediction Using Molecular and Machine Learning Approaches in Large Vessel Occlusion

by
Madhusmita Rout
1,
April Vaughan
2,
Evgeny V. Sidorov
2,* and
Dharambir K. Sanghera
1,3,4,5,6,*
1
Department of Pediatrics, College of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
2
Department of Neurology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
3
Department of Pharmaceutical Sciences, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
4
Department of Physiology, College of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
5
Oklahoma Center for Neuroscience, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
6
Harold Hamm Diabetes Center, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104, USA
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2024, 13(19), 5917; https://doi.org/10.3390/jcm13195917
Submission received: 5 September 2024 / Revised: 25 September 2024 / Accepted: 1 October 2024 / Published: 3 October 2024
(This article belongs to the Special Issue Stroke Diagnosis and Outcome Prediction)
Browse Figures
Versions Notes

Abstract

:
Introduction: Predicting stroke outcomes in acute ischemic stroke (AIS) can be challenging, especially for patients with large vessel occlusion (LVO). Available tools such as infarct volume and the National Institute of Health Stroke Scale (NIHSS) have shown limited accuracy in predicting outcomes for this specific patient population. The present study aimed to confirm whether sudden metabolic changes due to blood-brain barrier (BBB) disruption during LVO reflect differences in circulating metabolites and RNA between small and large core strokes. The second objective was to evaluate whether integrating molecular markers with existing neurological and imaging tools can enhance outcome predictions in LVO strokes. Methods: The infarction volume in patients was measured using magnetic resonance diffusion-weighted images, and the 90-day stroke outcome was defined by a modified Rankin Scale (mRS). Differential expression patterns of miRNAs were identified by RNA sequencing of serum-driven exosomes. Nuclear magnetic resonance (NMR) spectroscopy was used to identify metabolites associated with AIS with small and large infarctions. Results: We identified 41 miRNAs and 11 metabolites to be significantly associated with infarct volume in a multivariate regression analysis after adjusting for the confounders. Eight miRNAs and ketone bodies correlated significantly with infarct volume, NIHSS (severity), and mRS (outcome). Through integrative analysis of clinical, radiological, and omics data using machine learning, our study identified 11 top features for predicting stroke outcomes with an accuracy of 0.81 and AUC of 0.91. Conclusions: Our study provides a future framework for advancing stroke therapeutics by incorporating molecular markers into the existing neurological and imaging tools to improve predictive efficacy and enhance patient outcomes.

1. Introduction

Stroke is the third-leading cause of death and disability worldwide, and acute ischemic stroke (AIS) accounts for 85% of all strokes. According to estimates from the Global Burden of Diseases, Injuries, and Risk Factors Study, AIS mortality led to approximately 6.55 million deaths worldwide from 1990 to 2019 [1]. Large vessel occlusions (LVOs) in the AIS occur due to sudden interruption of blood flow in the large artery due to atherothrombosis or cardiac embolism, which leads to rapid death of neurons and neurological deterioration [2,3,4]. LVOs are identified based on occlusion in the internal carotid artery, the first and second segments of the middle cerebral artery, the basilar artery, the vertebral artery, the first and second segments of the posterior cerebral artery, and the first segment of the anterior cerebral artery [5]. Patients with LVOs constitute a significant subpopulation of AIS, and their annual incidence in the United States is around 80,000 [6].
At present, stroke is diagnosed based on clinical assessments and radiological exams by magnetic resonance imaging (MRI) and computed tomography (CT), along with currently available treatment therapies, including reperfusion and endovascular thrombectomy (EVT). The severity of a stroke is typically assessed using the National Institute of Health Stroke Scale (NIHSS) scores, while the outcome is evaluated with the modified Rankin Scale (mRS) [7]. The infarct volume is an important determinant of long-term consequences in the AIS, leading to extensive neuronal damage and subsequent neuronal deficit. Even though the infarct volume strongly correlates with clinical outcomes, it cannot explain why some patients after large infarctions demonstrate significant improvement while others develop poor outcomes with small strokes. There is currently limited understanding of the connection between infarct volume, severity, and outcomes [8].
Recent improvements in high-throughput metabolomic technologies (liquid-chromatography and mass spectroscopy (LCMS) and nuclear magnetic resonance (NMR)) and transcriptomics have made it possible to characterize molecular changes in multiple tissues of a large number of patients at the same time. Noncoding regions of the genome encoding RNA, microRNAs (miRNAs), play an important role in regulating the metabolic homeostasis of different diseases by influencing the regulation of target genes [9]. MiRNAs have been implicated in the disruption of the blood–brain barrier (BBB) and induce cellular changes in AIS, including mitochondrial dysfunction, cytokine-mediated cytotoxicity, oxidative stress, glial cell activation, and inflammation [10]. Thus far, the majority of biomarker discoveries in AIS have mainly been carried out in rodents utilizing middle cerebral artery occlusion (MCAO) models. However, their predictability in human strokes has proven to be challenging [11,12]. There is a need to identify chemical/biological markers to improve stroke prognostics and accurately predict the stroke outcomes following EVT or thrombolysis. During the acute phase of the LVO, the aberrant expression of the extracellular transcriptome and alterations in metabolite species in the brain regulate neuroinflammation, causing BBB interruption. The present study aimed to confirm whether sudden metabolic changes due to BBB disruption reflect differences in circulating metabolites and RNA between small and large core strokes. The second objective was to utilize machine learning as a predictive tool to integrate clinical, radiological imaging, transcriptomic (miRNA), and metabolomic data to determine functional outcomes in patients with LVO strokes.

2. Experimental Procedure

2.1. Study Participants

The study comprised 85 LVO patients who were part of our ongoing Metabolome in Ischemic Stroke Study (MISS), which began in 2017 for recruiting patients from the Comprehensive Stroke Center at the Oklahoma University Medical Center. The MISS is designed to identify biomarker predictors of AIS and other strokes using metabolomics and other omics technologies in longitudinal and cross-sectional study cohorts [13,14,15,16]. The blood samples were collected during acute (<72 h of stroke) and follow-up phases (within 3 months). All patients used for this cross-sectional study underwent MRI with diffusion-weighted imaging (DWI) sequence during admission for AIS definition. Hemorrhagic stroke was classified based on subarachnoid and intracerebral hemorrhage. Patients were evaluated upon admission with the NIHSS (0–42 points reflecting the least to the most severe neurological deficit) and then with the Miami Emergency Neurological Deficit examination to monitor for signs of neurological deterioration during hospital admission at the acute stage [17,18]. The inclusion criteria comprised (i) ischemic strokes with LVO with a large ischemic core (>40 mL), (ii) a small ischemic core (<40 mL) admitted within 24 h of symptom onset, and (iii) patients with atrial fibrillation, hypertension, obesity, and previous stroke. We excluded those patients with (i) hemorrhagic conversion of ischemic stroke and formation of hematoma defined by the European Cooperative Acute Stroke Study [19], (ii) fluctuating neurological exams during the acute stage, (iii) systemic infection at presentation or during admission with fever >38.0 °C, elevated white blood cells, diagnosis of pneumonia, or urinary tract infection, (iv) renal disease with a glomerular filtration rate (GFR) <45; and (v) recurrent stroke in the chronic stage that reduces confounding due to these conditions as described [15]. Of the 85 LVO cases, 14% had atrial fibrillation, 64% had hypertension, 40% were obese (BMI > 30), 45% were diabetic, 19% had coronary artery disease (CAD), and 24% had a previous history of AIS. The University of Oklahoma Health Sciences Center’s Institutional Review Board approved all study protocols. Figure 1 depicts the workflow of the study.

2.2. Neuroimaging, Infarct Volume and Clinical Outcome

Two independent investigators measured stroke (ischemic core) volume using MR diffusion-weighted images 12 to 48 h after admission using MR segmentation software [15]. As described in our earlier study, this software utilizes an interactive computer-aided detection (CAD) scheme. The user identifies and draws the stroke boundary on each of the identified stroke slices, and the CAD scheme computes the resulting infarction volume [15]. The infarct volume cut-off was used to categorize small (≤40; N = 42) and large (>40; N = 43) infarct volumes. The mRS assesses the clinical outcome of AIS patients after the stroke [20,21]. It grades the patient’s disability from 0 (no symptoms) to 6 (death), thus capturing the whole spectrum of functional states 30 and 90 days following stroke [20,22,23].

2.3. Metabolome Analysis

Aliquots of frozen serum from MISS patients were shipped on dry ice to Nightingale Health Ltd., Helsinki, Finland (https://nightingalehealth.com/), to measure lipidomics profiling using targeted high-throughput NMR. This platform provided simultaneous quantification of 250 metabolic measures, routine lipids, lipoprotein subclass profiling with lipid concentrations within 14 subclasses, esterified fatty acid composition, and various low molecular metabolites, including amino acids, ketone bodies, and gluconeogenesis-related metabolites in molar concentration units as described previously [24,25].

2.4. Exosome Extraction and RNA Isolation

The blood sample collected in anti-coagulant-free tubes was centrifuged for 10 min at 1258× g (4 °C), and the supernatant (serum) was transferred into a fresh tube and stored at –80 °C until further analysis. Exosomal miRNA quantification and differential expression analysis were performed using serum specimens. Following the manufacturer’s protocol, the Serum/Plasma Exosome Kit (Qiagen Inc., Chatsworth, CA, USA) was used. We first characterized the particle size morphology of the exosomes extracted from the sera of patients using transmission electron microscopy and via nanoparticle tracking analysis (NTA) using Nanosight 300 as described earlier [26]. Further, the exosome pellet was resuspended in a resuspension buffer to extract the total RNA. RNA was extracted and purified using the miRNeasy mini kit (Qiagen Inc., Chatsworth, CA, USA) and ExoQuick Exosome Isolation and RNA Purification Kit (System Biosciences, Palo Alto, CA, USA) according to the manufacturer’s protocol.

2.5. miRNA Sequencing Quality Control (QC) and Analysis

RNA concentrations were estimated and validated for their quantity and integrity using Agilent’s autoanalyzer. RNA quality was determined using a tape station. The nextGen sequencing libraries were created using QIAGEN’s QIAseq miRNA library kit, which incorporates Unique Molecular Indices, and the QC of libraries was performed using Agilent’s bioanalyzer to determine library size and quality. Sequencing was performed on an Illumina sequencer using single-end, 75 bp reads to a depth of 12–14 million reads at the Genomics Core facility of the Oklahoma Medical Research Foundation (OMRF). Sequence reads were pre-processed to remove adapters and were aligned to the reference human genome (hg19) using the Bowtie aligner as described [27,28]. The reads aligning to each known mature miRNA were counted using Bioconductor packages for next-generation sequencing data analysis [29] based on miRNA definitions in the miRBase database [30]. Sequencing analysis was performed using QIAGEN’s CLC workbench.

2.6. Statistical Analysis

The clinical and demographic characteristics of the study participants were analyzed using the phenotypic mean for continuous variables and percentages for categorical variables. We characterized and prioritized differentially expressed miRNAs correlating with the acute stroke using the normalized reads and a regression approach (R package DESeq2). A series of standardization steps were performed to make the expression phenotypes comparable across individuals and across miRNAs. To minimize the influence of overall signal levels, abundance values of all retained miRNAs were standardized within individual z-scores, followed by linear regression against individual-specific average log-transformed raw signals and their squared value. Before the statistical analysis, the raw miRNA and metabolite profile data were normalized through log transformation, which resulted in some negative values and was scaled by MetaboAnalyst 5.0 (www.metaboanalyst.ca/). We first used two-tailed, unpaired Wilcoxon rank-sum tests with p < 0.05 to identify metabolites and miRNAs that differed significantly in small and large infarct volumes using MetaboAnalyst. A false discovery rate (FDR) was implemented to correct for the multiple comparisons. Multiple linear regression analyses were performed to assess the impact of individual miRNA and metabolite markers on infarct volume after adjusting for covariates such as age, gender, ethnicity, body mass index (BMI), smoking, etc. Pearson correlation analysis was performed to determine the linear association between miRNAs and metabolites and the association of mRS and NIHSS scores with miRNAs and metabolites. All analyses were performed using SVS version 8.9.1 (Golden Helix, Bozeman, MT, USA) and SPSS version 29. For a sample size of 43 patients with large-volume strokes and 42 patients with small-volume strokes, with a significance level (alpha) of 2.0 × 10−4 and 250 metabolites, our study has 80% power to detect a mean difference of 0.58 standard deviation (SD) between small- and large-volume groups. After assuming conservative multi-testing correction for ~2600 miRNA transcripts at a significant level of 1.92 × 10−5 for the same patients, our study will achieve ~75% power to detect a mean difference of 1 SD to identify miRNAs differing between groups using a two-sided two-sample equal-variance t-test.

2.7. Machine Learning

Machine learning was used as a predictive tool to integrate clinical, miRNAs, and metabolomics biomarker data differing among 68 stroke cases in good (mRS = 0–3) and bad (mRS = 3–6) stroke outcomes. A binary classification method was employed to determine stroke outcomes in Python (v3.11.5) using sci-kit-learn [31] and pandas libraries. In binary classification, class labels were determined by stroke cases with mRS (≤3; N = 37) characterized as a good outcome and those with mRS (>3; N = 31) characterized as a bad outcome. Chi-Square and F-regression tests were used for selecting features, and the Shapley Additive exPlanations (SHAP) framework from Slundberg (https://github.com/slundberg/shap accessed on 21 June 2024) was used to identify important features using an XGBoost model with a Python 3 kernel (v3.11.5). Summary plots depict the top influential biomarkers and how they influence the model prediction [32,33]. Among the 216 miRNAs and 250 metabolites, only those features that were highly dependent on the outcome were selected. The 11 selected features were scaled as initial QC to improve the performance and accuracy of machine learning algorithms. Scaling ensures all features contribute equally by converting them to a common scale to eliminate potential biases. The data file was split into 80% training and 20% testing partitions using a defined seed value. Seeds were only chosen such that the resulting training and testing accuracies were similar, preventing overfit to the small sample size. Tenfold cross-validation was used to evaluate the model and detect overfitting. Supervised learning classifiers such as Logistic Regression (LR), Linear Discriminant Analysis (LDA), K-Nearest Neighbors Classifier (KNN), Gaussian Naïve Bayes (NB), Support Vector Machine (SVM) and a deep learning method, Artificial Neural Network (ANN), were used to classify the stroke outcome. The model was evaluated using various classification metrics, such as accuracy, which is defined as the number of correct predictions divided by the total number of predictions. The receiver operator characteristics (ROC) curve is a measure of sensitivity and specificity at different decision thresholds, and the area under the curve (AUC) offers a single metric for summarizing the overall performance of the model.

3. Results

The clinical and demographic characteristics of the study participants are presented in Table 1. There were no significant differences in the average age, BMI, and diastolic blood pressure between large and small infarct volume cases. Plasma glucose was significantly elevated (162.9 ± 96.7) in cases with large infarcts compared to small infarcts (128.5 ± 59.3). Systolic blood pressure was significantly higher in large (158.1 ± 32.7) vs. small strokes (145.7 ± 21.9). The study included predominantly Caucasian participants (93% and 74%), with smaller percentages of African Americans (7% and 14%) in the large and small infarct volume groups, respectively. The prevalence of type 2 diabetes was 53% in individuals with large infarct volumes and 36% in those with small infarct volumes. Smoking was prevalent in patients with small infarcts (50%) compared to those with large infarcts (44%). The average infarct volume differed significantly between large (117.4 ± 69.1) and small (19.1 ± 14.1) strokes (p = 7.43 × 10−12). Notably, patients with a large infarct volume had significantly worse outcomes, as indicated by the average mRS of (4.2 ± 1.9) compared to the mRS of (2.6 ± 1.9) of those with a small infarct volume (p = 0.0004) (Table 1). Incidentally, our data revealed a stronger correlation between the infarct volume and mRS (r = 0.41; p = 8.2 × 10−4) than the infarct volume and NIHSS (r = 0.32; p = 0.004) (Figure 2).

3.1. Discovery of Significantly Altered Patterns of Circulating miRNAs Associated with Infarct Volume

Sequencing analysis of exosomal miRNAs identified 2632 total miRNAs in our samples. Of these, 41 miRNAs were differentially expressed between large and small infarct volumes. Most of the miRNAs were significantly negatively regulated in large infarcts except for a few, which were upregulated in large infarcts: miR-21-5p (t.stat = 3.5; p = 7.1 × 10−4), miR-101-3p (t.stat = 3.8; p = 2.5 × 10−4), and miR-16-5p (t.stat = 2.3; p = 0.03) (Supplementary Table S1; Figure 3). Selected miRNAs such as miR-144-5p (r = −0.67; p = 2.5 × 10−10), miR-486-5p (r = −0.44; p = 1.8 × 10−4), let-7f-5p (r = −0.45; p = 1.2 × 10−4) and miR-96-5p (r = −0.44; p = 1.3 × 10−4) showed a strong negative correlation with infarct volume (Figure 3).

3.2. Association of miRNAs with NIHSS and Stroke Outcome (mRS Score)

The miRNAs that were significantly associated with high infarct volume also correlated with NIHSS and stroke outcomes (mRS). The majority of the miRNAs showed a significant negative correlation with NIHSS scores. The top ones were miR-625-3p (r = −0.51; p = 9.8 × 10−7) and miR-486-5p (r = −0.41; p = 1.8 × 10−4). However, miR-16-5p (r = 0.52; p = 8.5 × 10−7) displayed a positive correlation with NIHSS (Supplementary Table S2). Additionally, eight miRNAs were negatively correlated with increased mRS (Supplementary Table S3).

3.3. Discovery of Significantly Altered Patterns of Serum Metabolites with Infarct Volume

Selected metabolites such as glucose (r = 0.42; p = 3.6 × 10−4) and acetoacetate (r = 0.33; p = 0.006) were positively associated with infarct volume (Supplementary Figure S1). Multivariate regression analysis for infarct volume identified 11 metabolites, mostly lipids, to be associated with stroke volume (Table 2). Glucose (Beta (SE): 0.51 (0.17); p = 3.2 × 10−3), free cholesterol to total lipids ratio in small HDL (S-HDL-FC%) (Beta (SE): 0.53 (0.18); p = 4.6 × 10−3) and ketone bodies such as acetoacetate (Beta (SE): 0.36 (0.16); p = 0.03) were increased in patients with a high stroke volume (Table 2).

3.4. Association of Metabolites with NIHSS and Stroke Outcomes (mRS)

Since NIHSS correlated with the infarction volume in this study, it was not surprising to find metabolites associated with the infarct volume also correlated with NIHSS. Metabolites of ketone bodies such as acetoacetate (r = 0.35; p = 5.2 × 10−3) and lipid S-HDL-FC% (r = 0.29; p = 0.02) showed a positive significant correlation with NIHSS, while glycerol (r = −0.31; p = 0.01) and albumin (r = −0.27; p = 0.03) correlated negatively with NIHSS (Supplementary Table S4). Interestingly, ketone bodies acetoacetate (r = 0.43; p = 1.2 × 10−3) and 3-hydroxybutyrate (r = 0.34; p = 0.01) were positively correlated with mRS, depicting poor clinical outcomes (Supplementary Table S5).

3.5. Association of Exosomal miRNAs with Metabolites

We identified 41 miRNAs that were significantly differentially regulated with respect to infarct volume with exceptions for miR-21-5p, miR-27a-3p, miR-29c-3p, miR-101-3p, and miR-16-5p (Figure 3 and Figure 4). In total, 21 metabolites showed a significant correlation with 41 miRNAs. Amino acids such as alanine (Ala), glycine (Gly), and histidine (His) correlated positively with miRNAs, whereas glycoprotein acetyls (GlycA), glutamine (Gln), and phenylalanine (Phe) were negatively correlated with the miRNAs that were down-regulated in large infarct strokes. Interestingly, the upregulated miRNAs (miR-101-3p, miR-16-5p, and miR-21-5p) in large infarct strokes correlated positively with ketone bodies, glucose, apolipoprotein B (ApoB) and phospholipids in small LDL (S-LDL-PL) (Supplementary Figure S2).
Overall, eight miRNAs (miR-485-3p, miR-625-3p, miR-324-5p, miR-139-3p, miR-140-3p, miR-127-3p, miR-501-3p and miR-144-5p) and one metabolite (acetoacetate) correlated strongly with infarct volume, NIHSS, and mRS (Figure 5). The increased expression of miR-21-5p, miR-29c-3p, and miR-27a-3p and increased levels of Gln, GlycA, glucose, and S-LDL-PL correlated only with infarct volume (Figure 5).

3.6. Machine Learning Analysis

The clinical, radiological, transcriptomic, and metabolic markers that showed the strongest association were utilized in the final analyses for predicting stroke outcomes. Feature selection using statistical tests such as Chi-square or F-regression tests identified nine metabolites and seven miRNAs that could be used as features for predicting stroke outcomes (Supplementary Table S6). Upon further analysis using SHAP, we identified 11 interactive features with a higher impact on the model. These included clinical parameters such as age, gender, infarct volume, and NIHSS, miRNAs including miR-625-3p, miR-23a-3p, miR-486-5p, and miR-199a-3p, and metabolites like acetoacetate, glucose, and total fatty acids (total-FA) as detailed in Figure 6A,B. The combined list of features was used in tenfold cross-validation. The ANN was the most powerful classifier, with an accuracy of 0.81 and an AUC of 0.91. It was followed by LDA with an accuracy of 0.79 and an AUC of 0.79, LR with an accuracy of 0.79 and an AUC of 0.77, and KNN with an accuracy of 0.71 and an AUC of 0.73. NB and SVM had similar accuracies and AUCs, both with an accuracy of 0.71 and an AUC of 0.71 for predicting stroke outcomes (Figure 6C).

4. Discussion

The disruption of the BBB leads to rapid changes in circulatory fluids. We hypothesize that these changes in metabolites and small extracellular RNAs in the bloodstream are associated with the size of the stroke. Measuring these changes might provide additional predictive value when combined with clinical and imaging tools to determine clinical outcomes. Our study showed that the infarct volume was a stronger determinant of stroke outcomes (r = 0.41; p = 8.2 × 10−4) than the NIHSS (r = 0.32; p = 0.004) (Figure 2). A strong correlation between NIHSS and infarct volume (r = 0.82; p < 0.0001) reported by Furlanis et al. [34] earlier could not be confirmed in our study (r = 0.32; p = 0.004). NIHSS and volume have been shown to predict stroke outcomes in AIS patients treated with intravenous thrombolysis [35]. Intravenous (IV) tissue plasminogen activator (rt-PA), a clot-dissolving drug, is effective in treating AIS and improving survival outcomes; however, the use of IVtPA for the majority of cases remained restricted to 4.5 h after the onset [21]. A recent study by Popa et al. showed that the door-to-CT time was the only factor significantly correlated with the death outcomes in comparison to door-to-physician time, NIHSS at admission, or after 24 h [36]. EVT is a mechanical procedure to remove a thrombus from the blocked brain artery during the AIS. In 2017, the use of EVT can be extended to 24 h in selected patients with small ischemic cores (<40mL). However, in patients with large ischemic cores, the application of EVT procedures showed higher rates of vascular and hemorrhagic complications leading to disabilities in recent clinical trials [37,38,39]. Infarct volume has been shown to discriminate patient outcomes in mild strokes in earlier studies [40], but no data are available on large core strokes. Therefore, given the limitations of existing technologies and stringency in the therapeutic window, the discovery of actionable molecular markers and their integration into the existing radiological and clinical protocols can improve the accurate prognostication of outcomes in patients with large ischemic cores.
Of the 41 identified small RNAs in patients with large infarct volumes, miR-144-5p was the most downregulated considerably (log2FC = −4.55; p = 1.1 × 10−5) in large strokes, and miR 144-5p and miR-625-3p showed a significant negative association with clinical outcome (mRS). Studies on animal models have shown miR-144 to have a protective role against neuroinflammation and oxidative stress by targeting the pathway of miR-451-14-3-3ζ-FoxO3 in intracerebral hemorrhage. [41]. Additionally, miR-144 downregulates PTEN and upregulates Akt in mice, supporting its protective role in AIS [42]. Decreased miR-625 expression was observed in laryngeal squamous cell carcinoma tissues, breast cancer, and melanoma, but some studies reported conflicting results [43,44,45,46]. Elevated serum levels of miR-21 were attributed to a 6.2-fold increase in stroke risk in patients with AIS [47]. MiR-21 was increased in patients with subarachnoid hemorrhage [48] and cerebral ischemia, leading to increased ERK-stimulated MMP9 levels [49].
Several circulating metabolites were differentially regulated with respect to infarct volume in AIS. Patients with large infarct sizes showed significantly increased levels of glucose and ketone bodies. Hyperglycemia and ketone bodies were highly prevalent in AIS in our patients from Oklahoma and the UK biobank [16,25]. In this study, glucose and ketone bodies were associated with poor outcomes supported by other published studies [50,51,52]. Urine ketone bodies were also positively correlated with recurrent stroke [53]. Multivariate regression and correlation analysis identified that eight miRNAs and one metabolite correlated with all clinical and neurological factors, including infarct volume, stroke severity (NIHSS), and outcomes (Figure 5). Through integrative analysis using machine learning, our study identified and confirmed the top 11 features to predict stroke outcomes by machine learning (Figure 6A,B). ANN was the most powerful classifier for predicting poor (mRS > 3) and good (mRS = 0 to 3) outcomes with an accuracy of (0.81) and AUC of (0.91). Notably, miR-199a-3p and miR-23a-3p identified by machine learning were not among the top miRNAs found in statistical analysis. MiR-199a-5p promotes functional recovery after cerebral ischemia by inhibiting Cav-1 to facilitate neurogenesis; it has been shown to be downregulated in rats with cerebral ischemia [54,55]. Similarly, miR-23a-3p lowers oxidative stress in cerebral ischemia [56]. Fatty acids, such as omega-6, mono- and poly-unsaturated fatty acids, and saturated fatty acids, were significantly decreased in AIS and inversely correlated with the infarct volume, as we reported earlier [25].
Besides molecular, clinical, and radiological factors, our interactive machine learning models identified gender to be a significant determinant of stroke outcomes. Male patients with small and large cores exhibited good outcomes compared to women (Figure 6). Sex differences have been shown to determine the treatment and outcomes of AIS [57]. Even though women had better survival rates after AIS, they also were more disabled and experienced poor quality of life in a study reported on 19000 patients [58]. Increased age was also a strong predictor of poor outcomes in our study (Figure 6), which is an important unmodifiable factor for poor recovery and neurological worsening in AIS. Age-related risk is modulated by gender and other risk factors in our study, as supported by earlier studies [59]. Machine learning and deep learning algorithms can identify outliers within the small patient datasets that may influence disease diagnosis [60,61]. To avoid overfitting, we detected and rectified the outliers through data scaling and tenfold cross-validation using supervised methods. However, independent validation using an external dataset would still be warranted.

Study Strengths and Limitations

The strength of this study includes a well-characterized and deeply phenotyped LVO patient cohort and the availability of radiological images, and infarct volume measures, which allowed the comparison of molecular markers with stroke severity and the outcomes. The discovery of actionable molecules implicated in the development of ischemic core and patient outcomes will be essential to better understanding AIS pathophysiology and improving treatment strategies. Machine learning approaches have been successfully used to integrate multi-omics data to predict hepatic carcinomas [62] and diabetes and classify stroke etiology [63]. This is the first human study reporting a joint evaluation using clinical, radiological, molecular, and artificial intelligence tools for predicting the functional outcomes of LVO strokes. The weaknesses include a relatively small sample size of LVO patients. We analyzed RNA sequencing data for miRNAs, which account for only 2% of the entire EV transcriptome. In addition, our study examined targeted metabolite profiles for 250 analytes using NMR. The use of global metabolomics (>3000 analytes) and their absolute quantifications would improve the detection of critical chemical features to predict stroke outcomes in large strokes. Even though the MISS patients were recruited from multiethnic populations, because of Oklahoma demographics, there is an over-representation of European Americans in this study, which limits the generalizability. Replication and validation of these results in other ethnic datasets would improve generalizability.

5. Conclusions

In summary, no currently existing prediction models can guide stroke prognostication in LVOs with large core strokes. This study lays the groundwork for further research to incorporate new molecular characteristics into the current determinants, aiming to enhance their predictive accuracy. Ultimately, our work provides a framework for advancing stroke therapeutics and care by adding molecular markers to improve prognostication based on clinical assessments, radiological exams, and artificial intelligence to enhance patient care and functional outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm13195917/s1: Figure S1: Scatter plots depicting the correlation of selected serum metabolites associated with infarct volume; Figure S2: Heat maps depicting the correlation of significant miRNAs with serum metabolites; Table S1: Association of serum exosome miRNAs with stroke infarct volume; Table S2: Correlation of serum exosome miRNAs with NIHSS; Table S3: Correlation of serum exosome miRNAs with stroke outcome; Table S4: Correlation of serum metabolites with NIHSS; Table S5: Correlation of serum metabolites with stroke outcome; Table S6: Feature selection of metabolites for machine learning based on stroke outcome.

Author Contributions

Methodology, E.V.S. and D.K.S.; Formal analysis, M.R.; Resources, A.V., E.V.S. and D.K.S.; Writing—original draft, M.R.; Writing—review and editing, D.K.S.; Supervision, E.V.S. and D.K.S.; Project administration, D.K.S.; Funding acquisition, D.K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was co-funded by grants from the College of Medicine Alumni Association (COMAA), the Presbyterian Health Foundation, the Leinbach Foundation, and Dr. Geoffrey Altshuler Endowment funds from the Children’s Health Foundation of the University of Oklahoma Health Sciences Center. The National Institute of Health (NIH/NIDDK) R01DK118427 grant also partly supported this study.

Institutional Review Board Statement

All study protocols were approved by the University of Oklahoma Health Sciences Center’s Institutional Review Board (IRB# 5349; Approval date 28 July 2024).

Informed Consent Statement

All participants provided written informed consent.

Data Availability Statement

The data used and/or analyzed during the current study will be available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank all the participants of MISS who made this study possible.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Feigin, V.L.; Stark, B.A.; Johnson, C.O.; Roth, G.A.; Bisignano, C.; Abady, G.G.; Abbasifard, M.; Abbasi-Kangevari, M.; Abd-Allah, F.; Abedi, V.; et al. Global, regional, and national burden of stroke and its risk factors, 1990-2019: A systematic analysis for the Global Burden of Disease Study 2019. Lancet Neurol. 2021, 20, 795–820. [Google Scholar] [CrossRef] [PubMed]
  2. Qin, H.; Wang, P.; Zhang, R.; Yu, M.; Zhang, G.; Liu, G.; Wang, Y. Stroke History is an Independent Risk Factor for Poor Prognosis in Ischemic Stroke Patients: Results from a Large Nationwide Stroke Registry. Curr. Neurovasc. Res. 2020, 17, 487–494. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, Y.H.; Lei, L.; Bao, Y.P.; Zhang, L. An Integrated Metabolomic Screening Platform Discovers the Potential Biomarkers of Ischemic Stroke and Reveals the Protective Effect and Mechanism of Folic Acid. Front. Mol. Biosci. 2022, 9, 783793. [Google Scholar] [CrossRef] [PubMed]
  4. Virani, S.S.; Alonso, A.; Benjamin, E.J.; Bittencourt, M.S.; Callaway, C.W.; Carson, A.P.; Chamberlain, A.M.; Chang, A.R.; Cheng, S.; Delling, F.N.; et al. Heart Disease and Stroke Statistics-2020 Update: A Report From the American Heart Association. Circulation 2020, 141, e139–e596. [Google Scholar] [CrossRef]
  5. Turc, G.; Maier, B.; Naggara, O.; Seners, P.; Isabel, C.; Tisserand, M.; Raynouard, I.; Edjlali, M.; Calvet, D.; Baron, J.C.; et al. Clinical Scales Do Not Reliably Identify Acute Ischemic Stroke Patients With Large-Artery Occlusion. Stroke 2016, 47, 1466–1472. [Google Scholar] [CrossRef]
  6. Rai, A.T.; Seldon, A.E.; Boo, S.; Link, P.S.; Domico, J.R.; Tarabishy, A.R.; Lucke-Wold, N.; Carpenter, J.S. A population-based incidence of acute large vessel occlusions and thrombectomy eligible patients indicates significant potential for growth of endovascular stroke therapy in the USA. J. Neurointerv. Surg. 2017, 9, 722–726. [Google Scholar] [CrossRef]
  7. Hiraga, A.; Yamaoka, T.; Sakai, Y.; Osakabe, Y.; Suzuki, A.; Hirose, N. Relationship between outcome in acute stroke patients and multiple stroke related scores obtained after onset of stroke. J. Phys. Ther. Sci. 2018, 30, 1310–1314. [Google Scholar] [CrossRef]
  8. Menezes, N.M.; Ay, H.; Wang Zhu, M.; Lopez, C.J.; Singhal, A.B.; Karonen, J.O.; Aronen, H.J.; Liu, Y.; Nuutinen, J.; Koroshetz, W.J.; et al. The real estate factor: Quantifying the impact of infarct location on stroke severity. Stroke 2007, 38, 194–197. [Google Scholar] [CrossRef]
  9. Rottiers, V.; Naar, A.M. MicroRNAs in metabolism and metabolic disorders. Nat. Rev. Mol. Cell Biol. 2012, 13, 239–250. [Google Scholar] [CrossRef]
  10. Vasudeva, K.; Munshi, A. miRNA dysregulation in ischaemic stroke: Focus on diagnosis, prognosis, therapeutic and protective biomarkers. Eur. J. Neurosci. 2020, 52, 3610–3627. [Google Scholar] [CrossRef]
  11. Gu, R.F.; Fang, T.; Nelson, A.; Gyoneva, S.; Gao, B.; Hedde, J.; Henry, K.; Peterson, E.; Burkly, L.C.; Wei, R. Proteomic Characterization of the Dynamics of Ischemic Stroke in Mice. J. Proteome Res. 2021, 20, 3689–3700. [Google Scholar] [CrossRef] [PubMed]
  12. Ruan, J.; Yao, Y. Behavioral tests in rodent models of stroke. Brain Hemorrhages 2020, 1, 171–184. [Google Scholar] [CrossRef] [PubMed]
  13. Sidorov, E.; Bejar, C.; Xu, C.; Ray, B.; Reddivari, L.; Chainakul, J.; Vanamala, J.K.P.; Sanghera, D.K. Potential Metabolite Biomarkers for Acute Versus Chronic Stage of Ischemic Stroke: A Pilot Study. J. Stroke Cerebrovasc. Dis. 2020, 29, 104618. [Google Scholar] [CrossRef] [PubMed]
  14. Sidorov, E.; Sanghera, D.K.; Vanamala, J.K.P. Biomarker for Ischemic Stroke Using Metabolome: A Clinician Perspective. J. Stroke 2019, 21, 31–41. [Google Scholar] [CrossRef]
  15. Sidorov, E.V.; Bejar, C.; Xu, C.; Ray, B.; Gordon, D.; Chainakul, J.; Sanghera, D.K. Novel Metabolites as Potential Indicators of Ischemic Infarction Volume: A Pilot Study. Transl. Stroke Res. 2021, 12, 778–784. [Google Scholar] [CrossRef]
  16. Sidorov, E.V.; Xu, C.; Garcia-Ramiu, J.; Blair, A.; Ortiz-Garcia, J.; Gordon, D.; Chainakul, J.; Sanghera, D.K. Global Metabolomic Profiling Reveals Disrupted Lipid and Amino Acid Metabolism Between the Acute and Chronic Stages of Ischemic Stroke. J. Stroke Cerebrovasc. Dis. 2022, 31, 106320. [Google Scholar] [CrossRef]
  17. LaCombe, D.M.; Gordon, D.L.; Issenberg, S.B.; Vega, A.; Brocato, C.; Siegel, M.; Reynolds, P.S. Stroke on the mend. JEMS 2000, 25, 32–41. [Google Scholar]
  18. Spilker, J.; Kongable, G.; Barch, C.; Braimah, J.; Brattina, P.; Daley, S.; Donnarumma, R.; Rapp, K.; Sailor, S. Using the NIH Stroke Scale to assess stroke patients. The NINDS rt-PA Stroke Study Group. J. Neurosci. Nurs. 1997, 29, 384–392. [Google Scholar] [CrossRef]
  19. del Zoppo, G.J.; Poeck, K.; Pessin, M.S.; Wolpert, S.M.; Furlan, A.J.; Ferbert, A.; Alberts, M.J.; Zivin, J.A.; Wechsler, L.; Busse, O.; et al. Recombinant tissue plasminogen activator in acute thrombotic and embolic stroke. Ann. Neurol. 1992, 32, 78–86. [Google Scholar] [CrossRef]
  20. Banks, J.L.; Marotta, C.A. Outcomes validity and reliability of the modified Rankin scale: Implications for stroke clinical trials: A literature review and synthesis. Stroke 2007, 38, 1091–1096. [Google Scholar] [CrossRef]
  21. Powers, W.J.; Rabinstein, A.A.; Ackerson, T.; Adeoye, O.M.; Bambakidis, N.C.; Becker, K.; Biller, J.; Brown, M.; Demaerschalk, B.M.; Hoh, B.; et al. Guidelines for the Early Management of Patients With Acute Ischemic Stroke: 2019 Update to the 2018 Guidelines for the Early Management of Acute Ischemic Stroke: A Guideline for Healthcare Professionals From the American Heart Association/American Stroke Association. Stroke A J. Cereb. Circ. 2019, 50, e344–e418. [Google Scholar] [CrossRef]
  22. Nunn, A.; Bath, P.M.; Gray, L.J. Analysis of the Modified Rankin Scale in Randomised Controlled Trials of Acute Ischaemic Stroke: A Systematic Review. Stroke Res. Treat. 2016, 2016, 9482876. [Google Scholar] [CrossRef] [PubMed]
  23. Saver, J.L.; Filip, B.; Hamilton, S.; Yanes, A.; Craig, S.; Cho, M.; Conwit, R.; Starkman, S.; Investigators, F.-M. Coordinators Improving the reliability of stroke disability grading in clinical trials and clinical practice: The Rankin Focused Assessment (RFA). Stroke 2010, 41, 992–995. [Google Scholar] [CrossRef] [PubMed]
  24. Onderwater, G.L.J.; Ligthart, L.; Bot, M.; Demirkan, A.; Fu, J.; van der Kallen, C.J.H.; Vijfhuizen, L.S.; Pool, R.; Liu, J.; Vanmolkot, F.H.M.; et al. Large-scale plasma metabolome analysis reveals alterations in HDL metabolism in migraine. Neurology 2019, 92, e1899–e1911. [Google Scholar] [CrossRef] [PubMed]
  25. Rout, M.; Vaughan, A.; Blair, A.; Stavrakis, S.; Sidorov, E.V.; Sanghera, D.K. Discovery and validation of circulating stroke metabolites by NMR-based analyses using patients from the MISS and UK Biobank. Neurochem. Int. 2023, 169, 105588. [Google Scholar] [CrossRef]
  26. Rout, M.; Malone-Perez, M.W.; Park, G.; Lerner, M.; Kimble Frazer, J.; Apple, B.; Vaughn, A.; Payton, M.; Stavrakis, S.; Sidorov, E.; et al. Contribution of circulating Mfge8 to human T2DM and cardiovascular disease. Gene 2024, 927, 148712. [Google Scholar] [CrossRef]
  27. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, R25. [Google Scholar] [CrossRef]
  28. Sanghera, D.K.; Hopkins, R.; Malone-Perez, M.W.; Bejar, C.; Tan, C.; Mussa, H.; Whitby, P.; Fowler, B.; Rao, C.V.; Fung, K.A.; et al. Targeted sequencing of candidate genes of dyslipidemia in Punjabi Sikhs: Population-specific rare variants in GCKR promote ectopic fat deposition. PLoS ONE 2019, 14, e0211661. [Google Scholar] [CrossRef]
  29. Huber, W.; Carey, V.J.; Gentleman, R.; Anders, S.; Carlson, M.; Carvalho, B.S.; Bravo, H.C.; Davis, S.; Gatto, L.; Girke, T.; et al. Orchestrating high-throughput genomic analysis with Bioconductor. Nat. Methods 2015, 12, 115–121. [Google Scholar] [CrossRef]
  30. Ambros, V.; Bartel, B.; Bartel, D.P.; Burge, C.B.; Carrington, J.C.; Chen, X.; Dreyfuss, G.; Eddy, S.R.; Griffiths-Jones, S.; Marshall, M.; et al. A uniform system for microRNA annotation. RNA 2003, 9, 277–279. [Google Scholar] [CrossRef]
  31. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  32. Chen, T.G.C. XGBoost: A scalable tree boosting system. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; ACM: New York, NY, USA, 2016; pp. 785–794. [Google Scholar]
  33. Lundberg, S.M.; Erion, G.G.; Lee, S.I. Consistent individualized feature attribution for tree ensembles. arXiv 2018, arXiv:1802.03888. [Google Scholar]
  34. Furlanis, G.; Ajcevic, M.; Stragapede, L.; Lugnan, C.; Ridolfi, M.; Caruso, P.; Naccarato, M.; Ukmar, M.; Manganotti, P. Ischemic Volume and Neurological Deficit: Correlation of Computed Tomography Perfusion with the National Institutes of Health Stroke Scale Score in Acute Ischemic Stroke. J. Stroke Cerebrovasc. Dis. 2018, 27, 2200–2207. [Google Scholar] [CrossRef] [PubMed]
  35. Sanak, D.; Herzig, R.; Zapletalova, J.; Horak, D.; Kral, M.; Skoloudik, D.; Bartkova, A.; Veverka, T.; Herman, M.; Kanovsky, P. Predictors of good clinical outcome in acute stroke patients treated with intravenous thrombolysis. Acta Neurol. Scand. 2011, 123, 339–344. [Google Scholar] [CrossRef]
  36. Popa, D.; Iancu, A.; Petrica, A.; Buleu, F.; Williams, C.G.; Sutoi, D.; Trebuian, C.; Tudor, A.; Mederle, O.A. Emergency Department Time Targets for Interhospital Transfer of Patients with Acute Ischemic Stroke. J. Pers. Med. 2023, 14, 13. [Google Scholar] [CrossRef]
  37. Yoshimura, S.; Sakai, N.; Yamagami, H.; Uchida, K.; Beppu, M.; Toyoda, K.; Matsumaru, Y.; Matsumoto, Y.; Kimura, K.; Takeuchi, M.; et al. Endovascular Therapy for Acute Stroke with a Large Ischemic Region. N. Engl. J. Med. 2022, 386, 1303–1313. [Google Scholar] [CrossRef]
  38. Huo, X.; Ma, G.; Tong, X.; Zhang, X.; Pan, Y.; Nguyen, T.N.; Yuan, G.; Han, H.; Chen, W.; Wei, M.; et al. Trial of Endovascular Therapy for Acute Ischemic Stroke with Large Infarct. N. Engl. J. Med. 2023, 388, 1272–1283. [Google Scholar] [CrossRef]
  39. Sarraj, A.; Hassan, A.E.; Abraham, M.G.; Ortega-Gutierrez, S.; Kasner, S.E.; Hussain, M.S.; Chen, M.; Blackburn, S.; Sitton, C.W.; Churilov, L.; et al. Trial of Endovascular Thrombectomy for Large Ischemic Strokes. N. Engl. J. Med. 2023, 388, 1259–1271. [Google Scholar] [CrossRef]
  40. Vagal, A.S.; Sucharew, H.; Prabhakaran, S.; Khatri, P.; Jovin, T.; Michel, P.; Wintermark, M. Final infarct volume discriminates outcome in mild strokes. Neuroradiol. J. 2015, 28, 404–408. [Google Scholar] [CrossRef]
  41. Wang, X.; Hong, Y.; Wu, L.; Duan, X.; Hu, Y.; Sun, Y.; Wei, Y.; Dong, Z.; Wu, C.; Yu, D.; et al. Deletion of MicroRNA-144/451 Cluster Aggravated Brain Injury in Intracerebral Hemorrhage Mice by Targeting 14-3-3zeta. Front Neurol. 2020, 11, 551411. [Google Scholar] [CrossRef]
  42. Zhong, S.J.; Cui, M.M.; Gao, Y.T.; Cao, X.Y.; Chen, B.; Wen, X.R. MicroRNA-144 promotes remote limb ischemic preconditioning-mediated neuroprotection against ischemic stroke via PTEN/Akt pathway. Acta Neurol. Belg. 2021, 121, 95–106. [Google Scholar] [CrossRef] [PubMed]
  43. Zhou, W.B.; Zhong, C.N.; Luo, X.P.; Zhang, Y.Y.; Zhang, G.Y.; Zhou, D.X.; Liu, L.P. miR-625 suppresses cell proliferation and migration by targeting HMGA1 in breast cancer. Biochem. Biophys. Res. Commun. 2016, 470, 838–844. [Google Scholar] [CrossRef] [PubMed]
  44. Fang, W.; Fan, Y.; Fa, Z.; Xu, J.; Yu, H.; Li, P.; Gu, J. microRNA-625 inhibits tumorigenicity by suppressing proliferation, migration and invasion in malignant melanoma. Oncotarget 2017, 8, 13253–13263. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, J.; Zhang, J.; Zhang, J.; Qiu, W.; Xu, S.; Yu, Q.; Liu, C.; Wang, Y.; Lu, A.; Zhang, J.; et al. MicroRNA-625 inhibits the proliferation and increases the chemosensitivity of glioma by directly targeting AKT2. Am. J. Cancer Res. 2017, 7, 1835–1849. [Google Scholar]
  46. Zhang, M.; Xiong, F.; Zhang, S.; Guo, W.; He, Y. Crucial Roles of miR-625 in Human Cancer. Front. Med. 2022, 9, 845094. [Google Scholar] [CrossRef]
  47. Tsai, P.C.; Liao, Y.C.; Wang, Y.S.; Lin, H.F.; Lin, R.T.; Juo, S.H. Serum microRNA-21 and microRNA-221 as potential biomarkers for cerebrovascular disease. J. Vasc. Res. 2013, 50, 346–354. [Google Scholar] [CrossRef]
  48. Bache, S.; Rasmussen, R.; Rossing, M.; Laigaard, F.P.; Nielsen, F.C.; Moller, K. MicroRNA Changes in Cerebrospinal Fluid After Subarachnoid Hemorrhage. Stroke 2017, 48, 2391–2398. [Google Scholar] [CrossRef]
  49. Deng, X.; Zhong, Y.; Gu, L.; Shen, W.; Guo, J. MiR-21 involve in ERK-mediated upregulation of MMP9 in the rat hippocampus following cerebral ischemia. Brain Res. Bull. 2013, 94, 56–62. [Google Scholar] [CrossRef]
  50. Chamorro, A.; Brown, S.; Amaro, S.; Hill, M.D.; Muir, K.W.; Dippel, D.W.J.; van Zwam, W.; Butcher, K.; Ford, G.A.; den Hertog, H.M.; et al. Glucose Modifies the Effect of Endovascular Thrombectomy in Patients With Acute Stroke. Stroke 2019, 50, 690–696. [Google Scholar] [CrossRef]
  51. Osei, E.; Fonville, S.; Zandbergen, A.A.M.; Koudstaal, P.J.; Dippel, D.W.J.; den Hertog, H.M. Impaired fasting glucose is associated with unfavorable outcome in ischemic stroke patients treated with intravenous alteplase. J. Neurol. 2018, 265, 1426–1431. [Google Scholar] [CrossRef]
  52. Desilles, J.P.; Meseguer, E.; Labreuche, J.; Lapergue, B.; Sirimarco, G.; Gonzalez-Valcarcel, J.; Lavallee, P.; Cabrejo, L.; Guidoux, C.; Klein, I.; et al. Diabetes mellitus, admission glucose, and outcomes after stroke thrombolysis: A registry and systematic review. Stroke 2013, 44, 1915–1923. [Google Scholar] [CrossRef] [PubMed]
  53. Dai, L.; Wang, A.; Gu, H.; Zhang, Y.; Zuo, Y.; Meng, X.; Chen, P.; Tian, X.; Li, H.; Wang, Y. Urinary ketone bodies and stroke recurrence in patient with acute ischemic stroke or TIA. J. Clin. Neurosci. 2023, 117, 79–83. [Google Scholar] [CrossRef] [PubMed]
  54. Jin, H.Q.; Jiang, W.F.; Zheng, X.T.; Li, L.; Fang, Y.; Yang, Y.; Hu, X.W.; Chu, L.S. MiR-199a-5p enhances neuronal differentiation of neural stem cells and promotes neurogenesis by targeting Cav-1 after cerebral ischemia. CNS Neurosci. Ther. 2023, 29, 3967–3979. [Google Scholar] [CrossRef] [PubMed]
  55. Li, M.; Luan, L.; Liu, Q.; Liu, Y.; Lan, X.; Li, Z.; Liu, W. MiRNA-199a-5p Protects Against Cerebral Ischemic Injury by Down-Regulating DDR1 in Rats. World Neurosurg. 2019, 131, e486–e494. [Google Scholar] [CrossRef]
  56. Zhao, H.; Tao, Z.; Wang, R.; Liu, P.; Yan, F.; Li, J.; Zhang, C.; Ji, X.; Luo, Y. MicroRNA-23a-3p attenuates oxidative stress injury in a mouse model of focal cerebral ischemia-reperfusion. Brain Res. 2014, 1592, 65–72. [Google Scholar] [CrossRef]
  57. Ouyang, M.; Shajahan, S.; Liu, X.; Sun, L.; Carcel, C.; Harris, K.; Anderson, C.S.; Woodward, M.; Wang, X. Sex differences in the utilization and outcomes of endovascular treatment after acute ischemic stroke: A systematic review and meta-analysis. Front. Glob. Womens Health 2022, 3, 1032592. [Google Scholar] [CrossRef]
  58. Carcel, C.; Wang, X.; Sandset, E.C.; Delcourt, C.; Arima, H.; Lindley, R.; Hackett, M.L.; Lavados, P.; Robinson, T.G.; Munoz Venturelli, P.; et al. Sex differences in treatment and outcome after stroke: Pooled analysis including 19,000 participants. Neurology 2019, 93, e2170–e2180. [Google Scholar] [CrossRef]
  59. Ohya, Y.; Matsuo, R.; Sato, N.; Irie, F.; Wakisaka, Y.; Ago, T.; Kamouchi, M.; Kitazono, T.; Investigators for Fukuoka Stroke, R. Modification of the effects of age on clinical outcomes through management of lifestyle-related factors in patients with acute ischemic stroke. J. Neurol. Sci. 2023, 446, 120589. [Google Scholar] [CrossRef]
  60. Kokol, P.; Kokol, M.; Zagoranski, S. Machine learning on small size samples: A synthetic knowledge synthesis. Sci. Prog. 2022, 105, 368504211029777. [Google Scholar] [CrossRef]
  61. Terrades, O.R.; Berenguel, A.; Gil, D. A Flexible Outlier Detector Based on a Topology Given by Graph Communities. Big Data Res. 2022, 29, 100332. [Google Scholar] [CrossRef]
  62. Chaudhary, K.; Poirion, O.B.; Lu, L.; Garmire, L.X. Deep Learning-Based Multi-Omics Integration Robustly Predicts Survival in Liver Cancer. Clin. Cancer Res. 2018, 24, 1248–1259. [Google Scholar] [CrossRef]
  63. Labarga, A.; Martinez-Gonzalez, J.; Barajas, M. Integrative Multi-Omics Analysis for Etiology Classification and Biomarker Discovery in Stroke: Advancing towards Precision Medicine. Biology 2024, 13, 338. [Google Scholar] [CrossRef]
Jcm 13 05917 g001
Figure 1. Flow chart of the study.
Figure 1. Flow chart of the study.
Jcm 13 05917 g001
Jcm 13 05917 g002
Figure 2. Correlation analysis in the scatter plots shows correlations between (A) NIHSS with infarct volume and (B) mRS with infarct volume in the log-transformed data values of NIHSS and mRS. NIHSS: National Institute of Health Stroke Scale; mRS: modified Rankin Score.
Figure 2. Correlation analysis in the scatter plots shows correlations between (A) NIHSS with infarct volume and (B) mRS with infarct volume in the log-transformed data values of NIHSS and mRS. NIHSS: National Institute of Health Stroke Scale; mRS: modified Rankin Score.
Jcm 13 05917 g002
Jcm 13 05917 g003
Figure 3. Scatter plots of correlation analysis depicting the miRNAs (miR-144-5p, miR-486-5p, miR-96-5p and let-7f-5p) whose expression was downregulated in patients with large infarcts and miRNAs (miR-101-3p and miR-21-5p), which were upregulated with increased infarct volume.
Figure 3. Scatter plots of correlation analysis depicting the miRNAs (miR-144-5p, miR-486-5p, miR-96-5p and let-7f-5p) whose expression was downregulated in patients with large infarcts and miRNAs (miR-101-3p and miR-21-5p), which were upregulated with increased infarct volume.
Jcm 13 05917 g003
Jcm 13 05917 g004
Figure 4. Heat maps depicting the differential regulation of significant miRNAs in large and small infarct volumes.
Figure 4. Heat maps depicting the differential regulation of significant miRNAs in large and small infarct volumes.
Jcm 13 05917 g004
Jcm 13 05917 g005
Figure 5. Venn diagrams showing the most relevant (A) miRNAs and (B) metabolites associated with infarct volume, NIHSS, and mRS. NIHSS: National Institute of Health Stroke Scale; mRS: modified Rankin Score; Gln: glutamine; GlycA: glycoprotein acetyls; S-LDL-PL: phospholipids in small LDL; L-HDL-FC%: free cholesterol to total lipids ratio in large HDL; S-LDL-C%: cholesterol to total lipids ratio in small LDL; M-LDL-C%: cholesterol to total lipids ratio in medium LDL; IDL-TG%: triglycerides to total lipids ratio in IDL; L-LDL-TG%: triglycerides to total lipids ratio in large LDL; L-HDL-TG%: triglycerides to total lipids ratio in large HDL; L-HDL-C: cholesterol in large HDL; L-HDL-PL: phospholipids in large HDL; L-HDL-CE: cholesteryl esters in large HDL; L-HDL-L: total lipids in large HDL: XS-VLDL-TG%: triglycerides to total lipids ratio in very small VLDL; bOHbutyrate: 3-Hydroxybutyrate. Upward arrow- increased metabolite level/miRNA expression; Downward arrow- decreased metabolite level or miRNA expression.
Figure 5. Venn diagrams showing the most relevant (A) miRNAs and (B) metabolites associated with infarct volume, NIHSS, and mRS. NIHSS: National Institute of Health Stroke Scale; mRS: modified Rankin Score; Gln: glutamine; GlycA: glycoprotein acetyls; S-LDL-PL: phospholipids in small LDL; L-HDL-FC%: free cholesterol to total lipids ratio in large HDL; S-LDL-C%: cholesterol to total lipids ratio in small LDL; M-LDL-C%: cholesterol to total lipids ratio in medium LDL; IDL-TG%: triglycerides to total lipids ratio in IDL; L-LDL-TG%: triglycerides to total lipids ratio in large LDL; L-HDL-TG%: triglycerides to total lipids ratio in large HDL; L-HDL-C: cholesterol in large HDL; L-HDL-PL: phospholipids in large HDL; L-HDL-CE: cholesteryl esters in large HDL; L-HDL-L: total lipids in large HDL: XS-VLDL-TG%: triglycerides to total lipids ratio in very small VLDL; bOHbutyrate: 3-Hydroxybutyrate. Upward arrow- increased metabolite level/miRNA expression; Downward arrow- decreased metabolite level or miRNA expression.
Jcm 13 05917 g005
Jcm 13 05917 g006
Figure 6. (A) The most important predictive parameters using binary classification with mRS, with the absolute value of a feature being high (red) or low (blue), depicting bad or good stroke outcomes. (B) The VIP plot shows the most influential biomarkers and their impact on stroke outcomes. (C) Machine learning analysis for clinical, miRNAs, and metabolite features using binary classification and artificial neural networks was the most informative model showing high sensitivity and specificity in AUC. NIHSS: National Institute of Health Stroke Scale; FA: fatty acids; mRS: modified Rankin Scale.
Figure 6. (A) The most important predictive parameters using binary classification with mRS, with the absolute value of a feature being high (red) or low (blue), depicting bad or good stroke outcomes. (B) The VIP plot shows the most influential biomarkers and their impact on stroke outcomes. (C) Machine learning analysis for clinical, miRNAs, and metabolite features using binary classification and artificial neural networks was the most informative model showing high sensitivity and specificity in AUC. NIHSS: National Institute of Health Stroke Scale; FA: fatty acids; mRS: modified Rankin Scale.
Jcm 13 05917 g006
Table 1. Clinical characteristics of study participants.
Table 1. Clinical characteristics of study participants.
TRAITLarge Infarct Volume
(N = 43)		Small Infarct Volume
(N = 42)		p-Value
Age (years)64.9 ± 12.862.2 ± 16.40.40
BMI (kg/m2)29.5 ± 7.029.6 ± 8.90.94
Glucose (mg/dL)162.9 ± 96.7128.5 ± 59.30.05
Gender/Male (%)51400.33
Ethnicity (%)Caucasians9374-
African American714-
Others-12-
Systolic BP (mmHg)158.1 ± 32.7145.7 ± 21.90.04
Diastolic BP (mmHg)83.7 ± 19.087.0 ± 15.10.40
Type 2 diabetes (%)53360.12
Smoking (%)44500.66
Stroke outcome4.2 ± 1.92.6 ± 1.90.0004
NIHSS12.6 ± 9.49.5 ± 8.00.11
Infarct Volume117.4 ± 69.119.1 ± 14.17.43 × 10−12
Values are displayed as mean±SD; BMI: body mass index; BP: blood pressure; NIHSS: National Institute of Health Stroke Scale.
Table 2. Association of serum metabolites with infarct volume.
Table 2. Association of serum metabolites with infarct volume.
MetabolitesBiomarker Namet.statp-Valuelog2(FC)BetaSEp-Value
GlucoseGlucose2.242.84 ×10−20.890.510.173.23 ×10−3
S-HDL-FC%Free cholesterol to total lipids ratio in small HDL2.054.47 ×10−20.930.530.184.61 ×10−3
S-LDL-PLPhospholipids in small LDL2.679.68 ×10−30.980.340.197.74 ×10−2
GlnGlutamine2.163.43 ×10−21.970.630.379.14 ×10−2
AcetoacetateAcetoacetate2.094.09 ×10−21.180.360.162.97 ×10−2
GlycAGlycoprotein acetyls1.441.54 ×10−10.720.480.181.15 ×10−2
XL-HDL-PConcentration of very large HDL particles1.411.63 ×10−10.460.640.237.24 ×10−3
XL-HDL-CECholesteryl esters in very large HDL1.381.72 ×10−10.480.730.266.96 ×10−3
XL-HDL-CCholesterol in very large HDL1.242.20 ×10−10.480.730.281.12 ×10−2
XL-HDL-LTotal lipids in very large HDL1.102.77 ×10−10.440.690.291.97 ×10−2
XL-HDL-FC%Free cholesterol to total lipids ratio in very large HDL−1.292.02 ×10−1−0.70−0.750.291.21 ×10−2
SE: standard error; FC: fold change; data were analyzed using Metaboanalyst 5.0 and multivariate linear regression analysis adjusting for covariates of age, gender, and BMI.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rout, M.; Vaughan, A.; Sidorov, E.V.; Sanghera, D.K. Improving Stroke Outcome Prediction Using Molecular and Machine Learning Approaches in Large Vessel Occlusion. J. Clin. Med. 2024, 13, 5917. https://doi.org/10.3390/jcm13195917

AMA Style

Rout M, Vaughan A, Sidorov EV, Sanghera DK. Improving Stroke Outcome Prediction Using Molecular and Machine Learning Approaches in Large Vessel Occlusion. Journal of Clinical Medicine. 2024; 13(19):5917. https://doi.org/10.3390/jcm13195917

Chicago/Turabian Style

Rout, Madhusmita, April Vaughan, Evgeny V. Sidorov, and Dharambir K. Sanghera. 2024. "Improving Stroke Outcome Prediction Using Molecular and Machine Learning Approaches in Large Vessel Occlusion" Journal of Clinical Medicine 13, no. 19: 5917. https://doi.org/10.3390/jcm13195917

APA Style

Rout, M., Vaughan, A., Sidorov, E. V., & Sanghera, D. K. (2024). Improving Stroke Outcome Prediction Using Molecular and Machine Learning Approaches in Large Vessel Occlusion. Journal of Clinical Medicine, 13(19), 5917. https://doi.org/10.3390/jcm13195917

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop