PubMed2026
BACKGROUND: Increasing evidence suggests that Porphyromonas gingivalis (Pg) is associated with oral squamous cell carcinoma (OSCC) development and progression. This study aimed to identify Pg-associated genes with prognostic relevance in OSCC through integrated bioinformatics analysis. METHODS: OSCC-related differentially expressed genes (DEGs) were identified from the The Cancer Genome Atlas (TCGA)-OSCC cohort and intersected with Pg supernatant-associated DEGs from GSE192887. Raw count data were analyzed with DESeq2, whereas transcripts per million (TPM)-transformed expression values were used for downstream visualization and model construction. Weighted gene co-expression network analysis (WGCNA), univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariable Cox modeling were used to develop a seven-gene prognostic signature, which was externally evaluated in GSE41613. Additional analyses examined treatment-associated expression changes in the seven model genes, pairwise correlations among the model genes, and correlations between Pg supernatant-associated differentially expressed gene (PgSDEG)-derived module eigengenes and immune-cell fractions. Quantitative reverse-transcription polymerase chain reaction (qRT-PCR) was performed in eight paired OSCC and adjacent non-tumor tissues and in supplemented-brain heart infusion (BHI) vehicle-control and Pg culture-supernatant-treated HOK, HSC-3, and CAL-27 cells. RESULTS: A prognostic signature comprising CXCL8, GAST, HBQ1, PADI3, STC1, TEX19, and TMEM92 was established. The signature showed limited-to-moderate discrimination in the TCGA training cohort, with 1-, 3-, and 5-year areas under the curve (AUCs) of 0.68, 0.69, and 0.69, respectively, and limited discrimination in the GSE41613 external cohort (AUCs: 0.66, 0.67, and 0.61). Kaplan-Meier analysis showed poorer survival in the high-risk group in both cohorts. The GSE192887 analysis showed significant treatment-associated expression changes in all seven genes after Pg culture-supernatant exposure. In paired tissues, CXCL8 and TMEM92 were significantly higher in OSCC tissues, whereas STC1 was not significant after Holm correction. In CAL-27 cells, CXCL8, STC1, and TMEM92 increased significantly after culture-supernatant treatment, whereas the corresponding comparisons were not significant in HOK or HSC-3 cells after adjustment. CONCLUSIONS: This study developed a seven-gene Pg-associated prognostic signature for OSCC and provided complementary transcriptomic, immune-correlation, tissue, and cell-based evidence that placed the signature in biological context. The model showed limited-to-moderate discrimination and is not ready for clinical use. The enrichment, gene-correlation, and immune-correlation findings are hypothesis-generating rather than mechanistic evidence. Further independent validation and dedicated functional studies are required.
Oral squamous cell carcinoma (OSCC)↗