Next Article in Journal
EPR-Effect Enhancers Strongly Potentiate Tumor-Targeted Delivery of Nanomedicines to Advanced Cancers: Further Extension to Enhancement of the Therapeutic Effect
Next Article in Special Issue
Exploring the Oral Microbiome in Rheumatic Diseases, State of Art and Future Prospective in Personalized Medicine with an AI Approach
Previous Article in Journal
Improving the Pharmacological Properties of Ciclopirox for Its Use in Congenital Erythropoietic Porphyria
Previous Article in Special Issue
Rifaximin as a Potential Treatment for IgA Nephropathy in a Humanized Mice Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 6
10.3390/jpm11060486
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Different Associations between Tonsil Microbiome, Chronic Tonsillitis, and Intermittent Hypoxemia among Obstructive Sleep Apnea Children of Different Weight Status: A Pilot Case-Control Study

by
Hai-Hua Chuang
1,2,3,4,
Jen-Fu Hsu
2,5,
Li-Pang Chuang
2,6,
Cheng-Hsun Chiu
2,5,
Yen-Lin Huang
2,7,
Hsueh-Yu Li
2,8,
Ning-Hung Chen
2,6,
Yu-Shu Huang
2,9,
Chun-Wei Chuang
10,
Chung-Guei Huang
10,11,
Hsin-Chih Lai
10,11 and
Li-Ang Lee
2,8,*
1
Department of Family Medicine, Chang Gung Memorial Hospital, Taipei Branch and Linkou Main Branch, Taoyuan 33305, Taiwan
2
College of Medicine, Chang Gung University, Taoyuan 33302, Taiwan
3
Department of Industrial Engineering and Management, National Taipei University of Technology, Taipei 10608, Taiwan
4
Obesity Institute, Genomic Medicine Institute, Geisinger, Danville, PA 17822, USA
5
Department of Pediatrics, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
6
Department of Pulmonary and Critical Care Medicine, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
7
Department of Pathology, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
8
Department of Otorhinolaryngology-Head and Neck Surgery, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
9
Department of Child Psychiatry, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
10
Department of Laboratory Medicine, Chang Gung Memorial Hospital, Linkou Main Branch, Taoyuan 33305, Taiwan
11
Department of Medical Biotechnology and Laboratory Science, Graduate Institute of Biomedical Sciences, Chang Gung University, Taoyuan 33302, Taiwan
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(6), 486; https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11060486
Submission received: 4 May 2021 / Revised: 20 May 2021 / Accepted: 26 May 2021 / Published: 28 May 2021
(This article belongs to the Special Issue Targeting the Microbiome for Disease Diagnosis and Therapy: New Frontiers for Personalized Medicine)
Browse Figures
Review Reports Versions Notes

Abstract

:
The tonsil microbiome is associated with chronic tonsillitis and obstructive sleep apnea (OSA) in children, and the gut microbiome is associated with host weight status. In this study, we hypothesized that weight status may be associated with clinical profiles and the tonsil microbiome in children with OSA. We prospectively enrolled 33 non-healthy-weight (cases) and 33 healthy-weight (controls) pediatric OSA patients matched by the proportion of chronic tonsillitis. Differences in the tonsil microbiome between the non-healthy-weight and healthy-weight subgroups and relationships between the tonsil microbiome and clinical variables were investigated. Non-healthy weight was associated with significant intermittent hypoxemia (oxygen desaturation index, mean blood saturation (SpO2), and minimal SpO2) and higher systolic blood pressure percentile, but was not related to the tonsil microbiome. However, chronic tonsillitis was related to Acidobacteria in the non-healthy-weight subgroup, and oxygen desaturation index was associated with Bacteroidetes in the healthy-weight subgroup. In post hoc analysis, the children with mean SpO2 ≤ 97% had reduced α and β diversities and a higher abundance of Bacteroidetes than those with mean SpO2 > 97%. These preliminary findings are novel and provide insights into future research to understand the pathogenesis of the disease and develop personalized treatments for pediatric OSA.

1. Introduction

Obstructive sleep apnea (OSA) is a chronic disorder characterized by intermittent partial or complete upper airway obstruction during sleep. The prevalence of pediatric OSA is estimated to be 1–4%, with adenotonsillar hypertrophy and overweight/obesity being the two most important risk factors [1,2,3]. Pediatric OSA is of great clinical significance since evidence has shown a wide range of detrimental long-term effects associated with the condition [1]. For example, children with OSA show higher risks of neurobehavioral impairment [4], metabolic alterations [5], and cardiovascular dysfunction [6].
The role of microbiota in the development and aggravation of OSA has gained increasing attention. Previous studies have reported that the gut microbiota is involved in the pathogenesis of OSA [7,8], obesity [9,10], and hypertension [11,12]. For example, the transplant of fecal microbiota has been shown to elicit sleep disturbance [8], obesity [13], and hypertension [12] in animal models. The gut microbiota has also been associated with intermittent hypoxia and systemic inflammation [7,10], which are both well-documented manifestations of OSA [14,15,16]. More recent studies have suggested that, in addition to the gut microbiota, OSA is linked to alterations in various other microbiomes in the human body such as the nasal cavity [17], adenoids [18], tonsils [19], oropharynx [20], oral cavity [21], lungs [22], and urine [21]. In children who snore, the adenotonsillar microbiome has been shown to interact with the regional mucosal immune system such as interleukin-8 and heat shock protein 27 [23].
Tonsil size is one of the most important predictors for apnea-hypopnea index (AHI) in preschoolers and school-age children [24]. Therefore, the influence of the tonsil microbiome may be significant on OSA in young children. Two main methods are used to detect bacterial communities on tonsils: swab cultures and culture-free molecular tests based on 16S ribosomal RNA or ribosomal DNA sequencing [25]. Notably, molecular tests enable metagenomic studies to better detect slow-growing, uncultivable, and rare bacteria [26]. The advent of metagenomics has led to an increase in investigations on human microbiota. However, studies on the tonsil microbiomes in pediatric OSA patients and their relationships with patient characteristics, disease severity, and hypertension are still lacking. To the best of our knowledge, the clinical significance of tonsil microbiota in children with OSA has not been comprehensively elucidated.
We hypothesized that the tonsil microbiome may be associated with the weight status and anthropometrics of pediatric OSA patients. Furthermore, the relationships between the tonsil microbiome, OSA severity, intermittent hypoxemia, and hypertension may differ across patients with various demographic and clinical parameters. Therefore, among a cohort of children with OSA, the first aim of this study was to investigate differences in the tonsil microbiome between non-healthy-weight and healthy-weight subgroups. The second aim was to perform post hoc analysis to understand the correlations between the tonsil microbiome and other variables of interest, including OSA severity, intermittent hypoxemia, and hypertension.

2. Materials and Methods

2.1. Ethical Considerations

This was a prospective case-control study. Consecutive pediatric patients referred to the Department of Otolaryngology at Chang Gung Memorial Hospital (Linkou Main Branch, Taoyuan, Taiwan) for adenotonsillectomy between 1 March 2017 and 31 January 2019 were recruited. The Institutional Review Board of Chang Gung Medical Foundation approved this study (201507279A3), and all procedures were conducted in compliance with the Declaration of Helsinki 1975. Written informed consent was obtained from all parents and participants ≥6 years of age.

2.2. Patient Selection and Grouping

