Gut microbiota in obstructive sleep apnea–hypopnea syndrome: disease-related dysbiosis and metabolic comorbidities

Obstructive sleep apnea–hypopnea syndrome (OSAHS) is characterized by repeated collapse episodes in the upper airway during sleep, resulting in intermittent hypoxia (IH), together with symptoms of sleep fragmentation (SF) and daytime sleepiness. OSAHS is a systemic disorder; it is highly prevalent in Asia and is estimated to occur in 3.7–97.3% of the population [1]. Further, it is one of the most common sleep disorders and is associated with a growing list of OSAHS-induced end-organ morbidities requiring expertise in respiratory medicine, cardiology, neurology, and endocrinology. OSAHS-mediated oxidative stress and endothelial dysfunction are associated with the pathogenesis of atherosclerosis and cardiovascular morbidity [2]. In this context, IH/reoxygenation and SF are the key features of OSAHS-enhanced inflammation and the oxidative stress cascade [3,4]. By activating multiple systemic inflammatory mediators found in immune cells, OSAHS leads to an enhanced, sustained proinflammatory state that contributes to multiorgan morbidities [2], thereby suggesting a linkage with the host’s gut microbiota.

The gut microbiome plays a crucial role in modulating the risk for several chronic diseases, including obesity, metabolic abnormalities, cardiovascular disease, and inflammatory bowel disease (IBD) [5]. Delicate balance in the composition of the gut microbiome is the key to maintain intestinal immunity and whole-body homeostasis [6]. Gram-positive Firmicutes (F) and Gram-negative Bacteroidetes (B) are predominant in the gut microbiome [5]. Imbalance in the Firmicutes/Bacteroidetes (F/B) ratio has been closely related to inflammatory and immune function disorders, obesity, and metabolic disease [7].

Experiments conducted in an animal model of IH mimicking OSAHS and in an animal model of IH + high-fat diet (HFD) have demonstrated alterations in the diversity of gut microbiota. Previous studies have shown an increase in the F/B ratio, which is in turn associated with an increase in the relative abundance of lactate-producing bacteria and a decrease in the relative abundance of short-chain fatty acid (SCFA)-producing bacteria [3,6,8,9]. In addition, chronic SF alters feeding behaviors, thereby promoting obesity and metabolic abnormalities. The host gut microbiome changes, increasing intestinal permeability as well as systemic and adipose tissue inflammation. These changes are accompanied by insulin resistance [4,10]. Although these findings relate to rodent paradigms, we aimed to determine whether these characteristics are also present in humans with OSAHS-related metabolic comorbidities.

Enterotype is another means to investigate the gut microbiota. Major enterotype bacteria may be subdivided into three types: Bacteroides (enterotype 1), Ruminococcus (enterotype 2), and Prevotella (enterotype 3) [11]. Dietary habits affect enterotypes, thereby emphasizing the key role of diet in the composition of the gut microbiome [12]. The balance of the three enterotypes is maintained by various enzymes and is not associated with ethnicity, gender, age, or body mass index (BMI) [11].

Alteration in the gut microbiota has been reported in IH and SF mimicking OSAHS rodent paradigms. However, only few studies have directly investigated the gut microbiota in patients with OSAHS, who often present with metabolic comorbidities. Hence, here, using taxa and enterotype data, we aimed to investigate gut microbial dysbiosis in individuals with varying scores on the apnea–hypopnea index (AHI). This classification differs from IH modeling paradigm, which considers only two groups. In patients with OSAHS, homocysteine (HCY) levels are associated with an increased risk of cardiovascular events and metabolic abnormalities [13]. Therefore, we assessed the correlation between the gut microbiota and HCY levels in relation to OSAHS-related metabolic comorbidities.

We recruited 113 subjects and examined them by a full night of polysomnography (PSG; SOMNOscreen™ plus PSG⁺, SOMNOmedics GmbH, Randersacker, Germany), conducted by technologists in a sleep laboratory from 10 p.m. to 8 a.m. at the Department of Pulmonary and Critical Care Medicine. The Institutional Review Board of the Second Affiliated Hospital of Fujian Medical University approved the present study (IRB No. 2017-78). We collected fasting blood and fecal samples the following morning.

