Study of the inflammatory activating process in the early stage of Fusobacterium nucleatum infected PDLSCs

Human PDLSCs from healthy and young volunteers were successfully isolated and cultured as described in the Materials and Methods. The cultured PDLSCs exhibited a spindle-shaped fibroblast-like morphology (Fig. S1a). In the multidifferentiation assay, Alizarin Red–positive mineralized matrix and Oil Red O–positive lipid droplets were observed (Fig. S1b, c). Immunophenotypic analysis showed that PDLSCs expressed MSC-specific surface markers, but not hematopoietic or endothelial cell-specific markers (Fig. S1d).

F. nucleatum inhibits proliferation, and facilitates apoptosis, ferroptosis, and inflammatory cytokine production in PDLSCs

To determine the pathogenic effect of F. nucleatum on PDLSCs, we first evaluated the viability of PDLSCs exposed to different MOIs of F. nucleatum. The results showed that F. nucleatum significantly inhibited the proliferation of PDLSCs in a time- and dose-dependent manner (P < 0.001) (Fig. 1a–c). Proliferation was almost blocked at MOIs of 200 and 400. Next, we evaluated the apoptosis ratios of F. nucleatum-infected PDLSCs. As shown in Fig. 1d and Fig. S2, F. nucleatum significantly increased the apoptosis rate of PDLSCs in a dose- and time-dependent manner (P < 0.05). Notably, early apoptosis mainly occurred at 6 h, while late apoptotic cells accounted for a substantial portion of total apoptotic cells at 24 h and 48 h.

Fig. 1figure 1

F. nucleatum inhibits proliferation, facilitates apoptosis, and ferroptosis in PDLSCs. a Cell-counting assay of PDLSCs cocultured with or without F. nucleatum. b, c Cell proliferation rate of PDLSCs detected by EdU assay. Scale bar: 100 μm. d Cell apoptosis examined using Annexin V/PI staining. Annexin V-/PI- represents live cells, Annexin V + /PI- early apoptosis, Annexin V + /PI + late apoptosis, and Annexin V-/PI + necrosis. e Intracellular Fe2+ detected by FerroOrange. Scale bar: 20 μm. f Confocal images of JC-1 in PDLSCs. Scale bar: 20 μm. g Quantitative assessment of FerroOrange and JC-1 fluorescence. h The ratio of JC-1 monomers/aggregates. Data were expressed as mean ± SD. (n = 3) (*P < 0.05; **P < 0.01; ***P < 0.001, compared with the control group)

Ferroptosis, which is a novel necrotic cell death pathway, is triggered by iron overload.21 Perturbations in iron homeostasis are major pathogenic strategies for bacterial infection.22 To investigate whether F. nucleatum induced ferroptosis in PDLSCs, we compared intracellular free iron levels between normal and F. nucleatum infected PDLSCs. The fluorescence intensity of Fe2+ was significantly enhanced in the F. nucleatum-infected group (P < 0.01) (Fig. 1e, g). Iron overload leads to mitochondrial dysfunction, which mainly manifests as mitochondrial membrane potential (MitoMP) depolarization.23,24 Therefore, we evaluated the intracellular MitoMP of PDLSCs using JC-1 fluorescent dye and observed that F. nucleatum treatment reduced the fluorescence intensity of JC-1 aggregates and enhanced the green fluorescence of JC-1 monomers (Fig. 1f-g). Quantitative analysis showed that the JC-1 monomer/aggregate intensity ratio was increased after F. nucleatum infection, indicating that F. nucleatum-induced iron overload may impair mitochondrial function in PDLSCs (Fig. 1h). These results first showed that F. nucleatum treatment could increase the intracellular labile iron levels and promote MitoMP depolarization in host cells.

We next evaluated whether F. nucleatum could trigger inflammatory responses in PDLSCs by qRT-PCR and ELISA. As shown in Fig. 2a, the gene expression levels of IL-1β, IL-6, and IL-8 increased with increasing stimulation time and peaked at 6 h in a dose-dependent manner, and they regressed with the duration of stimulation. At the protein level, IL-1β, IL-6, and IL-8 were also consistently increased with increasing stimulation time in the early stage and reached a maximum level at 12 h (Fig. 2b). These results suggest the potential immunomodulatory effect of PDLSCs under the stimulation of periodontal pathogens.

