- 1Faculty of Applied Science, Universiti Teknologi MARA Shah Alam, Shah Alam, Selangor, Malaysia
- 2Department of Food Research, Zhejiang Business College, Hangzhou, Zhejiang, China
- 3Zhejiang Institute for Food and Drug Control, Hangzhou, China
- 4Integrative Pharmacogenomics Institute (iPROMISE), Puncak Alam, Bandar Puncak Alam, Selangor, Malaysia
Background: Extracellular vesicles (EVs) of probiotics are an important pathway for probiotics to exert probiotic effects. EVs of Clostridium butyricum (C. butyricum) has great potential in immunotherapy. The aim of this study is to investigate the regulatory effect of C. butyricum EVs on the pyroptosis-related gene network in LPS-stimulated macrophages and explore its potential molecular mechanisms.
Methods: Differential centrifugation was used to isolate C. butyricum EVs, and EVs were characterized by transmission electron microscopy, NTA, and ZETA potential. PI staining was employed to confirm the stimulatory effect of LPS on RAW264.7 cells. Transcriptome sequencing and qPCR were performed to validate the regulatory effects of EVs on pyroptosis-related genes (PRGs) in LPS-stimulated RAW264.7 cells.
Results: C. butyricum EVs exhibited a bilayer membrane structure with a median particle size of 144.4 ± 85.7 nm and a ZETA potential of -37.56 ± 1.08 mV. PI staining showed that LPS induced cell death in RAW264.7 macrophages, while EV treatment significantly affected the expression of pyroptosis-related genes IL-18 and IL-33 (q < 0.05). Transcriptome analysis revealed that C. butyricum EVs induce inflammation and cytokine secretion in RAW264.7 macrophages through the JAK-STAT signaling pathway and Cytokine receptor interaction. EVs pretreatment (10/20 μg/mL) regulated the expression of 14 PRGs in a concentration-related manner, including significant upregulation of IL-6 and P2rx7 (q < 0.01) and inhibition of ZBP1 and Tnf (q < 0.05). The gene interaction network showed a positive correlation (r=0.82) between IL-6 and P2rx7, and a negative correlation (r=-0.75) between ZBP1 and P2rx7.
Conclusion: C. butyricum EVs remodel the pyroptosis-related regulatory network by modulating key genes such as IL-6, P2rx7, and ZBP1, suggesting their potential role in fine-tuning LPS-induced inflammatory responses in macrophages.
1 Introduction
Butyric acid-producing Clostridium butyricum (C. butyricum) is a type of intestinal probiotic that can regulate the structure of the host’s intestinal microbiota and the body’s immune function, and is of great significance for maintaining intestinal health (1, 2). C. butyricium itself has been shown to significantly affect host metabolism and gut microbiota structure (3). Studies have shown that C. butyricum can produce extracellular vesicles (EVs), which can regulate the host’s cellular immune function through EVs, regulate intestinal homeostasis and improve colitis (4). While C. butyricum EVs have been shown to be non-toxic and have a higher yield than Lactobacillus rhamnosus (5), their precise immunomodulatory mechanisms, particularly their impact on the transcriptional programs of immune cells, remain to be fully characterized. One important step is to define the gene regulatory networks influenced by these EVs. Specifically, the transcriptional mechanisms through which C. butyricum EVs interact with macrophages to reshape the expression of inflammation- and pyroptosis-related genes are not yet clear. Elucidating this gene-level interaction is essential for understanding their immunoregulatory effects and future therapeutic applications.
Macrophages are the core components of the innate immune system (6), and their functional activation is central to regulating inflammatory responses. A key aspect of macrophage activation is the secretion of pro-inflammatory cytokines such as TNF-α and IL-1β (7). However, excessive immune responses can lead to tissue damage and negatively impact the body. Excessive secretion of inflammatory cytokines can cause pyroptosis, necroptosis, and even cell apoptosis. Studies have shown that IFN can induce ZBP1, causing inflammatory cell death in macrophages (8). Inflammatory cell death of macrophages, especially pyroptosis, plays an important role in the occurrence of various diseases (9, 10).
Studies have shown that the extracellular vesicles of intestinal commensal bacteria can regulate cellular inflammatory responses, affect the expression of inflammatory cytokines, and also affect LPS induced inflammatory responses (11–13). Some intestinal bacterial extracellular vesicles weaken LPS induced NLRP3 inflammasome activation by affecting the NLRP3-Caspase 1-ASC signaling pathway, thereby reducing inflammatory responses (14). Apoptosis is closely related to inflammatory response, and NLRP3 inflammasome is an important link in cell apoptosis (15). Regulating the NLRP3/Caspase 1 signaling pathway is an effective method to inhibit cell apoptosis (16). EVs from beneficial intestinal bacteria, such as C. butyricum, may possess the capacity to influence these processes. Preliminary evidence suggests that C. butyricum EVs can regulate inflammatory cytokines and mitigate LPS-induced inflammatory responses, warranting further investigation into their potential role in modulating inflammatory cell death. This study aims to explore the potential mechanism of the immunomodulatory effect of C. butyricum EV on RAW264.7 macrophages, and investigate their regulatory impact on the pyroptosis-related gene network in LPS-stimulated RAW264.7 macrophages.
2 Materials and methods
2.1 Materials
The following materials were utilized: RAW264.7 macrophages (CTCC-001-0048) sourced from Zhejiang Meisen Cell Technology Co., Ltd. (China); C. butyricum strain CBM588 procured from Miyarisan (Japan); The C. butyricum CBM588 (MIYAIRI 588) is a well-documented and commercially used probiotic strain, employed in commercial pharmaceutical preparations such as MIYA-BM TABLETS. Its selection is based on its established safety profile and proven efficacy in modulating host biology, which lends significant relevance to our investigation into host-inflammatory responses. Reinforced Clostridial Medium (RCM) from Haibo Biotechnology Co., Ltd. (China); and Eagle’s Minimum Essential Medium (EMEM, 1X) from Gibco (USA). Fetal bovine serum (FBS; heat-inactivated) and lipopolysaccharide (LPS) were supplied by Shanghai Datsil Bio-Tech Co., Ltd. (China) and Beyotime (China), respectively. Molecular reagents included a reverse transcription kit (Vazyme, China), AceQ qPCR SYBR Green Master Mix (Vazyme, China), and a total RNA extraction kit (TIANGEN, China). DEPC-treated water was acquired from Sangon Biotech (China). DL-Dithiothreitol (DTT) and Tris were purchased from Solarbio (China). Trypsin was obtained from Promega (USA). Iodoacetamide (IAM) was sourced from Aladdin (China). Phenylmethanesulfonyl fluoride (PMSF) was bought from Xiya Reagent (China). Tetraethylammonium bromide (TEAB) and Urea were procured from Sigma (USA). Ethylene Diamine Tetraacetic Acid (EDTA), Xylene brilliant cyaninG (G250), Sodium dodecyl sulfate (SDS), Thiourea, and acetone were supplied by Sinopharm (China). PCR primers and the BCA Protein Assay Kit were obtained from Beyotime (China). Pyridine iodide (PI) ws obtained from Biosharp (Beijing, China).
Instrumentation encompassed a Nano particle tracking analyzer (ZetaView, Particle Metrix, Germany) equipped with software version 8.05.14 SP7, ultracentrifuges (Avanti JXN26 and Optima XPN100, Beckman, USA), and an FEI Tecnai transmission electron microscope (USA). Consumables such as fluorescence quantitative PCR plates, polystyrene microspheres, and filter membranes were procured from Axygen (USA), Applied-Microspheres (Germany), and Merck Millipore (Ireland), respectively. The mass spectrometer (Orbitrap™ Astral™) was purchased from Thermo Fisher Scientific (USA). The ultrasonic cell disruptor (JY92-11N) was sourced from SCIENTZ (China). The benchtop high-speed refrigerated centrifuge (TGL-20M) was obtained from CENCE (China). The vacuum freeze dryer (CV 600) was acquired from Beijing JM Technology Co., Ltd (China). The electrophoresis power supply (DYY-6C) and electrophoresis tank (DYCZ-24DN) were purchased from LIUYI (China). The microplate reader (Cmax Plus) was obtained from Molecular Devices (USA).
2.2 Bacterial culture protocol
C. butyricum colonies were initially streaked onto solid Reinforced Clostridium Medium (RCM) and incubated anaerobically at 37 °C overnight. A single colony was transferred to liquid RCM medium for 24–36 hr under anaerobic conditions. Subsequently, 1%–2% (v/v) of this culture was inoculated into 2 L MCP liquid medium and incubated anaerobically at 37 °C for 24 hr. The bacterial suspension was adjusted to OD600 = 0.5 using fresh RCM medium and further incubated anaerobically at 37 °C for 24 hr prior to extracellular vesicle (EV) extraction.
2.3 EVs isolation and characterization
EVs were purified from C. butyricum cultures following Morishita et al. (5). Briefly, 100 mL supernatant was sequentially centrifuged (4,000g, 10 min; 10,000g, 30 min; 4 °C) to remove cellular debris, filtered through a 0.45 μm membrane, and ultracentrifuged (100,000g, 70 min; 4 °C). Pelleted EVs were washed twice with PBS via ultracentrifugation (100,000g, 70 min; 4 °C) and quantified via Bradford assay.
Transmission electron microscopy (TEM) was utilized to corroborate the morphological integrity of the isolated EVs. A small volume (15 μL) of the EVs suspension was adsorbed onto a carbon-coated copper grid for one minute. Any residual liquid was subsequently removed using filter paper. Negative staining was carried out by applying 15 μL of a 2% (v/v) uranyl acetate solution to the grid for one minute, after which the stain was blotted off. The grids were air-dried completely prior to imaging. TEM observations were conducted to assess the ultrastructural features of the EVs.
Nanoparticle tracking analysis (NTA) was employed to determine the size distribution and concentration of C. butyricum EVs. Following purification, EVs samples were diluted to a concentration of approximately 5 × 107 particles/mL. For each measurement, 500 μL of the diluted sample was introduced into the nanoparticle tracking analyzer. All samples were analyzed by capturing data from five randomly selected visual fields. The instrument settings were configured as follows: camera sensitivity was set to 16, the recording duration per video was 60 seconds, and the detection threshold was fixed at 7. Remaining parameters were maintained at their default values. Data analysis was performed using the instrument’s native software (version 8.04.02). Zeta potential was also determined via nanoparticle tracking analysis (ZetaView) with five fields analyzed per sample.
All EVs preparations were subjected to sterility testing by incubation in bacterial culture medium at 37 °C for 48 hours prior to cellular treatments. Only batches confirmed to be free of microbial contamination were used in subsequent experiments.
2.4 Macrophage culture, treatment, and assessment of pyroptosis
RAW264.7 cells were seeded in 12-well plates at a density of 1×106 cells per well and maintained in F12K medium supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin-streptomycin, and 5% CO2 at 37 °C. Five experimental groups were established (Table 1). After the respective treatments, cellular morphological changes in all groups were directly observed and recorded using an inverted phase-contrast microscope, with particular attention to pyroptosis-related features such as cell swelling and membrane blebbing.
To assess plasma membrane integrity, Propidium Iodide (PI) staining was performed. The cell culture supernatant from each treatment group was collected. Adherent cells were gently washed once with phosphate-buffered saline (PBS). The supernatant and the wash PBS were combined and centrifuged to pellet any detached cells. The adherent cells were then trypsinized and pooled with the pelleted, detached cells. The combined cell pellet was washed twice with cold PBS and resuspended. A PI staining solution was added to the cell suspension at a final concentration of 0.5 mg/mL, followed by incubation at 4 °C in the dark for 15 minutes. After staining, the cells were immediately observed and imaged using a fluorescence microscope with excitation and emission wavelengths of approximately 488 nm and 630 nm, respectively. Cells exhibiting red fluorescence were considered PI-positive, indicating a loss of plasma membrane integrity.
2.5 Total RNA extraction, cDNA construction and transcriptome analysis
RNA extraction followed the manufacturer’s protocol using a Total RNA Kit I. Cells were lysed with Trizol, and RNA was purified through chloroform extraction and isopropanol precipitation. Purified RNA was dissolved in DEPC-treated water. cDNA synthesis was conducted following a denaturation step at 65 °C for 5 minutes, reverse transcription at 42 °C for 30 minutes, and enzyme inactivation at 70 °C for 15 minutes.
Transcriptome libraries were sequenced on an Illumina NovaSeq 6000 platform at Shanghai Biozeron Co., Ltd (Shanghai, China). Use Trimmatic software to filter out adapters and low-quality reads (reads with a length less than 35 NT, containing n bases, and a Q value<20) to obtain high-quality clean data. Then, the transcriptome was assembled by de novo splicing using Trinity software, and redundant sequences were removed using rseqc (RNA SEQ data quality control) software to obtain the UniGene of the transcriptome. UniGene with a length exceeding 200 bp will be compared and annotated in multiple public databases. For example, nucleic acid sequence database (NT), conserved domain database (CDD), non-redundant protein database (NR), protein family alignment and hidden Markov model database (PFAM), and homologous protein cluster database (COG).
2.6 q-PCR
Gene expression analysis of target genes was performed using q-PCR on a CFX96 real-time system (Bio-Rad, USA). Primer sequences and reaction conditions are detailed in Tables 2–4.
2.7 Bioinformatics analysis
The bioinformatics analysis of the transcriptomic data was performed using a series of computational tools and pipelines (17). To explore the role of pyroptosis-related genes (PRGs), we integrated our DEGs with a list of 121 PRGs previously reported in the literature (17–19). The gene difference analysis between groups were analyzed using edgeR software, with significant conditions of abs (logFC)>=1 and padj<0.05. According to the intersection genes (the intersection genes of differential genes and pyroptosis genes), the calculation of differential genes in sample grouping was conducted using limma software to explore the expression differences between high and low expression genes. A Pvalue less than 0.05 and | logFC |>0.5 were considered as differences. Use the clusterProfiler package to enrich upregulated differentially expressed genes and downregulated differentially expressed genes separately, using algorithms ORA and GSEA, with the database used for ORA being Gene Ontology (GO), Encyclopedia of Genes and Genomes (KEGG), Reactome. The functional databases used by GSEA are HALLMRARK, GOBP, and KEGG, and the signatures of the functions come from the MSIGDB database. Visualize the enrichment results using the enrichplot package.
2.8 Data analysis
All experiments were repeated three times and all data were presented as mean ± standard deviation (n=3). Statistical significance was assessed using independent t-tests, with p<0.05 considered significant.
3 Results
3.1 Characteristics of C. butyricum-derived EVs
After centrifugation, a pale yellow precipitate was found at the bottom of the centrifuge tube, which was completely dissolved after adding sterile water. Transmission Electron Microscope images showed that C. butyricum-derived EVs exhibited a distinct double-layered vesicle structure (Figure 1A, consistent with previous research reports (4). Nanoparticle tracking analysis (NTA) indicated that the median particle size of C. butyricum EVs was 144.4 ± 85.7 nm (Figure 1B), with a particle concentration of 3.9 × 1011 particles/mL. Zeta Potential measurements showed that the zeta potential of C. butyricum EVs was -37.56 ± 1.08 mV (Figure 1C). The protein concentration of the EVs preparation, as determined by the BCA assay, was 119 μg/mL. Consequently, the particle-to-protein ratio was calculated to be approximately 3.28 × 109 particles/μg.
Figure 1. Characteristics of C. butyricum EVs. (A) TEM image of C. butyricum-derived EVs, (B) NTA result, (C) Zeta Potential result.
3.2 Observation of cell morphology and evaluation of cell death by PI staining
To evaluate the cell death status induced by LPS, this study conducted inverted phase contrast microscopy morphological observation and propidium iodide (PI) fluorescence staining detection on five experimental groups.
The PI staining results (Figure 2) showed that compared with the negative control group and the simple EVs treatment group, the LPS treatment group exhibited high-intensity PI fluorescence signals, indicating that LPS stimulation led to significant loss of membrane integrity in RAW264.7 cells. However, no significant differences were observed in PI fluorescence intensity between the LPS treatment group, the EVs (10 μg/mL) + LPS treatment group, and EVs (20 μg/mL) + LPS treatment group.
Figure 2. RAW264.7 macrophage PI cell staining. (A) Control, (B) EVs only, (C) LPS only, (D) EVs (10 μg/mL, 12 h) + LPS (1 μg/mL, 12 h), (E) EVs (20 μg/mL, 12 h) + LPS (1 μg/mL, 12 h).
In terms of cell morphology (Supplementary Figure S1), under an inverted phase contrast microscope, it was observed that the cells in the LPS treated group exhibited significant morphological features similar to pyroptosis, such as swelling and membrane rupture. Similar phenomena were also observed in the EVs+LPS treatment group. It is worth noting that in the pure EVs treatment group, although some changes in cell morphology can be observed, the PI fluorescence signal intensity is much lower than that of the LPS treatment group and its combined treatment group.
These results indicate that LPS treatment can induce significant cell death in RAW264.7 cells, but the specific mode of death and the role of C. butyricum EVs in this process need further validation.
3.3 Inflammatory cytokines and pyroptosis-related genes
We detected the expression of four inflammatory cytokine genes, IL-1β, IL-12, IL-18, and IL-33 (Figure 3), among which IL-1β, IL-18, and IL-33 are pyroptosis-related genes (20, 21). The results showed that C. butyricum EVs modulated the LPS-induced expression of these genes, including a non-significant reduction in IL-1β (q>0.05) and a significant inhibition of IL-18 (q<0.05). Under the treatment of 10 μg/mL EVs, the expression of IL-33 induced by LPS significantly increased (q<0.001), while under the treatment of 20 μg/mL EVs, it was significantly decreased (q<0.05). Treatment with C. butyricum EVs also reduced the expression of IL-12 induced by LPS, but the results were not significant.
Figure 3. Inflammatory cytokines and pyroptosis markers mRNA expression in RAW264.7 macrophages. (A) IL-18, (B) IL-12, (C) IL-1β, (D) IL-33. *P < 0.05, **P < 0.01, ***P < 0.0001.
3.4 Identification of differentially expressed genes and key pathways in RAW264.7 cells treated with EVs
We compared the transcriptome of cells treated with 10 μg/mL C. butyricum EVs with that of the control group. According to Figures 4A, B, a total of 5871 DEGs were identified (P<0.05, | logFC |>1), with 3410 genes upregulated and 2461 genes downregulated. Figures 4C, D illustrate the KEGG pathways and GO annotations enriched by upregulated genes, while Figures 4E, F depict those for downregulated genes.
Figure 4. Identification of DEGs and Key Pathways in RAW264.7 Cells Treated with EVs. (A) The volcano plot of DEGs, (B) The heatmap of DEGs, (C) The KEGG pathways enriched by upregulated genes, (D) The GO annotations enriched by upregulated genes, (E) The KEGG pathways enriched by downregulated genes, (F) The GO annotations enriched by downregulated genes, (G) The GSEA enrichment analysis results for upregulated genes, (H) The GSEA enrichment analysis results for downregulated genes.
The enrichment analysis revealed that the upregulated gene pathways primarily involve “Cytokine-cytokine receptor interaction,” “ECM-receptor interaction,” “JAK-STAT signaling pathway,” “Hematopoietic cell lineage,” and “Viral protein interaction with cytokine and cytokine receptor.” In contrast, the major downregulated pathways include “Alcoholism,” “DNA replication,” “Biosynthesis of cofactors,” “Cell cycle,” “Fanconi anemia pathway,” “Base excision repair,” “Aminoacyl-tRNA biosynthesis,” and “Purine metabolism.”
In the biological processes (BP) category of GO annotation results, upregulated genes were predominantly enriched in “leukocyte migration,” “cell chemotaxis,” “positive regulation of cytokine production,” “cell-substrate adhesion,” “regulation of cell-cell adhesion,” “negative regulation of response to external stimulus,” and “ERK1 and ERK2 cascade.” Downregulated genes showed significant enrichment in “ncRNA metabolic process,” “ncRNA processing,” “ribonucleoprotein complex biogenesis,” “ribosome biogenesis,” “rRNA processing,” “rRNA metabolic process,” and “tRNA metabolic process.”
Regarding cellular components (CC) in GO annotations, the top five significantly enriched terms for upregulated genes included “receptor complex,” “collagen-containing extracellular matrix,” “apical part of cell,” “apical plasma membrane,” and “cluster of actin-based cell projections.” For downregulated genes, the top five significantly enriched terms were “chromosomal region,” “chromosome,” “centromeric region,” “condensed chromosome,” and “kinetochore” or “spindle.”
In the molecular function (MF) category of GO annotations, the top five significantly enriched terms for upregulated genes comprised “receptor ligand activity,” “cytokine activity,” “cytokine receptor binding,” “peptide receptor activity,” and “cell adhesion molecule binding.” Meanwhile, “catalytic activity acting on RNA,” “catalytic activity acting on a tRNA,” “catalytic activity acting on DNA,” “transferase activity transferring one-carbon groups,” and “ligase activity” were the most significantly enriched terms for downregulated genes.
Figures 4G, H represent the GSEA enrichment analysis results for upregulated and downregulated genes, respectively. According to the GSEA analysis, the functions of differentially expressed genes mainly encompass “inositol phosphate metabolism” (Enrichment Score: -0.6442, P Value: 0.0000), “regulation of vacuole organization” (Enrichment Score: -0.6694, P Value: 0.0001), “autophagosome organization” (Enrichment Score: -0.5536, P Value: 0.0002), “protein localization to extracellular region” (Enrichment Score: 0.5170, P Value: 0.0002), and “regulation of granulocyte chemotaxis” (Enrichment Score: 0.6939, P Value: 0.0003).
The results indicated that after treatment with C. butyricum EVs, the upregulated genes in RAW264.7 cells were significantly enriched in immune-modulating pathways such as “Cytokine-cytokine receptor interaction” (e.g., IL-6, TNF signaling), “JAK-STAT signaling pathway,” and “ECM-receptor interaction.” The MF of these genes were concentrated in “receptor ligand activity” and “cytokine activity,” while BP involved “leukocyte migration,” “cell chemotaxis,” and “ERK1/2 signaling cascade.” In contrast, downregulated genes primarily participated in fundamental metabolic processes including “DNA replication,” “cell cycle,” and “purine metabolism,” and were enriched in chromosome-associated structures within CC (such as “chromosomal region,” “kinetochore”). Further GSEA analysis revealed that EVs treatment significantly inhibited pathways like “phosphatidylinositol metabolism” and “vacuole organization” (P<0.001), while activating “protein extracellular localization” and “regulation of granulocyte chemotaxis”. These findings suggest that C. butyricum EVs modulate macrophage inflammatory responses and functional homeostasis by reshaping immune-related gene networks (such as cytokine signaling) and suppressing cell proliferation-related pathways.
3.5 Identification of DEGs and key pathways in LPS-stimulated RAW264.7 cells treated with EVs
We pre-treated LPS-stimulated RAW264.7 macrophages with C. butyricum EVs at concentrations of 10 μg/mL and 20 μg/mL, respectively. By comparing the transcriptomes of cells treated with 10 μg/mL of C. butyricum EVs to those of the control group, 2648 DEGs were identified based on the screening criteria (P < 0.05 and |logFC| > 1), of which 2007 genes were upregulated and 641 genes were downregulated (Figures 5A, B). Figures 5C–E present the ORA enrichment analysis results for the BP, CC, and MF categories of GO for upregulated genes. Figure 3F displays the KEGG pathway analysis results for upregulated genes.
Figure 5. The Impact of Pre-treatment with 10 μg/mL C. butyricum EVs on the Transcriptome of Macrophages Following LPS Stimulation. (A) The volcano plot of DEGs, (B) The heatmap of DEGs, (C) The Biological Process (BP) categories of GO for upregulated genes, (D) The Cellular Component (CC) categories of GO for upregulated genes, (E) The Molecular Function (MF) categories of GO for upregulated genes, (F) The KEGG pathways enriched by upregulated genes.
According to the GO annotation results, the upregulated DEGs were primarily enriched in BP such as “cellular process involved in reproduction in multicellular organisms,” “regulation of peptidase activity,” “germ cell development,” “regulation of vasculature development,” “muscle contraction,” and “extracellular structure organization.” In the CC category, significantly enriched terms included “collagen-containing extracellular matrix,” “receptor complex,” “secretory granule,” “apical part of cell,” “synaptic membrane,” and “cell projection membrane.” Notably, in the MF category, the upregulated genes were significantly enriched in terms like “channel activity,” “passive transmembrane transporter activity,” “receptor ligand activity,” “metal ion transmembrane transporter activity,” and “ion channel activity.”
The KEGG pathway analysis highlighted several important pathways, including “Herpes simplex virus 1 infection,” “Systemic lupus erythematosus,” “Alcoholism,” “Axon guidance,” “ECM-receptor interaction,” and “Amoebiasis.”
The results of pre-treatment with 20 μg/mL of C. butyricum EVs were similar to those obtained with 10 μg/mL C. butyricum EVs. By comparing the transcriptomes of cells treated with 20 μg/mL of C. butyricum EVs to those of the control group, a total of 2950 DEGs were identified, among which 2280 genes were upregulated and 670 genes were downregulated. The volcano plot and heatmap of these DEGs are shown in Figures 6A, B, respectively. Figures 6C–E present the ORA enrichment analysis results for the BP, CC, and MF categories of GO for upregulated genes. Figure 4F displays the KEGG pathway analysis results for upregulated genes.
Figure 6. The Impact of Pre-treatment with 20 μg/mL C. butyricum EVs on the Transcriptome of Macrophages Following LPS Stimulation. (A) The volcano plot of DEGs, (B) The heatmap of DEGs, (C) The Biological Process (BP) categories of GO for upregulated genes, (D) The Cellular Component (CC) categories of GO for upregulated genes, (E) The Molecular Function (MF) categories of GO for upregulated genes, (F) The KEGG pathways enriched by upregulated genes.
3.6 Expression patterns of differentially expressed PRGs in RAW264.7 cells
The qPCR results showed that C. butyricum EVs could regulate the pyroptosis-related genes IL-1β, IL-33 and IL-18. To further understand the potential mechanism by which C. butyricum EVs modulate the pyroptosis-related gene network, we analyzed the expression patterns of differentially expressed PRGs in RAW264.7 cells. We conducted an intersection analysis between a panel of PRGs and the DEGs identified from transcriptome sequencing after EVs pre-treatment (Figure 7A). A total of 14 differentially expressed PRGs were identified. The PRGs following treatment with 10 μg/mL C. butyricum EVs included Gsdma, Gsdmc2, Gsdmc4, Gzmn, IL-6, Pjvk, ZBP1, P2rx7, Nlrp1a,Tgfb2, and Timp2. The PRGs following treatment with 20 μg/mL C. butyricum EVs included Gsdma, Gsdmc4, Gzmn, IL-6, Nlrp6, Pjvk, Tnf, Zbp1, P2rx7, Pelp1, and Tgfb2. Figure 7B illustrates the expression differences of these 14 PRGs between the EVs groups and the LPS control group. Among them, Gsdma, Gsdmc4, IL-6, Nlrp6, P2rx7, Pelp1, and Timp2 were found to be upregulated in the EVs-treated groups, while Gsdmc2, Gzmn, Nlrp1a, Pjvk, Tnf, Zbp1, and Tgfb2 were downregulated. Moreover, the expression changes of these PRGs were concentration-responsive, with the extent of upregulation or downregulation being significantly altered at the higher EV concentration (20 µg/mL) compared to the lower dose (10 µg/mL). The heatmap in Figure 7C represents the correlation of expression levels among the intersecting genes. Correlation analysis indicated that among the PRGs, ZBP1 showed the strongest negative correlation with P2rx7, whereas P2rx7 exhibited the strongest positive correlation with IL-6, IL-6 with Tgfb2, and IL-6 with Timp2.
Figure 7. Expression Patterns of Differentially Expressed PRGs in RAW264.7 Cells. (A) The Venn diagram of DEGs and Pyroptosis, (B) Differential expression of PRGs, (C) The heatmap of PRGs gene expression correlation, (D) GO enrichment analysis of PRGs following treatment with 10 μg/mL C. butyricum EVs, (E) GO enrichment analysis of PRGs following treatment with 20 μg/mL C. butyricum EVs, (F) KEGG pathway enrichment analysis of PRGs following treatment with 10 μg/mL C. butyricum EVs, (G) KEGG pathway enrichment analysis of PRGs following treatment with 20 μg/mL C. butyricum EVs, (H) Enrichment analysis results of PRGs in the Reactome database following treatment with 10 μg/mL C. butyricum EVs, (I) Enrichment analysis results of PRGs in the Reactome database following treatment with 20 μg/mL C. butyricum EVs.
GO enrichment analysis suggested that these PRGs are primarily involved in “pyroptosis,” “defense response to bacterium,” “positive regulation of interleukin-1 beta production,” and “positive regulation of interleukin-1 production” (Figures 7D, E). KEGG pathway enrichment analysis indicated that these PRGs are mainly associated with the “NOD-like receptor signaling pathway,” “Malaria,” “Inflammatory bowel disease,” and “Intestinal immune network for IgA production” (Figures 7F, G). Additionally, Figures 5H, I show the enrichment analysis results for the intersecting genes in the Reactome database, indicating that PRGs are most highly enriched in “Inflammasomes,” “Nucleotide-binding domain,” and “leucine-rich repeat-containing receptor (NLR) signaling pathways.” Collectively, these bioinformatic analyses suggest that the regulatory effects of C. butyricum EVs are functionally linked to core pyroptosis and inflammatory pathways. These findings are consistent with a concentration-dependent remodeling of the pyroptosis-related transcriptional network by C. butyricum EVs, wherein the shift from 10 to 20 µg/mL differentially modulated the expression of key genes, including IL-6, P2rx7, and ZBP1.
3.7 Biological process analysis of PRGs
To gain deeper insights into the pathways associated with the identified PRGs, GSEA analysis was conducted on the 14 PRGs. Among them, 7 genes (Gsdmc4, Gzmn, IL-6, P2rx7, Tgfb2, Timp2, and ZBP1) were analyzed in the 10 μ g/mL C. butyricum EVs treatment group. Ten genes (Gsdma, Gsdmc4, Gzmn, IL-6, Nlrp6, Pjvk, Tnf, Zbp1, P2rx7, and Tgfb2) were analyzed from the 20 μg/mL C. butyricum EVs treatment group, while Pelp1, Gsdmc2 and Nlrp1a did not obtain GSEA results. Figure 8 shows the top five GSEA results of the KEGG pathway based on these genes, with Gsdmc4, Gzmn, IL-6, P2rx7, Tgfb2, Timp2, and ZBP1 from the 10 μ g/mL C. butyricum EVs treatment group, and Gsdma, Nlrp6, Pjvk, and Tnf from the 20 μ g/mL C. butyricum EVs treatment group. According to the GSEA analysis, the primary functions related to these PRGs include the “GNRH signaling pathway” (Enrichment Score: 0.5728, P Value: 0.0000), “Prostate cancer” (Enrichment Score: 0.5268, P Value: 0.0002), “Cytokine-cytokine receptor interaction” (Enrichment Score: 0.3200, P Value: 0.0002), “Pathogenic Escherichia coli infection” (Enrichment Score: 0.6012, P Value: 0.0003), and “Type II diabetes mellitus” (Enrichment Score: 0.6124, P Value: 0.0003).
Figure 8. Related biological process of differentially expressed PRGs. (A) Gsdmc4, (B) Gzmn, (C) IL-6, (D) P2rx7, (E) Tgfb2, (F) Timp2, (G) ZBP1, (H) Gsdma, (I) Nlrp6, (J) Pjvk, (K) Tnf.
3.8 qPCR validation of PRGs
Quantitative PCR (qPCR) was employed to validate the expression of 14 differentially expressed PRGs (Gsdma, Gsdmc2, Gsdmc4, Gzmn, IL-6, Pjvk, ZBP1, P2rx7, Nlrp1a, Tgfb2, Timp2, Nlrp6, Pelp1, Tnf) in LPS-induced RAW264.7 macrophage models following EVs treatment. Figure 9 demonstrated that EVs significantly inhibited the expression of ZBP1 (q < 0.001), Gzmn (q < 0.05) and Tnf (q<0.01). Moreover, EVs significantly promoted the expression of IL-6 (q < 0.01), P2rx7 (q < 0.001), Tgfb2 (q < 0.01) and Pelp1(q<0.01). Additionally, EVs appeared to inhibit the expression of Gsdma, Timp2, and Pjvk, as well as promote the expression of Gsdmc4, Gsdmc2 and Nlrp6, however these effects were not statistically significant (q > 0.05).
Figure 9. PRGs mRNA expression in RAW264.7 macrophages. (A) ZBP1, (B) IL-6, (C) P2rx7, (D) Tgfb2, (E) Pelp1, (F) Tnf, (G) Gzmn, (H) Gsdmc2, (I) Timp2, (J) Nlrp6, (K) Nlrp1a, (L) Gsdmc4, (M) Gsdma, (N) Pjvk. *P < 0.05, **P < 0.01, ***P < 0.0001.
4 Discussion
This study elucidates the modulatory effects of C. butyricum EVs on the transcriptional network of PRGs and associated signaling pathways in LPS-stimulated RAW264.7 macrophages. Through transcriptome analysis and qPCR validation, we characterized the physicochemical properties of the EVs, delineated their impact on the macrophage transcriptome, and provided a theoretical foundation for understanding the gene-regulatory potential of probiotic-derived EVs in the context of inflammatory responses.
Research has shown that extracellular vesicles are an important pathway for the interaction between gut microbiota and host cells. Bacterial extracellular vesicles can transmit various information molecules to host cells, triggering immune responses (22). This study used nanoparticle tracking analysis (NTA) to show that the median particle size of EVs was 144.4 ± 85.7 nm, with a ZETA potential of -37.56 ± 1.08 mV, which was similar to that of previous studies (5). The particle size distribution was smaller than that of EVs derived from Lactobacillus casei (202.9 ± 88.9 nm) (23). The particle-to-protein ratio of our EV preparation was calculated to be 3.28 × 109 particles/μg. A ratio below 1.5 × 109 particles/μg is typically considered indicative of impurity (Whittaker et al., 2020). Therefore, our value suggests a relatively high purity of the isolated EVs, though it also implies the potential co-isolation of a minor amount of non-vesicular protein.
Studies have shown that C. butyricum can regulate the host’s cellular immune function through EVs, regulate intestinal homeostasis and improve colitis (4, 24). We found that C. butyricum EVs modulated the expression of inflammatory cytokines in RAW264.7 macrophages, exhibiting distinct effects under basal conditions and upon LPS stimulation. Transcriptome analysis showed that EVs treatment alone significantly upregulated the JAK-STAT signaling pathway in RAW264.7 cells. Research has shown a direct correlation between the JAK-STAT pathway and inflammatory response (25, 26). The JAK-STAT pathway is the main downstream signaling pathway for key inflammatory factors such as IL-6. After IL-6 binds to receptors, it activates STAT3 through JAK1/2, which can promote the expression of inflammatory cytokines (27). QPCR validation confirmed that treatment of RAW264.7 cells with C. butyricum EVs upregulated the expression level of IL-6, which may contribute to the observed upregulation of the JAK-STAT signaling pathway. Under LPS stimulation, EVs pretreatment further upregulated the expression of IL-6 (q<0.01), while inhibiting the expression of Tnf (q<0.01). This pattern of gene regulation is consistent with previous studies indicating that IL-6 can suppress LPS-induced Tnf expression, suggesting a potential anti-inflammatory mechanism (28).
Pyroptosis is an inflammatory programmed cell death (29). During the process of cell pyroptosis, it is generally accompanied by the formation of cell membrane pores, membrane rupture, cell swelling, and release of cellular contents. The occurrence of cell pyroptosis is closely related to inflammatory cytokines, and has become a hot topic in programmed cell death (30). Cellular pyroptosis often interacts with apoptosis and necroptosis, leading to inflammatory cell death. Together, the three are synthesized into PANoptosis (31). During the process of cell pyroptosis, the cell membrane ruptures and cell contents are released through GSDM mediation, thus cell pyroptosis is also defined as GSDM mediated programmed cell death (32). The GSDM family (GSDMA, GSDMB, GSDMC, GSDMD, GSDME, PVJK/DFNB59, etc.) is the executor of cell pyroptosis (33). Most of the proteins in the GSDM family can form oligomers, insert into the cell membrane or mitochondrial membrane, form pores, cause cell rupture and pyroptosis, and release IL-1 β, IL-18, and DAMPs, thereby triggering a strong inflammatory response (34). It is generally believed that Gsdma is related to cellular mitochondrial homeostasis (35), but recent studies have shown that Gsdma can also mediate cell apoptosis. Studies have shown that severe starvation can induce cell pyroptosis through phosphorylation induced Gsdma (36). African swine fever virus (ASFV) can activate caspase-3 and caspase-4, cleave GSDMA to regulate cell pyroptosis (37), and SpeB from group A Streptococcus can cleave and activate Gsdma to mediate cell pyroptosis (38). Gsdmc can mediate cell pyroptosis, and the caspase-8/GSDMC pathway is an important pathway for cell pyroptosis. GSDMC is specifically cleaved by caspase-8 under TNF - α treatment, producing a GSDMC-N-terminal domain that forms pores on the cell membrane and induces pyroptosis (39, 40). Gsdmc2 and Gsdmc4 are mainly expressed in the stomach, large intestine, small intestine, bladder, and prostate of mice (32). Research has found that overexpression of Gsdmc2 in HEK293 cells can trigger cell pyroptosis and lytic cell death (41). Pjvk is also a member of the GSDM protein family, but it does not have the ability to punch holes on the cell membrane, and its mechanism of action in cell pyroptosis has not been reported (42). Our study found that after treating LPS induced RAW264.7 macrophages with C. butyricum EVs, although the expression of Gsdmc2 and Gsdmc4 was upregulated and the expression of Gsdma and Pjvk was downregulated. However, the expression levels of these genes were low and the observed changes were not statistically significant (q>0.05). Therefore, the regulatory contribution of these particular GSDM family members to the overall transcriptional network alterations induced by C. butyricum EVs appears to be limited.
P2rx7(P2X7R) is believed to have the function of an inflammatory regulator, which can affect the regulation of the body’s immune system, influence the occurrence of inflammation and cancer development, and is one of the key genes in the body’s immune regulation (43). P2X7R is an adenosine triphosphate (ATP) gated ion channel that can be activated by ATP and BzATP, among others. The P2X7R-NEK7-NLRP3 axis is an important pathway that triggers cell pyroptosis. Research has shown that P2X7R mediated potassium efflux can induce NLRP3 inflammasome assembly and activation, thereby inducing pyroptosis of prostate epithelial cells (44). During LPS induced apoptosis of mouse liver cells, LPS activates NLRP3 through P2X7R to activate macrophages and induce hepatocyte pyroptosis (45). Meanwhile, studies have shown that the ATP mediated P2X7R paracrine mechanism is a P2X7 receptor dependent membrane repair mechanism. P2X7R achieves isolation and detachment of damaged plasma membrane segments by inducing Ca2+influx, foaming the cell membrane, and protecting the cell body from further damage by limiting elevated intracellular Ca2+. Studies have shown that the expression of P2X7R can enhance the resistance of HEK cells to bacterial pore forming toxin attacks (46). Transcriptome analysis and qPCR results showed that C. butyricum EVs significantly increased the expression of P2rx7 in LPS induced RAW264.7 macrophages (q<0.001). This upregulation suggests that the effect of EVs is not mediated by suppressing P2rx7. Given the dual role of P2X7R in promoting pyroptosis and membrane repair, the functional consequence of its increased expression warrants further investigation to determine if it contributes to a pro-pyroptotic signal or a protective membrane repair response in this context.
Pelp1 is a scaffold protein that can act as a co regulatory factor for multiple transcription factors and nuclear receptors (47). Recent studies have found that Pelp1 also plays an important role in cell pyroptosis, and inhibiting the expression of Pelp1 can reduce cell pyroptosis. For example, miR-195-5p alleviates pyroptosis in GC-1 cells by inhibiting Pelp1 expression during OGD/R, and this protective effect is blocked when miR-195-5p is downregulated (48). Metformin induces GSDMD mediated pyroptosis in esophageal squamous cell carcinoma cells by targeting the miR-497/Pelp1 axis (49). We found through qPCR that C. butyricum EVs significantly increased the expression of Pelp1 in LPS induced RAW264.7 macrophages (q<0.001). This upregulation suggests that the modulation of the pyroptosis-related gene network by C. butyricum EVs is not mediated through the suppression of Pelp1 expression.
ZBP1 is an important protein that regulates cell death and plays a regulatory role in cell pyroptosis and necroptosis (50, 51). ZBP1 is an important regulatory factor for cell pyroptosis, and studies have found that ZBP1 can induce pyroptosis in acute liver injury through the PGAM5/ROS pathway cluster (52). ZBP1 is an important cell apoptosis regulator in the occurrence and progression of systemic lupus erythematosus (53). ZBP1 activates NLRP3 inflammasome through the RIPK1-RIPK3-caspase-8 axis (54), and ZBP1-NLRP3 inflammasome promotes the maturation of pro-inflammatory cytokines and GSDMD by activating caspase-1. GSDMD is cleaved by caspase-1, forming a pore on the cell membrane that releases cytokines and causes cell pyroptosis (55). Gzmn is a granule enzyme that induce cell apoptosis and pyroptosis, but the specific mechanism is not yet clear. Research has found that GZMA can cleave GSDMB, causing pores in the cell membrane and leading to cell pyroptosis (56). Our study found that C. butyricum EVs significantly inhibited the expression of ZBP1 and Gzmin in LPS induced RAW264.7 macrophages (q<0.05). The downregulation of these key pyroptosis-related genes suggests that ZBP1 and Gzmn may represent important nodes in the transcriptional network through which C. butyricum EVs modulate pyroptosis-related signaling.
Research has shown that bacterial EVs can predict diseases, regulate host immune responses, and serve as drug delivery carriers (57). The production of C. butyricum EVs is relatively high, almost free of toxic substances, and has great potential for medical applications. C. butyricum EVs can regulate the structure of mouse gut microbiota, regulate intestinal homeostasis, and improve acute colitis (4). C. butyricum EVs can also alleviate LPS induced acute lung injury (58). However, due to the influence of bacterial culture conditions on the types and characteristics of bacterial EVs’ contents, the medical ability of bacterial EVs varies under different culture conditions (59). Therefore, the stability of the efficacy of bacterial EVs is an important factor limiting their application. Our study elucidates the transcriptional regulatory effects of C. butyricum EVs on pyroptosis-related genes in RAW264.7 macrophages in vitro, which has significant limitations.
While this study provides novel insights into the transcriptional regulatory effects of C. butyricum EVs, it is important to acknowledge its limitations. The primary limitation lies in the lack of direct functional evidence to confirm the execution of pyroptosis, such as measuring GSDMD cleavage or LDH release. Furthermore, the characterization of EVs purity and the confirmation of membrane-enclosed cargo were limited by the lack of specific assays, such as Western blot analysis for EVs-specific markers and protease/nuclease protection assays. Additionally, the particle-to-protein ratio observed in our EVs preparations is consistent with the potential co-isolation of non-vesicular components, a recognized limitation of the differential centrifugation method employed. An additional limitation is the narrow concentration range of EVs (10 and 20 µg/mL) tested, which, while showing bioactivity, precludes a definitive dose-response analysis and the identification of a potential cytotoxic threshold. Consequently, the optimal effective concentration range of C. butyricum EVs remains to be fully elucidated. At the same time, the discovery of mouse macrophage lineage models needs to be validated in more complex in vivo systems. Future research should focus on employing more stringent EVs purification methods, incorporating the aforementioned functional assays to validate pyroptosis, and performing a comprehensive dose-response assessment using a broader concentration gradient (e.g., up to 50 µg/mL) to determine the IC50 and assess potential cytotoxicity. Ultimately evaluating the therapeutic potential of these EVs in relevant animal models of inflammatory disease.
5 Conclusion
This study successfully isolated and identified extracellular vesicles (EVs) of C. butyricum, with a median particle size of 144.4 ± 85.7 nm and a ZETA potential of -37.56 ± 1.08 mV. Transcriptome analysis revealed that EVs treatment significantly regulated 3410 genes in RAW264.7 macrophages, mainly enriched in signaling pathways such as JAK-STAT signaling pathway and Cytokine-cytokine receptor interaction. EVs pretreatment significantly regulated 2648 genes in LPS-stimulated RAW264.7 macrophages, and 14 PRGs (Gsdma, Gsdmc2, Gsdmc4, Gzmn, IL-6, Pjvk, ZBP1, P2rx7, Nlrp1a, Tgfb2, Timp2, Nlrp6, Pelp1, Tnf) were identified among them. Among them, in the 10 μg/mL EVs treatment group, the expression of IL-6, P2rx7, Tgfb2, and Timp2 was upregulated (q<0.01), while the expression of ZBP1, Gzmn, Tnf, and Gsdmc2 was inhibited (q<0.05); The 20 μg/mL EVs treatment group further upregulated NLRP6 and inhibited Pjvk (q<0.01). Gene interaction analysis showed that IL-6 was significantly positively correlated with P2rx7(r=0.82), while ZBP1 was negatively correlated with P2rx7 (r=-0.75). These data indicate that C. butyricum EVs modulate inflammatory and cytokine-related pathways in RAW264.7 macrophages through pathways such as JAK-STAT signaling pathway and Cytokine-cytokine receptor interaction and remodel the expression network of key pyroptosis-related genes. This transcriptional reprogramming suggests a potential role for C. butyricum EVs in fine-tuning the macrophage response to LPS stimulation. This study provides novel insights into the gene-regulatory function of probiotic-derived EVs.
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 below: PRJNA1308677 (SRA).
Author contributions
ZQ: Formal analysis, Funding acquisition, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. HK: Formal analysis, Investigation, Writing – review & editing. MA: Conceptualization, Supervision, Writing – review & editing. MY: Conceptualization, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research was supported by Zhejiang Provincial Natural Science Foundation of China under Grant (No. LTGN24C200009); the Zhejiang Provincial Department of Education under Grant (No. Y202353218).
Acknowledgments
The authors would like to acknowledge the support of Faculty of Applied Sciences, Universiti Teknologi MARA and Department of Food Research, Zhejiang, Zhejiang Business College for providing the research facilities.
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.
The reviewer WT declared a shared affiliation with the author HK to the handling editor at the time of the review.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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/fimmu.2025.1686347/full#supplementary-material
References
1. Chen Q, Ma L, Wen Y, Lyu W, Yu M, Yang H, et al. The effect of clostridium butyricum-derived lipoteichoic acid on lipopolysaccharide-stimulated porcine intestinal epithelial cells. Veterinary Med Sci. (2025) 11:e70157. doi: 10.1002/vms3.70157
2. Zhang C, Hou T, Yu Q, Wang J, Ni M, Zi Y, et al. Clostridium butyricum improves the intestinal health of goats by regulating the intestinal microbial community. Front Microbiol. (2022) 13:991266. doi: 10.3389/fmicb.2022.991266
3. Ma L, Tao S, Song T, Lyu W, Li Y, Wang W, et al. Clostridium butyricum and carbohydrate active enzymes contribute to the reduced fat deposition in pigs. iMeta. (2024) 3:e160. doi: 10.1002/imt2.160
4. Ma L, Shen Q, Lyu W, Lv L, Wang W, Yu M, et al. Clostridium butyricum and its derived extracellular vesicles modulate gut homeostasis and ameliorate acute experimental colitis. Microbiol Spectr. (2022) 10:e0136822. doi: 10.1128/spectrum.01368-22
5. Morishita M, Horita M, Higuchi A, Marui M, Katsumi H, and Yamamoto A. Characterizing different probiotic-derived extracellular vesicles as a novel adjuvant for immunotherapy. Mol pharmaceutics. (2021) 18:1080–92. doi: 10.1021/acs.molpharmaceut.0c01011
6. Marrufo AM and Flores-Mireles AL. Macrophage fate: to kill or not to kill? Infection Immun. (2024) 92:e0047623. doi: 10.1128/iai.00476-23
7. Shapouri-Moghaddam A, Mohammadian S, Vazini H, Taghadosi M, Esmaeili SA, Mardani F, et al. Macrophage plasticity, polarization, and function in health and disease. J Cell Physiol. (2018) 233:6425–40. doi: 10.1002/jcp.26429
8. Karki R, Lee S, Mall R, Pandian N, Wang Y, Sharma BR, et al. ZBP1-dependent inflammatory cell death, PANoptosis, and cytokine storm disrupt IFN therapeutic efficacy during coronavirus infection. Sci Immunol. (2022) 7:eabo6294. doi: 10.1126/sciimmunol.abo6294
9. Hu XM, Zheng S, Zhang Q, Wan X, Li J, Mao R, et al. PANoptosis signaling enables broad immune response in psoriasis: From pathogenesis to new therapeutic strategies. Comput Struct Biotechnol J. (2023) 23:64–76. doi: 10.1016/j.csbj.2023.11.049
10. Chen H, Xie X, Xiao H, Liang W, Lin ZJ, Lin B, et al. A pilot study about the role of PANoptosis-based genes in atherosclerosis development. J Inflammation Res. (2023) 16:6283–99. doi: 10.2147/JIR.S442260
11. Kim W, Lee EJ, Bae IH, Myoung K, Kim ST, Park PJ, et al. Lactobacillus plantarum-derived extracellular vesicles induce anti-inflammatory M2 macrophage polarization in vitro. J extracellular vesicles. (2020) 9:1793514. doi: 10.1080/20013078.2020.1793514
12. Fonseca S, Carvalho AL, Miquel-Clopés A, Jones EJ, Juodeikis R, Stentz R, et al. Extracellular vesicles produced by the human gut commensal bacterium Bacteroides thetaiotaomicron elicit anti-inflammatory responses from innate immune cells. Front Microbiol. (2022) 13:1050271. doi: 10.3389/fmicb.2022.1050271
13. Zhao S, Xiang J, Abedin M, Wang J, Zhang Z, Zhang Z, et al. Characterization and anti-inflammatory effects of Akkermansia muciniphila-derived extracellular vesicles. Microorganisms. (2025) 13:464. doi: 10.3390/microorganisms13020464
14. Zhang H, Zhao Y, Li D, Li H, Wang Z, Zhang L, et al. Anti-inflammatory Effects of Membrane Vesicles from Eubacterium rectale via the NLRP3 Signal Pathway. Probiotics antimicrobial Proteins. (2024) 17:4841–50. doi: 10.1007/s12602-024-10432-y
15. Pandey A, Li Z, Gautam M, Ghosh A, and Man SM. Molecular mechanisms of emerging inflammasome complexes and their activation and signaling in inflammation and pyroptosis. Immunol Rev. (2025) 329:e13406. doi: 10.1111/imr.13406
16. Wang L, Xie X, Chen Q, Chen Y, Xu X, and Liang T. Puerarin reduces diabetic nephropathy-induced podocyte pyroptosis by modulating the SIRT1/NLRP3/caspase-1 pathway. Mol Cell Endocrinol. (2025) 595:112409. doi: 10.1016/j.mce.2024.112409
17. Zhu H, Wu J, Li C, Zeng Z, He T, Liu X, et al. Transcriptome analysis reveals the mechanism of pyroptosis-related genes in septic cardiomyopathy. PeerJ. (2023) 11:e16214. doi: 10.7717/peerj.16214
18. Song W, Ren J, Xiang R, Kong C, and Fu T. Identification of pyroptosis-related subtypes, the development of a prognosis model, and characterization of tumor microenvironment infiltration in colorectal cancer. Oncoimmunology. (2021) 10:1987636. doi: 10.1080/2162402X.2021.1987636
19. Zeng R, Huang S, Qiu X, Zhuo Z, Wu H, Jiang L, et al. Predicting the prognosis of esophageal adenocarcinoma by a pyroptosis-related gene signature. Front Pharmacol. (2021) 12:767187. doi: 10.3389/fphar.2021.767187
20. Zhan P, Feng Z, Huang X, Xu H, Xu S, and Wang S. Shared pyroptosis pathways and crosstalk genes underpin inflammatory links between periodontitis and atherosclerosis. Immunobiology. (2025) 230:152880. doi: 10.1016/j.imbio.2025.152880
21. Yang L, Zhang Y, Chai Z, Zhou Y, Li Z, and Wei Y. Regulation of pyroptosis by NF-κB signaling. Front Cell Death. (2025) 3:1503799. doi: 10.3389/fceld.2024.1503799
22. Kim JH, Lee J, Park J, and Gho YS. Gram-negative and Gram-positive bacterial extracellular vesicles. Semin Cell Dev Biol. (2015) 40:97–104. doi: 10.1016/j.semcdb.2015.02.006
23. Bäuerl C, Coll-Marqués JM, Tarazona-González C, and Pérez-Martínez G. Lactobacillus casei extracellular vesicles stimulate EGFR pathway likely due to the presence of proteins P40 and P75 bound to their surface. Sci Rep. (2020) 10:19237. doi: 10.1038/s41598-020-75930-9
24. Ma L, Lyu W, Song Y, Chen K, Lv L, Yang H, et al. Anti-Inflammatory Effect of Clostridium butyricum-Derived Extracellular Vesicles in Ulcerative Colitis: Impact on Host microRNAs Expressions and Gut Microbiome Profiles. Mol Nutr Food Res. (2023) 67:e2200884. doi: 10.1002/mnfr.202200884
25. Yang I, Jeong NH, Choi YA, Kwon TK, Lee S, Khang D, et al. Monotropein mitigates atopic dermatitis-like skin inflammation through JAK/STAT signaling pathway inhibition. Biomedicine pharmacotherapy = Biomedecine pharmacotherapie. (2024) 176:116911. doi: 10.1016/j.biopha.2024.116911
26. Samra S, Bergerson JRE, Freeman AF, and Turvey SE. JAK-STAT signaling pathway, immunodeficiency, inflammation, immune dysregulation, and inborn errors of immunity. J Allergy Clin Immunol. (2024) 155:357–67. doi: 10.1016/j.jaci.2024.09.020
27. Yeung YT, Aziz F, Guerrero-Castilla A, and Arguelles S. Signaling pathways in inflammation and anti-inflammatory therapies. Curr Pharm design. (2018) 24:1449–84. doi: 10.2174/1381612824666180327165604
28. Petersen A and Pedersen B. The role of IL-6 in mediating the anti inflammatory. J Physiol Pharmacol. (2006) 57:43–51.
29. Rao Z, Zhu Y, Yang P, Chen Z, Xia Y, Qiao C, et al. Pyroptosis in inflammatory diseases and cancer. Theranostics. (2022) 12:4310–29. doi: 10.7150/thno.71086
30. Vasudevan SO, Behl B, and Rathinam VA. Pyroptosis-induced inflammation and tissue damage. Semin Immunol. (2023) 69:101781. doi: 10.1016/j.smim.2023.101781
31. Chen W, Gullett JM, Tweedell RE, and Kanneganti TD. Innate immune inflammatory cell death: PANoptosis and PANoptosomes in host defense and disease. Eur J Immunol. (2023) 53:e2250235. doi: 10.1002/eji.202250235
32. Zuo Y, Chen L, Gu H, He X, Ye Z, Wang Z, et al. GSDMD-mediated pyroptosis: a critical mechanism of diabetic nephropathy. Expert Rev Mol Med. (2021) 23:e23. doi: 10.1017/erm.2021.27
33. Shi J, Zhao Y, Wang K, Shi X, Wang Y, Huang H, et al. Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death. Nature. (2015) 526:660–5. doi: 10.1038/nature15514
34. Wu C, Lu W, Zhang Y, Zhang G, Shi X, Hisada Y, et al. Inflammasome activation triggers blood clotting and host death through pyroptosis. Immunity. (2019) 50:1401–11.e4. doi: 10.1016/j.immuni.2019.04.003
35. Zou J, Zheng Y, Huang Y, Tang D, Kang R, and Chen R. The versatile gasdermin family: their function and roles in diseases. Front Immunol. (2021) 12:751533. doi: 10.3389/fimmu.2021.751533
36. Li X, Li X, Xiang C, Cao J, Guo J, Zhu S, et al. Starvation-induced phosphorylation activates gasdermin A to initiate pyroptosis. Cell Rep. (2024) 43:114728. doi: 10.1016/j.celrep.2024.114728
37. Li S, Song J, Liu J, Zhou S, Zhao G, Li T, et al. African swine fever virus infection regulates pyroptosis by cleaving gasdermin A via active caspase-3 and caspase-4. J Biol Chem. (2024) 300:107307. doi: 10.1016/j.jbc.2024.107307
38. Yin H, Zhao Z, Yan Y, Yuan Y, Qu W, Wang H, et al. Seneca Valley virus circumvents Gasdermin A-mediated inflammation by targeting the pore-formation domain for cleavage. mBio. (2024) 15:e0168024. doi: 10.1128/mbio.01680-24
39. Feng Y, Li M, Yangzhong X, Zhang X, Zu A, Hou Y, et al. Pyroptosis in inflammation-related respiratory disease. J Physiol Biochem. (2022) 78:721–37. doi: 10.1007/s13105-022-00909-1
40. Hou J, Zhao R, Xia W, Chang CW, You Y, Hsu JM, et al. PD-L1-mediated gasdermin C expression switches apoptosis to pyroptosis in cancer cells and facilitates tumour necrosis. Nat Cell Biol. (2020) 22:1264–75. doi: 10.1038/s41556-020-0575-z
41. Xi R, Montague J, Lin X, Lu C, Lei W, Tanaka K, et al. Up-regulation of gasdermin C in mouse small intestine is associated with lytic cell death in enterocytes in worm-induced type 2 immunity. Proc Natl Acad Sci United States America. (2021) 118:e2026307118. doi: 10.1073/pnas.2026307118
42. Tang D, Zheng Y, Wang G, Sheng C, Liu Z, Wang B, et al. HPV18 E6 inhibits α-ketoglutarate-induced pyroptosis of esophageal squamous cell carcinoma cells via the P53/MDH1/ROS/GSDMC pathway. FEBS Open Bio. (2023) 13:1522–35. doi: 10.1002/2211-5463.13666
43. Yu Q, Wang X, Li X, Bai X, Zhao R, and Peng X. Purinergic P2X7R as a potential target for pancreatic cancer. Clin Trans Oncol. (2023) 25:2297–305. doi: 10.1007/s12094-023-03123-7
44. Chen L, Liu Y, Yue S, Wang H, Chen J, Ma W, et al. P2X7R modulates NEK7-NLRP3 interaction to exacerbate experimental autoimmune prostatitis via GSDMD-mediated prostate epithelial cell pyroptosis. Int J Biol Sci. (2024) 20:3393–411. doi: 10.7150/ijbs.94704
45. Luo L, Fang Y, Yuan Q, Liao J, and Zhang Z. LPS activated macrophages induced hepatocyte pyroptosis via P2X7R activation of NLRP3 in mice. Iranian J immunology: IJI. (2022) 19:4. doi: 10.22034/IJI.2022.90579.2016
46. Schoenauer R, Atanassoff AP, Wolfmeier H, Pelegrin P, Babiychuk EB, and Draeger A. P2X7 receptors mediate resistance to toxin-induced cell lysis. Biochim Biophys Acta (BBA)-Molecular Cell Res. (2014) 1843:915–22. doi: 10.1016/j.bbamcr.2014.01.024
47. Sareddy GR and Vadlamudi RK. PELP1: Structure, biological function and clinical significance. Gene. (2016) 585:128–34. doi: 10.1016/j.gene.2016.03.017
48. He KX, Xu L, Ning JZ, and Cheng F. MiR-195-5p is involved in testicular ischemia/reperfusion injury by directly targeting PELP1 and regulating spermatogonia pyroptosis. Int Immunopharmacol. (2023) 121:110427. doi: 10.1016/j.intimp.2023.110427
49. Wang L, Li K, Lin X, Yao Z, Wang S, Xiong X, et al. Metformin induces human esophageal carcinoma cell pyroptosis by targeting the miR-497/PELP1 axis. Cancer Lett. (2019) 450:22–31. doi: 10.1016/j.canlet.2019.02.014
50. Banoth B, Tuladhar S, Karki R, Sharma BR, Briard B, Kesavardhana S, et al. ZBP1 promotes fungi-induced inflammasome activation and pyroptosis, apoptosis, and necroptosis (PANoptosis). J Biol Chem. (2020) 295:18276–83. doi: 10.1074/jbc.RA120.015924
51. Oh S and Lee S. Recent advances in ZBP1-derived PANoptosis against viral infections. Front Immunol. (2023) 14:1148727. doi: 10.3389/fimmu.2023.1148727
52. Yan S, Yu L, Chen Z, Xie D, Huang Z, and Ouyang S. ZBP1 promotes hepatocyte pyroptosis in acute liver injury by regulating the PGAM5/ROS pathway. Ann Hepatol. (2024) 29:101475. doi: 10.1016/j.aohep.2024.101475
53. Huang Y, Yang DD, Li XY, Fang DL, and Zhou WJ. ZBP1 is a significant pyroptosis regulator for systemic lupus erythematosus. Ann Trans Med. (2021) 9:1773. doi: 10.21037/atm-21-6193
54. Kuriakose T, Man SM, Malireddi RK, Karki R, Kesavardhana S, Place DE, et al. ZBP1/DAI is an innate sensor of influenza virus triggering the NLRP3 inflammasome and programmed cell death pathways. Sci Immunol. (2016) 1:aag2045. doi: 10.1126/sciimmunol.aag2045
55. Kayagaki N, Stowe IB, Lee BL, O’Rourke K, Anderson K, Warming S, et al. Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling. Nature. (2015) 526:666–71. doi: 10.1038/nature15541
56. Zhou Z, He H, Wang K, Shi X, Wang Y, Su Y, et al. Granzyme A from cytotoxic lymphocytes cleaves GSDMB to trigger pyroptosis in target cells. Sci (New York N.Y.). (2020) 368:eaaz7548. doi: 10.1126/science.aaz7548
57. Liu C, Yazdani N, Moran CS, Salomon C, Seneviratne CJ, Ivanovski S, et al. Unveiling clinical applications of bacterial extracellular vesicles as natural nanomaterials in disease diagnosis and therapeutics. Acta biomaterialia. (2024) 180:18–45. doi: 10.1016/j.actbio.2024.04.022
58. Liu L, Zheng C, Xu Z, Wang Z, Zhong Y, He Z, et al. Intranasal administration of Clostridium butyricum and its derived extracellular vesicles alleviate LPS-induced acute lung injury. Microbiol Spectr. (2024) 12:e0210824. doi: 10.1128/spectrum.02108-24
59. Müller L, Kuhn T, Koch M, and Fuhrmann G. Stimulation of probiotic bacteria induces release of membrane vesicles with augmented anti-inflammatory activity. ACS Appl Bio materials. (2021) 4:3739–48. doi: 10.1021/acsabm.0c01136
Keywords: Clostridium butyricum, extracellular vesicles, pyroptosis, macrophages, bioinformatics
Citation: Zhang Q, He K, Aazmi MS and Yahya MFZR (2026) Clostridium butyricum extracellular vesicles remodel the transcriptional network of pyroptosis-related genes in LPS-stimulated macrophages. Front. Immunol. 16:1686347. doi: 10.3389/fimmu.2025.1686347
Received: 15 August 2025; Accepted: 24 November 2025; Revised: 19 November 2025;
Published: 05 January 2026.
Edited by:
Yingping Xiao, Zhejiang Academy of Agricultural Sciences, ChinaReviewed by:
Jianping Quan, South China Agricultural University, ChinaZhang Yibo, Liaoning Medical University, China
Wenjie Tang, Zhejiang University, China
Copyright © 2026 Zhang, He, Aazmi and Yahya. 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: Mohd Fakharul Zaman Raja Yahya, ZmFraGFydWx6YW1hbkB1aXRtLmVkdS5teQ==
Kailun He3