Front. Med. Frontiers in Medicine Front. Med. 2296-858X Frontiers Media S.A. 10.3389/fmed.2022.813459 Medicine Original Research The Discovery, Validation, and Function of Hypoxia-Related Gene Biomarkers for Obstructive Sleep Apnea Wu Xiaofeng Pan Zhou Liu Wei Zha Shiqian Song Yan Zhang Qingfeng Hu Ke * Department of Respiratory and Critical Care Medicine, Renmin Hospital of Wuhan University, Wuhan, China

Edited by: Pedro Pereira Rodrigues, Universidade do Porto, Portugal

Reviewed by: Marta Drummond, CRI Sono e VNI Centro Hospitalar e Universitário São João, Portugal; Silvia V. Conde, New University of Lisbon, Portugal

 *Correspondence: Ke Hu huke-rmhospital@163.com

This article was submitted to Translational Medicine, a section of the journal Frontiers in Medicine

 17 03 2022 2022 9 813459 11 11 2021 10 02 2022 Copyright © 2022 Wu, Pan, Liu, Zha, Song, Zhang and Hu. 2022 Wu, Pan, Liu, Zha, Song, Zhang and Hu

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

While there is emerging evidence that hypoxia critically contributes to the pathobiology of obstructive sleep apnea (OSA), the diagnostic value of measuring hypoxia or its surrogates in OSA remains unclear. Here we investigated the diagnostic value of hypoxia-related genes and explored their potential molecular mechanisms of action in OSA. Expression data from OSA and control subjects were downloaded from the Gene Expression Omnibus database. Differentially-expressed genes (DEGs) between OSA and control subjects were identified using the limma R package and their biological functions investigated with the clusterProfiler R package. Hypoxia-related DEGs in OSA were obtained by overlapping DEGs with hypoxia-related genes. The diagnostic value of hypoxia-related DEGs in OSA was evaluated by receiver operating curve (ROC) analysis. Random forest (RF) and lasso machine learning algorithms were used to construct diagnostic models to distinguish OSA from control. Geneset enrichment analysis (GSEA) was performed to explore pathways related to key hypoxia-related genes in OSA. Sixty-three genes associated with hypoxia, transcriptional regulation, and inflammation were identified as differentially expressed between OSA and control samples. By intersecting these with known hypoxia-related genes, 17 hypoxia-related DEGs related to OSA were identified. Protein-protein interaction network analysis showed that 16 hypoxia-related genes interacted, and their diagnostic value was further explored. The 16 hypoxia-related genes accurately predicted OSA with AUCs >0.7. A lasso model constructed using AREG, ATF3, ZFP36, and DUSP1 had a better performance and accuracy in classifying OSA and control samples compared with an RF model as assessed by multiple metrics. Moreover, GSEA revealed that AREG, ATF3, ZFP36, and DUSP1 may regulate OSA via inflammation and contribute to OSA-related cancer risk. Here we constructed a reliable diagnostic model for OSA based on hypoxia-related genes. Furthermore, these transcriptional changes may contribute to the etiology, pathogenesis, and sequelae of OSA.

 obstructive sleep apnea hypoxia machine learning diagnosis AREG ATF3 ZFP36 DUSP1 Introduction

Obstructive sleep apnea (OSA), characterized by partial or total collapse of the upper airway during sleep, is a major and underappreciated public health problem. Untreated OSA patients are at significantly increased risk of multiple diseases and traffic accidents, and the economic burden of OSA increases with disease severity (1, 2). The prevalence of OSA is increasing as the global obesity epidemic worsens, and OSA currently affects 9% of women and 24% of men (3). The gold standard clinical diagnosis of OSA is based on polysomnography (PSG), but due to the need for specialist equipment, long-term continuous monitoring, and specialist input and training, PSG is not widely carried out in the community. Although other methods of risk stratification and screening for OSA exist, such as the Berlin questionnaire, expiratory NO detection and face recognition technology, these approaches are not recommended for the diagnosis and treatment of OSA (4). Furthermore, while home sleep apnea testing has high sensitivity [79% (95% CI, 71–86%)] and specificity [79% (95% CI, 63–89%)], it still takes all night and is costly (5).Therefore, new, low-cost, easily applicable diagnostic methods for OSA are required to replace the need for PSG.

Hypoxia is best defined as an imbalance between oxygen supply and demand, which varies widely depending on the organ or tissue in the body (6). The intermittent hypoxia (IH) caused by sleep apnea characterizes OSA (7). Similar to ischemia-reperfusion injury, chronic intermittent hypoxia (CIH) can lead to abundant reactive oxygen species (ROS) formation through repeated cycles of systemic hypoxia/reoxygenation (8) and increase blood pressure by regulating sympathetic activity (9). In both animal and human experiments, nerve damage and the slow autonomic response system caused by IH in OSA result in an ability to restore cardiovascular regulation (10). In addition, OSA-related IH is a risk factor for the development and progression of lung cancer. Thus, investigating the molecular basis of hypoxia and its effects in OSA could improve our understanding of the pathogenesis of OSA and related diseases (11–13).