Fig. 2figure 2

F. nucleatum promotes the production of inflammatory cytokines in PDLSCs. a The mRNA expression detected by qRT-PCR. b Protein levels measured by ELISA. Data were expressed as mean ± SD. (n = 3) (*P < 0.05; **P < 0.01; ***P < 0.001, compared with the control group)

FadA activates NF-κB and MAPK signaling pathways by interacting with PEBP1

Fusobacterium adhesin A (FadA) has been reported to be one of the most important adhesins and virulence factors of F. nucleatum.16 To explore the molecular mechanism of F. nucleatum infection, we investigated the pathogenic effect of FadA on PDLSCs. We obtained recombinant histidine (His)-tagged FadA through an E. coli expression system (Fig. S3a–b). After the addition of 0.5 mg·mL−1 FadA protein, the mRNA levels of IL1β, IL6, and IL8 were significantly increased compared with the controls (P < 0.001) (Fig. 3a). At the protein level, the IL1β level increased at 1 h, while IL6 and IL8 increased at 3 h after FadA stimulation (Fig. 3b).

Fig. 3figure 3

F. nucleatum adhesin FadA promotes inflammatory response via interaction with PEBP1. a The mRNA expression detected by qRT-PCR. b Protein levels measured by ELISA. c Co-IP analysis showed the interaction between FadA and PEBP1. d SPR analysis of PEBP1 binding to FadA. PEBP1 at various concentrations was used as analyte to detect binding to FadA immobilized on the sensor chip. e Effects of FadA on the phosphorylation of Raf1, IKK, PEBP1, p38, JNK, ERK, and NF-κB p65 detected by Western blot analysis. f Quantitative analysis of altered protein expression of these proteins used ImageJ. g Schematic diagram of the FadA-PEBP1-Raf1-MAPK pathway and FadA-PEBP1-IKK-NF-κB pathway. Created with BiorRender.com. Data were expressed as mean ± SD (n = 3) (*P < 0.05; **P < 0.01; ***P < 0.001, compared with the control group)

Next, we isolated FadA-binding proteins by a His pull-down assay, and identified all pull-down proteins by mass spectrometry. Among the candidates (Table S3), a cytoplasmic protein phosphatidylethanolamine-binding protein 1 (PEBP1), which is also known as Raf1 kinase inhibitory protein (RKIP),25,26,27 was proven to be co-immunoprecipitated with FadA by Co-IP assay (Fig. 3c). The direct binding of PEBP1 to FadA was further confirmed by surface plasmon resonance (SPR) analysis (Fig. S3c, d and Fig. 3d).

To study whether FadA induces the inflammatory response by interacting with PEBP1, we first evaluated the phosphorylation state of PEBP1. Figure 3e shows the binding of FadA-PEBP1 phosphorylated PEBP1 at S153. As the devitalization of PEBP1 could activate Raf1 and IKK, we hypothesized that FadA promoted the production of proinflammatory cytokines by activating the NF-κB and MAPK signaling pathways by binding to PEBP1. Western blot analysis showed that Raf1 and IKK were both significantly activated, and ERK-JNK-p38 MAPKs and NF-κB-p65 were subsequently significantly activated in PDLSCs (P < 0.05) (Fig. 3e–g). These findings suggest that FadA acts as a pathogenic effector of F. nucleatum and can initiate intracellular immune signal transduction in PDLSCs.

F. nucleatum infection induces dynamic gene activation in PDLSCs