All subjects underwent PSG (performed with a computerized polysomnographic system), which included electrocardiography, electroencephalography, electromyography, and electrooculography. After one night of the examination, AHI was calculated as the total number of episodes of apnea (continuous cessation of airflow for at least 10 s) and hypopnea (reduction in airflow for ≥10 s with oxygen desaturation of ≥4%) divided by the total duration of sleep events, according to the diagnostic criteria of the American Academy of Sleep Medicine. AHI ≤ 15 events/h was defined as non-OSAHS (control group), 5 < AHI ≤ 15 as mild OSAHS (Group 1), 15 < AHI ≤ 30 as moderate OSAHS (Group 2), and AHI ≥ 30 as severe OSAHS (Group 3).

We determined the level of fasting HCY serum with the Automatic Biochemical Analyzer (TBA-120 FR, Toshiba, Japan) and the HCY assay kit (Yong He Sun Biotech. Ltd., Hunan, China), utilizing the enzymatic cycling method.

Interleukin (IL)-6 (IL-6) and tumor necrosis factor α (TNF-α) were assayed by BD Human Enhanced Sensitivity Cytometric Bead Array Kit (BD Biosciences, New Jersey, U.S.A.) as described previously [14]. The standard coefficient of determination (r²) was greater than 0.995 (detail see in Supplementary Figure S1).

All fresh fecal samples were collected and stored in Microbiome Test Kit (G-BIO Biotech, Inc., Hangzhou, China). Magnetic bead isolation was done to extract genomic DNA using a TIANamp stool DNA kit (TIANGEN Biotech Co., Ltd., Beijing, China), according to the manufacturer’s instructions. The concentration of extracted DNA was determined by a Nanodrop ND-1000 spectrophotometer (Thermo Electron Corporation, U.S.A.), and DNA quality was confirmed using 1.0% agarose gel electrophoresis with 0.5 mg/ml ethidium bromide.

Isolated fecal DNA was used as a template to amplify the V3 and V4 hypervariable regions of the bacterial 16S rRNA gene. The V3 and V4 regions were PCR-amplified (forward primer, 5′-ACTCCTACGGGAGGCAGCAG-3′; reverse primer, 5′-GGACTACHVGGGTWTCTAAT-3′) [15]. The 16S target-specific sequence contained adaptor sequences permitting uniform amplification of a highly complex library ready for downstream next-generation sequencing with Illumina MiSeq (Illumina, U.S.A.). Negative DNA extraction controls (lysis buffer and kit reagents only) were amplified and sequenced as contamination controls. The amplicons were normalized, pooled, and sequenced on the Illumina MiSeq platform, using a V3 reagent kit with 2 × 300 cycles per sample, and with imported and prepared routine data (samsheet) run in the MiSeq sequence program. After sequencing, Q30 scores were ≥70%, the percentage of clusters passing filter (i.e. cluster PF) was ≥80%, and there were at least 30000 clean tags. Finally, image analysis and base calling were conducted with the MiSeq Control Software.

Based on the Quantitative Insights into Microbial Ecology bio-informatic pipeline for performing taxonomy assignment by the operational taxonomic unit (OTU) method, we used 113 sequence data to analyze the fecal microbiota taxa. We predicted the bacterial metabolic function using the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States bioinformatics software package and the Kyoto Encyclopedia of Genes and Genomes (KEGG).

We discriminated OSAHS patients from controls with an area under the receiver operating characteristic curve (ROC-AUC) with a patient disease index (PDI) cut-off, wherein ROC-AUC and PDI were calculated using a previous study [16].

Data were expressed as the mean ± standard deviation (S.D.). We analyzed differences in gut microbiota using the Kruskal–Wallis test, as appropriate, and performed principal coordinate analysis (PCoA) on the basis of the Bray–Curtis distance function, using R statistics. We performed other analyses using statistically with SPSS version 19.0 (SPSS Inc., Chicago, IL, U.S.A.) by one-way ANOVA, followed by Scheffe post hoc analyses, we considered a two-sided P<0.05 to be statistically significant. We evaluated the correlation coefficients between gut microbiota and HCY level using Spearman correlation.

