
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)13315-4
10.1016/j.heliyon.2024.e37284
e37284
Research Article
Gender-specific microbial signatures in saliva: Unveiling the association between the oral microbiome and the pathogenesis of glioma
Qin Hao haoqin@bjmu.edu.cn
a⁎1
Liu Jie b1
Li Yang-Yang c
Xu Ya-Lan d
Yan Yi-Fang yanyifangyanyf@hotmail.com
efgh⁎⁎
a State Key Laboratory of Molecular Oncology, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
b Department of Medical Records, Air Force Medical Center, PLA, Air Force Medical University, Beijing, China
c Medical Center for Human Reproduction, Beijing Chao-Yang Hospital, Capital Medical University, Beijing, China
d Department of Obstetrics and Gynecology, Peking Union Medical College Hospital, Beijing, China
e State Key Laboratory of Female Fertility Promotion, Center for Reproductive Medicine, Department of Obstetrics and Gynecology, Peking University Third Hospital, Beijing, China
f National Clinical Research Center for Obstetrics and Gynecology (Peking University Third Hospital), Beijing, China
g Key Laboratory of Assisted Reproduction (Peking University), Ministry of Education, Beijing, China
h Beijing Key Laboratory of Reproductive Endocrinology and Assisted Reproductive Technology, Beijing, China
⁎ Corresponding author. State Key Laboratory of Molecular Oncology, Cancer Institute Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, No. 17 Panjiayuan Nanli, Chaoyang District, Beijing, 100021, China. haoqin@bjmu.edu.cn
⁎⁎ Corresponding author. State Key Laboratory of Female Fertility Promotion, Center for Reproductive Medicine, Department of Obstetrics and Gynecology, Peking University Third Hospital, No. 49 North Huayuan Road, Haidian District, Beijing, 100191, China. yanyifangyanyf@hotmail.com
1 These authors contributed equally to this work.

31 8 2024
15 9 2024
31 8 2024
10 17 e3728421 3 2024
29 8 2024
30 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
The intricate interplay between the human oral microbiome and systemic health is increasingly being recognized, particularly in the context of central nervous system pathologies such as glioblastoma. In this study, we aimed to elucidate gender-specific differences in the salivary microbiome of glioma patients by utilizing 16S rRNA sequencing data from publicly available salivary microbiome datasets. We conducted comprehensive bioinformatics analysis, encompassing quality control, noise reduction, species classification, and microbial community composition analysis at various taxonomic levels. Machine learning algorithms were employed to identify microbial signatures associated with glioma. When compared to healthy controls, our analysis revealed distinct differences in the salivary microbiota of glioma patients. Notably, the genera Leptotrichia and Atopobium exhibited significant variations in abundance between genders. Leptotrichia was prevalent in healthy females but exhibited a reduced abundance in female glioma patients. In contrast, Atopobium was more abundant in male glioma patients. These findings suggest that hormonal influences might play a role in shaping the salivary microbiome and its association with glioma. We utilized a combination of LASSO-logistic regression and random forest models for feature selection, and identified key microbial features that differentiated glioma patients from healthy controls. We developed a diagnostic model with high predictive accuracy and area under the curve and principal component analysis metrics confirmed its robustness. The analysis of microbial markers, including Atopobium and Leptotrichia, highlighted the potential of the salivary microbiota as a non-invasive biomarker for the diagnosis and prognosis of glioma. Our findings highlight significant gender-specific disparities in the salivary microbiome of patients with glioma, offering new insights into the pathogenesis of glioma and paving the way for innovative diagnostic and therapeutic strategies. The use of saliva as a diagnostic fluid, given its ease of collection and non-invasive nature, holds immense promise for monitoring systemic health and the trajectory of disease. Future research should focus on investigating the underlying mechanisms by which the salivary microbiome influences the development of glioma and identifying potential microbiome-targeted therapies to enhance the management of glioma.

Keywords

Salivary flora
Glioma
Machine learning
Disease markers
Risk model
==== Body
pmc1 Introduction

A diverse array of microorganisms, predominantly bacteria, establishes a critical ecological balance in the human oral ecosystem. The integrity of this microbial consortium is intricately linked to the etiology of multiple pathologies. Therefore, delineating the symbiotic interactions between the host and its resident oral microbiota is imperative if we are to unravel the complexities underlying the pathogenesis of oral disease [1,2]. Daily salivary output, which fluctuates between 0.75 and 1.5 L, can exert significant impact on the gastrointestinal microbiota via the ingestion of saliva and accompanying nutrients. Perturbations in the equilibrium of the oral microbiota have been implicated in a broad spectrum of systemic diseases, thereby enhancing individual susceptibility to these disorders [3]. Emerging evidence has also identified a bidirectional pathway between salivary components and the brain, indicating that substances in the saliva may reach the brain through sublingual routes. This fascinating discovery opens new avenues for utilizing saliva in the diagnosis of brain diseases [4]. Given its ease of collection and non-invasive nature, saliva is posited as a paramount diagnostic fluid with immense potential and benefits for monitoring systemic health and the trajectory of diseases. Consequently, saliva is heralded as an indispensable biomarker for the diagnosis, prognostication, and therapeutic guidance of various illnesses, including neurological disorders [5]. Glioblastoma multiforme stands as the most common primary malignant tumor of the central nervous system in adults, and accounts for approximately 58.4 % of all gliomas [6]. This disease is associated with a poor five-year survival rate and high rates of mortality, with a median survival time of 14–15 months despite the application of various treatment modalities, including surgery, chemotherapy, and radiotherapy [7]. Consequently, the discovery of effective therapeutic targets for glioblastoma is imperative to ameliorate treatment outcomes for patients. Compounding the complexity of glioblastoma multiforme treatment is the inherent heterogeneity observed between different patients and within tumors themselves, in addition to the presence of multiple types of tumor stem cells and the variability of the tumor microenvironment. Collectively, these factors constitute significant obstacles to the diagnosis and treatment of glioblastoma multiforme, severely hindering therapeutic advancements [8]. Recent studies have highlighted the promising role of salivary biomarkers for monitoring therapeutic responses and disease progression in patients with glioma. For instance, advancements in understanding the potential utility of saliva-derived exosomes before and after treatment have been reported for patients with glioma. These findings suggest that salivary biomarkers could play a pivotal role in the management of glioma [9]. Furthermore, the complex interplay between the salivary microbiota and severity of glioma has gained significant attention. For example, research by Wen and coworkers identified a profound association between the composition of the salivary microbiota and the grade of glioma. Specifically, five distinct microbial taxa were identified as markers that possessed the ability to differentiate patients with high-grade gliomas from healthy controls. Notably, the abundance of Patescibacteria was inversely correlated with the malignancy of gliomas, thus serving as a negative risk factor [10]. This discovery highlighted the potential of the salivary microbiota for the diagnosis and prognosis of glioma, thus emphasizing the critical role of microbial research in understanding the pathogenesis of glioma.

