Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 12 June 2024
Sec. Stroke
This article is part of the Research Topic Omics-Based Approaches in Stroke Research View all 10 articles

Circulating miRNA profiles and the risk of hemorrhagic transformation after thrombolytic treatment of acute ischemic stroke: a pilot study

  • 1Department of Adult Neurology, Faculty of Medicine, Medical University of Gdańsk, Gdańsk, Poland
  • 2Department of Adult Neurology, University Clinical Center, Gdańsk, Poland
  • 3Brain Diseases Centre, Medical University of Gdańsk, Gdańsk, Poland
  • 4Laboratory for Regenerative Biotechnology, Department of Biotechnology and Microbiology, Gdańsk University of Technology, Gdańsk, Poland
  • 5Department of Biotechnology and Microbiology, Gdańsk University of Technology, Gdańsk, Poland
  • 6BioTechMed Center, Gdańsk University of Technology, Gdańsk, Poland

Background: Hemorrhagic transformation (HT) in acute ischemic stroke is likely to occur in patients treated with intravenous thrombolysis (IVT) and may lead to neurological deterioration and symptomatic intracranial hemorrhage (sICH). Despite the complex inclusion and exclusion criteria for IVT and some useful tools to stratify HT risk, sICH still occurs in approximately 6% of patients because some of the risk factors for this complication remain unknown.

Objective: This study aimed to explore whether there are any differences in circulating microRNA (miRNA) profiles between patients who develop HT after thrombolysis and those who do not.

Methods: Using qPCR, we quantified the expression of 84 miRNAs in plasma samples collected prior to thrombolytic treatment from 10 individuals who eventually developed HT and 10 patients who did not. For miRNAs that were downregulated (fold change (FC) <0.67) or upregulated (FC >1.5) with p < 0.10, we investigated the tissue specificity and performed KEGG pathway annotation using bioinformatics tools. Owing to the small patient sample size, instead of multivariate analysis with all major known HT risk factors, we matched the results with the admission NIHSS scores only.

Results: We observed trends towards downregulation of miR-1-3p, miR-133a-3p, miR-133b and miR-376c-3p, and upregulation of miR-7-5p, miR-17-3p, and miR-296-5p. Previously, the upregulated miR-7-5p was found to be highly expressed in the brain, whereas miR-1, miR-133a-3p and miR-133b appeared to be specific to the muscles and myocardium.

Conclusion: miRNA profiles tend to differ between patients who develop HT and those who do not, suggesting that miRNA profiling, likely in association with other omics approaches, may increase the current power of tools predicting thrombolysis-associated sICH in acute ischemic stroke patients. This study represents a free hypothesis-approach pilot study as a continuation from our previous work. Herein, we showed that applying mathematical analyses to extract information from raw big data may result in the identification of new pathophysiological pathways and may complete standard design works.

1 Introduction

Ischemic stroke was found to have an incidence of 7.6 million individuals worldwide in 2019, resulting in 63.48 million disability-adjusted life years (DALYs) and 3.29 million deaths. Ischemic stroke is a devastating neurological condition characterized by brain tissue damage caused by sudden obstruction of blood flow in the cerebral arteries (1, 2). Treatment in the acute phase aims to restore blood flow through intravenous thrombolysis and mechanical thrombectomy. The former method, which is used in up to 25% of patients, involves the administration of tissue-type plasminogen activator (rtPA), which promotes the formation of plasmin, a proteolytic enzyme. Plasmin breaks the crosslinks between fibrin molecules, leading to thrombus dissolution and restoration of blood flow (3, 4).

Hemorrhagic transformation (HT), which involves the extravasation of blood across a disrupted blood–brain barrier into the brain parenchyma, is one of the most common complications of ischemic stroke (5). According to the European Cooperative Acute Stroke Study (ECASS), HT can be categorized based on its intensity and radiological features into small petechial hemorrhagic infarction (HI1), confluent petechial hemorrhagic infarction (HI2), small parenchymal hemorrhage (PH1) (<30% infarct, mild mass effect), and large parenchymal hemorrhage (PH2, >30% infarct, marked mass effect) (6). Depending on its severity, HT may remain asymptomatic; however, if it is sufficiently large to exert a mass effect on brain tissue outside the infarct, it may cause neurological deterioration (7). Autopsy studies revealed hemorrhagic transformations in 18–42% of patients with acute ischemic stroke, and clinical assessment indicated symptomatic intracerebral hemorrhage after intravenous thrombolysis in approximately 6% of patients (8, 9).

Several studies aim at pinpointing reliable predictors of hemorrhagic transformation. The established clinical risk factors include baseline National Institutes of Health Stroke Scale (NIHSS) score, systolic and diastolic blood pressure, atrial fibrillation, antiplatelets use, age, and time from onset to treatment and hyperglycemia among others (10, 11). Radiological determinants of increased risk of hemorrhagic transformation include a large infarct size, early ischemic changes visible on computed tomography (CT), and absent or poor collaterals (10, 12). Among identified blood biomarkers, matrix metalloproteinase-9 (MMP-9), ferritin, and cellular fibronectin (c-Fn), as well as the neutrophil-to-lymphocyte ratio (NLR) and high-density lipoprotein (HDL), have been extensively studied across multiple experiments (13, 14).

