Effect of obstructive sleep apnea on cerebrovascular compliance and cerebral small vessel disease

Reduced cerebrovascular compliance is the major mechanism of cerebral small vessel disease (SVD). Obstructive sleep apnea (OSA) also promotes SVD development, but the underlying mechanism was not elucidated. We investigated the association among OSA, cerebrovascular compliance, and SVD parameters. This study retrospectively included individuals ≥ 50 years of age, underwent overnight polysomnographic (PSG) for the evaluation of OSA, and performed MRI and transcranial Doppler (TCD) within 12 months of interval without a neurological event between the evaluations. TCD parameters for the cerebrovascular compliance included middle cerebral artery pulsatility index (MCA PI) and mean MCA resistance index ratio (MRIR). SVD parameters included white matter hyperintensity (WMH) volume, number of lacunes, enlarged perivascular space (ePVS) score, and the presence of microbleeds or lacunes. Ninety-seven individuals (60.8% male, mean age 70.0±10.5 years) were included. MRIR was associated with higher respiratory distress index (B = 0.003; 95% confidence interval [CI] 0.001–0.005; P = 0.021), while MCA PI was not associated with any of the PSG markers for OSA severity. Apnea-hypopnea index was associated with the log-transformed total WMH volume (B = 0.008; 95% confidence interval [CI] 0.001–0.016; P = 0.020), subcortical WMH volume (B = 0.015; 95% CI 0.007–0.022; P<0.001), total ePVS score (B = 0.024; 95% CI 0.003–0.045; P = 0.026), and centrum semiovale ePVS score (B = 0.026; 95% CI 0.004–0.048; P = 0.019), and oxygen-desaturation index with periventricular WMH volume, independently from age, MCA PI, and MRIR. This study concluded that OSA is associated with reduced cerebrovascular compliance and also with SVD independently from cerebrovascular compliance. Underlying pathomechanistic link might be region specific.