After PSG assessment, we enrolled 93 patients with OSAHS [Group 1 (n=40), Group 2 (n=23), and Group 3 (n=30)] and 20 controls (Table 1). Body weight (P=0.006), BMI (P=0.010), and hip circumference (P=0.025) were significantly higher in the Group 3 than in the control group. Waist circumference [P<0.001 (Control vs. Group 1), P=0.005 (Control vs. Group 2), P<0.001 (Control vs. Group 3)], and waist-to-hip ratio [P<0.001 (Control vs. Group 1), P<0.001 (Control vs. Group 2), P<0.001 (Control vs. Group 3)] were significantly lower among the control group than in the OSAHS groups (Table 1).

Our PSG data indicated that AHI [P=0.014 (Control vs. Group 1), P<0.001 (Control vs. Group 2), P<0.001 (Control vs. Group 3)] and hypopnea index [P=0.036 (Control vs. Group 1), P<0.001 (Control vs. Group 2), P<0.001 (Control vs. Group 3)] were significantly higher in the OSAHS groups than in the control group. The lowest oxygen saturation in the OSAHS groups was lower than that in the control group [P=0.022 (Control vs. Group 1), P<0.001 (Control vs. Group 2), P<0.001 (Control vs. Group 3)]. Mean oxygen saturation was significantly lower in the Group 3 than in the control group (P<0.001; Table 1).

We examined the mean community diversity indices (Supplementary Figure S2) and the structure of the gut microbiome (Supplementary Figure S3). At the phylum level, we found no significant differences in relative abundance, including the F/B ratio (Figure 1A). The relative abundances of the following genera significantly differed among groups: Faecalibacterium (P=0.044), Megamonas (P=0.046), Ruminococcaceae (P=0.048), Clostridiales (P=0.021), Alistipes (P=0.024), Bifidobacterium (P=0.037), Dialister (P=0.007), Oscillibacter (P=0.008), Erysipelotrichaceae (P=0.025), Anaerotruncus (P=0.043), Cronobacter (P=0.006), Aquabacterium (P=0.033), Gammaproteobacteria (P=0.038), Synergistaceae (P=0.032), Betaproteobacteria (P=0.020), and Micrococcus (P=0.007) (Figure 1B; for details, see Supplementary information). We found no significant differences between the control and OSAHS groups in terms of the relative abundances of acetate-, butyrate-, propionate-, or lactate-producing bacteria (Figure 1C).

IL-6 of the Group 3 was significantly higher than the control group (P=0.006). However, there was not significantly different in TNF-α among the control group and OSAHS groups (Figure 1D).

According to KEGG, there were significant differences in the fecal microbiome among these study groups, covering 14 pathways involved in arginine and proline metabolism [P=0.016 (Control vs. Group 3)]; alanine, aspartate, and glutamate metabolism [P=0.049 (Control vs. Group 3)]; galactose metabolism [P=0.002 (Control vs. Group 1), P=0.003 (Control vs. Group 3)]; phosphotransferase system [P=0.037 (Control vs. Group 1), P=0.028 (Control vs. Group 3)]; propanoate metabolism [P=0.027 (Group 1 vs. Group 2)]; C5-branched dibasic acid metabolism [P=0.029 (Control vs. Group 1)]; insulin signaling pathway [P=0.012 (Control vs. Group 2)]; Vibrio cholerae pathogenic cycle [P=0.007 (Control vs. Group 1), P=0.006 (Control vs. Group 3)]; tropane, piperidine, and pyridine alkaloid iosynthesis [P=0.007 (Control vs. Group 3), P=0.031 (Group 2 vs. Group 3)]; primary immunodeficiency [P=0.031 (Group 2 vs. Group 3)]; linoleic acid metabolism [P=0.027 (Control vs. Group 3)]; prostate cancer [P=0.015 (Control vs. Group 1)]; flavone and flavonol biosynthesis [P<0.001 (Control vs. Group 1), P=0.027 (Control vs. Group 3)]; and bladder cancer [P=0.045 (Control vs. Group 2)] (Figure 2A).