Recent advances in artificial intelligence (AI) and omics have fostered their application in the search for novel HT biomarkers and predictive models. Machine learning methods have been used to develop predictive models based on clinical data and laboratory test results (15). In our previous study, we explored a hypothesis-free approach using MS proteomic data to identify new biomarkers (16). In that study, 15 proteins detected in the blood collected prior to rtPA treatment were unique to patients who developed HT.

MicroRNAs (miRNAs) are small non-coding RNA molecules composed of approximately 22 nucleotides that are known for their regulatory roles in various biological processes, mainly through the post-transcriptional regulation of gene expression (17). Their stability and detectability in various tissues, including blood, have attracted significant attention in the last decade, leading to their exploration as potential diagnostic and prognostic biomarkers, particularly in oncology (18). Circulating miRNAs have also emerged as valuable tools in stroke medicine. Numerous studies have identified miRNAs as diagnostic markers for ischemic stroke, with hsa-let-7e-5p, hsa-miR-124-3p, hsa-miR-17-5p, and hsa-miR-185-5p showing consistent differential expression (19). Furthermore, the combination of miR-124-3p, miR-125b-5p, and miR-192-5p expression has been shown to predict the extent of neurological deterioration in ischemic stroke patients treated with rtPA (20). In another study, miR-21-5p, miR-206, and miR-3123 were implicated in predicting the risk of hemorrhagic transformation in patients with cardioembolic stroke (21). Additionally, the assessment of RNA markers, including miRNA-23a, miRNA-193a, miRNA-128, miRNA-99a, miRNA-let-7a, miRNA-494, miRNA-424, and the long non-coding (lnc)RNA H19, has been shown to improve the prediction of symptomatic intracranial hemorrhage (sICH) after rtPA (22).

The findings of the above studies suggest that quantitative miRNA and proteomic data may increase the current power of the tools for predicting thrombolysis-associated sICH in patients with acute ischemic stroke [as we showed in our previous study (16)]. However, the main objective of the presented studies is to demonstrate a methodology for and the feasibility of such an approach. This pilot study only aimed to identify potential miRNAs indicative of an increased risk of HT occurrence.

2 Methods

2.1 Study population

This research is a continuation of our previous study on biomarkers of rtPA-associated intracranial hemorrhage in acute ischemic stroke (16). Participants were recruited between March 2019 and April 2022 from among patients at the Stroke Comprehensive Center of the Department of Adult Neurology, University Clinical Center, Medical University of Gdansk.

Before recruitment, approval was obtained from The Independent Bioethics Committee for Scientific Research at the Medical University of Gdańsk in Poland. This study was conducted in accordance with the principles outlined in the World Medical Association Declaration of Helsinki. All subjects were informed of the purpose and course of the study and signed an informed consent form.

For patients with acute stroke, a pretreatment MRI protocol (described in Section 2.2) was performed, in addition to the collection of two blood samples: plasma in an EDTA-coated tube and serum in a silica clot activator (Becton Dickinson Vacutainer). However, only plasma samples were used for further analyses. Upon confirmation of ischemic stroke diagnosis and patient adherence to thrombolysis inclusion criteria according to the stroke guidelines, alteplase was administered intravenously at a standard dose of 0.9 mg/kg. Subsequently, patients were assigned to either the HT or control group, based on follow-up MRI scans on day 5–9 after the treatment. In the initial cohort of 94 patients admitted for acute stroke, 56 received intravenous rtPA, met the standard inclusion and exclusion criteria, and had no contraindications for MRI. Due to limited funding during this phase of the study, 10 patients with HT were selected for further analysis and matched with 10 non-HT controls based on comparable risk factors for hemorrhagic transformation (NIHSS score, Oxfordshire Community Stroke Project classification) and similar miRNA profiles (age and sex).

The Propensity Score Matching method, which incorporates the k-nearest neighbor (k-NN) algorithm without replacement, was used to select a well-matched control group (non-HT) with a distribution of baseline characteristics similar to that of the HT group. The analysis was conducted using the MatchIt library in R software (23).

2.2 Brain imaging

Pretreatment MRI was performed using a Siemens Healthcare GOBrain application on a 1.5 T MRI-scanner (Siemens Magnetom Aera), comprising axial T2-weighted fluid attenuation inversion recovery, axial diffusion-weighted imaging (b-values 0,800 s/mm) with apparent diffusion coefficient maps, axial T2*-weighted, and sagittal T1-weighted images. These scans were used to assess the infarct volume and were analyzed by radiologists using the syngo.via software. The follow-up MRI scan was taken using the same 1.5 T Siemens device, but the protocol was expanded compared to the initial one and additionally included the following sequences: axial susceptibility-weighted imaging, sagittal T2-weighted, axial diffusion tensor imaging, and 3D axial T1-weighted sequences.

