Impact Factor 2.889 | CiteScore 3.1

More on impact ›

Original Research ARTICLE

Front. Neurol., 05 November 2020 | https://doi.org/10.3389/fneur.2020.575376

Machine Learning Analysis of the Cerebrovascular Thrombi Proteome in Human Ischemic Stroke: An Exploratory Study

Cyril Dargazanli1,2*, Emma Zub1, Jeremy Deverdun3, Mathilde Decourcelle4, Frédéric de Bock1, Julien Labreuche5, Pierre-Henri Lefèvre2, Grégory Gascou2, Imad Derraz2, Carlos Riquelme Bareiro2, Federico Cagnazzo2, Alain Bonafé2, Philippe Marin1, Vincent Costalat1,2* and Nicola Marchi1*
  • 1Institut de Génomique Fonctionnelle, Univ. Montpellier, UMR 5203 CNRS - U 1191 INSERM, Montpellier, France
  • 2Diagnostic and Interventional Neuroradiology Department, Gui de Chauliac Hospital, Montpellier, France
  • 3I2FH, Institut d'Imagerie Fonctionnelle Humaine, Gui de Chauliac Hospital, Montpellier, France
  • 4BioCampus Montpellier, CNRS, INSERM, Université de Montpellier, Montpellier, France
  • 5Santé Publique: Epidémiologie et Qualité des Soins, CHU Lille, University of Lille, Lille, France

Objective: Mechanical retrieval of thrombotic material from acute ischemic stroke patients provides a unique entry point for translational research investigations. Here, we resolved the proteomes of cardioembolic and atherothrombotic cerebrovascular human thrombi and applied an artificial intelligence routine to examine protein signatures between the two selected groups.

Methods: We specifically used n = 32 cardioembolic and n = 28 atherothrombotic diagnosed thrombi from patients suffering from acute stroke and treated by mechanical thrombectomy. Thrombi proteins were successfully separated by gel-electrophoresis. For each thrombi, peptide samples were analyzed by nano-flow liquid chromatography coupled to tandem mass spectrometry (nano-LC-MS/MS) to obtain specific proteomes. Relative protein quantification was performed using a label-free LFQ algorithm and all dataset were analyzed using a support-vector-machine (SVM) learning method. Data are available via ProteomeXchange with identifier PXD020398. Clinical data were also analyzed using SVM, alone or in combination with the proteomes.

Results: A total of 2,455 proteins were identified by nano-LC-MS/MS in the samples analyzed, with 438 proteins constantly detected in all samples. SVM analysis of LFQ proteomic data delivered combinations of three proteins achieving a maximum of 88.3% for correct classification of the cardioembolic and atherothrombotic samples in our cohort. The coagulation factor XIII appeared in all of the SVM protein trios, associating with cardioembolic thrombi. A combined SVM analysis of the LFQ proteome and clinical data did not deliver a better discriminatory score as compared to the proteome only.

Conclusion: Our results advance the portrayal of the human cerebrovascular thrombi proteome. The exploratory SVM analysis outlined sets of proteins for a proof-of-principle characterization of our cohort cardioembolic and atherothrombotic samples. The integrated analysis proposed herein could be further developed and retested on a larger patients population to better understand stroke origin and the associated cerebrovascular pathophysiology.

Introduction

Stroke is a major public health burden and the second most common cause of death worldwide (13). Currently, the incomplete molecular understanding of stroke pathophysiology negatively impacts patients' management, follow-up, and secondary prevention (3, 4). A recent consensus indicates that examinations of patients' intracranial thrombi could help unveil novel disease mechanisms (5). Studying the intracranial thrombi composition could advance our knowledge of the molecular mechanisms of local cerebrovascular cell damage in this disease setting (69).

Mechanical thrombectomy (MT) is a standard of care for patients presenting with acute ischemic stroke (AIS) due to large vessel occlusion (LVO) (10). MT allows the retrieval of cerebral thrombi from brain arteries, enabling subsequent samples storage and analysis. A few studies have analyzed the histological composition of intracranial thrombi (11, 12), describing architecture or reporting the presence of fibrin and leucocytes (13). However, an in depth characterization of the thrombi molecular components is currently lacking (11).

Here, we performed a quantitative proteomic analysis of intracranial thrombi retrieved using MT from a cohort of n = 32 cardioembolic and n = 28 atherothrombotic diagnosed AIS patients. We resolved the thrombi proteomes for our cohort samples and next applied a support-vector machine (SVM) learning approach to estimate whether specific sets of proteins, alone or in combination with available clinical data, could help differentiate the cardioembolic from atherothrombotic origin in our selected population.

Materials and Methods

Inclusion Criteria