Neither lactate- nor SCFA-producing bacteria had any correlation with the HCY level, apart from Lactobacillus, which is lactate- and acetate-producing bacterium, with a positive correlation with the HCY level (P=0.041; Table 2).

To select predictive features, we discriminated OSAHS patients from controls with ROC-AUC 0.789 using different taxa (Figure 2B). The best predictive point was the false positive rate of 0.14 and the true positive rate of 0.60, with a PDI cut-off of −0.00738 (P=0.00003; data not shown).

In Figure 3A, we present the distribution of fecal taxa in the control group and OSAHS groups. For enterotype 1, the relative abundances of the genera Lachnospiraceae, Erysipelotrichaceae, Bifidobacterium, Dialister, Burkholderiales, Oscillibacter, Acidaminococcus, Phascolarctobacterium, and Cronobacter were significantly different among four groups (Figure 3B; detail see in Supplementary information). For enterotype 2, the relative abundance of the genera Sutterella, Raoultella, Collinsella, Bacteroidetes, Methylobacterium, Prevotellaceae, and Victivallis were significant different among four groups (Figure 3C; for detail see Supplementary information). For enterotype 3, Howardella, Synergistes, and Brevibacterium were higher in the control group than in the OSAHS groups (Figure 3D; for detail see Supplementary information). After conducting further stratification analysis utilizing the enterotypes, the estimated OSAHS odds ratios were much higher in Ruminococcus enterotype (3.65) and Prevotella enterotype (2.15) for OSAHS risky (Table 3).

We compared the differences in the gut microbiome between the control and OSAHS groups. The number of candidate taxa associated with gut microbial dysbiosis in OSAHS groups were 2 at the class level, 3 at the order level, 10 at the family level (details in Supplementary Figure S4), and 16 at the genera level. Decrease in the relative abundance of SCFA-producing bacteria and increased levels of proinflammatory cytokines were also observed. These observations were supported by functional analyses of the microbiota. Lactobacillus levels correlated with HCY levels. Our validation cohort showed an ROC-AUC value of 0.789 for these samples, confirming that gut microbiota information can be used to identify patients. The Ruminococcus enterotype was associated with an increased risk for OSAHS. We contend that these findings illustrate the pathophysiological role played by changes in the gut microbiota in OSAHS. These changes might be associated with OSAHS-related metabolic comorbidities.

IH and SF may play pathophysiological roles in changes to the gut microbiota associated with OSAHS. Gut microbial dysbiosis has been demonstrated in several OSAHS-mimicking animal paradigms [3,4,8,9]. However, OSAHS is a systemic and complex disorder. IH, IH + HFD, or SF alone cannot be used to diagnose a patient with OSAHS. Contrary to the results reported previously for animal models [3,10], our results did not show significant differences among Chinese OSAHS patients in terms of Chao richness, Shannon diversity, and Simpson diversity (Supplementary Figure S2). Observations regarding alterations in the gut microbiome in the IH rodent model suggest that although the gut epithelium is remarkably resistant to hypoxia, the absorptive and barrier functions that regulate the intestinal epithelium are sensitive to the level of oxygen inside the gut. By increasing permeability and bacterial translocation and decreasing tight junction integrity, IH/reoxygenation may directly impair cellular function [17,18]. The increasing abundance of anaerobic bacteria noted in the IH model is due to the emergence of anoxic environment in the intestine, which is beneficial for the growth of obligate anaerobic bacteria and endogenous lipopolysaccharide (LPS) production by Gram-negative bacteria [3,8]. Severe circulating endotoxemia was maintained throughout prolonged normoxic recovery after IH exposure in mice [8], suggesting that IH alters the gut microbiome and endotoxin levels.