Recent studies have indicated significant gender differences in the incidence and prognosis of glioma, with males exhibiting a higher incidence rate and females experiencing better outcomes; these differences are probably related to differences in treatment responses that are medicated by sex hormones [11]. However, there is a significant lack of research focusing on differences in the salivary microbiome between male and female patients with glioma. This gap in knowledge was the prime motivation for our investigation into the salivary microbial communities of glioma patients, which utilized bioinformatics analysis to uncover potential sex-specific microbial signatures that could influence the pathology and treatment outcomes of glioma. Our bioinformatics approach included quality control, noise reduction, and species classification of sequencing data, followed by a comprehensive analysis of microbial community composition at various taxonomic levels. In addition, we employed machine learning algorithms for feature selection within the salivary microbiome in an attempt to identify signature microbial communities associated with glioma. Our findings, including the identification of salivary microbial taxa that exhibited significant differences in abundance between genders, highlight the complexity of the interaction between the salivary microbiome and glioma. These insights could pave the way for innovative diagnostic and therapeutic strategies by exploiting the non-invasive nature of saliva collection to enhance the management of glioma.

2 Materials and methods

2.1 Bioinformatics analysis of salivary microbial communities

Salivary microbiota profiling was conducted using 16S rRNA sequencing with BioProject accession number PRJNA750937, as described in the previous study by Wen et al. [10]. The raw sequencing sequences were subjected to quality control using Fastp software [12]. Following this, noise reduction and other processing steps were performed on the sequences using the DADA2 plugin within the Qiime2 workflow [13]. The sequences processed with DADA2 were named amplicon sequence variants (ASVs) and stored in qza format. To obtain species classification information for each ASV, we utilized the Greengenes2 reference database (version 13.8) [14]. Based on the principles of machine learning, the ASVs were taxonomically annotated using the classify-sklearn method within the q2-feature-classifier plugin. Furthermore, the community composition of each sample was statistically analyzed at various taxonomic levels, including kingdom, phylum, class, order, family, genus, and species, thereby yielding taxonomic annotation. A dilution curve can directly reflect the appropriateness of sequencing data volume and indirectly indicates the species richness within samples. By utilizing the “qiime2R″ and “microeco” tools in R software for graphical representation, we constructed a dilution curve by simulating the resampling process. This allowed for the observation of trends in species variation and the assessment of species richness within the environment. Next, the curve was used to compare the richness of species across samples with varying volumes of sequencing data, thus serving to determine the adequacy of the sequencing data volume for each sample.

2.2 Screening for salivary microbiome signature communities, model validation and evaluation

At the genus level, original ASVs were randomly divided into training and validation sets in a 1:1 ratio. Feature selection within the training set was conducted by adopting two approaches. First, the “microeco” package (version 1.2.0) was employed to select characteristic microbial communities via a random forest model algorithm [15]. Second, the “glmnet” package (version 4.1-7) in R, utilizing a 10-fold cross-validation with the family set to binomial, was used to identify model microbial communities and coefficients based on differential microbes [16]. The feature microbes identified by the two algorithms were then compared and the intersecting components were retained as optimized feature communities for further analysis.

2.3 Identification and analysis of salivary microbial diversity by community bar charts, linear discriminant analysis (LDA) effect size (LEfSe), and network analysis

