Skip to main content

Changes in the nasopharyngeal and oropharyngeal microbiota in pediatric obstructive sleep apnea before and after surgery: a prospective study

Abstract

Objective

To explore the changes and potential mechanisms of microbiome in different parts of the upper airway in the development of pediatric OSA and observe the impact of surgical intervention on oral microbiome for pediatric OSA.

Methods

Before adeno-tonsillectomy, we collected throat swab samples from different parts of the oropharynx and nasopharynx of 30 OSA patients and 10 non-OSA patients and collected throat swab samples from the oropharynx of the above patients one month after the adeno-tonsillectomy. The 16 S rRNA V3–V4 region was sequenced to identify the microbial communities. The correlation analysis was conducted based on clinical characteristics.

Results

There was a significant difference of alpha diversity in different parts of the upper airway of pediatric OSA, but this difference was not found in children with non-OSA. Beta diversity was significantly different between non-OSA and pediatric OSA. At the genus level, the composition of flora in different parts is different between non-OSA and pediatric OSA. The correlation analysis revealed that the relative abundance of Neisseria was significantly correlated with obstructive apnea hypopnea index. Furthermore, the functional prediction revealed that pathways related to cell proliferation and material metabolism were significantly different between non-OSA and pediatric OSA. Besides, the adeno-tonsillectomy has minimal impact on oral microbiota composition in short term.

Conclusion

The changes in upper airway microbiome are highly associated with pediatric OSA. The relative abundance of some bacteria was significantly different between OSA and non-OSA. These bacteria have the potential to become new diagnostic and early warning biomarkers.

Peer Review reports

Introduction

Pediatric obstructive sleep apnea (OSA) is characterized by the recurrent incidence of either partial or complete obstruction of the upper airway during sleep, which disturbs the normal sleep quality and sleep architecture, leads to long-term sleep hypoxia and hypercapnia in children, and causes a series of pathophysiological changes, including growth retardation, behavioral cognitive impairment, memory and intellectual decline in children, mainly manifested as snoring, suffocation, mouth opening breathing with the high incidence rate among children(1.2%-5.7%) [1, 2]. To some extent, OSA is more harmful to children than adults. Hence, the timely detection and intervention of pediatric OSA holds substantial importance in improving the prognosis.

Different from OSA in adult, the main cause of upper airway obstruction in children is adenoid and/or tonsillar hypertrophy [3]. The adenoids and tonsils, as integral components of the immune system, exert a significant influence on the body’s immune response to infectious pathogens within the upper respiratory tract [4]. There are many folds and recesses on the surface of adenoids and tonsils, which are occupied by symbiotic microorganisms shortly after birth [5]. These microorganisms subsequently affect the development of innate mucosal immune response, which may be an important factor in protecting or inducing individuals to produce mucosal immune response during infection [5]. It is reported that most patients with chronic tonsillitis have a layer of “bacterial biofilm” attached to the surface of the tonsils [6]. This biofilm is formed by bacteria adhering to the surface of the body to produce various polysaccharide matrices, proteins, and other complex substances that aggregate with each other, causing various bacteria to adhere to form microbiota colonies, resulting in membrane like substances.

Prior research has established a significant association between the oral microbiota and the development of various systemic ailments, including cardiovascular disease, rheumatoid arthritis, and complications arising from HIV [7]. Presently, the investigation of the upper respiratory microbiota in pediatric patients with OSA primarily centers on the microbiota in the vicinity of the tonsils and adenoids. Certain scholars have posited the “pathogen reservoir hypothesis,” which proposes that the adenoids serve as the origin of pathogens in the tonsils and middle ear; however, additional substantiation of this hypothesis is required [8, 9]. Furthermore, the “bacterial interference theory” has been proposed by certain individuals, which posits that oral symbiotic bacteria impede the growth of pathogens by engaging in resource competition or generating antagonistic nutrients [10]. Recent investigation also demonstrated a significant correlation between adenoid hypertrophy and alterations in the microbiota composition of the upper airways compared to healthy subjects [11]. As a result, the alterations in upper airway microbiota may be closely related to the presence of pediatric OSA, and in-depth research on the correlation between upper airway microbiota and pediatric OSA has important reference value.

In this study, we collected swab samples from different parts of the participants’ nasopharynx and oropharynx. In addition, we also collected oral swab samples from participants after adeno-tonsillectomy. 16 S rRNA pyrosequencing was used to investigate the diversity of upper airway microbiome in pediatric OSA. The objective of this study is (1) to examine the microbial community characteristics of various regions within the upper airway in pediatric patients with OSA; (2) to study whether there are differences in the composition and diversity of microbial communities between pediatric patients with OSA and pediatric patients with non-OSA; (3) to observe the impact of adeno-tonsillectomy on the oral microbiome of pediatric patients with OSA.

Methods

Study participant