All of the participants underwent comprehensive history-taking, physical examinations, and standard in-lab polysomnography (PSG). The protocol was previously published [15]. The inclusion criteria were: (1) age 5–12 years, and (2) AHI ≥ 5.0 events/h or AHI ≥ 2.0 events/h plus at least one morbidity (such as elevated blood pressure (BP), daytime sleepiness, learning problems, growth failure, or enuresis) [14,27,28]. Patients with craniofacial, neuromuscular, or chronic inflammatory disorders (such as atopic dermatitis, asthma, or autoimmune disease) were excluded [14,15]. The subjects were further divided into two subgroups according to body mass index (BMI) z-score: “non-healthy-weight” (≤−2.0 kg/m2 and ≥1.0 kg/m2) group, and “healthy-weight” (>−2.0 kg/m2 and <1.0 kg/m2) group [29]. Both groups were matched by the proportion of chronic tonsillitis. Chronic tonsillitis was defined as symptoms of tonsillitis that persisted for a period longer than three months [30]. Patients with acute inflammation, such as rhinosinusitis, tonsillitis, gastrointestinal infection, or other conditions that needed antibiotic treatment did not undergo surgery after the diseases diminished for at least 2 weeks [15]. Subjective OSA symptoms (evaluated using the Chinese version of the OSA-18 questionnaire [31,32]), tonsil size (rated using the Brodsky grading scale [33]), the adenoidal-nasopharyngeal ratio (ANR) (measured using lateral radiography of the nasopharynx [34]), and allergic rhinitis were recorded. Figure 1 shows the flow diagram of the study.

2.3. Polysomnography Variables

We assessed OSA severity variables (AHI, respiratory disturbance index (RDI), oxygen desaturation index (ODI), mean pulse oxygen saturation (SpO2), and minimal SpO2) by standard full-night, in-lab PSG, according to the 2012 American Academy of Sleep Medicine Manual [35]. Briefly, the AHI was defined as the sum of all apneas (≥90% decrease in airflow for a duration of ≥2 breaths) plus hypopneas (≥50% decrease in airflow and either ≥ 3% desaturation or electroencephalographic arousal, for a duration of ≥2 breaths) divided by the number of hours of total sleep time. The patients were then categorized as having either severe (obstructive AHI ≥ 10.0 events/h) or non-severe (obstructive AHI ≥ 2.0 events/h to < 10.0 events/h) OSA [36]. The RDI was defined as the average number of respiratory disturbances (obstructive apneas, hypopneas, and respiratory-event-related arousals) per hour. The ODI was calculated as the average number of respiratory events with a 3% drop in SpO2 per hour. Furthermore, sleep stages were scored from electroencephalography records according to conventional criteria. PSG scoring was performed by a technician blinded to the clinical status of the children. Details of the PSG protocol were described previously [15,37].

2.4. Nocturnal Blood Pressure

Before the PSG exam, nocturnal BP was measured three times with a standard sphygmomanometer between 10:00 and 11:00 PM. The detailed procedure of measuring BP is described elsewhere [38]. Age, sex, and height-corrected percentiles of systolic BP (SBP) and diastolic BP (DBP) were recorded for each child [39]. Pediatric hypertension was defined as an average clinic SBP and/or DBP ≥ 95th percentile [39].

2.5. Tonsil Microbiota