Patients with suspected ischemic stroke secondary to an LVO were prospectively recruited at a high-volume, comprehensive stroke center in France. Patients were required to present imaging evidence of occlusion of the internal carotid artery (ICA, cervical or intracranial part), the M1 or M2 branches of the middle cerebral artery (MCA), the basilar artery, or a tandem atheromatous occlusion defined by the occlusion of both cervical carotid artery and intracranial artery (carotid artery or MCA). Use of intravenous thrombolysis (IVT) treatment was allowed and administrated according to current guidelines (10). Stroke cause was defined by a stroke neurologist blinded to the proteomics analysis, according to the TOAST (Trial of ORG 10172 in Acute Stroke Treatment) (14) classification, after an exhaustive in-hospital workup (15) including at least computed tomography and magnetic resonance imaging, duplex sonography of the cervical arteries, blood coagulation tests, long-term electrocardiography, and transthoracic or transesophageal echocardiography. Stroke etiology was defined as “atherothrombotic tandem” when CT angiography and MR angiography demonstrated >50% stenosis or occlusion of the cervical carotid artery with associated intracranial ICA or MCA occlusion ipsilateral to the symptomatic hemisphere, in addition to exclusion of potential sources of cardiac embolism. Stroke etiology was defined as “cardioembolic” when at least one cardiac source for an embolus was identified after a complete cardiological work-up including Holter monitoring and echocardiography, in the absence of any stenosis of ipsilateral large extracranial arteries or atherosclerosis, excluding atrial fibrillation with non-cardioembolic strokes.

Exclusion criteria for the present study were: (1) failure of thrombus retrieval (failure of catheterization, patients with spontaneous reperfusion at the beginning of the procedure), (2) patients non-suitable for MT with a pre-stroke modified Rankin Scale (mRS) score of >3; (3) patients with non-atheromatous or non-cardioembolic tandem occlusions (intimal dysplasia/web, dissection), (4) patients having had MT but with a thromboembolic material unsuitable for proteomic analyses (mainly due to insufficient material amounts retrieved), (5) patients with no clear etiology or “undefined etiology” (defined as at least two possible etiologies found after a complete clinical, laboratory, and imaging work-up).

The study was approved by the local ethics committee, with the patients providing written informed consent in acute phase whenever possible. Otherwise, the consent form was signed by the patient's relatives.

Patient Characteristics

Patient demographics, vascular risk factors, imaging data, vital signs before treatment, severity of ischemic stroke, and clinical outcomes were collected with a structured questionnaire. Age, sex, cardiovascular risk factors (hypertension, dyslipidemia, diabetes mellitus, and smoking habits), time of symptom onset, National Institutes of Health Stroke Scale (NIHSS) at baseline, use of IVT, and its time from symptom onset were collected (see Table 1). The Alberta Stroke Program Early CT Score (ASPECT) on diffusion-weighted magnetic resonance or CT imaging was assessed by a neuroradiologist.

TABLE 1
www.frontiersin.org

Table 1. Patients data, Treatment Characteristics, Complications, and Outcomes according to stroke etiology.

Endovascular Procedure

All patients were treated in a dedicated neuroangiography suite under general anesthesia or conscious sedation, after evaluation by the anesthesiology team.

Most of the procedures were performed using the Trevo® device (Stryker, Kalamazoo, Michigan) or the Solitaire FR™ device (Medtronic, Dublin, Ireland) via the femoral artery approach. A balloon catheter was positioned in the ICA to allow flow arrest during thrombus retrieval. The stent retriever was delivered through a microcatheter and deployed across the thrombus. A distal aspiration during the stent retrieval was performed, according to the SAVE technique (16). A control angiogram was obtained to assess recanalization and reperfusion. This sequence was repeated until mTICI 2b or mTICI 2C/3 flow (defined as successful reperfusion) was established (17). The “retrograde approach” (also known as the distal-to-proximal or intracranial-first approach), aiming to recanalize the distal and symptomatic intracranial occlusion before addressing the cervical carotid lesion, was generally chosen for tandem occlusions. The interventional neuroradiologist used another thrombectomy device in the case of reperfusion failure (mTICI <2b) with the first stent retriever. Reperfusion results were reported by using the mTICI score (18). Peri-procedural complications [embolization in a new territory (defined as an angiographic occlusion in a previously unaffected vascular territory observed on the angiogram after clot removal), arterial dissection or perforation, vasospasm, and subarachnoid hemorrhage] were recorded.

Follow-Up and Outcome

All patients underwent cross-sectional imaging (computed tomography or magnetic resonance imaging) within a range of 18–24 h after the procedure. Intracranial hemorrhage was classified according to the ECASS (European Cooperative Acute Stroke Study) criteria (19). Symptomatic intracranial hemorrhage was defined as any intracerebral hemorrhage with an increase of at least four NIHSS points within 24 h, or resulting in death. The mRS at 90 days was assessed by trained research nurses unaware of the study group assignments, during face-to-face interviews, or via telephone conversations with the patients, their relatives, or their general practitioners.

Collection and Processing of Intracranial Thrombi

In the angiography room, after retrieval (Figure 1E), thrombi were immediately frozen at −80°C in a dedicated transportable azote freezer (Voyager, Air Liquide). In the laboratory, samples were prepared for mass spectrometry analysis. After initial mashing in a glass potter at 4°C in RIPA buffer, thrombi were further dissolved using an ultrasonic liquid processor (10 applications of 1 second each at 4°C; Vibra-cell VCX130PB, VWR) and then centrifuged (Eppendorf 5427R) at 1,200 RPM for 10 min at 4°C. Protein concentration was assessed by a bicinchoninic acid (BCA) assay. Protein extracts (20 μg) were separated by SDS-PAGE using a short (2 cm) migration. Single pieces of gel including separated proteins except hemoglobin were excised for each sample and proteins were digested in-gel using Trypsin (Trypsin Gold, Promega), as previous described (20).