Therefore, here we evaluated hypoxia-related gene expression in OSA to (i) explore the diagnostic value of hypoxia-related genes in OSA patients; and (ii) investigate their possible contribution to OSA pathobiology using in silico analyses.

 Materials and Methods Data Sources

The transcriptomic profiles of 34 OSA and eight control subjects were retrieved from GSE135917 in the Gene Expression Omnibus (GEO) database (www.ncbi.nlm.nih.gov/geo/). In this public dataset, the diagnosis of OSA patients were based on AHI≥30 events/h or RDI≥5 events/h, and the control subjects had normal AHI or RDI. Thirty-four OSA subjects, including 16 males and 18 females, had a mean age of 51.3 years with a mean BMI of 40.7 kg/m2. When diagnosis, 10 OSA patients were enrolled by RDI and their average RDI was 19.2 events/h, while the other 24 OSA patients were recruited by AHI with the average of 41.5 events/h. Besides, eight OSA subjects had diabetes mellitus, 12 subjects had hypertension disease, and two subjects had heart disease. Control subjects, including one male and seven females, had a mean age of 54.5 years with a mean BMI of 35.2 kg/m2. Among them, three control subjects had hypertension disease. Two hundred hypoxia-related genes were obtained from the HALLMARK_HYPOXIA gene set in the MSigDB database (www.gsea-msigdb.org/gsea/msigdb/cards/HALLMARK_HYPOXIA.html). The GSE135917 dataset was used to detect hypoxia-related genes involved in OSA, construct a diagnostic model, and identify potential diagnostic biomarkers for OSA. Ten OSA and eight control subjects from GSE38792 were used to validated the OSA diagnostic model.

 Identification of Hypoxia-Related Genes Involved in OSA

Differentially-expressed genes (DEGs) between OSA and control subjects in GSE135917 were identified using the limma package in R with the threshold |log2 FC| >1 and p < 0.05. The functional enrichment of DEGs in biological pathways and processes {according to the Gene Ontology [GO; biological process (BP), cellular component (CC), and molecular function (MF)] and KEGG databases} was analyzed using clusterProfiler in R. Thereafter, hypoxia-related genes involved in OSA were identified by overlapping DEGs with the 200 hypoxia-related genes identified in MSigDB. Interactions between hypoxia-related genes involved in OSA were investigated by inputting them into the STRING database (www.string-db.org/) to construct a protein-protein interaction (PPI) network, and genes within the PPI network were evaluated by correlation analysis.

 Identification of Potential Diagnostic Biomarkers for OSA

