- 1Department of Radiology, The First Affiliated Hospital of Nanchang University, Jiangxi, China
- 2Big Data Center, The Second Affiliated Hospital of Nanchang University, Jiangxi, China
- 3Department of PET Center, The First Affiliated Hospital of Nanchang University, Jiangxi, China
Obstructive sleep apnea (OSA) is the most common sleep disorder worldwide. Previous studies have shown that OSA patients are often accompanied by cognitive function loss, and the underlying neurophysiological mechanism is still unclear. This study aimed to determine whether there are differences in regional homogeneity (Reho) and functional connectivity (FC) across the brain between OSA patients with MCI (OSA-MCI) and those without MCI (OSA-nMCI) and whether such differences can be used to distinguish the two groups. Resting state magnetic resonance data were collected from 48 OSA-MCI patients and 47 OSA-nMCI patients. The brain regions with significant differences in Reho and FC between the two groups were identified, and the Reho and FC features were combined with machine learning methods for classification. Compared with OSA-nMCI patients, OSA-MCI patients showed significantly lower Reho in bilateral lingual gyrus and left superior temporal gyrus. OSA-MCI patients also showed significantly lower FC between the bilateral lingual gyrus and bilateral cuneus, left superior temporal gyrus and left middle temporal gyrus, middle frontal gyrus, and bilateral posterior cingulate/calcarine/cerebellar anterior lobe. Based on Reho and FC features, logistic regression classification accuracy was 0.87; sensitivity, 0.70; specificity, 0.89; and area under the curve, 0.85. Correlation analysis showed that MoCA scale score in OSA patients was significant positive correlation sleep efficiency and negatively correlation with neck circumference. In conclusion, our results showed that the OSA-MCI group showed decreased Reho and FC in specific brain regions compared with the OSA-nMCI group, which may help to understand the underlying neuroimaging mechanism of OSA leading to cognitive dysfunction and may serve as a potential biomarker to distinguish whether OSA is accompanied by cognitive impairment.
Introduction
Obstructive sleep apnea (OSA) is a highly common sleep disorder characterized by repeated partial collapse and obstruction of the upper respiratory tract, resulting in intermittent hypoxia, hypercapnia, and sleep fragmentation. OSA affects up to more than 30% of older adults (Peppard et al., 2013), with an incidence rate of ~14% in men and 5% in women (Heinzer et al., 2016). Long-term intermittent hypoxia is not only associated with various complications (hypertension, cardiovascular damage, and chronic kidney disease; Yayan et al., 2017), but can also lead to concomitant cognitive dysfunction, anxiety, and even Alzheimer’s disease (Daulatzai, 2015). Studies have also shown that OSA is associated with a variety of cognitive impairments, including deficits in attention, memory, executive function, visuospatial function, and language ability (Olaithe et al., 2018). However, the neuropathological mechanisms underlying cognitive dysfunction in OSA are not fully understood.
Neuroimaging studies to explain these cognitive deficits have revealed functional and structural changes in multiple brain regions in OSA patients (Kumar et al., 2012; Xia et al., 2016). Previous structural magnetic resonance imaging (sMRI) research has shown decreased white matter integrity and volume reduction (Lee et al., 2019) and regional cortical thinning in OSA patients (Joo et al., 2013). A resting-state function MRI (rs-fMRI) study to evaluate brain activity in the spontaneous state also showed that OSA patients had regional homogeneity changes (Santarnecchi et al., 2013) and abnormal global and regional functional connectivity (Park et al., 2016).It involves the default mode network (DMN), saliency network, and central executive network (Prilipko et al., 2011; Khazaie et al., 2017; Chen L. et al., 2018). Changes in the function and structure of these areas are associated with cognition, emotion, autonomic nerves, and sensory control (Chen L. T. et al., 2018; Yu et al., 2019; Kong et al., 2021). However, most previous studies have been limited to comparisons between OSA patients and healthy individuals, and functional differences between OSA patients with and without mild cognitive impairment (MCI) have not been assessed. Thus, how to account for the heterogeneity of concomitant cognitive impairment in the OSA population and the possible mechanisms linking OSA to MCI need to be clarified.
Rs-fMRI is a non-invasive technique that has become a hotspot in neuroimaging research owing to its high spatial and temporal resolution. Regional homogeneity (Reho) is a data-driven local measure of spontaneous neural activity that does not require the onset time of stimulation or prior knowledge and has good test–retest reliability (Xia et al., 2018). Functional connectivity (FC) is an fMRI data processing method used to explore the functional interactions between anatomically separated brain regions (Luo et al., 2021). Reho and FC are considered to be complementary in examining the synchronization of local and remote brain activity (Sun et al., 2018),and the simple computations and reliable representations involved have led to their use in exploring brain diseases with underlying functional alterations, such as epilepsy (Liu et al., 2021), Alzheimer’s (Ibrahim et al., 2021), and depression (Chen et al., 2021). Biomarkers obtained by various resting-state MRI techniques [functional connectivity, Reho, amplitude of low-frequency fluctuation (ALFF)] are widely used in disease classification (Hojjati et al., 2017, 2019; Ju et al., 2019; Wang et al., 2021). Long et al. used ALFF, regional uniformity, and gray matter density as classification features to improve the efficiency of MCI diagnosis by distinguishing between MCI patients and healthy controls. The classification model obtained high accuracy (sensitivity, 93.1%; specificity, 100%, area under the curve, 0.97; Long et al., 2016).
This study aimed to determine whether there are differences in regional consistency and FC across the brain between OSA patients with MCI (OSA-MCIs) and those without MCI (OSA-nMCIs) and whether such differences can be used to distinguish the two groups. We hypothesized that compare with OSA-nMCI, OSA-MCI would showed abnormal spontaneous functional activity and connectivity change. To test this hypothesis, the voxel-level Reho method was used to explore local spontaneous brain activity in OSA-MCI patients, while exploring functional connectivity with the whole brain using significantly different Reho brain regions as seed points. Then, logistic regression was used to evaluate whether Reho and FC values could be used as neuroimaging markers to distinguish whether MCI was associated with OSA.
Materials and methods
As a descriptive cross-sectional study, we conducted an function analysis using Gpower software. The specific settings were: Effect size d = 0.8, α = 0.05, 1 − β = 0.8. The total sample size was required to be 52. In this study, we recruited a total of 95 OSA patients diagnosed in the sleep monitoring room of the Respiratory Department of the First Affiliated Hospital of Nanchang University between August 2017 and August 2022. Diagnosis was jointly determined by experienced respiratory physicians according to the American Society of Sleep Medicine 2017 Clinical Practice Guidelines for adult obstructive sleep apnea (Kapur et al., 2017). The inclusion criterion was male with apnea-hypopnea index (AHI) > 15/h. The exclusion criteria were as follows: (1) sleep disorders other than OSA (e.g., insomnia, somnolence, and circadian dysrhythmicity sleep disorder); (2) respiratory diseases, cardiovascular diseases, diabetes mellitus, hypothyroidism, central nervous system diseases, metabolic diseases, and tumors; (3) alcohol or illicit drug abuse or current use of psychotropic substances; (4) contraindications for MRI, such as claustrophobia and various internal stents; and (5) image artifacts (motion, metal). All participants were right-handed, native Chinese speakers, and aged 20–60 years.
This study was approved by the Medical Ethics Committee of the First Affiliated Hospital of Nanchang University [2020(94)] and was conducted according to the principles of the Declaration of Helsinki. Written informed consent was obtained from all participants.
Polysomnography and neuropsychological assessment
All subjects were instructed to abstain from alcohol, coffee, and hypnotics prior to polysomnography monitoring. All participants underwent an overnight polysomnography (from 10 p.m. to 6 a.m. the next day) using a Respironics LE series physiological monitoring system (Alice 5 LE, Respironics, Orlando, FL, United States). Polysomnography monitoring included standard electrocardiogram, electroophthalmogram, electromyogram, electrocardiogram, oral and nasal airflow, chest and abdomen breathing movement, and snoring. Blood oxygen saturation (SaO2), sleep latency, total sleep time, sleep efficiency, sleep stage, arousal and respiratory events were also recorded.
Obstructive apnea was described as a sustained 90% reduction in airflow for ≥10 s with significant dyspnea. Hypopnea was defined as a ≥30% decrease in airflow for ≥10 s accompanied by an oxygen saturation of ≥4%. AHI was determined as the sum of apnea and hypopnea events per hour during sleep and was classified as mild (5–15), moderate (15–30), or severe (≥30; Berry et al., 2017). Cognitive function was assessed using the 11-item Montreal Cognitive Assessment (MoCA) scale (Chinese version) with the area under the ROC (AUC) of 0.930 (95%CI: 0.894; 0.965), a sensitivity of 92%, and a specificity of 85% (Hu et al., 2013). All MoCA scale assessments were performed by one senior experienced psychologists. Briefly, the MoCA scale examines eight cognitive domains, including executive functioning, language, attention, computation, abstraction, naming, memory, and orientation. A total MoCA score of <26 indicates cognitive impairment (Chen et al., 2016).
MRI data acquistion
All patients underwent MRI using a 3.0 Tesla MRI scanner with an 8-channel phased array head coil (Siemens, Munich, Germany) at our hospital. Foam pads and earplugs were used to reduce head movement and noise. Before the scan, the patients were instructed to close their eyes, stay awake, and not engage in specific thinking. First, a routine MRI scan was performed. The conventional T1-weighted imaging parameters were as follows: repetition time (TR) = 250 ms, echo time (TE) = 2.46 ms, thickness = 5 mm, gap = 1.5 mm, field of view (FOV) = 230 × 230 mm. For T2-weighted imaging, they were: TR = 4,000 ms, TE = 113 ms, thickness = 5 mm, gap = 1.5 mm, FOV = 230 × 230 mm, slice = 19. Then, high-resolution T1-weighted MRI images of brain structures were obtained using sagittal brain volume sequences (TR = 1,900 ms, TE = 2.26 ms, thickness = 1.0 mm, gap = 0.5 mm, FOV = 250 × 250 mm, Matrix = 256 × 256, turning angle = 9°, slice = 176). Finally, the rs-fMRI data were collected in the echo plane imaging sequence (TR = 2,000 ms, TE = 30 ms, flip angle = 90°, thickness = 4.0 mm, gap = 1.2 mm, FOV = 230 × 230 mm2, matrix size = 64 × 64, slice = 30). A total of 240 rs-fMRI images were obtained. Two experienced radiologists read the images to exclude macroscopic lesions (demyelinating encephalopathy, brain tumors) and motion artifacts.
Data preprocess
Imaging data were examined using MRIcro software1 to discard suboptimal data, such as the presence of deflected or missing images, etc. Data processing is based on Data Processing & Analysis for Brain Imaging 6.0 (DPABI6.0, Chinese Academy of Sciences, Beijing, China, http://rfmri.org/dpabi), which is based on statistical parameter mapping (SPM12, http://www.fil.ion.ucl.ac.uk/spm/software/spm12/), and in MATLAB2018b (Math Works, Natick, MA, United States). Firstly, the file format was converted from Digital Imaging and Communications in Medicine to Neuroimaging Informatics Technology Initiative. Secondly, The first 10 time points are removed to ensure signal stability. Then, time layer correction and 3D head motion correction were performed. Patients whose frame displacement exceeded 2.5 standard deviations were excluded (Van Dijk et al., 2012). Structural images were co-registered with functional images for each subject using a linear transformation. Thus, the new segmentation in SPM12 was used to segment the structural images of all subjects into grey matter, white matter and CSF. Then, the image space was normalized to the Montreal Neurological Institute (MNI) template and resampled to 3 × 3 × 3 mm voxels. Finally, the Friston 24 parameters, white matter signals, and cerebrospinal fluid signals were regressed from the time series of all voxels using linear regression, after filtering with a temporal filter (0.01–0.08 Hz). The data preprocessing is described in detail in our previous study (Li et al., 2021).
The Reho value was generated using Kenndall’s consistency coefficient to generate a separate Reho map (in a given voxel and its 26 nearest time series; Zang et al., 2004). The calculation formula is as follows:
where Reho is the consistency coefficient of Kenndall in a given voxel, and the range is 0–1; K is the voxel number of the time series in the measurement cluster, k = 27; n is the rank number; and Ri is the total rank of the ith time point, where = (n + 1)K/2 is the average value of Ri. When a given cluster is consistent with its neighbors in the time series, the consistency coefficient is close to 1. Then, to reduce the global effect of subject variability, we divided the Reho plot by the Reho value of the global tie and spatially smoothed with a gaussian kernel (half-height full-width = 6 mm). Reho plots were generated using Fisher’s r-to-z normalization, and szReho plots were used for final statistical analyses.
Functional connectivity analysis
Functional connectivity analysis was calculated using DPABI 6.0. The brain areas that exhibited significantly different Reho values between the OSA-MCI and OSA-nMCI groups were set as regions of interest (ROIs). First, the time series of each region (i.e., the average of the fMRI sequences of all voxels in that region) were extracted. Then, Pearson’s correlation coefficients between the time series and all other voxels in the brain were calculated to obtain a whole-brain FC map of the subject. Finally, Fish R-to-z changes were used to transform the correlation maps into Z-value maps.
Logistic regression
Logistic regression (LR) is a multivariable technique that requires a functional relationship between a number of predictor variables and a single output to predict the likelihood of a target variable. Reho and FC values, which were significantly different between OSA-MCI and OSA-nMCI groups, were used as features to construct LR models. First, we calculated the Pearson correlation coefficient between the characteristics of each group, with a 0.75 cutoff for the related absolute threshold. Given that the characteristics of correlation is greater than the threshold after comparing the mean absolute correlation, the mean absolute correlation variable was deleted. To weaken the multicollinearity at the expense of the smaller variable, the most valuable predictor variable was retained. Mesh optimization was used to optimize the hyperparameters, and leave-one-out cross validation and permutation test (5,000 times) were used to verify the model performance. The accuracy, sensitivity and specificity of the model were calculated, and the predictive performance of the model was evaluated according to the receiver operating characteristic curve and the area under the curve. The heterogeneity of the test results was evaluated by Kappa.
Statistics
For demographic and clinical evaluation data, we first processed the data using SPSS 23.0 software and tested the normality of the data by Kolmogorov–Smirnov. Then, two-sample t-test was performed on the data conforming to normal distribution, and Mann–Whitney U test was performed on the data not normally distributed. p < 0.05 was considered statistically significant.
Voxel-based comparisons of all Reho maps and functional connections were performed using DPABI6.0. Two-sample t-tests were performed to analyze differences in the Reho and FC maps among the groups. Multiple comparisons were corrected using the Gaussian random field method (voxel level, p < 0.005; cluster level, p < 0.05).
Furthermore, Pearson correlation analysis was performed to evaluate the relationship between clinical variables, MoCA scale scores and the abnormal signal value in Reho and FC with significant group differences. All statistical analyses were performed using SPSS 23.0 with a statistical significance level of p < 0.05, and the analyses were corrected for multiple comparisons using the Bonferroni correction.
Results
Demographic and clinical characteristics
There were significant differences in neck circumference (p = 0.018) and MoCA scale score (p < 0.001) between the OSA-MCI and OSA-nMCI groups. Meanwhile, there were no between-group significant differences in age, education, BMI, waist circumference, AHI, LSaO2, MSaO2, SaO2 < 90%, sleep efficiency, oxygen index reduction, and AI (all p > 0.05; Table 1).
Group differences in Reho maps
The group means of the Reho plots for OSA-MCI and OSA-nMCI are shown in Figure 1. Two-sample t-test was performed to compare Reho data between OSA-MCI group and OSA-nMCI group. Differences between groups are shown in Table 2 and Figure 2. Compared with OSA-nMCI, OSA-MCI group showed significant abnormalities in bilateral lingual gyrus and left superior temporal gyrus.
Figure 2. Significantly different Reho was observed between the OSA-MCI and OSA-nMCI groups. Reho values decreased significantly in the MCI group compared to the nMCI group (blue). B, bilateral; L, left; STG, superior temporal gyrus.
Group differences in functional connectivity maps
Patients with OSA-MCI showed decreased FC between left superior temporal gyrus (STG) and left middle temporal gyrus (MTG), left middle frontal gyrus (MFG), bilateral posterior cingutate (PC) /Calcarine/left cerebellum anterior lobe (CAL), and decreased FC between bilateral lingual gyrus and bilateral Cuneus (Table 3; Figure 3).
Figure 3. Significantly different FC values seeded as ROI in B-Lingual (A) and L-STG (B). In figure (A), red represents the lingual gyrus and yellow represents the cuneus; In figure (B), red represents the left superior temporal gyrus, green represents the left middle frontal gyrus, blue represents the left middle temporal gyrus, and yellow represents the bilateral posterior cingulate/calcarine/left cerebellum anterior lobe. The FCs between significantly different regions and ROIs in OSA-MCI group was significantly reduced compare with OSA-nMCI group. L, left; B, bilateral; ROI, region of interest; FC, functional connectivity; STG, superior temporal gyrus; MTG, middle temporal gyrus; MFG, middle frontal gyrus; PC, posterior cingutate; CAL, cerebellum anterior lobe.
Correlational analysis
In patients with OSA, the MoCA scale scores showed a significant positive correlation with sleep efficiency (r = 0.272, p = 0.008) and a negative correlation with neck circumference (r = −0.243, p = 0.018; Figure 4). No signifificant correlation was found between clinical variables, MoCA scale scores and the abnormal Reho and FC value in OSA-MCI or OSA-nMCI groups.
Figure 4. The MoCA scale scores was negatively correlated with neck circumference (A) and positively correlated with sleep efficiency (B).
Logistic regression
The screened features (Reho of bilateral lingual gyrus, left STG, FC between left STG and B PC/B Calcarine/L CAL) were optimized with hyperparameters and left cross validation to obtain the performance of LG machine learning model, as shown in Figure 4. The AUC was 0.85, the accuracy was 0.87, the sensitivity was 0.70, the specificity was 0.89, the Kappa coefficient was 0.60, and the p-value was 0.0014 after 5,000 permutation tests (Figure 5).
Figure 5. ROC curve of classifier based on altered Reho and FC values, AUC is 0.85. LR, logistic regression; ROC, receiver operating characteristic.
Discussion
In this study, voxel-based Reho and secondary FC were used to explore the abnormal brain activity in OSA-MCI and OSA-nMCI groups. Consistent with our hypothesis, compared with OSA-nMCI group, OSA-MCI group showed significantly lower Reho value in the bilateral lingual gyri and left STG. Decreased FC value was observed between left STG and MTG, left MFG, bilateral PC, calcarine cortex, left ACL, and between bilateral lingual gyrus and bilateral cuneus,which DMN, visual network (VN), and cerebellar networks were mainly involved. These findings indicating that the functional changes of such brain networks may be related to mild cognitive impairment. Secondly, based on the Reho and FC characteristics of significantly different brain regions, we can effectively distinguish whether OSA patients are accompanied by MCI, and provide more evidence for early clinical intervention.
Our study found that in OSA-MCI group, the Reho value of left STG was decreased, the FC value between left STG and left MTG, MFG, and bilateral PC was decreased, and these abnormal brain functional regions were mainly involved in DMN. The DMN is a group of brain regions, including the anterior and medial frontal cortex, bilateral temporal lobes, precuneus, and lateral parietal cortex, that are associated with social behavior, emotional control, memory, learning, and task execution (Raichle, 2015). Ya-ting Chang et al. showed that reduced functional connectivity in the DMN was associated with nocturnal hypoxemia rather than sleep rhythms (Chang et al., 2020). Our previous studies have shown that decreased FC between the DMN and the central executive network and decreased functional connectivity between brain networks are associated with delayed memory (Li et al., 2016). DMN was further divided into aDMN and pDMN (Yang et al., 2017; Wang et al., 2018), aDMN is related to social cognitive functions (emotional self-referential processing and inferential mental states of others), and pDMN is related to cognitive processes (temporal episodic memory and thinking). Zhang et al. (2013) found that FC in the anterior DMN of OSA patients decreased, while FC in the posterior DMN increased compensatively. Similar to previous studies, we found changes in FC values in multiple brain functional areas in the MCI group, which may be related to long-term intermittent hypoxia in OSA patients. In addition, compared with the nMCI group, the FC changes in multiple brain regions (aDMN and pDMN) were decreased, indicating that MCI patients may have more obvious brain function changes, and the possible disruption of compensatory mechanism in the DMN may be the potential cause of MCI.
In our study, we found that the Reho value of bilateral lingual gyrus and FC of left STG and calcarine were decreased in MCI group. The lingual gyrus is mainly involved in the processing of visual memory (Menon and Uddin, 2010) and is an important region of visual perception (Yang et al., 2015). The calcarine is a component of the primary visual cortex and is associated with spatial working memory response time (Valenzuela et al., 2015). Both lingual gyrus and calcarine cortex are involved in the visual network, and the functional changes of visual network are related to MCI (Amaefule et al., 2021). In patients with Alzheimer’s disease, task-state and resting-state MRI studies found that decreased FC of visual network function was significantly associated with cognitive impairment, and changes in visual network integrity were associated with the progression of MCI (Huang et al., 2021). Yang et al. (2022) showed that compared with healthy people, OSA patients took longer to process visual information, which reflected impaired visual perception from the perspective of electrophysiology, and suggested that impaired visual perception may be the potential mechanism of cognitive impairment. In our study, compared with the OSA-nMCI group, the function of brain areas related to visual network was impaired in the MCI group, which may be related to hypoxemia and oxidative stress (Archer and Pepin, 2013). However, reduced Reho in the bilateral lingual gyrus and FC in the left STG and talus were not associated with cognitive scales. Therefore, the potential causal relationship between changes in the visual network and cognitive dysfunction needs further investigation.
The cerebellum is a key module not only in the motor control system, but also in cognitive and emotional processing (Stoodley et al., 2012). Ping Xiao’s study showed that cerebellar gray matter thinning and decreased blood flow in OSA patients and found that structural changes in cerebellum and changes in blood circulation were associated with cognitive impairment (Xiao et al., 2022). Several other studies have shown that the cerebellum is vulnerable to hypoxia and ischemia, and that sleep deprivation is also an important factor affecting cerebellar function (Gazes et al., 2012; Kim et al., 2016). In a study of healthy people, small-world properties of the cerebellum and brain-cerebellum functional coupling were shown (Chen et al., 2022). A resting-state function-based study found that intermittent hypoxia resulted in impaired cerebellar network integration and cerebellar functional connectivity and was associated with cognitive impairment (Park et al., 2022). In this study, our results showed that FC decreased in the anterior cerebellar lobe and left STG, and this decreased functional connectivity in the cerebellum may be related to hypoxia or ischemia. Therefore, we hypothesized that hypoxia or sleep disturbance may lead to the disruption of functional connectivity between cerebellum and brain, which may be one of the causes of cognitive impairment.
In addition, we found that MoCA scale scores were associated with neck circumference and sleep efficiency. Increased neck circumference is associated with obesity, which is one of the risk factors for OSA. Obesity plays an important role in the pathogenesis of apnea syndrome by altering the caliber of the upper airway during breathing in subjects and thus in the pathogenesis of apnea syndrome (Kim et al., 2014). An arterial MRI study showed that the posterior palatal airway was narrower in obese subjects during respiratory arousal and was associated with OSA severity (Feng et al., 2018). Our results showed that the range of neck circumference was significantly different between the OSA-MCI and OSA-nMCI groups. Therefore, we hypothesized that the increase of neck circumference may imply changes in airway anatomy, which may be related to the severity of OSA, although no significance difference was found in BMI or AHI indices between the two groups. Characteristics of sleep disorders include decreased sleep duration and quality, reduced sleep efficiency, fragmented sleeping and sleepiness all day (Casagrande et al., 2022). Previous studies have also shown that sleep disorders are associated with declined cognitive level (Bubu et al., 2020), particularly impairements showed in attention, memory and executive. Our study showed a positive correlation between sleep efficiency and cognitive performance, suggesting that the decline in sleep efficiency may be associated with cognitive impairment.
At last, based on Reho and FC features, LR was used to distinguish OSA-MCIs and OSA-nMCIs. Khatri and Kwon (2022) showed that an efficient classification model was established based on various rs-fMRI and structural MRI, combined with machine learning, indicating that rs-fMRI can be used as a neuroimaging marker to solve prediction and classification. In this study, LR models with excellent classification accuracy (AUC = 0.85, accuracy rate 0.87, sensitivity 0.70, specificity 0.89) were constructed, which suggests that LR may be an effective tool to identify OSA-MCI from OSA patients.
Limitation
First, our subjects were male with severe OSA, and the results may not be generalizable to patients with mild OSA or women with OSA. Second, we only analyzed abnormal functional change and ignore potential morphological and microstructural changes, which may help us to deepen our understanding of OSA with cognitive impairment. Therefore in the future, we will combine other fMRI features (e.g., cortical thickness, white matter fiber) to improve the classification efficiency of the model. Finally, the sample size was small. Larger samples will be needed in future studies.
Conclusion
Our results show that OSA patients with MCI have spontaneous brain activity changes and decreased functional connectivity in multiple brain regions, mainly related to DMN, VN, and the cerebellar network. These provide additional information about the underlying neural mechanisms of OSA-related cognitive impairment. At the same time, our results demonstrate that an effective machine learning method can efficiently distinguish whether OSA patients are accompanied by MCI and provide potential imaging markers for clinical treatment.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Ethics Committee of the First Affiliated Hospital of Nanchang University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
YS and XL wrote, reviewed, and revised the manuscript. DP guided and designed the MRI experiment. HL analyzed the resting-state fMRI data. YS and HL analyzed and discussed the ideas of the paper. PY analyzed machine learning. WD, KL, YZ, and WX collected resting-state fMRI data and applied for the ethics approval. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (grant nos. 81860307 and 81560285); the Natural Science Foundation Project of Jiangxi, China (grant nos. 20202BABL216036, 20181ACB20023, and 20171BAB205070).
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.
Footnotes
References
Amaefule, C. O., Dyrba, M., Wolfsgruber, S., Polcher, A., Schneider, A., Fliessbach, K., et al. (2021). Association between composite scores of domain-specific cognitive functions and regional patterns of atrophy and functional connectivity in the Alzheimer's disease spectrum. Neuroimage Clin. 29:102533. doi: 10.1016/j.nicl.2020.102533
Archer, E. L., and Pepin, S. (2013). Obstructive sleep apnea and nonarteritic anterior ischemic optic neuropathy: evidence for an association. J. Clin. Sleep Med. 09, 613–618. doi: 10.5664/jcsm.2766
Berry, R. B., Brooks, R., Gamaldo, C., Harding, S. M., Lloyd, R. M., Quan, S. F., et al. (2017). AASM scoring manual updates for 2017 (version 2.4). J. Clin. Sleep Med. 13, 665–666. doi: 10.5664/jcsm.6576
Bubu, O. M., Andrade, A. G., Umasabor-Bubu, O. Q., Hogan, M. M., Turner, A. D., de Leon, M. J., et al. (2020). Obstructive sleep apnea, cognition and Alzheimer's disease: a systematic review integrating three decades of multidisciplinary research. Sleep Med. Rev. 50:101250. doi: 10.1016/j.smrv.2019.101250
Casagrande, M., Forte, G., Favieri, F., and Corbo, I. (2022). Sleep quality and aging: a systematic review on healthy older people, mild cognitive impairment and Alzheimer's disease. Int. J. Environ. Res. Public Health 19:8457. doi: 10.3390/ijerph19148457
Chang, Y. T., Chen, Y. C., Chen, Y. L., Hsu, S. W., Yang, F. Y., Lee, C. C., et al. (2020). Functional connectivity in default mode network correlates with severity of hypoxemia in obstructive sleep apnea. Brain Behav. 10:e01889. doi: 10.1002/brb3.1889
Chen, L., Fan, X., Li, H., Ye, C., Yu, H., Gong, H., et al. (2018). Topological reorganization of the default mode network in severe male obstructive sleep apnea. Front. Neurol. 9:363. doi: 10.3389/fneur.2018.00363
Chen, L. T., Fan, X. L., Li, H. J., Ye, C. L., Yu, H. H., Xin, H. Z., et al. (2018). Aberrant brain functional connectome in patients with obstructive sleep apnea. Neuropsychiatr. Dis. Treat. 14, 1059–1070. doi: 10.2147/NDT.S161085
Chen, L. M., Li, X. J., Xu, K., Guo, C. L., Zhang, G. L., Han, M., et al. (2021). Functional MRI-based study on neuromechanism of trans-auricular vagus nerve stimulation for treatment-resistant depression. Zhen Ci Yan Jiu 46, 869–874. doi: 10.13702/j.1000-0607.20210241
Chen, K., Xu, Y., Chu, A., Ding, D., Liang, X., Nasreddine, Z. S., et al. (2016). Validation of the Chinese version of Montreal cognitive assessment basic for screening mild cognitive impairment. J. Am. Geriatr. Soc. 64, e285–e290. doi: 10.1111/jgs.14530
Chen, Z., Zhang, R., Huo, H., Liu, P., Zhang, C., and Feng, T. (2022). Functional connectome of human cerebellum. NeuroImage 251:119015. doi: 10.1016/j.neuroimage.2022.119015
Daulatzai, M. A. (2015). Evidence of neurodegeneration in obstructive sleep apnea: relationship between obstructive sleep apnea and cognitive dysfunction in the elderly. J. Neurosci. Res. 93, 1778–1794. doi: 10.1002/jnr.23634
Feng, Y., Keenan, B. T., Wang, S., Leinwand, S., Wiemken, A., Pack, A. I., et al. (2018). Dynamic upper airway imaging during wakefulness in obese subjects with and without sleep apnea. Am. J. Respir. Crit. Care Med. 198, 1435–1443. doi: 10.1164/rccm.201711-2171OC
Gazes, Y., Rakitin, B. C., Steffener, J., Habeck, C., Lisanby, S. H., Butterfield, B., et al. (2012). Dual-tasking alleviated sleep deprivation disruption in visuomotor tracking: an fMRI study. Brain Cogn. 78, 248–256. doi: 10.1016/j.bandc.2012.01.004
Heinzer, R., Marti-Soler, H., and Haba-Rubio, J. (2016). Prevalence of sleep apnoea syndrome in the middle to old age general population. Lancet Respir. Med. 4, e5–e6. doi: 10.1016/S2213-2600(16)00006-0
Hojjati, S. H., Ebrahimzadeh, A., and Babajani-Feremi, A. (2019). Identification of the early stage of Alzheimer's disease using structural MRI and resting-state fMRI. Front. Neurol. 10:904. doi: 10.3389/fneur.2019.00904
Hojjati, S. H., Ebrahimzadeh, A., Khazaee, A., and Babajani-Feremi, A. (2017). Predicting conversion from MCI to AD using resting-state fMRI, graph theoretical approach and SVM. J. Neurosci. Methods 282, 69–80. doi: 10.1016/j.jneumeth.2017.03.006
Hu, J. B., Zhou, W. H., Hu, S. H., Huang, M. L., Wei, N., Qi, H. L., et al. (2013). Cross-cultural difference and validation of the Chinese version of Montreal cognitive assessment in older adults residing in eastern China: preliminary findings. Arch. Gerontol. Geriatr. 56, 38–43. doi: 10.1016/j.archger.2012.05.008
Huang, J., Beach, P., Bozoki, A., and Zhu, D. C. (2021). Alzheimer's disease progressively reduces visual functional network connectivity. J. Alzheimers Dis. Rep. 5, 549–562. doi: 10.3233/ADR-210017
Ibrahim, B., Suppiah, S., Ibrahim, N., Mohamad, M., Hassan, H. A., Nasser, N. S., et al. (2021). Diagnostic power of resting-state fMRI for detection of network connectivity in Alzheimer's disease and mild cognitive impairment: a systematic review. Hum. Brain Mapp. 42, 2941–2968. doi: 10.1002/hbm.25369
Joo, E. Y., Jeon, S., Kim, S. T., Lee, J. M., and Hong, S. B. (2013). Localized cortical thinning in patients with obstructive sleep apnea syndrome. Sleep 36, 1153–1162. doi: 10.5665/sleep.2876
Ju, R., Hu, C., Zhou, P., and Li, Q. (2019). Early diagnosis of Alzheimer's disease based on resting-state brain networks and deep learning. IEEE/ACM Trans. Comput. Biol. Bioinform. 16, 244–257. doi: 10.1109/TCBB.2017.2776910
Kapur, V. K., Auckley, D. H., Chowdhuri, S., Kuhlmann, D. C., Mehra, R., Ramar, K., et al. (2017). Clinical practice guideline for diagnostic testing for adult obstructive sleep apnea: an American Academy of sleep medicine clinical practice guideline. J. Clin. Sleep Med. 13, 479–504. doi: 10.5664/jcsm.6506
Khatri, U., and Kwon, G. R. (2022). Alzheimer's disease diagnosis and biomarker analysis using resting-state functional MRI functional brain network with multi-measures features and hippocampal subfield and amygdala volume of structural MRI. Front. Aging Neurosci. 14:818871. doi: 10.3389/fnagi.2022.818871
Khazaie, H., Veronese, M., Noori, K., Emamian, F., Zarei, M., Ashkan, K., et al. (2017). Functional reorganization in obstructive sleep apnoea and insomnia: a systematic review of the resting-state fMRI. Neurosci. Biobehav. Rev. 77, 219–231. doi: 10.1016/j.neubiorev.2017.03.013
Kim, H., Joo, E., Suh, S., Kim, J. H., Kim, S. T., and Hong, S. B. (2016). Effects of long-term treatment on brain volume in patients with obstructive sleep apnea syndrome. Hum. Brain Mapp. 37, 395–409. doi: 10.1002/hbm.23038
Kim, A. M., Keenan, B. T., Jackson, N., Chan, E. L., Staley, B., Poptani, H., et al. (2014). Tongue fat and its relationship to obstructive sleep apnea. Sleep 37, 1639–1648. doi: 10.5665/sleep.4072
Kong, L., Li, H., Shu, Y., Liu, X., Li, P., Li, K., et al. (2021). Aberrant resting-state functional brain connectivity of insular subregions in obstructive sleep apnea. Front. Neurosci. 15:765775. doi: 10.3389/fnins.2021.765775
Kumar, R., Chavez, A. S., Macey, P. M., Woo, M. A., Yan-Go, F. L., and Harper, R. M. (2012). Altered global and regional brain mean diffusivity in patients with obstructive sleep apnea. J. Neurosci. Res. 90, 2043–2052. doi: 10.1002/jnr.23083
Lee, M. H., Yun, C. H., Min, A., Hwang, Y. H., Lee, S. K., Kim, D. Y., et al. (2019). Altered structural brain network resulting from white matter injury in obstructive sleep apnea. Sleep 42:zsz120. doi: 10.1093/sleep/zsz120
Li, H., Li, L., Kong, L., Li, P., Zeng, Y., Li, K., et al. (2021). Frequency specific regional homogeneity alterations and cognitive function in obstructive sleep apnea before and after short-term continuous positive airway pressure treatment. Nat. Sci. Sleep 13, 2221–2238. doi: 10.2147/NSS.S344842
Li, H. J., Nie, X., Gong, H. H., Zhang, W., Nie, S., and Peng, D. C. (2016). Abnormal resting-state functional connectivity within the default mode network subregions in male patients with obstructive sleep apnea. Neuropsychiatr. Dis. Treat. 12, 203–212. doi: 10.2147/NDT.S97449
Liu, W., Yue, Q., Gong, Q., Zhou, D., and Wu, X. (2021). Regional and remote connectivity patterns in focal extratemporal lobe epilepsy. Ann. Transl. Med. 9:1128. doi: 10.21037/atm-21-1374
Long, Z., Jing, B., Yan, H., Dong, J., Liu, H., Mo, X., et al. (2016). A support vector machine-based method to identify mild cognitive impairment with multi-level characteristics of magnetic resonance imaging. Neuroscience 331, 169–176. doi: 10.1016/j.neuroscience.2016.06.025
Luo, B., Lu, Y., Qiu, C., Dong, W., Xue, C., Liu, D., et al. (2021). Altered regional homogeneity and functional connectivity during microlesion period after deep brain stimulation in Parkinson's disease. Parkinsons Dis. 2021:2711365. doi: 10.1155/2021/2711365
Menon, V., and Uddin, L. Q. (2010). Saliency, switching, attention and control: a network model of insula function. Brain Struct. Funct. 214, 655–667. doi: 10.1007/s00429-010-0262-0
Olaithe, M., Bucks, R. S., Hillman, D. R., and Eastwood, P. R. (2018). Cognitive deficits in obstructive sleep apnea: insights from a meta-review and comparison with deficits observed in COPD, insomnia, and sleep deprivation. Sleep Med. Rev. 38, 39–49. doi: 10.1016/j.smrv.2017.03.005
Park, H. R., Cha, J., Joo, E. Y., and Kim, H. (2022). Altered cerebrocerebellar functional connectivity in patients with obstructive sleep apnea and its association with cognitive function. Sleep 45:zsab209. doi: 10.1093/sleep/zsab209
Park, B., Palomares, J. A., Woo, M. A., Kang, D. W., Macey, P. M., Yan-Go, F. L., et al. (2016). Aberrant insular functional network integrity in patients with obstructive sleep apnea. Sleep 39, 989–1000. doi: 10.5665/sleep.5738
Peppard, P. E., Young, T., Barnet, J. H., Palta, M., Hagen, E. W., and Hla, K. M. (2013). Increased prevalence of sleep-disordered breathing in adults. Am. J. Epidemiol. 177, 1006–1014. doi: 10.1093/aje/kws342
Prilipko, O., Huynh, N., Schwartz, S., Tantrakul, V., Kim, J. H., Peralta, A. R., et al. (2011). Task positive and default mode networks during a parametric working memory task in obstructive sleep apnea patients and healthy controls. Sleep 34, 293–301. doi: 10.1093/sleep/34.3.293
Raichle, M. E. (2015). The brain's default mode network. Annu. Rev. Neurosci. 38, 433–447. doi: 10.1146/annurevneuro-071013-014030
Santarnecchi, E., Sicilia, I., Richiardi, J., Vatti, G., Polizzotto, N. R., Marino, D., et al. (2013). Altered cortical and subcortical local coherence in obstructive sleep apnea: a functional magnetic resonance imaging study. J. Sleep Res. 22, 337–347. doi: 10.1111/jsr.12006
Stoodley, C. J., Valera, E. M., and Schmahmann, J. D. (2012). Functional topography of the cerebellum for motor and cognitive tasks: an fMRI study. NeuroImage 59, 1560–1570. doi: 10.1016/j.neuroimage.2011.08.065
Sun, Q., Fan, W., Ye, J., and Han, P. (2018). Abnormal regional homogeneity and functional connectivity of baseline brain activity in hepatitis B virus-related cirrhosis with and without minimal hepatic encephalopathy. Front. Hum. Neurosci. 12:245. doi: 10.3389/fnhum.2018.00245
Valenzuela, M. J., Turner, A. J., Kochan, N. A., Wen, W., Suo, C., Hallock, H., et al. (2015). Posterior compensatory network in cognitively intact elders with hippocampal atrophy. Hippocampus 25, 581–593. doi: 10.1002/hipo.22395
Van Dijk, K. R., Sabuncu, M. R., and Buckner, R. L. (2012). The influence of head motion on intrinsic functional connectivity MRI. NeuroImage 59, 431–438. doi: 10.1016/j.neuroimage.2011.07.044
Wang, L., Feng, Q., Wang, M., Zhu, T., Yu, E., Niu, J., et al. (2021). An effective brain imaging biomarker for AD and aMCI: ALFF in Slow-5 frequency band. Curr. Alzheimer Res. 18, 45–55. doi: 10.2174/1567205018666210324130502
Wang, C., Pan, Y., Liu, Y., Xu, K., Hao, L., Huang, F., et al. (2018). Aberrant default mode network in amnestic mild cognitive impairment: a meta-analysis of independent component analysis studies. Neurol. Sci. 39, 919–931. doi: 10.1007/s10072-018-3306-5
Xia, W., Chen, Y. C., Luo, Y., Zhang, D. F., Chen, H., Ma, J., et al. (2018). Decreased spontaneous brain activity and functional connectivity in type 1 diabetic patients without microvascular complications. Cell. Physiol. Biochem. 51, 2694–2703. doi: 10.1159/000495960
Xia, Y., Fu, Y., Xu, H., Guan, J., Yi, H., and Yin, S. (2016). Changes in cerebral metabolites in obstructive sleep apnea: a systemic review and meta-analysis. Sci. Rep. 6:28712. doi: 10.1038/srep28712
Xiao, P., Hua, K., Chen, F., Yin, Y., Wang, J., Fu, X., et al. (2022). Abnormal cerebral blood flow and volumetric brain Morphometry in patients with obstructive sleep apnea. Front. Neurosci. 16:934166. doi: 10.3389/fnins.2022.934166
Yang, Y. L., Deng, H. X., Xing, G. Y., Xia, X. L., and Li, H. F. (2015). Brain functional network connectivity based on a visual task: visual information processing-related brain regions are significantly activated in the task state. Neural Regen. Res. 10, 298–307. doi: 10.4103/1673-5374.152386
Yang, C., Wang, C., Chen, X., Xiao, B., Fu, N., Ren, B., et al. (2022). Event-related potential assessment of visual perception abnormality in patients with obstructive sleep apnea: a preliminary study. Front. Hum. Neurosci. 16:895826. doi: 10.3389/fnhum.2022.895826
Yang, H., Wang, C., Zhang, Y., Xia, L., Feng, Z., Li, D., et al. (2017). Disrupted causal connectivity anchored in the posterior cingulate cortex in amnestic mild cognitive impairment. Front. Neurol. 8:10. doi: 10.3389/fneur.2017.00010
Yayan, J., Rasche, K., and Vlachou, A. (2017). Obstructive sleep apnea and chronic kidney disease. Adv. Exp. Med. Biol. 1022, 11–18. doi: 10.1007/5584_2017_35
Yu, H., Chen, L., Li, H., Xin, H., Zhang, J., Wei, Z., et al. (2019). Abnormal resting-state functional connectivity of amygdala subregions in patients with obstructive sleep apnea. Neuropsychiatr. Dis. Treat. 15, 977–987. doi: 10.2147/NDT.S191441
Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. NeuroImage 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030
Keywords: obstructive sleep apnea, mild cognitive impairment, rs-fMRI, regional homogeneity, functional connectivity
Citation: Shu Y, Liu X, Yu P, Li H, Duan W, Wei Z, Li K, Xie W, Zeng Y and Peng D (2022) Inherent regional brain activity changes in male obstructive sleep apnea with mild cognitive impairment: A resting-state magnetic resonance study. Front. Aging Neurosci. 14:1022628. doi: 10.3389/fnagi.2022.1022628
Edited by:
Joel Ramirez, University of Toronto, CanadaReviewed by:
Haiqing Huang, University of Pittsburgh, United StatesSy Duong-Quy, Lam Dong Medical College, Vietnam
Copyright © 2022 Shu, Liu, Yu, Li, Duan, Wei, Li, Xie, Zeng and Peng. 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.
*Correspondence: Dechang Peng, cGVuZ2RjZG9jdG9yQDE2My5jb20=
†These authors have contributed equally to this work