Introduction
Cerebral small vessel disease (SVD) is a very prevalent age-related phenomenon of the brain [1,2]. There is no approved treatment for SVD, although a recent randomized study has shown that remote ischemic conditioning reduces the recurrence of ischemic stroke and is an effective therapy for patients with SVD [3]. Reduced compliance in the cerebral arterioles resulting from chronic dysregulated vascular remodeling has been regarded the fundamental pathomechanism of SVD progression, which provokes the blood-brain-barrier (BBB) damage, intermittent hypoperfusion, and impaired clearance of the cerebral waste and ultimately leads to the chronic inflammation and subclinical ischemia in the brain parenchyma [2,[4][5][6][7].
Obstructive sleep apnea (OSA) is a very prevalent sleep related breathing disorder in elderly population. The reported incidence of OSA is from 5.6% to 60% in people over 65 while SVD presents an increased one from 6% at age 60 to 28% at age 80 [8]. Aging is associated with an increased risk and severity of OSA, and OSA has also been shown to be related to olfactory disorders, the origin of which could be attributed to disorders of the central nervous system [9].
OSA and SVD share common risk factors such as increasing age, hypertension, diabetes, and obesity, and OSA is associated with the progression of SVD [10][11][12]. OSA is implicated in the endothelial dysfunction and decreased vascular compliance [7,9,13,14]. A recent study suggested that endothelial dysfunction in OSA may be related to the consequent generation of reactive oxygen species and pro-inflammatory molecules, resulting in microvascular damage in OSA patients [15]. However, the effect of OSA on SVD progression is independent from those cerebrovascular risk factors and its underlying mechanisms might be distinct from the dysregulated vascular remodeling [16,17]. The candidate pathophysiologic link between OSA and SVD includes abrupt increment in the intrathoracic pressure that interferes with the adequate venous return and cardiac output, frequent arousal, intermittent brain hypoxia, and provocation of arrhythmia [13,[18][19][20]. Frequent arousal and intermittent hypoxia can recurrently activate the sympathetic nervous system, provoke oxidative stress, and induce parenchymal inflammation [13,18,[21][22][23]. Also, altered intrathoracic pressure and frequent arousal might impair the activation of cerebral waste clearance system during sleep [13,19,20,[24][25][26]. Although highly interrelated, discrimination of the major mechanism underlying the progression of SVD in patients with OSA might be important to predict how the adequate treatment of OSA could also modify the progression of SVD [27,28].
Transcranial Doppler sonography (TCD) is a widely used noninvasive method to measure the blood flow in the brain [29]. Some parameters of TCD are useful indicators for the cerebrovascular compliance [30]. Pulsatility index (PI) evaluates the stiffness of the vessel and also the resistance of the distal arterial bed [31,32]. Recently, it is reported that a reduction of resistance index (RI) along the middle cerebral artery (MCA RI ratio, MRIR) might be an indirect marker for the compliance of the cerebral small vessels [30].
In this study, we aimed to determine the pathomechanism linking OSA, altered cerebrovascular compliance, and SVD development. In this regard, we investigated whether the polysomnographic (PSG) parameters for the OSA severity was associated with the TCD parameters for the cerebrovascular compliance. Then, we evaluated the associations of the OSA severity factors and the vascular compliance markers with the parameters of cerebral SVD, such as white matter hyperintensity (WMH) volume, enlarged perivascular space (ePVS) score, and the presence of cerebral microbleeds (CMB) or lacunes.

Study subjects and clinical data
The design of this study was approved by the institutional review board of Seoul National University Boramae Hospital. Informed consent was waived by the IRB, as this study did not intervene with the evaluation or treatment process of the patients, and the patients were anonymized during the study process. This retrospective cohort study initially included all consecutive individuals who were � 50 years of age, visited Seoul National University Boramae Hospital between March 2008 and August 2019, underwent overnight PSG evaluation to evaluate for the presumptive diagnosis of OSA, and performed MRI and TCD within 12 months of interval. Among the initially included 175 individuals, the final study population was defined according to the following criteria: (1) no significant (� 50%) stenosis in the internal carotid arteries (ICAs) or the MCAs at the initial magnetic resonance angiography (MRA) or TCD; (2) no history of stroke other than an old lacunar stroke (> 90 days); (3) able to carry out daily activities independently without a neurological deficit; and (4) adequate temporal window for a TCD evaluation. According to the criteria, 15 patients with significant stenosis in the ICAs or the MCAs, 45 with a stroke history other than an old lacunar stroke, 11 with fixed neurological deficit or incapable of independent daily living, and 7 with inadequate temporal window were sequentially excluded and the remaining 97 individuals were included for the study analysis. Clinical profiles including demographics, history of lacunar infarction, and presence of hypertension, diabetes mellitus, hyperlipidemia, congestive heart failure, and atrial fibrillation were obtained from the patients' medical record. The design of this study was approved by the institutional review board of Seoul National University Boramae Hospital.

Transcranial Doppler evaluation
Intracranial arteries were sonographically evaluated using a 2-MHz pulsed-wave and rangegated TCD probe (Digital PMD 100 or ST3 Digital PMD 150; Spencer Technologies; Redmond, WA, USA) with a transmission power level of 100 mW/cm 2 , pulse repetition frequency of 8000 Hz, filter frequency of 200 Hz, gain of 4 decibels, and range of 200 decibels. The TCD evaluation protocol was standardized for every patient and was conducted by two skilled sonographers with 17 and 10 years of experience, respectively. The flow in the MCA, the anterior cerebral artery, and the internal carotid artery (ICA) was evaluated to exclude stenosis of the anterior cerebral arteries. Peak systolic, minimal diastolic, and mean flow velocities (PSV, MDV, and MFV, cm/s), pulsatility index (PI), and resistance index (RI) were obtained at five points along the M1 portion of the MCA with 5mm intervals insonation depths from the temporal windows [29,33]. MRIR was designated as the mean ratio of RI along M1 (RI at the shallowest [the most distal] point/ RI at the deepest [the most proximal] point) in each hemisphere [30]. The mean of PI in both MCAs was defined as the average of PIs btained at five insonation points in each hemisphere [25,30].
immediately preceding basal value, and arousal index as the mean number of arousals per hour [35,36]. A significant OSA (moderate to severe degree) was defined as an AHI value of �15/h [34].
The volumetric analysis of WMH was performed by a neurologist (WJL) blinded to other data, according to the protocols previously described. In brief, FLAIR sequences were registered into an offline workstation. A freeware NeuRoi (Nottingham university, Nottingham, UK) was used to identify the WMH lesion, designate its boundaries, calculate the volume [25,30,37]. The periventricular and subcortical WMH volume was separately measured and normalized by the total intracranial volume [30,38]. The reproducibility of the volumetric analysis was established as good in the previous report [30].
ePVS was designated as small and sharply demarcated ovoid or linear shaped lesions which appeared as hyperintensity in T2 weighted images and hypointensity in T1 weighted images. ePVS was rated in basal ganglia (BG) and centrum semiovale (CS) using the Wardlaw scale of which score range is 0-4 (0 for no, 1 for 1-10, 2 for 11-20, 3 for 21-40, and 4 for >40 ePVS) and total score range is 0-8 [7,39]. CMB was defined as a 2-10 millimeter punctate hypointensity in GRE/SWI images and lacune was designated as a hyperintense lesion with central hypointensity in FLAIR images [7]. SPSS 25.0 (SPSS Inc., Chicago, IL, USA) was used for all statistical analyses. Pearson coefficients was used to measure correlations between continuous variables, Spearman's Rho to measure correlations between ordinal variables, paired t test to compare mean values, and chisquare test to compare the frequencies. Age, sex, and variables with P values <0.15 in univariate analyses were included in regression analyses using backward elimination. Linear regression analyses were performed to evaluate the factors associated with WMH volume parameters. WMH volumes were log-transformed to obtain a normal distribution. Multi-collinearity between variables was estimated using the variance inflation factor (VIF). A scatterplot of the standardized predicted values and the standardized residuals was used to check for the linearity of the regression model. Ordinal regression analyses were performed to evaluate the factors associated with ePVS scores. For all analyses, P values < 0.05 were considered as statistically significant.

Results
Ninety-seven patients, including 59 males (60.8%) and 38 females (39.2%), with a mean age of 70.0±10.5 years (range, 50-88 years) were included. The median interval between the evaluations were 82 days (interquartile range, IQR 49-141 days, range -186-236 days) for the PSG and the MRI, 62 days (IQR 28-123 days, range -283-361 days) for the PSG and the TCD, and 17 days (IQR 10-23 days, range -356-72 days) for the MRI and TCD. The indications for the MRI and the TCD analysis were for follow-up evaluation of a pre-existing cerebral WMH in 16 (16.5%) patients or of an old lacunar infarction in 9 (9.3%) patients, and investigations for nonspecific neurologic symptoms such as headache, dizziness, or subjective memory impairment in 72 (74.2%) patients.
In the univariate analyses for the MCA PI, higher age, shorter total sleep time, lower sleep efficiency, and higher MCA MFV were associated with higher MCA PI, while none of the OSA parameters or other clinical parameters showed a significant association with MCA PI (S1 and S2 Tables). Following multivariate linear regression analysis revealed that higher age was the only parameter significantly associated with MCA PI (B coefficient = 0.007; 95% confidence interval [CI] 0.004-0.009; P<0.001, Table 2).
In univariate analyses for ePVS, total ePVS score was associated with higher age, AHI, RDI, MRIR, and MCA PI, and the presence of diabetes mellitus (S1 and S2 Tables). Following multivariate ordinal regression analysis revealed that total ePVS score was associated with age

PLOS ONE
Association of OSA with cerebrovascular compliance and cerebral SVD 2.551; P = 0.001), but only marginally with RDI (B = 0.022; 95% CI -0.001-0.045; P = 0.060, Tables 4 and S1 and S2 for the univariate analyses). None of the clinical, PSG, or TCD parameters was associated with the presence of CMB or lacunes. Multivariate analysis was not performed due to the small numbers of the patients with CMB or lacunes (S3 Table).

Discussion
The purpose of the current study was to investigate whether the OSA severity was associated with the cerebrovascular compliance and to evaluate the associations of the OSA severity factors with the SVD parameters. This study found that RDI was significantly associated with the TCD marker for cerebral microvascular compliance, MRIR, but not with another marker for cerebrovascular remodeling, MCA PI. AHI and those cerebral vascular remodeling parameters were independently associated with the cerebral WMH volume and ePVS scores. Additionally, AHI was significantly associated with the subcortical SVD parameters such as subcortical WMH volume and CS ePVS score, whereas ODI was associated with WMH volume in deep structures. These findings are concordant to that of a recent meta-analysis, which reported that moderate to severe OSA is positively associated with WMH and silent brain infarction, but not with CMBs [40]. The major finding of the current study is that OSA is associated with cerebral SVD independently from the reduced cerebrovascular compliance. This finding suggests that OSA influences the pathogenesis of SVD by multiple mechanisms and some of them might be distinct from the chronic cerebrovascular remodeling. Considering that chronic aging related cerebrovascular remodeling process has progressive and irreversible feature [2,41,42], the finding of this study also indicates that the effect of OSA on the SVD might be reversible in some extent. Previous studies reported that continuous positive airway pressure (CPAP) treatment reversed the altered diffusion tensor image parameters for the white matter integrity and the voxelbased morphometry based gray matter volume parameters in patients with OSA, and explained the reversibility by the improvement of the early osmotic changes in the cells induced by mild ischemia and the inflammation that impairs the cell membrane integrity [27,28]. Additionally, surgical treatment for OSA by relocation pharyngoplasty, a variant of uvulopalatopharyngoplasty after, also improved the level of high sensitivity C-reactive protein and reduced cardiovascular risk in subjects with OSA [43,44].
Some plausible pathophysiologic link might explain the association between OSA and SVD independent from the vascular remodeling process. Recent findings that sleep activates the perivascular system function, which is a major pathway for the clearance of cerebral wastes to maintain brain homeostasis of which dysfunction is associated with the pathogenesis of SVD [19,45,46]. Sleep induced amplification of the very low frequency oscillations in the cerebral flow velocity might mediate the activation of this system, and disturbed sleep is associated with the impaired regulation of the major cerebral wastes such as β-amyloid or tau proteins [19,25]. Adequate regulation of the cerebral blood flow pulsation and the CSF recirculation is required for the functioning of the cerebral waste clearance system [37,47]. Therefore, respiratory event related abrupt intrathoracic pressure increment and the disruption of the sleep structures by frequent arousals might significantly interfere with the cerebral waste clearance during sleep [13,26]. Additionally, intermittent brain hypoxia provokes oxidative stress and activates inflammatory cascade, which results in the disruption of the BBB integrity and the glial and neuronal cell damage [18,24]. Respiratory event during sleep also provokes paroxysmal cardiac arrhythmia that is associated with the pathogenesis of SVD [20,48]. As those pathologic effects can be efficiently prevented by CPAP treatment, OSA might be considered as another modifiable risk factor for the cerebral SVD progression, independent from the conventional vascular risk factors (Fig 1).
Another major finding of the current study is that the effect of AHI on the SVD was different between the brain regions. AHI was significantly associated with the subcortical SVD parameters. However, the SVD parameters of periventricular or basal ganglia was not associated with AHI but rather with ODI, and the conventional cerebrovascular compliance factors such as the presence of hypertension, MCA PI, and MRIR. This might be due to the difference in the dominant pathomechanism of SVD between the subcortical and deep brain structures. As periventricular and basal ganglia have less production of cerebral wastes and are supplied the highly compliant perforating arterioles ramifying directly and perpendicularly from the M1 of MCA, those regions might be less dependent to the sleep-related cerebral waste clearance system function and more susceptible to the conventional cerebrovascular remodeling or the hypoxic damage resulting from the respiratory events during sleep, measured by ODI [7,35,37,38,47,49].
Additionally, this study also observed that RDI was associated with higher MRIR, a marker for a reduced compliance of the perforating arterioles ramifying from M1 of MCA [30]. Considering that RDI reflects the burden of sympathetic activation by sleep related respiratory events, and is associated with the development of hypertension and coronary artery disease [50,51], this association indicates that sympathetic activation during sleep has a major role in the chronic reduction of the cerebral microvascular compliance.
Several limitations in this study need to be addressed. First, due to the retrospective design, the clinical profiles of the study population and the indications for the MRI and the TCD evaluations were heterogeneous. For instance, a history of lacunar infarction or subjective memory impairment might have implicated in the severity of cerebral WMH. Second, the intervals between the PSG, MRI, and TCD were unstandardized. Therefore, the WMH burden at the time of MRI evaluation might be different from that at the time when PSG or the TCD was The red lines with arrowheads denotes the pathophysiologic mechanisms dependent on the cerebral vascular remodeling process and the blue lines denotes the mechanisms independent from the cerebral vascular remodeling process. SNS: sympathetic nerve system, BBB: blood brain barrier, and SVD: small vessel disease.
https://doi.org/10.1371/journal.pone.0259469.g001 performed, although the patients were stable with neurological status during the study intervals. Third, due to the relatively small study population, some factors with actual pathophysiologic relevance might have not reached a statistical significance. Fourth, due to the crosssectional design, the findings of the current study did not ensure the causative relationship among the OSA, vascular remodeling, and the progression of WMH. Future prospective studies with larger study populations, standardized follow-up protocols, and evaluation for the effect of CPAP treatment are warranted to clarify the clinical significance of OSA in the pathomechanism of cerebral WMH and the possible therapeutic effect.
Supporting information S1