- 1The Second Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangdong Provincial Hospital of Chinese Medicine, Guangzhou, China
- 2The Second Clinical College of Guangzhou University of Chinese Medicine, Guangzhou, China
Introduction: In this study, we aimed to investigate the association between gut microbiota and high on-treatment platelet reactivity (HTPR) in patients with acute ischemic stroke (AIS).
Methods: We enrolled a total of 48 AIS patients, including 19 HTPR patients and 29 non-high on-treatment platelet reactivity (NHTPR) patients, along with 10 healthy controls. Clinical and laboratory data, as well as stool samples, were collected from all participants. The composition and function of gut microbiota were assessed using 16S rRNA sequencing. Differences in the gut microbiota between the two groups were analyzed, and a diagnostic model based on the gut microbiota was established using random forest model.
Results: HTPR patients exhibited a decreased microbial richness compared to NHTPR patients. Additionally, the relative abundance of unidentified_Clostridia and Ralstonia was lower in HTPR patients. Significant differences in biological functions, such as toxoplasmosis, were observed between the two groups. The combination of Ralstonia, unidentified-Clostridia, Mailhella, Anaerofustis, and Aggregatibacter showed excellent predictive ability for HTPR occurrence (AUC=0.896). When comparing AIS patients with healthy controls, alterations in the microbiota structure were observed in AIS patients, with imbalances in short-chain fatty acid-producing bacteria and pathogenic bacteria. Significant differences in biological functions, such as oxidative phosphorylation, were noted between the two groups. The combination of Alloprevotella, Terrisporobacter, Streptococcus, Proteus, and unidentified_Bacteria exhibited strong predictive power for AIS occurrence (AUC=0.994).
Conclusions: This study is the first to uncover the microbial characteristics of HTPR in AIS patients and demonstrate the predictive potential of specific bacterial combinations for HTPR occurrence.
1 Introduction
Stroke has emerged as the leading cause of disability and mortality among Chinese adults (Hu et al., 2020). In recent years, the incidence and prevalence of acute ischemic stroke (AIS) have continued to rise, even after standardizing for age (Wang et al., 2022). Antiplatelet therapy, particularly aspirin and clopidogrel, is fundamental in secondary prevention of AIS. However, despite persistent adherence to guideline-based secondary stroke prevention, approximately 8.3% of AIS patients still experience recurrent stroke within 12 months (Pan et al., 2021). The reason may be that antiplatelet drugs do not play the expected effect, which is known as high on-treatment platelet reactivity (HTPR) (Fiolaki et al., 2017; Shah et al., 2022). HTPR has been significantly associated with unfavorable outcomes following AIS (Zhou et al., 2022). Therefore, it is important to identify factors associated with HTPR in AIS patients.
Previous study revealed that trimethylamine N-oxide (TMAO), a gut microbial metabolite, increased platelet reactivity by promoting intracellular Ca2+ release, thereby contributing to thrombosis in mice; and specific gut microbial taxa were related to thrombosis in mice, for example, microbial taxa characterized by a high choline diet were significantly positively associated with enhanced thrombosis (Zhu et al., 2016), which highlighted the significant role of gut microbiota in the thrombotic process. In patients with ST-segment elevation myocardial infarction (STEMI) receiving ticagrelor treatment, it was observed that HTPR patients exhibited higher microbial richness and diversity compared to non-high on-treatment platelet reactivity (NHTPR) patients, and fecal transplantation experiments demonstrated that gut microbiota dysbiosis might be an important mechanism for the ticagrelor of HTPR (Zhang et al., 2022). Nevertheless, no previous studies have reported the association between gut microbiota and HTPR in AIS patients, leaving the microbial characteristics of HTPR in AIS patients largely unexplored.
Hence, this study aimed to investigate the microbial characteristics of HTPR in AIS patients, elucidate the potential functions of the gut microbiota in HTPR patients with AIS, and establish a predictive model for HTPR based on the gut microbiota.
2 Materials and methods
2.1 Study patients and data collection
All participants, including AIS patients and healthy controls, were recruited from the Brain Center of Guangdong Provincial Hospital of Traditional Chinese Medicine between February 2022 and January 2023. AIS diagnosis was made based on the Chinese guidelines for the diagnosis and treatment of acute ischemic stroke 2018 (Neurology and Society, 2018). The inclusion criteria for AIS patients were as follows: (1) age ≥ 18 years, (2) administration of clopidogrel (300 mg/day) or aspirin (300 mg/day), or a combination of clopidogrel (300 mg/day) and aspirin (300 mg/day) for at least 12 hours prior to the platelet function test; or administration of aspirin (100 mg/day) or clopidogrel (75 mg/day), or a combination of clopidogrel (75 mg/day) and aspirin (100 mg/day) for at least 5 days before the platelet function test, (3) informed consent obtained. The exclusion criteria for AIS patients were as follows: (1) cardiogenic stroke, (2) use of drugs that affect platelet aggregation, such as traditional non-steroidal antiplatelet drugs, ticlopidine, warfarin, heparin, and low molecular weight heparin, within 1 week before enrollment; use of antibiotics, probiotics, or glucocorticoids within 3 months before enrollment, (3) platelet count < 100×10^9/L or > 450×10^9/L, (4) myelodysplastic syndromes, (5) history of bleeding within the last 3 months or major surgical procedures within the last 4 weeks, (6) history of gastrointestinal diseases such as inflammatory bowel disease or gastrointestinal tumor; or use of alcohol, laxatives, or a history of drug abuse, (7) contraindication to clopidogrel or aspirin, (8) liver or kidney function damage (≥ three times the upper normal limit), (9) pregnancy or lactation, (10) participation in other clinical trials within the last 3 months. The inclusion criteria for healthy controls were as follows: (1) age ≥ 18 years, (2) no history of acute cardio-cerebrovascular events or myocardial infarction within the past 6 months, (3) informed consent obtained. The exclusion criteria for healthy controls were as follows: (1) use of antibiotics, probiotics, or glucocorticoids within 3 months before enrollment, (2) history of gastrointestinal diseases such as inflammatory bowel disease or gastrointestinal tumor, (3) use of alcohol, laxatives, or a history of drug abuse, (4) pregnancy or lactation, (5) participation in other clinical trials within the last 3 months.
Clinical data, including age, gender, medical history, medication history, as well as laboratory data, such as blood lipid levels and blood glucose levels, were collected.
2.2 Sample size calculation
The sample size was calculated using PASS software. According to the previous study, we used PD whole tree diversity index as the evaluating index. Based on Zhang’s study (Zhang et al., 2022), the average values of PD whole tree diversity index in HTPR and NHTPR groups were approximately16 and 13, respectively. The sample sizes of the HRPR and NHTPR groups should be 20 and 30, respectively (Power=0.8, significance level=0.05, mean difference=3, SD=3, allocation ratio=1.5).
2.3 Platelet function assessment
Peripheral venous blood samples were collected in 3.2% trisodium citrate for platelet function assessment. Two methods, namely light transmission aggregometry (LTA) and thromboelastography (TEG), were used in this study to evaluate platelet function.
LTA: The blood samples were centrifuged at 150 g/min for 10 minutes at room temperature to obtain platelet-rich plasma (PRP), while platelet-poor plasma (PPP) was obtained by centrifuging PRP at 3000 g/min for 10 minutes. ADP (150 μmol/L) or arachidonic acid (5 mg/mL) was added to the PRP as inducers. The decrease in plasma turbidity and increase in light transmittance were measured. The platelet aggregator recorded the dynamic change in light transmittance and generated the platelet aggregation curve. In this study, a maximum platelet aggregation rate of ≥50% when ADP (150 μmol/L) was used as the inducer was defined as HTPR on clopidogrel (Jia et al., 2021). Similarly, a maximum platelet aggregation rate of ≥20% when arachidonic acid (5 mg/mL) was used as the inducer was defined as HTPR on aspirin (Jia et al., 2021).
TEG: According to manufacturers’ instructions, TEG used four channels to assess platelet function. The maximum amplitude (MA), which represents the maximal clot strength, was used to evaluate the maximum platelet function stimulated by the inducer. MAADP represented the ADP-induced maximal clot strength, MAAA represented the arachidonic acid-induced maximal clot strength, MAfibrin represented the maximal clot strength induced by fibrin, and MAthrombin represented the maximal clot strength induced by thrombin. The percentage of platelet inhibition induced by ADP/AA was calculated using the formula: ADP/AA%= [(MAthrombin – MAADP/AA)/(MAthrombin - MAfibrin)] ×100%. In this study, ADP% < 30% or MAADP > 47 mm was considered as HTPR on clopidogrel, while AA% < 50% was considered as HTPR on aspirin (Yang et al., 2020).
AIS patients exhibiting HTPR on clopidogrel or aspirin were classified as the HTPR group, while AIS patients without HTPR were classified as the NHTPR group.
2.4 Stool sample collection, DNA extraction and 16S RNA gene amplicon and sequencing
Stool samples (200 mg) were collected and stored in the Fecal Microbial Collection and Preservation Kit (ML-001A, Shenzhen Dayun Gene Technology Co., Ltd.) and saved in an -80°C refrigerator within 3 days (Chen et al., 2022; Luo et al., 2022). Total DNA was extracted following the manufacturer’s instructions using CTAB method. To further ensure the quality of extracted total genomic DNA, the integrity and concentration of the extract were evaluated using 1% agarose gel electrophoresis. The 16S V4 region of the bacterial 16S rRNA gene was amplified using specific primers (515F - 806R). Subsequently, the amplified products were evaluated by 2% agarose gel electrophoresis. PCR products was mixed in equidensity ratios. Then, mixture PCR products was purified with Qiagen Gel Extraction Kit (Qiagen, Germany). Sequencing libraries were generated using the TruSeq DNA PCR-Free Sample Preparation Kit (Illumina, USA) and index codes were added. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. At last, the library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated. Raw Tags were filtered by Qiime (Version 1.9.1, http://qiime.org/scripts/split_libraries_fastq.html) (Caporaso et al., 2010). High-quality sequences were clustered into Operational Taxonomic Units (OTUs) with similarity≥97% by Uparse software (Uparse v7.0.1001, http://drive5.com/uparse/) (Edgar, 2013). Then OTUs were classified into the statistics of each level (phylum, class, order, family, genus, species), and eventually an OTU table was created. Qiime (Version 1.9.1) was used to analyze Alpha diversity and Beta diversity. Linear discriminant analysis (LDA) effect size (LEfSe) was used to find taxa with significant differences (LDA > 2.0).
2.5 Statistical analyses
Statistical analysis was performed using SPSS 18.0 software (Statistical Package for Social Sciences, Chicago, USA). Categorical variables were presented as numbers and percentages and analyzed using the chi-square test. Continuous variables with a normal distribution were expressed as mean ± standard deviation (SD) and analyzed using the Student’s t-test. Continuous variables without a normal distribution were described as median (interquartile range (IQR)) and analyzed using the non-parametric Wilcoxon test. A generalized linear model (GLM) was employed to model the microbiota that was significantly different between the two groups after controlling for possible confounding factors. Spearman correlation analysis was used to identify associations between genera with significant differences in LEfSe analysis and pathways with significant differences in the predictive function analysis.
3 Results
3.1 Baseline characteristics
A total of 463 patients were initially screened, there were 53 AIS patients signing informed consent after strict screening according to inclusion and exclusion criteria, 5 AIS patients without stool sample were excluded, and ultimately, 48 AIS patients were included in the study which were divided into 19 HTPR patients and 29 NHTPR patients. Additionally, 10 healthy controls were enrolled. The flow chart illustrating the participant selection process is presented in Figure 1. Table 1 displays the baseline characteristics of the HTPR group and NHTPR group. The NHTPR group had a higher median BMI compared to the HTPR group (24.22 vs. 21.59, P=0.045). However, no significant differences were observed in other baseline data between the HTPR and NHTPR groups. Table 2 presents the baseline characteristics of AIS patients and healthy controls. AIS patients were older (65.92 ± 10.29 vs. 53.00 ± 4.00, P<0.001) and had higher blood sugar levels (7.01 vs. 5.17, P=0.018) compared to the healthy controls. No statistically significant differences in other baseline data were observed between the AIS and healthy control groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g001.jpg)
Figure 1 Flow chart. AIS: acute ischemic stroke, PLT, platelet count; HTPR, high on-treatment platelet reactivity; NHTPR, non-high on-treatment platelet reactivity.
3.2 Comparison of gut microbiota diversity
To assess microbial richness, α-diversity analyses were performed using the Shannon index, Simpson index, Chao1 index, and Ace index. β-diversity analyses, employing principal coordinate analysis (PCoA), were conducted to compare differences in microbiota structure. The results indicated no statistical differences in the Simpson index and Shannon index between the HTPR and NHTPR groups (Wilcoxon rank-sum test, P=0.558 for the Simpson index, Figure 2A; P=0.476 for the Shannon index, Figure 2B). However, the Chao1 index and Ace index of the NHTPR group were significantly higher than those of the HTPR group (Wilcoxon rank-sum test, P=0.019 for the Chao1 index, Figure 2C; P=0.022 for the Ace index, Figure 2D). These findings suggested that compared to the NHTPR group, the HTPR group exhibited decreased microbial richness. PCoA analysis indicated a similarity in microbial structures between the HTPR and NHTPR groups (P=0.063, Figure 2E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g002.jpg)
Figure 2 Comparison of the microbial communities of the HTPR and NHTPR groups. Box plots depict differences in the gut microbiota diversity indices between the HTPR and NHTPR groups according to the Simpson index (A), Shannon index (B), Chao1 index (C), and Ace index (D); each box plot represents the median, interquartile range, minimum, and maximum values. PCoA shows the gut microbiota composition in the HTPR and NHTPR groups (E); PCoA, principal coordinates analysis.
No statistically significant differences were observed in the Simpson index, Shannon index, Chao1 index, and Ace index between the AIS and healthy control groups (Wilcoxon rank-sum test, P=0.628 for the Simpson index, Figure 3A; P=0.658 for the Shannon index, Figure 3B; P=0.879 for the Chao1 index, Figure 3C; P=0.723 for the Ace index, Figure 3D), indicating similar microbial richness between the AIS and healthy control groups. PCoA analysis demonstrated a significant difference in microbial structure between the AIS and healthy control groups (P<0.001, Figure 3E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g003.jpg)
Figure 3 Comparison of the microbial communities of the AIS and healthy control groups. Box plots depict differences in the gut microbiota diversity indices between the AIS and healthy control groups according to the Simpson index (A), Shannon index (B), Chao1 index (C), and Ace index (D); each box plot represents the median, interquartile range, minimum, and maximum values. PCoA shows the gut microbiota composition in the AIS and healthy control groups (E); PCoA, principal coordinates analysis.
3.3 LEfSe analysis
LEfSe analysis was employed to compare the gut microbiota composition among different groups. The LDA score cutoff of 2.0 was applied to identify taxonomic differences of significance between the groups. The distribution of LDA scores is depicted in Figures 4A, 5A, while the cladogram based on the LEfSe method is presented in Figure 4B and Figure 5B. Here, we considered the differences at the genus level. Figure 4A reveals five bacterial species at the genus level that differ between HTPR patients and NHTPR patients, namely Floricoccus, Ralstonia, Bombella, unidentified-Clostridia, Aggregatibacter. After controlling for potential confounding factors (such as BMI), GLM was employed to assess the significant differences in the identified genera between the two groups. Notably, unidentified_Clostridia and Ralstonia exhibited substantial differences (P < 0.05), suggesting that these genera were associated with HTPR (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g004.jpg)
Figure 4 LEfSe analysis revealed taxonomic differences of gut microbiota in the HTPR and NHTPR groups. (A) The bar plots showed significant differences in gut microbiota between the HTPR group and NHTPR group. The LDA scores (log10) >2 and p < 0.05 are listed. (B) Cladogram indicated the phylogenetic distribution of the gut microbiota associated with the HTPR and NHTPR groups. Red indicates communities with increased relative abundance in the HTPR group, green indicates communities with increased relative abundance in the NHTPR group, and yellow nodes indicate microbial communities that do not play an important role in both groups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g005.jpg)
Figure 5 LEfSe analysis revealed taxonomic differences of gut microbiota in the AIS and healthy control groups. (A) The bar plots showed significant differences in gut microbiota between the AIS group and healthy control group. The LDA scores (log10) >2 and p < 0.05 are listed. (B) Cladogram indicated the phylogenetic distribution of the gut microbiota associated with the AIS and healthy control groups. Red indicates communities with increased relative abundance in the AIS group, green indicates communities with increased relative abundance in the healthy control group, and yellow nodes indicate microbial communities that do not play an important role in both groups.
Figure 5A displays 28 bacterial species at the genus level that demonstrate differences between the AIS and healthy control groups. Following adjustments for potential confounders (such as age and blood glucose), the primary distinctions between the two groups were linked to Parabacteroides, Streptococcus, Blautia, Lachnoclostridium, and unidentified_Bacteria (P < 0.05), indicating that these genera were related to AIS (Table 4).
3.4 Predictive function analysis
To explore the functional characteristics of the gut microbiota in HTPR patients, we employed the PICRUSt2 analysis method, utilizing the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and 16S rRNA gene sequence data for functional prediction. Figure 6A presents the main significantly different metabolic pathways at level 2 between the HTPR and NHTPR groups. Compared to the NHTPR group, the HTPR group exhibited reduced relative abundance in nucleotide metabolism, infectious disease-viral, and sensory system pathways. At level 3, there were 47 metabolic pathways displaying differences between the HTPR and NHTPR groups. Figure 6B illustrates the top 30 metabolic pathways with significant differences, all of which exhibited lower relative abundance in the HTPR group compared to the NHTPR group. Notably, insect hormone biosynthesis displayed the highest relative abundance, followed by toxoplasmosis. Spearman’s correlation analysis was performed between the pathways with significant differences, the top 10 relative abundance at level 3, and the genera exhibiting significant differences. Ralstonia demonstrated a moderate correlation with the above 10 pathways, with the highest correlation observed with toxoplasmosis. On the other hand, unidentified_Clostridia displayed a weak correlation with the above 10 pathways, with the lowest correlation observed with insect hormone biosynthesis (Figure 6C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g006.jpg)
Figure 6 Functional predictions of the gut microbiota and heatmap of Spearman’ s correlation analysis in the HTPR and NHTPR groups. (A) The main metabolic pathways with significant differences between the HTPR and NHTPR groups at level 2. (B) The top 30 metabolic pathways with significant differences between the HTPR and NHTPR groups at level 3. (C) Heatmap of Spearman’ s correlation analysis between the pathways with significant differences and the top 10 relative abundance at level 3 and the genera with significant differences. *p < 0.05, **p < 0.01.
At level 2, there were 7 metabolic pathways exhibiting significant differences between the AIS and healthy control groups (Figure 7A). Among these, metabolism of cofactors and vitamins displayed the highest relative abundance, followed by energy metabolism. At level 3, there were 69 metabolic pathways showing differences between the AIS and healthy control groups. Figure 7B depicts the top 30 metabolic pathways with significant differences, with pyrimidine metabolism displaying the highest relative abundance, followed by oxidative phosphorylation. Spearman’s correlation analysis was conducted between the pathways with significant differences and the top 10 relative abundance at level 3 and the genera displaying significant differences. Figure 7C reveals that Lachnoclostridium exhibits a positive correlation with oxidative phosphorylation, transcription machinery, and vancomycin resistance, and a negative correlation with glutathione metabolism; similarly, Blautia demonstrats a positive correlation with oxidative phosphorylation and a negative correlation with glutathione metabolism; Parabacteroides displays a positive correlation with protein processing in the endoplasmic reticulum and a negative correlation with riboflavin metabolism; unidentified_Bacteria exhibits a positive correlation with glutathione metabolism and a negative correlation with riboflavin metabolism, oxidative phosphorylation, transcription machinery, and pyrimidine metabolism; Streptococcus displays a positive correlation with glutathione metabolism and a negative correlation with protein processing in the endoplasmic reticulum, oxidative phosphorylation, transcription machinery, nucleotide excision repair, vancomycin resistance, riboflavin metabolism, tuberculosis, and legionellosis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g007.jpg)
Figure 7 Functional predictions of the gut microbiota and heatmap of Spearman’ s correlation analysis in the AIS and healthy control groups. (A) The main metabolic pathways with significant differences between the AIS and healthy control groups at level 2. (B) The top 30 metabolic pathways with significant differences between the AIS and healthy control groups at level 3. (C) Heatmap of Spearman’ s correlation analysis between the pathways with significant differences and the top 10 relative abundance at level 3 and the genera with significant differences. *p < 0.05, **p < 0.01.
3.5 Random forest predictive models
To further investigate the potential of gut microbiota in diagnosing HTPR, we employed random forest to develop a prediction model using the bacterial species with significant importance at the genus level as input. The relative importance of each genus in the predictive model was evaluated using mean decrease Gini. Among the 30 genera analyzed, unidentified-Clostridia ranked first in terms of importance based on mean decrease Gini (Figure 8A). Notably, Ralstonia, unidentified-Clostridia, Mailhella, Anaerofustis, and Aggregatibacter showed significant associations with the occurrence of HTPR, with areas under the receiver operating characteristic (ROC) curve (AUC) of 0.750, 0.716, 0.705, 0.683, and 0.676, respectively. To enhance diagnostic accuracy, we constructed an ROC curve incorporating these five bacterial species, resulting in an AUC of 0.896 (95% confidence interval: 0.806–0.986) (Figure 8B). Hence, gut microbiota holds promise as a potential diagnostic tool for HTPR.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g008.jpg)
Figure 8 The predictive model based on genus-level abundance taxa using a random forest model in the HTPR group. (A) The relative importance of each genus was assessed by the mean decrease gini. (B) ROC curve constructed by random forest in gut microbiota. AUC is shown in the bottom right corner of the figure with 95% confidence intervals in parentheses. ROC, receiver operating characteristic; AUC, area under the ROC curve.
Similarly, we developed a prediction model using random forest and the bacterial species with significant importance at the genus level to predict the occurrence of AIS. Among the 30 genera analyzed, Alloprevotella ranked first in terms of mean decrease Gini (Figure 9A). Alloprevotella, Terrisporobacter, Streptococcus, Proteus, and unidentified_Bacteria were significantly associated with the occurrence of AIS, exhibiting AUCs of 0.928, 0.884, 0.879, 0.858, and 0.856, respectively. To improve diagnostic performance, we constructed a ROC curve incorporating these five bacterial species, yielding an AUC of 0.994 (95% confidence interval: 0.980–1.000) (Figure 9B). These findings suggest that gut microbiota has the potential to serve as a valuable diagnostic tool for AIS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1257317/fcimb-13-1257317-HTML/image_m/fcimb-13-1257317-g009.jpg)
Figure 9 The predictive model based on genus- level abundance taxa using a random forest model in the AIS group. (A) The relative importance of each genus was assessed by the mean decrease gini. (B) ROC curve constructed by random forest in gut microbiota. AUC is shown in the bottom right corner of the figure with 95% confidence intervals in parentheses. ROC, receiver operating characteristic; AUC, area under the ROC curve.
4 Discussion
In our study, we observed alterations in the gut microbiota and functions in patients with HTPR and AIS. Specific combinations of bacterial genera demonstrated predictive potential for HTPR and AIS occurrence. Here are the key findings regarding the gut microbiota characteristics in HTPR and AIS patients: HTPR patients exhibited a decreased microbial richness compared to NHTPR patients; additionally, the relative abundance of unidentified_Clostridia and Ralstonia was lower in HTPR patients; significant differences in biological functions, such as toxoplasmosis, were observed between the two groups; the combination of Ralstonia, unidentified-Clostridia, Mailhella, Anaerofustis, and Aggregatibacter showed excellent predictive ability for HTPR occurrence; when comparing AIS patients with healthy controls, alterations in the microbiota structure were observed in AIS patients, with imbalances in short-chain fatty acid-producing bacteria and pathogenic bacteria; significant differences in biological functions, such as oxidative phosphorylation, were noted between the two groups; the combination of Alloprevotella, Terrisporobacter, Streptococcus, Proteus, and unidentified_Bacteria exhibited strong predictive power for AIS occurrence.
4.1 The microbial richness decreased and the bacterial taxonomic composition changed in HTPR patients
So far, there is only one study reported by Zhang X et al. investigating the relationship between gut microbiota and HTPR in patients with STEMI after ticagrelor treatment (Zhang et al., 2022). There is no published study on the relationship between the gut microbiota and HTPR in AIS patients. Zhang X et al. (Zhang et al., 2022) explored the mechanism of the poor response of ticagrelor in patients with STEMI and found that HTPR patients had increased richness and altered microbiota structure compared with NHTPR patients, while the results of our study showed that the microbial richness of the HTPR patients decreased compared with NHTPR patients, and the microbial structures between the HTPR and NHPTR groups were similar. As we know, the diversity of gut microbiota is affected by many factors such as age, disease, diet, and lifestyle (Lozupone et al., 2012), which may be the reason for the different results of the two studies.
Furthermore, we observed a decrease in the relative abundance of unidentified_Clostridia and Ralstonia in HTPR patients. Unidentified_Clostridia have been reported to inhibit the growth of harmful bacteria, promote the growth of beneficial bacteria, enhance immunity, and support nutrient digestion and absorption (Liu et al., 2021). Therefore, the decrease in unidentified_Clostridia could disrupt the balance between beneficial and harmful bacteria, which might have a negative effect on health. Moreover, the absence of unidentified_Clostridia might influence the digestion and absorption of antiplatelet drugs, potentially reducing their bioavailability and efficacy, which could contribute to the occurrence of HTPR. Ralstonia, a prevalent pro-inflammatory opportunistic pathogen, has been previously found to be more abundant in fecal samples and sigmoid mucosa of Parkinson’s disease patients compared to controls (Keshavarzian et al., 2015). However, the role of Ralstonia in HTPR patients remains unclear and further investigations are warranted to determine its significance.
Functional predictions using PICRUSt2 analysis revealed differences in biological functions, such as toxoplasmosis, between the HTPR and NHTPR groups. Spearman’s correlation analyses indicated associations between Ralstonia and unidentified_Clostridia and specific pathways at level 3. For example, Ralstonia showed a moderate correlation with toxoplasmosis, while unidentified_Clostridia exhibited a weak correlation with insect hormone biosynthesis. These findings suggest that gut microbiota may influence the occurrence and progression of HTPR through these specific pathways. However, further studies are required to elucidate the potential mechanisms by which unidentified_Clostridia and Ralstonia impact HTPR through these pathways.
4.2 The combination of special bacteria could predict the occurrence of HTPR
In the study reported by Zhang et al. (Zhang et al., 2022), the occurrence of HTPR was predicted using random forest with the relative importance order of 21 genera as input, yielding the AUC of 0.807. Similarly, in our study, we successfully distinguished HTPR from NHTPR using specific bacterial genera. The AUCs of Ralstonia, unidentified-Clostridia, Mailhella, Anaerofustis, and Aggregatibacter were 0.750, 0.716, 0.705, 0.683, and 0.676, respectively, and the AUC of the combination of these five bacterial species was 0.896, highlighting the significant role of gut microbiota in HTPR diagnosis. However, it is worth noting that the diagnostic value of specific gut microbiota can vary across different studies. Therefore, large-scale multicenter, cross-ethnic, and cross-regional investigations are necessary to thoroughly evaluate the diagnostic potential of gut microbiota in HTPR.
4.3 The microbial structure and taxonomic composition changed in AIS patients
In our study, we did not observe a statistically significant difference in α-diversity analysis between the AIS and healthy control groups. However, β-diversity analysis revealed significant differences in the microbial structure between the AIS and healthy control groups. The relationship between α and β diversity of gut microbiota and AIS is currently a matter of debate. A previous study has shown that Chao1, observed species, and PD whole tree in IS or TIA group were significantly higher than those in control group, the Shannon index showed a similar trend, but the difference was not statistically significant, the β-diversity analysis showed that the microbial structure of IS or TIA group was significantly different from that of control group (Yin et al., 2015). Li et al. reported that there were no differences in Shannon index, Simpson index, ACE index, Chao1 index, and the bacterial community structures between ischemic stroke and healthy control groups (Li et al., 2019). The diversity of gut microbiota is influenced by various factors, which may be the reason for the differences in diversity analysis among different studies.
Moreover, we observed increased relative abundance of Parabacteroides, Streptococcus, and unidentified_Bacteria, and decreased relative abundance of Lachnoclostridium and Blautia in AIS patients. Parabacteroides is known to be an opportunistic pathogen, and its increased abundance has been associated with sepsis in ICU patients (Agudelo-Ochoa et al., 2020). Streptococcus can cause various infections, including impetigo, pharyngitis, and necrotizing fasciitis (Li et al., 2023). Mulla et al. has reported that Streptococcus may be related to interleukin-6 (IL6), tumor necrosis factor (TNF)-alpha, and other cytokines (Xiong et al., 2023). The increased abundance of Parabacteroides and Streptococcus suggests that opportunistic pathogens may proliferate in AIS patients, consistent with the pattern of post-stroke gut microbiota dysbiosis observed in previous studies (Yin et al., 2015; Sun et al., 2022). It has been reported that unidentified_Bacteria may be involved in acetic oxidation from butyrate (a short-chain fatty acid) metabolism (Yi et al., 2020). Lachnoclostridium and Blautia are known SCFA-producing bacteria (Kang et al., 2021; Zhao et al., 2021), and Blautia has anti-inflammatory effects and contributes to the recovery of intestinal mucosal injury (Zhao et al., 2021). The differential trends observed in these SCFA-producing bacteria among AIS patients warrant further investigation to elucidate their roles in AIS.
PICRUSt2 analysis revealed differences in biological functions between the AIS and healthy control groups. Spearman’s correlation analyses indicated that Parabacteroides, Streptococcus, unidentified_Bacteria, Lachnoclostridium and Blautia were associated with specific pathways at level 3. Future studies are required to validate the potential impact of these specific bacteria on AIS through these pathways.
4.4 The combination of special bacteria could predict the occurrence of AIS
In this study, the occurrence of AIS could be effectively diagnosed by special bacterial combination at the genus level. The AUCs of Alloprevotella, Terrisporobacter, Streptococcus, Proteus, unidentified_Bacteria were 0.928, 0.884, 0.879, 0.858, and 0.856, respectively. Remarkably, the combination of these five bacterial species yielded an impressive AUC of 0.994.
At present, there are some studies focusing on the gut microbiota of AIS or specific AIS populations, such as the study reported by Chen et al. (Chen et al., 2022). Chen’s study was to investigate the characteristic gut microbiota of poor outcome of AIS with hyperlipidemia (POAH) patients. They found that in POAH patients, the relative abundance of Enterococcaceae and Enterococcus were increased, and the relative abundance of Lachnospiraceae, Faecalibacterium, Rothia and Butyricicoccus were decreased. And they established a ROC model based on the combination of characteristic microbiota (Unclassified-f-Lachnospiraceae, Enterococcus, Faecalibacterium, Lachnospiraceae-UCG-010, and norank-f-Lachnospiraceae) to predict POAH (AUC = 0.694, 95% CI 0.618-0.770).
However, our study is the first to investigate the association between the gut microbiota and HTPR in AIS patients. In this study, we found that the relative abundance of unidentified_Clostridia and Ralstonia was lower in HTPR patients, and the combination of Ralstonia, unidentified-Clostridia, Mailhella, Anaerofustis, and Aggregatibacter showed excellent predictive ability for HTPR occurrence (AUC=0.896). And we observed increased relative abundance of Parabacteroides, Streptococcus, and unidentified_Bacteria, and decreased relative abundance of Lachnoclostridium and Blautia in AIS patients, the combination of Alloprevotella, Terrisporobacter, Streptococcus, Proteus, and unidentified_Bacteria exhibited strong predictive power for AIS occurrence (AUC=0.994).
In this study, several limitations should be mentioned. Firstly, this was a small-scale, single-center observational study, which may introduce inherent biases. Secondly, the use of 16S rRNA sequencing to analyze changes in gut microbiota may provide limited information compared to metagenomic sequencing, which offers more comprehensive insights into bacterial species and functions. Despite these limitations, our study provides novel insights into the microbial characteristics of HTPR patients with AIS, enhancing our understanding of the association between gut microbiota and HTPR in AIS patients. These findings also offer a theoretical basis for developing personalized medication strategies for HTPR from the perspective of gut microbiota.
5 Conclusion
In summary, this study is the first to uncover the microbial characteristics of HTPR in AIS patients and demonstrate the predictive potential of specific bacterial combinations for HTPR occurrence.
Data availability statement
The data presented in the study are deposited in the National Center for Biotechnology Information repository, and the accession number is PRJNA995199.
Ethics statement
This study was approved by the Ethics Committee of the Guangdong Provincial Hospital of Traditional Chinese Medicine (Grant No. YE2022-022-01) and conducted in accordance with the principles of the Declaration of Helsinki. Written informed consent was obtained from all participants.
Author contributions
ZL: Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. HO: Data curation, Writing – review & editing. GC: Data curation, Writing – review & editing. XL: Data curation, Writing – review & editing. HC: Data curation, Writing – review & editing. YZhan: Investigation, Writing – review & editing. LP: Investigation, Writing – review & editing. CD: Investigation, Writing – review & editing. ZZ: Methodology, Writing – review & editing. LW: Methodology, Writing – review & editing. HX: Methodology, Writing – review & editing. MZ: Methodology, Supervision, Validation, Writing – review & editing. YZhao: Methodology, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the “State Key Laboratory of Dampness Syndrome of Chinese Medicine (SZ2021ZZ46 and SZ2022KF22)”, “Guangzhou Science and Technology Project (202002020034 and 202201020506)”, “Guangdong Province Key Laboratory of Chinese Medicine Emergency Special Project (2019KT1340)”, “Guangdong Provincial Bureau of Traditional Chinese Medicine research project (20225030)”, “National Administration of Traditional Chinese Medicine (0102016401)”, “NATCM’s Project of High-level Construction of Key TCM Disciplines (zyyzdxk-2023154)”. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Acknowledgments
We sincerely thank all participants who volunteered to participate in this study.
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.
References
Agudelo-Ochoa, G. M., Valdés-Duque, B. E., Giraldo-Giraldo, N. A., Jaillier-Ramírez, A. M., Giraldo-Villa, A., Acevedo-Castaño, I., et al. (2020). Gut microbiota profiles in critically ill patients, potential biomarkers and risk variables for sepsis. Gut. Microbes 12 (1), 1707610. doi: 10.1080/19490976.2019.1707610
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). Qiime allows analysis of high-throughput community sequencing data. Nat. Methods. 7 (5), 335–336. doi: 10.1038/nmeth.f.303
Chen, J., Chi, B., Ma, J., Zhang, J., Gu, Q., Xie, H., et al. (2022). Gut microbiota signature as predictors of adverse outcomes after acute ischemic stroke in patients with hyperlipidemia. Front. Cell. Infect. Microbiol. 12. doi: 10.3389/fcimb.2022.1073113
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods. 10 (10), 996–998. doi: 10.1038/nmeth.2604
Fiolaki, A., Katsanos, A. H., Kyritsis, A. P., Papadaki, S., Kosmidou, M., Moschonas, I. C., et al. (2017). High on treatment platelet reactivity to aspirin and clopidogrel in ischemic stroke: a systematic review and meta-analysis. J. Neurol. Sci. 376, 112–116. doi: 10.1016/j.jns.2017.03.010
Hu, S., Cui, B., Mlynash, M., Zhang, X., Mehta, K. M., Lansberg, M. G. (2020). Stroke epidemiology and stroke policies in China from 1980 to 2017: a systematic review and meta-analysis. Int. J. Stroke. 15 (1), 18–28. doi: 10.1177/1747493019873562
Jia, W., Jia, Q., Zhang, Y., Zhao, X., Wang, Y. (2021). Effect of prediabetes on asprin or clopidogrel resistance in patients with recent ischemic stroke/tia. Neurol. Sci. 42 (7), 2829–2835. doi: 10.1007/s10072-020-04881-w
Kang, L., Li, P., Wang, D., Wang, T., Hao, D., Qu, X. (2021). Alterations in intestinal microbiota diversity, composition, and function in patients with sarcopenia. Sci. Rep. 11 (1), 4628. doi: 10.1038/s41598-021-84031-0
Keshavarzian, A., Green, S. J., Engen, P. A., Voigt, R. M., Naqib, A., Forsyth, C. B., et al. (2015). Colonic bacterial composition in parkinson’s disease. Mov. Disord. 30 (10), 1351–1360. doi: 10.1002/mds.26307
Li, N., Wang, X., Sun, C., Wu, X., Lu, M., Si, Y., et al. (2019). Change of intestinal microbiota in cerebral ischemic stroke patients. BMC Microbiol. 19 (1), 191. doi: 10.1186/s12866-019-1552-1
Li, Z., Liang, H., Hu, Y., Lu, L., Zheng, C., Fan, Y., et al. (2023). Gut bacterial profiles in parkinson’s disease: a systematic review. CNS Neurosci. Ther. 29 (1), 140–157. doi: 10.1111/cns.13990
Liu, Y., Zhou, M., Yang, M., Jin, C., Song, Y., Chen, J., et al. (2021). Pulsatilla chinensis saponins ameliorate inflammation and dss-induced ulcerative colitis in rats by regulating the composition and diversity of intestinal flora. Front. Cell. Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.728929
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K., Knight, R. (2012). Diversity, stability and resilience of the human gut microbiota. Nature 489 (7415), 220–230. doi: 10.1038/nature11550
Luo, J., Chen, Y., Tang, G., Li, Z., Yang, X., Shang, X., et al. (2022). Gut microbiota composition reflects disease progression, severity and outcome, and dysfunctional immune responses in patients with hypertensive intracerebral hemorrhage. Front. Immunol. 13. doi: 10.3389/fimmu.2022.869846
Neurology, C. S. O., Society, C. S. (2018). Chinese guidelines for diagnosis and treatment of acute ischemic stroke 2018. Chin. J. Neurol. 51 (09), 666–682. doi: 10.3760/cma.j.issn.1006-7876.2018.09.004
Pan, Y., Li, Z., Li, J., Jin, A., Lin, J., Jing, J., et al. (2021). Residual risk and its risk factors for ischemic stroke with adherence to guideline-based secondary stroke prevention. J. Stroke. 23 (1), 51–60. doi: 10.5853/jos.2020.03391
Shah, J., Liu, S., Yu, W. (2022). Contemporary antiplatelet therapy for secondary stroke prevention: a narrative review of current literature and guidelines. Stroke. Vasc. Neurol. 7 (5), 406–414. doi: 10.1136/svn-2021-001166
Sun, H., Gu, M., Li, Z., Chen, X., Zhou, J. (2022). Gut microbiota dysbiosis in acute ischemic stroke associated with 3-month unfavorable outcome. Front. Neurol. 12. doi: 10.3389/fneur.2021.799222
Wang, Y., Li, Z., Gu, H., Zhai, Y., Zhou, Q., Jiang, Y., et al. (2022). China stroke statistics: an update on the 2019 report from the national center for healthcare quality management in neurological diseases, China national clinical research center for neurological diseases, the chinese stroke association, national center for chronic and non-communicable disease control and prevention, chinese center for disease control and prevention and institute for global neuroscience and stroke collaborations. Stroke. Vasc. Neurol. 7 (5), 415–450. doi: 10.1136/svn-2021-001374
Xiong, Z., Wang, Q., Pei, S., Zhu, Z. (2023). The causal role of intestinal microbiome in development of pre-eclampsia. Funct. Integr. Genomics 23 (2), 127. doi: 10.1007/s10142-023-01054-8
Yang, Y., Chen, W., Pan, Y., Yan, H., Meng, X., Liu, L., et al. (2020). Effect of ticagrelor versus clopidogrel on platelet reactivity measured by thrombelastography in patients with minor stroke or tia. Aging (Albany. Ny.). 12 (20), 20085–20094. doi: 10.18632/aging.103452
Yi, Y., Wang, H., Chen, Y., Gou, M., Xia, Z., Hu, B., et al. (2020). Identification of novel butyrate- and acetate-oxidizing bacteria in butyrate-fed mesophilic anaerobic chemostats by dna-based stable isotope probing. Microb. Ecol. 79 (2), 285–298. doi: 10.1007/s00248-019-01400-z
Yin, J., Liao, S. X., He, Y., Wang, S., Xia, G. H., Liu, F. T., et al. (2015). Dysbiosis of gut microbiota with reduced trimethylamine-n-oxide level in patients with large-artery atherosclerotic stroke or transient ischemic attack. J. Am. Heart Assoc. 4 (11), e002699. doi: 10.1161/JAHA.115.002699
Zhang, X., Zhang, X., Tong, F., Cai, Y., Zhang, Y., Song, H., et al. (2022). Gut microbiota induces high platelet response in patients with st segment elevation myocardial infarction after ticagrelor treatment. Elife 11, e70240. doi: 10.7554/eLife.70240
Zhao, T., Zhan, L., Zhou, W., Chen, W., Luo, J., Zhang, L., et al. (2021). The effects of erchen decoction on gut microbiota and lipid metabolism disorders in zucker diabetic fatty rats. Front. Pharmacol. 12. doi: 10.3389/fphar.2021.647529
Zhou, K., Yu, S., Li, J., Tan, Y., Xing, S., Chen, Y., et al. (2022). High on-treatment platelet reactivity is associated with poor outcomes after ischemic stroke: a meta-analysis. Acta Neurol. Scand. 146 (3), 205–224. doi: 10.1111/ane.13655
Keywords: acute ischemic stroke, high on-treatment platelet reactivity, gut microbiota, 16S rRNA sequencing, random forest model
Citation: Lou Z, Ouyang H, Chen G, Li X, Chen H, Zhan Y, Peng L, Du C, Zheng Z, Wen L, Xu H, Zhao M and Zhao Y (2024) Gut microbiota as predictors of the occurrence of high on-treatment platelet reactivity in acute ischemic stroke patients. Front. Cell. Infect. Microbiol. 13:1257317. doi: 10.3389/fcimb.2023.1257317
Received: 12 July 2023; Accepted: 07 December 2023;
Published: 04 January 2024.
Edited by:
Julio Villena, CONICET Centro de Referencia para Lactobacilos (CERELA), ArgentinaReviewed by:
Michiyo Matsumoto-Nakano, Okayama University, JapanJustyna Agier, Medical University of Lodz, Poland
Copyright © 2024 Lou, Ouyang, Chen, Li, Chen, Zhan, Peng, Du, Zheng, Wen, Xu, Zhao and Zhao. 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: Yuanqi Zhao, zyq2022@gzucm.edu.cn; Min Zhao, 376479808@qq.com