In addition, SF-induced sleep perturbations led to profound alterations in the gut microbiome, such as increases in the circulating levels of IL-6, neutrophil gelatinase-associated lipocalin, and LPS-binding protein, particularly in patients with obesity [4,19]. These observations suggest the induction of inflammatory processes, perhaps owing to the leakage of microbial metabolites into the circulation, which is similar to our findings in patients with severe OSAHS. SF-induced sleep disturbances also promote the innate immune response, thereby leading to low-grade systemic and adipose tissue inflammation and ultimately resulting in metabolic homeostasis disruption [4]. SF treatment altered the feeding behaviors, eventually promoting obesity and metabolic abnormalities. These findings indicate that SF changes the host’s gut microbiota, increases intestinal permeability, and promotes inflammation in adipose tissue and throughout the body. These changes are accompanied by insulin resistance [4]. Taken together, IH and SF may induce alterations in the intestinal epithelial barrier and increase intestinal permeability, leading to local and systemic inflammatory responses and multi-metabolic abnormalities, such as obesity, hypertension, and diabetes [20,21].

The prevalence of OSAHS in patients with obesity/metabolic syndrome is estimated at 45–60% [22,23]. In our study, patients with OSAHS had a higher waist circumference and waist-to-hip ratio than patients in the control group; in particular, patients with severe OSAHS had higher body weight, BMI, and IL-6 levels. We also observed an increased risk for metabolic abnormalities with increasing age and AHI in males with OSAHS, as reported previously [22]. OSAHS may alter the expression of intestinal epithelial barrier markers, thereby increasing intestinal permeability. Persistently increased permeability may be an important contributor to the development of metabolic complications [21]. Gut permeability is further increased during the development of obesity and metabolic disorders [19]. Obesity involves low-grade inflammation and elevations in the F/B ratio [5]. Although we noted that patients with OSAHS have a higher BMI and greater risk for metabolic abnormalities, we cannot draw a firm conclusion from our data. However, functional analysis of fecal samples collected from patients with OSAHS indicated effects on pathways involved in the down-regulation of amino acid metabolism (arginine and proline metabolism; alanine, aspartate, and glutamate metabolism); the elevation of galactose and propanoate metabolism; the down-regulation of C5-branched dibasic acid metabolism, linoleic acid metabolism, and flavone and flavonol biosynthesis; and the down-regulation of the insulin signaling pathway. These factors are connected to the energy metabolism that promotes obesity, metabolic abnormalities, and inflammation [4,10], indirectly corresponding to the primary immune disruption in our OSAHS groups, as revealed by functional analysis.

The disruption of gut permeability is initially a cascading factor; however, SCFA may regulate gut integrity [24,25]. Our OSAHS groups had lower fecal counts of Faecalibacterium, Megasphaera, Ruminococcaceae, Dialister, and Oscillibacter, which are lactate- or SCFA-producing bacteria [26]. Microbial fermentation of dietary fiber generates SCFA. In some bacterial species, lactate can be converted into butyrate. SCFA-producing bacteria also produce acetate, propionate, and butyrate. A decrease in SCFA production results in intestinal barrier dysfunction [27,28]. SCFA has a direct anti-inflammatory effect on the gut, contributes to mucin synthesis, decreases bacterial translocation, maintains gut integrity, and mitigates inflammation in the intestine [8,24,25]. Moreover, SCFAs participate in immunity, adipogenesis, insulin sensitivity, and oxidative stress [26]. Faecalibacterium, Oscillibacter, and Howardella genera were highly enriched in controls. The levels of Faecalibacterium and Oscillibacter, which contribute to anti-inflammatory activity within the gut [29], were negatively correlated with IBD [25]. We speculate that the relative abundance of anti-inflammation-associated bacteria among the patients with OSAHS at least partly mitigated the degree of inflammation. We also detected an overgrowth of proinflammatory enteric pathogens (the class Gammaproteobacteria and the family Enterobacteriaceae; Supplementary information), such as common Gram-negative bacteria that promote LPS production and thereby increase inflammation [30].

