Identification of novel biomarkers in obstructive sleep apnea via integrated bioinformatics analysis and experimental validation

Abstract

Background Obstructive sleep apnea (OSA) is a complex and multi-gene inherited disease caused by both genetic and environmental factors. However, due to the high cost of diagnosis and complex operation, its clinical application is limited. This study aims to explore potential target genes associated with OSA and establish a corresponding diagnostic model. Methods This study used microarray datasets from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) related to OSA and perform functional annotation and pathway analysis. The study employed multi-scale embedded gene co-expression network analysis (MEGENA) combined with least absolute shrinkage and selection operator (LASSO) regression analysis to select hub genes and construct a diagnostic model for OSA. In addition, the study conducted correlation analysis between hub genes and OSA-related genes, immunoinfiltration, gene set enrichment analysis (GSEA), miRNA network analysis, and identified potential transcription factors (TFs) and targeted drugs for hub genes. Finally, the study used chronic intermittent hypoxia (CIH) mouse model to simulate OSA hypoxic conditions and verify the expression of hub genes in CIH mice. Results In this study, a total of 401 upregulated genes and 275 downregulated genes were identified, and enrichment analysis revealed that these differentially expressed genes may be associated with pathways such as vasculature development, cellular response to cytokine stimulus, and negative regulation of cell population proliferation. Through MEGENA combined with LASSO regression, seven OSA hub genes were identified, including C12orf54, FOS, GPR1, OR9A4, MYO5B, RAB39B, and KLHL4. The diagnostic model constructed based on these genes showed strong stability. The expression levels of hub genes were significantly correlated with the expression levels of OSA-related genes and mainly acted on pathways such as the JAK/STAT signaling pathway and the cytosolic DNA-sensing pathway. Drug-target predictions for hub genes were made using the Connectivity Map (CMap) database and the Drug-Gene Interaction database (Dgidb), which identified targeted therapeutic drugs for the hub genes. In vivo experiments showed that the hub genes were all decreasing in the OSA mouse model. Conclusions This study identified novel biomarkers for OSA and established a reliable diagnostic model. The transcriptional changes identified may help to reveal the pathogenesis, mechanisms, and sequelae of OSA.