This study was designed and performed based on the Helsinki Declaration and the prospective specimen collection [12]. Ethical clearance was approved by the Ethical Review Board of the Zhongnan Hospital of Wuhan University (No.2,021,056). The written consent of every participant was obtained prior to the commencement of the study. All samples were collected exclusively from the Zhongnan Hospital of Wuhan University. Each participant exhibited symptoms of snoring and underwent standard polysomnography. All participants were diagnosed with adenoid and tonsillar hypertrophy through nasopharyngeal endoscopy. Severity of sleep apnea was categorized by Society of Otorhinolaryngology Head and Neck Surgery in Chinese Medical Association criteria [13]: non-OSA (obstructive apnea hypopnea index [OAHI], < 1 events/h), mild OSA (OAHI, 1–5 events/h), moderate OSA (OAHI, 6–10 events/h), and severe OSA (OAHI > 10 events/h). Participants were excluded if they met any of the following criteria: (1) There are infected lesions in the oral cavity; (2) People or their family members with acute respiratory infectious diseases; (3)People with nose, pharynx, throat diseases and chronic respiratory diseases that affect sleep; (4)Patients with dental caries and periodontal disease; (5)Abnormal preoperative biochemical examination results. The participants were categorized into four groups, namely the non-OSA group, mild OSA group, moderate OSA group, and severe OSA group, based on the results obtained from polysomnography.

Sample collection and DNA extraction

All subjects were forbidden to brush their teeth, rinse their mouth and eat the preceding night before collection. The trained professionals collected samples during adenotonsillectomy. Sterile throat swab shall be used to wipe the mucosa of the patient’s tongue base (oropharynx site)(Part A), surface of palatine tonsil (Part B), palatine tonsillar capsule after tonsillectomy (Part C) and adenoid (nasopharynx site) (Part D). The detailed collection location is shown in supplementary Figure S1. At the same time, it shall be noted that the throat swab shall not contact the mucosa of both sides of the mouth and cheek, hard palate teeth and other parts. Similarly, during a follow-up one month after adeno-tonsillectomy, swab samples from the participants’ tongue base were collected. After collection, put the sample into the liquid nitrogen tank for quick freezing and quickly transfer it to the laboratory (Fig. 1).

Fig. 1
figure 1

Overview of study design. Data analysis of the microbiome was performed to explore their correlation and relevance to clinical parameters

After placing swab samples in TENS buffer (5 M Sodium Chloride, 10% SDS, Triton X-1000, Tris-HCl, EDTA), 10% SDS and 20 mg/mL proteinase K, the samples were incubated overnight at 55°C. Proteins were eliminated through phenol/chloroform/isoamyl alcohol extractions, and DNA was subsequently precipitated using isopropanol on the following day. After being washed in 75% ethanol, the DNA was then resuspended in TE buffer. The quantification of DNA was conducted using the Qubit Fluorometer, employing the Qubit dsDNA BR Assay kit (Invitrogen, USA), while the quality assessment involved running an aliquot on a 1% agarose gel. The amplification of the variable regions V3–V4 of the bacterial 16S rRNA gene was achieved using degenerate PCR primers, namely 341F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’). Both the forward and reverse primers were modified with Illumina adapter, pad, and linker sequences. PCR enrichment was conducted in a 50 µL reaction mixture containing 30ng of template, fusion PCR primer, and PCR master mix. The PCR cycling conditions consisted of an initial denaturation step at 94 °C for 3 min, followed by 30 cycles of denaturation at 94 °C for 30 s, annealing at 56 °C for 45 s, extension at 72 °C for 45 s, and a final extension step at 72 °C for 10 min. The resulting PCR products were purified using AmpureXP beads and eluted in Elution buffer. The quality of the libraries was assessed using the Agilent 2100 bioanalyzer (Agilent, USA). The sequencing process on the Illumina MiSeq platform (BGI, Shenzhen, China) utilized validated libraries and adhered to the standard pipelines of Illumina, resulting in the production of 2 × 300 bp paired-end reads.

Bacterial 16  rDNA sequencing and bioinformatic analysis

The raw reads underwent filtration to eliminate adaptors, low-quality, and ambiguous bases. Subsequently, the paired-end reads were merged with tags using the Fast Length Adjustment of Short reads program (FLASH, v1.2.11) [14] to obtain the tags. These tags were then clustered into operational taxonomic units (OTUs) using the UPARSE software (v7.0.1090) [15] with a 97% cutoff value. Additionally, chimera sequences were identified by comparing them with the Gold database using UCHIME (v4.2.40) [16]. Subsequently, the taxonomic classification of OTU representative sequences was conducted using the Ribosomal Database Project (RDP) Classifier v.2.2, employing a minimum confidence threshold of 0.6. The classifier was trained on the Greengenes database v201305 by QIIME v1.8.0 [17]. To obtain the OTU abundance statistics table for each sample, the USEARCH_global [18] was employed to compare all Tags back to OTU. At the OTU level, MOTHUR (v1.31.2) [19] and QIIME (v1.8.0) [17] were utilized to estimate alpha and beta diversity, respectively. Principal Coordinate Analysis (PCoA) was conducted using QIIME (v1.8.0) [17]. Sample clustering was performed by QIIME (v1.8.0) [15] using the UPGMA method. LEfSe clustering and LDA analysis were carried out using LEfSe (https://huttenhower.sph.harvard.edu/galaxy/). COG and KEGG functions were predicted using the PICRUSt software [20]. A barplot representing different classification levels was generated using GraphPad Prism (version 9.4.1). The raw data of this study were deposited in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with the accession numbers PRJNA996776.

Statistical analysis

The Mann-Whitney test was employed to assess the presence of a significant difference between the two groups in relation to continuous variables, while the Kruskal-Wallis test was utilized to evaluate the existence of a significant difference among multiple groups. In order to examine any disparities between groups in terms of frequency distribution, the chi-square test was employed. The dependent variables in this study were the clinical variables, while the independent variables consisted of the microbiota characteristics associated with OSA. Multiple linear regression analysis was conducted to examine the internal correlation between these variables. The statistical analysis was performed using GraphPad Prism (version 9.4.1), and p < 0.05 was considered statistically significant.

Results

Characteristics of the participants

This study included a cohort of 40 participants, comprising 10 individuals diagnosed with non-OSA, 10 individuals with mild OSA, 10 individuals with moderate OSA, and 10 individuals with severe OSA. The average age of the patients was 7.1 ± 2.9 years old, including 10 female and 30 male. Except for OAHI, micro-arousal index (MAI) and lowest blood oxygen saturation (LSaO2), there were no statistically significant differences observed among the four groups with regards to age, sex, BMI, mean blood oxygen saturation (MBOS), and white blood cell count (WBCC) (P > 0.05). The detailed demographic and clinical characteristics of the four groups are presented in Table 1. Following the screening process, 160 samples from 40 patients were subjected to 16 S rRNA MiSeq sequencing in accordance with established standards.

Table 1 Demographics and Clinical Characteristics

Microbiome diversity analysis

The differences in microbiome between children without OSA and children with various severity of OSA were evaluated. Analysis of alpha diversity in OSA groups showed that part B had significantly decreased shannon index and significantly increased simpson index compared with other parts (Fig. 2), while part C had significantly increased chao index and ace index compared with other parts (Fig. S2). However, the statistically significant difference of alpha diversity index among four parts in non-OSA group was different from OSA groups (Figure S2). Among non-OSA group and various severity of OSA groups, alpha diversity had no statistically significant difference in same parts of upper airway(P > 0.05) (Figure S3).

In order to visually represent the variation in microbiome composition across samples, beta diversity was assessed using Principal Coordinate Analysis (PCoA). PCoA of unweighted UniFrac (Fig. 3) showed that the samples of non-OSA group and various severity of OSA group were significantly separated, indicating that the beta diversity of overall upper airway microbial composition was different between non-OSA group and various severity of OSA group. On the contrary, there is no significant difference in beta diversity among various severity of OSA group (Figure S4).

Fig. 2
figure 2

The microbial alpha diversity comparison among four parts in OSA group (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001)

Fig. 3
figure 3

Comparisons of beta diversity by PCoA between non-OSA group and various severity of OSA group. (A) Mild OSA vs. Non-OSA (R2 = 0.026, P < 0.01); (B) Moderate OSA vs. Non-OSA (R2 = 0.035, P < 0.01); (C) Severe OSA vs. Non-OSA (R2 = 0.024, P < 0.05)

Composition and comparison of the Microbiome

Based on the annotation of OTUs, we calculated the composition of microbiomes and plotted the microbial relative abundance on different levels in four parts of non-OSA group and various severity of OSA group. The composition of upper airway microbiome was different among different parts. At the genus level in part A, the proportion of Prevotella and Veillonella was up to 44-57% in non-OSA group, mild OSA group and moderate OSA group (Fig. 4), while in severe OSA group, the proportion of Prevotella, Neisseria and Veillonella was the highest (Fig. 4). At the genus level in part B, the proportion of Prevotella and Fusobacterium was up to 45-63% in non-OSA group and various severity of OSA group (Fig. 4). At the genus level in part C, the proportion of Prevotella, Haemophilus and Fusobacterium were the highest in non-OSA group and mild OSA group (Fig. 4), and the proportion of Prevotella, Haemophilus and Streptococcus were the highest in moderate OSA group and severe OSA group (Fig. 4). At the genus level in part D, the proportion of Prevotella and Haemophilus were the highest in non-OSA group, and the proportion of Prevotella, Moraxella and Haemophilus were the highest in mild OSA group, and the proportion of Prevotella and Fusobacterium were the highest in moderate OSA group, and the proportion of Prevotella and Others were the highest in severe OSA group (Fig. 4).

At the genus level, we selected the 10 species with the highest abundance to compare the relative abundance differences of each species in each group. In non-OSA group, the relative abundance of Fusobacterium, Veillonella and Treponema were significantly different among the four parts (Fig. 5A). In mild OSA group, the relative abundance of Fusobacterium, Veillonella, Treponema and Neisseria were significantly different among the four parts (Fig. 5B). In moderate OSA group, the relative abundance of Veillonella was significantly different among the four parts (Fig. 5C). In severe OSA group, the relative abundance of Veillonella and Corynebacterium were significantly different among the four parts (Fig. 5D). Overall, the relative abundance of Haemophilus, Neisseria and Alloprevotella were significantly different between OSA group and non-OSA group (Fig. 5E), and the relative abundance of Haemophilus and Neisseria were significantly different between non-OSA group and various severity of OSA group (Fig. 5F). In part A and part D, there was no significant difference in the relative abundance of bacterium between non-OSA group and various severity of OSA group (Figure S5). In part B, the relative abundance of Haemophilus, Treponema and Neisseria were significantly different between non-OSA group and various severity of OSA group (Figure S5). In part C, the relative abundance of Prevotella was significantly different between non-OSA group and various severity of OSA group (Figure S5).

Fig. 4
figure 4

The microbial relative abundance at the genus level in four parts. (A) Non-OSA; (B) Mild OSA; (C) Moderate OSA; (D) Severe OSA

Fig. 5
figure 5

Comparison of key species differences.(A) Comparison of four parts in non-OSA group; (B) Comparison of four parts in mild OSA group; (C) Comparison of four parts in moderate OSA group; (D) Comparison of four parts in severe OSA group; (E) OSA vs. Non-OSA; (F) Non-OSA vs. mild OSA vs. moderate OSA vs. severe OSA (*P < 0.05, **P < 0.01, ***P < 0.001)

Phylogenetic characteristics of microbial communities

Using Lefse analysis and linear discriminant analysis (LDA) Score, we conducted a characterization of the upper airway microbiota in both the non-OSA group and the various severity levels of the OSA group. The phylogenetic distribution of upper airway microbiomes associated with non-OSA and the different severity levels of OSA is visually represented in the cladogram (Figure S6). The histograms depicting the LDA scores were computed, revealing significant difference in bacterial composition between the non-OSA group and the different severity levels of the OSA group (Figure S7). In Figure S7A, 22 kinds of biomarkers were significantly enriched in non-OSA group (P < 0.05; LDA score > 2), and 25 kinds of biomarkers were significantly enriched in mild OSA group (P < 0.05; LDA score > 2). In Figure S7B, 13 kinds of biomarkers were significantly enriched in non-OSA group (P < 0.05; LDA score > 2), and 16 kinds of biomarkers were significantly enriched in moderate OSA group (P < 0.05; LDA score > 2). In Figure S7C, 5 kinds of biomarkers were significantly enriched in non-OSA group (P < 0.05; LDA score > 2), and 30 kinds of biomarkers were significantly enriched in severe OSA group (P < 0.05; LDA score > 2).

Multivariate analysis for microbiota signatures

The clinical variables (Age, BMI, OAHI, MAI, MBOS and WBCC) were taken as dependent variables respectively, and the relative abundance of the different microflora was taken as independent variables. We mainly analyzed the correlation between dependent variable and independent variable in the tonsil area (Tonsil recess and Tonsillar capsule after tonsillectomy). As shown in Table S1, the relative abundance of various microflora was significantly correlated with clinical characteristics. However, only the relative abundance of Neisseria was significantly correlated with OAHI (P < 0.05).

Function analysis

To investigate the potential mechanisms underlying the upper airway microbiome in pediatric OSA, we conducted an analysis utilizing Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG). There were 25 functional pathwaysIn in COG analysis (Figure S8A) and 28 functional pathways in KEGG analysis (Figure S8B). What’s more, we compared the functional pathways with significant difference between non-OSA group and various severity of OSA group. In COG analysis, differential functional pathways between non-OSA group and various severity of OSA group were mainly related to cellular processes and information storage. In KEGG analysis, differential functional pathways between non-OSA group and various severity of OSA group were mainly related to metabolism. The detailed information of functional pathways with significant differences were showed in Table S2.

Postoperative oral microbiota characteristics

A total of 38 patients completed the collection of oral microbiota samples in tongue base (oropharynx site) after surgery. One month after the patient underwent adeno-tonsillectomy, there was no difference in alpha diversity of microbiota in the tongue base of each group (P>0.05), but the difference in beta diversity of microbiota still existed among the groups (P<0.001) (Figure S9). Comparison of preoperative and postoperative results, there was no difference in alpha diversity of oral microbiota (P>0.05) (Figure S10), while there was a difference in beta diversity of oral microbiota (P<0.001) (Fig. 6A). To evaluate the differences more intuitively between groups, we next conducted an analysis of genus abundance of predominant bacterial. The differential microbiota among the postoperative groups were Fusobacterium, Veillonella, Streptococcus, Neisseria, and Treponema (Fig. 6B). In addition, there was no significant difference in the relative abundance of predominant bacteria between the preoperative and postoperative groups in the non OSA group, mild OSA group and moderate OSA group (Fig. 6C-E). Only in the severe OSA group, there was a significant difference in the relative abundance of Moraxella, Treponema, and Porphyromonas in the predominant oral microbiota between preoperative and postoperative severe OSA patients (Fig. 6F).

Fig. 6
figure 6

Postoperative oral microbiota characteristics. (A) Comparisons of beta diversity by PCoA among preoperative and postoperative groups; (B) Analysis of genus abundance of predominant bacteria among postoperative groups; (C-F) Analysis of genus abundance of predominant bacteria between preoperative and postoperative groups; C) non-OSA; (D) Mild OSA; (E) Moderate OSA; (F) Severe OSA

Discussion

Pediatric OSA is characterized by a range of pathophysiological alterations resulting from the recurrent incidence of partial or complete obstruction in the upper airway during sleep, thereby disrupting the respiratory function and sleep patterns in children [2]. OSA is the most serious disease in pediatric sleep disordered breathing. Due to its elevated prevalence and consequential long-term complications, OSA has garnered increasing attention from both parents and society [1, 2]. The oral cavity is an important part of the upper respiratory tract of the human body. It stores a large number of biological information substances, such as human oral microorganisms, DNA, RNA, proteins, metabolites, etc., which can reflect the physiological and pathological conditions of the body locally and generally [21, 22]. Compared with blood and samples from other parts of the body, the collection of oral swabs has the advantages of non-invasive, convenient collection, convenient storage and transportation, and is conducive to large-scale screening of certain diseases [21]. The research on oral biomarkers has attracted more and more attention, and has been widely used to study the pathogenesis, drug monitoring and efficacy evaluation of various diseases, such as oral diseases (caries, periodontal disease, oral squamous cell carcinoma, etc.) [23, 24] and systemic diseases (autoimmune diseases, cardiovascular diseases, HIV, etc.) [25,26,27]. It is the complex interaction between the microbial flora and the host microenvironment that maintains the dynamic balance of the microecology of the upper airways. Therefore, studying the changes in the microbiota of the upper airwaysis is crucial in gaining a better comprehension of the pathogenesis of pediatric OSA.

Based on our research findings, a notable variation in alpha diversity was observed across distinct regions of the upper airway among children diagnosed with OSA, especially around the tonsils, but this difference was not found in children with non-OSA. In addition, the PCoA of upper airway microbiome can significantly distinguish patients with non-OSA and patients with OSA of different severity. The aforementioned findings suggested that there may be alterations in the upper airway microbiome of children diagnosed with OSA. In OSA patients, the fluctuation in upper airway pressure would significantly reduce airflow, which affects the moisture and oxygen concentration in the upper respiratory tract and could probably lead to disorders of the upper airway microbiome in patients with OSA [28].

In the previous study, Jervis et al. [29] analyzed the middle ear effusion, nasopharyngeal swab and adenoid biopsy samples of 11 children with secretory otitis media in Australia, and found that there were significant differences between the middle ear effusion, nasopharyngeal swab and adenoid flora, and Haemophilus, Streptococcus and Moraxella existed in all samples. Ren et al. [10] analyzed the microflora in adenoids of 67 children who had undergone surgical removal of adenoids, and found that the microflora of adenoids were very different from that of other parts of the human body. The microflora of adenoids was also different among different individuals, but the core flora was similar, mainly including four phyla: Firmicutes, Proteobacteria, Fusobacteria and Actinobacteria. Similarly, our study also found that the microbial composition of different parts was different. Prevotella was highly enriched in each part, while Corynebacterium was highly enriched only in the region of nasopharynx site of severe OSA. In addition, the relative abundance of Prevotella, Fusobacterium, Veillonella, Treponema, Neisseria and Alloprevotella between OSA group and non-OSA group was significantly different, which indicated that the flora of upper airway is diverse, and the oral-colonized bacteria in different parts are quite different. Pediatric OSA probably is the result of multiple symbiotic or pathogenic bacteria and their immune responses, rather than a single bacterium on the surface of tissues [30]. These symbiotic bacteria exist in the oral tissues of OSA patients for a long time, which may be the reason for the previous failure of antibiotic treatment for chronic tonsillitis. It also shows that our understanding of the pathogenesis of OSA in children is not deep enough.

Currently, there is no study on the change of upper airway microbiome with the development of OSA and the upper airway microbiome reflecting the severity of symptoms is also unclear. A previous large-scale study on nasal microbiome of adult OSA demonstrated that there was a significant and positive correlation observed between the abundance of Streptococcus in the nasal mucosa and the apnea-hypnea index in adult individuals with OSA [31]. The study by Marazzato et al. suggested that adenoid hypertrophy is correlated with alterations in the nasal microbiota composition when compared to individuals in a healthy state [11]. Given that there may be involvement of multiple bacteria in the occurrence and progression of pediatric OSA, we conducted a correlation analysis with the relative abundance of significantly different bacteria between groups as the independent variable and clinical characteristics as the dependent variable. It was found that in the tonsil area (Tonsil recess and tonsillar capsule), only the relative abundance of Neisseria had a significant correlation with OAHI. With the increase of OAHI, the relative abundance of Neisseria also increased.

Neisseria is a group of gram-negative diplococci with oxidase positive. Neisseria gonorrhoeae and Neisseria meningitidis in Neisseria are common pathogenic Neisseria causing human diseases, while other strains belong to normal flora or conditional pathogens parasitizing human upper airways [32]. According to report, certain strains of Neisseria that are not gonococci or meningococcal have the potential to cause primary or secondary respiratory tract infections, thereby posing a significant challenge to accurate clinical and laboratory diagnosis and treatment [33]. Additionally, it is noteworthy that symbiotic Neisseria exhibits immunosuppressive properties that can enhance the expression of host immunosuppressive molecules and impede the release of immune inducible factors into the host environment [34, 35]. The research conducted by Zhu et al. demonstrates that symbiotic Neisseria can prompt dendritic cells to generate immunosuppressive IL-10, which consequently restrains the proliferation of T cells upon antigenic stimulation [36]. It is reported that the bacteria that cause tonsillitis not only reside on the surface of the tonsils, but also in the deep tissues of the tonsils. Therefore, the culture results on the surface of the tonsils may not always indicate the true pathogen, leading to sometimes unsuccessful antibiotic treatment in clinical practice [37]. Consequently, it is reasonable to postulate that the proliferation of Neisseria within the tonsil area of pediatric patients with OSA may be intricately associated with the hypertrophy of tonsils. Nevertheless, the causative association between these factors remains uncertain and warrants further investigation for validation.

