2. Materials and Methods
2.1. Study subjects
Patients enrolled in the study were urgently admitted to University Hospital Centre “Zagreb“ Neurology Department and University Hospital Centre “Sestre milosrdnice” Radiology Department in Zagreb, Croatia, receiving mechanical thrombectomy (MT) with or without bridging intravenous thrombolytic therapy (IVT) between. They had been selected based on our University Centre stroke guidelines which are mainly based on European Stroke Organization (ESO) mechanical thrombectomy guidelines published in 2019 [
19] and ESO thrombolysis guidelines published in 2021 [
20] and ESO thrombolysis recommendation before MT published in 2022 [
21].
All included patients (N = 14) were adults (> 18 years), had a focal neurologic deficit defined as at least 4 points on National Institute of Health Stroke Scale (NIHSS) and duration of symptoms up to 24 hours. The non-contrast computed tomography (CT) was used to exclude cerebral haemorrhage or other causative intracranial pathology other than ischaemia. CT angiography confirmed intracranial large vessel occlusion of the anterior or posterior circulation (internal carotid artery [ICA], carotid T, M1, M2 segments of the middle cerebral artery [MCA], A1 segment of the anterior cerebral artery [ACA], P1 segment of the posterior cerebral artery [PCA], basilar artery [BA], vertebral artery [VA]) and an Alberta Stroke Program Early CT Score (ASPECT) between 6 and 10. Patients that arrived later than 6 hours from symptom onset and underwent MT had CT perfusion imaging demonstrating core/penumbra mismatch defined as core volume less than 70 mL, critically hypoperfused volume/core ratio larger than 1, 2 and mismatch volume larger than 10 mL. Patients that received bridging IVT had symptom duration up to 4 - 5 hours or perfusion mismatch on CT perfusion imaging according to the aforementioned criteria. All the patients that received IVT were administered with alteplase (0,9 mg/kg).
Types of strokes were defined as cardioembolism, large-artery atherosclerosis or other causes based on the Trial of ORG 10172 in Acute Stroke Treatment (TOAST) classification [
22]. Mechanical thrombectomy was performed using all three main techniques: stent retriever technique, contact aspiration and combined technique.
Participants in the control group (N = 14) were selected from Department of Neurosurgery, University Hospital Centre “Sestre milosrdnice”, Zagreb, Croatia. The baseline assessments were conducted for all participants, including demographic information, medical history, and relevant clinical measures ensuring comparability with the stroke group as closely as possible.
In the control group, exclusion criteria included history of stroke and any other cerebrovascular events or neurological disorders (other than stroke), severe cardiovascular diseases, confirmed chronic inflammatory conditions, confirmed severe renal or hepatic diseases, uncontrolled hypertension as well as pregnancy and intake of medications that could interfere with the study outcomes.
2.2. Clinical and Biochemical Measurements
For clinical laboratory and targeted metabolomics analysis, venous blood samples were collected after 12 h fasting for control group and immediately before IVT and MT for patient group. 10 mL of venous blood was collected and centrifuged (750 × g for 15 min after 30 min at room temperature) to separate serum from cellular components. The serum aliquots were stored in glass tubes at -80°C until measurements were performed.
Biochemical parameters were measured at the Biochemistry Laboratory of University Hospital Centre “Zagreb” and University Hospital Centre “Sestre milosrdnice” in Zagreb, Croatia. Biochemical analyses of serums samples in clinical setting were involved measurement of blood glucose, lipid profile (total cholesterol [tCH], triglycerides [TG], HDL-cholesterol [HDL-CH], LDL-cholesterol [LDL-CH]) as well as coagulation profile (prothrombin time [PT], activated partial thromboplastin time [aPTT]), international normalized ratio [INR]) by standard laboratory procedures [
23,
24]. Weight (kg) and height (m) measurements were done in a standardized manner. Body weight (BW) was measured by standard medical balance, accurate to 0.1 kg. Body height (BH) was measured using the stadiometer with a precision of 0.1 cm. Body mass index (BMI) was calculated by standard formula BMI = BW (kg)/BH (m
2). Systolic and diastolic blood pressure were measured using sphygmomanometer.
2.3. Preparation of Stroke Thrombi for Targeted Metabolomics Analysis
After MT procedure, each thrombus sample was washed with cold saline solution and stored in glass tubes at -80°C until metabolomics analysis. Collected stroke thrombi were weighted on ice and 3-fold volume of ice-cold isopropanol was used as extraction solvent (
Table 1).
The samples were vortexed and sonicated on ice 3 times for 30 seconds at max amplitude. Samples were spin-down between cycles and held on ice for 1 min. After third sonication cycle, samples were vortexed for 30 sec and centrifuged at 10 000 x g for 5 min at 4°C. The supernatants were used for further analysis on AbsoluteIDQ p400 kit (Biocrates Life Science AG, Innsbruck, Austria) with mass spectrometry (Q Exactive Plus hybrid quadrupole - Orbitrap; Thermo Fisher Scientific, Bremen, Germany).
Serums and stroke thrombi samples were transported in dry ice to the Laboratory of Proteomics, Clinic for Internal Diseases, Faculty of Veterinary Medicine in Zagreb, Croatia, where targeted metabolomics analysis was performed.
2.4. Targeted Metabolomics Analysis
Concentration of 408 endogenous metabolites divided into 11 analytical classes, including amino acids (N = 21), biogenic amines (N = 21), monosaharides comprising the sum of hexoses, including glucose (N = 1), phosphatidylcholines (N = 172), lysophosphatidylcholines (N = 24), sphingomyelins (N = 31), ceramides (N = 9), diglycerides (N = 18), triglycerides (N = 42), cholesteryl esters (N = 14) and acylcarnitines (N = 55) were measured in patient serums, control serums and stroke thrombi samples with Biocrates AbsoluteIDQ p400 kit using mass spectrometry (MS) according to the manufacturer’s protocol. Metabolites were quantified using the combination of liquid chromatography-mass spectrometry (LC-MS/MS) and flow-injection analysis-mass spectrometry (FIA-MS/MS). The LC-MS/MS analysis allowed the quantification of amino acids and biogenic amines while the FIA-MS/MS was used to quantify acylcarnitines, glycerophospholipids, glycerides, hexoses, cholesteryl esters, and sphingolipids.
Sample preparation and metabolite extraction were performed according to the manufacturer instructions provided with the kit and described in detail previously [
25]. In short, the 10 μL of each sample, calibration standards, blank samples (phosphate-buffered saline; BDH PROLABO, Lutterworth, UK), and three quality control samples (low, QC1; medium, QC2 and high, QC3) were prepared on the specific 96-well plate for protein removal, internal standard normalisation and derivatization. The samples were extracted using 5 mM ammonium acetate (Sigma-Aldrich, St. Louis, MO, USA) in methanol (Honeywell, Charlotte, NC, USA). Finally, extracts were diluted with the LC-MS-grade water prior to LC-MS/MS analysis and with FIA mobile phase (made by mixing 290 mL MeOH and a 10 mL ampule Biocrates FIA mobile phase additive, provided with the kit) prior to FIA-MS/MS analysis.
Mass spectrometry analysis was performed on a Dionex Ultimate 3000 UHPLC system (Thermo Fisher Scientific, Germering, Germany) coupled to a Q Exactive Plus hybrid quadrupole - Orbitrap mass spectrometer using an electrospray ionization source. A Thermo p400 HR UHPLC column provided with the kit was used to separate metabolites from the samples. The column temperature was maintained at 50°C. Mobile phase A contained 0.2 % formic acid in H2O, and B contained 0.2 % formic acid in acetonitrile. The injection volume was 5 μL. The total run time was 5.81 min, and the gradient change of 0 to 95 % of mobile phase B over 4 min. The flow rate was 0.8 mL/min.
For the FIA-MS/MS analysis, metabolites were eluted using the FIA mobile phase at a flow rate of 0.05 mL/min for the first 1.6 min, then the flow rate increased to 0.2 mL/min for 1.2 min and then decreased back to 0.05 mL/min for the end of the sequence.
The mass spectrometer was operated in positive and negative ion modes for both LC-MS/MS and FIA-MS/MS according to the instructions from Biocrates kit.
The MetIDQ software package Version Boron, an integral part of the Absolute IDQ p400 kit, was used for the data processing, quality assessment, and data export. The quantification of the LC-MS metabolites was processed via XCalibur Quan 4.1 software (Thermo Fisher Scientific, Waltham, MA, USA) based on a seven-point calibration curve and isotope labelled internal standards for most analytes, while the FIA-MS/MS analysis used a single-point calibrator with the representative internal standards. All reagents used in this analysis were of LC-MS grade and purchased from Merck (Darmstadt, Germany), Honeywell (Charlotte, NC, USA) and Sigma-Aldrich (St. Louis, MO, USA).
Blank samples were used for the calculation of the limits of detection (LOD). In terms of quantification, if the compounds were quantified with restriction, then the calibration curves had expected coefficients of determination (R2) < 0.99 according to the manufacturer guidelines. When specific standards were not commercially available and verification of the accuracy was not possible by the manufacturer, the measuring was performed “semi-quantitatively”.
2.5. Statistical Analyses
The MedCalc 22.019 (Frank Schoonjans, Mariakerke, Belgium) was used for descriptive statistical analysis and for the comparison of clinical and biochemical data between the study groups using Student’s t-test or Mann-Whitney test, depending on distribution normality.
The MetaboAnalyst v.4.0 software (
http://www.metaboanalyst.ca, accessed on 15 November 2023) was used for univariate and multivariate statistical analysis of metabolites. Metabolites with 50% missing values were removed and replace by LoDs (1/5 of the minimum positive value of each variable). The concentrations of each metabolite were normalised by median, log transformed, and Pareto scaled. The overall differences in the metabolomics profile between control and patient serum samples and between patient serums and stroke thrombi were analysed by Student’s t-test, a partial least square–discriminant analysis (PLS-DA), the variable importance on projection (VIP) and hierarchical clustering analysis (HCA). For all statistical analysis a
p-value ≤ 0.05 was considered to be statistically significant.
4. Discussion
Ischemic stroke (IS) followed by cerebral ischemia causes the cascade of highly complex and interrelated pathophysiological processes such as energy failure and excitotoxicity, oxidative stress, blood-brain barrier dysfunction, microvascular injury, post-ischemic inflammation, and finally, death of neurons, glia, and endothelial cells [
26]. Due to the impaired integrity of the blood-brain barrier after cerebral ischemia, many of the post-ischemic metabolites are more easily found in the plasma or serum of the patient [
16]. Therefore, metabolomics represents a key tool for analysing end products of cellular metabolism in various pathologies [
8]. Metabolomic research of stroke thrombi in ischemic stroke is very limited. To our knowledge, only several studies have been performed as non-target metabolomic analysis and using solely stroke thrombi, being mainly focused on single group of metabolites such as glucose and sorbitol [
17] or glycerophospholipids in very small number of stroke thrombi samples [
18]. In this study, we have performed comparative targeted metabolomic analysis of 14 stroke thrombi, corresponding patient serums and healthy control serums and identified specific subset of stroke-related metabolites.
According to both multivariate and univariate analysis, the concentrations of some metabolites were statistically significant between patient and control serums, while some were statistically significant between patient serums and stroke thrombi.
Glutamate concentration was significantly increased in patient serums compared to control serums (
Figure 2 and
Figure 3), with the highest VIP score (VIP > 2.5,
Figure 1b), which is in accordance to reported studies. Glutamate is excitatory amino acid and neurotransmitter having important role in mediating neuronal damage during cerebral ischemia [
13]. It is considered the main contributor to the ischemic brain tissue excitotoxicity as a result of energy failure [
27] and higher concentrations of glutamate in the blood and cerebrospinal fluid have been reported and associated with poor clinical outcome and neurological impairment after stroke [
15,
28].
Contrary to glutamate,
serotonin concentration was significantly decreased in patient serums compared to control serums (
Figure 2 and
Figure 3), also having high VIP score (VIP > 2,
Figure 1b). This is in accordance to some early reports indicating decreased plasma serotonin levels in patients with acute ischemic stroke [
29]. On the other hand, some studies reported elevated post-stroke plasma serotonin levels compared to healthy subjects [
30] while some reported no significant differences in serotonin levels [
31]. The suggested reason could be in the dynamics of the platelet serotonin (PS) content in one- day period after stroke, since the patients with low PS content had more severe status than patients with values greater than normal [
31]. Serotonin is well known as excitatory glutamatergic neurotransmitter that modulates neural activity and regulates mood and behaviour, but it also very important regulator of platelet aggregation and cardiovascular function [
32]. Abnormal serotoninergic mechanism may lead to a pro-thrombotic state, while downregulation may increase the risk of bleeding [
33]. Many studies have shown that selective serotonin reuptake inhibitors (SSRIs) can improve clinical outcome from ischemic stroke [
34] possibly due to the effecting neuronal cell survival, synaptic plasticity and neuronal connections [
35].
Phenylalanine concentration was also increased in patients
vs. controls serums in this study (
Figure 2 and
Figure 3) with VIP > 2 (
Figure 1b), which is also in accordance with related studies, where it is suggested to be a compensatory response to high neurotoxic concentrations of glutamate, due to phenylalanine inhibition of excitatory glutamatergic synaptic transmission [
27].
Excessive release of
aspartate, another excitatory amino acid significantly elevated in patient serums (
Figure 2 and
Figure 3), is believed to be as a result of membrane depolarization after cerebral ischemia and related to excitotoxicity and oxidative stress [
36].
In microglial cells,
ornithine is a key substrate for polyamines biosynthesis, which are involved in regulation of inflammatory reactions and cell renewal processes after cerebral ischemia [
37], but the main metabolic fate for ornithine in the brain is most probably conversion to glutamate and GABA via ornithine aminotransferase [
38]. Therefore, increased ornithine concentration in patient serums (
Figure 2 and
Figure 3) suggests both involvement in neuronal damage and cell renewal after cerebral ischemia.
Methionine sulfoxide, a cell oxidative stress product, was significantly elevated in stroke thrombi compared to patients serums (
Figure 5 and
Figure 6; VIP > 1.5,
Figure 4b), as previously reported to be elevated in stroke patients’ plasma and associated with increased risk of ischemic stroke [
39,
40].
Beside several significantly altered amino acids that are involved in ischemic metabolic serum profile,
glucose concentrations have also been significantly increased in patient serums showed by both standard laboratory blood tests (
Table 2) and metabolomics serum analysis (
Figure 2 and
Figure 3; VIP > 1.5,
Figure 1b), as previously reported in many studies and considered to be an predictor of poor stroke outcome [
41]. Hyperglycaemia was reported to be present in almost 40 % of stroke patients [
42] and it is suggested to promote ischemic injury by several proposed mechanisms, including acidosis due to anaerobic metabolism of glucose to lactic acid, superoxide and nitric oxide production, enhanced glucose-sodium exchange, abnormal protein glycosylation and advanced glycation products [
43].
The majority of statistically significant metabolites in this study were lipids - phospholipids, sphingolipids and triacylglycerols, representing 7/13 metabolites between patient and control serums (
Figure 3) and 14/15 metabolites between stroke thrombi and patient serums (
Figure 6). The only lipid metabolite that was significantly elevated in patient serums
vs. controls was
ether phosphatidylcholine (PC-O [42:4]) (
Figure 2 and
Figure 3; VIP > 2,
Figure 1b), possibly having C22:0 fatty alcohol linked via ether linkage at the position C-1 and eicosatetranoic acid (C20:4) at the position C-2, which categorizes this ether phospholipid as plasmalogen [
44]. Ether phospholipids, including plasmalogens, have different biological functions such as modulation of membrane trafficking, cell signalling, oxidative status and storage of polyunsaturated fatty acids (PUFA) and their compositions and their contents are altered in the plasma of patients suffering from different brain disorders [
45]. Interestingly, all other ether phospholipids in this study (PC-O [34:0], PC-O [32:0], PC-O [34:1], PC-O [36:2],
Figure 5 and
Figure 6) were characterized in significantly higher concentrations in stroke thrombi compared to patient serum, with
PC-O (34:0) having the highest VIP score (VIP > 2;
Figure 4b), pointing out their potential as a stroke biomarker of oxidative stress and brain injury. On the other hand, “classical” phospholipid species such as
phosphatidyl choline (PC [32:0], PC [36:5], PC [38:7], PC [40:9]; VIP > 1.5,
Figure 4b) were significantly increased in patient serum samples compared to stroke thrombi (
Figure 5 and
Figure 6), moreover, all having very high abundance of PUFA content. The only exception was
PC (32:0) as the only phosphatidyl choline species that was significantly decreased in patient serums compared to stroke thrombi (
Figure 5 and
Figure 6). Interestingly, the fatty acids composition is identical as in PC-O (32:0), which is also much more expressed in thrombi than in patient serum samples. The only
lysophosphatidyl choline species characterized as statistically significant and discriminatory in this study was
LPC (18:2), which was lower in patient serum compared to healthy control serums (
Figure 3; VIP > 1.5,
Figure 1b). Choline-containing phospholipids are acetylcholine precursors and they have been proposed as adjuvant therapy in acute stroke due to their anti-inflammatory effects and choline’s involvement in membrane synthesis, which is suggested as beneficial for reducing cell injury in ischemic brain [
46].
Clinical laboratory analysis of all lipid parameters (TG, tCH, LDL-CH and HDL-CH) did not show significant differences between patient and control serums (
Table 2), despite their common association with high risk of atherosclerosis and subsequently ischemic stroke [
47]. In accordance to this, in our metabolomics study there were no significant differences detected in cholesterol or cholesterol ester species between these two groups. On the other hand,
triacylglycerol species (TG [48:1] and [48:2], VIP > 2; TG [50:3], [53:3] and [56:6], VIP > 1.5,
Figure 1b) were significantly decreased in patient serums compared to control serums (
Figure 2 and
Figure 3), while TG (52:4) and TG (54:3) (VIP > 1.5,
Figure 4b) were decreased in stroke thrombi compared to patient serums (
Figure 5 and
Figure 6). This is in accordance with previous studies, characterizing TG as bioenergetic compounds of lipid droplets that are also present in neural cells, which were reported to have low abundance in plasma patients having stroke recurrence after transient ischemic attack, possibly due to low accumulation or formation of cerebral lipid droplets as a product of a defective ischemia-associated stress response [
48].
From all sphingolipid metabolites characterized in this study,
sphingomyelins (SM [40:1], [42:1] and [44:1]) and
ceramide Cer (42:2) were detected as significantly higher in stroke thrombi than in patient serums (
Figure 5 and
Figure 6; VIP > 1.5,
Figure 4b). Taking into consideration only univariate analysis, from 26 detected SM species between patient and control serums, 11 SM species were significantly increased in patient serums (
Figure 2). Sphingolipids, beside their well-known role as structural membrane components, have important functions as bioactive, signaling molecules involved in multiple cellular processes and mechanisms related to ischemia-related stress response [
49], inflammation [
50], and also serving as second messengers during platelet activation in coagulation process [
51]. In thrombin-stimulated platelets, active sphingomyelinase enzyme hydrolyses sphingomyelin to phosphocholine and ceramide, which is involved in signal transduction events during platelet activation [
52]. Variations in sphingolipid metabolites plasma levels have also been associated with numerus metabolic and vascular diseases [
53]. Long-chain ceramides are involved in apoptotic pathways and inflammation related to cerebral ischemia [
54] and their elevated serum levels have been reported as predictors of risk and severity of the ischemic stroke [
55,
56] as well as related to adverse cardiovascular risk and events [
57]. Rise in ceramide levels is a result of inflammation and tissue damage as well as a product of a strong increase in acidic and/or neutral sphingomyelinase activity (ASM/NSM), producing ceramide from sphingomyelin, which was reported to occur in focal cerebral ischemia [
58]. The ASM / ceramide system is critically involved in ischemic stroke pathogenesis and studies have reported improvement in ischemic brain injury after ASM inhibition and lowering ceramide content [
59]. Sphingomyelins, as major constituents of lipid rafts, are involved in signaling cascades and their deficiency has shown to suppress the inflammation induced by cerebral ischemia damage after ischemic stroke in mice [
60] while their higher serum content was related to higher risk of myocardial infarction [
61]. Besides ceramides, sphingomyelin SM (44:2) was reported as a potential stroke biomarker [
62], SM (38:1) and Cer (34:1) were reported as biomarkers of cerebral microvascular disease [
63], while SM (32:1) levels were reported to be inversely related to incident ischemic stroke in large multi-cohort metabolomics study [
64]. The reason why some particular SM or Cer species are elevated or decreased in different pathologies remains unclear.
4.1. Study Limitations
Our study has several limitations. The research was conducted using a relatively small number of subjects, which may influence the statistical power of the study and affect generalizability of study results due to the heterogeneity of aetiologies, comorbidities, clinical presentations and demographics of ischemic stroke patients. Since the patients’ blood was not collected after overnight fasting, potential nutritional influence on metabolite profile cannot be excluded. Furthermore, other stroke risk factors such as physical activity, smoking and alcohol consummation were not evaluated. Age differences between the analysed groups may also influence the accuracy of the compared metabolic profiles. Therefore, this research could be validated and confirmed in a larger and more diverse study.
Author Contributions
Conceptualization, D.F. and I.K.; methodology, D.F, I.K., I.R., D.R.M, J.K., D.O., V.K., M.J.S., R.K. and K.S.; validation, I.K., D.F., I.R. and K.S.; formal analysis, D.F., I.K. and I.R.; investigation, D.F., I.K., I.R, K.S., and R.K.; resources, T.S., K.R., M.J.S., M.R.B., R.K., Z.P and V.M.; data curation, I.R., I.K. and D.F.; writing - original draft preparation, D.F, I.K., I.R. and K.S.; writing - review and editing, D.F, I.K., Z.P, M.R.B. and V.M.; visualization, I.K, D.F. and I.R.; supervision, D.F. and I.K.; project administration, D.F., I.K. and V.M.; funding acquisition, V.M, D.F., I.K. All authors have read and agreed to the published version of the manuscript.