2.3 miRNA quantification

Plasma samples previously stored at deep-freezing temperatures (−80°C) were used for miRNA quantification. Briefly, plasma miRNAs were isolated using the miRNeasy Serum/Plasma Advanced Kit (Cat. no. 217204; Qiagen, Hilden, Germany). cDNA was synthesized using a miRCURY LNA RT Kit (Cat. no. 339340; Qiagen, Hilden, Germany). Each sample contained miRNA corresponding to 16 μL of plasma sample, 1 x buffer, 1 x miRCURY RT Enzyme mix, synthetic RNA spike-ins. All cDNA samples were synthesized from a master mix of reagents. Samples were incubated for 60 min at 42°C followed by incubation for 5 min at 95°C. The efficiency of miRNA isolation and cDNA synthesis was evaluated by PCR analysis of exogenous synthetic miRNAs (Cat. no. 339390; Qiagen, Hilden, Germany). DNA isolation and cDNA synthesis were successful in all samples.

The quantification of 84 miRNAs was performed using the Serum/Plasma Small Focus miRCURY LNA Panel (cat. no. 339325; Qiagen), which comprised 84 human miRNAs curated by the manufacturer based on their highest expression in serum and plasma samples obtained from healthy individuals and those afflicted with various conditions. The LC96 Roche platform was used. The same panel was used to determine the occurrence of hemolysis. The determinant that indicated an increased risk of hemolysis exceeded the value of dCq = 7, which was the difference between the value of miR-23a, which is constant in the blood, and the value of miR-451a, which is specific to red blood cells. The miRNA analyses were performed at the Department of Molecular Biotechnology and Microbiology of Gdansk University of Technology.

2.4 Statistical analysis and results exploration

Normalized miRNA expression was calculated as 2-(Cq(miRNA) - Cq(mean of references)). The fold change was calculated for each miRNA as the quotient of the mean normalized expression in the HT and non-HT groups.

The Yuen–Welch t-test permutation was applied to compare the trimmed means of the two study groups. The number of permutations was set to 1,000 with a random seed of 123 in the R program. A 10%-trimmed mean value was established. To address the issue of multiple hypothesis testing, p-values were adjusted using the Benjamini-Hochberg correction method. All statistical analyses were performed using R software version 3.6.3.

A multifactor analysis allowed us to control for confounding factors. However, the sample size in the present study was too small to achieve statistical validity; (24) therefore, we limited our single-factor analysis.

To assess miRNA tissue/organ specificity, we utilized the TissueAtlas (24). Subsequently, the significant miRNAs were subjected to functional analysis using DIANA mirPath v.3 (25). To identify miRNA target genes, we selected TarBase v.7.0, a database of experimentally validated miRNA-gene interactions (26). The tool was set to annotate target genes to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, with results merging set to a pathway union (27).

3 Results

3.1 Patients’ characteristic

Of 75 patients with acute ischemic stroke treated with intravenous thrombolysis, 25 were excluded because of incomplete data. Among the remaining 50 patients, HT was detected in 10. All patients with HT experienced stroke in the anterior circulation territory and none had lacunar strokes. As this was only the pilot phase of the study, HT was not differentiated into symptomatic or asymptomatic. Ten patients with anterior circulation and non-lacunar stroke who developed HT were matched with 10 patients without HT, as described in the Methods section. The clinical and demographic characteristics of the patients selected for further investigation are presented in Table 1. The complete dataset is shown in Supplementary Table S1.

Table 1
www.frontiersin.org

Table 1. Demographical and clinical characteristics of patients selected for miRNA analysis.

3.2 miRNA isolation and cDNA synthesis

To estimate the RNA isolation efficiency, UniSp2, UniSp4, and UniSp5 spike-ins were added to the plasma samples, which allowed for the comparison of RNA isolation between samples after RT-qPCR. RNA isolation efficiency remained stable across all samples, with UniSp2 presenting a minimum quantification cycle (Cq) value of 22.01 and a maximum Cq value of 19.47 for samples 108 and 160, respectively. In addition, UniSp4 (present at a 100-fold lower concentration than UniSp2) and UniSp5 (present at a 100-fold lower concentration than UniSp4) were detected in all samples. Similarly, cDNA synthesis efficiency was assessed using the UniSp6 spike-in, which was amplified with the min Cq value of 19.11 and max Cq value of 18.56 among all samples.

3.3 Hemolysis in samples

Based on the analysis of the differences in the Cq values of miR-23a and miR-451a, it was noticed that in six samples, the probability of hemolysis increased (the dCq value was higher than 7). In four samples we observed the hemolysis marker at values between 7.26 and 7.76 and in two samples at values of 8.43 and 8.68 (the exact values of dCq for each sample are presented in the Figure 1). However, since low levels of hemolysis are difficult or impossible to eliminate in clinical practice and we cannot determine the degree of hemolysis based on these results, while our samples were not highly hemolyzed, we decided to continue the research using these samples.