At present, the gene regulation process in oral cells in the early stage of F. nucleatum infection is unknown. We performed time course RNA-seq analysis of PDLSCs under F. nucleatum infection for 1 h, 3 h, 6 h, and 12 h. Principal component analysis (PCA) showed that the transcriptomes of the control groups exhibited a stable gene expression state, while the F. nucleatum-infected groups changed continuously in a particular pattern (Fig. S4). The gene expression profiles of the control group and experimental group at 1, 3, 6, and 12 h were compared by DESeq2, and 25, 271, 495, and 619 differentially expressed genes (DEGs) were identified at each time point, respectively (Fig. 4a). Among these genes, 18, 235, 415, and 495 were upregulated, while 7, 36, 80, and 124 were downregulated at each time point (Fig. 4b). Notably, the Venn diagram showed that 4 DEGs were consistently upregulated in the F. nucleatum-stimulated group across the four time points (Fig. 4a and Fig. S5). CXCL1 and CXCL2 are two vital neutrophil chemoattractants. The continuous upregulation of these two chemokines indicated the key role of PDLSCs in recruiting immune cells during F. nucleatum infection.

Fig. 4figure 4

Dynamic gene expression patterns of PDLSCs under F. nucleatum Infection. a Venn diagram depicted the extent of overlap between DEGs at different time points in PDLSCs. b The histograms of upregulated (red) and downregulated (blue) DEGs at each time point. c Mfuzz analysis clustered genes into different classes. Green font represents proliferation-related genes, blue font apoptosis-related genes, orange font ferroptosis-related genes, and red font inflammation-related genes

To reveal the gene expression patterns related to the cytological phenotypes that are altered by F. nucleatum, we clustered all the genes into 30 expression patterns in control group and 20 clusters in F. nucleatum group by Mfuzz analysis (Fig. S6). As shown in Fig. 4c, we found several proliferation-inhibition genes that showed a continuous upregulation trend after F. nucleatum infection, including IFIT3, ING1, and EIF2AK2. Apoptosis-related gene such as IFI27 and PML showed a consistent increasing trend. Notably, some genes related to iron metabolism (such as TFRC and TF) were continuously upregulated while those associated with ROS detoxification (such as GPX4) were gradually decreasing, suggest that F. nucleatum could induce ferroptosis by aggravating intracellular iron overload and inhibiting lipid hydroperoxide detoxification. With respect to inflammation-related DEGs, various proinflammatory cytokines were classified into cluster 20 in F. nucleatum group (such as IL-1β, IL-6, and IL-8), which showed a gradual increasing trend. The expression levels of these genes gradually increased during infection. Some of the genes were validated by qRT-PCR, and the expression levels were consistent with the RNA-Seq results (Fig. S7).

Inflammatory genes are expressed sequentially in response to F. nucleatum infection

To explore the intracellular cascade induced by F. nucleatum-infection, we further analyzed the coexpressed genes at two adjacent time points. As shown in Fig. 5a, 9 members of the CXC chemokine subfamily, 11 members of the CC chemokine subfamily, and some proinflammatory cytokines formed a coexpression network. Interestingly, all of the chemokines in this network are inflammatory chemokines and are mainly involved in the recruitment of leukocytes to inflamed tissues.28 In addition, CCL11 and CCL20 have homeostatic functions and act as dual-function chemokines. The inflammatory chemokines CXCL1, CXCL2, CCL3L1, and CCL3L3 and inflammatory cytokine TNF were first released in response to F. nucleatum stimulation. With increasing duration of infection, the types and expression levels of chemokines increased gradually, which suggesting the potential of PDLSCs to recruit immune cells.

Fig. 5figure 5

Gene regulation process in response to F. nucleatum infection in PDLSCs. a Coexpression network of DEGs showed the intracellular cascade induced by F. nucleatum-infection. (The figure only shows the results of inflammation-related genes.) b Transcription factors NFKB1 and NFKB2 up-regulated multiple DEGs

Additionally, several transcription factors were differentially expressed at 3 h and regulated a series of target genes among DEGs at 3 h, 6 h, and 12 h. Notably, NFKB1 and NFKB2, two central activators of genes involved in inflammation and immune function, were significantly upregulated at 3 h after F. nucleatum-infection and showed sustained activation at 6 h and 12 h. The increases in NFKB1 and NFKB2 sequentially regulated downstream target genes of the NF-κB signaling pathway, such as MAP3K8, NFKBIA, and REL (Fig. 5b).