The diagnostic value of genes within the PPI network was evaluated by area under the receiver operating characteristic (AUC-ROC) curves. Hypoxia-related genes involved in OSA with AUCs >0.7 were further explored as potential diagnostic biomarkers for OSA. Moreover, micro-(mi)RNA-diagnostic biomarker and transcription factor (TF)-diagnostic biomarker regulatory networks were constructed using the mirnet database (https://www.mirnet.ca). Furthermore, geneset enrichment analysis (GSEA) was performed using KEGG genesets as reference to reveal potential mechanisms by which the identified diagnostic biomarkers might regulate OSA. Reference genesets were downloaded from the MSigDB database (http://software.broadinstitute.org/gsea/msigdb/index.jsp).

 Identification of Diagnostic Biomarkers and Development of an OSA Diagnostic Model

Machine learning methods including random forest (RF) (14) and lasso regression (15) were applied to construct a diagnostic model based on the OSA diagnostic biomarkers identified in the GSE135917 dataset using the randomForest and glmnet R packages, respectively. ROC curves and their derived sensitivities, specificities, and AUCs were constructed to evaluate the accuracy of the established diagnostic models. To compare the performance of the RF and lasso algorithms, accuracy, error, precision, recall, F1-score, and Kolmogorov-Smirnov (K-S) values were calculated using the following formulae: accuracy = (TP + TN)/(TP + FP + TN + FN); error = 1-Accuracy; precision = TP/(TP + FP); recall = TP/(TP + FN); F1-score = 2 × precision × recall/(precision + recall); K-S = max[TP/(TP + FN)-FP/(FP + TN)], where TP is true positive, FN is false negative, FP is false positive, and TN is true negative. The diagnostic model showing the best performance for classifying OSA and control subjects was selected and validated in the GSE38792 dataset.

 Results DEGs Associated With Hypoxia and Inflammation in OSA

Sixty-three DEGs between OSA and control subjects were identified in GSE135917 (Supplementary Table 1). Among them, 11 genes were upregulated and 52 were downregulated in OSA subjects relative to control (Figure 1A). The expression levels of the top ten upregulated (EIF1AY, USP9Y, UTY, ERAP2, IMPAD1, PPBP, HBA2, HBB, TM4SF19, EGFL6) and downregulated (PKGD1L1, LINC00917, IL6, MIR21, DUSP1, CD69, PTGS2, FOS, FOSB, EGR1) genes in OSA are shown in the heatmap in Figure 1B. Functional enrichment analysis showed that those DEGs were significantly enriched in 245 BP, 2 CC, 36 MF (Supplementary Table 2), and 53 KEGG pathways (Supplementary Table 3). The top GO terms were associated with hypoxia and transcription such as fat cell differentiation, response to oxidative stress, response to glucocorticoids, transcription regulator complex, DNA-binding transcription activator activity, and DNA-binding transcription factor binding (Figure 1C). Interestingly, DEGs were also involved in inflammation-which has complex interactions with hypoxia-such as IL-17 signaling, TNF signaling, and MAPK signaling (Figure 1D). In addition, there was enrichment for diseases closely associated with OSA such as rheumatoid arthritis and non-alcoholic fatty liver disease (Figure 1D).

Identification and functional analysis of differentially expressed genes between OSA and control subjects. (A) Volcano plot of the DEGs, including 11 upregulated genes in red and 52 downregulated genes in green, between OSA and control subjects. (B) Heatmap showing the expression of the top 10 upregulated genes and top 10 downregulated genes in OSA compared with control subjects. (C) Bar chart showing the GO enrichment analysis of DEGs, including biological processes, molecular functions, and cellular components. (D) Bar chart showing KEGG enrichment analysis of DEGs.

 Hypoxia-Related Genes Involved in OSA

In view of the identified DEGs and their functional enrichment, we hypothesized that hypoxia plays an important role in OSA. Thus, we extracted 17 hypoxia-related genes involved in OSA by intersecting the 63 DEGs with 200 hypoxia-related genes in MSigDB (Figure 2A); all 17 genes were significantly downregulated in OSA subjects compared with control subjects (Supplementary Figure 1). PPI analysis showed that FOS, IL6, ATF3, JUN, DUSP1, EGR1, PTGS2, FOSB, CXCL2, NR4A1, SOCS3, ZFP36, NR4A2, KLF4, AREG, and CD69 interacted with each other (Figure 2B), with FOS and IL6 probable hub nodes, followed by ATF3 and JUN, since they had the most interactions with hypoxia-related genes involved in OSA (Figure 2C). Moreover, the expression levels of the 16 genes in the PPI network were strongly positively correlated (≥0.85, Figure 2D), further demonstrating their positive interactions.

Identification and characterization of hypoxia-related DEGs in OSA. (A) Venn diagram of overlapping hypoxia-related genes and DEGs. (B) PPI network of hypoxia-related DEGs. (C) Number of adjacent nodes of hypoxia-related DEGs in the PPI network. (D) Heatmap of correlations between hypoxia-related DEGs in the PPI network.

 Hypoxia-Related Genes are Accurate Biomarkers for the Diagnosis of OSA

We next explored the diagnostic value of the 16 hypoxia-related genes in the PPI network. The AUCs of FOS, IL6, ATF3, JUN, DUSP1, EGR1, PTGS2, FOSB, CXCL2, NR4A1, SOCS3, ZFP36, NR4A2, KLF4, AREG, and CD69 were all >0.7 (Figure 3), indicating that they might be useful diagnostic biomarkers for OSA. Next, we investigated which miRNAs and TFs might regulate their expression. Using mirnet, we detected 290 miRNA-diagnostic biomarker pairs, which allowed us to construct an miRNA-mRNA network composed of 16 hypoxia-related genes and 85 miRNAs (Figure 4A). We also identified 198 TF-diagnostic biomarker pairs, which produced a TF-target network composed of 15 diagnostic biomarkers and 118 TFs (Figure 4B).

ROC curves showing the diagnostic value of hypoxia-related DEGs in the PPI network.

Regulation network of hypoxia-related DEGs in the PPI network. (A) Construction of the miRNA-mRNA network. Hypoxia-related DEGs are represented by magenta triangles, and the corresponding miRNAs are represented by blue triangles. (B) Construction of the TF-mRNA network. Hypoxia-related DEGs are represented by magenta triangles, and the TFs targeting hypoxia-related DEGs are represented by green diamonds.

 Lasso Diagnostic Model for OSA Based on Hypoxia-Related Genes Involved in OSA

We next used the 16 hypoxia-related genes in the PPI network to construct OSA diagnostic models. The confusion matrix of the RF model showed correct classification of 30 OSA patients and four control subjects (Figure 5A) with an AUC of 0.667 (Figure 5B). Similarly, using lasso, AREG, ATF3, ZFP36, and DUSP1 were identified as a predictive gene signature (Figures 6A,B), and the diagnostic OSA model based on these four genes had an AUC of 0.842 (Figure 6C). The confusion matrix of the lasso model also demonstrated its superior performance over the RF model, with 32 OSA patients and five control subjects correctly identified (Figure 6D). Moreover, as shown in Table 1, the lasso model had higher precision and K-S values, indicating that it was better able to distinguish OSA from control subjects than the RF model. Therefore, the lasso model was selected as the OSA diagnostic model of choice, with validation in the external dataset showing an AUC of 0.787 (Figure 6E).

Establishment of a random forest model based on hypoxia-related DEGs in OSA. (A) Confusion matrix of GSE135917 using the random forest model. (B) The diagnostic performance of the random forest model by ROC curve analysis.

Establishment and verification of the lasso model. (A,B) The lasso algorithm was used to screen the gene signature in OSA. (C) The diagnostic performance of the lasso model by ROC curve analysis in GSE135917. (D) Confusion matrix of GSE135917 using the lasso model. (E) Verification of the diagnostic performance of the lasso model by ROC curve analysis in GSE38792.

The diagnostic accuracy of the lasso and random forest models.

Model Precision Recall Fl-score Accuracy Error K-S
Random forest 0.88 0.88 0.79 0.81 0.19 0.38
Logistic regression 0.91 0.94 0.91 0.88 0.12 0.535
GSEA Reveals a Potential Role for Hypoxia-Related Genes in OSA

Next, we explored the role of AREG, ATF3, ZFP36, and DUSP1 in OSA by GSEA. According to the expression of each gene, OSA patients were divided into low- and high-expression groups. Interestingly, AREG, ATF3, ZFP36, and DUSP1 were all associated with inflammation and cancer, specifically enrichment of adipocytokine signaling, bladder cancer, epithelial cell signaling in Helicobacter pylori infection, and JAK/STAT signaling pathways in patients with high AREG expression (Figure 7A). ATF3 expression was related to glycosaminoglycan degradation, pathways in cancer, prostate cancer, and TGF-β signaling (Figure 7B), and DUSP1 was associated with adipocytokine signaling, graft vs. host disease, pathways in cancer, and small cell lung cancer (Figure 7C). Finally, ZFP36 was associated with cytokine-cytokine receptor interactions, glycosphingolipid biosynthesis (lacto and neolecto series), hypertrophic cardiomyopathy, and pathways in cancer (Figure 7D). REG, ATF3, ZFP36, and DUSP1 therefore have both common and distinct roles in OSA.

GSEA analysis of functions related to AREG (A), ATF3 (B), DUSP1 (C), and ZFP36 (D).

 Discussion

The diagnostic value of hypoxia or its surrogates in OSA remains unclear. Here we investigated the diagnostic value of hypoxia-related genes in OSA and analyzed potential mechanisms by which the identified key hypoxia-related genes might be involved in OSA pathobiology.

Sixty-three genes were differentially expressed between control and OSA samples, which were significantly enriched in hypoxia-related pathways such as oxidative stress. Oxidative stress refers to an imbalance between anti-oxidants and pro-oxidants, resulting in the disruption of redox signaling and/or molecular damage as a consequence of ROS and reactive nitrogen species (RNS) overproduction (16). In OSA, the upper airways of patients collapse repeatedly during sleep, resulting in intermittent hypoxia and reoxygenation (17). This hypoxia/reoxygenation process is an important mechanism of target organ injury in OSA and it produces a large amount of ROS, leading to oxidative stress (18). Our results further demonstrate an important role for hypoxia in OSA. In addition, the identified genes were significantly enriched in inflammation-related signaling pathways such as IL-17 signaling, TNF signaling, and MAPK signaling. IL-17 is mainly secreted by Th17 cells, where it binds to cell surface receptors to activate downstream cellular signal transduction pathways and induce the secretion of pro-inflammatory cytokines (19). IL-17 also participates in the proliferation of neutrophils and dendritic cells (DC), immune cell maturation and chemotaxis, and synergizes with some cytokines to activate T cells, thereby amplifying the immune response and inflammatory destruction of target organs and contributing to several autoimmune and infectious diseases (20). In their analysis of adenoid hypertrophy of children with OSA, Sade et al. found a negative correlation between Th17/T cell subsets and OSA symptom severity (21). During long-term chronic hypoxia in adult OSA, Th17/Tregs become imbalanced, and Th17 and associated inflammatory cytokines may be involved in the occurrence and development of OSA (22). Therefore, IL-17 signaling may play an important role in the pathogenesis of OSA and its complication.

TNF-α plays an important role in sleep regulation, and excessive daytime drowsiness, one of the main consequences of OSA, is accompanied by an increase in TNF-α levels (23–25). In addition, TNF-α is a crucial player in immune and inflammatory responses (26), which may be involved in IH-driven OSA-related cardio-metabolic disease (27, 28). The MAPK signaling pathway plays an important role in central injury and regulatory mechanisms induced by CIH (29). Activated by hypoxia and other stimuli, MAPK enters the nucleus to regulate transcription and participate in cell proliferation, development, apoptosis, as well as the recognition, transmission, and amplification of a variety of biochemical response signals between cells (30), indicating that MAPK may mediate the OSA-related changes in cellular physiology caused by IH.

By intersecting DEGs with hypoxia-related genes, we obtained 17 hypoxia-related genes. Among them, we found that the expression of IL6 was significantly down-regulated in OSA samples compared with control samples, which is inconsistent with previous finding in OSA (31). According to the BMI in GSE135917, both OSA and control samples were obese and under chronic inflammation, which may affect the expressions of IL6. Besides, we also found that some of OSA and control individuals were taking medicines, such as duloxetine, insulin and vitamin D, which may either elevate or reduce the expression of IL6 (32–34). Those characteristics of OSA and control samples used in the current study may cause the conflicted results with previous finding. We subsequently constructed a diagnostic model that accurately distinguished between OSA and control samples based on four hypoxia-related genes: AREG, ATF3, ZFP36, and DUSP1. AREG, a member of the epidermal growth factor family (35), is an autocrine growth factor and mitogen for astrocytes, Schwann cells, and fibroblasts (36). It has been reported that hypoxia promotes AREG expression in intestinal epithelial cells in a cAMP response element-dependent manner (37). Furthermore, Yoji et al. detected significant increases in AREG transcript in intermittent hypoxia-treated rat vascular smooth muscle cells (38). Its role in OSA remains unknown but, tor the first time, our study found that AREG1 is a potential diagnostic biomarker with high accuracy (AUC 0.856).

The expression of ATF3, a member of the CREB/ATF family, can be induced by hypoxia (39, 40), and numerous studies have reported its important role in hypoxia-related diseases (41, 42) but not in OSA. Jiang et al. demonstrated that hypoxia or ATF3 overexpression impaired mitochondrial function in differentiated 3T3-L1 fibroblast-like cells, and ATF3 deletion partially recovered hypoxia-mediated mitochondrial dysfunction, suggesting that ATF3 may play a role in adipocyte hypoxia-mediated mitochondrial dysfunction in obesity (39). However, its exact role in OSA requires further in vivo and in vitro investigation. ZFP36 plays a significant role in regulating immune responses and inflammatory diseases by inhibiting the production of various inflammatory cytokines such as TNF-α in macrophages (43). Recent data have indicated that inflammation may play an important role in OSA (44). Zhang et al. found that the genotype distribution of single nucleotide polymorphisms rs251864 and rs17879933 in ZFP36 were significantly different between moderate-to-severe OSA patients and control controls, but they were not directly associated with sleep apnea parameters, suggesting that ATF may be involved in regulating OSA in cooperation with other factors (45). It has also been reported that overnight IH in patients with OSA induces expression of DUSP1, which may mediate increases in manganese superoxide dismutase (MnSOD) expression and activity, contributing significantly to neutralizing the effects of the ROS elicited by OSA (46). Collectively, AREG, ATF3, ZFP36, and DUSP1 are promising biomarkers for the diagnosis of OSA and are possibly functional in OSA pathogenesis. However, the detailed mechanisms of how they participate in OSA need to be established in future studies.

Finally, to obtain further insights into potentially relevant mechanisms in OSA pathogenesis, we used GSEA to analyze the signaling pathways related to AREG, ATF3, ZFP36, and DUSP1, and found that all four markers were related to inflammation and cancer-related pathways. For example, AREG and DUSP2 were associated with adipocytokine signaling. There is increasing evidence that IH mediates at least some of its detrimental effects through adipose tissue inflammation and dysfunction (47), and adipocytokines secreted by adipocytes such as leptin, chemerin, resistin, adiponectin, and omentin 1 are closely associated with OSA (48). ATF3 was related to TGF-β signaling, which has been reported to be necessary for nuclear accumulation of the ATF3/cJun transcription complex and induction of pro-inflammatory responses (49). Compared with control subjects, Lin et al. found that OSA patients had higher TGF-β concentrations in exhaled breath condensates. These reports suggest that ATF3 may function in OSA via TGF-β signaling-mediated inflammation. Similarly, cytokine-cytokine interactions were significantly enriched in the high ZPF36-expressing group. Elevated levels of several cytokines, such as IL6 and TNF-α, are a common feature of OSA (50), suggesting that their interactions may also be important in OSA. Our results also reinforce the hypothesis that OSA may be considered a systemic inflammatory disease, with an imbalance in pro/anti-inflammatory cytokines contributing to its development and the progression of comorbidities. Of note, AREG, ATF3, DUSP1, and ZFP36 were all related to cancers or pathways in cancer, and numerous studies have reported a role for these genes in different tumors (51–54). OSA is a recognized risk factor for cancer, mainly through hypoxia (55), but the exact mechanisms remain unclear. Combined with findings from our current study, we speculate that AREG, ATF3, DUSP1, and ZFP36 may serve as a bridge between cancer and OSA.

This study helps to clarify the relationship between OSA hypopnea syndrome and hypoxia genes. However, it also has some limitations. First, this research is based on a retrospective analysis of existing databases, so prospective studies are needed to verify the in silico results. Second, our goal was to find a practical diagnostic model, but exactly how to use and measure these biomarkers in clinical practice requires further exploration. And it is worth mentioning that the OSA samples used in the current study were obese, and obesity is closely related to hypoxia. Therefore, the relationships between hypoxia-related genes in the diagnostic model and OSA may be due to obesity. Further collecting samples from subgroups of lean and obese OSA patients are needed to validate the clinical use of the current diagnostic model. Third, although OSA and related hypoxia genes were identified, their exact biological functions in OSA require direct expenormalrimental verification.

In conclusion, this study explored the diagnostic value of hypoxia-related genes in OSA, constructed a diagnostic model that can help to distinguish OSA from subjects, and unveiled potential molecular mechanisms of key hypoxia-related genes in OSA. Our findings may pave the way for accurately diagnosing OSA without the need for PSG.

 Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: the OSA-related data employed in this research were obtained from the GEO database GSE135917 (https://www.ncbi.nlm.nih.gov/geo/). Hypoxia-related genes were obtained from HALLMARK_HYPOXIA from the MSigDB database (http://www.gsea-msigdb.org/gsea/msigdb/cards/HALLMARK_HYPOXIA.html). The PPI network was from the STRING database (www.string-db.org/).

 Author Contributions

KH supervised and conceptualized the study. XW and ZP performed most of the bioinformatics analysis. WL, SZ, and QZ participated in the bioinformatics and statistical analyses. XW and KH wrote the manuscript. All authors read and approved the final manuscript.

 Funding

This study was supported by the National Natural Science Foundation of China (No. 81970082).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2022.813459/full#supplementary-material

The expressions of hypoxia-related DEGs in OSA and control subjects.

Identified DEGs between OSA and control subjects in GSE135917.

GO enrichment analysis of DEGs identified in GSE135917.

KEGG pathway enrichment analysis of DEGs identified in GSE135917.

 References Lyons MM Bhatt NY Pack AI Magalang UJ. Global burden of sleep-disordered breathing and its implications. Respirology. (2020) 25:690–702. 10.1111/resp.1383832436658 Léger D Stepnowsky C. The economic and societal burden of excessive daytime sleepiness in patients with obstructive sleep apnea. Sleep Med Rev. (2020) 51:101275. 10.1016/j.smrv.2020.10127532169792 Sarkissian L Kitipornchai L Cistulli P Mackay SG. An update on the current management of adult obstructive sleep apnoea. Aust J Gen Pract. (2019) 48:182–6. 10.31128/AJGP-10-18-474931256486 Kapur VK Auckley DH Chowdhuri S Kuhlmann DC Mehra R Ramar K . Clinical practice guideline for diagnostic testing for adult obstructive sleep apnea: an American academy of sleep medicine clinical practice guideline. J Clin Sleep Med. (2017) 13:479–504. 10.5664/jcsm.650633524841 Gottlieb DJ Punjabi NM. Diagnosis and management of obstructive sleep apnea: a review. JAMA. (2020) 323:1389–400. 10.1001/jama.2020.351432286648 McGarry T Biniecka M Veale DJ Fearon U. Hypoxia, oxidative stress and inflammation. Free Radic Biol Med. (2018) 125:15–24. 10.1016/j.freeradbiomed.2018.03.04229601945 Cori JM O'Donoghue FJ Jordan AS. Sleeping tongue: current perspectives of genioglossus control in healthy individuals and patients with obstructive sleep apnea. Nat Sci Sleep. (2018) 10:169–79. 10.2147/NSS.S14329629942169 Zhao YS An JR Yang S Guan P Yu FY Li W . Hydrogen and oxygen mixture to improve cardiac dysfunction and myocardial pathological changes induced by intermittent hypoxia in rats. Oxid Med Cell Longev. (2019) 2019:7415212. 10.1155/2019/741521230984338 Dewan NA Nieto FJ Somers VK. Intermittent hypoxemia and OSA: implications for comorbidities. Chest. (2015) 147:266–74. 10.1378/chest.14-050025560865 Romero D Jane R. Non-linear HRV analysis to quantify the effects of intermittent hypoxia using an OSA rat model. Annu Int Conf IEEE Eng Med Biol Soc. (2019) 2019:4994–7. 10.1109/EMBC.2019.885763631946981 Wang W He L Ouyang C Chen C Xu X Ye X. Key common genes in obstructive sleep apnea and lung cancer are associated with prognosis of lung cancer patients. Int J Gen Med. (2021) 14:5381–96. 10.2147/IJGM.S33068134526807 Kang HS Kwon HY Kim IK Ban WH Kim SW Kang HH . Intermittent hypoxia exacerbates tumor progression in a mouse model of lung cancer. Sci Rep. (2020) 10:1854. 10.1038/s41598-020-58906-732024881 Pérez-Warnisher MT Cabezas E Troncoso MF Gómez T Melchor R Pinillos EJ . Sleep disordered breathing and nocturnal hypoxemia are very prevalent in a lung cancer screening population and may condition lung cancer screening findings: results of the prospective sleep apnea in lung cancer screening (SAILS) study. Sleep Med. (2019) 54:181–6. 10.1016/j.sleep.2018.10.02030580192 Breiman L. Random forests, machine learning. J Clin Micorbiol. (2001) 45:199–228. 10.1023/A:1010933404324 Vasquez MM Hu C Roe DJ Chen Z Halonen M Guerra S. Least absolute shrinkage and selection operator type methods for the identification of serum biomarkers of overweight and obesity: simulation and application. BMC Med Res Methodol. (2016) 16:154. 10.1186/s12874-016-0254-827842498 Debevec T Millet GP Pialoux V. Hypoxia-induced oxidative stress modulation with physical activity. Front Physiol. (2017) 8:84. 10.3389/fphys.2017.0008428243207 Younes M. Role of respiratory control mechanisms in the pathogenesis of obstructive sleep disorders. J Appl Physiol. (1985) 105:1389–405. 10.1152/japplphysiol.90408.200818787092 Maniaci A Iannella G Cocuzza S Vicini C Magliulo G Ferlito S . Oxidative stress and inflammation biomarker expression in obstructive sleep apnea patients. J Clin Med. (2021) 10:277. 10.3390/jcm1002027733451164 Komiyama Y Nakae S Matsuki T Nambu A Ishigame H Kakuta S . IL-17 plays an important role in the development of experimental autoimmune encephalomyelitis. J Immunol. (2006) 177:566–73. 10.4049/jimmunol.177.1.56635135958 Ahluwalia B Magnusson MK Öhman L. Mucosal immune system of the gastrointestinal tract: maintaining balance between the good and the bad. Scand J Gastroenterol. (2017) 52:1185–93. 10.1080/00365521.2017.134917328697651 Sade K Fishman G Kivity S DeRowe A Langier S. Expression of Th17 and Treg lymphocyte subsets in hypertrophied adenoids of children and its clinical significance. Immunol Invest. (2011) 40:657–66. 10.3109/08820139.2011.57542621542720 Liu Y Gu YC Tang Y Lü XL Chen H Wu QZ. (Role and mechanism of Th17/Treg in obstructive sleep apnea hypopnea syndrome). Zhonghua Yi Xue Za Zhi. (2016) 96:3577–81. 10.3760/cma.j.issn.0376-2491.2016.44.00927916079 Bhushan B Guleria R Misra A Luthra K Vikram NK. TNF-alpha gene polymorphism and TNF-alpha levels in obese Asian Indians with obstructive sleep apnea. Respir Med. (2009) 103:386–92. 10.1016/j.rmed.2008.10.00119022640 de Lima FF Mazzotti DR Tufik S Bittencourt L. The role inflammatory response genes in obstructive sleep apnea syndrome: a review. Sleep Breath. (2016) 20:331–8. 10.1007/s11325-015-1226-726201496 Varvarigou V Dahabreh IJ Malhotra A Kales SN. A review of genetic association studies of obstructive sleep apnea: field synopsis and meta-analysis. Sleep. (2011) 34:1461–8. 10.5665/sleep.137622043116 Wang L Liang Q Lin A Chen X Wu Y Zhang B . Puerarin increases survival and protects against organ injury by suppressing NF-κB/JNK signaling in experimental sepsis. Front Pharmacol. (2020) 11:560. 10.3389/fphar.2020.0056032457606 Ryan S. Adipose tissue inflammation by intermittent hypoxia: mechanistic link between obstructive sleep apnoea and metabolic dysfunction. J Physiol. (2017) 595:2423–30. 10.1113/JP27331227901270 Kent BD Ryan S McNicholas WT. Obstructive sleep apnea and inflammation: relationship to cardiovascular co-morbidity. Respir Physiol Neurobiol. (2011) 178:475–81. 10.1016/j.resp.2011.03.01521439407 Wang B Li W Jin H Nie X Shen H Li E . Curcumin attenuates chronic intermittent hypoxia-induced brain injuries by inhibiting AQP4 and p38 MAPK pathway. Respir Physiol Neurobiol. (2018) 255:50–7. 10.1016/j.resp.2018.05.00629758366 Craig EA Stevens MV Vaillancourt RR Camenisch TD. MAP3Ks as central regulators of cell fate during development. Dev Dyn. (2008) 237:3102–14. 10.1002/dvdy.2175018855897 Wali SO Al-Mughales J Alhejaili F Manzar MD Alsallum F Almojaddidi H . The utility of proinflammatory markers in patients with obstructive sleep apnea. Sleep Breath. (2021) 25:545–53. 10.1007/s11325-020-02149-332705528 Khan J Alghamdi H Anwer MM Eliav E Ziccardi V. Role of collagen conduit with duloxetine and/or pregabalin in the management of partial peripheral nerve injury. J Oral Maxillofac Surg. (2016) 74:1120–30. 10.1016/j.joms.2016.01.02226878362 Qiao L Li Y Sun S. Insulin exacerbates inflammation in fibroblast-like synoviocytes. Inflammation. (2020) 43:916–36. 10.1007/s10753-020-01178-031981062 Lakkireddy M Gadiga SG Malathi RD Karra ML Raju ISSVPM Ragini . Impact of daily high dose oral vitamin D therapy on the inflammatory markers in patients with COVID 19 disease. Sci Rep. (2021) 11:10641. 10.1038/s41598-021-97181-y34462608 Zaiss D Gause WC Osborne LC Artis D. Emerging functions of amphiregulin in orchestrating immunity, inflammation, and tissue repair. Immunity. (2015) 42:216–26. 10.1016/j.immuni.2015.01.02025692699 Berasain C Avila MA. Amphiregulin. Semin Cell Dev Biol. (2014) 28:31–41. 10.1016/j.semcdb.2014.01.00524463227 O'Reilly SM Leonard MO Kieran N Comerford KM Cummins E Pouliot M . Hypoxia induces epithelial amphiregulin gene expression in a CREB-dependent manner. Am J Physiol Cell Physiol. (2006) 290:C592–600. 10.1152/ajpcell.00278.200516207795 Kyotani Y Ota H Itaya-Hironaka A Yamauchi A Sakuramoto-Tsuchida S Zhao J . Intermittent hypoxia induces the proliferation of rat vascular smooth muscle cell with the increases in epidermal growth factor family and erbB2 receptor. Exp Cell Res. (2013) 319:3042–50. 10.1016/j.yexcr.2013.08.01423968588 Jang MK Son Y Jung MH. ATF3 plays a role in adipocyte hypoxia-mediated mitochondria dysfunction in obesity. Biochem Biophys Res Commun. (2013) 431:421–7. 10.1016/j.bbrc.2012.12.15423333392 Chen SC Liu YC Shyu KG Wang DL. Acute hypoxia to endothelial cells induces activating transcription factor 3 (ATF3) expression that is mediated via nitric oxide. Atherosclerosis. (2008) 201:281–8. 10.1016/j.atherosclerosis.2008.02.01418377912 Hai T Wolford CC Chang YS. ATF3, a hub of the cellular adaptive-response network, in the pathogenesis of diseases: is modulation of inflammation a unifying component? Gene Expr. (2010) 15:1–11. 10.3727/105221610X1281968655501521061913 Shi Z Xu L Xie H Ouyang R Ke Y Zhou R . Attenuation of intermittent hypoxia-induced apoptosis and fibrosis in pulmonary tissues via suppression of ER stress activation. BMC Pulm Med. (2020) 20:92. 10.1186/s12890-020-1123-032299413 Makita S Takatori H Nakajima H. Post-transcriptional regulation of immune responses and inflammatory diseases by RNA-Binding ZFP36 family proteins. Front Immunol. (2021) 12:711633. 10.3389/fimmu.2021.71163334276705 Racanelli AC Kikkers SA Choi A Cloonan SM. Autophagy and inflammation in chronic respiratory disease. Autophagy. (2018) 14:221–32. 10.1080/15548627.2017.138982329130366 Zhang Y Li NF Abulikemu S Zhang DL Wang YC Kong JQ . Relationship between zinc finger protein 36 (ZFP36) gene polymorphisms and obstructive sleep apnea. Genet Mol Res. (2015) 14:6733–43. 10.4238/2015.June.18.1726125882 Hoffmann MS Singh P Wolk R Narkiewicz K Somers VK. Obstructive sleep apnea and intermittent hypoxia increase expression of dual specificity phosphatase 1. Atheroscerosis. (2013) 231:378–83. 10.1016/j.atherosclerosis.2013.09.03324267255 Ryan S Arnaud C Fitzpatrick SF Gaucher J Tamisier R Pépin JL. Adipose tissue as a key player in obstructive sleep apnoea. Eur Respir Rev. (2019) 28:190006. 10.1183/16000617.0006-201931243096 Xu X Xu J. Effects of different obesity-related adipokines on the occurrence of obstructive sleep apnea. Endocr J. (2020) 67:485–500. 10.1507/endocrj.EJ20-003632307359 Eiselein L Nyunt T Lamé MW Ng KF Wilson DW Rutledge JC . TGRL lipolysis products induce stress protein atf3 via the TGF-β receptor pathway in human aortic endothelial cells. PLoS ONE. (2015) 10:e145523. 10.1371/journal.pone.014552326709509 Cherniack EP Cherniack NS. Obstructive sleep apnea, metabolic syndrome, and age: will geriatricians be caught asleep on the job? Aging Clin Exp Res. (2010) 22:1–7. 10.1007/BF0332480820305362 Wang L Wang L Zhang H Lu J Zhang Z Wu H . AREG mediates the epithelial-mesenchymal transition in pancreatic cancer cells via the EGFR/ERK/NF-κB signalling pathway. Oncol Rep. (2020) 43:1558–68. 10.3892/or.2020.752332323797 Avraham S Korin B Aviram S Shechter D Shaked Y Aronheim A. ATF3 and JDP2 deficiency in cancer associated fibroblasts promotes tumor growth via SDF-1 transcription. Oncogene. (2019) 38:3812–23. 10.1038/s41388-019-0692-y30670778 Teng F Xu Z Chen J Zheng G Zheng G Lv H . DUSP1 induces apatinib resistance by activating the MAPK pathway in gastric cancer. Oncol Rep. (2018) 40:1203–22. 10.3892/or.2018.652029956792 Xia W Liu Y Du Y Cheng T Hu X Li X. MicroRNA-423 drug resistance and proliferation of breast cancer cells by targeting ZFP36. Onco Targets Ther. (2020) 13:769–82. 10.2147/OTT.S21774532158228 Marrone O Bonsignore MR. Obstructive sleep apnea and cancer: a complex relationship. Curr Opin Pulm Med. (2020) 26:657–67. 10.1097/MCP.000000000000072932925366