Figure 1
www.frontiersin.org

Figure 1. Graph showing the difference in the Cq values of miR-23a and miR-451a (dCq) in the analyzed samples. The red line indicates the value of dCq above which there is an increased risk of hemolysis.

3.4 miRNA differential expression

The miRNAs hsa-miR-30c-5p, hsa-miR-103a-3p, and hsa-miR-23a-3p were selected as reference miRNAs to calculate the relative expression of plasma miRNAs. Normalized results of miRNA quantification are provided in Supplementary Table S2. When the mean Cq of the reference miRNAs was calculated, sample no. 47 was found to be an outlier. The mean Cq value for the reference miRNAs was 34.67, while the lowest mean Cq value among the other samples was 30.03. Thus, the levels of miRNAs in this sample were notably lower, which resulted in the lack of detection of several plasma miRNAs. This led us to exclude sample 47 from further analyses.

Of the 84 assessed miRNAs, 12 were excluded from further analyses due to missing data in more than 50% of the samples and one was excluded due to non-specific amplification. In our study, “missing data,” referred to situations in which the concentration of genetic material was so low that we were unable to detect it. Following normalization of miRNA levels, differential expression was calculated. The results are shown in Figure 2. Thirteen miRNAs reached a Permutation Yuen-Welch t-test p-value <0.10, with five having a p-value <0.05. However, after Benjamini and Hochberg correction for multiple testing problems, none of the analyzed miRNAs had a p-value <0.05.

Figure 2
www.frontiersin.org

Figure 2. Plot of log2(FC) against –log10(p-value) for all analyzed miRNAs. The values presented on the plot were calculated using the p-values obtained after performing Permutation Yuen-Welch t test. A blue line represents the level of p-value = 0.1. Red lines indicate the range of FC < 0, 67 or > 1, 5.

Regarding the magnitude of difference between the HT and non-HT groups, miRNAs with a fold change (FC) in relative expression (FC) <0.67 were considered downregulated in HT, while miRNAs with FC >1.5 were deemed upregulated. Owing to the nature of the pilot study, the level of significance was set to 0.10. Using the commonly encountered value of 0.05 could potentially result in the pilot study being underpowered (28). The results of the statistical analysis of all miRNAs with a complete set of data are provided in Supplementary Table S3. Sources for the analysis are avaliable in online repository: https://github.com/barbara-kolodziej-gumed/Circulating-miRNA-profiles-and-the-risk-of-HT-after-rtPA.

The analysis included two previously described stroke biomarkers, hsa-miR-124-3p and hsa-miR-17-5p. The concentration of hsa-miR-124-3p was so low in 12 samples that quantitative determination was not possible. For hsa-miR-17-5p, no statistically significant differences were found between the HT and non-HT groups.

MiRNAs with p-values <0.1 and FC <0.67 or > 1.5 are listed in Table 2. Specifically, hsa-miR-133b, hsa-miR-376c-3p, hsa-miR-133a-3p, and hsa-miR-1-3p were downregulated, whereas hsa-miR-17-3p, hsa-miR-296-5p and hsa-miR-7-5p were upregulated in the ICH group. The distributions of their values in the HT and non-HT groups are shown in Figure 3.

Table 2
www.frontiersin.org

Table 2. miRNAs with p < 0.05 and FC >1.5 or < 0.67.

Figure 3
www.frontiersin.org

Figure 3. Box plots of miRNA expression in HT and non-HT group. Only results with p < 0.1 and FC >1.5 or < 0.67 are presented.

Next, the TissueAtlas was used to determine whether these seven miRNAs exhibited specificity to any organ or tissue. Specifically, miR-133b, miR-1-3p, and miR133a-5p appeared to be predominantly present in muscles and myocardium, miR-7-5p was found to be abundant in the brain, miR-376c-3p in the brain, bone, and tunica albuginea, and miR-17-3p in muscle and thyroid tissues, whereas miR-296-5p levels were uniform across all the investigated tissues (Supplementary Table S4).

Subsequently, the DIANA miR path tool was used to assess the potential functional roles of the seven miRNAs. The KEGG pathways significantly associated with the investigated miRNAs are shown in Figure 4. These pathways include those associated with cell adhesion molecules, fatty acid synthesis, and metabolism, as well as those associated with signal transduction, amino acid metabolism, and pathways explored in carcinogenesis.

Figure 4
www.frontiersin.org

Figure 4. KEGG pathways annotated to the 7 investigated miRNAs.

4 Discussion

This pilot study highlighted the differences in circulating miRNA profiles between patients who developed HT and those who did not. We hypothesized that miRNA profiling and other omic approaches may increase the predictive power of thrombolysis-associated sICH in acute ischemic stroke; however, this requires further investigation. This study was a free hypothesis-approach pilot study based on previous work from our group, in which we showed that applying mathematical analyses to extract information from raw big data may result in the identification of new pathophysiological pathways and may complete standard design work.