In KEGG analysis, we found that, compared with non OSA group, “Phosphonate and phosphinate metabolism”, “Tyrosine metabolism”, “Glutathione metabolism” and “Tryptophan metabolism” significantly increased in severe OSA group. This result is different from that of zhang et al., which indicated that the function of adenoid flora in pediatric OSA patients was characterized by decreased amino acid metabolism [38]. Previous study demonstrated that Nitrotyrosine levels were positively correlated with hypoxemia severity in OSA, and intracellular Glutathione level significantly increased for OSA patients [39, 40]. What’s more, some people hypothesized the increase in Tryptophan metabolism probably is important in OSA-related cardiovascular diseases and might be one of the players in the relationship between cancers and OSA [41]. All the above results indicate that the research on upper airway microbiome should not only focus on the change of microbial composition, but also explore the metabolic pathway of upper airway microecological function genes to enable people to better understand the pathogenesis of OSA.

To facilitate a more comprehensive analysis of the impact of surgical interventions on the oral microbiome of pediatric patients with OSA, oral swab samples were collected one month post adeno-tonsillectomy. Prior research has identified substantial alterations in the diversity of gut microbiome among patients following interventions such as surgical procedures or medication administration [42, 43]. What is more, Lin et al. reported that oral microbial diversity were significantly decreased in postoperative thyroid cancer patients [44]. However, in our analysis results, unexpectedly, there was no significant alterations between preoperative and postoperative groups in alpha diversity and relative abundance of predominant oral microflora, except for those with severe OSA. These findings indicate that while surgery may enhance the ventilatory function, it may not significantly affect the homeostasis of the upper airway microenvironment in short term. A prior randomized controlled trial demonstrated that the utilization of antibiotics at varying follow-up intervals resulted in significant alterations to the intestinal microbiota of the subjects, while the composition of oral microbiome remained relatively stable [45]. This suggests that oral microbiome may possess robust tolerance and resistance to numerous unfavorable conditions, which may partially elucidate the limited fluctuations in oral microbiome after the operation [46]. Additionally, it may also be caused by the transfer of bacterial communities in other parts of the upper airway. Therefore, the enduring influence of oral microbiota on pediatric patients with OSA may potentially contribute to the recurrence of adenoidectomy in certain children. Nonetheless, our findings can only serve as a point of reference, given the dearth of pertinent research on the impact of adeno-tonsillectomy on the microbiome of the upper airways.

Although we try to avoid the design defects of previous studies, our research still has some limitations. Firstly, 16 S rDNA sequencing rather than metagenome sequencing limits the further interpretation of species precision and function analysis in this study. Secondly, although strict inclusion and exclusion criteria were set, individual differences were inevitable, and there were always confounding factors such as genetics, diet, environment, etc. Third, the subjects only included children with OSA or non-OSA. It is necessary to further evaluate the impact of upper airway microbiome on adult OSA and include healthy people as controls. Fourth, this was a single center study. Most of the samples were obtained from nearby regions, and the sample size was comparatively limited. Consequently, it is imperative to conduct a multi-center study encompassing a substantial sample size to corroborate these findings. Finally, this study only provided evidence of correlation rather than causality. It is necessary to further experiment to verify whether candidate bacteria and their functional pathways related to pediatric OSA. Despite these limitations, our study offers valuable insights for future research exploring floral disturbances in the upper airway of pediatric patients with OSA.

Conclusion

In summary, our study demonstrated that changes in upper airway microbiota, especially in the tonsil area, are highly associated with pediatric OSA. The relative abundance of some bacteria were significantly different among OSA groups and control group. The impact of adeno-tonsillectomy on oral microbiota composition is relatively small in the short term. The associations between upper airway microbiota and OSA require further investigations to evaluate causal relationships.

Data availability

Raw data in this study have been deposited into the Sequence Read Archive (SRA) and are publicly available under the BioProject accession number of PRJNA100021.