Data analysis was conducted using the Majorbio Cloud Platform (https://www.majorbio.com). Microbiological sequencing data from saliva, sourced from the BioProject database under the accession number PRJNA750937, were acquired as reported previously by Wen et al. [10]. The present study included a total of 47 saliva microbiome samples: 10 healthy males, 14 healthy females, 13 males with high-grade gliomas, and 10 females with high-grade gliomas. Profiling of the microbiota was conducted by 16S rRNA sequencing, as described by Wen et al. [10]. By utilizing LEfSe analysis (http://huttenhower.sph.harvard.edu/LEfSe) [17], we subsequently identified salivary microbial taxa that exhibited significant differential abundances across groups, from domain to the genus level (P < 0.05, LDA score >2.0). We also used the Kruskal-Wallis rank-sum test, a non-parametric factorial approach, to identify microbial features that exhibited significant differences in abundance. Next, we used the LDA method to evaluate the magnitude of the effect size relating to each of the identified features. Inter-species correlations were calculated by adopting a network analysis approach, which allowed us to construct a correlation network diagram. This analysis was performed at the genus taxonomic level and used Spearman's rank correlation coefficient to quantify the relationships between species at the genus level. Only correlations with an absolute value of ≥0.5 and a significance level of P < 0.05 were considered significant. In the resultant network diagrams, the thickness of the lines is proportional to the strength of the correlation between species, whereas the abundance of lines associated with a species represented the level of connectivity to other species within the network. In addition, we performed community bar chart analysis in R software (version 3.3.1) to identify dominant species across various taxonomic levels (phylum, class, order, family, genus, and species). This analysis was performed by representing each species based on its relative abundance (percentage contribution).

3 Results

3.1 PLS-DA and LEfSe analysis and the detection of variations in the salivary microbiota of patients with glioma and healthy controls

First, we employed the PLS-DA method to group samples so that we could intuitively investigate the distribution of salivary microbiota in patients with glioma compared to healthy controls. This analysis revealed the differences and similarities in the composition of the salivary microbiome (Fig. 1A). Subsequently, we conducted multi-level differential analysis using the LEfSe method (LDA threshold = 2.0) to investigate differences in the salivary microbiota between patients with glioma and healthy controls across various taxonomic levels, from domain to genus (Fig. 1B). Moreover, community composition analysis and bar graph presentations were used to illustrate trends in microbial variation between different groups, in which the y-axis represents the relative proportion of microbial genera in both groups, and the length of the bars indicates the abundance of species (Fig. 1C). These bar graphs also display LDA scores for different species, with higher LDA scores signifying a greater impact of species abundance on the inter-group differences. The results indicate that salivary bacteria, such as f_Porphyromonadaceae and g_Porphyromonas, play a key role in healthy individuals, while bacteria such as f_Lactobacillaceae and g_Atopobium are more abundant in patients with glioma (Fig. 1D). Next, we used network analysis to detect correlations between microbial species in healthy individuals and patients with glioma, and to analyze interactions among dominant species. Our analysis showed that the dominant species in healthy individuals were g_Prevotella_7, g_Prevotella, g_Streptococcus, g_Porphyromonas, and g_Veillonella. In patients with glioma, the dominant species were g_Haemophilus, g_Gemella, g_Neisseria, g_Prevotella_7, g_Granulicatella, g_Porphyromonas, g_Actinomyces, g_Veillonella, and g_Prevotella (Fig. 1E). Subsequently, we analyzed rarefaction curves by simulating the resampling process; this analysis revealed that curves became more level as the volume of sequencing data increased. This trend suggested that the sequencing data was approaching saturation, thus providing a more accurate reflection of species richness within samples (Supplementary Figure 1).

3.2 Gender-specific differences in the salivary microbiome between patients with glioma and healthy controls, with unique changes in microbial abundance

By applying a PLS-DA regression model, we detected notable differences in compositions of the salivary microbiome between patients with glioma (for both genders) and healthy controls (Fig. 2A and B). These differences were highlighted by community bar plots which depicted the top 30 abundant genera in terms of their relative abundance, with “others” representing species with lower abundance. These results indicated that the Streptococcus, Prevotella_7, and Veillonella genera predominated in the salivary microbiome. Specifically, the proportion of Streptococcus was higher in female glioma patients when compared to healthy females, whereas the Prevotella_7 genus was more abundant in male glioma patients than in healthy males, thus indicating distinct microbial compositions in the saliva of patients with glioma when compared to healthy individuals, with unique changes in abundance detected between male and female patients (Fig. 2C and D). In addition, we investigated microbial changes at the phylum, class, order, family, and species levels. Results indicated that no significant change in the abundance of Gammaproteobacteria was detected between female glioma patients and healthy females, while the abundance of this microorganism was reduced in male glioma patients when compared to healthy males. In contrast, there was an increased abundance of the Prevotellaceae family in male glioma patients when compared to healthy males; in contrast, no significant difference was detected between female glioma patients and healthy females (Supplementary Figure 2). In addition, Circos software (version 0.67-7, http://circos.ca/) was used to generate Circos plots to illustrate the distribution of salivary microbial compositions between different genders of glioma patients and healthy individuals. These plots confirmed our previous bar plots, showing similar compositions of dominant microbes across different groups. Streptococcus was the dominant genus, followed by Prevotella_7 as the second most dominant. In the salivary microbiome of healthy males and male glioma patients, the abundance of Prevotella_7 increased from 36 % to 64 %, but remained consistent at 48 % and 52 % in healthy females and female glioma patients, respectively (Fig. 2E and F). Similar patterns were detected across various taxonomic levels, including phylum, class, order, family, and species. Taking the family level as an example, the relative abundance of Prevotellaceae in the salivary microbiota of healthy women was 48 %, compared to 52 % in women with glioma. In contrast, the relative abundance of Prevotellaceae was 39 % in healthy men but was more abundant (61 %) in men with glioma (Supplementary Figure 3). In addition, our analyses identified gender-specific variations in the salivary microbiota of healthy individuals. These differences involved specific bacterial groups, including Eubacterium nodatum, Solobacterium, and certain unclassified members within the family (Supplementary Figure 4). These findings demonstrated significant differences in the compositions of the salivary microbiome between male and female glioma patients when compared to healthy individuals, thus indicating gender-specific changes in microbial abundance.Fig. 1 Comprehensive analysis and visualization of the salivary microbiota in patients with glioma compared to healthy controls, illustrating genus-level proportions, phylogenetic relationships, and interaction networks. (A) PLS-DA analysis of salivary microbiota characteristics in samples grouped by category (HC, healthy control; HGG, high-grade gliomas). (B) Evolutionary relationships among the salivary microbiota of patients with glioma and healthy controls, as determined by LEfSe analysis. The selected microbial representatives, comprising d: class Gracilibacteria; n: family Eggerthellaceae; h: order Flavobacteriales; k: family Aerococcaceae; p: family Lactobacillaceae; t: family Porphyromonadaceae; and r: family Paludibacteraceae. (C) Proportional distribution of genus-level salivary microbiota within samples. (D) Statistical depiction of salivary microbiota LDA scores featured in the results arising from LEfSe analysis. (E) Interaction networks among different genera, with positive and negative associations represented in red and green, respectively. The size of each node reflects the abundance of the species, and the thickness of the lines indicates the strength of the correlation, with thicker lines denoting higher levels of interconnectivity among species.)

Fig. 1

Fig. 2 Gender-specific salivary microbiome profiles in glioma patients in comparison to healthy individuals. PLS-DA analyses contrasting salivary microbiome samples between healthy females and glioma patients (A), and between healthy males and glioma patients (B). Community bar charts depicting the top 30 microbial species by abundance in the saliva of healthy females in comparison to patients with glioma (C), and in healthy males in comparison with patients with glioma, with species abundance depicted along the y-axis against sample names on the x-axis (D). Circos plots illustrating the distribution of microbial species in the saliva across genders in both healthy individuals and glioma patients, in which the inner ribbon colors represent different groups and their lengths indicate the proportion of species distribution within the samples (E–F).

Fig. 2

3.3 The impact of gender on disparities in the salivary microbiota of patients with glioma, as determined by intergroup significance testing

Next, we investigated how gender influenced disparities in the salivary microbiota between glioma patients and healthy individuals by performing rigorous intergroup significance testing. This analysis aimed to identify bacteria at the genus level that exhibited significant variances, thereby revealing the specific mechanisms by which gender exerts impact on microbial composition. Our findings indicated a significant reduction in the abundance of Leptotrichia, Selenomonas, Abiotrophia, and Lachnoanaerobaculum genera in female patients with glioma when compared to the healthy control group, a phenomenon that may relate to the unique physiological or immune responses specific to female patients with glioma (Fig. 3A and B). Conversely, in male patients with glioma, we detected a reduction in the abundance of Haemophilus, Granulicatella, and Capnocytophaga genera, with a significant increase in Atopobium, potentially reflecting distinct adjustments in the structure of the salivary microbiota in male glioma patients (Fig. 3C and D). Bar chart analysis further revealed differences in LDA scores among species within the salivary microbiomes (LDA threshold = 2.0). These findings indicated a marked predominance of bacterial species such as g_Cryptobacterium in female patients with glioma. In contrast, species such as p_Patescibacteria and c_Fusobacteriia were more prevalent in the salivary microbiomes of healthy females. We also found that c_Coriobacteriia and o_Coriobacteriales played a pivotal role in the salivary microbiomes of male patients with glioma. Furthermore, the presence of o_Enterobacterales and f_Pasteurellaceae was more significant in the salivary microbiomes of healthy males (Supplementary Figure 5). Next, network analysis was used to elucidate the complex correlations between the salivary microbiota of healthy individuals of different genders and those of patients with glioma, thus providing comprehensive insights into microbial interactions in the context of glioma (Supplementary Figure 6). Significant differences in microbial populations could have potential associations with the onset and progression of glioma, thus suggesting that gender plays a crucial role in the variation of salivary microbiome composition in patients with glioma.Fig. 3 Bar graphs depicting gender-specific differences in salivary microbiome abundance between patients with glioma and healthy individuals. These bar graphs compare the relative abundance of salivary microbial genera between female patients with glioma and healthy females, represented as percentages of relative abundance (A). Statistical analysis of microbial differences within the female cohort, including g_Leptotrichia, g_Selenomonas, g_Abiotrophia, and g_Lachnoanaerobaculum (B). Bar chart comparing the relative abundance of salivary microbial genera between male patients with glioma and healthy male controls, expressed as a percentage (C). Statistical analysis of microbial differences illustrating disparities in the relative abundance of salivary microbial genera between male patients with glioma and healthy males, emphasizing the importance of g_Haemophilus, g_Granulicatella, g_Capnocytophaga, and g_Atopobium (D). *P < 0.05 and **P < 0.01.

Fig. 3

3.4 Selection and analysis of key microbial features for the construction and evaluation of a diagnostic model

Next, we performed a systematic feature selection process using two sequential algorithms: LASSO-logistic regression and random forest analysis. First, we used the “glmnet” package in R with 10-fold cross-validation under a binomial family setting to determine model microbial communities and coefficients based on differential microbes. The LASSO-logistic regression algorithm successfully identified eight significant microbial features (Fig. 4A and B). Next, we used the “microeco” package to select characteristic microbial communities using a random forest model; this identified twelve additional key microbial features. The left panel of Fig. 4C shows a MeanDecreaseGini plot, which was used to analyze the contribution of each species to model classification, with higher values indicating better classification performance. The right panel shows a box plot illustrating differences in relative abundance. The features identified by both algorithms were compared, and overlapping features were selected as optimal microbial communities for further analysis. Subsequently, we performed rigorous Venn diagram analysis to extract intersections of the microbial features selected by the two algorithms. This meticulous process ultimately identified five key microbial features that were deemed to play pivotal roles in the construction of a diagnostic model. The formula used to calculate the Risk Score within the diagnostic model was derived as follows: Risk Score = g_Atopobium × 0.61733369 + g_Bacteroides × 2.93668876 + g_Capnocytophaga × (−1.79325740) + g_Leptotrichia × (−0.07011431) + g_Cardiobacterium × (−0.64509608). This formula represents a critical component of the diagnostic model and provides a quantitative assessment of risk based on the presence and abundance of specific microbial features. To rigorously assess the reliability and stability of our diagnostic model, we conducted several key analyses. First, we calculated the area under the curve (AUC) metric, as depicted by Fig. 4D. The AUC is a critical measure used to evaluate the performance of binary classification models, and provides an indicator of a model's ability to distinguish between different classes. In addition, we used principal component analysis (PCA) to provide a more rigorous analysis of data structure (Fig. 4E). PCA was used to visualize the separation of data points in a multidimensional space, thereby validating the effectiveness of the model. Collectively, these analyses highlight the robustness and predictive power of our diagnostic model, indicating its potential utility for evaluating the risk of glioblastoma.Fig. 4 Enhanced model accuracy derived from LASSO regularization and RandomForests feature selection, as demonstrated by advanced AUC and PCA metrics. The coefficient of variation and distribution under LASSO regularization, highlighting the precision of our methodology (A). The optimized selection of lambda.min based on likelihood deviation within the LASSO coefficient distribution (B). The MeanDecreaseGini chart (left), generated by RandomForests feature selection analysis, was used to quantify the influence of various species on the classification accuracy of the model, in which higher MeanDecreaseGini values correlate with enhanced performance. Conversely, the right section of (C) features a boxplot showcasing relative abundance differences, thus offering insights into species impact. The models exceptional AUC (D) and PCA (E) outcomes, confirmed the predictive reliability and dimensional reduction capability of the model, thus emphasizing the efficacy of the method.

Fig. 4

4 Discussion

The advancing field of microbiome research highlights the intricate interplay between microbial communities and human health. Recent investigations have highlighted the potential impact of the oral microbiota on a variety of central nervous system diseases, including Alzheimer's disease [18] and ischemic stroke [19]. In oncological research, dysbiosis of the oral microbiome has been closely associated with the development of oral squamous cell carcinoma [20], head and neck cancer [21], and glioma [10]. Oral microbes can influence the gut microbiota of a host via the oral-gut microbiome axis, thus leading to dysbiosis and the onset of disease [22]. A previous study identified imbalances in the gut microbiome of a mouse model of glioma and patients with glioma when compared to healthy controls, including alterations in β-diversity [23]. Our comprehensive analysis of the saliva microbiome of patients with glioma, focusing particularly on gender-specific differences, is critically important if we are to enhance our understanding of the pathogenesis of glioma and its potential diagnostic and therapeutic implications. In our study, we elucidated the impact of the salivary microbiome on the risk and progression of glioma, highlighting the intricate interplay between microbial dysbiosis and alterations in host-microbiome interactions. Notably, the presence of gender-specific microbial signatures highlighted the potentially pivotal role of biological gender in modulating these interactions. Given the observed disparities in the incidence and prognosis of glioma between genders, we hypothesize that this variation may partly be attributed to differences in the composition of the salivary microbiome between males and females. These findings are in alignment with recent studies indicating that male patients with brain gliomas are more susceptible to this disease and have poorer outcomes when compared to their female counterparts, exhibiting distinct gender disparities [24]. Moreover, the ability of the gut microbiome to regulate hormone levels by altering reabsorption rates, as well as modulate endocrine signaling through the biotransformation of endogenous hormones [25], highlights the significance of this research. Our study accentuates the potential value of gender-dependent modulation by the salivary microbiome on the risk and therapeutic response to brain glioma disease.

Moreover, the findings of this study identified significant differences in the salivary microbiome between male and female patients with gliomas when compared to healthy individuals. Specifically, the predominant differential bacterial genera in male patients with glioma included Haemophilus, Granulicatella, Capnocytophaga, and Atopobium. In female patients with glioma, the genera predominantly consisted of Leptotrichia, Selenomonas, Abiotrophia, and Lachnoanaerobaculum. Notably, Leptotrichia has been reported to be the most abundant bacterial genus in the oral gingival plaques of pregnant women, coinciding with the period of elevated estrogen levels during human pregnancy [26,27]. Our analysis discovered a significant increase in Leptotrichia in the healthy female population when compared to those with gliomas, leading us to hypothesize that a positive correlation between Leptotrichia and estrogen levels might be associated with the incidence of gliomas in females. The elevation in estrogen levels mediated by Leptotrichia may potentially reduce the incidence rate of gliomas in women; however, whether this association is causal remains unclear. Furthermore, we identified a significantly higher abundance of the Atopobium genus in the salivary microbiome of male patients with glioma when compared to healthy controls. Intriguingly, the abundance of Atopobium has been shown to differ significantly in vaginal samples from patients with endometriosis when compared to healthy individuals, with almost no presence in the former population [28], a condition that is characterized by high estrogen dependency and secretion as its endocrine hallmark [29]. This suggests that the antagonistic relationship between Atopobium and estrogen might explain the higher salivary abundance of Atopobium in male patients with glioma when compared to healthy individuals. The observed differences in salivary microbial abundance provide us with valuable insights into the etiology of glioma. Estrogen-related bacterial genera may directly or indirectly modulate the microenvironment of the brain, thus influencing the development and progression of glioma. Future identification of specific microbial characteristics will further elucidate the potential mechanisms by which the oral microbiome may impact glioma biology.

Recent applications of machine learning algorithms, such as PLS-DA and LEfSe, have facilitated detailed analyses of the salivary microbiome, uncovering significant disparities between patients with glioma and healthy individuals. These computational approaches have been instrumental in identifying key salivary microbial markers for glioma, employing algorithms such as LASSO-logistic regression and random forest analyses. Notably, LASSO-logistic regression identified eight characteristic microbes, while random forest revealed twelve. The diagnostic model was subsequently refined to include only the overlapping microbial markers, thus marking the clear advantages of utilizing algorithms for the selection of diagnostic or prognostic factors in recent studies [[30], [31], [32]]. Moreover, the high predictive accuracy of this diagnostic model, validated by PCA, highlights its potential utility in clinical diagnostics, offering a non-invasive, saliva-based tool for the detection of glioma. This enhances our comprehension of the microbial landscape associated with glioma. Our findings also highlight the potential therapeutic implications of targeting the oral microbiome when treating glioma. Given the observed microbial dysbiosis in patients with glioma, strategies aimed at modulating the oral microbiome could represent novel therapeutic avenues. For instance, the application of probiotics to modulate the microbial equilibrium and influence glioma-associated biological activities has already been investigated. One previous study highlighted the capacity of a probiotic blend to suppress intracranial tumor growth in a mouse model of glioma by altering the gut microbiome, facilitated by the restoration of intestinal barrier integrity [33]. Furthermore, our discovery of gender-specific microbial features in the salivary microbiome indicates significant gender differences in glioma, suggesting that personalized microbiome interventions that consider individual patient characteristics (including gender) could enhance therapeutic outcomes. Future research should aim to elucidate the mechanisms underlying observed gender differences in the impact of the salivary microbiome on glioma and investigate the therapeutic potential of modulating the microbiome in the treatment of glioma.

5 Conclusion

Our comprehensive analysis demonstrates the critical role of the salivary microbiome in the pathogenesis of glioma and potential diagnostic applications, particularly highlighting gender-specific differences. Our research successfully identified distinctive microbial taxa in the saliva of male and female patients with glioma, distinguishing them from healthy controls. Notably, genera such as Atopobium and Leptotrichia showed significant variations in abundance linked to gender, thus suggesting their influence on glioma susceptibility and hormone interactions. These findings emphasize the potential of saliva as a non-invasive diagnostic medium that is capable of reflecting underlying disease processes via microbial signatures. The integration of machine learning techniques has further refined the identification of key microbial markers, thus enhancing the diagnostic precision for glioma. These findings emphasize the promise of saliva-based diagnostics in clinical settings, offering insights into disease progression and response to treatment. Future research should investigate the specific mechanisms by which the salivary microbiota can influence the development of glioma and investigate targeted therapeutic interventions that could modulate these microbial communities.

Ethics approval and consent to participate

Not applicable.

Funding

This work was supported by 10.13039/501100004826 Beijing Natural Science Foundation (L232075 ), the 10.13039/501100001809 National Natural Science Foundation of China (Nos. 82304151 and 32100647 ), Supported by the 10.13039/501100012226 Fundamental Research Funds for the Central Universities (3332023033 ), Funding from National Clinical Research Center for Obstetrics and Gynecology (10.13039/501100009399 Peking University Third Hospital ) (No. BYSYSZKF2023014 ), Key Clinical Project of Peking University Third Hospital (No. BYSYZD2023036 ), and the Cancer Hospital of Chinese Academy of Medical Sciences-Shenzhen Hospital Cooperation Fund (CFA202202023).

Data availability statement

The data used to support the findings of this study are available from the corresponding authors on reasonable request.

CRediT authorship contribution statement

Hao Qin: Writing – original draft, Supervision, Methodology, Funding acquisition, Conceptualization. Jie Liu: Writing – review & editing, Methodology, Investigation, Data curation. Yang-Yang Li: Investigation. Ya-Lan Xu: Writing – review & editing. Yi-Fang Yan: Writing – original draft, Funding acquisition, Conceptualization.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A Supplementary data

The following are the supplementary data to this article:Supplementary Fig. 1 Dilution curves provide a direct indicator of the appropriateness of sequencing data volume. A trend towards a plateau in the latter part of a curve signifies a rational level of sequencing data, beyond which only incremental discoveries of new species are anticipated.

Supplementary Fig. 1

Supplementary Fig. 2 Bar charts depicting the differences in microbial composition within saliva samples from healthy controls and patients with glioma, for both females and males. These charts highlight the top 15 microbial species in terms of abundance across various taxonomic levels, including phylum, class, order, family, and species, along with their respective proportions.

Supplementary Fig. 2

Supplementary Fig. 3 Circos diagram illustrating the distribution of microbial species in saliva samples from healthy male or female control groups compared to those of patients with glioma. Microorganisms are detailed across five taxonomic levels: phylum, class, order, family, and species. The length of each segment within the diagram proportionately reflects the distribution of microbial species within the saliva samples.

Supplementary Fig. 3

Supplementary Fig. 4 Bar charts illustrating gender differences in the abundance of the salivary microbiome in healthy individuals. These charts compare the relative abundances of salivary microbes in healthy male and female populations, expressed as percentages (A). Statistical analysis of microbial diversity revealing significant differences in the relative abundance of saliva microbes between healthy males and females, including Eubacterium nodatum group, Solobacterium, and unclassified members of the family Pasteurellaceae (B). *P < 0.05 and **P < 0.01.

Supplementary Fig. 4

Supplementary Fig. 5 Bar graphs depicting LDA scores for differences in microbial species in saliva between healthy individuals of different genders and patients with glioma. An increase in LDA score indicates a greater influence of the abundance of microbial species on the differentiation effect.

Supplementary Fig. 5

Supplementary Fig. 6 Network analysis depicting correlations among microbial species in the salivary microbiome of healthy individuals (both female and male) compared to patients with glioma. The analysis focused on the top 50 species by overall abundance. Node sizes represent the abundance of microbial species; red lines indicate positive correlations, green lines indicate negative correlations, and a higher number of connections per node denotes closer inter-species relationships.

Supplementary Fig. 6

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.heliyon.2024.e37284.
==== Refs
References

1 Wade W.G. Characterisation of the human oral microbiome J. Oral Biosci. 55 2013 143 148 10.1016/j.job.2013.06.001
2 Keijser B.J.F. Zaura E. Huse S.M. van der Vossen J.M.B.M. Schuren F.H.J. Montijn R.C. ten Cate J.M. Crielaard W. Pyrosequencing analysis of the Oral Microflora of healthy adults J. Dent. Res. 87 2008 1016 1020 10.1177/154405910808701104 18946007
3 Sedghi L. DiMassa V. Harrington A. Lynch S.V. Kapila Y.L. The oral microbiome: role of key organisms and complex networks in oral health and disease Periodontol. 2000 87 2021 107 131 10.1111/prd.12393 34463991
4 Zürcher C. Humpel C. Saliva: a challenging human fluid to diagnose brain disorders with a focus on Alzheimer's disease Neural regeneration research 18 2023 2606 2610 10.4103/1673-5374.373675 37449596
5 Kumar P. Gupta S. Das B.C. Saliva as a potential non-invasive liquid biopsy for early and easy diagnosis/prognosis of head and neck cancer Translational oncology 40 2024 101827 10.1016/j.tranon.2023.101827
6 Ostrom Q.T. Cioffi G. Waite K. Kruchko C. Barnholtz-Sloan J.S. CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2014-2018 Neuro Oncol. 23 2021 iii1 iii105 10.1093/neuonc/noab200 34608945
7 Sampson J.H. Gunn M.D. Fecci P.E. Ashley D.M. Brain immunology and immunotherapy in brain tumours Nat. Rev. Cancer 20 2020 12 25 10.1038/s41568-019-0224-7 31806885
8 Eisenbarth D. Wang Y.A. Glioblastoma heterogeneity at single cell resolution Oncogene 42 2023 2155 2165 10.1038/s41388-023-02738-y 37277603
9 Müller Bark J. Trevisan França de Lima L. Zhang X. Broszczak D. Leo P.J. Jeffree R.L. Chua B. Day B.W. Punyadeera C. Proteome profiling of salivary small extracellular vesicles in glioblastoma patients Cancer 129 2023 2836 2847 10.1002/cncr.34888 37254878
10 Wen Y. Feng L. Wang H. Zhou H. Li Q. Zhang W. Wang M. Li Y. Luan X. Jiang Z. Chen L. Zhou J. Association between oral microbiota and human brain glioma grade: a case-control study Front. Microbiol. 12 2021 746568 10.3389/fmicb.2021.746568
11 Carrano A. Juarez J.J. Incontri D. Ibarra A. Guerrero Cazares H. Sex-specific differences in glioblastoma Cells 10 2021 10.3390/cells10071783
12 Chen S. Zhou Y. Chen Y. Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor Bioinformatics 34 2018 i884 i890 10.1093/bioinformatics/bty560 30423086
13 Bolyen E. Rideout J.R. Dillon M.R. Bokulich N.A. Abnet C.C. Al-Ghalith G.A. Alexander H. Alm E.J. Arumugam M. Asnicar F. Bai Y. Bisanz J.E. Bittinger K. Brejnrod A. Brislawn C.J. Brown C.T. Callahan B.J. Caraballo-Rodríguez A.M. Chase J. Cope E.K. Da Silva R. Diener C. Dorrestein P.C. Douglas G.M. Durall D.M. Duvallet C. Edwardson C.F. Ernst M. Estaki M. Fouquier J. Gauglitz J.M. Gibbons S.M. Gibson D.L. Gonzalez A. Gorlick K. Guo J. Hillmann B. Holmes S. Holste H. Huttenhower C. Huttley G.A. Janssen S. Jarmusch A.K. Jiang L. Kaehler B.D. Kang K.B. Keefe C.R. Keim P. Kelley S.T. Knights D. Koester I. Kosciolek T. Kreps J. Langille M.G.I. Lee J. Ley R. Liu Y.X. Loftfield E. Lozupone C. Maher M. Marotz C. Martin B.D. McDonald D. McIver L.J. Melnik A.V. Metcalf J.L. Morgan S.C. Morton J.T. Naimey A.T. Navas-Molina J.A. Nothias L.F. Orchanian S.B. Pearson T. Peoples S.L. Petras D. Preuss M.L. Pruesse E. Rasmussen L.B. Rivers A. Robeson M.S. 2nd Rosenthal P. Segata N. Shaffer M. Shiffer A. Sinha R. Song S.J. Spear J.R. Swafford A.D. Thompson L.R. Torres P.J. Trinh P. Tripathi A. Turnbaugh P.J. Ul-Hasan S. van der Hooft J.J.J. Vargas F. Vázquez-Baeza Y. Vogtmann E. von Hippel M. Walters W. Wan Y. Wang M. Warren J. Weber K.C. Williamson C.H.D. Willis A.D. Xu Z.Z. Zaneveld J.R. Zhang Y. Zhu Q. Knight R. Caporaso J.G. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2 Nat. Biotechnol. 37 2019 852 857 10.1038/s41587-019-0209-9 31341288
14 McDonald D. Jiang Y. Balaban M. Cantrell K. Zhu Q. Gonzalez A. Morton J.T. Nicolaou G. Parks D.H. Karst S.M. Albertsen M. Hugenholtz P. DeSantis T. Song S.J. Bartko A. Havulinna A.S. Jousilahti P. Cheng S. Inouye M. Niiranen T. Jain M. Salomaa V. Lahti L. Mirarab S. Knight R. Greengenes2 unifies microbial data in a single reference tree Nat. Biotechnol. 2023 10.1038/s41587-023-01845-1
15 Liu C. Cui Y. Li X. Yao M. microeco: an R package for data mining in microbial community ecology FEMS Microbiol. Ecol. 97 2021 10.1093/femsec/fiaa255
16 Jiang B. Sun P. Luo B. GLMNet: graph learning-matching convolutional networks for feature matching Pattern Recogn. 121 2022 108167 10.1016/j.patcog.2021.108167
17 Segata N. Izard J. Waldron L. Gevers D. Miropolsky L. Garrett W.S. Huttenhower C. Metagenomic biomarker discovery and explanation Genome Biol. 12 2011 R60 10.1186/gb-2011-12-6-r60 21702898
18 Wan J. Fan H. Oral microbiome and Alzheimer's disease Microorganisms 11 2023 2550 10.3390/microorganisms11102550 37894208
19 Sun W. Huang S. Yang X. Luo Y. Liu L. Wu D. The oral microbiome of patients with ischemic stroke predicts their severity and prognosis Front. Immunol. 14 2023 1171898 10.3389/fimmu.2023.1171898
20 Zhang L. Liu Y. Zheng H.J. Zhang C.P. The oral microbiota may have influence on oral cancer Front. Cell. Infect. Microbiol. 9 2019 476 10.3389/fcimb.2019.00476 32010645
21 Wu Z. Han Y. Wan Y. Hua X. Chill S.S. Teshome K. Zhou W. Liu J. Wu D. Hutchinson A. Jones K. Dagnall C.L. Hicks B.D. Liao L. Hallen-Adams H. Shi J. Abnet C.C. Sinha R. Chaturvedi A. Vogtmann E. Oral microbiome and risk of incident head and neck cancer: a nested case-control study Oral Oncol. 137 2023 106305 10.1016/j.oraloncology.2022.106305
22 Lu Y. Li Z. Peng X. Regulatory effects of oral microbe on intestinal microbiota and the illness Front. Cell. Infect. Microbiol. 13 2023 1093967 10.3389/fcimb.2023.1093967
23 Patrizz A. Dono A. Zorofchian S. Hines G. Takayasu T. Husein N. Otani Y. Arevalo O. Choi H.A. Savarraj J. Tandon N. Ganesh B.P. Kaur B. McCullough L.D. Ballester L.Y. Esquenazi Y. Glioma and temozolomide induced alterations in gut microbiome Sci. Rep. 10 2020 21002 10.1038/s41598-020-77919-w
24 Tournant F. Immunological sex differences in glioblastoma Nat. Can. (Ott.) 4 2023 10.1038/s43018-023-00660-6 1636-1636
25 Cox L.M. Abou-El-Hassan H. Maghzi A.H. Vincentini J. Weiner H.L. The sex-specific interaction of the microbiome in neurodegenerative diseases Brain Res. 1724 2019 146385 10.1016/j.brainres.2019.146385
26 Condon J.C. Kyathanahalli C. Anamthathmakula P. Jeyasuria P. Estrogen/estrogen receptor action and the pregnant myometrium Current Opinion in Physiology 13 2020 135 140 10.1016/j.cophys.2019.10.017
27 Zhang Y. Wu Z. Li L. Wang X. Fan W. Zhao J. Characterizing the supragingival microbiome of healthy pregnant women Front. Cell. Infect. Microbiol. 12 2022 10.3389/fcimb.2022.1016523
28 Ata B. Yildiz S. Turkgeldi E. Brocal V.P. Dinleyici E.C. Moya A. Urman B. The endobiota study: comparison of vaginal, cervical and gut microbiota between women with stage 3/4 endometriosis and healthy controls Sci. Rep. 2019 2204 10.1038/s41598-019-39700-6 30778155
29 Chantalat E. Valera M.C. Vaysse C. Noirrit E. Rusidze M. Weyl A. Vergriete K. Buscail E. Lluel P. Fontaine C. Arnal J.F. Lenfant F. Estrogen receptors and endometriosis Int. J. Mol. Sci. 21 2020 10.3390/ijms21082815
30 Wu Z. Chen H. Ke S. Mo L. Qiu M. Zhu G. Zhu W. Liu L. Identifying potential biomarkers of idiopathic pulmonary fibrosis through machine learning analysis Sci. Rep. 13 2023 16559 10.1038/s41598-023-43834-z
31 Li Z. Qin Y. Liu X. Chen J. Tang A. Yan S. Zhang G. Identification of predictors for neurological outcome after cardiac arrest in peripheral blood mononuclear cells through integrated bioinformatics analysis and machine learning Funct. Integr. Genom. 23 2023 83 10.1007/s10142-023-01016-0
32 Lin D. Zhao W. Yang J. Wang H. Zhang H. Integrative analysis of biomarkers and mechanisms in adamantinomatous craniopharyngioma Front. Genet. 13 2022 10.3389/fgene.2022.830793
33 Wang L. Li S. Fan H. Han M. Xie J. Du J. Peng F. Bifidobacterium lactis combined with Lactobacillus plantarum inhibit glioma growth in mice through modulating PI3K/AKT pathway and gut microbiota Front. Microbiol. 13 2022 986837 10.3389/fmicb.2022.986837
