MiRNA expression profiles in healthy OSAHS and OSAHS with arterial hypertension: potential diagnostic and early warning markers
Respiratory Research volume 19, Article number: 194 (2018)
Obstructive sleep apnea-hypopnea syndrome (OSAHS) is prone to being complicated with various cardiovascular, cerebrovascular and metabolic conditions. OSAHS, due to its multifactorial nature, entails individualized and comprehensive treatment. So far, no well-established diagnostic criteria for the disease are available. In recent years, miRNA has been shown to be a sensitive biomarker suggestive of the progression and prognosis of many diseases. In this study, we examined some serum miRNAs in healthy OSAHS (OSAHS patients without complication) and OSAHS with arterial hypertension, with an attempt to understand the potential effects on the disease, improve the diagnosis of OSAHS and find OSAHS-related diagnostic markers.
Against various diagnostic criteria, participants were divided into three groups: healthy OSAHS, OSAHS with arterial hypertension and healthy controls. Their serum miRNA profiles were assessed by microarray technology, and then differentially expressed miRNAs were verified by quantitative real-time PCR (qRT-PCR). The receiver operating characteristic (ROC) curves of miRNAs were constructed and the areas under the curve (AUC) were calculated. Meanwhile, the miRNAs were subjected to logistic regression analysis. The target genes were bioinformatically assessed, their functions and signaling pathways further determined and eventually an miRNA-gene network was established.
Analysis with the miRNA array exhibited that, compared with the control group, 12 differentially expressed miRNAs were found in healthy OSAHS, and 33 were found in OSAHS with arterial hypertension. The expression of miR-126-3p, let-7d-5p, miR-7641 and miR-1233-5p, miR-320b, miR-145-5p, miR-107, miR-26a-5p were validated by using qRT-PCR. Bioinformatics analysis predicted that the potential target genes of these miRNAs might be involved in metabolism, and the regulation of endothelial cells and nervous system. Moreover, the ROC analysis showed that the using miR-145-5p and let-7d-5p in combination can identify the healthy OSAHS, presence of miR-126-3p, miR-26a-5p and miR-107 was well indicative of OSAHS with arterial hypertension.
A cluster of dysregulation miRNAs have been found to be involved in the development of OSAHS patients. Moreover, these miRNAs might be used to be potential diagnostic and early warning markers.
OSAHS is a complex, multifactorial disorder known to affect, with varying degrees, millions of people worldwide . The clinical manifestations and therapeutic effects vary considerably with the stage of the disease. Therefore, early diagnosis and treatment are crucial to treatment efficacy and the prevention of complications. The apnea hypopnea index (AHI), the sum of apneas and hypopneas per hour of sleep, has been extensively used as a criterion for the diagnosis and assessment of OSAHS . Nonetheless, the incidence of OSAHS-related complications tend to vary substantially even in patients with similar AHI. Obviously, efficacy of the treatment vary with patients even with similar AHI. In fact, AHI, as a single parameter, is not a reliable measure of the severity of OSAHS [3,4,5]. Ideally, any treatment decision about OSAHS should be based on stratification of OSAHS symptoms. This “ideal OSAHS decision” should not be made against AHI alone. MiRNAs serve as important regulators in multicellular organisms and can regulate protein expression by blocking mRNA translation . By coordinating the expression of multiple genes, miRNAs establish a broad network that regulates various biological processes, such as cellular differentiation, proliferation, apoptosis and metabolism. In humans, miRNAs are ubiquitous in tissues or body fluids and are stably expressed in blood [7, 8]. Changes in miRNAs may precede the appearance of physical symptoms of certain diseases [9, 10]. The altered expression of miRNAs have been found in a wide array of pathologies, including cancers, diabetes, hypoxia and hypertension [11,12,13,14]. Many studies have demonstrated that miRNAs played pivotal parts in the maintenance of metabolic homeostasis and stable blood pressure [15, 16]. Therefore, in-depth study on their roles can help us better understand the mechanisms of diseases and provide new targets for their treatment.
As we know, OSAHS involves multiple organs and systems, such as cardiovascular systems. About one-third of hypertensive patients have OSAHS, and about half of OSAHS patients have hypertension . A multitude of studies have demonstrated the association between OSAHS and hypertension, especially the resistant hypertension [18, 19]. Individuals suffering from both OSAHS and hypertension tend to have increased cardiovascular risk.
In this study, we examined some diagnostic evidence and identified a set of warning signs for the progression of the disease and compared healthy OSAHS and OSAHS with arterial hypertension, with an attempt to provide diagnostic evidence and early warning indicators for a better diagnosis of OSAHS.
This study was approved by the Ethics Committee of the Union Hospital of Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China, and written informed consent was obtained from all participants.
The subjects were those who visited the Department of Otorhinolaryngology, the Union Hospital of Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China, due to snoring or apnea from March 2017 to November 2017. Healthy controls were recruited from the Health Screening Center of the same hospital and didn’t have any medical diseases or snoring.
All participants received overnight polysomnography (PSG), OSAHS patients and normal subjects were confirmed on the basis of American Academy of Sleep Medicine (AASM) Guidelines. Meanwhile, all participants underwent physical checkups, blood pressure measurement and tests of blood biochemistry.
The criteria for exclusion included: (1) presence of other sleep disorders (International Classification of Sleep Disorders (ICSD-II) diagnosis) and history of having received treatment for sleep related breathing disorders; (2) other diseases such as central nervous system diseases, cardiopathy, diabetes, renal disease, thyroid disease, cancer, ongoing infections; (3) old age (> 60 year-old), young age (< 24 year-old), or being pregnant.
They were divided into 3 groups, each group having 20 subjects: (1) healthy controls (AHI < 5, without hypertension, hyperglycemia or dyslipidemia), (2) healthy OSAHS (AHI > 5, without hypertension, hyperglycemia or dyslipidemia) and (3) OSAHS with arterial hypertension (AHI > 5, with hypertension without hyperglycemia or dyslipidemia). In each group, 5 subjects were subjected to microarray and the other 15 subjects received qRT-PCR for validation.
Blood sample collection
Blood samples were collected from all the participants at AM 8:00 after receiving full-night in-laboratory polysomnography and an overnight fasting. For miRNA studies, blood was harvested into an inert separation gel vacuum procoagulant tube. Within 60 min of blood collection, the blood was centrifuged at 3000 g for 10 mins to obtain serum. The supernatant was transferred to RNase-free Eppendorf tubes and then stored at − 80 °C until RNA extraction.
MiRNA isolation and expression analysis by microarray
Total RNA, including miRNA, was extracted from 15 frozen serum by using QIAGEN RNeasy Mini Kit (217,004; QIAGEN, Germany) according to the manufacturer’s instructions. RNA concentration was determined on an ND-2000 spectrophotometer (Nanodrop™, Thermo Fisher Scientific, USA) at 260 nm, 280 nm and RNA integrity was determined by RNA gel electrophoresis. Small isolated RNAs was processed with the FlashTag®Biotin HSR RNA labeling kit (Affymetrix, USA) and was subsequently subjected to hybridization in the GeneChip® Hybridization Oven 645 (Affymetrix, USA). Washing and scanning were respectively conducted on the GeneChip® Fluidics Station 450 (Affymetrix, USA) and the GeneChip® Scanner 3000 7G (Affymetrix, USA) by following instructions. The data were analyzed against Affymetrix Expression Console-1_4_0 (EC 5.0) by using Robust Multichip Average (RMA) as normalization method.
qRT-PCR verification of miRNA array results
Of the significantly differentially expressed miRNAs (there was at least a two-fold level change between the two groups and p value < 0.05) as found by the microarray, the miRNAs found in both healthy OSAHS and in OSAHS with arterial hypertension and in those that have been reportedly associated with mental disorder, endothelial dysfunction, oxidative stress or angiogenesis, were selected. Those miRNAs were verified in 45 blood samples by qRT-PCR. Total RNA was extracted according to the instructions of miRNeasy serum kit (217,184; QIAGEN, Germany). Moreover, due to the unavailability of stable internal control in the serum, C. elegans miR-39 miRNA mimic, the miRNeasy serum/plasma spike-in control (219,610; QIAGEN, Germany), was used to monitor miRNA purification and amplification. cDNA was synthesized from total RNA using the miScript II RT kit (218,161; QIAGEN, Germany) and qRT-PCR was performed in a Lightcycle 480 RT-PCR system (Roche Diagnostics Ltd., Rotkreuz, Switzerland) using the miScript SYBR Green PCR kit (218,073; QIAGEN, Germany). The thermal cycle profile was as follows: an initial activation of 15 min at 95 °C, and 45 cycles of denaturation (15 s at 94 °C), annealing (30 s at 55 °C) and extension (30 s at 70 °C). The specificity and identity of PCR products were evaluated by melting curve analysis, with cel-miR-39 used as a normalization standard for miRNAs. The relative levels of each miRNAs were determined by 2−[Ct(miRNA)-Ct(cel-miR-39)]. Primer sequences are listed in Table 1.
Bioinformatic analysis of proven miRNAs
We identified the biological processes that might be affected by these differentially expressed miRNAs by employing genomic enrichment analysis of regulated target genes. On the basis of the two groups of miRNAs obtained from the analysis, we used the miRanda and TargetScan databases to predict the target genes of the differentially expressed miRNAs. Next, the function of the differentially expressed genes were analyzed by using GO analysis . Meanwhile, KEGG pathways analysis was applied to predict the significant signaling pathways of the target genes, according to the microarray data [21, 22]. Then, we obtained the intersection target genes identified by GO and KEGG prediction. Besides, we constructed an miRNA-gene-network according to the targeting relationship between the miRNAs and these genes.
Normally distributed data were presented as means ± SEM, non-normally distributed data were expressed as medians (the first and third quartile values), and qualitative variables were given as frequencies and percentages. Possible differences in demographic and clinical variables among the groups were evaluated by one-way analysis of variance (ANOVA) and Kruskal-Wallis. For microarray analysis, we adopted the Random variance model modified t-test to filter the differentially expressed miRNAs between patients and healthy controls. We distinguished the different miRNAs by using the significance analysis and FDR analysis. Fisher’s exact test and χ2 test were employed to classify the GO categories and signaling pathways, and the false discovery rate (FDR) was calculated to calibrate the p value. For the results of qRT-PCR, independent sample Student t-test and Mann-Whitney U tests were performed for statistical analysis. To further determine the prediction values of the selected miRNAs, ROC curves were constructed among the three groups. The AUC was calculated, and the optimum sensitivity and specificity were determined by Youden index. Logistic regression models were generated to assess the value of the differentially expressed miRNAs in combination in the identification of OSAHS patients. Variables were included in the model if they exhibited statistically significant contributions as shown by the likelihood ratio test. The AUC was used to determine the maximum differentiating ability of the model. Variables which had no statistically significant contributions to the model but could increase its AUC were also selected. All analyses were performed using GraphPad Prism (version 5.0; GraphPad Software) and statistical software package SPSS (version 22.0, SPSS Inc., Chicago, IL, USA). All tests were two-sided and a p value < 0.05 was considered to be statistically significant.
Demographic data and clinical characteristics of participants
All participants enrolled were Chinese. Table 2 details the clinical characteristics, covering body mass index (BMI), PSG, and biochemical data. The age, BMI, gender, smoking and drinking were all matched among the three subgroups.
Differential expression of miRNAs
The microarray revealed that miRNAs experienced changes in healthy OSAHS and their counterparts with arterial hypertension as compared to controls. Compared to the healthy controls, 12 miRNAs were deregulated (P < 0.05) in healthy OSAHS. Among them, 6 were up-regulated and 6 were down-regulated (P < 0.05) (Fig. 1a). In OSAHS with arterial hypertension, 33 miRNAs were deregulated (P < 0.05), compared with the healthy control group. Of them, 20 were up-regulated and 13 down-regulated (P < 0.05) (Fig. 1b). Of note, 4 miRNAs were similarly deregulated in healthy OSAHS and OSAHS with arterial hypertension (Fig. 1c). The results were submitted to NCBI’s Gene Expression Omnibus (GEO).
qRT-PCR identification of differentially regulated miRNAs
As aforementioned, we chose 6 miRNAs for qRT-PCR validation. Of them, 4 were overlapping miRNAs from healthy OSAHS and OSAHS with arterial hypertension (miR-126-3p, let-7d-5p, miR-7641 and miR-1233-5p), two were from OSAHS with arterial hypertension (miR-107, miR-26a-5p). Despite their low fold-change of 1.45 and 1.2, miR-145-5p and miR-320b were also chosen from healthy OSAHS, because, according to the selection criteria, only overlapping miRNAs would be selected in healthy OSAHS. Meanwhile, these two miRNAs had been considered to be potential biomarkers indicative of cardiovascular conditions and hypoxia. Thus, 8 miRNAs were selected for validation by qRT-PCR. Only miR-126-3p was verified to be down-regulated in both OSAHS groups (Fig. 2a, b) (P < 0.05). Moreover, let-7d-5p was down-regulated only in healthy OSAHS (Fig. 2c, d) (P < 0.05). Significant down-regulation of other 4 miRNAs (miR-320b, miR-145-5p, miR-107, miR-26a-5p) was also confirmed (Fig. 2e, f, g, h) (P < 0.05). The expression levels of miR-7641 and miR-1233-5p were too low to be appreciable in the three groups (data not show).
Microarray-based bioinformatical prediction
From the miRanda and TargetScan databases, we retrieved 959 mRNAs corresponding to miRNAs we verified in healthy OSAHS and 537 mRNAs corresponding to miRNAs we verified in OSAHS with arterial hypertension. Next, these target genes were subjected to GO analysis and KEGG pathway analysis. In the healthy OSAHS group, GO analysis revealed these target genes involved 26 GOs, which were related to nervous systems, cardiovascular system, cancer and signal transduction (P < 0.01) (Fig. 3a). At the same time, 28 GOs were demonstrated to be regulated in the OSAHS with arterial hypertension (P < 0.01) (Fig. 3b). These GOs were implicated in metabolism, hormone secretion, vascular endothelial growth, signal transduction. Similarly, the pathway analysis showed that 25 pathways were regulated in the healthy OSAHS and 28 pathways in the OSAHS with arterial hypertension (P < 0.05) (Fig. 4a, b). Furthermore, the miRNA-gene-network consisting of the core miRNAs and key target genes revealed the relationship between miRNA and target genes (Fig. 5a, b).
Diagnostic accuracy of the selected miRNAs
The sensitivity and specificity of these serum miRNAs as potential diagnostic and early warning markers in OSAHS patients were investigated by ROC analysis. The AUC of miR-126-3p in the two groups was 0.822 (95% confidence interval (CI) = 0.671–0.973; sensitivity: 86.7%; specificity: 66.7%; P = 0.0027), 0.871 (95% CI = 0.747–0.995; sensitivity: 86.7%; specificity: 73.3%; P < 0.001) (Fig. 6a). The AUC for let-7d-5p was 0.778 (95% CI = 0.598–0.958; sensitivity: 66.7%; specificity: 86.7%; P < 0.001) (Fig. 6b). The AUC for miR-320b was 0.751 (95% CI = 0.570–0.932; sensitivity: 86.7%; specificity: 66.7%; P = 0.019) (Fig. 6c). The AUC for miR-145-5p was 0.849 (95% CI = 0.711–0.986; sensitivity: 66.7%; specificity: 93.3%; P = 0.0011) (Fig. 6d). The AUC for miR-107 was 0.733 (95% CI = 0.552–0.915; sensitivity: 86.7%; specificity: 60%; P = 0.029) (Fig. 6e). The AUC for miR-26a-5p was 0.729 (95% CI = 0.543–0.914; sensitivity: 100%; specificity: 53.3%; P = 0.032) (Fig. 6f). Logistic regression analysis was used to evaluate the differentiating ability of the miRNAs in combination. A first multivariate model (model 1) (Table 3) was employed to identify the healthy OSAHS by using the expression levels of both miR-145-5p and let-7d-5p. Both of these miRNAs contributed significantly to the model as shown by the likelihood ratio test. Model 1 exhibited good differentiating ability (AUC = 0.951; 95% CI: 0.856 to 1.00; sensitivity: 93%; specificity: 100%; P < 0.001) (Fig. 7a). Model 2 only included miR-126-3p, given that it was not a miRNA specific to OSAHS with arterial hypertension, and among the subset of differentially expressed miRNAs, miR-26a-5p and miR-107 could increase the AUC of model, they were included in the final model (model 3) (Table 4). Figure 7b shows the ROC of the model 3 (AUC = 0.969; 95% CI: 0.915 to 1.00; sensitivity: 86.7%; specificity: 100%; P < 0.001). These models were shown to be able to well differentiate OSAHS patients from healthy controls.
As a serious poorly-understood clinical condition, OSAHS is associated with many diurnal symptoms and has been a subject of active studies [23,24,25]. So far, progress has been limited in spite of the continuous research efforts . To find better diagnostic and therapeutic strategies, we need to fully understand its underlying mechanisms. Serum miRNA, due to its unique pattern of the expression, can serve a fingerprint for the evaluation of prognosis of various diseases. Previous studies prompted us to speculate that, by identifying dysregulated miRNAs, it is possible to find specific biomarkers of different subsets of OSAHS. In this study, we tried to find representative expression signatures in peripheral blood that might be associated with OSAHS, with an attempt to identify the early warning signs and prognostic markers of healthy OSAHS and OSAHS with arterial hypertension. Moreover, we theorized that the combination of several biomarkers may help achieve more accurate assessment of OSAHS.
In this study, we analyzed the differences in serum miRNA expression in healthy controls, healthy OSAHS and OSAHS with arterial hypertension using microarray technique, qRT-PCR and statistical analysis. Our data showed six serum miRNAs had significant changes in OSAHS patients as compared to healthy controls. A single cluster of miRNAs appears to be able to specifically identify healthy OSAHS and OSAHS with arterial hypertension. The specificity and sensitivity of the six miRNAs were not good in predicting OSAHS. We used a logistic regression analysis to select those miRNAs that were associated with the two subtypes of OSAHS. The results showed that the combination of miR-145-5p and let-7d-5p and the combination of miR-126-3p, miR-26a-5p and miR-107 was more valuable, in terms of specificity and sensitivity, in differentiating the OSAHS subtypes.
Among the dysregulated miRNAs, the miR-126-3p was found to be consistently down-regulated in healthy OSAHS and in OSAHS with arterial hypertension, and in both types of OSAHS, miR-126-3p had good predictive value, indicating that it, as a marker, was involved in the pathogenesis of both healthy OSAHS and OSAHS with arterial hypertension. MiR-126, one of the most abundant miRNAs in endothelial cells , has been found to regulate multiple processes, such as promotion of endothelial cell differentiation and maturation in embryonic vasculogenesis, inhibition of angiogenesis and proliferation of mature endothelial cells [28,29,30]. A study showed that miR-126 was down-regulated in cells under hypoxic conditions both in vitro and in vivo . Meanwhile, mounting evidence indicates that endothelial function is impaired in OSAHS patients [32,33,34]. In addition, Kontaraki et al. found that expression of miR-126 was down-regulated in hypertensive patients compared to healthy controls . Endothelial dysfunction and vascular remodeling are responses to the organ damage under hypertensive and hypoxic states. As an endothelium-specific factor, miR-126 could protect organs or tissues from hypoxia/reoxygenation (H/R)-induced injury  and inflammation, by enhancing endothelial repair during vascular remodeling . Our results indicated that this miRNA might be involved in the specific regulation of endothelial function in OSAHS patients.
Interestingly, among specific miRNAs in the two groups, most of them were hypoxamiRs (hypoxia-associated miRNAs), a specific group of miRNAs that has been proved to be associated with hypoxia . In OSAHS with arterial hypertension, miR-26a-5p, miR-107 were confirmed to be down-regulated. Several mechanisms have been suggested for the pathogenesis of hypertension, including increased activity of the sympathetic nervous system, overactivation of the renin-angiotensin-aldosterone system, dysfunction of the vascular endothelium, impaired platelet function, thrombogenesis, vascular smooth muscle and cardiac hypertrophy, and altered angiogenesis [14, 39]. These two miRNAs were reportedly dysregulated by hypoxia and were involved in cardiovascular disease, suggesting that they are involved in the pathogenesis of OSAHS with arterial hypertension. In fact, miR-26 has been consistently reported to be a contributor to endothelial dysfunction and impaired angiogenesis. One study found that miR-26a might exert an anti-apoptotic effect in endothelium by inhibiting TRPC6-induced calcium overload, and its expression was reduced in ApoE−/− mice with atherosclerosis induced by high-fat diet . Moreover, in another similar study, inhibiting miR-26a expression could promote apoptosis of vascular smooth muscle cells . These studies demonstrated that, in endothelial cells, miR-26a possessed significant anti-apoptotic effect. Icli et al. showed that miR-26a acts as an anti-angiogenic factor, and the reduced miR-26a level could promote angiogenesis in endothelial cells . Both miR-26a and miR-107 contributed to cell cycle arrest during oxygen deprivation that resulted in impaired cell proliferation . Indeed, miR-107 blocked hypoxic signaling by suppressing HIF-1β expression, and mediated p53 regulation of hypoxic signaling and tumor angiogenesis . In particular, miR-107 was a lipid-modulated miRNA that plays an important role in metabolic diseases. MiR-107 has been shown as a key regulator of insulin sensitivity , and was an important modulator in hepatic lipid metabolism . A recent study exhibited that miR-107 could regulate circadian rhythm in cultured cells by targeting CLOCK genes. Moreover, the effect of miR-107 on lipid metabolism could be partially mediated by modifying the circadian system . In addition, the miR-107-mediated regulation of intestinal microbiota and proinflammatory cytokines has been considered to be important for the maintenance of intestinal homeostasis . Presumably, miR-107 might regulate a good many events involved in the development of OSAHS with arterial hypertension.
Meanwhile, let-7d-5p, miR-320b and miR-145-5p were validated to be the specific miRNAs in healthy OSAHS. It was reported that let-7 miRNA was one of the highly-expressed miRNA in the human brain . Specifically, let-7 miRNAs was the first known human miRNA to play a vital role in brain development and could regulate neuronal differentiation/maturation, and nerve degeneration/regeneration [50,51,52]. More importantly, let-7 family members respond to hypoxia in a cell-specific manner . Let-7d-5p was previously reported to be involved in the development of psychiatric diseases or degenerative disease of the central nervous system [53, 54]. Consistent with previous reports, our study showed that let-7d-5p miRNA was down-regulated in healthy OSAHS. Meanwhile, a large proportion of adults with OSAHS were found to suffer from one or more neurocognitive impairment(s), including excessive daytime sleepiness, fatigue, depressed mood, impaired memory, and/or poor concentration . Accordingly, decreased serum let-7d-5p in OSAHS patients might account for the hypoxia-induced neural injury.
Among the miRNAs investigated in this study, cardiovascular-related miRNAs are of great value in the research of OSAHS. It has been reported that miR-320b was down-regulated during acute myocardial infarction and was intimately related to plaque stability [56, 57]. Previous studies suggested that miR-320b participates in a wide array of biological activities, including ischemia-reperfusion injury, angiogenesis [58, 59]. On the other hand, miR-320b exerts a significant effect on endothelial cells. During ischemia, the down-regulated miR-320b level could increase the release of endothelial vasoactive factors, such as VEGF, ET-1 and FN , thereby increasing the risk of atherosclerosis and ischemic cerebrovascular diseases. The coincidence between our finding and previous results suggests that miR-320b may cause cardiovascular and cerebrovascular disease in patients with OSAHS through recurrent hypoxia and reoxygenation. A previous study suggested miR-145 might be a regulator in vascular smooth muscle cells (VSMCs) both in vivo and in vitro, and might be a major miRNA that efficiently drove VSMC differentiation from multipotent stem cells . A recent study exhibited that miR-145 was significantly suppressed and Smad3 was subsequently increased in hypoxia-treated VSMCs. Moreover, miR-145/Smad3 signaling pathway might promote OSAHS-induced aortic remodeling, which might be initiated by inflammation and oxidative stress . In addition, hypoxia was found to aggravate inflammatory response and miR-145-5p might play an anti-inflammatory role to protect cells from ischemic and hypoxic injury . It is well-known that, during sleep, intermittent hypoxia can cause sympathetic activation, evoke systemic inflammation, oxidative stresses and vascular dysfunction [64,65,66]. Meanwhile, we observed that the circulating miR-145-5p experienced significant changes, suggesting that healthy OSAHS are more vulnerable to cardiovascular diseases.
Microarray tends to yield a long string of target genes. GO and KEGG pathway analyses are generally employed to know their molecular functions and biological processes of these genes [21, 67, 68]. In this study, we used GO analysis and found that the miRNAs we singled out were correlated with nervous system differentiation/development, cancer development, metabolism, vascular regulation, and signal transduction. In addition, KEGG pathway analysis revealed that these miRNAs were involved in metabolism, circadian rhythm, receptor interaction, HIF-1 signaling pathway, neurotrophin signaling pathway, calcium signaling pathway. In previous studies conducted in OSAHS patients, anatomical and neuroimaging findings showed that the patients had nervous system damage. Moreover, hypoxia-inducible factor (HIF), as a key regulator of hypoxia response in cells and tissues, plays an indispensable role in hypoxia. Meanwhile, OSAHS was shown to be an independent risk factor for a variety of metabolic diseases. These may explain the link between the above-mentioned biological processes and OSAHS. Besides, according to the miRNA-gene-network, these genes targeted by the miRNAs were involved in hypoxia, nerve growth, metabolism, vascular endothelial differentiation. Especially, insulin-like growth factor-1/insulin-like growth factor 1 receptor (IGF1/IGF1R) were found to partake in nerve regeneration , regulation of the viability of schwann cells [69, 70] and the axon growth of motor neurons . Another key gene was the angiomotin (AMOT), which plays an important part in angiogenesis . More importantly, an OSAHS genome-wide gene expression array study showed that AMOT was involved in OSAHS-related excessive daytime sleepiness by regulating a variety of endothelial cell functions .
However, several limitations of the present study should be mentioned. First, the sample size of our study was relatively small. We sub-grouped OSAHS patients in terms of freedom of complications and presence of arterial hypertension only, which could substantially limit number of available subjects. Many adult OSAHS patients often had more than one complication, such as diabetes, hypertension, dyslipidemia or cardio-cerebrovascular diseases. On the other hand, absence of complications and the single complication in our series increased the power of our study. Numerous therapeutic studies have shown that early OSAHS damage were reversible . Study of these special types of OSAHS help us better understand the progression of the diseases. Another limitation was that in comparing controls and OSAHS with arterial hypertension, let-7d-5p was down-regulated in the training set, but was significantly up-regulated in the validation set. The results of qRT-PCR and microarray were obviously inconsistent. We believe this inconsistency had something to do the small sample size of microarray, indicating that the results of microarray still needed to be further validated in a larger samples. In addition, during the subsequent verification, the miRNAs we selected were interpreted primarily on the basis of microarray results and previous findings, a number of miRNAs with low expression levels might have been left out. Regardless of these limitations, our study provides direct evidence that dysregulated miRNAs are important contributors to OSAHS progression. Moreover, the target genes of these six miRNAs and the potential molecular mechanisms may provide useful information for the treatment of OSAHS.
Our study demonstrated that, the miRNA expression was dysregulated in OSAHS patients, and some miRNAs were differentially expressed in healthy OSAHS and OSAHS with arterial hypertension. These miRNAs might be used as biomarkers for stratifying OSAHS patients, achieving early diagnosis and improving diagnostic accuracy.
American Academy of Sleep Medicine
Apnea hypopnea index
One-way analysis of variance
Area under the curve
Body mass index
False discovery rate
International Classification of Sleep Disorders
Insulin-like growth factor-1/insulin-like growth factor 1 receptor
Obstructive sleep apnea-hypopnea syndrome
Quantitative real-time PCR
Receiver operating characteristic
Vascular smooth muscle cell
Peppard PE, Young T, Barnet JH, et al. Increased prevalence of sleep-disordered breathing in adults. Am J Epidemiol. 2013;177(9):1006–14.
The Report of an American Academy of Sleep Medicine Task Force. Sleep-related breathing disorders in adults: recommendations for syndrome definition and measurement techniques in clinical research. Sleep. 1999;22(5):667–89.
Redline S, Sanders M. Hypopnea, a floating metric: implications for prevalence, morbidity estimates, and case finding. Sleep. 1997;20(12):1209–17.
Stepnowsky CJ, Berry C, Dimsdale JE. The effect of measurement unreliability on sleep and respiratory variables. Sleep. 2004;27(5):990–5.
Ruehland WR, Rochford PD, O’Donoghue FJ, et al. The new AASM criteria for scoring hypopneas: impact on the apnea hypopnea index. Sleep. 2009;32(2):150–7.
Bartel DP. miRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
Mitchell PS, Parkin RK, Kroh EM, et al. Circulating miRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. 2008;105(30):10513–8.
Weber JA, Baxter DH, Zhang S, et al. The microRNA spectrum in 12 body fluids. Clin Chem. 2010;56(11):1733–41.
Urbich C, Kuehbacher A, Dimmeler S. Role of miRNAs in vascular diseases, inflammation, and angiogenesis. Cardiovasc Res. 2008;79(4):581–8.
Creemers EE, Tijsen AJ, Pinto YM. Circulating miRNAs: novel biomarkers and extracellular communicators in cardiovascular disease? Circ Res. 2012;110(3):483–95.
Chen X, Ba Y, Ma L, et al. Characterization of miRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18(10):997–1006.
Yang Z, Chen H, Si H, et al. Serum miR-23a, a potential biomarker for diagnosis of pre-diabetes and type 2 diabetes. Acta Diabetol. 2014;51(5):823–31.
Kulshreshtha R, Davuluri RV, Calin GA, et al. A microRNA component of the hypoxic response. Cell Death Differ. 2008;15(4):667–71.
Batkai S, Thum T. miRNAs in hypertension: mechanisms and therapeutic targets. Curr Hypertens Rep. 2012;14(1):79–87.
Karolina DS, Tavintharan S, Armugam A, et al. Circulating miRNA profiles in patients with metabolic syndrome. J Clin Endocrinol Metab. 2012;97(12):E2271–6.
Romaine SP, Charchar FJ, Samani NJ, et al. Circulating miRNAs and hypertension--from new insights into blood pressure regulation to biomarkers of cardiovascular risk. Curr Opin Pharmacol. 2016;27:1–7.
Somers VK, White DP, Amin R, et al. Sleep apnea and cardiovascular disease: an American Heart Association/american college of Cardiology Foundation scientific statement from the American Heart Association Council for high blood pressure research professional education committee, council on clinical cardiology, stroke council, and council on cardiovascular nursing. In collaboration with the National Heart, Lung, and Blood Institute National Center on sleep disorders research (National Institutes of Health). Circulation. 2008;118(10):1080–111.
Pedrosa RP, Drager LF, Gonzaga CC, et al. Obstructive sleep apnea: the most common secondary cause of hypertension associated with resistant hypertension. Hypertension. 2011;58(5):811–7.
Goncalves SC, Martinez D, Gus M, et al. Obstructive sleep apnea and resistant hypertension: a case-control study. Chest. 2007;132(6):1858–62.
Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.
Draghici S, Khatri P, Tarca AL, et al. A systems biology approach for pathway level analysis. Genome Res. 2007;17(10):1537–45.
Kanehisa M, Furumichi M, Tanabe M, et al. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.
Maspero C, Giannini L, Galbiati G, Rosso G, Farronato G. Obstructive sleep apnea syndrome: a literature review. Minerva Stomatol. 2015;64(2):97–109.
Jackson ML, Howard ME, Barnes M. Cognition and daytime functioning in sleep-related breathing disorders. Prog Brain Res. 2011;190:53–68.
Myers KA, Mrkobrada M, Simel DL. Does this patient have obstructive sleep apnea?: the rational clinical examination systematic review. JAMA. 2013;310(7):731–41.
Dempsey JA, Veasey SC, Morgan BJ, O'Donnell CP. Pathophysiology of sleep apnea. Physiol Rev. 2010;90(1):47–112.
Chistiakov DA, Orekhov AN, Bobryshev YV. The role of miR-126 in embryonic angiogenesis, adult vascular homeostasis, and vascular repair and its alterations in atherosclerotic disease. J Mol Cell Cardiol. 2016;97:47–55.
van Solingen C, Seghers L, Bijkerk R, et al. Antagomir-mediated silencing of endothelial cell specific microRNA-126 impairs ischemia-induced angiogenesis. J Cell Mol Med. 2009;13(8A):1577–85.
van Solingen C, Bijkerk R, de Boer HC, et al. The role of microRNA-126 in vascular homeostasis. Curr Vasc Pharmacol. 2015;13(3):341–51.
Fish JE, Santoro MM, Morton SU, et al. miR-126 regulates angiogenic signaling and vascular integrity. Dev Cell. 2008;15(2):272–84.
Ye P, Liu J, He F, et al. Hypoxia-induced deregulation of miR-126 and its regulative effect on VEGF and MMP-9 expression. Int J Med Sci. 2014;11(1):17–23.
Kato M, Roberts-Thomson P, Phillips BG, et al. Impairment of endothelium-dependent vasodilation of resistance vessels in patients with obstructive sleep apnea. Circulation. 2000;102(21):2607–10.
Budhiraja R, Parthasarathy S, Quan SF. Endothelial dysfunction in obstructive sleep apnea. J Clin Sleep Med. 2007;3(4):409–15.
Jelic S, Padeletti M, Kawut SM, et al. Inflammation, oxidative stress, and repair capacity of the vascular endothelium in obstructive sleep apnea. Circulation. 2008;117(17):2270–8.
Kontaraki JE, Marketou ME, Zacharis EA, et al. MicroRNA-9 and microRNA-126 expression levels in patients with essential hypertension: potential markers of target-organ damage. J Am Soc Hypertens. 2014;8(6):368–75.
Yang HH, Chen Y, Gao CY, et al. Protective effects of MicroRNA-126 on human cardiac microvascular endothelial cells against hypoxia/Reoxygenation-induced injury and inflammatory response by activating PI3K/Akt/eNOS signaling pathway. Cell Physiol Biochem. 2017;42(2):506–18.
Wei Y, Schober A, Weber C. Pathogenic arterial remodeling: the good and bad of miRNAs. Am J Physiol Heart Circ Physiol. 2013;304(8):H1050–9.
Azzouzi HE, Leptidis S, Doevendans PA, et al. HypoxamiRs: regulators of cardiac hypoxia and energy metabolism. Trends Endocrinol Metab. 2015;26(9):502–8.
Nadar SK, Tayebjee MH, Messerli F, et al. Target organ damage in hypertension: pathophysiology and implications for drug therapy. Curr Pharm Des. 2006;12(13):1581–92.
Zhang Y, Qin W, Zhang L, et al. MicroRNA-26a prevents endothelial cell apoptosis by directly targeting TRPC6 in the setting of atherosclerosis. Sci Rep. 2015;5:9401.
Leeper NJ, Raiesdana A, Kojima Y, et al. MicroRNA-26a is a novel regulator of vascular smooth muscle cell function. J Cell Physiol. 2011;226(4):1035–43.
Icli B, Wara AK, Moslehi J, et al. MicroRNA-26a regulates pathological and physiological angiogenesis by targeting BMP/SMAD1 signaling. Circ Res. 2013;113(11):1231–41.
Kulshreshtha R, Ferracin M, Wojcik SE, et al. A microRNA signature of hypoxia. Mol Cell Biol. 2007;27(5):1859–67.
Yamakuchi M, Lotterman CD, Bao C, et al. P53-induced microRNA-107 inhibits HIF-1 and tumor angiogenesis. Proc Natl Acad Sci U S A. 2010;107(14):6334–9.
Trajkovski M, Hausser J, Soutschek J, et al. miRNAs 103 and 107 regulate insulin sensitivity. Nature. 2011;474(7353):649–53.
Bhatia H, Pattnaik BR, Datta M. Inhibition of mitochondrial beta-oxidation by miR-107 promotes hepatic lipid accumulation and impairs glucose tolerance in vivo. Int J Obes. 2016;40(5):861–9.
Daimiel-Ruiz L, Klett-Mingo M, Konstantinidou V, et al. Dietary lipids modulate the expression of miR-107, an miRNA that regulates the circadian system. Mol Nutr Food Res. 2015;59(3):552–65.
Xue X, Cao AT, Cao X, et al. Downregulation of microRNA-107 in intestinal CD11c(+) myeloid cells in response to microbiota and proinflammatory cytokines increases IL-23p19 expression. Eur J Immunol. 2014;44(3):673–82.
Anacker AM, Beery AK. Life in groups: the roles of oxytocin in mammalian sociality. Front Behav Neurosci. 2013;7:185.
Roush S, Slack FJ. The let-7 family of miRNAs. Trends Cell Biol. 2008;18(10):505–16.
Shao NY, Hu HY, Yan Z, et al. Comprehensive survey of human brain microRNA by deep sequencing. BMC Genomics. 2010;11:409.
Leinders M, Doppler K, Klein T, et al. Increased cutaneous miR-let-7d expression correlates with small nerve fiber pathology in patients with fibromyalgia syndrome. Pain. 2016;157(11):2493–503.
Maffioletti E, Cattaneo A, Rosso G, et al. Peripheral whole blood microRNA alterations in major depression and bipolar disorder. J Affect Disord. 2016;200:250–8.
Tan L, Yu JT, Tan MS, et al. Genome-wide serum microRNA expression profiling identifies serum biomarkers for Alzheimer's disease. J Alzheimers Dis. 2014;40(4):1017–27.
Lal C, Strange C, Bachman D. Neurocognitive impairment in obstructive sleep apnea. Chest. 2012;141(6):1601–10.
Huang S, Chen M, Li L, et al. Circulating miRNAs and the occurrence of acute myocardial infarction in Chinese populations. Circ Cardiovasc Genet. 2014;7(2):189–98.
Gidlof O, van der Brug M, Ohman J, et al. Platelets activated during myocardial infarction release functional miRNA, which can be taken up by endothelial cells and regulate ICAM1 expression. Blood. 2013;121(19):3908–17 S1-S26.
Feng B, Chakrabarti S. miR-320 regulates glucose-induced gene expression in diabetes. ISRN Endocrinol. 2012;2012:549875.
Wu YY, Chen YL, Jao YC, et al. miR-320 regulates tumor angiogenesis driven by vascular endothelial cells in oral cancer by silencing neuropilin 1. Angiogenesis. 2014;17(1):247–60.
Zhang R, Qin Y, Zhu G, et al. Low serum miR-320b expression as a novel indicator of carotid atherosclerosis. J Clin Neurosci. 2016;33:252–8.
Cordes KR, Sheehy NT, White MP, et al. miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature. 2009;460(7256):705–10.
Yu C, Liu Y, Sun L, et al. Chronic obstructive sleep apnea promotes aortic remodeling in canines through miR-145/Smad3 signaling pathway. Oncotarget. 2017;8(23):37705–16.
Yuan M, Zhang L, You F, et al. MiR-145-5p regulates hypoxia-induced inflammatory response and apoptosis in cardiomyocytes by targeting CD40. Mol Cell Biochem. 2017;431(1–2):123–31.
Garpestad E, Parker JA, Katayama H, et al. Decrease in ventricular stroke volume at apnea termination is independent of oxygen desaturation. J Appl Physiol (1985). 1994;77(4):1602–8.
Gaisl T, Bratton DJ, Kohler M. The impact of obstructive sleep apnoea on the aorta. Eur Respir J. 2015;46(2):532–44.
Kohler M, Stradling JR. Mechanisms of vascular damage in obstructive sleep apnea. Nat Rev Cardiol. 2010;7(12):677–85.
Khatri P, Draghici S, Ostermeier GC, et al. Profiling gene expression using onto-express. Genomics. 2002;79(2):266–70.
Hashimoto K, Goto S, Kawano S, et al. KEGG as a glycome informatics resource. Glycobiology. 2006;16(5):63R–70R.
Lewis ME, Neff NT, Contreras PC, et al. Insulin-like growth factor-I: potential for treatment of motor neuronal disorders. Exp Neurol. 1993;124(1):73–88.
Syroid DE, Zorick TS, Arbet-Engels C, et al. A role for insulin-like growth factor-I in the regulation of Schwann cell survival. J Neurosci. 1999;19(6):2059–68.
Ozdinler PH, Macklis JD. IGF-I specifically enhances axon outgrowth of corticospinal motor neurons. Nat Neurosci. 2006;9(11):1371–81.
Bratt A, Birot O, Sinha I, et al. Angiomotin regulates endothelial cell-cell junctions and cell motility. J Biol Chem. 2005;280(41):34859–69.
Chen YC, Chen KD, Su MC, et al. Genome-wide gene expression array identifies novel genes related to disease severity and excessive daytime sleepiness in patients with obstructive sleep apnea. PLoS One. 2017;12(5):e176575.
Patil SP, Schneider H, Schwartz AR, et al. Adult obstructive sleep apnea: pathophysiology and diagnosis. Chest. 2007;132(1):325–37.
We are indebted to all the subjects who took part in this study, and the staff of the Sleep Center of Union Hospital of Tongji Medical College, Huazhong University of Science and Technology.
This work is supported by the National Natural Science Foundation of China (81570903).
Availability of data and materials
Microarray data can be found under the GEO accession number GSE112093 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE112093).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Union Hospital of Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China (2015-S1001), and written informed consent was obtained from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, X., Niu, X., Xiao, Y. et al. MiRNA expression profiles in healthy OSAHS and OSAHS with arterial hypertension: potential diagnostic and early warning markers. Respir Res 19, 194 (2018). https://doi.org/10.1186/s12931-018-0894-9