In addition, Lactobacillus and Escherichia are positively correlated with IBD [25]. Altered levels of microbial butyrate and lactate may contribute to OSAHS-induced hypertension [6,9]. Lactate bacteria in the human gastrointestinal tract, including Bifidobacteria, Lactobacilli, Streptococci, and Enterococci, predominantly produce l- and/or d-lactate, which can also be produced by strict anaerobes, such as Eubacterium spp., Selenomonas, Megasphaera, and Veillonella that can convert lactate into acetate and propionate [27]. l-lactate may enter the gut from host tissues. In short-bowel syndrome, d-lactate accumulation may lead to serious outcomes, including neurotoxicity and cardiac arrhythmia. Under normal conditions, lactate is seldom detected in human feces or gut content as a major fermentation product of mixed anaerobic communities [27]. However, plasma lactate concentrations are positively associated with increased blood pressure [31]. Our data reveal that HCY levels are positively correlated with the levels of Lactobacillus, which is a lactate bacterium. In patients with OSAHS, HCY levels are associated with an increased risk for metabolic abnormalities and hypertension [13]. Lactobacillus plantarum has been demonstrated to reduce HCY levels [32] and to reduce the risk for cardiovascular disease [33]. This positive effect may be responsible for the production of SCFA through the bacterial fermentation of fiber [33]. However, the results of our functional analyses of the microbiome did not reflect information about the depletion of vitamins B₆ and B₁₂, folates, or methionine in relation to HCY metabolism. Instead, arginine participates in the synthesis of nitric oxide (NO), which declines with increased blood pressure. Arginine levels decreased in our patients with OSAHS, 40–67% of whom suffer from hypertension. Arginine, methionine, HCY, and vitamins influence NO synthesis, eventually affecting the endothelial function that controls vascular homeostasis [34]. Furthermore, diet seemed to play a key role in the synthesis of HCY and NO as well as the composition of the gut microbiota. We suggest that Lactobacillus plays important roles in OSAHS-related metabolic comorbidities; however, this hypothesis requires additional evidence.

On the other side, enterotype analysis has been proposed as a useful method to understand human gut microbial communities, including Bacteroides, Ruminococcus, and Prevotella enterotypes, irrespective of ethnicity, gender, age, or BMI [11]. The IH rodent model may be classified as enterotypes Bacteroides and Prevotella [3]. Using this model, researchers found an increase in the abundance of Prevotella and Desulfovibrio; however, this trend was not observed in our human data described above. Two types of mucus-degrading bacteria deplete the mucosa on the epithelial layer of the gut when exposed to IH, resulting in altered intestinal permeability [3]. However, we conducted stratification analysis using various enterotypes to reveal that the Ruminococcus enterotype posed the highest risk for OSAHS. Although the Prevotella enterotype belongs to the phyla Bacteroidetes, the accompanying endotoxins render the host prone to low-grade inflammation, triggering the progression of OSAHS. The Prevotella enterotype is associated with diets high in carbohydrates (fiber) [12]. Therefore, we speculated that could promote microbial fermentation and produce SCFA to protect the host.

In contrast, the Bacteroides enterotype is associated with Western-style diets, including the consumption of high amounts of protein and fat [12]. According to our findings, this enterotype comprises the putative anti-inflammation-associated bacteria Lachnospiraceae, Phascolarctobacterium, Dialister, Oscillibacter, Acidaminococcus, and Bifidobacterium [26,35,36] to a greater degree than the pathobiont Burkholderiales. Their gut microbiota protects the intestinal mucosa, maintains intestinal permeability, and attenuates the inflammatory response, ultimately reducing the risk for OSAHS. Patients with autism spectrum disorder with the Ruminococcus enterotype have high levels of the enriched bacteria Sutterella spp. and Ruminococcus spp. [37]. Patients with OSAHS also have increased putative pathobionts (e.g. Raoultella spp., Methylobacterium), which can trigger infection [38,39] and decrease the relative abundance of Collinsella in smokers with IBD [40], as observed in our patients. However, the non-OSAHS patients in this enterotype may maintain constant health, leading to the speculation that Prevotellaceae, Bacteroidetes, and Victivallis spp. are putative SCFA-producing bacteria [41,42] and occur to a lesser degree in patients with OSAHS than in normal controls.

There are certain limitations to our study. The first limitation is the small sample size, particularly for the Ruminococcus enterotype; the present data must therefore be conservatively interpreted. Second, BMI should not affect the enterotype theory [11]. A future prospective study should include controls and OSAHS groups with comparable BMI. Third, although the gut bacterial profiles are influenced by dietary habits [5,12], we did not collect or analyze data regarding the participants’ dietary habits; this will need to be corrected in the future because long-term dietary differences could cause differences in a patient’s microbiome and confer different health risks. Fourth, we did not comprehensively analyze lactate and endotoxins in blood; thus, even greater caution is required when interpreting alterations in the relative abundance of microbiota taxa with regard to their effects on the pathogenesis of OSAHS.

IH and SF may cause gut microbial dysbiosis. Changes in the gut microbiome are associated with reduced SCFA production and increased pathogen levels, and those changes induce alterations in the levels of intestinal epithelial barrier markers and increase intestinal permeability, leading to local and systemic inflammatory responses and metabolic comorbidities. Our findings show that changes in the gut microbiota have a pathophysiological role in OSAHS and suggest that changes in the gut microbiome may be associated with the pathophysiology of metabolic comorbidities in patients with OSAHS.

• Gut microbial dysbiosis is found in varying degrees in patients with OSAHS. We performed enterotypes stratification analysis to identify a panel of gut microbiome biomarkers that could be used as a non-invasive test with which to accurately diagnose OSAHS.

• Decreases in the relative abundance of SCFA-producing bacteria and elevated pathogens, accompanied by elevated levels of proinflammatory cytokines, are associated with the pathophysiology of OSAHS-related metabolic comorbidities.

• Lactobacillus levels were correlated with HCY levels. Elucidating the interactions between microbiomes and body homeostasis may shed some new insight into the pathophysiology of metabolic comorbidities in OSAHS.

We thank all the participants and their families who took part in the present study. We thank the assistance provided by Huan Wu (G-BIO Biotech, Inc., Hangzhou, China) for performing the bioinformatics analysis.

The authors declare that they have no financial or personal relationships with others who might inappropriately influence the results or interpretation in this manuscript.

We thank the Fujian Provincial Health and Family Planning Commission, China [grant number 2018-CX-36]; Fujian Province Science and Technology Project, China [grant numbers 2018J01290 and 2015J01445]; and Quanzhou Science and Technology Project, China [grant numbers 2018Z107, 2018N003S, 2017Z016, 2013Z55]; and the Young people training project from Fujian Province Health System [grant number 2014JQNJC20]. This article was also subsidized by Strait Exchanges Postdoctoral Funding Schemes of Fujian Province, China, and academic funding of the Second Affiliated Hospital of Fujian Medical University [grant number BSH001].

Conception and design: C.-Y.K., H.-P.Z., and Y.-M.Z. Acquisition of data: C.-Y.K., Q.-Q.L., H.-Z.S., J.-M.F., J.-H.Y., A.-K.H., and Y.-Q.L. Analysis and interpretation of data: C.-Y.K., D.C., H.-Z.S., H.-P.Z., and Y.-M.Z. Drafting/revising of the article: C.-Y.K., D.C., H.-P.Z., and Y.-M.Z.

AHI

apnea–hypopnea index

BMI

body mass index

F/B

Firmicutes/Bacteroidetes

HCY

homocysteine

HFD

high-fat diet

IBD

inflammatory bowel disease

IH

intermittent hypoxia

IL-6

interleukin-6

KEGG

Kyoto Encyclopedia of Genes and Genomes

LPS

lipopolysaccharide

NO

nitric oxide

OSAHS

obstructive sleep apnea–hypopnea syndrome

PSG

polysomnography

ROC-AUC

area under the receiver operating characteristic curve

SCFA

short-chain fatty acid

SF

sleep fragmentation

TNF-α

tumor necrosis factor α