Abstract
Introduction:
Diabetes mellitus (DM), a globally prevalent chronic metabolic disorder characterized by persistent hyperglycemia, has been increasingly linked to dysbiosis of the oral microbiome. However, the relationship between the virome, a crucial component of the oral microbiome, and DM remains poorly understood.
Methods:
To explore the characteristics of the oral virome in DM patients, we analyze the oral viral communities of 45 DM patients and 40 healthy controls (HC) using a publicly available metagenomic dataset.
Results:
Our analysis revealed no significant differences in a-diversity between DM patients and HC. However, Podovirus was enriched in DM patients, whereas Microviridae was more prevalent in HC. A total of 1,131 virus signal was identified, primarily belonging to the Siphovirus and Myovirus taxa. Notably, HC-enriched vOTUs exhibited broader host tropism, predominantly infecting Prevotella, Fusobacterium, and Gemella, whereas DM-enriched vOTUs showed narrower specificity for Pauljensenia and Veillonella. Cross-kingdom network analysis suggested that certain viruses (HMP_1157.k81_309051) may have potential links to the development of DM, and the bacteria genus F0040 might play a significant role in maintaining oral health. Additionally, the random forest model based on viral markers effectively distinguished between HC and DM patients (AUC =90.8%), significantly outperforming the bacterial model.
Discussion:
This indicates that these unique viral markers could serve as potential targets for DM intervention. Taken together, our findings reveal distinct alterations in the oral virome of DM patients and highlight its promise as a novel diagnostic and therapeutic target in metabolic disease research.
Introduction
Diabetes mellitus (DM) not only compromises patients’ quality of life but also induces long-term damage, dysfunction, and failure in ocular, renal, neural, cardiac, and vascular systems (). With a global prevalence exceeding 820 million adults and projections estimating a rise to 1.31 billion by 2050, DM has emerged as a major chronic disease threatening public health (). While its pathogenesis is influenced by both environmental and genetic factors, recent studies highlight the significant role of the gut microbiota in diabetes onset and progression.
However, the oral microbiota, as the second-largest microbial community in the body, also plays a crucial role in health maintenance (Sabharwal et al., 2019; ). In diabetic patients, the oral cavity shows increased levels of Lactobacillus, Corynebacterium, and Pseudomonas, while Porphyromonas, Treponema, Prevotella, and Parvimonas are reduced (). Additionally, the prevalence of Candida albicans-induced oral candidiasis in diabetic patients is as high as 30%, further underscoring microbial dysbiosis (). Although research traditionally focused on oral bacteria and fungi, yet oral viruses, despite being a key component of the oral microbiome, have received limited attention. Investigating the role of oral viruses in DM pathogenesis may provide novel insights into disease mechanisms and therapeutic opportunities.
The oral virome encompasses the entire viral community in the oral cavity, including bacteriophages, eukaryotic viruses, and free viral genetic material (). Metagenomic analyses indicate that bacteriophages constitute the dominant viral population in healthy oral ecosystems (; ). Under specific conditions, Mechanistically, these viral agents might promote pro-inflammatory cytokine production (e.g., IFN-γ, IL - 12, and TNF-γ) and oxidative stress, initiating inflammatory cascades that worsen systemic inflammation in diabetic individuals and potentially impair pancreatic islet function (Wang et al., 2018; ; ). Beyond direct inflammatory effects, oral viruses may contribute to diabetes pathogenesis by altering microbial communities in both oral and gastrointestinal environments. Currently, research on the oral virome in DM patients remains limited. Considering the documented alterations in oral bacteriome composition among DM patients and the intricate virome-bacteriome interplay, comprehensive investigation of the oral virome in DM holds substantial scientific merit.
Recent technological advancements, particularly in high-throughput sequencing, have significantly enhanced our ability to characterize the complexity of the human virome. In light of growing interest in the DM virome, we analyzed publicly available deep shotgun metagenomic sequencing data from 45 DM patients and 40 HC to characterize the oral virome in DM (). Through comprehensive reanalysis, we recovered over 35,000 non-redundant viral sequences and performed an in-depth examination of the oral virome. Our findings reveal substantial differences in oral viral community composition between DM patients and HC, with 1,131 non-redundant viruses showing significant abundance variations between the two groups. Additionally, our findings highlight complex interactions between oral viruses and bacterial populations, with the classification potential of these signatures evaluated in distinguishing DM patients from HC. This research provides new insights into potential mechanisms underlying DM pathogenesis and may help identify novel diagnostic or therapeutic targets.
Materials and methods
Participants and data summary
In this study, we reanalyzed 85 oral metagenomic sequencing samples obtained from a publicly available cohort (Supplementary Table S1), which were downloaded from the European Bioinformatics Institute database under the accession code PRJNA289586. All diagnosed DM cases had a duration exceeding 5 years and were successfully managed with insulin replacement therapy, while the HC group exhibited no evidence of imminent DM risk based on autoantibody status.
Raw reads were preprocessed using fastp (v0.20.1) with the following parameters: -l 90 -q 20 -u 30 -y –trim_poly_g, to eliminate low-quality sequences. Quality-filtered reads were aligned to the human reference genome GRCh38 using bowtie2 (v2.4.1), and host-derived contaminants were subsequently identified and removed using default parameters.
Oral virome database
To construct the oral virome reference database, we aggregated 2,792 publicly available oral metagenomic datasets from 16 global studies, encompassing samples collected from 12 distinct anatomical sites. These datasets were uniformly processed via a standardized bioinformatics pipeline, yielding 9.3 Tb of high-quality non-human metagenomic data. This pipeline culminated in an oral virome database (OVD) comprising 48,425 non-redundant vOTUs. The pipeline involved several key steps: (1) raw reads were quality-filtered and decontaminated from human sequences; (2) clean reads were assembled into contigs using SPAdes; (3) viral sequences were identified using CheckV (), DeepVirFinder (), and VIBRANT (), and decontaminated by removing those with high BUSCO ratios; and (4) viral sequences were clustered into viral operational taxonomic units (vOTUs) at 95% nucleotide identity across 75% of the sequence, with the longest sequence in each cluster serving as the representative. Specific data processing procedures were referenced from previous studies ().
Metagenome-based oral virome profiling and analyses
To profile the oral viral community in fecal metagenomes of both DM patients and HCs, we utilized Bowtie 2 (v2.4.1) () to align high-quality reads of all samples into the OVD. We set a nucleotide similarity threshold of 95% to define viral “species-level” taxa, following an established phylogenetic threshold (). By doing so, the abundance profile of vOTUs for each oral sample was produced by aggregating the number of reads mapped to each vOTU, and the relative abundances were generated by dividing it by the total mapped reads in that sample. We also produced a relative abundance profile at the viral family level by aggregating the relative abundances of vOTUs assigned to the same family.
We measured oral viral community richness using the observed number of vOTUs per family, while Shannon and Simpson indices were applied to assess virome diversity across all samples. These measures were computed using the vegan package () in the R platform, with a consistent number of reads (20 million) for each sample.
Oral bacteriome profiling
The oral microbial genome catalog was constructed from 3,569 prokaryotic species, encompassing over 50,000 metagenome-assembled genomes (MAGs), as reported in a previous study (Zhu et al., 2022). To analyze the bacterial communities in the oral metagenomes of DM patients and HC, high-quality reads from all samples were aligned against these MAGs using Bowtie 2 (v2.4.1). Phylum- and genus-level abundances were calculated through hierarchical aggregation of the proportions of constituent species.
Construction of cross-kingdom co-occurrence network
Bacterial-vOTU interactions were analyzed using Spearman’s rank correlation with a stringent threshold (| r | > 0.8). Significant associations were visualized as interaction networks through Cytoscape v3.8.2 (Su et al., 2014).
Statistical analysis
All statistical analyses in this study were conducted on the R (v4.2.1) platform. Principal coordinates analysis (PCoA) of the Bray-Curtis distance was carried out and visualized using the vegan package (). Permutational multivariate analysis of variance (PERMANOVA) was realized with the adonis function of the vegan package, with an adonis p-value generated from 1,000 permutations. We employed the Student’s t-test and Wilcoxon rank-sum test to assess statistical differences in the diversity and taxonomic levels between the two cohorts, respectively.
Comparison of functional differences
In view of the functional differences between DM-enriched and HS-enriched vOTUs, the occurrence frequency of each KO in each group was first statistically analyzed. To test whether there is a significant difference in KO frequency between the two groups, Fisher’s exact test is performed through the fisher.test function. If the corrected p value is less than 0.05, it is determined that there is a significant functional difference of this KO between the two groups.
Random forest modeling
To identify differential viral signatures distinguishing DM patients from HCs, we trained a Random Forest (RF) classifier using the randomForest package with 1,000 decision trees. RF is an ensemble learning algorithm that builds multiple decision trees and outputs either the majority vote (for classification) or the average prediction (for regression). In this study, the RF model was trained on high-dimensional, sparse vOTU abundance data, leveraging its robustness to such data structures. Feature importance was assessed via mean decrease in accuracy and Gini index, enabling implicit feature selection. To prevent overfitting, bootstrap aggregation (bagging) and random feature selection were employed. Model performance was evaluated using leave-one-out cross-validation (LOOCV), which iteratively trains the model on all samples except one and tests on the held-out sample, ensuring a rigorous assessment of predictive accuracy. We considered alternative models (for example, support vector machine, LASSO regression) but did not adopt them; RF is given priority due to its outstanding performance in dealing with nonlinear interactions and noise metagenomic characteristics.
Results
Diversity and structure of the oral virome in DM patients
To assess the impact of disease status on oral virome richness and evenness, we employed rarefaction curve analysis to quantify viral richness, measured as the number of observed vOTUs. Consistent with prior studies, no significant difference in viral richness was observed between DM patients and HC at equivalent sequencing depths (Figure 1A). Similarly, α-diversity metrics—including the Shannon, Simpson, and richness indices—showed no significant differences between the groups (Wilcoxon rank-sum test, p > 0.05; Figures 1B–D). For β-diversity analysis, we utilized Bray-Curtis dissimilarities at the vOTU level, incorporating PCoA and PERMANOVA. A significant separation between the DM and HC groups was observed (PERMANOVA: R² = 0.019, p = 0.042; Figure 1E). Taxonomic profiling at the family level revealed that annotated viral families constituted 29.31% of the total viral families. Specifically, the oral virome of DM patients was dominated by Siphovirus (mean relative abundance: 20.4%), Myovirus (7.2%), and Podovirus (0.86%). In contrast, HC individuals exhibited dominance of Siphovirus (20.16%), Myovirus (6.34%) and Podovirus (0.21%; Figure 1F). Notably, Podovirus (p = 0.00009) and Microviridae (p = 0.025) showed statistically significant abundance differences between the cohorts (Wilcoxon rank-sum test; Figure 1F; Supplementary Table S2). Podovirus was enriched in DM patients, whereas Microviridae was more prevalent in HC groups. These findings suggest that dysbiosis in the oral virome may be linked to diabetes pathogenesis (Figure 1G).
Figure 1
Identification of oral viral signatures associated with DM
Comparative analysis of the oral virome between DM patients and HC identified 1,131 vOTUs with significantly differential abundance (Wilcoxon rank-sum test, p < 0.05; Figure 2A; Supplementary Tables S3). Among these, 361 vOTUs were enriched in the DM patients, primarily affiliated with Siphovirus (37.1%), Myovirus (8%) and Herelleviridae (0.6%). In contrast, 770 HC-enriched vOTUs spanned seven families, dominated by Siphovirus (16.4%), Myovirus (8.8%) and Quimbyviridae (Figures 2B, C). We constructed a virus-host relationship based on genome homology alignment and CRISPR spacer sequence matching using established methods. The analysis revealed that DM-enriched vOTUs and HC-enriched vOTUs were primarily hosted by Streptococcus (>40%), followed by Prevotella, Pauljensenia, and others. However, we observed that HC-enriched vOTUs exhibited a stronger association with Streptococcus, Prevotella, Fusobacterium and Gemella species, while DM-enriched vOTUs were more closely linked to Pauljensenia and Veillonella bacteria. Notably, significant differences in vOTU signals were identified among DM patients (Figure 2D). Immediately after, we conducted a functional analysis of the differentiated votus and found that the enrichment patterns of different functional pathways in the DM group and the HC group were significantly different (Figure 2E; Supplementary Table S4). The proportion of K01447 et al. in DM-enriched was close to 20%, which was much higher than that in HC-enriched (about 7%), indicating that this functional pathway was highly enriched in the diabetes group and might be closely related to the pathophysiological process of diabetes. The enrichment of K17733 et al. was more prominent in the HC group, reflecting that this functional pathway was more active in the healthy control group and perhaps an important pathway for maintaining normal physiological states.
Figure 2
Oral cross-kingdom co-occurrence network analysis in DM patients
Using an analytical approach analogous to our virome methodology, we identified 195 differential bacterial species between DM patients and HCs. To investigate cross-kingdom interactions between DM-associated oral viral and bacterial signatures, we performed Spearman’s rank correlation analysis on 1,131 viral signatures and 195 bacterial signatures (Supplementary Table S5). Applying a stringent threshold (|r| > 0.8, p < 0.05), we identified 189 significant virus-bacteria associations in DM patients (Figure 3A) and 209 in HC (Figure 3B). The two networks exhibited distinct interaction patterns. Specifically, the DM network was predominantly composed of positive correlations, indicative of a relatively unstable network structure (). Further analysis of network hub nodes revealed that, in the DM network, a specific viral species (HMP_1157.k81_309051) played a central role, influencing multiple bacterial taxa (e.g., Lancefieldella, Pauljensenia, Figure 3C). In contrast, the HC network highlighted the potential importance of the bacterial genus F0040 (Figure 3D), which appeared to play a key role in maintaining healthy oral interactions. Notably, this relationship was decrease in the DM network, suggesting that F0040 may be crucial for sustaining healthy oral microbial interactions.
Figure 3
Diagnostic potential of oral virome and bacteriome in DM
To evaluate the diagnostic potential of oral microbiome features for DM classification, we trained random forest models using three distinct feature sets: bacteriome, virome, and a combined (bacteriome-virome) approach. The virome model demonstrated the highest diagnostic accuracy, achieving an area under the receiver operating characteristic curve (AUC) of 90.8% (95% CI: 84.1 – 97.5%). In comparison, the bacteriome model yielded a lower AUC of 80.1% (95% CI: 70.5 – 89.8%), while the combined model exhibited performance comparable to the virome model (AUC = 90.4%, 95% CI: 83.6 – 97.2%) (Figure 4A). Furthermore, by analyzing feature importance across these models, we consistently observed that viral features held greater significance in the diagnostic framework (Figures 4B–D). Collectively, these findings highlight the potential of oral viral signatures as valuable biomarkers for DM diagnosis, offering a promising avenue for future microbiome-based diagnostic strategies.
Figure 4
Discussion
The global rise in DM, driven by aging populations and environmental factors, has established it as a leading endocrine-metabolic disorder worldwide. While a strong bidirectional relationship between oral microbial dysbiosis and DM pathogenesis has been documented (Sanz et al., 2018; Rohani, 2019), the role of the virome in DM remains poorly characterized. In this study, we analyzed the oral virome of a cohort comprising 85 individuals (DM patients vs. HC group) to evaluate viral diversity, taxonomic composition, and cross-kingdom network shifts. Specifically, although no significant differences in viral α-diversity were observed between groups, a large number of viral signatures revealed significant alterations in the oral virome of DM patients, demonstrating strong predictive performance in disease diagnosis. Furthermore, cross-kingdom network analysis demonstrated unique microbial interactions in the oral microbiome of DM patients.
At the family level, compositional profiling revealed significant enrichment of Podovirus in the DM group, which aligns with previous findings in severe early childhood caries (). In diabetes patients, oral inflammation (e.g., periodontitis) leads to suboptimal glycemic control, a higher incidence of diabetes-related complications, and conversely, poor blood sugar management exacerbates oral microbial pathogenicity, forming a vicious cycle (; ). By contrast, Microviridae was more abundant in the HC group. Studies have shown that Microviridae are small phages that infect Gram-negative commensal bacteria. They may protect microenvironment homeostasis by modulating bacterial population balance, suppressing overgrowth of opportunistic pathogens, and sustaining the metabolic activity of short-chain fatty acid-producing bacteria (Ríos-Covián et al., 2016; ; ). The distinct enrichment patterns of Podovirus and Microviridae suggest their divergent ecological roles in the diabetic oral niche. Podovirus primarily infecting Gram-positive bacteria, likely target key oral pathogens such as Streptococcus and Actinomyces ()– consistent with our observations of DM-enriched Streptococcus and Pauljensenia (Actinomyces) phages. In diabetic individuals, chronic hyperglycemia may promote lytic activation of these phages, exacerbating tissue damage through bacterial lysis-induced inflammation (; Xu et al., 2022). Conversely, Microviridae’s preference for Gram-negative commensals (e.g., Prevotella, Fusobacterium) positions them as potential homeostasis modulators (). Their reduction in DM could permit overgrowth of LPS-producing bacteria, impairing epithelial barrier function and fueling metabolic endotoxemia (). These findings highlight the intricate interplay between viral communities, oral and systemic health, and metabolic regulation, emphasizing the need for further research to elucidate these complex interactions.
We identified numerous DM-associated vOTUs and observed that DM-enriched viruses were predominantly composed of Streptococcus, Pauljensenia, and Veillonella phages, whereas HC-enriched vOTUs were dominated by Prevotella, Fusobacterium, and Gemella. It is plausible that the association between Gemella-infecting bacteriophages—enriched in healthy controls—and Gemella itself may be linked to butyric acid production. This connection, although speculative, suggests a potential role for these phages in supporting intestinal barrier integrity and limiting LPS translocation through modulation of butyrate-producing pathways (; ; Zhang and Wang, 2023). K01447 (N-acetylmuramyl-L-alanine amidase) was significantly enriched in the DM group. This enzyme hydrolyzes bacterial cell wall peptidoglycan to release muramyl-dipeptide (MDP), which activates the host’s innate immune receptor NOD2 and drives chronic inflammatory responses (). Importantly, current evidence suggests that this process may operate synergistically with Veillonella phage-mediated excessive Th17 activation. This potential synergy could contribute to IL - 17A/IL-22 releases, which may further exacerbate insulin resistance and systemic inflammation (; Talukdar et al., 2012; Smith et al., 2013). K17733 (peptidoglycan LD-endopeptidase CwlK) was highly active in the HC group. It functions to maintain the structural stability of the flora by precisely regulating peptidoglycan hydrolysis.
Microbial interaction network analysis in this study further revealed that the restructuring of the diabetic oral ecosystem is characterized not only by shifts in viral abundance but also by a dynamic imbalance in virus–bacteria synergistic relationships. Although the total number of co-abundance relationships was similar between the HC and DM groups, the feature of viral nodes occupying key hub positions was particularly prominent in the DM group. In the bacterial-viral correlation analysis of the oral microbiome, the virus HMP_0806.k81_47237_1 maintained high connectivity (degree > 20) in both the HC and DM groups, indicating its critical role in maintaining oral homeostasis. By contrast, the association strength of HMP_1157.k81_309051 increased significantly in DM, suggesting that under specific environmental pressures (e.g., DM microenvironment), this virus may exacerbate microbial dysbiosis and host inflammation. Notably, F0040 (bacterial genus) played a central role as a core species in the HC network, yet its potential function in the oral cavity remains underexplored in diabetes-related studies, despite its observed associations with Veillonella and Streptococcus in other research contexts (). Moreover, in the construction of diagnostic models, viral community features outperformed bacterial features in predictive accuracy. This challenges the conventional bacteria-centric paradigm in microbiome-based diagnostics and highlights the potential of the oral virome as a novel biomarker and early warning indicator for diabetes. These findings offer new insights and directions for the early detection and intervention of diabetes through virome-based approaches. It is worth noting that the current interaction network is constructed based on Spearman correlations, which capture statistical associations but do not imply causality. Although this method is widely used to reveal co-abundance patterns in microbiome research, it does not reflect the directionality or underlying mechanisms of virus–bacteria interactions. Therefore, future studies should consider incorporating longitudinal sampling and causal inference techniques—such as Granger causality, structural equation modeling, or dynamic Bayesian networks—to more precisely elucidate the direction and impact of these microbial interactions under disease conditions.
However, this study has limitations that warrant consideration. First, the moderate sample size (n = 85) may constrain the statistical power and generalizability of our findings. Although compelling, these results necessitate validation in larger, multi-center cohorts to confirm their robustness. Second, the absence of longitudinal or prospective data limits our ability to infer causal relationships between oral virome dynamics and DM pathogenesis. Third, critical confounding variables—including periodontal status, oral hygiene habits, medication use (e.g., antibiotics or antidiabetic drugs), and detailed glycemic control indicators (e.g., HbA1c levels)—were unavailable in the original public dataset (PRJNA289586). This limitation was unavoidable due to the metadata constraints of the source data, and may bias the observed associations between virome alterations and diabetes. Additionally, current taxonomic ambiguities in oral virus classification hinder mechanistic exploration. Thus, integrating animal models is essential to elucidate potential pathogenic mechanisms linking oral virome to DM.
Conclusion
In summary, through metagenomic sequencing of oral samples from DM patients, we characterized the oral virome alterations in DM, identifying 1,131 DM-associated viral signatures and disrupted cross-kingdom interactions. Notably, virus-based classifiers achieved high accuracy (AUC = 90.8%) in predicting DM status, highlighting translational potential. Critically, our finding that DM-associated virome alterations center on bacteriophage families (Podovirus/Microviridae) with opposing impacts on bacterial dysregulation reveals an underappreciated viral layer in diabetes pathology. Nevertheless, the moderate cohort size (n = 85) and lack of external validation require caution in generalizing these findings. Future studies must verify these results in independent, expanded cohorts and address causal links through longitudinal designs. Our work provides a foundation for mechanistic and clinical studies on oral virome-DM interplay.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
YiZ: Data curation, Formal Analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing. YuZ: Data curation, Formal Analysis, Investigation, Software, Visualization, Writing – review & editing. GX: Data curation, Formal Analysis, Writing – review & editing. TM: Writing – original draft. MW: Data curation, Writing – review & editing. CH: Writing – original draft. HY: Investigation, Writing – review & editing. YZh: Writing – original draft. SY: Investigation, Writing – review & editing. QY: Conceptualization, Writing – review & editing, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization. SL: Conceptualization, Writing – review & editing. CC: Conceptualization, Funding acquisition, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by National Natural Science Foundation of China (82370563), the China Postdoctoral Science Foundation (2024MD754009), Young Top Talents Program of the Liaoning Revitalization Talents Initiative (XLYC2403207), Liaoning Provincial Natural Science Foundation Joint Project (2023-MSLH-032), Basic Research Project of Liaoning Provincial Department of Education for Universities (LJ212410161043), Dalian Municipal Guidance Program for the Life and Health Sector (2024ZDJH01PT130), and the Key Laboratory of Guizhou Provincial Education Department (Guizhou Education Technology [2023] No.017).
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.
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/fcimb.2025.1607798/full#supplementary-material
References
1
AghoE. T.OwotadeF. J.KolawoleB. A.OyetolaE. O.AdedejiT. A. (2021). Salivary inflammatory biomarkers and glycated haemoglobin among patients with type 2 diabetic mellitus. BMC Oral. Health21, 101. doi: 10.1186/s12903-021-01453-y
2
American Diabetes Association. (2011). Diagnosis and classification of diabetes mellitus. Diabetes Care34 (Suppl 1), S62–S69. doi: 10.2337/dc11-S062
3
ArpaiaN.CampbellC.FanX.DikiyS.van der VeekenJ.deRoosP.et al. (2013). Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature504, 451–455. doi: 10.1038/nature12726
4
BakerJ. L.Mark WelchJ. L.KauffmanK. M.McLeanJ. S.HeX. (2024). The oral microbiome: diversity, biogeography and human health. Nat. Rev. Microbiol.22, 89–104. doi: 10.1038/s41579-023-00963-6
5
ByndlossM. X.OlsanE. E.Rivera-ChávezF.TiffanyC. R.CevallosS. A.LokkenK. L.et al. (2017). Microbiota-activated PPAR-γ signaling inhibits dysbiotic Enterobacteriaceae expansion. Science357, 570–575. doi: 10.1126/science.aam9949
6
ChenX.ZouT.ZengQ.ChenY.ZhangC.JiangS.et al. (2024). Metagenomic analysis reveals ecological and functional signatures of oral phageome associated with severe early childhood caries. J. Dentistry146, 105059. doi: 10.1016/j.jdent.2024.105059
7
CoyteK. Z.SchluterJ.FosterK. R. (2015). The ecology of the microbiome: networks, competition, and stability. Science. 350, 663–666. doi: 10.1126/science.aad2602
8
DiaoJ.ZhangY.ZhangX.JiaS.LeiY.MaB.et al. (2025). Integrated analysis of oral microbiome and metabolome in T2DM patients with varying glycemic status. Oral Dis. 31 (7), 2308–2322. doi: 10.1111/odi.15220
9
DixonP. (2003). VEGAN, a package of R functions for community ecology. J. vegetation Sci.14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
10
FanG.CaoF.KuangT.YiH.ZhaoC.WangL.et al. (2023). Alterations in the gut virome are associated with type 2 diabetes and diabetic nephropathy. Gut Microbes15, 2226925. doi: 10.1080/19490976.2023.2226925
11
GBD 2021 Diabetes Collaborators. (2023). Global, regional, and national burden of diabetes from 1990 to 2021, with projections of prevalence to 2050: a systematic analysis for the Global Burden of Disease Study 2021. Lancet402 (10397), 203–234. doi: 10.1016/S0140-6736(23)01301-6
12
GanesanS. M.JoshiV.FellowsM.DabdoubS. M.NagarajaH. N.O’DonnellB.et al. (2017). A tale of two risks: smoking, diabetes and the subgingival microbiome. Isme J.11, 2075–2089. doi: 10.1038/ismej.2017.73
13
GiriB.DeyS.DasT.SarkarM.BanerjeeJ.DashS. K. (2018). Chronic hyperglycemia mediated physiological alteration and metabolic distortion leads to organ dysfunction, infection, cancer progression and other pathophysiological consequences: An update on glucose toxicity. BioMed. Pharmacother.107, 306–328. doi: 10.1016/j.biopha.2018.07.157
14
GogokhiaL.BuhrkeK.BellR.HoffmanB.BrownD. G.Hanke-GogokhiaC.et al. (2019). Expansion of bacteriophages is linked to aggravated intestinal inflammation and colitis. Cell Host Microbe25, 285–299.e288. doi: 10.1016/j.chom.2019.01.008
15
GrazianiF.GennaiS.SoliniA.PetriniM. (2018). A systematic review and meta-analysis of epidemiologic observational evidence on the effect of periodontitis on diabetes An update of the EFP-AAP review. J Clin. Periodonto. 45, 167–187. doi: 10.1111/jcpe.12837
16
GregoryA. C.ZayedA. A.Conceição-NetoN.TempertonB.BolducB.AlbertiA.et al. (2019). Marine DNA viral macro-and microdiversity from pole to pole. Cell177, 1109–1123.e1114. doi: 10.1016/j.cell.2019.03.040
17
Heintz-BuschartA.MayP.LacznyC. C.LebrunL. A.BelloraC.KrishnaA.et al. (2016). Integrated multi-omics of the human gut microbiome in a case study of familial type 1 diabetes. Nat. Microbiol.2, 16180. doi: 10.1038/nmicrobiol.2016.180
18
KarinaD.HeldayaniI.HidayatW. (2025). Oral opportunistic infection induced by stress and silent type 2 diabetes mellitus in young adult patient: A case report. Int. Med. Case Rep. J.18, 59–66. doi: 10.2147/IMCRJ.S488127
19
KornT.MitsdoerfferM.CroxfordA. L.AwasthiA.DardalhonV. A.GalileosG.et al. (2008). IL - 6 controls Th17 immunity in vivo by inhibiting the conversion of conventional T cells into Foxp3+ regulatory T cells. Proc. Natl. Acad. Sci. U.S.A.105, 18460–18465. doi: 10.1073/pnas.0809850105
20
KövérZ.Johansen NordskagV.BánÁ.GajdácsM.UrbánE. (2023). The role of Actinomyces spp. and related organisms in cervicofacial infections: Pathomechanism, diagnosis and therapeutic aspects. Anaerobe82, 102767. doi: 10.1016/j.anaerobe.2023.102767
21
KrupovicM.ForterreP. (2011). Microviridae goes temperate: microvirus-related proviruses reside in the genomes of bacteroidetes. PloS One6 (5), e19893. doi: 10.1371/journal.pone.0019893
22
LangmeadB. (2012). Salzberg SLJNm: Fast gapped-read alignment with Bowtie 2. Nat. Methods9, 357–359. doi: 10.1038/nmeth.1923
23
LekunberriI.SubiratsJ.BorregoC. M.BalcázarJ. L. (2017). Exploring the contribution of bacteriophages to antibiotic resistance. Environ. pollut.220, 981–984. doi: 10.1016/j.envpol.2016.11.059
24
LiS.GuoR.ZhangY.LiP.ChenF.WangX.et al. (2022). A catalog of 48,425 nonredundant viruses from oral metagenomes expands the horizon of the human oral virome. iScience25, 104418. doi: 10.1016/j.isci.2022.104418
25
LiC.ZhangY.YanQ.GuoR.ChenC.LiS.et al. (2023). Alterations in the gut virome in patients with ankylosing spondylitis. Front. Immunol.14, 1154380. doi: 10.3389/fimmu.2023.1154380
26
LiangG.BushmanF. D. (2021). The human virome: assembly, composition and host interactions. Nat. Rev. Microbiol.19, 514–527. doi: 10.1038/s41579-021-00536-5
27
LiuX.TongX.JieZ.ZhuJ.TianL.SunQ.et al. (2023). Sex differences in the oral microbiome, host traits, and their causal relationships. iScience26 (1), 105839. doi: 10.1016/j.isci.2022.105839
28
LvL. J.WenJ. Y.ZhangY.GuoR. C.LiH.YiZ. T.et al. (2024). Deep metagenomic characterization of the gut virome in pregnant women with preeclampsia. mSphere9, e0067623. doi: 10.1128/msphere.00676-23
29
MorrisonA. G.SarkarS.UmarS.LeeS. T. M.ThomasS. M. (2023). The contribution of the human oral microbiome to oral disease: A review. Microorganisms11 (2), 318. doi: 10.3390/microorganisms11020318
30
NayfachS.CamargoA. P.SchulzF.Eloe-FadroshE.RouxS.KyrpidesN. C. (2021). CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nat. Biotechnol.39, 578–585. doi: 10.1038/s41587-020-00774-7
31
NegroniA.PierdomenicoM.CucchiaraS.StronatiL. (2018). NOD2 and inflammation: current insights. J. Inflammation Res.11, 49–60. doi: 10.2147/JIR.S137606
32
RenJ.SongK.DengC.AhlgrenN. A.FuhrmanJ. A.LiY.et al. (2020). Identifying viruses from metagenomic data using deep learning. Quant Biol.8, 64–77. doi: 10.1007/s40484-019-0187-4
33
Ríos-CoviánD.Ruas-MadiedoP.MargollesA.GueimondeM.de Los Reyes-GavilánC. G.SalazarN. (2016). Intestinal short chain fatty acids and their link with diet and human health. Front. Microbiol.7, 185. doi: 10.3389/fmicb.2016.00185
34
RohaniB. (2019). Oral manifestations in patients with diabetes mellitus. World J. Diabetes10, 485–489. doi: 10.4239/wjd.v10.i9.485
35
SabharwalA.GanleyK.MiecznikowskiJ. C.HaaseE. M.BarnesV.ScannapiecoF. A. (2019). The salivary microbiome of diabetic and non-diabetic adults with periodontal disease. J. Periodontol90, 26–34. doi: 10.1002/JPER.18-0167
36
SanzM.CerielloA.BuysschaertM.ChappleI.DemmerR. T.GrazianiF.et al. (2018). Scientific evidence on the links between periodontal diseases and diabetes: Consensus report and guidelines of the joint workshop on periodontal diseases and diabetes by the International Diabetes Federation and the European Federation of Periodontology. J. Clin. Periodontol45, 138–149. doi: 10.1111/jcpe.12808
37
SmithP. M.HowittM. R.PanikovN.MichaudM.GalliniC. A.BohloolyY. M.et al. (2013). The microbial metabolites, short-chain fatty acids, regulate colonic Treg cell homeostasis. Science341, 569–573. doi: 10.1126/science.1241165
38
SuG.MorrisJ. H.DemchakB.BaderG. D. (2014). Biological network exploration with Cytoscape 3. Curr. Protoc. Bioinf.47, 8.13.1–8.13.24. doi: 10.1002/0471250953.bi0813s47
39
TalukdarS.OhD. Y.BandyopadhyayG.LiD.XuJ.McNelisJ.et al. (2012). Neutrophils mediate insulin resistance in mice fed a high-fat diet through secreted elastase. Nat. Med.18, 1407–1412. doi: 10.1038/nm.2885
40
WangJ.ZhengJ.ShiW.DuN.XuX.ZhangY.et al. (2018). Dysbiosis of maternal and neonatal microbiota associated with gestational diabetes mellitus. Gut67, 1614–1625. doi: 10.1136/gutjnl-2018-315988
41
XuH. M.XuW. M.ZhangL. (2022). Current Status of Phage Therapy against Infectious Diseases and Potential Application beyond Infectious Diseases. Int. J. Of Clin. Pract.2022, 4913146. doi: 10.1155/2022/4913146
42
ZhangY.WangR. (2023). The human gut phageome: composition, development, and alterations in disease. Front. Microbiol.14, 1213625. doi: 10.3389/fmicb.2023.1213625
43
ZhuJ.TianL.ChenP.HanM.SongL.TongX.et al. (2022). Over 50,000 metagenomically assembled draft genomes for the human oral microbiome reveal new taxa. Genom. Proteom. Bioinform.20, 246–259. doi: 10.1016/j.gpb.2021.05.001
Summary
Keywords
oral viruses, diabetes mellitus, viral diversity, viral function, host
Citation
Zhang Y, Zhang Y, Xing G, Mei T, Wang M, Huang C, Yi H, Zhan Y, Yang S, Yan Q, Li S and Chen C (2025) Characterization of the oral virome in patients with diabetes mellitus. Front. Cell. Infect. Microbiol. 15:1607798. doi: 10.3389/fcimb.2025.1607798
Received
08 April 2025
Accepted
17 August 2025
Published
03 September 2025
Volume
15 - 2025
Edited by
Velmurugan Ganesan, KMCH Research Foundation, India
Reviewed by
Dóra Tombácz, University of Szeged, Hungary
Monica Viveros-Rogel, National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico
Updates
Copyright
© 2025 Zhang, Zhang, Xing, Mei, Wang, Huang, Yi, Zhan, Yang, Yan, Li and Chen.
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: Changming Chen, dychangming_chen@163.com; Shenghui Li, lsh2@qq.com; Qiulong Yan, qiulongy1988@163.com
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.