Despite the low number of patients investigated in this study and the restricted miRNA panel, to the best of our knowledge, this study still reveals the broadest panel of miRNAs that might play a role in hemorrhagic transformation associated with the thrombolytic treatment of acute ischemic stroke.

The miRNAs miRNA-1, −133a and −133b, which were downregulated in HT patients, are known to be present in abundance in both skeletal and cardiac muscle (29). They have further been reported to be elevated in the circulation following a myocardial infarction (30). Notably, miR-133b alone has been found to be a risk factor for cardiovascular disease (31). Another miRNA downregulated in HT, miR-376c, has been extensively studied as a cancer modulator, with evidence supporting a role in both cancer suppression and progression (32). In neonatal hypoxic–ischemic encephalopathy, miR-376c demonstrates protective effects against oxygen–glucose deprivation-induced cell injury (33). One of the upregulated miRNAs, miR-296-5p, has previously been found to be present at high levels in deep venous thrombosis patients. Interestingly, in a mouse model, its elevation inhibited deep venous thrombosis by suppressing S100A4 expression, which influences the release of prothrombosis-related factors (34). Therefore, the observed increase in the level of this miRNA in patients with HT after intravenous thrombolysis in our study is of particular interest. MiR-17-3p is a member of the miR-17/92 cluster and plays a role in cell cycle regulation, proliferation, and apoptosis. It is frequently dysregulated in cardiovascular, immune, and neurodegenerative diseases (35). Notably, in a mouse model, miR-17-3p overexpression suppressed cardiomyocyte apoptosis induced by oxygen–glucose deprivation and reperfusion (36).

The last upregulated miRNA, miR-7-5p, is highly expressed in the brain and is involved in cerebral cortex development. It has also been shown to play a role in neurodegenerative diseases, neuroinflammation, and mental disorders (37). In rodent models, miR-7-5p expression was found to be downregulated in the brain following focal ischemia. Interestingly, administration of miR-7 mimic resulted in smaller infarct volume and better functional recovery, supposedly by repressing α-synuclein (38). Among the imRNAs excluded from the differential expression analysis due to low signal intensity, hsa-miR-124-3p was identified as potentially important. This miRNA is downregulated in ischemic stroke (39), which could explain its undetectable expression in our samples.

The detected trends in miRNA expression aligned with patterns observed in various cardiovascular diseases and the pathophysiological mechanisms associated with ischemic stroke. Implementing our proposed methodology in a large-scale study could facilitate the identification of significant biomarkers indicative of the risk of hemorrhagic transformation. Broadening the spectrum of analyzed miRNAs through either an expanded qPCR panel or small RNA sequencing would further enrich the study and potentially deepen our understanding of the pathophysiology of hemorrhagic transformation.

4.1 Limitations

In this miRNA analysis, owing to the small sample size, a multiple regression model that included the baseline characteristics of the patients used in the matching process was not calculated. The small size of the recruited study group resulted in none of the adjusted p-values being <0.1. Therefore, based on the obtained results, it is not possible to definitively determine which miRNAs may serve as markers of HT after rtPA treatment for stroke. This article is a pilot study and primarily serves as a description of miRNA measurement methods in clinical practice. It should also be noted that the quality of the results may also be influenced by hemolysis, which was observed in six samples. However, in clinical practice, the risks of suboptimal collection and handling of biological materials must be considered. Therefore, these samples were included in our analysis.

5 Conclusion

The quantitative determination of miRNA expression in HT and non-HT patients serves as a complement to our previous analysis (16) and may potentially expand the predictive power of currently used calculators for thrombolysis-associated HT in acute ischemic stroke based on clinical and neuroimaging data; however, further studies are needed to verify this idea. We found differences in the blood concentrations of seven miRNAs, and further analysis of the available miRNA databases revealed that miR-296-5p was most directly linked to the occurrence of HT. However, owing to the small sample size, these results may only be considered as a pilot for any major study. Overall, this study shows that analyzing datasets for a sufficiently large number of patients using a hypothesis-free approach could be used to identify complex pathophysiological associations that may be unachievable in standard-design studies.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

The studies involving humans were approved by The Independent Bioethics Committee for Scientific Research at the Medical University of Gdańsk, Poland. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MS: Formal analysis, Investigation, Project administration, Resources, Visualization, Writing – original draft. AW: Formal analysis, Writing – original draft. PS: Investigation, Writing – review & editing. BKo: Project administration, Visualization, Writing – original draft. AS-S: Writing – original draft, Formal analysis, Investigation, Methodology. BKa: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing. BJ: Writing – review & editing, Data curation.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was Funded by Siemens Helthineers (research grant for the Medical University of Gdansk, PI—BKa, 2018; Contract Number C00229230), by the Excellence Initiative – Research University program at the Medical University of Gdansk and by the Medical University of Gdansk.

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