FIGURE 1
www.frontiersin.org

Figure 1. Cerebral angiography showing a right MCA occlusion before (A: anteroposterior projection, B: lateral projection) and after recanalization (C: anteroposterior projection, D: lateral projection) achieved by mechanical thrombectomy. (E) Clot engaged by the Trevo® stent-retriever. (F) Illustration of thrombi protein separation on 4–15% polyacrylamide gels and Hemoglobin depletion (*) prior to in-gel protein digestion by trypsin.

Mass Spectrometry

The resulting peptide samples were analyzed online using Q-Exactive HF mass spectrometer coupled with an Ultimate 3000 RSLC (Thermo Fisher Scientific) fitted with a stainless-steel emitter (Thermo Fisher Scientific). Desalting and pre-concentration of samples were performed online on a Pepmap® pre-column (0.3 mm × 10 mm, Thermo Fisher Scientific). A gradient consisting of 2–40% of buffer B in 123 min, then 90% of buffer B during 5 min (A: 0.1% formic acid in water; B: 0.1% formic acid 80% ACN) at 300 nL/min was used to elute peptides from the capillary reverse-phase column (0.075 mm × 500 mm, Pepmap® C18, Thermo Fisher Scientific). Spectra were acquired with Xcalibur software (v4.1 Thermo Fisher Scientific). MS/MS analyses were performed in a data-dependent mode with standard settings. MS data analysis was performed using the MaxQuant software with default settings (v1.5.5.1) (21). All MS/MS spectra were searched using the Andromeda search engine (22) against the UniProtKB Reference proteome UP000005640 database for Homo sapiens (release 2019_03, https://www.uniprot.org/) and the contaminant database in MaxQuant. Default search parameters were applied. Oxidation (Met) and Acetylation (N-term) were used as variable modifications and Carbamidomethyl (Cys) was used as fixed modification. FDR was set to 1% for peptides and proteins. A representative protein ID in each protein group was automatically selected using an in-house bioinformatics tool (Leading v3.2). First, proteins with the most numerous identified peptides were isolated in a “match group” (proteins from the “Protein IDs” column with the maximum number of “peptides counts”). For the match groups where more than one protein ID were present after filtering, the best annotated protein in UniProtKB [reviewed entries rather than automatic ones, highest evidence for protein existence, most annotated protein according to the number of Gene Ontology Annotations (UniProtKB-GOA, made on 20190416)] was defined as the “leading” protein. Graphical representation and statistical analysis of MS/MS data were performed using Perseus (v1.6.1.1). Label-free quantification (MaxQuant LFQ) was used to highlight proteins differentially expressed between samples (23).

The mass spectrometry proteomics data have been uploaded to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD020398 (24).

Data Analysis

Descriptive Analysis

Data in Table 1 are presented as median (range) for quantitative variables, and percentage (count) for categorical variables. Baseline and treatment characteristics, complications and outcomes were compared according to stroke etiology using Chi-Square or Fisher's exact tests for categorical variables and the Mann-Whitney U-test for quantitative variables. No statistical comparisons were done for categorical variables with frequency <5. Statistical testing was done at the 2-tailed α level of 0.05. Data were analyzed using the SAS package, release 9.4 (SAS Institute, Cary, NC).

A support-vector machine (SVM) approach was implemented using MATLAB (r2018a, MathWorks, Natick, MA, USA). The SVM algorithm analyzes and learns from the dataset (Supplementary Table 2) to identify the hyperplanes for the best segregation of data according to a known discriminatory characteristic (25). In our work, the relatively small sample size prevents from achieving a correct validation step and SVM was used as a statistical tool to examine whether hyperplanes exist splitting our two groups. Here, we specifically tested whether samples segregation is attainable using combinations of up to 3 proteins (trios) from those commonly detected in all samples. Each possible combinations of three proteins from the data set in Supplementary Table 2 was tested (n = 13,908,836), the corresponding X/Y/Z hyperplanes were defined by the SVM (see Figure 3A), and the percentage of correct sample classification was obtained. The protein combinations achieving the best discriminatory score for our populations were retained. SVM analysis was also performed using clinical data in Table 1.

Results

Clinical Data, Peripheral Blood and Thrombi Characteristics

Baseline clinical data, treatment characteristics, early complications and outcomes are provided in Table 1. In the selected population, subjects suffering from atherothrombotic stroke were younger (67.5 vs. 79.5 years old, p = 0.005), presented no cardiac failure (0 vs. 18.8%, p = 0.047), no significant atrial fibrillation (3.6 vs. 62.5%, p < 0.001), and displayed higher systolic and diastolic blood pressure at admission (152 and 90 mmHg vs. 136 and 80 mmHg, p = 0.006 and 0.033). M1 occlusions were more frequent in the cardioembolic group (56.3 vs. 7.1%, p < 0.001). Groin puncture to reperfusion time was longer in the atherothrombotic group, which included 85.7% of tandem occlusions (72 vs. 40 min., p = 0.002). Complete blood count at admission indicated that platelet levels were higher in the atherothrombotic group (250 × 109/L vs. 203 × 109/L, p = 0.011; Table 1). Weight of the retrieved thrombi was 31.2 mg for the cardioembolic group (range 5.8–206.2 mg) and 36 mg for the atherothrombotic group (range 3.2–136.2; p = 0.85). Total protein concentrations were 11.20 μg/μl (5.3–22.1) and 11.1 μg/μl (4–26.5; p = 0.82) for the cardioembolic and atherothrombotic groups, respectively.

Analysis of the Intracranial Human Thrombi Proteome

All thrombus samples were individually processed by SDS-page chromatography and the hemoglobin band excised (Figures 1 A–F). Mass spectrometry analysis identified a total of 2,455 proteins in the samples analyzed. The complete list of all proteins detected in each sample is provided in Supplementary Table 1. A total of 438 proteins were commonly present in all the samples analyzed (Supplementary Table 2). Analysis of ClueGO annotations of the thrombi proteome, according to UniProtKB or EBI GOA databases, showed protein clusters for key biological pathways including metabolic processes, cytokines production, and cell proliferation, activation, or motility (Figure 2A). Indicating an inflammatory track are proteins associated with leukocytes activation, migration, and cell adhesion (Figure 2B; high definition zoom-in). This dataset constitutes the largest human thrombus proteome available and a shared library for the investigation of the molecular mechanisms of thrombus formation and ischemic stroke pathophysiology.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Graphic representation of the proteome fingerprint and protein clustering according to cellular functions. (B) List of specific cellular processes color-coded to panel A (high resolution zoom-in).

Exploring the Use of Support-Vector-Machine Learning to Analyse the Thrombi Proteome

The proteomic LFQ data obtained from our samples cohort were analyzed using a SVM routine to mathematically examine potential signatures existing between the cardioembolic and atherothrombotic proteomes. The SVM algorithm does not handle missing data across samples and the analysis was performed using the proteins commonly detected in all thrombi (438 proteins; Supplementary Table 2). In our SVM study we specifically aimed at identifying small set of discriminatory elements, here up to 3 proteins (see Methods). As a result, proteins trios were found by SVM providing a 88.3% accuracy of correct classification of our two sample groups. Proteins and their biological functions are detailed in Table 2. Factor XIII, which catalyzes the last step of the coagulation cascade by crosslinking fibrin fibers, was present in all combinations. Figure 3A shows an illustration of the SVM hyperplane classification for the cardioembolic and atherothrombotic samples according to the protein trio Eukaryotic translation initiation factor 2 subunit 3, Ras GTPase-activating-like protein IQGAP2, and Coagulation factor XIII. Using this specific setting, four and three patients were misclassified (light green squares in Figure 3A) as cardioembolic and atherothrombotic, respectively. In univariate analysis (Wilcoxon test), the coagulation Factor XIII, the Eukaryotic translation initiation factor 2 subunit 3, and the Myosin light chain kinase levels were significantly different between the cardioembolic and atherothrombotic groups, with respective p-values of 0.002, 0.04, and 0.01 (see Table 2). These results have a dual value, suggesting potential molecular differences between cardioembolic and atherothrombotic thrombi while supporting the notion of protein biomarkers to understand clot origin.

TABLE 2
www.frontiersin.org

Table 2. SVM protein trios allowing 88.3% accuracy of correct classification of cardioembolic and atherothrombotic thrombi.

FIGURE 3
www.frontiersin.org

Figure 3. (A) Example of SVM plot and classification for the combination of Eukaryotic translation initiation factor 2 subunit 3, Ras GTPase-activating-like protein IQGAP2, and coagulation Factor XIII. Groups A: atherothrombotic; Group C: cardioembolic. Light green squares and red circles indicate labeling errors for misclassified patients. (B) Volcano plot showing proteins significantly (p < 0.05) enriched using LFQ values log2 transformed (red dots; see Table 3 for protein details). Student's t-test is performed by using Perseus algorithms. Blue dots indicate the SVM-identified proteins (see Table 2 for details).

Integrating SVM Analyses of Clinical Data and Thrombi Proteome

In an attempt to identify additional SVM differentiation factors, we performed an analysis using patients clinical data (Table 1; age, sex, history of cardiac failure or atrial fibrillation, previous antithrombotic medication, glycemia, weight and BMI, thrombus weight and global protein concentration, hemoglobin, leucocytes, and platelet rate). SVM identified history of cardiac failure and atrial fibrillation as variables differentiating the two population with an 81.36% accuracy. This result is obvious considering our study design and because history of cardiac failure was used as one of the criteria to diagnose etiology at enrollement (see Methods). Cardiac failure and atrial fibrillation are two known risk factors linked to cardioembolic stroke (3). Interestingly, when atrial fibrillation was excluded from the SVM analysis, patient age and thrombus protein concentration provided a differentiation level of 74.58% within our sample cohorts. The latter results indicate thrombus total protein concentration as a new SVM analytical variable. Addition of a third variable did not improve the SVM score (not shown). We do acknowledge that combining the protein trio 1 (see Table 2), history of cardiac failure, and protein concentration we obtained a SVM score of 96.6%.

Testing Proteome Using LFQ Statistics

The selected SVM method tests all combinations of three inter-dependent proteins, obtaining solutions for data clusterization that are not executable using LFQ and standard statistics (26). Thus, a Student's T-test (Perseus algorithms) analysis on the proteins (log2 transformed) detected in all samples did not deliver significant difference between the studied cardioembolic and atherothrombotic populations. Furthermore, we applied a conventional method where proteomes (Supplementary Table 1) are filtered to include proteins with at least 50% of valid LFQ values. By using this approach, Student's t-test identified four proteins (PHB, SLC25A11, ATP5A1, and APOE; see Table 3) that display an abundance in cardioembolic as compared to atherothrombotic thrombi (volcano plot in Figure 3B). However, LFQ T-test difference was low (x-axis = −1.2; red dots in Figure 3B) with the crucial caveat that, because of method design, these proteins were undetectable in an elevated number of samples (Supplementary Table 1), therefore impeding group discrimination. These results support the relevance and the efficiency of SVM to analyze the proteome thombi dataset in our experimental settings.

TABLE 3
www.frontiersin.org

Table 3. LFQ (log2 transformed) of single proteins enriched in cardioembolic as compared to atherothrombotic thrombi.

Discussion

Our study advances the knowledge of the human cerebrovascular thombi composition by delivering the largest proteome dataset available to date. We focused on the protemic analysis of cardioembolic and atherothrombotic thrombi and we applied a support-vector machine learning routine in an exploratory, proof-of-concept, attempt to identify protein candidates segragating the two selected populations. Our research supports the general notion that direct analysis of the thrombi material could unveil, in the future, new disease players and candidate biomarkers potentially aiding stroke diagnosis. The SVM method used herein was set to identify combinations of protein trios (Table 2) in the intracranial thrombi, and it allowed for an 88.3% correct classification of our selected cardioembolic and atherothrombotic populations (Table 1). We here underscore that histological, cellular (e.g., red blood cells, platelets, white blood cells), and molecular (omics) analyses should be all integrated to obtain a complete and multi-level depiction of the thrombi structure and biology.

Understanding the composition of the human clot was previously attempted in two studies, although limited in sample size or lacking SVM analysis (12, 27) A first proteomic investigation correlated 2 inflammation-associated proteins (integrin alpha-M and mitochondrial superoxide dismutase) to high blood LDL (27). Mitochondrial superoxide dismutase was previously associated to unstable carotid plaques (28). These proteins were detected in our study, although without significant differences between cardioembolic and atherothrombotic thrombi. A second study analyzed 4 thrombi, with ~1,600 proteins identified (12). An earlier investigation, focused on human coronary thrombi in patients with ST-segment elevation in acute myocardial infarction, identified 708 proteins. The implication of platelet activation during the formation of thrombus causing acute coronary syndrome was suggested (29).

Combining Mass-Spectrometry With SVM Analysis: Initial Feasibility and Proposed Applicability to Human Ischemic Stroke

An innovative aspect of the presented study is the methodological combination of large-scale proteomic tools and machine learning models or algorithms to define and potentially categorize the thrombi proteomes (3). In our patients' cohort, the fibrin stabilizing or coagulation Factor XIII (FXIII) was identified by SVM as one potential differentiating element between the cardioembolic and atherothrombotic thrombi analyzed (Table 2). FXIII is a key enzyme in the coagulation cascade that allows the cross-linking of fibrin chains with subsequent increase of mechanical clot strength and resistance to fibrinolysis (30). FXIII was also reported in embolized thrombi from the cardiac left atrial appendage in atrial fibrillation patients (31).

Interestingly, it has been recently shown that FXIII levels are higher after myocardial injury and that FXIII harbors an important role in cardiac healing and remodeling (32). Moreover, valine-to-leucine (V34L) single-nucleotide polymorphism (SNP), which is associated with higher levels of FXIIIa, appears to be associated with a lower risk of pathological thrombosis in ischemic heart disease (33, 34). Importantly, atrial fibrillation or atrial cardiopathies that share a common mechanism of thrombus formation in the left atrial appendage should be identified as soon as possible after stroke occurrence to initiate anticoagulation therapy (35). Our SVM learning analysis also identified proteins involved in the cellular cytoskeleton assembly (Table 2), namely the myosin light chain kinase and F-actin-capping protein. In general, the large scale proteomic analysis of human clots here executed discloses pathways and molecular players of clot-endothelium interplay and local inflammation related to cerebrovascular damage (Figure 2). The latter is important because cerebrovascular breakdown contributes to the development of central nervous system disease (68, 36), in this case potentially enabling post-stroke sequelae.

Study Limitations and Prospectives

To further explore the utility of the protein candidates here discovered (Table 2) a validation step using an independent, and larger sample population will be necessary to define reproducibility and accuracy parameters (e.g., sensitivity, specificity, positive ad negative predictive values). Our SVM analysis, due to a relatively small sample size, only allowed accuracy estimation. A compelling question is whether our integrated proteomic-SVM method could be next used to examine specific signatures in case of cryptogenic stroke. We are aware that the proteins here identified may be not helpful in a population of cryptogenic stroke that includes etiologies other that the two studied here. We are aware that an efficient transition from SVM proteome analysis to clinical laboratory tools (e.g., Elisa) could be challenging and time consuming. (12, 27). The latter will be possible only when definitive molecular candidate(s) will be confirmed in larger populations with results replicated across stroke centers. Nevertheless, our study provide a proof of principle model that could be further developed and applied. Our proteome results (Supplementary Tables 1, 2) are shared and available to be re-analyzed using more advanced or alternative SVM methods.

We here recognize that the cohort used in the present study is heterogeneous in respect to age and blood platelet levels. Although blood platelet levels have been associated to stroke outcome (37), it is unknown whether a correlation with stroke etiology exists. One study showed that high platelet content of intracranial thrombi associates with large artery atherosclerosis. However, the authors did not study the correlation between blood platelet content and stroke cause (38). Another possible limitation of our approach concerns the retrieved material that may not represent the entire thrombus, although the analyses presented here were performed on the largest portion of clots retrieved at one pass of the thrombectomy device. IVT may also alter the samples, although this effect is likely to be limited due to the short time between IVT and thrombus extraction and processing. Finally, pre-stroke antithrombotic therapy may alter thrombus proteome composition (39).

Conclusions

In summary, quantitative proteomics and SVM analysis can be feasibly combined to examine the variation of intracranial human thrombi proteomes. If further developed and tested on larger cohorts, these methods have the potential to discover precise and novel pathophysiological players and biomarkers, with the ideal goal of aiding cerebrovascular stroke diagnosis and secondary prevention.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [1] partner repository with the dataset identifier PXD020398.

Ethics Statement

The studies involving human participants were reviewed and approved by Comité de Protection des Personnes ≪Sud-Méditerranée IV≫, Centre Hospitalier Universitaire de Montpellier, hôpital Saint-Eloi, 34295 Montpellier Cedex 5. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

CD, JD, PM, VC, and NM: conception and design of the study, analysis of data, and drafting of the manuscript. EZ, MD, FB, and JL: acquisition and analysis of data, drafting of the manuscript, and figures. PH-L, GG, ID, CR, FC, and AB: acquisition of data. CD and VC: emergency surgery interventions, samples collection and patients’ approval. All authors contributed to the article and approved the submitted version.

Funding

Funds from Stryker Neurovascular were used to performed this study. Stryker was not involved in study design, monitoring, data collection, statistical analysis or interpretation of results.

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.

Acknowledgments

Mass spectrometry experiments were carried out using facilities of the Functional Proteomics Platform of the Proteomics Pole of Montpellier. We would like to thank Leonie Runtz for initial testing. We also thank Marine Blaquiere (IGF) for her technical involvement.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.575376/full#supplementary-material

Supplementary Table 1. Complete list of proteins and mass-spectrometry data for each sample.

Supplementary Table 2. List of proteins commonly present in all samples.

References

1. Amarenco P, Bogousslavsky J, Caplan LR, Donnan GA, Hennerici MG. Classification of stroke subtypes. Cerebrovasc Dis Basel Switz. (2009) 27:493–501. doi: 10.1159/000210432

CrossRef Full Text | Google Scholar

2. Ornello R, Degan D, Tiseo C, Carmine CD, Perciballi L, Pistoia F, et al. Distribution and temporal trends from 1993 to 2015 of ischemic stroke subtypes: a systematic review and meta-analysis. Stroke. (2018) 49:814–9. doi: 10.1161/STROKEAHA.117.020031

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Yaghi S, Bernstein RA, Passman R, Okin PM, Furie KL. Cryptogenic stroke: research and practice. Circ Res. (2017) 120:527–40. doi: 10.1161/CIRCRESAHA.116.308447

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Jickling GC, Sharp FR. Biomarker panels in ischemic stroke. Stroke J Cereb Circ. (2015) 46:915–20. doi: 10.1161/STROKEAHA.114.005604

PubMed Abstract | CrossRef Full Text | Google Scholar

5. De Meyer SF, Andersson T, Baxter B, Bendszus M, Brouwer P, Brinjikji W, et al. Analyses of thrombi in acute ischemic stroke: a consensus statement on current knowledge and future directions. Int J Stroke. (2017) 12:606–14. doi: 10.1177/1747493017709671

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sweeney MD, Zhao Z, Montagne A, Nelson AR, Zlokovic BV. Blood-brain barrier: from physiology to disease and back. Physiol Rev. (2019) 99:21–78. doi: 10.1152/physrev.00050.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Librizzi L, de Cutis M, Janigro D, Runtz L, de Bock F, Barbier EL, et al. Cerebrovascular heterogeneity and neuronal excitability. Neurosci Lett. (2018) 667:75–83. doi: 10.1016/j.neulet.2017.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Giannoni P, Badaut J, Dargazanli C, De Maudave AF, Klement W, Costalat V, et al. The pericyte-glia interface at the blood-brain barrier. Clin Sci Lond Engl 1979. (2018) 132:361–74. doi: 10.1042/CS20171634

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nation DA, Sweeney MD, Montagne A, Sagare AP, D'Orazio LM, Pachicano M, et al. Blood-brain barrier breakdown is an early biomarker of human cognitive dysfunction. Nat Med. (2019) 25:270–6. doi: 10.1038/s41591-018-0297-y

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Furie KL, Jayaraman MV. 2018 guidelines for the early management of patients with acute ischemic stroke. Stroke. (2018) 49:509–10. doi: 10.1161/STROKEAHA.118.020176

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Brinjikji W, Duffy S, Burrows A, Hacke W, Liebeskind D, Majoie CBLM, et al. Correlation of imaging and histopathology of thrombi in acute ischemic stroke with etiology and outcome: a systematic review. J NeuroInterventional Surg. (2017) 9:529–34. doi: 10.1136/neurintsurg-2016-012391

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Muñoz R, Santamaría E, Rubio I, Ausín K, Ostolaza A, Labarga A, et al. Mass spectrometry-based proteomic profiling of thrombotic material obtained by endovascular thrombectomy in patients with ischemic stroke. Int J Mol Sci. (2018) 19:498. doi: 10.3390/ijms19020498

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Dargazanli C, Rigau V, Omer E. High CD3+ cells in intracranial thrombi represent a biomarker of atherothrombotic stroke. PLoS ONE. (2016) 11:e0154945. doi: 10.1371/journal.pone.0154945

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Adams HP, Bendixen BH, Kappelle LJ, Biller J, Love BB, Gordon DL, et al. Classification of subtype of acute ischemic stroke. Definitions for use in a multicenter clinical trial. TOAST. trial of org 10172 in acute stroke treatment. Stroke. (1993) 24:35–41. doi: 10.1161/01.STR.24.1.35

PubMed Abstract | CrossRef Full Text | Google Scholar

15. McMahon NE, Bangee M, Benedetto V, Bray EP, Georgiou RF, Gibson JME, et al. Etiologic Workup in Cases of Cryptogenic Stroke. Stroke. (2020) 51:1419–27. doi: 10.1161/STROKEAHA.119.027123

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Maus V, Behme D, Kabbasch C, Borggrefe J, Tsogkas I, Nikoubashman O, et al. Maximizing first-pass complete reperfusion with save. Clin Neuroradiol. (2018) 28:327–38. doi: 10.1007/s00062-017-0566-z

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dargazanli C, Fahed R, Blanc R, Gory B, Labreuche J, Duhamel A, et al. Modified thrombolysis in cerebral infarction 2c/thrombolysis in cerebral infarction 3 reperfusion should be the aim of mechanical thrombectomy: insights from the aster trial (contact aspiration versus stent retriever for successful revascularization). Stroke. (2018) 49:1189–96. doi: 10.1161/STROKEAHA.118.020700

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Fugate JE, Klunder AM, Kallmes DF. What is meant by “tici”? Am J Neuroradiol. (2013) 34:1792–7. doi: 10.3174/ajnr.A3496

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Hacke W, Kaste M, Fieschi C, Toni D, Lesaffre E, von Kummer R, et al. Intravenous thrombolysis with recombinant tissue plasminogen activator for acute hemispheric stroke. The European cooperative acute stroke study (ECASS). JAMA. (1995) 274:1017–25. doi: 10.1001/jama.1995.03530130023023

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Shevchenko A, Tomas H, Havlis J, Olsen JV, Mann M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. (2006) 1:2856–60. doi: 10.1038/nprot.2006.468

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. (2008) 26:1367–72. doi: 10.1038/nbt.1511

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M. Andromeda: a peptide search engine integrated into the MaxQuant environment. J Proteome Res. (2011) 10:1794–805. doi: 10.1021/pr101065j

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, et al. the perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods. (2016) 13:731–40. doi: 10.1038/nmeth.3901

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. (2019) 47:D442–50. doi: 10.1093/nar/gky1106

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gholami R, Fakhari N. Support vector machine: principles, parameters, applications. In: Samui P, Sekhar S, Balas VE, editors. Handbook of Neural Computation. Academic Press (2017). p. 515–35. doi: 10.1016/B978-0-12-811318-9.00027-2

CrossRef Full Text | Google Scholar

26. Tyanova S, Albrechtsen R, Kronqvist P, Cox J, Mann M, Geiger T. Proteomic maps of breast cancer subtypes. Nat Commun. (2016) 7:10259. doi: 10.1038/ncomms10259

PubMed Abstract | CrossRef Full Text

27. Rao NM, Capri J, Cohn W, Abdaljaleel M, Restrepo L, Gornbein JA, et al. Peptide composition of stroke causing emboli correlate with serum markers of atherosclerosis and inflammation. Front Neurol. (2017) 8:427. doi: 10.3389/fneur.2017.00427

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lepedda AJ, Cigliano A, Cherchi GM, Spirito R, Maggioni M, Carta F, et al. A proteomic approach to differentiate histologically classified stable and unstable plaques from human carotid arteries. Atherosclerosis. (2009) 203:112–8. doi: 10.1016/j.atherosclerosis.2008.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Alonso-Orgaz S, Moreno-Luna R, López JA, Gil-Dones F, Padial LR, Moreu J, et al. Proteomic characterization of human coronary thrombus in patients with ST-segment elevation acute myocardial infarction. J Proteomics. (2014) 109:368–81. doi: 10.1016/j.jprot.2014.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Muszbek L, Yee VC, Hevessy Z. Blood coagulation factor XIII: structure and function. Thromb Res. (1999) 94:271–305. doi: 10.1016/S0049-3848(99)00023-7

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Gosk-Bierska I, McBane RD, Wu Y, Mruk J, Tafur A, McLeod T, et al. Platelet factor XIII gene expression and embolic propensity in atrial fibrillation. Thromb Haemost. (2011) 106:75–82. doi: 10.1160/TH10-11-0765

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Frey A, Gassenmaier T, Hofmann U, Schmitt D, Fette G, Marx A, et al. Coagulation factor XIII activity predicts left ventricular remodelling after acute myocardial infarction. ESC Heart Fail. 7:2354–64. doi: 10.1002/ehf2.12774

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Bagoly Z, Koncz Z, Hársfalvi J, Muszbek L. Factor XIII, clot structure, thrombosis. Thromb Res. (2012) 129:382–7. doi: 10.1016/j.thromres.2011.11.040

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Wartiovaara U, Mikkola H, Szôke G, Haramura G, Kárpáti L, Balogh I, et al. Effect of Val34Leu polymorphism on the activation of the coagulation factor XIII-A. Thromb Haemost. (2000) 84:595–600. doi: 10.1055/s-0037-1614073

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kamel H, Okin PM, Elkind MSV, Iadecola C. Atrial fibrillation and mechanisms of stroke: time for a new model. Stroke. (2016) 47:895–900. doi: 10.1161/STROKEAHA.115.012004

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Klement W, Blaquiere M, Zub E, deBock F, Boux F, Barbier E, et al. A pericyte-glia scarring develops at the leaky capillaries in the hippocampus during seizure activity. Epilepsia. (2019) 60:1399–411. doi: 10.1111/epi.16019

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Yang M, Pan Y, Li Z, Yan H, Zhao X, Liu L, et al. Platelet count predicts adverse clinical outcomes after ischemic stroke or TIA: subgroup analysis of CNSR II. Front Neurol. (2019) 10:370. doi: 10.3389/fneur.2019.00370

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Fitzgerald S, Dai D, Wang S, Douglas A, Kadirvel R, Layton KF, et al. Platelet-rich emboli in cerebral large vessel occlusion are associated with a large artery atherosclerosis source. Stroke. (2019) 50:1907–10. doi: 10.1161/STROKEAHA.118.024543

CrossRef Full Text | Google Scholar

39. Marcone S, Dervin F, Fitzgerald DJ. Proteomic signatures of antiplatelet drugs: new approaches to exploring drug effects. J Thromb Haemost. (2015) 13:S323–31. doi: 10.1111/jth.12943

PubMed Abstract | CrossRef Full Text | Google Scholar

40. UniProt Consortium T. UniProt: the universal protein knowledgebase. Nucleic Acids Res. (2018) 46:2699. doi: 10.1093/nar/gky092

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: stroke, thrombus, cerebrovascular, mechanical thrombectomy, proteome, support vector machine learning, neuroradiology

Citation: Dargazanli C, Zub E, Deverdun J, Decourcelle M, de Bock F, Labreuche J, Lefèvre P-H, Gascou G, Derraz I, Riquelme Bareiro C, Cagnazzo F, Bonafé A, Marin P, Costalat V and Marchi N (2020) Machine Learning Analysis of the Cerebrovascular Thrombi Proteome in Human Ischemic Stroke: An Exploratory Study. Front. Neurol. 11:575376. doi: 10.3389/fneur.2020.575376

Received: 23 June 2020; Accepted: 30 September 2020;
Published: 05 November 2020.

Edited by:

Bruce Campbell, The University of Melbourne, Australia

Reviewed by:

Simon F. De Meyer, KU Leuven, Belgium
Hidetoshi Kasuya, Tokyo Women's Medical University Medical Center East, Japan
Miriam Priglinger-Coorey, Royal North Shore Hospital, Australia

Copyright © 2020 Dargazanli, Zub, Deverdun, Decourcelle, de Bock, Labreuche, Lefèvre, Gascou, Derraz, Riquelme Bareiro, Cagnazzo, Bonafé, Marin, Costalat and Marchi. 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: Nicola Marchi, nicola.marchi@igf.cnrs.fr
Vincent Costalat, v-costalat@chu-montpellier.fr
Cyril Dargazanli, c-dargazanli@chu-montpellier.fr

These authors have contributed equally to this work