To improve our understanding of the biological functions of DEGs, we further performed Metascape analysis and displayed the top 20 enriched clusters in Fig. 6a. The regulation of cytokine production, the MAPK cascade, the apoptotic signaling pathway and the negative regulation of cell proliferation were highlighted in the network. Consistent with these results, GO analysis of the DEGs indicated that they were involved in the inflammatory response process (GO:0006954) and chemokine-mediated signaling pathway (GO:0070098) during F. nucleatum infection. Immune-related processes were changed in the initial phase of infection. With increasing duration of infection, genes were enriched in the apoptotic process (GO:0006915) and the negative regulation of cell proliferation (GO:0008285) (Fig. 6b and Fig. S8).

Fig. 6figure 6

GO and KEGG pathway enrichment analysis of differentially expressed genes. a Network of enriched terms colored by cluster identity. b GO biological process analysis of DEGs. c KEGG pathway analysis of DEGs of 1 h, 3 h, 6 h, and 12 h

The top 10 abundant KEGG pathways are displayed in Fig. 6c, which shows that F. nucleatum stimulation mainly resulted in changes in immune-related pathways. DEGs at 1 h were mainly enriched in the cytokine receptor pathway, while at the three other time points, the IL-17, TNF and NF-κB signaling pathways were prominent. Pathview analysis showed the activated genes in the TNF, IL-17, and NF-κB signaling pathways. Most of the significantly altered genes were associated with survival and inflammation (Fig. S9). These findings suggested that the reactions of PDLSCs to F. nucleatum infection involved the recognition of bacterial surface epitopes by host receptors at 1 h, followed by activation of the host defense system within 12 h after infection. These results collectively indicated that F. nucleatum could induce an inflammatory response in PDLSCs associated with activation of the NF-κB signaling pathway and the production of inflammatory chemokines.

Screening miRNAs and transcription factors and constructing the gene regulatory network of PDLSCs

Next, we sought to identify the miRNAs and their target mRNAs that were specifically expressed during F. nucleatum infection as described in the Materials and Methods. The differentially expressed miRNAs at each time point are listed in Table S4. After matching the differentially expressed miRNAs with their predicted target genes, we constructed networks containing 3 upregulated miRNAs and 1 downregulated miRNA associated with a total of 22 target genes at 6 h, and 3 downregulated miRNAs associated with 23 target genes at 12 h (Fig. 7a). Notably, target genes of miR-4257 were significantly enriched in cysteine-type endopeptidase activity involved in apoptotic process; target genes of miR-4696 were significantly enriched in iron ion homeostasis and positive regulation of MAPK cascade (Fig. S10).

Fig. 7figure 7

Construction of miRNA network and transcription factors network. a miRNA-mRNA regulatory network. b K-means method clustered the DEGs into 5 modules. c Gene expression trends in 5 modules. d The functional linkages between the modules. Red arrows represent positive regulation, and blue arrows represent negative regulation. e Regulatory networks of module 1 and module 5

Next, we clustered the DEGs by the similarity of expression patterns to investigate gene regulation during F. nucleatum infection. As shown in Fig. 7b, all DEGs were grouped into 5 clusters, which were named Module 1- Module 5 (M1-M5): (1) genes in Module 1 were continuously downregulated during F. nucleatum infection; (2) genes in Module 2 were upregulated in the early stages of F. nucleatum infection and downregulated at the following time points; (3) genes in Module 3 was gradually up-regulated within 6 h after infection; (4) genes in Module 4 were upregulated after 1 h; and (5) and genes in Module 5 were downregulated between 1 h and 3 h and recovered after 3 h (Fig. 7c). To reveal the overall regulatory relationships of the F. nucleatum-infected PDLSCs, we constructed a regulatory network between modules using high-dimensional ordinary differential equations,29 as shown in Fig. 7d.

At the functional level, GO analysis confirmed that genes in M1 were involved in signal transduction, genes in M5 were enriched in activating G-protein coupled receptor signaling pathway, and genes in M4 were involved in the defense response to virus and the innate immune response. To decipher the inter-module regulatory relationships, we integrated the regulatory linkages between the TFs and their target genes, and constructed regulatory networks (Fig. 7e and Fig. S11). Consistent with previous results, NFKB1 was predicted in all of the modules. Taken together, these results indicate the pivotal role of NFKB1 and revealed the gene regulation process at different time points of F. nucleatum infection.