SUPPLEMENTARY TABLE S1 | Clinical and neuroimaging data of 20 stroke patients with and without hemorrhagic transformation.

SUPPLEMENTARY TABLE S2 | Normalized results of miRNA quantification.

SUPPLEMENTARY TABLE S3 | MicroRNA comparison between study groups.

SUPPLEMENTARY TABLE S4 | Expression of selected miRNAs in various organs according to the TissueAtlas.

References

1. Feigin, VL, Stark, BA, Johnson, CO, Roth, GA, Bisignano, C, Abady, GG, 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. doi: 10.1016/S1474-4422(21)00252-0

PubMed Abstract | Crossref Full Text | Google Scholar

2. Kuriakose, D, and Xiao, Z. Pathophysiology and treatment of stroke: present status and future perspectives. Int J Mol Sci. (2020) 21:1–24. doi: 10.3390/ijms21207609

Crossref Full Text | Google Scholar

3. Bivard, A, Lin, L, and Parsons, MW. Review of stroke Thrombolytics. J Stroke. (2013) 15:90. doi: 10.5853/jos.2013.15.2.90

Crossref Full Text | Google Scholar

4. Van Wijngaarden, JDH, Dirks, M, Huijsman, R, Niessen, LW, Fabbricotti, IN, and Dippel, DWJ. Hospital rates of thrombolysis for acute ischemic stroke: the influence of organizational culture. Stroke. (2009) 40:3390–2. doi: 10.1161/STROKEAHA.109.559492

Crossref Full Text | Google Scholar

5. Spronk, E, Sykes, G, Falcione, S, Munsterman, D, Joy, T, Kamtchum-Tatuene, J, et al. Hemorrhagic transformation in ischemic stroke and the role of inflammation. Front Neurol. (2021) 12:661955. doi: 10.3389/fneur.2021.661955

PubMed Abstract | Crossref Full Text | Google Scholar

6. Fiorelli, M, Bastianello, S, Von Kummer, R, Del Zoppo, GJ, Larrue, V, Lesaffre, E, et al. Hemorrhagic transformation within 36 hours of a cerebral infarct relationships with early clinical deterioration and 3-month outcome in the European cooperative acute stroke study I (ECASS I) cohort. Stroke. (1999) 30:2280–4. doi: 10.1161/01.STR.30.11.2280

Crossref Full Text | Google Scholar

7. Mazya, M, Egido, JA, Ford, GA, Lees, KR, Mikulik, R, Toni, D, et al. Predicting the risk of symptomatic intracerebral hemorrhage in ischemic stroke treated with intravenous alteplase: safe implementation of treatments in stroke (SITS) symptomatic intracerebral hemorrhage risk score. Stroke. (2012) 43:1524–31. doi: 10.1161/STROKEAHA.111.644815

Crossref Full Text | Google Scholar

8. Hong, JM, Kim, DS, and Kim, M. Hemorrhagic transformation after ischemic stroke: mechanisms and management. Front Neurol. (2021) 12:703258. doi: 10.3389/fneur.2021.703258

PubMed Abstract | Crossref Full Text | Google Scholar

9. Charbonnier, G, Bonnet, L, Biondi, A, and Moulin, T. Intracranial bleeding after reperfusion therapy in acute ischemic stroke. Front Neurol. (2021) 11:629920. doi: 10.3389/fneur.2020.629920

PubMed Abstract | Crossref Full Text | Google Scholar

10. Van Kranendonk, KR, Treurniet, KM, Boers, AMM, Berkhemer, OA, Van Den Berg, LA, Chalos, V, et al. Clinical and imaging markers associated with hemorrhagic transformation in patients with acute ischemic stroke. Stroke. (2019) 50:2037–43. doi: 10.1161/STROKEAHA.118.024255

PubMed Abstract | Crossref Full Text | Google Scholar

11. Hao, Y, Zhou, H, Pan, C, Xie, G, Hu, J, Zhang, B, et al. Prediction factors and clinical significance of different types of hemorrhagic transformation after intravenous thrombolysis. Eur J Med Res. (2023) 28:509. doi: 10.1186/s40001-023-01503-x

PubMed Abstract | Crossref Full Text | Google Scholar

12. Whiteley, WN, Slot, KB, Fernandes, P, Sandercock, P, and Wardlaw, J. Risk factors for intracranial hemorrhage in acute ischemic stroke patients treated with recombinant tissue plasminogen activator. Stroke. (2012) 43:2904–9. doi: 10.1161/STROKEAHA.112.665331

Crossref Full Text | Google Scholar

13. Krishnamoorthy, S, Singh, G, Jose, KJ, Soman, B, Foerch, C, Kimberly, WT, et al. Biomarkers in the prediction of hemorrhagic transformation in acute stroke: a systematic review and Meta-analysis. Cerebrovasc Dis. (2022) 51:235–47. doi: 10.1159/000518570