References

  1. Mitchell RB, Archer SM, Ishman SL, Rosenfeld RM, Coles S, Finestone SA, Friedman NR, Giordano T, Hildrew DM, Kim TW, et al. Clinical practice Guideline: Tonsillectomy in Children (Update)-Executive Summary. Otolaryngol Head Neck Surg. 2019;160(2):187–205.

    Article  PubMed  Google Scholar 

  2. Marcus CL, Brooks LJ, Draper KA, Gozal D, Halbower AC, Jones J, Schechter MS, Sheldon SH, Spruyt K, Ward SD, et al. Diagnosis and management of childhood obstructive sleep apnea syndrome. Pediatrics. 2012;130(3):576–84.

    Article  PubMed  Google Scholar 

  3. Tan H-L, Gozal D, Kheirandish-Gozal L. Obstructive sleep apnea in children: a critical update. Nat Sci Sleep. 2013;5:109–23.

    PubMed  PubMed Central  Google Scholar 

  4. Sigurdardottir SL, Thorleifsdottir RH, Valdimarsson H, Johnston A. The role of the palatine tonsils in the pathogenesis and treatment of psoriasis. Br J Dermatol. 2013;168(2):237–42.

    Article  CAS  PubMed  Google Scholar 

  5. Winther B, Gross BC, Hendley JO, Early SV. Location of bacterial biofilm in the mucus overlying the adenoid by light microscopy. Arch Otolaryngol Head Neck Surg. 2009;135(12):1239–45.

    Article  PubMed  Google Scholar 

  6. Kania RE, Lamers GEM, Vonk MJ, Huy PTB, Hiemstra PS, Bloemberg GV, Grote JJ. Demonstration of bacterial cells and glycocalyx in biofilms on human tonsils. Arch Otolaryngol Head Neck Surg. 2007;133(2):115–21.

    Article  PubMed  Google Scholar 

  7. Peng X, Cheng L, You Y, Tang C, Ren B, Li Y, Xu X, Zhou X. Oral microbiota in human systematic diseases. Int J Oral Sci. 2022;14(1):14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Nistico L, Kreft R, Gieseke A, Coticchia JM, Burrows A, Khampang P, Liu Y, Kerschner JE, Post JC, Lonergan S, et al. Adenoid reservoir for pathogenic biofilm bacteria. J Clin Microbiol. 2011;49(4):1411–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Stępińska M, Olszewska-Sosińska O, Lau-Dworak M, Zielnik-Jurkiewicz B, Trafny EA. Identification of intracellular bacteria in adenoid and tonsil tissue specimens: the efficiency of culture versus fluorescent in situ hybridization (FISH). Curr Microbiol. 2014;68(1):21–9.

    Article  PubMed  Google Scholar 

  10. Ren T, Glatt DU, Nguyen TN, Allen EK, Early SV, Sale M, Winther B, Wu M. 16S rRNA survey revealed complex bacterial communities and evidence of bacterial interference on human adenoids. Environ Microbiol. 2013;15(2):535–47.

    Article  CAS  PubMed  Google Scholar 

  11. Marazzato M, Zicari AM, Aleandri M, Conte AL, Longhi C, Vitanza L, Bolognino V, Zagaglia C, De Castro G, Brindisi G, et al. 16S Metagenomics reveals dysbiosis of nasal core microbiota in children with chronic nasal inflammation: role of adenoid hypertrophy and allergic Rhinitis. Front Cell Infect Microbiol. 2020;10:458.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. World Medical Association. Declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA. 2013;310(20):2191–4.

    Article  Google Scholar 

  13. Chinese guideline for the diagnosis and treatment of childhood obstructive sleep apnea. (2020). World J Otorhinolaryngol Head Neck Surg 2021, 7(3):201–220.

  14. Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.

    Article  CAS  PubMed  Google Scholar 

  16. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.

    Article  CAS  PubMed  Google Scholar 

  19. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wilkinson TJ, Huws SA, Edwards JE, Kingston-Smith AH, Siu-Ting K, Hughes M, Rubino F, Friedersdorff M, Creevey CJ. CowPI: a Rumen Microbiome Focussed Version of the PICRUSt Functional Inference Software. Front Microbiol. 2018;9:1095.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Liu J, Duan Y. Saliva: a potential media for disease diagnostics and monitoring. Oral Oncol. 2012;48(7):569–77.

    Article  PubMed  Google Scholar 

  22. Javaid MA, Ahmed AS, Durand R, Tran SD. Saliva as a diagnostic tool for oral and systemic diseases. J Oral Biol Craniofac Res. 2016;6(1):66–75.

    Article  PubMed  Google Scholar 

  23. Chen W-P, Chang S-H, Tang C-Y, Liou M-L, Tsai S-JJ, Lin Y-L. Composition analysis and feature selection of the oral Microbiota Associated with Periodontal Disease. Biomed Res Int. 2018;2018:3130607.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Mok SF, Karuthan C, Cheah YK, Ngeow WC, Rosnah Z, Yap SF, Ong HKA. The oral microbiome community variations associated with normal, potentially malignant disorders and malignant lesions of the oral cavity. Malays J Pathol 2017, 39(1).

  25. Farrugia C, Stafford GP, Murdoch C. Outer membrane vesicles increase vascular permeability. J Dent Res. 2020;99(13):1494–501.

    Article  CAS  PubMed  Google Scholar 

  26. Liu X, Tian K, Ma X, Wang S, Luo C, Du Q. Analysis of subgingival microbiome of periodontal disease and rheumatoid arthritis in Chinese: a case-control study. Saudi J Biol Sci. 2020;27(7):1835–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Li J, Chang S, Guo H, Ji Y, Jiang H, Ruan L, Du M. Altered salivary microbiome in the early stage of HIV infections among Young Chinese men who have sex with men (MSM). Pathogens 2020, 9(11).

  28. Yang W, Shao L, Heizhati M, Wu T, Yao X, Wang Y, Wang L, Li N. Oropharyngeal Microbiome in Obstructive Sleep Apnea: decreased diversity and abundance. J Clin Sleep Med. 2019;15(12):1777–88.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Jervis-Bardy J, Rogers GB, Morris PS, Smith-Vaughan HC, Nosworthy E, Leong LEX, Smith RJ, Weyrich LS, De Haan J, Carney AS, et al. The microbiome of otitis media with effusion in indigenous Australian children. Int J Pediatr Otorhinolaryngol. 2015;79(9):1548–55.

    Article  PubMed  Google Scholar 

  30. Swidsinski A, Göktas O, Bessler C, Loening-Baucke V, Hale LP, Andree H, Weizenegger M, Hölzl M, Scherer H, Lochs H. Spatial organisation of microbiota in quiescent adenoiditis and tonsillitis. J Clin Pathol. 2007;60(3):253–60.

    Article  CAS  PubMed  Google Scholar 

  31. Wu BG, Sulaiman I, Wang J, Shen N, Clemente JC, Li Y, Laumbach RJ, Lu S-E, Udasin I, Le-Hoang O et al. Severe obstructive sleep apnea is Associated with alterations in the nasal microbiome and an increase in inflammation. Am J Respir Crit Care Med 2019, 199(1).

  32. Maiden MC. Population genomics: diversity and virulence in the Neisseria. Curr Opin Microbiol. 2008;11(5):467–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Baerentsen R, Tang CM, Exley RM. Et Tu, Neisseria? Conflicts of interest between Neisseria species. Front Cell Infect Microbiol. 2022;12:913292.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. John CM, Liu M, Phillips NJ, Yang Z, Funk CR, Zimmerman LI, Griffiss JM, Stein DC, Jarvis GA. Lack of lipid A pyrophosphorylation and functional lptA reduces inflammation by Neisseria commensals. Infect Immun. 2012;80(11):4014–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Escobar A, Rodas PI, Acuña-Castillo C. Macrophage-Neisseria gonorrhoeae interactions: a Better Understanding of Pathogen Mechanisms of Immunomodulation. Front Immunol. 2018;9:3044.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Zhu W, Cardenas-Alvarez MX, Tomberg J, Little MB, Duncan JA, Nicholas RA. Commensal Neisseria species share immune suppressive mechanisms with Neisseria gonorrhoeae. PLoS ONE. 2023;18(4):e0284062.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Brook I, Yocum P, Friedman EM. Aerobic and anaerobic bacteria in tonsils of children with recurrent tonsillitis. Ann Otol Rhinol Laryngol. 1981;90(3 Pt 1):261–3.

    Article  CAS  PubMed  Google Scholar 

  38. Zhang X, Wang S, Xu H, Yi H, Guan J, Yin S. Metabolomics and microbiome profiling as biomarkers in obstructive sleep apnoea: a comprehensive review. Eur Respir Rev 2021, 30(160).

  39. Tahrani AA, Ali A, Raymond NT, Begum S, Dubb K, Mughal S, Jose B, Piya MK, Barnett AH, Stevens MJ. Obstructive sleep apnea and diabetic neuropathy: a novel association in patients with type 2 diabetes. Am J Respir Crit Care Med. 2012;186(5):434–41.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Boneberg R, Pardun A, Hannemann L, Hildebrandt O, Koehler U, Kinscherf R, Hildebrandt W. High plasma cystine levels are Associated with blood pressure and reversed by CPAP in patients with obstructive sleep apnea. J Clin Med 2021, 10(7).

  41. Boulet L, Flore P, Le Gouellec A, Toussaint B, Pépin JL, Faure P. Is tryptophan metabolism involved in sleep apnea-related cardiovascular co-morbidities and cancer progression? Med Hypotheses. 2015;85(4):415–23.

    Article  CAS  PubMed  Google Scholar 

  42. Weersma RK, Zhernakova A, Fu J. Interaction between drugs and the gut microbiome. Gut. 2020;69(8):1510–9.

    Article  CAS  PubMed  Google Scholar 

  43. Kim YJ, Womble JT, Gunsch CK, Ingram JL. The Gut/Lung Microbiome Axis in obesity, Asthma, and bariatric surgery: a Literature Review. Obes (Silver Spring). 2021;29(4):636–44.

    Article  Google Scholar 

  44. Lin B, Zhao F, Liu Y, Wu X, Feng J, Jin X, Yan W, Guo X, Shi S, Li Z, et al. Randomized Clinical Trial: Probiotics alleviated oral-gut microbiota dysbiosis and thyroid hormone withdrawal-related complications in thyroid Cancer patients before Radioiodine Therapy following Thyroidectomy. Front Endocrinol (Lausanne). 2022;13:834674.

    Article  PubMed  Google Scholar 

  45. Zaura E, Brandt BW, Teixeira de Mattos MJ, Buijs MJ, Caspers MPM, Rashid M-U, Weintraub A, Nord CE, Savell A, Hu Y, et al. Same exposure but two radically different responses to antibiotics: resilience of the salivary Microbiome versus Long-Term Microbial shifts in feces. mBio. 2015;6(6):e01693–e01615.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Høiby N, Bjarnsholt T, Givskov M, Molin S, Ciofu O. Antibiotic resistance of bacterial biofilms. Int J Antimicrob Agents. 2010;35(4):322–32.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Thank you to all participants who have contributed to this study.

Funding

This project is supported by the National Natural Science Foundation of China (No. 82071033).

Author information

Authors and Affiliations

Authors

Contributions

L.C.F and A.T had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. A.T participated in study conception. L.C.F, W.Y.C recruited patients. L.C.F and X.C drafted and revised the paper. All authors approved the final draft of the manuscript for publication.

Corresponding author

Correspondence to Xiong Chen.

Ethics declarations

Ethical approval

This study was approved by the Ethics Committee of Zhongnan Hospital of Wuhan University (No.2021056). Informed consent was obtained from all subjects involved in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fang, L., Tuohuti, A., Cai, W. et al. Changes in the nasopharyngeal and oropharyngeal microbiota in pediatric obstructive sleep apnea before and after surgery: a prospective study. BMC Microbiol 24, 79 (2024). https://doi.org/10.1186/s12866-024-03230-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03230-7

Keywords