Identification of the therapeutic targets to attenuate the negative effects of F. nucleatum infection

We used cogena30 to perform coexpression analysis and divided the DEGs into three clusters (Fig. 8a). KEGG pathway enrichment analysis of the coexpressed genes showed that genes in clusters 1 and 2 were highly enriched in immune-related pathways, while genes in cluster 3 were enriched in calcium signaling pathway and inositol phosphate metabolism (Fig. 8b). Considering the major pathological changes in F. nucleatum-infected PDLSCs, we further performed drug repositioning analysis of clusters 1 and 2 to identify potential drug candidates. The list of drug candidates targeting the coexpressed genes in clusters 1 and 2 is shown in Fig. 8c. Pathway enrichment analysis of the target genes of each candidate drug was performed to narrow the field of candidates (Fig. S12). Based on the enrichment results, we ultimately selected six drugs and assessed the therapeutic value of these six candidates.

Fig. 8figure 8

Drug repositioning based on transcriptome revealed candidate drugs. a K-means method clustered the DEGs of 1 h, 3 h, 6 h and 12 h into 3 clusters. b KEGG pathway analysis for coexpressed genes generated by cogena. c Computational drug repositioning for the coexpressed genes. Drug candidates target DEGs in (1) cluster 1 and (2) cluster 2

Cytotoxicity assays helped us to select the drug concentrations that cells could tolerate for the follow-up experiments (Fig. S13a). After 12 h of F. nucleatum infection, PDLSCs had significantly elevated mRNA expression of IL1β, IL6 and IL8 (P < 0.05), and all six candidates significantly reduced the expression level of these inflammatory genes (P < 0.05) (Fig. S13b). These results validated the efficacy of our predicted agents.

Among the drugs investigated, piperlongumine and fisetin exhibited the best attenuating effects, which prompted us to further assess their effects on FadA-induced inflammation. Similarly, piperlongumine and fisetin significantly decreased the FadA-induced proinflammatory cytokine production (P < 0.001) (Fig. 9a, b). Considering the ferroptotic effects of F. nucleatum on PDLSCs, we further examined the effects of piperlongumine and fisetin on ferroptosis. As shown in Fig. 9c–f, piperlongumine and fisetin reversed this trend, reduced the level of intracellular Fe2+, and ameliorated the impairment in mitochondrial function. As the IL-17 signaling pathway and NF-κB signaling pathway were enriched by KEGG analysis, we used molecular docking to simulate the binding of piperlongumine or fisetin with key protein targets of the IL-17 and NF-κB signaling pathways. The predicted hub targets of piperlongumine are displayed in the 3D results in Fig. 9g. To validate these prediction results, we next evaluated the effects of piperlongumine on the F. nucleatum-induced expression of downstream markers. Western blot analysis showed that F. nucleatum infection increased the phosphorylation of IKK, p65 and p38, and piperlongumine treatment significantly inhibited the F. nucleatum-induced activation of IKK, p65 and p38 (Fig. 9h–i).

Fig. 9figure 9

Experimental validation of the candidate drugs. a QRT-PCR and b ELISA showed piperlongumine and fisetin suppressed the expression of inflammatory cytokines. c Piperlongumine and fisetin reversed the increased intracellular Fe2+ induced by F. nucleatum stimulation. Scale bar: 20 μm. d Effect of piperlongumine and fisetin on the mitochondrial membrane potential. Scale bar: 20 μm. e Quantitative assessment of FerroOrange and JC-1 fluorescence. f The ratio of JC-1 monomers/aggregates. g Molecular docking between piperlongumine and main targets. h Effects of piperlongumine on the phosphorylation of IKK, p38, and NF-κB p65 detected by Western blot analysis. i Quantitative analysis of altered protein expression of these proteins used ImageJ. Data were expressed as mean ± SD. (n = 3) (*P < 0.05; **P < 0.01; ***P < 0.001, compared with the control group; #P < 0.05; ##P < 0.01; ###P < 0.001)

留言 (0)

沒有登入
gif