PubMed Abstract | Crossref Full Text | Google Scholar

14. Liu, J, Wang, Y, Jin, Y, Guo, W, Song, Q, Wei, C, et al. Prediction of hemorrhagic transformation after ischemic stroke: development and validation study of a novel multi-biomarker model. Front Aging Neurosci. (2021) 13:667934. doi: 10.3389/fnagi.2021.667934

Crossref Full Text | Google Scholar

15. Xu, Y, Li, X, Wu, D, Zhang, Z, and Jiang, A. Machine learning-based model for prediction of hemorrhage transformation in acute ischemic stroke after Alteplase. Front Neurol. (2022) 13:897903. doi: 10.3389/fneur.2022.897903

Crossref Full Text | Google Scholar

16. Karaszewski, B, Gójska-Grymajło, A, Czaplewska, P, Jabłoński, B, Lewandowska, AE, Ossowska, D, et al. SWATH-MS for prospective identification of protein blood biomarkers of rtPA-associated intracranial hemorrhage in acute ischemic stroke: a pilot study. Sci Rep. (2021) 11:18765. doi: 10.1038/s41598-021-97710-9

PubMed Abstract | Crossref Full Text | Google Scholar

17. O’Brien, J, Hayder, H, Zayed, Y, and Peng, C. Overview of microRNA biogenesis, mechanisms of actions, and circulation. Front Endocrinol. (2018) 9:402. doi: 10.3389/fendo.2018.00402

PubMed Abstract | Crossref Full Text | Google Scholar

18. Zhang, C, Sun, C, Zhao, Y, Wang, Q, Guo, J, Ye, B, et al. Overview of MicroRNAs as diagnostic and prognostic biomarkers for high-incidence cancers in 2021. Int J Mol Sci. (2022) 23:11389. doi: 10.3390/ijms231911389

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wang, Y, Su, X, Leung, GHD, Ren, B, Zhang, Q, Xiong, Z, et al. Circulating microRNAs as diagnostic biomarkers for ischemic stroke: evidence from comprehensive analysis and real-world validation. Int J Med Sci. (2023) 20:1009–23. doi: 10.7150/ijms.83963

PubMed Abstract | Crossref Full Text | Google Scholar

20. He, XW, Shi, YH, Liu, YS, Li, GF, Zhao, R, Hu, Y, et al. Increased plasma levels of miR-124-3p, miR-125b-5p and miR-192-5p are associated with outcomes in acute ischaemic stroke patients receiving thrombolysis. Atherosclerosis. (2019) 289:36–43. doi: 10.1016/j.atherosclerosis.2019.08.002

PubMed Abstract | Crossref Full Text | Google Scholar

21. Zheng, L, Xiong, Y, Liu, J, Yang, X, Wang, L, Zhang, S, et al. MMP-9-related microRNAs as prognostic markers for hemorrhagic transformation in Cardioembolic stroke patients. Front Neurol. (2019) 10:945. doi: 10.3389/fneur.2019.00945

PubMed Abstract | Crossref Full Text | Google Scholar

22. Han, Z, Li, L, Tao, Z, Wang, R, Zhao, H, Zheng, Y, et al. Neutrophilic noncoding RNAs predict outcomes of acute ischemic stroke patients treated with recombinant tissue plasminogen activator. Front Pharmacol. (2022) 13:1003806. doi: 10.3389/fphar.2022.1003806

PubMed Abstract | Crossref Full Text | Google Scholar

23. Stuart, E, King, G, Imai, K, and Ho, D. MatchIT: nonparametric preprocessing for parametric causal inference. J Stat Softw. (2011) 42:1–28. doi: 10.18637/jss.v042.i08

Crossref Full Text | Google Scholar

24. Harrell, FE . Regression modeling strategies: With applications to linear models, logistic and ordinal regression, and survival analysis. Cham, Switzerland: Springer International Publishing (2015).

Google Scholar

25. Keller, A, Gröger, L, Tschernig, T, Solomon, J, Laham, O, Schaum, N, et al. miRNATissueAtlas2: an update to the human miRNA tissue atlas. Nucleic Acids Res. (2022) 50:D211–21. doi: 10.1093/nar/gkab808

PubMed Abstract | Crossref Full Text | Google Scholar

26. Vlachos, IS, Zagganas, K, Paraskevopoulou, MD, Georgakilas, G, Karagkouni, D, Vergoulis, T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. (2015) 43:W460–6. doi: 10.1093/nar/gkv403

PubMed Abstract | Crossref Full Text | Google Scholar

27. Vlachos, IS, Paraskevopoulou, MD, Karagkouni, D, Georgakilas, G, Vergoulis, T, Kanellos, I, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. (2015) 43:D153–9. doi: 10.1093/nar/gku1215

Crossref Full Text | Google Scholar

28. Lee, EC, Whitehead, AL, Jacques, RM, and Julious, SA. The statistical interpretation of pilot trials: should significance thresholds be reconsidered? BMC Med Res Methodol. (2014) 14:41. doi: 10.1186/1471-2288-14-41