The current study chose tonsil as the primary site for investigation since previous research suggested that the adenoidal microbiome was compatible with the tonsillar microbiome at the phylum level [23]. Also, since bacterial colonies were mainly observed in the tonsil crypts rather than in the tonsil follicles, the superficial tonsils with crypts were used for molecular examinations.
Tonsils with crypts were excised using sterile scissors during adenotonsillectomy. The specimens were rinsed with normal saline to remove superficial debris several times after harvesting. Genomic DNA was immediately extracted from the superficially biopsied specimens (3 mm × 3 mm × 3 mm) using an EasyPrep Genomic DNA Extraction Kit (Biotools Co., Ltd., New Taipei, Taiwan). Tonsil tissue was treated with 4 μL of RNase A (100 mg/mL) for 5 min at room temperature followed by 20 μL of Proteinase K at 56 °C until completely lysed, and then 200 μL ethanol (96–100%) for 15 s [40]. The quality and quantity of genomic DNA were measured using a NanoPhotometer P360 system (Implen, Westlake Village, CA, USA). Polymerase chain reaction (PCR) was used to amplify the V3–V4 regions of the gene encoding for 16S rRNA in bacteria using composite primers, including the forward primer 5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTAYGGGRBGCASCAG-3′ and the reverse primer 5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGGACTACNNGGGTATCTAAT-3′ [41]. Amplicons were purified using a QiaQuick PCR Purification Kit (Qiagen, Hilden, Germany). PCR amplicons were sequenced using the Illumina HiSeq 2500 platform (Illumina, Inc., San Diego, CA, USA) following the manufacturer’s instructions to generate 250 bp paired-end reads.
All of the paired-end reads were assembled using FLASH software (version 1.2.7; http://0-ccb-jhu-edu.brum.beds.ac.uk/software/FLASH/ (accessed on 5 February 2019)) [42], and reads with a quality score < 20 were removed using QIIME software (version 1.7; http://qiime.org/ (accessed on 5 February 2019)) [43]. Sequences were chimera-checked using UCHIME software (http://drive5.com/uchime/ (accessed on 5 February 2019)) [44] and filtered from the dataset before the operational taxonomic unit (OTU) picking of 97% sequence identity using USEARCH (version 7) [45]. Taxonomy classification was annotated according to the Greengenes database (version 13.8; http://greengenes.secondgenome.com/ (accessed on 5 February 2019)) [46]. Multiple sequences were aligned using PyNAST software (version 1.2; https://pypi.org/project/pynast/ (accessed on 5 February 2019)) against the Greengenes core set database to identify the relationships between different OTUs [47]. We used Graphical Phylogenetic Analysis to visualize microbial genomes and metagenomes [48]. Detailed protocols of bioinformatics were described previously [41,49]. During data collection and analysis, the investigators were blinded to group allocation.

2.6. Sample Size Estimation

The sample size was estimated using primary outcome effects (BMI z-score) based on a priori study criteria [15] (healthy-weight group = 0.52 ± 1.12 and non-healthy-weight group = 1.53 ± 1.05). We used a two-tailed Wilcoxon–Mann–Whitney test to calculate the sample size (effect size = 0.93; type I error = 0.05; power = 0.95), which generated a sample size of 33 in each group.

2.7. Statistical Analysis

The D’Agostino and Pearson normality test showed that most variables had non-normal distribution. Therefore, descriptive statistics were expressed as the median, interquartile range (IQR), or frequency. Differences in variables of interest between specific (weight, OSA severity, BP) subgroups were determined using the Mann–Whitney U test, Kruskal–Wallis test, or chi-square test as appropriate. Analysis of similarity was performed to compare bacterial communities [50]. The α diversity (the diversity within each sample) of the tonsil sample was calculated using the observed richness based on the frequency of OTUs and genera in the sequence collections [51]. The β diversity (the number of species shared between two groups) was calculated using the weighted UniFrac measure [52]. Spearman’s correlation test was used to determine associations among major (>0.1% abundance and present in >90% of samples [49]) or minor phyla with the patient characteristics. Overall taxonomic or phylum-level abundances were included when determining the most discriminatory taxa between the two groups. Statistical significance was established at p < 0.05. p-values were corrected for multiple comparisons using the Benjamini–Hochberg method at 0.1 and reported as q-values when appropriate [53]. All statistical analyses were conducted using R software (versions 2.15.3 and 3.6.1, R Foundation for Statistical Computing, Vienna, Austria; http://www.r-project.org/ (accessed on 26 February 2021)) and Graph Pad Prism software (version 9.00; Graph Pad Software Inc., San Diego, CA, USA).

3. Results

3.1. Participants’ Characteristics

Figure 1 demonstrates the study flow diagram. Seventy-six Taiwanese children of Han ancestry with OSA were assessed for eligibility, 10 of whom were excluded from this study. Therefore, a total of 66 children with OSA (16 girls and 50 boys; median age, 6.5 years (IQR, 6.0–9.0); median BMI, 17.1 kg/m2 (IQR, 15.2–22.7); median AHI, 8.5 events/hour (IQR, 4.1–19.5); median SBP, 103 mmHg (IQR, 95–114); median DBP, 64 mmHg (IQR, 58–71)) were enrolled.
The children were further divided into two subgroups according to BMI z-score: non-healthy-weight subgroup (cases; n = 33), and healthy-weight subgroup (controls; n = 33). In the non-healthy-weight subgroup, 30 (91%) children had a BMI z-score ≥ 1.0 kg/m2 and three (9%) had a BMI z-score ≤ −2.0 kg/m2. All participants underwent adenotonsillectomy and were included for primary statistical analysis. The median time interval between PSG and adenotonsillectomy was 1 week (IQR: 6–19 weeks).

3.2. Differences in Participants’ Characteristics, PSG Variables, BP, and Tonsil Microbiome between the Different Weight Status Subgroups

3.2.1. Differences in Participants’ Characteristics, PSG Variables, and BP between the Different Weight Status Subgroups

As expected, there was no significant difference in the proportion of chronic tonsillitis between the two subgroups (Table 1). Furthermore, there were no statistically significant differences in the proportions of male sex, allergic rhinitis, tonsil size, ANR, OSA-18 score, AHI, RDI, N1 stage, N2 stage, N3 stage, rapid eye movement (REM) stage, and DBP percentile. Notably, the non-healthy-weight group had significantly higher age, BMI z-score, ODI, SBP, DBP, and SBP percentile, and lower mean SpO2 and minimal SpO2 than the healthy-weight group. The difference in the tonsil size between children with chronic tonsillitis (3 (IQR, 3–3)) and children without chronic tonsillitis (3 (IQR, 3–4)) did not reach a statistical significance (p = 0.13).
Furthermore, the differences in participants’ characteristics, PSG variables, and BP between the overweight and underweight subgroups were not statistically significant (Supplementary Table S1). These variables were also comparable across the overweight, underweight, and healthy-weight subgroups (all p > 0.05).

3.2.2. Differences in Tonsil Microbiome between the Different Weight Status Subgroups

After quality assessment, a total of 4,207,400 16S rRNA paired-end reads with an average of 105,185 ± 19,957 paired-end reads per sample passed the filters. Figure 2A,B show OTU trees of the non-healthy-weight subgroup (OTU = 9318) and healthy-weight subgroup (OTU = 9886). Figure 2C shows that both subgroups shared 6539 OTUs; otherwise, there were 2779 and 3347 deferential OTUs in the non-healthy-weight and healthy-weight subgroups, respectively. However, there were no significant differences in α diversity, β diversity, and relative abundances of the top 10 tonsil families between the non-healthy-weight and healthy-weight subgroups (Figure 2D–F; all p > 0.05). Furthermore, the α diversity, β diversity, and relative abundance of the top 10 tonsil families of the over-weight were comparable with those of the underweight subgroups (p = 0.064, 0.106, 0.492, respectively). Moreover, the differences in α diversity, β diversity, and relative abundances of the top 10 tonsil families across the overweight, underweight, and healthy-weight subgroups (p = 0.088, 0.119, 0.700, respectively).
Fifty-five phyla were identified from the tonsil samples. There were six major phyla (>0.1% abundance and present in >90% of the samples [49]), including Proteobacteria, Firmicutes, Bacteroidetes, Fusobacteria, Actinobacteria, and Epsilonbacteraeota, and 49 minor phyla. In descending order of median relative abundance, the 10 most common phyla were Proteobacteria, Firmicutes, Bacteroidetes, Fusobacteria, Actinobacteria, Epsilonbacteraeota, Patescibacteria, Cyanobacteria, Tenericutes, and Acidobacteria.
Figure 3A,B demonstrate the similar distributions of the relative abundances of these phyla in the non-healthy-weight and healthy-weight subgroups (p > 0.05). Furthermore, the relative abundances of these phyla in the overweight and underweight subgroups were comparable (p = 0.536). Additionally, the differences in the 10 most common phyla across the overweight, underweight, and healthy-weight subgroups were not statistically significant (p = 0.814).

3.2.3. Associations of Tonsil Phyla, Participants’ Characteristics, PSG Variables, and BPs in the Different Weight Status Subgroups

In the non-healthy-weight group, age was related to Cyanobacteria and Acidobacteria, chronic tonsillitis was correlated with Acidobacteria, and SBP percentile was associated with Firmicutes (Figure 3C). The positive relationship between chronic tonsillitis and Acidobacteria remained significant after applying the Benjamini–Hochberg method (r = 0.53, q = 0.015).
Although there were several weak associations between the PSG variables (AHI, RDI, ODI, mean SpO2, minimal SpO2), DBP percentile, and Bacteroidetes, ODI and Actinobacteria, stage 1 sleep and Firmicutes, stage 3 sleep and Tenericutes, and REM stage and Proteobacteria, only the significant positive association between ODI and Bacteroidetes persisted after applying the Benjamini–Hochberg method (r = 0.52, q = 0.020) (Figure 3D).

3.3. Post Hoc Analysis

After studying the differences between the two weight subgroups, we wondered whether other classifications may be associated with the tonsil microbiome. We performed median splits of the participants’ characteristics, PSG variables, and BPs. The α diversities of the tonsil microbiome were not associated with age ≥ 6 years, male sex, chronic tonsillitis, BMI z-score ≥ 1.00, AHI ≥ 9.0 events/h, RDI ≥ 9.0 events/h, ODI ≥ 5.0 events/h, minimal SpO2 ≤ 90%, N1 stage ≥ 10%, N2 stage ≥ 40%, N3 stage ≤ 26%, REM stage ≥ 19%, SBP percentile ≥ 70%, and DSBP percentile ≥ 70%. In addition, differences in the tonsil microbiome between various (moderate-to-severe OSA and mild OSA; severe OSA, and non-severe OSA; hypertension and non-hypertension) subgroups did not reach statistical significance.
Interestingly, both α (Figure 4A) and β (Figure 4B) diversity indices of the mean SpO2 ≤ 97% subgroup were significantly lower than those of the mean SpO2 > 97% subgroup (p = 0.030 and 0.0005, respectively). The relative abundances of the top 10 tonsil phyla in the mean SpO2 ≤ 97% subgroup were significantly different from those in the mean SpO2 > 97% subgroup (p = 0.014; analysis of similarity test) (Figure 4C). Notably, the relative abundance of Bacteroidetes in the mean SpO2 ≤ 97% subgroup was significantly higher than that in the mean SpO2 > 97% subgroup (q = 0.026) (Figure 4D).

4. Discussion

In the following paragraphs, we address several novel and interesting findings of this study regarding the relationships of the tonsil microbiome with weight status, OSA severity, and hypoxemia among a sample of pediatric OSA patients.
Some clinical parameters were significantly different between the non-healthy-weight and healthy-weight subgroups. The non-healthy-weight subgroup had more profound intermittent hypoxemia (ODI, mean SpO2, and minimal SpO2) and a higher SBP percentile, which is consistent with the well-known connections between obesity and manifestations of OSA. However, no significant association was observed between the tonsil microbiome and weight status. The data did not support our primary hypothesis. Notably, although weight status was not directly associated with the tonsil microbiome in the overall cohort, relationships between the tonsil microbiome and other clinical parameters differed between the patients with different weight statuses. Chronic tonsillitis was related to Acidobacteria in the non-healthy-weight subgroup, while ODI was associated with Bacteroidetes in the healthy-weight subgroup. In post hoc analysis, we found that the children with or without a mean SpO2 ≤ 97% had significantly different microbial profiles, especially with regards to Bacteroidetes. It seemed that, instead of weight status, hypoxemia status was the key differentiating factor for the tonsil microbiota among pediatric OSA patients.
Previous studies have suggested that the presence and severity of OSA are associated with microbial profiles [19,20,23]. Yang et al. used 16S ribosomal DNA sequencing to investigate the oropharyngeal microbiome and demonstrated that adults with OSA had less oropharyngeal species diversity and altered abundance compared with non-OSA controls, and that the relative abundance of Neisseria (a genus of Proteobacteria) increased with higher OSA severity [20]. In our results, neither Proteobacteria nor Neisseria was related to AHI (p = 0.905 and 0.246, respectively). However, the children with a mean SpO2 ≤ 97% had less tonsillar species diversity and altered abundance compared to those with a mean SpO2 > 97%. Again, this suggested that the degree of intermittent hypoxemia may be more influential on the tonsil microbiome than OSA severity.
Johnston et al. investigated the tonsillar crypt microbiota in children with recurrent tonsillitis [54] and OSA [19], and found that Fusobacteria, Proteobacteria, Bacteroidetes, and Firmicutes were major phyla of the tonsil specimens in both groups (α diversity: p = 0.66; β diversity: p = 0.52) [19]. In another study on the adenotonsillar microbiome of children who snored [23], Kim and colleagues also reported similar major phyla, namely Proteobacteria, Actinobacteria, Firmicutes, Fusobacteria, Bacteroidetes, and Tenericutes. They further demonstrated that α diversity indices were related to some patient characteristics such as sex, emotional stress, and interleukin-8, and that β diversity indices were related to heat shock protein 70. These findings suggest possible connections between demographic characteristics, clinical symptoms, regional mucosal immune environment, and tonsil microbiome.
Our data are compatible with the results of the previous studies concerning the general picture of the tonsil microbiome. Moreover, we found a positive correlation between chronic tonsillitis and Acidobacteria in the non-healthy-weight subgroup. To the best of our knowledge, this is the first study to suggest that Acidobacteria may be implicated in tonsillar infections among pediatric OSA patients. Acidobacteria are Gram-negative rod-shaped bacteria. The majority of Acidobacteria strains have been described as aerobes, and they are ubiquitous in soil [55]. An increasing number of studies have investigated Acidobacteria in humans. For example, Acidobacteria are reported to be the fifth most dominant phyla in the bronchoalveolar lavage of adults with OSA [22]. In addition, both Acidobacteria and obesity are associated with an increase in the fecal levels of valeric acid [56,57]. Furthermore, Acidobacteria are shown to be significantly enriched in patients with chronic endodontic infection [58]. These observations provide indirect evidence for possible connections between tonsillar Acidobacteria infection, OSA, and obesity, and future investigations are warranted to confirm the causality.
The other novel finding of our study is the significant association between intermittent hypoxemia and Bacteroidetes in the healthy-weight subgroup. Members of the phylum Bacteroidetes are Gram-negative, rod-shaped, anaerobic or aerobic bacteria, and they are commonly found in the oral cavity and gastrointestinal tract. The abundance of Bacteroidetes and the Bacteroidetes–Firmicutes ratio are reported to be decreased in obese individuals compared to lean individuals [59]. However, our data suggested that it was the healthy-weight pediatric OSA patients in whom ODI significantly interacted with Bacteroidetes. Bacteroidetes of the tonsils may thus be implicated in the pathophysiology of OSA.
Previous studies have suggested associations between OSA and alterations in the composition and diversity of fecal microbiota. In a murine model, intermittent hypoxia exposure led to a lower abundance of Bacteroidetes in the feces [60], and reintroduction of a normoxic environment did not reverse the negative alterations of the gut microbiota [61]. In a sample of adults with OSA, Ko et al. found that fecal Bacteroides (a major genus of Bacteroidetes) were not associated with AHI or intermittent hypoxemia [62]. However, the results were very different for the oral microbiome. Xu et al. reported a higher abundance of oral Bacteroidetes in children with OSA [21]. Consistent with this finding, we also found a higher abundance of Bacteroidetes in pediatric OSA patients with a lower mean SpO2. The discrepancy between fecal and oral microbiota as well as the differences between patients with and without profound hypoxia are very interesting and may be explained by the oxygen concentration level at different sites and the aerotolerant abilities of different bacteria.
The environments of sites in the human body impact which microorganisms can inhabit these sites. Unlike the environment of the intestine, which is extremely low in oxygen concentration [63], to survive in the oral cavity, bacteria need to overcome the challenge of atmospheric oxygen exposure. Moreover, mouth breathing is very common in children with OSA [64] and positively associated with AHI [65]. OSA-related mouth breathing further increases exposure of the tonsil microbiome to atmospheric oxygen. Bacteroides species are among the most aerotolerant anaerobes and are able to tolerate oxygen in room air for up to 3 days [66]. On the other hand, Fusobacterium species (a genus of Fusobacteria) and Clostridium species (a genus of Firmicutes) are less aerotolerant anaerobes. Intermittent hypoxemia (in terms of mean SpO2 ≤ 97%) further enhances the survival advantage of Bacteroidetes relative to other aerobes of the tonsils. Mouse models of chronic intermittent hypoxia would be helpful to further validate these inferences [60].
Several limitations should be addressed in this study. First, the patient number of the underweight subgroup was too small and insufficient to make a conclusion. Future studies with a larger sample size of each weight status subgroup are warranted to further understand how obesity or underweight may impact the tonsil microbiome and its interaction with OSA. Second, the study cohort was predominantly male and mostly Han in ethnicity, which may limit the generalizability of the results. Third, some of the children may have had co-existing chronic tonsillitis, which would interfere with the analysis of the tonsil microbiome. However, the proportions of sex and chronic tonsillitis were comparable in both weight subgroups to minimize confounding effects from baseline characteristics. Forth, the study was cross-sectional and thus unable to conclude the direction of associations or causal effects. The relationships between the microbiota and OSA need to be further explored. Also, future prospective investigations on the effects of OSA treatment on the tonsil microbiome with a larger sample size will be of interest.

5. Conclusions

The advent of metagenomics has led to an increase in investigations on human microbiota. The tonsil microbiome plays a role in pediatric OSA, and it seems to have different effects depending on weight status. We preliminarily found that chronic tonsillitis was related to Acidobacteria in children with OSA and non-healthy weight, and that ODI was associated with Bacteroidetes in the children with OSA and healthy weight. In addition, children with OSA with or without mean SpO2 ≤ 97% had significantly different microbial profiles, particularly with regards to Bacteroidetes. Future studies to investigate associations among alterations of the tonsil microbiome and exacerbations or reductions of OSA severity are warranted. Furthermore, this study also suggests the possibility of personalized treatment of pediatric OSA based on the tonsil microbiome.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/jpm11060486/s1, Table S1: Patient characteristics, polysomnography variables, and blood pressures of the over-weight and under-weight subgroups.

Author Contributions

Conceptualization, H.-H.C., J.-F.H., C.-G.H., H.-C.L., and L.-A.L.; methodology, L.-P.C., Y.-L.H., H.-Y.L., Y.-S.H., C.-W.C., C.-G.H., H.-C.L., and L.-A.L.; software, J.-F.H., N.-H.C., C.-W.C., C.-G.H., and H.-C.L.; validation, H.-H.C., C.-H.C., Y.-L.H., H.-Y.L., N.-H.C., Y.-S.H., C.-G.H., H.-C.L., and L.-A.L.; formal analysis, H.-H.C. and L.-A.L.; investigation, H.-H.C., L.-P.C., C.-H.C., Y.-L.H., H.-Y.L., Y.-S.H., C.-G.H., and L.-A.L.; resources, L.-P.C., C.-H.C., Y.-L.H., Y.-S.H., C.-G.H., H.-C.L., and L.-A.L.; data curation, J.-F.H., C.-H.C., H.-Y.L., N.-H.C., Y.-S.H., C.-G.H., and H.-C.L.; writing—original draft, H.-H.C., J.-F.H., L.-P.C., C.-W.C., and L.-A.L.; writing—review and editing, C.-H.C., Y.-L.H., H.-Y.L., N.-H.C., Y.-S.H., C.-G.H., and H.-C.L.; visualization, H.-H.C. and L.-A.L.; supervision, C.-H.C., H.-Y.L., N.-H.C., Y.-S.H., and H.-C.L.; project administration, L.-A.L.; funding acquisition, H.-H.C. and L.-A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chang Gung Medical Foundation, Taiwan, grant numbers CMRPG1J0041, EMRPD1I0411, and 1K0391 (H.-H.C.), and CMRPG3F1091, 3F1092, 3F1093, and 3J1701 (L.-A.L.). The APC was funded by the Chang Gung Medical Foundation, Taiwan.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Chang Gung Medical Foundation, Taoyuan, Taiwan (protocol code 201507279A3 and date of approval 15 April 2016).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to ethical restrictions.

Acknowledgments

The authors would like to thank Ruo-Chi Wang and Chung-Fang Hsiao (Department of Otorhinolaryngology, Head and Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan City, Taiwan) for their technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Marcus, C.L.; Brooks, L.J.; Draper, K.A.; Gozal, D.; Halbower, A.C.; Jones, J.; Schechter, M.S.; Sheldon, S.H.; Spruyt, K.; Ward, S.D.; et al. Diagnosis and management of childhood obstructive sleep apnea syndrome. Pediatrics 2012, 130, 576–584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Andersen, I.G.; Holm, J.C.; Homoe, P. Obstructive sleep apnea in children and adolescents with and without obesity. Eur. Arch. Otorhinolaryngol. 2019, 276, 871–878. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, J.; Zhao, Y.; Yang, W.; Shen, T.; Xue, P.; Yan, X.; Chen, D.; Qiao, Y.; Chen, M.; Ren, R.; et al. Correlations between obstructive sleep apnea and adenotonsillar hypertrophy in children of different weight status. Sci. Rep. 2019, 9, 11455. [Google Scholar] [CrossRef] [PubMed]
  4. Madaeva, I.; Berdina, O.; Polyakov, V.; Kolesnikov, S. Obstructive Sleep Apnea and Hypertension in Adolescents: Effect on Neurobehavioral and Cognitive Functioning. Can. Respir. J. 2016, 2016, 3950914. [Google Scholar] [CrossRef] [Green Version]
  5. Shalitin, S.; Deutsch, V.; Tauman, R. Hepcidin, soluble transferrin receptor and IL-6 levels in obese children and adolescents with and without type 2 diabetes mellitus/impaired glucose tolerance and their association with obstructive sleep apnea. J. Endocrinol. Investig. 2018, 41, 969–975. [Google Scholar] [CrossRef]
  6. Smith, D.F.; Hossain, M.M.; Hura, A.; Huang, G.; McConnell, K.; Ishman, S.L.; Amin, R.S. Inflammatory Milieu and Cardiovascular Homeostasis in Children With Obstructive Sleep Apnea. Sleep 2017, 40. [Google Scholar] [CrossRef] [Green Version]
  7. Tripathi, A.; Melnik, A.V.; Xue, J.; Poulsen, O.; Meehan, M.J.; Humphrey, G.; Jiang, L.; Ackermann, G.; McDonald, D.; Zhou, D.; et al. Intermittent Hypoxia and Hypercapnia, a Hallmark of Obstructive Sleep Apnea, Alters the Gut Microbiome and Metabolome. mSystems 2018, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Badran, M.; Khalyfa, A.; Ericsson, A.; Gozal, D. Fecal microbiota transplantation from mice exposed to chronic intermittent hypoxia elicits sleep disturbances in naive mice. Exp. Neurol. 2020, 334, 113439. [Google Scholar] [CrossRef]
  9. Turnbaugh, P.J.; Ley, R.E.; Mahowald, M.A.; Magrini, V.; Mardis, E.R.; Gordon, J.I. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 2006, 444, 1027–1031. [Google Scholar] [CrossRef]
  10. Cox, A.J.; West, N.P.; Cripps, A.W. Obesity, inflammation, and the gut microbiota. Lancet Diabetes Endocrinol. 2015, 3, 207–215. [Google Scholar] [CrossRef]
  11. Pevsner-Fischer, M.; Blacher, E.; Tatirovsky, E.; Ben-Dov, I.Z.; Elinav, E. The gut microbiome and hypertension. Curr. Opin. Nephrol. Hypertens. 2017, 26, 1–8. [Google Scholar] [CrossRef] [PubMed]
  12. Durgan, D.J.; Ganesh, B.P.; Cope, J.L.; Ajami, N.J.; Phillips, S.C.; Petrosino, J.F.; Hollister, E.B.; Bryan, R.M., Jr. Role of the Gut Microbiome in Obstructive Sleep Apnea-Induced Hypertension. Hypertension 2016, 67, 469–474. [Google Scholar] [CrossRef] [Green Version]
  13. Ridaura, V.K.; Faith, J.J.; Rey, F.E.; Cheng, J.; Duncan, A.E.; Kau, A.L.; Griffin, N.W.; Lombard, V.; Henrissat, B.; Bain, J.R.; et al. Gut microbiota from twins discordant for obesity modulate metabolism in mice. Science 2013, 341, 1241214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Huang, Y.S.; Guilleminault, C.; Hwang, F.M.; Cheng, C.; Lin, C.H.; Li, H.Y.; Lee, L.A. Inflammatory cytokines in pediatric obstructive sleep apnea. Medicine (Baltimore) 2016, 95, e4944. [Google Scholar] [CrossRef]
  15. Chuang, H.H.; Huang, C.G.; Chuang, L.P.; Huang, Y.S.; Chen, N.H.; Li, H.Y.; Fang, T.J.; Hsu, J.F.; Lai, H.C.; Chen, J.Y.; et al. Relationships Among and Predictive Values of Obesity, Inflammation Markers, and Disease Severity in Pediatric Patients with Obstructive Sleep Apnea Before and After Adenotonsillectomy. J. Clin. Med. 2020, 9, 579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Eltzschig, H.K.; Carmeliet, P. Hypoxia and inflammation. N. Engl. J. Med. 2011, 364, 656–665. [Google Scholar] [CrossRef] [Green Version]
  17. Wu, B.G.; Sulaiman, I.; Wang, J.; Shen, N.; Clemente, J.C.; Li, Y.; Laumbach, R.J.; 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, 99–109. [Google Scholar] [CrossRef]
  18. Dirain, C.O.; Silva, R.C.; Collins, W.O.; Antonelli, P.J. The Adenoid Microbiome in Recurrent Acute Otitis Media and Obstructive Sleep Apnea. J. Int. Adv. Otol. 2017, 13, 333–339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Johnston, J.; Hoggard, M.; Biswas, K.; Astudillo-Garcia, C.; Waldvogel-Thurlow, S.; Radcliff, F.J.; Mahadevan, M.; Douglas, R.G. The bacterial community and local lymphocyte response are markedly different in patients with recurrent tonsillitis compared to obstructive sleep apnoea. Int. J. Pediatric Otorhinolaryngol. 2018, 113, 281–288. [Google Scholar] [CrossRef]
  20. 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. JCSM Off. Publ. Am. Acad. Sleep Med. 2019, 15, 1777–1788. [Google Scholar] [CrossRef] [PubMed]
  21. Xu, H.; Li, X.; Zheng, X.; Xia, Y.; Fu, Y.; Li, X.; Qian, Y.; Zou, J.; Zhao, A.; Guan, J.; et al. Pediatric Obstructive Sleep Apnea is Associated With Changes in the Oral Microbiome and Urinary Metabolomics Profile: A Pilot Study. J. Clin. Sleep Med. 2018, 14, 1559–1567. [Google Scholar] [CrossRef]
  22. Lu, D.; Yao, X.; Abulimiti, A.; Cai, L.; Zhou, L.; Hong, J.; Li, N. Profiling of lung microbiota in the patients with obstructive sleep apnea. Medicine (Baltimore) 2018, 97, e11175. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, K.S.; Min, H.J. Correlation between adenotonsillar microbiome and clinical characteristics of pediatrics with snoring. Clin. Exp. Otorhinolaryngol. 2020. [Google Scholar] [CrossRef]
  24. Chuang, H.H.; Hsu, J.F.; Chuang, L.P.; Chen, N.H.; Huang, Y.S.; Li, H.Y.; Chen, J.Y.; Lee, L.A.; Huang, C.G. Differences in Anthropometric and Clinical Features among Preschoolers, School-Age Children, and Adolescents with Obstructive Sleep Apnea-A Hospital-Based Study in Taiwan. Int. J. Environ. Res. Public Health 2020, 17, 4663. [Google Scholar] [CrossRef]
  25. Johnston, J.J.; Douglas, R. Adenotonsillar microbiome: An update. Postgrad. Med. J. 2018, 94, 398–403. [Google Scholar] [CrossRef] [PubMed]
  26. Verma, D.; Garg, P.K.; Dubey, A.K. Insights into the human oral microbiome. Arch. Microbiol. 2018, 200, 525–540. [Google Scholar] [CrossRef] [PubMed]
  27. Baugh, R.F.; Archer, S.M.; Mitchell, R.B.; Rosenfeld, R.M.; Amin, R.; Burns, J.J.; Darrow, D.H.; Giordano, T.; Litman, R.S.; Li, K.K.; et al. Clinical practice guideline: Tonsillectomy in children. Otolaryngol. Head Neck Surg. 2011, 144, S1–S30. [Google Scholar] [CrossRef] [Green Version]
  28. Kaditis, A.; Kheirandish-Gozal, L.; Gozal, D. Algorithm for the diagnosis and treatment of pediatric OSA: A proposal of two pediatric sleep centers. Sleep Med. 2012, 13, 217–227. [Google Scholar] [CrossRef]
  29. de Onis, M.; Onyango, A.W.; Borghi, E.; Siyam, A.; Nishida, C.; Siekmann, J. Development of a WHO growth reference for school-aged children and adolescents. Bull. World Health Organ. 2007, 85, 660–667. [Google Scholar] [CrossRef]
  30. Burton, M.J.; Glasziou, P.P.; Chong, L.Y.; Venekamp, R.P. Tonsillectomy or adenotonsillectomy versus non-surgical treatment for chronic/recurrent acute tonsillitis. Cochrane Database Syst. Rev. 2014, CD001802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Franco, R.A., Jr.; Rosenfeld, R.M.; Rao, M. First place—Resident clinical science award 1999. Quality of life for children with obstructive sleep apnea. Otolaryngol. Head Neck Surg. Off. J. Am. Acad. Otolaryngol. Head Neck Surg. 2000, 123, 9–16. [Google Scholar] [CrossRef] [PubMed]
  32. Huang, Y.-S.; Hwang, F.-M.; Lin, C.-H.; Lee, L.-A.; Huang, P.-Y.; Chiu, S.-T. Clinical manifestations of pediatric obstructive sleep apnea syndrome: Clinical utility of the Chinese-version Obstructive Sleep Apnea Questionaire-18. Psychiatry Clin. Neurosci. 2015, 69, 752–762. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Brodsky, L. Modern Assessment of Tonsils and Adenoids. Pediatric Clin. N. Am. 1989, 36, 1551–1569. [Google Scholar] [CrossRef]
  34. Fujioka, M.; Young, L.W.; Girdany, B.R. Radiographic evaluation of adenoidal size in children: Adenoidal-nasopharyngeal ratio. AJR Am. J. Roentgenol. 1979, 133, 401–404. [Google Scholar] [CrossRef]
  35. Berry, R.B.; Budhiraja, R.; Gottlieb, D.J.; Gozal, D.; Iber, C.; Kapur, V.K.; Marcus, C.L.; Mehra, R.; Parthasarathy, S.; Quan, S.F.; et al. Rules for scoring respiratory events in sleep: Update of the 2007 AASM Manual for the Scoring of Sleep and Associated Events. Deliberations of the Sleep Apnea Definitions Task Force of the American Academy of Sleep Medicine. J. Clin. Sleep Med. 2012, 8, 597–619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Dehlink, E.; Tan, H.L. Update on paediatric obstructive sleep apnoea. J. Thorac. Dis. 2016, 8, 224–235. [Google Scholar] [CrossRef] [PubMed]
  37. Lu, C.T.; Li, H.Y.; Lee, G.S.; Huang, Y.S.; Huang, C.G.; Chen, N.H.; Lee, L.A. Snoring sound energy as a potential biomarker for disease severity and surgical response in childhood obstructive sleep apnoea: A pilot study. Clin. Otolaryngol. 2019, 44, 47–52. [Google Scholar] [CrossRef]
  38. The fourth report on the diagnosis, evaluation, and treatment of high blood pressure in children and adolescents. Pediatrics 2004, 114, 555–576. [CrossRef]
  39. Flynn, J.T.; Kaelber, D.C.; Baker-Smith, C.M.; Blowey, D.; Carroll, A.E.; Daniels, S.R.; de Ferranti, S.D.; Dionne, J.M.; Falkner, B.; Flinn, S.K.; et al. Clinical Practice Guideline for Screening and Management of High Blood Pressure in Children and Adolescents. cPediatrics 2017, 140, e20171904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Chang, C.C.; Huang, Y.S.; Lin, Y.M.; Lin, C.J.; Jeng, J.C.; Liu, S.M.; Ho, T.L.; Chang, R.T.; Changou, C.A.; Ho, C.C.; et al. The role of sentrin-specific protease 2 substrate recognition in TGF-beta-induced tumorigenesis. Sci. Rep. 2018, 8, 9786. [Google Scholar] [CrossRef]
  41. Wu, T.R.; Lin, C.S.; Chang, C.J.; Lin, T.L.; Martel, J.; Ko, Y.F.; Ojcius, D.M.; Lu, C.C.; Young, J.D.; Lai, H.C. Gut commensal Parabacteroides goldsteinii plays a predominant role in the anti-obesity effects of polysaccharides isolated from Hirsutella sinensis. Gut 2019, 68, 248–262. [Google Scholar] [CrossRef]
  42. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  43. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  44. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  46. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; DeSantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef] [PubMed]
  47. Caporaso, J.G.; Bittinger, K.; Bushman, F.D.; DeSantis, T.Z.; Andersen, G.L.; Knight, R. PyNAST: A flexible tool for aligning sequences to a template alignment. Bioinformatics 2010, 26, 266–267. [Google Scholar] [CrossRef] [Green Version]
  48. Asnicar, F.; Weingart, G.; Tickle, T.L.; Huttenhower, C.; Segata, N. Compact graphical representation of phylogenetic data and metadata with GraPhlAn. PeerJ 2015, 3, e1029. [Google Scholar] [CrossRef] [PubMed]
  49. Wu, I.W.; Lin, C.Y.; Chang, L.C.; Lee, C.C.; Chiu, C.Y.; Hsu, H.J.; Sun, C.Y.; Chen, Y.C.; Kuo, Y.L.; Yang, C.W.; et al. Gut Microbiota as Diagnostic Tools for Mirroring Disease Progression and Circulating Nephrotoxin Levels in Chronic Kidney Disease: Discovery and Validation Study. Int. J. Biol. Sci. 2020, 16, 420–434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Austral Ecol. 1993, 18, 117–143. [Google Scholar] [CrossRef]
  51. Schloss, P.D.; Gevers, D.; Westcott, S.L. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE 2011, 6, e27310. [Google Scholar] [CrossRef] [Green Version]
  52. Lozupone, C.A.; Hamady, M.; Kelley, S.T.; Knight, R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 2007, 73, 1576–1585. [Google Scholar] [CrossRef] [Green Version]
  53. Cornejo-Pareja, I.; Ruiz-Limon, P.; Gomez-Perez, A.M.; Molina-Vega, M.; Moreno-Indias, I.; Tinahones, F.J. Differential Microbial Pattern Description in Subjects with Autoimmune-Based Thyroid Diseases: A Pilot Study. J. Pers. Med. 2020, 10, 192. [Google Scholar] [CrossRef] [PubMed]
  54. Jensen, A.; Fago-Olsen, H.; Sorensen, C.H.; Kilian, M. Molecular mapping to species level of the tonsillar crypt microbiota associated with health and recurrent tonsillitis. PLoS ONE 2013, 8, e56418. [Google Scholar] [CrossRef] [Green Version]
  55. Eichorst, S.A.; Trojan, D.; Roux, S.; Herbold, C.; Rattei, T.; Woebken, D. Genomic insights into the Acidobacteria reveal strategies for their success in terrestrial environments. Environ. Microbiol. 2018, 20, 1041–1063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Liu, S.; Li, E.; Sun, Z.; Fu, D.; Duan, G.; Jiang, M.; Yu, Y.; Mei, L.; Yang, P.; Tang, Y.; et al. Altered gut microbiota and short chain fatty acids in Chinese children with autism spectrum disorder. Sci. Rep. 2019, 9, 287. [Google Scholar] [CrossRef] [PubMed]
  57. Tiihonen, K.; Ouwehand, A.C.; Rautonen, N. Effect of overweight on gastrointestinal microbiology and immunology: Correlation with blood biomarkers. Br. J. Nutr. 2010, 103, 1070–1078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Tzanetakis, G.N.; Azcarate-Peril, M.A.; Zachaki, S.; Panopoulos, P.; Kontakiotis, E.G.; Madianos, P.N.; Divaris, K. Comparison of Bacterial Community Composition of Primary and Persistent Endodontic Infections Using Pyrosequencing. J. Endod. 2015, 41, 1226–1233. [Google Scholar] [CrossRef] [Green Version]
  59. Ley, R.E.; Turnbaugh, P.J.; Klein, S.; Gordon, J.I. Microbial ecology: Human gut microbes associated with obesity. Nature 2006, 444, 1022–1023. [Google Scholar] [CrossRef] [PubMed]
  60. Moreno-Indias, I.; Torres, M.; Montserrat, J.M.; Sanchez-Alcoholado, L.; Cardona, F.; Tinahones, F.J.; Gozal, D.; Poroyko, V.A.; Navajas, D.; Queipo-Ortuno, M.I.; et al. Intermittent hypoxia alters gut microbiota diversity in a mouse model of sleep apnoea. Eur. Respir. J. 2015, 45, 1055–1065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Moreno-Indias, I.; Torres, M.; Sanchez-Alcoholado, L.; Cardona, F.; Almendros, I.; Gozal, D.; Montserrat, J.M.; Queipo-Ortuno, M.I.; Farre, R. Normoxic Recovery Mimicking Treatment of Sleep Apnea Does Not Reverse Intermittent Hypoxia-Induced Bacterial Dysbiosis and Low-Grade Endotoxemia in Mice. Sleep 2016, 39, 1891–1897. [Google Scholar] [CrossRef] [PubMed]
  62. Ko, C.Y.; Fan, J.M.; Hu, A.K.; Su, H.Z.; Yang, J.H.; Huang, L.M.; Yan, F.R.; Zhang, H.P.; Zeng, Y.M. Disruption of sleep architecture in Prevotella enterotype of patients with obstructive sleep apnea-hypopnea syndrome. Brain Behav. 2019, 9, e01287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Albenberg, L.; Esipova, T.V.; Judge, C.P.; Bittinger, K.; Chen, J.; Laughlin, A.; Grunberg, S.; Baldassano, R.N.; Lewis, J.D.; Li, H.; et al. Correlation between intraluminal oxygen gradient and radial partitioning of intestinal microbiota. Gastroenterology 2014, 147, 1055–1063.e8. [Google Scholar] [CrossRef] [Green Version]
  64. Li, H.Y.; Lee, L.A. Sleep-disordered breathing in children. Chang Gung Med. J. 2009, 32, 247–257. [Google Scholar] [PubMed]
  65. Lai, C.C.; Lin, P.W.; Lin, H.C.; Friedman, M.; Chang, H.W.; Salapatas, A.M.; Lin, M.C.; Wang, P.C. Clinical Predictors of Pediatric Obstructive Sleep Apnea Syndrome. Ann. Otol. Rhinol. Laryngol. 2018, 127, 608–613. [Google Scholar] [CrossRef] [PubMed]
  66. Tally, F.P.; Stewart, P.R.; Sutter, V.L.; Rosenblatt, J.E. Oxygen tolerance of fresh clinical anaerobic bacteria. J. Clin. Microbiol. 1975, 1, 161–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jpm 11 00486 g001 550
Figure 1. Flow diagram of the study. Seventy-six children with obstructive sleep apnea were assessed for eligibility. However, three did not meet the inclusion criteria, three met the exclusion criteria, and four withdrew consent. Therefore, a total of 66 children were recruited. The non-healthy-weight group included 33 children, and the healthy-weight group included 33 children. Both groups were matched by the proportion of chronic tonsillitis. All participants underwent adenotonsillectomy. Therefore, 66 participants were included in the primary analysis.
Figure 1. Flow diagram of the study. Seventy-six children with obstructive sleep apnea were assessed for eligibility. However, three did not meet the inclusion criteria, three met the exclusion criteria, and four withdrew consent. Therefore, a total of 66 children were recruited. The non-healthy-weight group included 33 children, and the healthy-weight group included 33 children. Both groups were matched by the proportion of chronic tonsillitis. All participants underwent adenotonsillectomy. Therefore, 66 participants were included in the primary analysis.
Jpm 11 00486 g001
Jpm 11 00486 g002 550
Figure 2. Tonsil microbiome in children with obstructive sleep apnea. (A) The operational taxonomic unit (OTU) tree of the non-healthy-weight subgroup (n = 33) included 9318 OTUs, assessed by Graphical Phylogenetic Analysis. (B) The OTU tree of the healthy-weight subgroup (n = 33) included 9886 OTUs. (C) A Venn diagram demonstrated that both subgroups shared 6539 OTUs; otherwise, the non-healthy-weight subgroup had 2779 deferential OTUs and the healthy-weight subgroup had 3347 deferential OTUs. (D) The α diversities of both subgroups were equal (p = 0.300; Mann–Whitney U test). (E) Furthermore, the β diversity of the non-healthy-weight subgroup was comparable to that of the healthy-weight group (p = 0.147; Mann–Whitney U test). (F) The relative abundances of the top 10 tonsil families in the non-healthy-weight subgroup were similar to those in the healthy-weight subgroup (p = 0.651; analysis of similarity test).
Figure 2. Tonsil microbiome in children with obstructive sleep apnea. (A) The operational taxonomic unit (OTU) tree of the non-healthy-weight subgroup (n = 33) included 9318 OTUs, assessed by Graphical Phylogenetic Analysis. (B) The OTU tree of the healthy-weight subgroup (n = 33) included 9886 OTUs. (C) A Venn diagram demonstrated that both subgroups shared 6539 OTUs; otherwise, the non-healthy-weight subgroup had 2779 deferential OTUs and the healthy-weight subgroup had 3347 deferential OTUs. (D) The α diversities of both subgroups were equal (p = 0.300; Mann–Whitney U test). (E) Furthermore, the β diversity of the non-healthy-weight subgroup was comparable to that of the healthy-weight group (p = 0.147; Mann–Whitney U test). (F) The relative abundances of the top 10 tonsil families in the non-healthy-weight subgroup were similar to those in the healthy-weight subgroup (p = 0.651; analysis of similarity test).
Jpm 11 00486 g002
Jpm 11 00486 g003 550
Figure 3. The top 10 tonsil phyla in both weight status subgroups. The relative abundances of the top 10 tonsil phyla in the non-healthy-weight subgroup (A) and the healthy-weight subgroup (B) (p = 0.651; analysis of similarity test). Chronic tonsillitis was significantly associated with a relative abundance of Acidobacteria in the non-healthy-weight subgroup (r = 0.53, q = 0.015) (C), whereas the oxygen desaturation index (ODI) was significantly associated with Bacteroidetes in the healthy-weight subgroup (r = 0.52, q = 0.020) (D). Abbreviations: AHI, apnea–hypopnea index; BMI, body mass index; RDI, respiratory disturbance index; REM, rapid eye movement; SpO2, pulse oxygen saturation.
Figure 3. The top 10 tonsil phyla in both weight status subgroups. The relative abundances of the top 10 tonsil phyla in the non-healthy-weight subgroup (A) and the healthy-weight subgroup (B) (p = 0.651; analysis of similarity test). Chronic tonsillitis was significantly associated with a relative abundance of Acidobacteria in the non-healthy-weight subgroup (r = 0.53, q = 0.015) (C), whereas the oxygen desaturation index (ODI) was significantly associated with Bacteroidetes in the healthy-weight subgroup (r = 0.52, q = 0.020) (D). Abbreviations: AHI, apnea–hypopnea index; BMI, body mass index; RDI, respiratory disturbance index; REM, rapid eye movement; SpO2, pulse oxygen saturation.
Jpm 11 00486 g003
Jpm 11 00486 g004 550
Figure 4. Comparison of the tonsil microbiome between the mean oxygen saturation (SpO2) ≤ 97% and > 97% subgroups. (A) The α diversity of the mean SpO2 ≤ 97% subgroup was significantly lower than that of the mean SpO2 > 97% subgroup (p = 0.030; Mann–Whitney U test). (B) The β diversity of the mean SpO2 ≤ 97% subgroup was significantly lower than that of the SpO2 > 97% subgroup (p = 0.0005; Mann–Whitney U test). (C) The relative abundances of the top 10 tonsil phyla in the mean SpO2 ≤ 97% subgroup were significantly different from those in the mean SpO2 > 97% subgroup (p = 0.014; analysis of similarity test). (D) Furthermore, the relative abundance of Bacteroidetes in the mean SpO2 ≤ 97% subgroup was significantly higher than that in the mean SpO2 > 97% subgroup (q = 0.026; Mann–Whitney U test).
Figure 4. Comparison of the tonsil microbiome between the mean oxygen saturation (SpO2) ≤ 97% and > 97% subgroups. (A) The α diversity of the mean SpO2 ≤ 97% subgroup was significantly lower than that of the mean SpO2 > 97% subgroup (p = 0.030; Mann–Whitney U test). (B) The β diversity of the mean SpO2 ≤ 97% subgroup was significantly lower than that of the SpO2 > 97% subgroup (p = 0.0005; Mann–Whitney U test). (C) The relative abundances of the top 10 tonsil phyla in the mean SpO2 ≤ 97% subgroup were significantly different from those in the mean SpO2 > 97% subgroup (p = 0.014; analysis of similarity test). (D) Furthermore, the relative abundance of Bacteroidetes in the mean SpO2 ≤ 97% subgroup was significantly higher than that in the mean SpO2 > 97% subgroup (q = 0.026; Mann–Whitney U test).
Jpm 11 00486 g004
Table
Table 1. Patient characteristics, polysomnography variables, and blood pressures of the different weight status subgroups.
Table 1. Patient characteristics, polysomnography variables, and blood pressures of the different weight status subgroups.
VariablesNon-Healthy-Weight SubgroupHealthy-Weight Subgroupp-Value 1
Patient Characteristics
Age (years)7.0 (6.0–8.0)6.0 (5.0–7.5)0.015 *
Male sex, n (%)28 (85)22 (67)0.150
Chronic tonsillitis6 (18)10 (30)0.389
Allergic rhinitis, n (%)22 (67%)25 (76%)0.587
BMI (kg/m2) z-score2.01 (1.46–2.38)−0.36 (−1.16–0.18)<0.001 *
Tonsil size3 (3–4)2 (3–4)0.461
ANR0.73 (0.62–0.83)0.81 (0.72–0.87)0.053
OSA-18 score80 (69–92)81 (70–91)0.928
Polysomnography variables
AHI (events/h)9.6 (5.0–25.2)5.4 (3.9–16.5)0.074
RDI (events/h)12.1 (5.3–27.6)6.1 (4.9–17.9)0.158
ODI (events/h)7.3 (3.6–22.3)3.2 (1.6–9.3)0.006 *
Mean SpO2 (%)97 (96–98)98 (97–98)0.030 *
Minimal SpO2 (%)89 (83–91)91 (88–93)0.022 *
N1 stage13 (6–21)9 (6–13)0.142
N2 stage38 (33–46)41 (36–44)0.807
N3 stage28 (23–30)28 (22–36)0.663
REM stage18 (13–22)21 (16–25)0.221
Blood pressure variables
Systolic BP, mmHg111 (100–121)98 (87–107)0.001 *
Diastolic BP, mmHg67 (61–76)60 (58–68)0.011 *
Systolic BP percentile (%)84 (55–91)48 (25–89)0.018 *
Diastolic BP percentile (%)75 (55–87)68 (50–77)0.174
Note: Data are summarized as median (interquartile range) or n (%) as appropriate. Abbreviations: AHI, apnea–hypopnea index; ANR, adenoid–nasopharyngeal ratio; BMI, body mass index; BP, blood pressure; ODI, oxygen desaturation index; OSA, obstructive sleep apnea; RDI, respiratory disturbance index; REM, rapid eye movement; SpO2, pulse oxygen saturation. 1 Data were compared using the Mann–Whitney U test for continuous variables, and the chi-square test for categorical variables. * Significant differences p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chuang, H.-H.; Hsu, J.-F.; Chuang, L.-P.; Chiu, C.-H.; Huang, Y.-L.; Li, H.-Y.; Chen, N.-H.; Huang, Y.-S.; Chuang, C.-W.; Huang, C.-G.; et al. Different Associations between Tonsil Microbiome, Chronic Tonsillitis, and Intermittent Hypoxemia among Obstructive Sleep Apnea Children of Different Weight Status: A Pilot Case-Control Study. J. Pers. Med. 2021, 11, 486. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11060486

AMA Style

Chuang H-H, Hsu J-F, Chuang L-P, Chiu C-H, Huang Y-L, Li H-Y, Chen N-H, Huang Y-S, Chuang C-W, Huang C-G, et al. Different Associations between Tonsil Microbiome, Chronic Tonsillitis, and Intermittent Hypoxemia among Obstructive Sleep Apnea Children of Different Weight Status: A Pilot Case-Control Study. Journal of Personalized Medicine. 2021; 11(6):486. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11060486

Chicago/Turabian Style

Chuang, Hai-Hua, Jen-Fu Hsu, Li-Pang Chuang, Cheng-Hsun Chiu, Yen-Lin Huang, Hsueh-Yu Li, Ning-Hung Chen, Yu-Shu Huang, Chun-Wei Chuang, Chung-Guei Huang, and et al. 2021. "Different Associations between Tonsil Microbiome, Chronic Tonsillitis, and Intermittent Hypoxemia among Obstructive Sleep Apnea Children of Different Weight Status: A Pilot Case-Control Study" Journal of Personalized Medicine 11, no. 6: 486. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11060486

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop