key: cord-0682104-oxazgpq9 authors: Chen, Xuehui; Chen, Yanlong; Feng, Mengqi; Huang, Xin; Li, Changtao; Han, Fang; Zhang, Qian; Gao, Xuemei title: Altered Salivary Microbiota in Patients with Obstructive Sleep Apnea Comorbid Hypertension date: 2022-04-08 journal: Nat Sci Sleep DOI: 10.2147/nss.s347630 sha: c7e16bcdf6fd7d858761067fc2cc4fc0be706ec1 doc_id: 682104 cord_uid: oxazgpq9 PURPOSE: Microorganisms contribute to the pathogenesis of obstructive sleep apnea (OSA)-associated hypertension (HTN), while more studies focus on intestinal microbiome. However, the relationship between oral microbiota and OSA-associated HTN has yet to be elucidated. This study aimed to identify differences in salivary microbiota between patients with OSA comorbid HTN compared with OSA patients, and furthermore evaluate the relationship between oral microbiome changes and increased blood pressure in patients with OSA. PATIENTS AND METHODS: This study collected salivary samples from 103 participants, including 27 healthy controls, 27 patients with OSA, 23 patients with HTN, and 26 patients with OSA comorbid HTN, to explore alterations of the oral microbiome using 16S rRNA gene V3-V4 high-throughput sequencing. And ultra-high-performance liquid chromatography was used for metabolomic analysis. RESULTS: Alpha- and beta-diversity analyses revealed a substantial difference in community structure and diversity in patients with OSA comorbid HTN compared with patients with OSA or HTN. The relative abundance of the genus Actinomyces was significantly decreased in patients with HTN compared with healthy controls, and those with OSA concomitant HTN compared with the patients in OSA, but was not significantly different between patients with OSA and healthy controls. Linear discriminant analysis effect size and variance analysis also indicated that the genera Haemophilus, Neisseria, and Lautropia were enriched in HTN. In addition, Oribacterium was an unique taxa in the OSA comorbid HTN group compared with the control group. Metabolomic analysis of saliva identified compounds associated with cardiovascular disease in patients with OSA comorbid HTN.2-hydroxyadenine, was significantly increased in the group of patients with OSA compared with controls, and L-carnitine was significantly decreased in patients with OSA comorbid HTN compared with OSA patients. CONCLUSION: This study highlighted noninvasive biomarkers for patients with OSA comorbid HTN. As the first study to find alterations of the salivary microbiome in patients with OSA comorbid HTN, it may provide a theoretical foundation for clinical diagnosis and treatment of this condition. Systemic hypertension, a very common chronic disease, is a leading risk factor for stroke, heart failure, atrial fibrillation, chronic kidney disease, and other life-threatening diseases. 1 It also coexists with obstructive sleep apnea (OSA), which is characterized by intermittent and repeated collapse of the upper airway, resulting in intermittent hypoxia (IH), hypercapnia, daytime sleepiness, and other consequences such as hypertension (HTN). 2 Approximately 50% of adult patients requirements of each group. Exclusion criteria for the study were: a) inappropriate disease history or serious disease(s)/conditions of all systems other than HTN-for example, serious cardiovascular disease (myocardial infarction, etc.), respiratory disease (chronic obstructive pulmonary disease, COVID-19, etc.), other sleep conditions (insomnia, etc.), serious neuropsychiatric diseases, serious endocrine system conditions (diabetes, etc.) or the presence of a malignant tumor; b) inappropriate oral signs, such as severe periodontitis, Sjogren's syndrome, oral ulcer, presence of a dental prosthesis in the mouth, or other oral diseases; c) inappropriate lifestyle habits and medical history, including vegan diet, periodontal treatment in the past 6 months, and use of antibiotics within the last 3 months. OSA was defined by standard American Academy of Sleep Medicine criteria: AHI≥5 was defined as OSA. 28 HTN was diagnosed by ambulatory blood pressure monitoring, defined as 24-hour systolic blood pressure (SBP) ≥130 mmHg or diastolic blood pressure (DBP) ≥80 mmHg or use of medications for HTN, which is equal to clinic SBP ≥140 mmHg or DBP ≥90 mmHg. The measurement of blood pressure accorded with guidelines. Participants sat in a chair relax for at least 5 min, avoided caffeine, exercise and did not smoke for at least 30 min. 29 Four groups were included in the study. There were 27 patients in OSA group who went through full-night hospitalization of polysomnography (PSG) in Sleep Center of Peking University People's Hospital. The blood pressure of these participants remained at normal levels. The complication group comprised 26 patients who had both OSA and HTN. The control group included 27 patients that did not exhibit positive findings by PSG or blood pressure monitoring. Finally, the HTN group contained 23 patients who came from the Cardiology Department of Beijing Haidian Hospital having been diagnosed by 24-h ambulatory blood pressure monitoring; subjects were recruited to this group when their clinic blood pressure was ≥140/90 mmHg and measured their real-time blood pressure and recorded the data as their blood pressure before saliva collection. Home sleep test indicated that they did not meet the diagnostic conditions of OSA. All subjects did not eat, drink, or brush their teeth for at least 2 h prior to saliva collection at 7:00-9:00 a.m., then rinsed their mouth with pure water, and rested for at least 10 min before saliva collection. A minimum of 3 mL unstimulated whole saliva samples were collected into sterile centrifugal tubes, then immediately placed on ice and processed within 2 h. Saliva samples were centrifuged at 12,000 g for 30 min at 4°C, and the resulting precipitate and supernatant were stored separately at −80°C. Genomic DNA was extracted using a QIAamp DNA mini kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The purity and concentration of extracted DNA was tested by a NanoDrop 8000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), and the molecular size of the DNA was estimated by 1% agarose gel electrophoresis. DNA samples were stored at −80°C until further use. Extracted DNA was used as template to amplify the V3-V4 variable region of bacterial 16S rRNA genes. The forward primer was 5′-ACTCCTACGGGAGGCAGCA-3′ and the reverse primer was 5′-GGACTACHVGGGTWTCTAAT-3′. The thermal cycle consisted of an initial denaturation at 98°C for 5 min, followed by 25 cycles of denaturation at 98°C for 30s, annealing at 53°C for 30s, and extension at 72°C for 45s, then a final extension for 5 min at 72°C. After purification using Vazyme VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China) and quantification by the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA), amplicons were pooled in equal amounts, and pairedend (2×250 bp) sequencing was performed using an Illumina Novaseq TM 6000 platform (Illumina, San Diego, CA, USA) at Shanghai Personal Biotechnology Co., Ltd (Shanghai, China). After trimming off the barcodes and primers, the raw sequence data were filtered, denoised, merged, and chimeras removed using the DADA2 plugin. Amplicon sequence variants (ASVs) were classified taxonomically according to the Human Oral Microbiome Database (HOMD) (V15.2). Metabolomic analysis was performed on random samples from each group, including six samples from the control group, six from the HTN group, six from the OSA group, and seven samples from the complication group. The supernatant of the salivary samples was thawed at 4°C and mixed in equal portions of 100 μL with 400 μL cold methanol/acetonitrile (1:1, V/V) to remove protein. The mixture was centrifuged for 15 min (14,000 g, 4°C) and the resulting supernatant was dried in a vacuum centrifuge. For liquid chromatography-mass spectrometry (LC-MS) analysis, the samples were redissolved in 100 μL acetonitrile/water (1:1, V/V) solvent. Analyses were performed using an ultra-high-performance liquid chromatography (UHPLC) system (1290 Infinity LC, Agilent Technologies) coupled to a quadrupole time-of-flight mass spectrometer (AB Sciex TripleTOF 6600). The samples were separated on the Agilent 1290 Infinity LC by hydrophilic interaction liquid chromatography (HILIC) column with a column temperature of 25°C, flow rate of 0.5 mL/min, injection volume of 2 μL, and mobile phase composition of 25 mM ammonium acetate, 25 mM ammonia, water, and acetonitrile. QC samples were inserted into the sample queue to monitor and evaluate the stability of the system and the reliability of the experimental data. One-way analysis of variance (ANOVA) in SPSS 22.0 (IBM, Armonk, NY, USA) was used to evaluate the differences of clinical indicators among the four experimental groups. Sequence data analyses were mainly performed using QIIME2 (2019.4) 30 and R packages (v3.2.0). Multivariate linear regression was performed to test associations between clinical covariates (age, BMI) and the salivary microbiome, with a q-value <0.10 used as significance cut-off. (MaAsLin, Galaxy Version 1.0.1). 16 Alpha diversity was measured with the indices of Chao1, Shannon, and Richness, respectively. Beta diversity analysis was performed to investigate the structural variation of microbial communities across samples using Bray-Curtis metrics and was visualized using principal coordinates analysis (PCoA). Comparisons of the similarities between different groups were conducted using Analysis of similarities (Anosim) based on Bray-Curtis metrics. Kruskal-Wallis tests in SPSS 22.0 were used to analyze the different taxa at the genus level, with p<0.05 considered statistically significant. Linear discriminant analysis (LDA) effect size (LEfSe) analysis was used to determine the organisms most likely to explain differences among the groups. A nested crossvalidation based on the random forest model was conducted to illustrate the diagnostic value of the salivary microbiome, using a receiver operating characteristic (ROC) curve to evaluate the diagnostic value of the salivary microbiome. Microbial functions were predicted by PICRUSt2 with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and STAMP software was used to identify the difference of pathways by White nonparametric t-test, with p<0.05 considered statistically significant. For metabolomic analyses, after normalizing the total peak intensity to processed data, the R package was used for multivariate data analysis, including Pareto-scaled principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA). When VIP value >1, metabolites were further tested at the univariate level by Student's t-test to measure the significance of each metabolite; p<0.05 was considered statistically significant. A total of 103 subjects were recruited to the study. There were 27 patients in the control group, 27 in the OSA group, 23 in the HTN group, and 26 in the complication group (OSA+HTN). General information for the different groups is shown in Table 1 . No significant differences were observed in terms of age, body mass index (BMI), and AHI between the OSA group and the complication group, However, SBP and DBP were significantly increased in the complication group compared with the OSA group (p<0.001 and p=0.002 for SBP and DBP, respectively). In addition, compared with the HTN group, the complication group had a considerably higher BMI and AHI (p=0.031 and p=0.001, respectively), despite no significant differences in age, SBP, or DBP. In consideration of age and BMI, were significantly differences among groups, multivariate linear regression analysis revealed that age and BMI did not affect the difference in salivary microbiomes among groups (q >0.10). The shape of the rarefaction curve (Supplementary Figure 1 ) demonstrated that the plateauing stage had been reached, indicating that the sequencing depth of all samples was reasonable and the sequencing results could reflect the microbial information in the saliva samples. Alpha diversity among the groups was described by evenness, richness, and diversity. The Chao1 index of the HTN group was considerably lower than that of the control group (p<0.001), whereas there was no significant change in Chao1 index in the OSA and complication groups compared with the control group (p=0.15 and p=0.42, respectively). The Chao1 index was substantially lower in the HTN group compared with the complication group (p=0.014), whereas the Richness index was significantly lower in HTN group compared with control and complication groups (p=0.001 and p=0.0087, respectively). The Shannon index did not reveal a significant variation in alpha diversity among the four groups (p=0.43) ( Figure 1A -C). Variation of salivary microbial community structure was investigated using PCoA based on the Bray-Curtis distance ( Figure 1D and E). Samples obtained from patients in the OSA and complication groups clearly clustered away from those obtained from the control group, highlighting a significant difference in the phylogenetic structures of the salivary microbial communities (Anosim, p=0.026 for control vs OSA group; p=0.014 for control vs complication group, respectively). In addition, there were significant differences among the control, HTN, and complication groups (Anosim, p=0.001 for control vs HTN group; p=0.021 for HTN vs complication group, respectively). The top 10 phyla in terms of abundance are presented in Figure 2A . Firmicutes, Actinobacteria, Proteobacteria, and Bacteroidetes were four components in the saliva microbiome with a relative abundance greater than 10%. At the genus level, Figure 2B shows the top 20 most abundant taxa. The genera with highest abundance in the HTN group were Haemophilus, Lautropia, Neisseria, Fusobacterium, and Aggregatibacter based on the heatmap. The genera and Peptostreptococcus had the highest abundance in the OSA group, while Porphyromonas, Leptotrichia, and Treponema exhibited a higher abundance in the complication group ( Figure 2C ). Compared with the control group, the abundances of Rothia, Granulicatella, and Peptostreptococcus in the OSA group were significantly increased (p=0.027, p=0.013, and p=0.014, respectively), while that of the genus Haemophilus was Compared with the control group, three genera, including Rothia, Neisseria, and Lautropia, were significantly increased in the HTN group (p=0.004, p=0.043, and p<0.001, respectively), while Actinomyces was significantly decreased (p=0.044). Furthermore, compared with the HTN group, the abundances of Lautropia and Haemophilus were significantly decreased in the complication group (p=0.01 and p=0.003, respectively) and that of Peptostreptococcus was significantly increased (p=0.019) (Figure 3 ). Detailed statistics represented in Supplementary Table 1 . The LEfSe method was applied to analyze more specific differences in microflora composition and look for potential biomarkers of the four experimental groups. A cladogram representative of salivary microbial structure and predominant bacteria showing the greatest differences in taxa among the four groups was produced based on an LDA threshold of 2 (Figure 4) . At the genus level, the predominant bacteria of the OSA group were Peptostreptococcaceae__XI__G_1_, Microbial functions were predicted by PICRUSt2 using the KEGG database and indicated that pathways of nitrogen metabolism, citrate cycle, bacterial secretion system, toluene degradation, phenylalanine metabolism, lipopolysaccharide biosynthesis, and fatty acid metabolism had significantly higher abundance in the HTN group compared with the control group. In addition, compared with control group, the pathways of metabolism of xenobiotics by cytochrome P450, polyketide sugar unit biosynthesis, phosphonate and phosphinate metabolism, and proteasome exhibited significantly higher abundance in the OSA group. Furthermore, the KEGG pathways of vasopressin-regulated water reabsorption, steroid biosynthesis, inositol phosphate metabolism, and lipoic acid metabolism presented significantly lower abundance in the complication group compared with the OSA group ( Supplementary Figures 2-6 ). To illustrate the diagnostic value of the salivary microbiome for the OSA and complication groups, random forest classifier models were constructed to identify OSA samples from control samples Solobacterium, and Veillonella were selected as optimal marker genera and the AUC value was 78.5% (95% CI, 65.85-91.13%) ( Figure 5A -C). PCA of the metabolomics samples revealed that the QC samples were closely clustered in the positive mode, indicating that the experimental results were reproducible. In addition, there was a tendency of separation among the groups. Data were interrogated based on the orthogonal partial least-squares discrimination analysis (OPLS-DA), which established intergroup separation between the OSA group and the control group ( Figure 6A) , and between the complication and control groups ( Figure 6B ). The discriminant metabolites that account for the intergroup separation, between the complication and HTN groups, and between the complication group and the OSA group, were also identified by OPLS-DA ( Figure 6C and D) . Four discriminant metabolites were measured in the OSA group relative to the control group, (2-hydroxyadenine, erucic acid, 7-methylxanthine, and pyruvaldehyde), while 26 discriminant metabolites were measured in the HTN group relative to the control group. The levels of 2-hydroxyadenine and 7-methylxanthine were significantly increased in the patients with OSA (p=0.004 and p=0.04, respectively) compared with the control group. Furthermore, there were 6 metabolites that were altered in the complication group compared with the OSA group, while four discriminant metabolites (L-carnitine, Tyr-Gln, N-(omega)-hydroxyarginine, and His-Ala) were measured in the complication group compared with the HTN group; all four of these metabolites were significantly increased in the complication group compared with HTN group.(p=0.038, p=0.039, p=0.045, and p=0.048, respectively). The changes in metabolites are shown in Supplementary Table 2 and Supplementary Figures 7-11 . In recent years, attention has increased towards OSA combined with HTN. 6, 31 This subtype of OSA is common and greatly impacts health. However, these patients have fewer typical OSA symptoms such as sleep disturbance and daytime sleepiness, so it may not be easily detected in clinical practice, and it is also somewhat inconvenient to identify in the case of troublesome polysomnography. Therefore, the identification of biomarkers for this subtype of OSA is desirable and may facilitate a more targeted treatment approach. Studies of the gut microbiome of patients with OSA have revealed that the microbiome is associated with OSA-related HTN, but there are limited studies examining oral microbiomes in patients with OSA. Saliva samples have been routinely taken as a representative average of the entire ecosystem of the oral cavity, therefore, the current study conducted salivary microbiome analysis in patients with OSA combined with HTN to confirm A B C the hypothesis that oral microbiomes might be associated with the distinct subtype of OSA combined with HTN. Alpha diversity of patients with OSA or HTN decreased compared with the control group, indicating a decrease in the richness of the microflora, which was similar to the results of previous studies. 16, 32 Simultaneously, beta-diversity analysis suggested there was significantly different community evolution structures among the different groups of patients. Ko et al found the abundance of the Porphyromonas was significantly higher in oral mucosal samples of OSA patients with or without HTN compared with controls. Furthermore, the abundance of the Syntrophomonas was significantly higher in patients with mild OSA plus HTN compared with that in patients with mild OSA without HTN. 27 We did not find Porphyromonas has significant different among groups probably because saliva samples were taken in our study. Moreover, individuals with HTN but without OSA were recruited in our study to verify whether the microbiological alterations were caused by HTN or OSA. In the current study, some taxa correlated with OSA-related HTN were defined. Rothia was relatively more abundant in all three disease groups (OSA, HTN, and complication) compared with the control group. study reported that Rothia was correlated with apnea-hypopnea of patients with severe OSA, 16 while another study showed that Rothia was present at a high level in the gut microbiome of people with HTN, 33 which was congruent with the current study. Therefore, Rothia may indicate the emergence of a disease state or unhealthy outcomes. In addition, in most of the oral microbial studies on diseases, an increased abundance of Oribacterium is associated with oral lesions. 34, 35 In the current study, Oribacterium was a specific taxa in the complication group compared with control group, and while the underlying mechanism for the OSA-related HTN was not identified, it might be that Oribacterium can be used as a biomarker of the complication group as it was proved in random forest model 2 to distinguish patients of the complication group from those in the HTN group. Furthermore, in the current study, the abundance of Actinomyces was significantly lower in the complication group compared with the OSA group, and as one of the most copious taxa among the nitrate-reducing bacteria, it was also decreased in the HTN group. Nitrate-reducing bacteria such as Actinomyces were found on the dorsum of the tongue, cleaning of the tongue and the use of chlorhexidine mouthwash prevents increased proportion of nitrate reducing species. 12 However, there was no significant difference between the OSA and control groups. Therefore, we speculate that the decrease of Actinomyces in patients with OSA is associated with a potential mechanism of OSA-related hypertension, and it may be possible to predict the tendency to develop HTN in patients with OSA by the reduction of Actinomyces in the salivary microbiome of these individuals. As the second largest biobank in the human body, microbes have attracted increasing attention regarding the occurrence and development of diseases. Various microorganisms exist in human tissues and organs, including the skin, digestive tract, oral cavity, and so on. The oral microbiome is related to nitrate metabolism and cardiovascular health. Approximately 1/4 of the nitrate intake returns to the mouth via an entero-salivary nitrate-nitrite-nitric oxide (NO) pathway. 13 Nitrate-reducing bacteria in oral cavity can reduce nitrate to nitrite, which is absorbed into the bloodstream through the digestive tract and converted into NO, 36,37 a multifunctional signaling molecule involved in cardiovascular homeostasis and an effective endogenous vasodilator. Blockage of the nitrate-nitrite-NO pathway by the use of antimicrobial mouthwash was paralleled by a small increase in SBP in treated men and women with HTN. 38 Findings from the current study are consistent with this theory, hence the development of the hypothesis that the decrease of Actinomyces, which leads to the decrease of NO production, is associated with the increase of blood pressure in patients in the HTN and complication groups. The relative abundance of Absconditabacteria_(SR1) in the oral flora of patients with inflammatory bowel disease (IBD) was higher than normal, and was positively correlated with the levels of inflammatory cytokines such as fecal calprotectin. 39 Absconditabacteria_(SR1) may be related to the inflammatory state caused by OSA, but there was no evidence of a possible mechanism of OSA-related HTN associated with Absconditabacteria_(SR1)_ [G-1] and Peptostreptococcaceae_[XI] [G-1] in the previous literature, therefore we regarded these taxa more likely as outcome indicators. Streptococcus and Prevotella are known to be related to the production of acetate, one of the short-chain fatty acids (SCFAs) that result from anaerobic fermentation of dietary fibers by gut microbiota and which has an important vasodilator effect. 40, 41 However, no significant differences of Streptococcus and Prevotella were identified among the four groups in the current study. In addition, Haemophilus, Neisseria, and Lautropia were speculated to be biomarkers of HTN. Haemophilus, another nitrate-reducing bacteria, 13 had a significantly lower abundance in the complication group compared with the HTN group, while simultaneously presented a significantly lower level in the OSA group compared with the control group. Consequently, we speculate that OSA might cause HTN by reducing the production of Haemophilus. Furthermore, the abundance of Neisseria in the HTN group was significantly higher than that in the control group, which was consistent with previous studies on intestinal microorganisms in patients with HTN. 18 Concurrently, in a study of oral microorganisms and salivary NO, Neisseria was more abundant in subjects with HTN compared with normal subjects, suggesting that the concentration of this bacteria was a biomarker for the development of HTN. 14 Moreover, Lautropia exhibited a significantly higher relative abundance in the HTN group compared with the control group in the current study. This was consistent with a previous study showing that Lautropia was higher in abundance in pharyngeal swabs from patients with pulmonary HTN that from control samples. 32 In OSCC cell lines, Lautropia mirabilis upregulated IL-6 and TNF-α, which are cytokines related to the inflammatory response. 42 All these genera were verified in LEfSe analysis in the current study. Bacterial metabolites may also be responsible for OSA-related HTN. Oral microorganisms may contribute to the occurrence of cardiovascular disease through direct invasion, immune responses, promotion of platelet aggregation, inflammatory responses, and oxidative stress responses. 43 OSA may induce related inflammatory factors such as C-reactive protein, TNF-α, and IL-6, 44 through oral bacterial pathogens, like Rothia, Haemophilus, and Actinomyces, which can promote the occurrence and development of HTN. DNA and lipopolysaccharide from oral bacterial pathogens activate nuclear factor-κB (NF-κB) and activator protein 1 (AP1), which lead to increased production of inflammatory factors and activation of neutrophils, leading to the production of numerous reactive oxygen species (ROS). 45 Intermittent hypoxia induced by OSA induces systemic inflammation and sympathetic hyperactivity by enhancing the production of ROS. The ROS-HIF-1-endothelium signaling pathway and the ROS-TLR-NF-κB signaling pathway are responsible for hypertension and cardiovascular dysfunction in patients with OSA. 46, 47 In addition, hypoxia induced by OSA can lead to mitochondrial dysfunction, which affects related enzyme systems, such as endothelial NO synthase production. 48 Previous study compared clinical chemistry measurements of enzymes, metabolites between patients with HTN and controls, 49 few studies used metabolomic to analysis metabolites in patients with HTN, especially OSA related HTN. Ultra-high-performance liquid chromatography was used for metabolomic analysis to screen for differential metabolites in our study, including indicators that are not routinely used in clinical. Metabolomic analysis found evidence of cardiovascular damage associated with OSA. As previously stated, OSA enhances ROS production. 2-Hydroxyadenine, a purine produced in cells by reactive oxygen, 50 was significantly increased in the group of patients with OSA in the current study. Meanwhile, L-carnitine was significantly decreased in the complication group compared with the OSA group, and this metabolite triggers cardioprotective effects through reduced inflammation, oxidative stress, and necrosis of cardiac myocytes, and can also reduce HTN. 51 The decrease of L-carnitine in the complication group may be responsible for the increase of blood pressure in these patients. This study does have some limitations. OSA may have an impact on multiple systems of the whole body. In the study, HTN was confirmed with 24-h ambulatory blood pressure monitoring, OSA was confirmed with full-night hospitalization PSG, and the medical history and complications of all subjects were carefully taken in the study. However, it is inevitable that other complications, although very minor, might also interfere with the data. Such as age and BMI, thought they were proved no influence on the comparison among groups they should be better matched in the future study. Moreover, salivary flora has the advantage of easy collection and may induce changes in the gut microbiome through swallowing, but the performance of bacteria in the mouth is not necessarily identical to their role(s) in the intestine. This study allows the screening different flora, but causal relationships between salivary flora and disease require further study by verification of CPAP therapy intervention or through in vivo experiments. In addition, this study does not represent the performance of other genders, ages, and regional races. This study demonstrates the characteristic changes in salivary microbiota of patients with OSA, patients with HTN, and patients with both OSA and HTN. Findings from the study led to the hypothesis that the decrease of Actinomyces may be associated with OSA-related HTN. To identify specific microbial markers for the different groups of patients, models were constructed to distinguish all-OSA patients from control subjects, patients with OSA and HTN (complication group) from patients with OSA, and patients with OSA and HTN (complication group) from patients with HTN. The diagnostic efficacy of these models was validated and highlights potential noninvasive diagnostic biomarkers for the subtype of OSA that is associated with HTN. Moreover, bacterial metabolites may also be responsible for OSA-related hypertension. OSA produced ROS can increased 2-hydroxyadenine indirectly. And the decrease of L-Carnitine may be associated with the hypertension in OSA patients. OSA, obstructive sleep apnea; HTN, hypertension; BMI, body mass index; AHI, apnea hypopnea index; SBP, systolic blood pressure; DBP, diastolic blood pressure. High blood pressure and cardiovascular disease Obstructive sleep apnea in adults Obstructive sleep apnea: the most common secondary cause of hypertension associated with resistant hypertension Relationship of daytime blood pressure and severity of obstructive sleep apnea among Chinese: a multi-center investigation in China The different clinical faces of obstructive sleep apnoea: a cluster analysis Phenotypes of comorbidity in OSAS patients: combining categorical principal component analysis with cluster analysis Pathophysiology of sleep apnea The role of the microbiome in human health and disease: an introduction for clinicians The microbiome and breast cancer: a review The influence of the gut microbiome on cancer, immunity, and cancer immunotherapy The human microbiome and obesity: moving beyond associations Current insights in microbiome shifts in Sjogren's syndrome and possible therapeutic interventions How periodontal disease and presence of nitric oxide reducing oral bacteria can affect blood pressure Association between hypertension, oral microbiome and salivary nitric oxide: a case-control study The microbiome in obstructive sleep apnea Severe obstructive sleep apnea is associated with alterations in the nasal microbiome and an increase in inflammation Oropharyngeal microbiome in obstructive sleep apnea: decreased diversity and abundance Role of the gut microbiome in obstructive sleep apnea-induced hypertension Prebiotics, probiotics, and acetate supplementation prevent hypertension in a model of obstructive sleep apnea Obstructive sleep apnea and systemic hypertension: gut dysbiosis as the mediator? Lactobacillus rhamnosus GG strain mitigated the development of obstructive sleep apnea-induced hypertension in a high salt diet via regulating TMAO level and CD4(+) T cell induced-type I inflammation The oral microbiota -a mechanistic role for systemic diseases Role of the oral microbiota in cancer evolution and progression Links between atherosclerotic and periodontal disease Elevated levels of circulating cell adhesion molecules in uncomplicated essential hypertension Oral bacteria affect the gut microbiome and intestinal immunity Analysis of oral microbiota in patients with obstructive sleep apnea-associated hypertension 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 NMA/PCNA Guideline for the prevention, detection, evaluation, and management of high blood pressure in adults: executive summary: a report of the American college of cardiology/American heart association task force on clinical practice guidelines Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2 Polysomnographic phenotypes and their cardiovascular implications in obstructive sleep apnoea Altered airway microbiota composition in patients with pulmonary hypertension Characteristics and variation of fecal bacterial communities and functions in isolated systolic and diastolic hypertensive patients Betel nut chewing, oral premalignant lesions, and the oral microbiome Dysbiosis of saliva microbiome in patients with oral lichen planus Inorganic nitrate supplementation lowers blood pressure in humans: role for nitrite-derived NO Nitrate, bacteria and human health Antibacterial mouthwash blunts oral nitrate reduction and increases blood pressure in treated hypertensive men and women High-throughput sequencing provides insights into oral microbiota dysbiosis in association with inflammatory bowel disease Dissecting the in vivo metabolic potential of two human gut acetogens Microbial short chain fatty acid metabolites lower blood pressure via endothelial G protein-coupled receptor 41 Screening of health-associated oral bacteria for anticancer properties in vitro A review of the association between oral bacterial flora and obstructive sleep apnea-hypopnea syndrome comorbid with cardiovascular disease Systemic inflammation: a key factor in the pathogenesis of cardiovascular complications in obstructive sleep apnoea syndrome? Oxidative stress, systemic inflammation, and severe periodontitis Co-upregulation of Toll-like receptors 2 and 6 on peripheral blood cells in patients with obstructive sleep apnea Endothelin contributes to the blood pressure rise triggered by hypoxia in severe obstructive sleep apnea Molecular mechanisms of cardiovascular disease in OSAHS: the oxidative stress link Dysbiosis of the salivary microbiome is associated with hypertension and correlated with metabolic syndrome biomarkers Mutagenicities of 8-hydroxyguanine and 2-hydroxyadenine produced by reactive oxygen species L-carnitine and heart disease Nature and Science of Sleep is an international, peer-reviewed, open access journal covering all aspects of sleep science and sleep medicine, including the neurophysiology and functions of sleep, the genetics of sleep, sleep and society, biological rhythms, dreaming, sleep disorders and therapy, and strategies to optimize healthy sleep. The manuscript management system is completely online and includes a very quick and fair peer-review system We thank Dr Bing Zhou, Yuhua Zuo, Xueli Zhang, Long Zhao, Jing Li (Sleep Center of Peking University People's Hospital), Dr Qingchun Ma (Cardiology Department of Peking University People's Hospital), Dr Liping Song, and Ping Han (Cardiology Department of Beijing Haidian Hospital) for help in sample collection and generous volunteer subjects who recruited in this study. We also thank all the participants in the trial. All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas: took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work. This work was supported partly by National Program for Multidisciplinary Cooperative Treatment on Major Diseases (PKUSSNMP-201902), National Natural Science Foundation of China (82170102), and Science Foundation of HaiDian Hospital (KYM2017014). The sequencing data from this study have been deposited in Genome Sequence Archive in National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA006039) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa. The authors declare that they have no conflicts of interest in this work.