PubMed Abstract | Crossref Full Text | Google Scholar

29. Townley-Tilson, WHD, Callis, TE, and Wang, D. MicroRNAs 1, 133, and 206: critical factors of skeletal and cardiac muscle development, function, and disease. Int J Biochem Cell Biol. (2010) 42:1252–5. doi: 10.1016/j.biocel.2009.03.002

PubMed Abstract | Crossref Full Text | Google Scholar

30. Dimmeler, S, and Zeiher, AM. Circulating microRNAs: novel biomarkers for cardiovascular diseases? Eur Heart J. (2010) 31:2705–7. doi: 10.1093/eurheartj/ehq221

Crossref Full Text | Google Scholar

31. Gacoń, J, Badacz, R, Stȩpień, E, Karch, I, Enguita, FJ, Zmudka, K, et al. Diagnostic and prognostic micro-RNAs in ischaemic stroke due to carotid artery stenosis and in acute coronary syndrome: a four-year prospective study. Kardiol Pol. (2018) 76:362–9. doi: 10.5603/KP.a2017.0243

PubMed Abstract | Crossref Full Text | Google Scholar

32. Zhao, F, Zhong, M, Pei, W, Tian, B, and Cai, Y. miR-376c-3p modulates the properties of breast cancer stem cells by targeting RAB2A. Exp Ther Med. (2020) 20:1. doi: 10.3892/etm.2020.9196

Crossref Full Text | Google Scholar

33. Zhang, H, Zhou, J, Zhang, M, Yi, Y, and He, B. Upregulation of miR-376c-3p alleviates oxygen–glucose deprivation-induced cell injury by targeting ING5. Cell Mol Biol Lett. (2019) 24:67. doi: 10.1186/s11658-019-0189-2

PubMed Abstract | Crossref Full Text | Google Scholar

34. Pan, Z, Zhang, Y, Li, C, Yin, Y, Liu, R, Zheng, G, et al. MiR-296-5p ameliorates deep venous thrombosis by inactivating S100A4. Exp Biol Med. (2021) 246:2259–68. doi: 10.1177/15353702211023034

PubMed Abstract | Crossref Full Text | Google Scholar

35. Mogilyansky, E, and Rigoutsos, I. The miR-17/92 cluster: a comprehensive update on its genomics, genetics, functions and increasingly important and numerous roles in health and disease. Cell Death Differ. (2013) 20:1603–14. doi: 10.1038/cdd.2013.125

PubMed Abstract | Crossref Full Text | Google Scholar

36. Shi, J, Bei, Y, Kong, X, Liu, X, Lei, Z, Xu, T, et al. miR-17-3p contributes to exercise-induced cardiac growth and protects against myocardial ischemia-reperfusion injury. Theranostics. (2017) 7:664–76. doi: 10.7150/thno.15162

PubMed Abstract | Crossref Full Text | Google Scholar

37. Zhao, J, Zhou, Y, Guo, M, Yue, D, Chen, C, Liang, G, et al. MicroRNA-7: expression and function in brain physiological and pathological processes. Cell Biosci. (2020) 10:77. doi: 10.1186/s13578-020-00436-w

PubMed Abstract | Crossref Full Text | Google Scholar

38. Kim, TH, Mehta, SL, Morris-Blanco, KC, Chokkalla, AK, Chelluboina, B, Lopez, M, et al. The microRNA miR-7a-5p ameliorates ischemic brain damage by repressing-synuclein. Sci Signal. (2018) 11:eaat4285. doi: 10.1126/scisignal.aat4285

PubMed Abstract | Crossref Full Text | Google Scholar

39. Kanehisa, M, Sato, Y, Kawashima, M, Furumichi, M, and Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. (2016) 44:D457–62. doi: 10.1093/nar/gkv1070

Crossref Full Text | Google Scholar

Keywords: microRNA, miR, ischemic stroke, hemorrhagic transformation, thrombolysis, biomarker

Citation: Stańczak M, Wyszomirski A, Słonimska P, Kołodziej B, Jabłoński B, Stanisławska-Sachadyn A and Karaszewski B (2024) Circulating miRNA profiles and the risk of hemorrhagic transformation after thrombolytic treatment of acute ischemic stroke: a pilot study. Front. Neurol. 15:1399345. doi: 10.3389/fneur.2024.1399345

Received: 11 March 2024; Accepted: 21 May 2024;
Published: 12 June 2024.

Edited by:

Shubham Misra, Yale University, United States

Reviewed by:

Marta Kopańska, University of Rzeszow, Poland
Yuhua Ji, Jinan University, China

Copyright © 2024 Stańczak, Wyszomirski, Słonimska, Kołodziej, Jabłoński, Stanisławska-Sachadyn and Karaszewski. 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: Bartosz Karaszewski, YmFydG9zekBrYXJhc3pld3NraS5vcmc=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.