Brain tissue integrity mapping in adults with obstructive sleep apnea using T1-weighted and T2-weighted images

Abstract

Background:

Obstructive sleep apnea (OSA) is accompanied by both gray and white matter differences in brain areas that regulate autonomic, cognitive, and mood functions, which are deficient in the condition. Such tissue changes have been examined through diffusion tensor and diffusion kurtosis imaging-based procedures. However, poor in-plane spatial resolution of these techniques precludes precise determination of the extent of tissue injury. Tissue texture maps derived from the ratio of T1-weighted and T2-weighted images can provide more adequate in-plane assessment of brain tissue differences.

Objectives:

To examine brain tissue integrity in recently diagnosed, treatment-naïve OSA subjects, relative to age- and sex-comparable control subjects using T1-weighted and T2-weighted images.

Design:

A cross-sectional study.

Methods:

We examined the extent of tissue changes in 106 OSA over 115 control subjects using high-resolution T1- and T2-weighted images collected from a 3.0-Tesla scanner (analysis of covariance; covariates: age, sex, body-mass-index, Pittsburgh sleep quality index, Epworth sleepiness scale, Beck Anxiety Inventory, and Beck Depression Inventory II; false discovery rate corrected; p < 0.01).

Results:

OSA subjects showed significantly lowered tissue integrity in several brain regions, including the frontal, cingulate and insular cortices, cingulum bundle, thalamus, corpus callosum, caudate and putamen, pons, temporal, occipital, and parietal sites, cerebellar peduncles, and medial medullary sites, compared with controls.

Conclusion:

OSA subjects show widespread lowered tissue integrity in autonomic, mood, and cognitive control sites over healthy controls. The pathological processes contributing to the alterations may include repetitive hypoxic and hypercarbic processes and excitotoxic injury, leading to altered brain tissue integrity in OSA.

Keywords

autonomic cognition hippocampus insula mood tissue injury

Introduction

Obstructive sleep apnea (OSA) is a condition characterized by persistent events of partial or complete upper airway obstruction, with continued diaphragmatic efforts to breathe during sleep. In addition to breathing issues, OSA subjects show severe deficits in autonomic functions, expressed principally as blood pressure and sweating disturbances, and cognitive deficits, along with depression and anxiety.¹ Untreated OSA can lead to severely reduced attention, memory, and cognitive functions, and consequently reduced quality of life.² Gray matter injury in OSA subjects is principally located in autonomic, respiratory, mood, and cognitive control sites, as observed in various studies, including observations with procedures such as voxel-based morphometry, diffusion tensor imaging, magnetic resonance spectroscopy, and manual volumetric-based magnetic resonance imaging (MRI).^3–7 Regional brain gray matter changes can be accompanied by subcortical and white matter alterations, as reported in other conditions.^7–9 However, the extent of gray and white matter microstructural differences that may result from inflammation or edema in OSA subjects is unclear.

Numerous MRI methods, including magnetization transfer imaging, diffusion tensor imaging-based radial diffusivity, and diffusion kurtosis imaging-based radial kurtosis procedures are sensitive to examine changes in gray and white matter tissue integrity.^10–13 Precise determination of brain tissue integrity is essential for protection from further injury, or to intervene with restorative procedures after initial damage; there is some evidence that aggressive intervention for OSA with breathing support can at least partially recoup from early injury.^14,15 Assessing the nature of tissue injury, and determining potential therapies for restoring the integrity may hasten that recovery.

T1- to T2-weighted ratio imaging was initially developed for brain myelin mapping at cortical sites,¹⁶ and later extended to accommodate data from different scanners for subject comparisons.¹⁷ Both T1- and T2-weighted images show different signal intensities in white matter regions compared with surrounding gray and deep-gray matter areas due to the presence of myelin, and white matter sites show higher signal intensity on T1-weighted than in T2-weighted images. However, when brain tissue integrity is compromised, signal intensities of both T1-weighted and T2-weighted images vary, but in opposite directions, and signal reductions on T1-weighted images are more prominent than T2-weighted, which leads to enhanced microstructural differences on T1-weighted to T2-weighted ratio maps. The T1-weighted to T2-weighted ratio map has been used to assess tissue integrity in some clinical conditions,⁸ and thus, the approach may be useful to examine tissue integrity in gray and white matter regions in OSA subjects.

Our aim was to examine brain tissue integrity in recently diagnosed, treatment-naïve OSA subjects, relative to age- and sex-comparable control subjects using T1-weighted and T2-weighted images. Such precise assessments are essential for determining potential therapies for restricting further or rescuing initial white matter damage, since early tissue changes are partially restorable with aggressive breathing intervention in OSA subjects.^14,15 We hypothesized that this technique would be sensitive to detect gray and white matter microstructural changes due to repetitive hypoxic processes that lead to deficit in autonomic, respiratory, cognitive, and mood functions in OSA.

Methods

Subjects

A total of 106 newly diagnosed, treatment-naïve OSA and 115 age- and sex-comparable control subjects were included in this study (Table 1). Control subjects were recruited through posting flyers at the University of California Los Angeles (UCLA) campus and West Los Angeles Area. No control subjects were taking any medication as described below that might alter brain tissue. OSA subjects were recruited from the Sleep Disorders Center at the UCLA Medical Center, were newly diagnosed via overnight polysomnography with at-least mild severity (apnea-hypopnea-index ⩾ 5 events/h), and were treatment-naïve for breathing issues. OSA subjects were not taking any mood regulating or cardiovascular medications, such as serotonin reuptake inhibitors, α-agonists, angiotensin-converting enzyme inhibitors, or vasodilators. All OSA and control subjects were compatible with the MRI scanner environment and were without history of stroke, heart failure, diagnosed brain conditions that may introduce brain changes, or metallic implants. Although all OSA subjects were newly diagnosed, of 106 OSA, 43 were taking multivitamins, 19 were on blood pressure medication, 18 were on cholesterol medication, 21 were taking pain medication, 7 were taking anti-diabetes medication, 5 were using proton pump inhibitors, 4 were taking anti-depression medication, and 6 were taking sleep medication. Of 106 OSA subjects, 7 subjects were diagnosed with Type 2 diabetes, 19 with hypertension, and 6 with thyroid issue. All control subjects included here were healthy, without taking any medication, and without any neurological or psychiatric condition.

Table 1.

Demographic, physiologic, neuropsychologic, sleep, and cognitive data of control and OSA subjects.

Variables	Control (mean ± SD)[n = 115]	OSA (mean ± SD)[n = 106]	p-values
Age (years)	47.4 ± 9.2	48.1 ± 9.6	0.578
Gender (Male: Female)	71:44	77:29	0.085
BMI (kg/m²)	25.1 ± 3.5	30.9 ± 5.8	<0.001
Heart rate (beats/min)	70.3 ± 10 (n = 64)	73.8 ± 11.8 (n = 56)	0.096
Systolic BP (mmHg)	117.3 ± 15.1 (n = 64)	126.9 ± 16.1 (n = 56)	0.001
Diastolic BP (mmHg)	74.8 ± 10.2 (n = 64)	79.4 ± 9.4 (n = 56)	0.011
PSQI	3.7 ± 2.4	8.9 ± 4.1	<0.001
ESS	4.9 ± 3.3	9.6 ± 4.7	<0.001
BDI-II	3.3 ± 4.2	8.1 ± 7.3	<0.001
BAI	3.2 ± 4.5	8.5 ± 10	<0.001
AHI (events/h)	–	35.8 ± 22.5	–
SaO₂ baseline (%)	–	94.7 ± 1.9 (n = 105)	–
SaO₂ nadir (%)	–	79.8 ± 8.6 (n = 105)	–
Global MoCA scores	26.9 ± 2.3 (n = 65)	26.3 ± 3.1 (n = 57)	0.235
MoCA: Visuospatial	4.5 ± 0.8 (n = 65)	4.1 ± 0.9 (n = 57)	0.011
MoCA: Naming	2.9 ± 0.3 (n = 65)	2.9 ± 0.2 (n = 57)	0.621
MoCA: Attention	5.4 ± 0.9 (n = 65)	5.4 ± 1.0 (n = 57)	0.749
MoCA: Language	2.6 ± 0.7 (n = 65)	2.5 ± 0.8 (n = 57)	0.248
MoCA: Abstraction	2.0 ± 0.2 (n = 65)	1.9 ± 0.5 (n = 57)	0.141
MoCA: Delayed recall	3.5 ± 1.3 (n = 65)	3.4 ± 0.5 (n = 57)	0.693
MoCA: Orientation	5.9 ± 0.2 (n = 65)	5.9 ± 0.2 (n = 57)	0.870
Trail A	29.4 ± 11.4	30.3 ± 14.4	0.606
Trail B	63.5 ± 24.9	68.4 ± 31.4	0.198

AHI, Apnea-Hypopnea-Index; BAI, Beck Anxiety Inventory; BDI-II, Beck Depression Inventory II; BMI, body mass index; ESS, Epworth sleepiness scale; MoCA, Montreal Cognitive Assessment; OSA, Obstructive sleep apnea; PSQI, Pittsburgh sleep quality index; SaO₂, Oxygen saturation; SD, standard deviation.

Assessment of sleep, mood, and anxiety

To investigate sleep quality and daytime sleepiness, two self-administered questionnaires were employed.¹⁸ In both OSA and control subjects, sleep quality was evaluated based on the Pittsburgh sleep quality index (PSQI), and daytime sleepiness based on the Epworth Sleepiness Scale (ESS). The Beck anxiety inventory (BAI) and Beck depression inventory II (BDI-II) questionnaires, which are self-administered questionnaires on symptoms severity (21 questions; each score ranged from 0 to 3), with each score ranging from 0 to 63, were used for assessment of anxiety¹⁹ and depressive²⁰ symptoms in OSA and control subjects.

Cognition assessment

The Montreal Cognitive Assessment (MoCA) test was performed for cognitive assessment in OSA and control subjects. This test includes fast screening of various cognitive domains, including attention and concentration, language, memory, executive functions, visuo-constructional skills, conceptual thinking, calculations, and orientation. Trail-making A & B tests, which are used for screening of dementia and cognitive functions, including executive function, visual search, mental flexibility, scanning and processing speed, were conducted in OSA and control subjects.

Magnetic resonance imaging

Brain images from OSA and control subjects were acquired from 3.0-Tesla MRI (Siemens, Magnetom Tim-Trio and Prisma, Erlangen, Germany) scanners. Foam pads were used on either side of the head to reduce motion-related artifacts during scanning. Two high-resolution T1-weighted scans were acquired with magnetization prepared rapid acquisition gradient echo pulse sequence in the sagittal plane [repetition time (TR) = 2200 ms; echo time (TE) = 2.4, 2.6 ms; inversion-time = 900 ms; flip-angle (FA) = 9^o; matrix size = 256 × 256, 320 × 320; field of view (FOV) = 230 × 230 mm²; slice thickness = 0.9, 0.72 mm). Proton density (PD) and T2 weighted images [TR = 10,000 ms; TE1, 2 = 12, 124 ms; FA = 130°] were acquired simultaneously using a dual-echo turbo spin-echo pulse sequence in the axial plane, with a 256 × 256 matrix size, 230 × 230 mm² FOV, 4.0, 3.5 mm slice thicknesses, and no inter-slice gap. Two control subjects were scanned twice, once on a Tim-Trio and the other on a Prisma MRI scanner.

Theoretical background ─ gray and white matter changes

Dividing the T1-weighted image by the aligned T2-weighted image mathematically cancels the signal intensity bias related to the sensitivity profile of the radio frequency receiver coils, which is assumed to be the same in both images.¹⁶ Therefore, calculating the ratio of T1-weighted/T2-weighted increases gray and white matter microstructural contrast. However, the contrast obtained through this approach is not completely bias-free due to the transmit field differences that remain.¹⁶ Previously, this technique was implemented for data acquired across the same scanner,¹⁶ and later advanced by introducing scaling factors to combine data across various scanners with same field strength.¹⁷ In case data need to be combined with different acquisition parameters a scaling factor can be introduced to lower bias so that the T1-weighted/T2-weighted intensity scaling is comparable across subjects.¹⁷ Thus, mathematically, the ratio of T1-weighted/T2-weighted images, could be modeled as:

\frac{T 1 - w}{T 2 - w} = \frac{α_{1} * x}{α_{2} * \frac{1}{x}} = \frac{α_{1}}{α_{2}} x^{2} = β x^{2}

(1)

where tissue integrity is represented by ‘x’, and variables α₁, α₂ are scaling factors related to bias field inhomogeneity of T1-weighted and T2-weighted images, respectively.

Data processing and analysis

We used various software tools for image visualization, data processing, and analyses that included the statistical parametric mapping package SPM12 (Wellcome Department of Cognitive Neurology, UK; http://www.fil.ion.ucl.ac.uk/spm), MRIcroN,²¹ and MATLAB-based custom routines (MATLAB, The MathWorks Inc., Natick, MA). For all OSA and control subjects, both T1- and T2-weighted images were visually examined to ensure that no serious anatomical defects were apparent. None of the subjects included in the analysis showed any serious brain pathology on visual examination of these structural brain images.

Calculation of T1-/T2-weighted ratio maps

Both high-resolution T1-weighted image volumes from each subject were realigned using SPM12 to remove any potential variation between scans and averaged. The realignment parameters for separation, smoothing were 2, 3 with a 7th degree B-spline interpolation. The averaged T1-weighted image and the single T2-weighted image were bias- corrected, using the intensity inhomogeneity correction tool implemented in SPM12 (smoothing, 60 mm; regularization, 0.00001 FWHM). After bias correction, the T1-weighted images were co-registered and re-sliced to the corresponding T2-weighted images, maintaining 2D resolution (normalized mutual information; interpolation, 7th degree B-Spline). Using bias-corrected T2-weighted and bias-corrected and re-sliced T1-weighted images, T1-/T2-weighted ratio maps, representing gray and white matter tissue integrity maps, were computed voxel-by-voxel, as described above (Figure 1).

Figure 1.

The graphical representation illustrating the steps that were used to calculate the T1w/T2w ratio maps. The tissue integrity is represented by ‘x’, and variables α₁, α₂ are scaling factors related to bias field inhomogeneity of T1-weighted and T2-weighted images, respectively.

Normalization and smoothing of T1-/T2-weighted ratio maps

T1-/T2-weighted ratio maps of OSA and control subjects were normalized to Montreal Neurological Institute space (1 × 1 × 1 mm³ voxel size), using a modified unified segmentation approach implemented in the SPM12 package. The warping parameters for x, y, z directions were obtained from the bias-corrected and re-sliced T1-weighted images via modified unified segmentation approach and applied to the corresponding ratio maps. The normalized ratio maps were smoothed using an isotropic Gaussian filter (kernel size, 8 mm).

Global brain mask

The averaged and bias-corrected T1-weighted images of each subject were partitioned into gray matter, white matter, and cerebrospinal fluid tissue types using SPM12. The gray and white matter probability maps of all subjects were normalized into MNI space, and averaged to create a global gray matter and a white matter probability map. The averaged gray and white matter probability maps were thresholded (probability threshold of 0.3), and combined to form a whole-brain brain mask.

Histogram calculation

T1-/T2-weighted ratio maps, calculated for the two subjects that were scanned across the two scanners, were first masked for skull and other non-brain tissue using the global brain mask. Histograms were calculated from these masked ratio maps, using the ‘hist’ function implemented in MATLAB. A total of 250 bins were used to estimate the frequency distribution of T1-weighted/T2-weighted ratio values ranging from 1 to 5.

Background images

Averaged high-resolution T1-weighted image of all OSA and control subjects were normalized to MNI space (1 × 1 × 1 mm³, voxel size), using the modified unified segmentation approach. The normalized images of all OSA and control subjects were averaged to create a whole-brain mean image (background image), which was used for structural identification.

Region-of-interest analyses

Region of interest (ROI) analyses were performed to calculate regional brain T1-weighted/T2-weighted ratio values to determine magnitude differences between groups. Using a neuromorphometric atlas (http://www.neuromorphometrics.com), ROIs were created for each brain site, and these ROIs were used to create anatomic-specific masks based on findings from group comparisons. In addition, manual masks were created for white matter sites that were not available in the neuromorphometric atlas using MRIcroN. These anatomic-specific masks were used to compute average values from individual smoothed T1-weighted/T2-weighted ratio maps of OSA and control subjects.

Statistical analyses

The SPM12 and the IBM statistical package for social sciences (SPSS v24) were employed for statistical analyses. Chi-square test and independent-sample t-tests were used to examine group differences in demographic, neuropsychologic, sleep, and cognitive data. The consistency of histograms between two scanners was examined with the Kolmogorov test, and the difference in shapes was assessed with the Wilcoxon rank-sum test.

The smoothed T1-/T2-weighted ratio maps were compared voxel-by-voxel between groups using analysis of covariance (ANCOVA), with age, sex, BMI, PSQI, ESS, BAI, and BDI-II included as covariates. Corrections for multiple comparisons (false discovery rate, p < 0.011) were performed with the REST toolbox.²² The statistical threshold (p < 0.011) was chosen to achieve clusters with maximum significant level corrected for multiple comparisons. The global brain mask was used to restrict the analysis to brain regions only, and sites with significant differences between groups were overlaid onto a background image for structural identification.

The mean ratio values, derived from ROI analyses, were compared between groups using ANCOVA (covariates, age, sex, BMI, PSQI, ESS, BAI, and BDI-II), and effect sizes were calculated. We considered a p value less than 0.05 as statistically significant.

To identify any associations between T1-/T2-weighted ratio values and cognitive, sleep, mood, and anxiety variables, we performed partial correlations between ratio values and BMI, PSQI, ESS, BDI-II, BAI, and MoCA in OSA subjects (covariates, age, and gender).

Results

Demographics and biophysical variables

No significant differences in age (p = 0.58) or gender (p = 0.09) emerged between groups (Table 1). However, OSA subjects showed higher BMI values compared with control subjects (p < 0.001; Table 1).

Physiology, sleep, depression, anxiety, and cognitive variables

No significant difference in heart rate (p = 0.096) appeared between OSA and controls. However, both systolic (p = 0.001) and diastolic (p = 0.01) blood pressure values were significantly elevated in OSA as compared with control subjects (Table 1). Sleep scores (PSQI, p < 0.001; ESS, p < 0.001), as well as BDI-II values (BDI-II, p = 0.045) were significantly higher in OSA compared with controls (Table 1). Although global MoCA scores did not differ significantly between OSA and control subjects, visuospatial scores were significantly lower in OSA subjects. In addition, Trail-A and Trail-B measures did not show any significant differences between OSA and control subjects (Table 1).

Histogram consistency and shape analysis

No significant differences were observed in shape (p = 0.56, Wilcoxon rank-sum test) or consistency (p = 0.46, Kolmogorov test) between the two histograms of T1-/T2-weighted maps of control subjects scanned across the two scanners (Figure 2), indicating that data collected from two different scanners could be pooled for the analysis. The two peaks represent the GM and WM distributions.

Figure 2.

T1-weighted/T2-weighted ratio histograms of two control subjects scanned on two different scanners (3.0-Tesla, Magnetom, Tim-Trio, and Prisma). T1-weighted/T2-weighted ratio histograms were plotted for descalped T1-weighted/T2-weighted ratio maps obtained from Tim-Trio (red) and Prisma (blue) scanners. The x-axis (intensity) represents the bins for T1-w/T2-w ratio values.

T1-/T2-weighted ratio differences

Ratio maps from an OSA and a control subject are shown, overlaid onto a background image (Figure 3). The OSA subject shows lowered brain tissue integrity across substantial areas of the whole-brain compared with control subjects [Figure 3(a) versus Figure 3(a)]. Multiple brain areas in OSA subjects showed lowered T1-/T2-weighted ratio values, indicating altered gray and white matter microstructural differences in those sites, compared with control subjects (Figure 4, p < 0.011, FDR corrected). No brain sites showed increased T1-/T2-weighted ratio values in OSA over controls. Brain sites in OSA subjects that showed lowered ratio values included the bilateral frontal cortices [Figure 4(a)], bilateral insular cortices [Figure 4(b)], posterior corpus callosum [Figure 4(c)], cerebellum and cerebellar peduncles [Figure 4(d)], anterior [Figure 4(e)], mid [Figure 4(f)], and posterior cingulate cortices [Figure 4(g)], bilateral amygdala [Figure 4(h)], fornix fibers [Figure 4(i)], left hippocampus [Figure 4(j)], cerebellar vermis [Figure 4(k)], bilateral occipital white matter [Figure 4(l)], parietal cortices [Figure 4(m)], caudate [Figure 4(n)] and putamen [Figure 4(o)], internal and external capsule [Figure 4(p)], thalamus [Figure 4(q)], basal forebrain [Figure 4(r)], ventral temporal white matter [Figure 4(s)], cerebellar cortices [Figure 4(t)], cerebral crus [Figure 4(u)], medial medulla [Figure 4(v)], basal and caudal pons [Figure 4(w)], bilateral anterior [Figure 4(z)], mid [Figure 4(a1)] and posterior [Figure 4(b1)] corona radiata. Remarkable was absence of differences in tissue integrity in major portions of the thalamus and the frontal cortex, although the ventral frontal cortex was severely affected.

Figure 3.

T1-w/T2-w ratio maps from an (a) OSA and (b) an age and gender-matched control subject. The OSA subject shows altered tissue integrity in large areas of the whole-brain relative to the control subject.

Figure 4.

Brain sites showing lowered T1-w/T2-w ratio values in OSA over control subjects (FDR corrected p < 0.011). These sites include the bilateral frontal cortices (a), bilateral insular cortices (b), posterior corpus callosum (c), cerebellum and cerebellar peduncles (d), anterior (e), mid (f), and posterior cingulate cortices (g), bilateral amygdala (h), fornix (i), hippocampus (j), vermis (k), bilateral occipital white matter (l), parietal cortices (m), caudate (n) and putamen (o), internal and external capsule (p), thalamus (q), basal forebrain (r), ventral temporal white matter (s), cerebellar cortices (t), cerebral crus (u), medial medulla (v), basal and caudal pons (w), anterior (x) and mid corpus callosum (y), bilateral anterior (z), mid (a1) and posterior (b1) corona radiata. Color bar represents t-statistic values (L = Left; R = Right).

ROI analyses

Regional brain T1-/T2-weighted ratio values from various sites of OSA and control subjects are summarized in Table 2. The regional brain mean ratio values were significantly lowered in OSA compared to control subjects, and showed predominantly large effect sizes between groups.

Table 2.

Regional T1-/T2-weighted ratio values of control and OSA subjects.

Brain region	Control(mean ± SD)N = 115	OSA(mean ± SD)N = 106	Cluster Size (voxels)	p values (effect size)
R Insula	0.38 ± 0.06	0.34 ± 0.06	286	<0.001 (0.66)
L Insula	0.48 ± 0.08	0.44 ± 0.08	556	0.001 (0.49)
R Putamen	0.66 ± 0.12	0.60 ± 0.12	3	0.002 (0.50)
L Putamen	0. 72 ± 0.13	0.65 ± 0.14	415	0.001 (0.51)
R Caudate	0.64 ± 0.13	0.61 ± 0.13	237	0.001 (0.23)
L Caudate	0.76 ± 0.14	0.68 ± 0.15	16	0.001 (0.55)
R Thalamus	0.79 ± 0.14	0.71 ± 0.14	473	0.001 (0.55)
L Thalamus	0.79 ± 0.14	0.71 ± 0.15	1183	0.001 (0.55)
L Hippocampus	0.48 ± 0.08	0.43 ± 0.08	734	<0.001 (0.62)
R Hippocampus	0.44 ± 0.07	0.40 ± 0.07	314	0.001 (0.57)
R Amygdala	0.43 ± 0.07	0.38 ± 0.07	198	0.001 (0.71)
L Amygdala	0.44 ± 0.07	0.40 ± 0.07	233	0.001 (0.57)
R Mid Cingulate	0.57 ± 0.10	0.52 ± 0.10	480	0.001 (0.49)
L Mid Cingulate	0.69 ± 0.12	0.62 ± 0.12	94	0.001 (0.58)
R Basal Forebrain	0.45 ± 0.08	0.41 ± 0.09	14	0.001 (0.46)
L Basal Forebrain	0.37 ± 0.06	0.33 ± 0.06	117	0.001 (0.66)
R Medial Orbital gyrus	0.44 ± 0.07	0.39 ± 0.08	1844	<0.001 (0.66)
L Medial Orbital gyrus	0.52 ± 0.08	0.46 ± 0.09	658	<0.001 (0.70)
Cerebellar Crus 1	0.42 ± 0.06	0.39 ± 0.06	1851	0.018 (0.49)
Midline Pons	0.76 ± 0.13	0.71 ± 0.13	3218	0.001 (0.38)
Cerebellar Vermis	0.40 ± 0.06	0.37 ± 0.06	2234	<0.001 (0.50)
Fornix	0.36 ± 0.07	0.33 ± 0.07	330	<0.001 (0.42)
Cerebellar Peduncles	0.85 ± 0.14	0.79 ± 0.15	3809	<0.001 (0.41)
Cerebellar Cortices	0.55 ± 0.08	0.51 ± 0.08	6302	<0.001 (0.50)
Medial Medulla	0.54 ± 0.09	0.49 ± 0.08	1773	<0.001 (0.58)

L, Left; OSA, obstructive sleep apnea; R, Right; SD, standard deviation.

Brain region labels correspond to Figure 2 and mean, SD, and p values are corrected for age, gender, BMI: body mass index; PSQI: Pittsburgh sleep quality index; ESS: Epworth sleepiness scale; BDI-II: Beck Depression Inventory II; and BAI: Beck Anxiety Inventory.

Correlations between T1-/T2-weighted ratios and sleep, mood, anxiety, and cognitive variables in OSA

BAI values showed significant negative correlations between T1-/T2-weighted ratio values of the left hippocampus and right amygdala in OSA subjects. No other variables showed any significant correlations with T1-/T2-weighted ratio values (Table 3).

Table 3.

Correlation r-values (p values) between T1-/T2-weighted ratio values and sleep, mood, anxiety, and cognitive scores.

Brain region	BMIn = 106	PSQIn = 106	ESSn = 106	BDI–IIn = 106	BAIn = 106	MoCAn = 57	AHI	SaO₂ delta
R Insula	–0.174 (0.204)	0.018 (0.897)	0.005 (0.971)	–0.101 (0.463)	–0.186 (0.174)	0.091 (0.508)	–0.019 (0.890)	–0.072 (0.603)
L Insula	–0.151 (0.272)	–0.050 (0.716)	–0.034 (0.803)	–0.185 (0.175)	–0.154 (0.260)	0.076 (0.583)	0.019 (0.892)	0.110 (0.426)
R Putamen	–0.076 (0.583)	–0.043 (0.755)	0.040 (0.775)	–0.196 (0.151)	–0.195 (0.155)	0.006 (0.964)	0.088 (0.525)	0.162 (0.238)
L Putamen	–0.183 (0.181)	–0.088 (0.521)	0.002 (0.990)	–0.204 (0.135)	–0.155 (0.259)	0.066 (0.630)	0.056 (0.684)	0.135 (0.324)
R Caudate	–0.172 (0.208)	–0.122 (0.373)	–0.001 (0.997)	–0.231 (0.090)	–0.186 (0.175)	0.015 (0.912)	0.116 (0.398)	0.159 (0.248)
L Caudate	–0.024 (0.238)	–0.016 (0.905)	–0.030 (0.829)	–0.212 (0.120)	–0.135 (0.327)	0.055 (0.689)	0.070 (0.612)	0.113 (0.412)
R Thalamus	–0.059 (0.670)	0.012 (0.933)	0.018 (0.896)	–0.171 (0.211)	–0.074 (0.590)	0.057 (0.679)	0.065 (0.636)	0.115 (0.402)
L Thalamus	0.010 (0.941)	–0.072 (0.600)	–0.090 (0.515)	–0.174 (0.203)	–0.236 (0.083)	0.009 (0.947)	0.042 (0.760)	0.154 (0.262)
R Hippocampus	–0.026 (0.851)	–0.054 (0.697)	–0.041 (0.764)	–0.150 (0.276)	–0.199 (0.145)	0.065 (0.637)	0.094 (0.494)	0.201 (0.141)
L Hippocampus	–0.019 (0.888)	–0.085 (0.539)	–0.129 (0.348)	–0.222 (0.104)	–0.281 (0.037)	–0.012 (0.932)	0.015 (0.913)	0.136 (0.324)
R Amygdala	–0.036 (0.792)	–0.086 (0.533)	–0.105 (0.447)	–0.245 (0.072)	–0.306 (0.023)	0.041 (0.764)	0.011 (0.939)	0.170 (0.213)
L Amygdala	–0.110 (0.422)	–0.080 (0.562)	–0.004 (0.977)	–0.178 (0.195)	–0.144 (0.295)	0.087 (0.527)	0.069 (0.618)	0.078 0.571
R Mid Cingulate	–0.066 (0.634)	–0.044 (0.747)	0.127 (0.356)	–0.162 (0.238)	–0.065 (0.638)	0.054 (0.693)	0.063 (0.649)	0.127 (0.356)
L Mid Cingulate	–0.101 (0.463)	–0.048 (0.729)	–0.080 (0.564)	–0.130 (0.343)	–0.203 (0.137)	0.025 (0.856)	0.009 (0.947)	0.130 (0.345)
R Basal Forebrain	–0.092 (0.505)	–0.023 (0.870)	0.014 (0.921)	–0.053 (0.703)	–0.209 (0.125)	0.050 (0.719)	0.039 (0.775)	0.150 (0.273)
L Basal Forebrain	–0.075 (0.585)	–0.074 (0.593)	0.054 (0.694)	–0.210 (0.123)	–0.207 (0.130)	0.240 (0.078)	0.147 (0.285)	0.127 (0.355)
R Medial Orbital gyrus	–0.043 (0.753)	–0.064 (0.644)	0.086 (0.530)	–0.199 (0.145)	–0.201 (0.140)	0.106 (0.442)	0.101 (0.465)	0.144 (0.296)
L Medial Orbital gyrus	0.030 (0.828)	0.044 (0.747)	0.036 (0.797)	–0.176 (0.199)	–0.120 (0.384)	0.152 (0.268)	0.139 (0.313)	0.210 (0.124)
Cerebellar Crus 1	–0.080 (0.559)	–0.074 (0.590)	–0.058 (0.672)	–0.242 (0.075)	–0.217 (0.112)	0.038 (0.784)	0.061 (0.659)	0.105 (0.443)
Midline Pons	–0.125 (0.365)	–0.129 (0.347)	–0.053 (0.702)	–0.163 (0.234)	–0.183 (0.182)	0.070 (0.610)	0.077 (0.574)	0.056 (0.685)
Cerebellar Vermis	–0.117 (0.395)	–0.031 (0.821)	0.021 (0.877)	–0.162 (0.236)	–0.139 (0.312)	0.012 (0.930)	0.045 0.743	0.124 0.368
Fornix	–0.053 (0.702)	–0.040 (0.770)	–0.070 (0.612)	–0.178 (0.193)	–0.161 (0.241)	0.048 (0.728)	0.078 (0.572)	0.108 (0.431)
Cerebellar Peduncles	–0.140 (0.308)	–0.137 (0.317)	–0.066 (0.635)	–0.075 (0.585)	–0.210 (0.123)	0.155 (0.257)	0.081 (0.558)	0.023 (0.870)
Cerebellar Cortices	–0.045 (0.744)	–0.003 (0.981)	–0.038 (0.785)	–0.178 (0.193)	–0.141 (0.305)	0.045 (0.746)	0.066 (0.634)	0.117 (0.397)
Medial Medulla	–0.125 (0.365)	–0.129 (0.347)	–0.053 (0.702)	–0.163 (0.234)	–0.183 (0.182)	0.070 (0.610)	0.077 (0.574)	0.056 (0.685)

AHI, Apnea-Hypopnea-Index; BAI, Beck anxiety inventory; BDI-II, Beck depression Inventory II; BMI, Body mass index; ESS, Epworth sleepiness scale; MoCA, Montreal cognitive assessment; PSQI, Pittsburgh sleep quality index; SaO₂, Oxygen saturation.

p values are corrected for age and gender.

Discussion

Overview

Tissue injury appeared in several white matter and subcortical gray matter regions in recently diagnosed, untreated OSA subjects, as indicated by significantly lowered T1-/T2-weighted ratio values. Damage appeared predominantly in the medial medulla, extending to the midline pons, regions which contribute significantly to breathing and sympathetic control. In addition, areas such as cerebellar cortex, hippocampus, and amygdala, regions important for muscle coordination, and mediating cognitive, memory, and autonomic aspects showed lowered tissue integrity, along with surrounding white matter. Fiber projections to the frontal cortex, as well as insular and cingulate cortices necessary for autonomic regulation, showed significantly lowered tissue integrity. Previous studies have shown several brain sites in OSA to be compromised, but those sites were not as widespread as shown here. The widespread white matter damage to cortical-subcortical circuitry may result from cellular damage from hypoxic or inflammatory processes accompanying the sleep disordered breathing condition, and could contribute to understanding neural dysfunctions that lead to lowered autonomic, cognitive, and mood functions in OSA subjects. Lowered T1-/T2-weighted ratio values preferentially emerged on the left side. The lateralization of injury follows a pattern in a number of OSA studies.^6,7,9 The mechanisms underlying the asymmetry are unclear, but may include relative perfusion differences between the left and right brain.

T1-weighted, T2-weighted images, and microstructural tissue differences

Conduction impairment from tissue damage can cause paralysis, sensory-motor dysfunction, cognitive deficits, mental retardation, and lead to compromised vital functions.²³ Important for OSA, slowing of conduction of signals to brainstem blood pressure and breathing arriving from cerebellar and forebrain sites can lead to delayed timing of activation of upper airway muscles relative to the diaphragm, accentuating airway obstruction. Moreover, late activation of transient blood pressure increases can lead to syncope during movement and slowed cognitive responses. We showed that the T1-weighted/T2-weighted ratio maps could be used to detect altered gray and white matter microstructural tissue differences in OSA. White matter structures, such as midline and surrounding occipital white matter,²⁴ extending to parietal and motor cortices, showed lowered ratio values in OSA subjects. This may be due to the increase in T2-weighted signal, especially in white matter areas, as a result of hypoxia accompanying the condition.²⁵ As the T2-weighted signal increases, the overall ratio values lower. In addition, frontal cortices also showed compromised tissue integrity, indicating that this procedure is robust in identifying various regions affected by microstructural tissue differences.

OSA findings and their implications

The underlying mechanisms contributing to gray and white matter microstructural differences in OSA may be a combination of hypoxic, hypercarbic, excitotoxic, and perfusion issues accompanying the syndrome. Upper airway blockage during sleep in OSA results in intermittent hypoxia, a process that reduces O₂ supply to brain tissue, and then rapidly restores oxygenation a repetitive sequence especially damaging to tissue.^26,27 The hypoxic processes initiate water movement from extracellular to intracellular spaces, leading to axonal and neuronal inflammation, including myelin swelling in acute conditions. The inflammatory responses can be observed as white matter differences, gray matter volume alterations, and metabolic deficits, and can be examined by various MRI and spectroscopy procedures in OSA, which are present in multiple brain areas.^5,9,28–30 In addition, white matter injury in the form of loss of fiber integrity has been shown by fraction anisotropy techniques³¹ in regions such as cingulate, corpus callosum, and multiple cortical projection areas. However, the findings here show widespread altered tissue integrity not only in the regions specific to autonomic, respiratory, mood, and cognitive roles, but also in surrounding white matter sites. Lowered T1-/T2-weighted ratio values appeared in frontal, insular, and cingulate cortices. Abnormalities in these sites can establish the basis for affective disorders, including depression and anxiety,^32,33 and may affect autonomic regulation.^34–36 White matter differences also extended to the cingulum bundle and cerebellar areas that control sympathetic and parasympathetic nervous systems.^37–39 The cerebellar and cingulum bundle findings may contribute to the distorted timing to autonomic challenges (e.g. delayed onset of heart rate changes)⁴⁰ and excessive sympathetic tone.^41,42

Areas implicated in normal emotional behavior have been closely related to the medial prefrontal cortex (MPFC), insula, amygdala, hippocampus, and ventromedial parts of the basal ganglia.^43,44 We found altered tissue integrity in white matter circuits between the orbital and medial prefrontal cortex, amygdala, hippocampus, continuing up to the fornix and mid cingulate. Conduction failure between these regions can compromise emotional and mood behavior in OSA subjects,⁴⁵ which is also evident from the clinical measures, as they differ between groups. Correlation analyses show significant associations between BAI and T1-/T2-weighted ratio values of the insula, hippocampus, amygdala, basal forebrain, and cerebellar vermis. The majority of these areas, including the amygdala, insula, hippocampus, and basal forebrain are involved in anxiety control, a common issue in OSA subjects (Table 3). In addition, tissue injury to these limbic pathways can also cause deficits in processing speed, working and episodic memory, and worsen the quality of life in OSA subjects, as readily observed in major depressive disorder cases.⁴⁶

Similarly altered tissue integrity was observed in basal forebrain and pathways projecting to the caudate, putamen, and hippocampus and can severely affect mental performance and learning.⁴⁷ The mid corpus callosum, along with surrounding white matter, leading to the cingulate cortices, showed impaired tissue integrity that can affect inter-hemispheric transfer of information at various visuomotor and cognitive levels.⁴⁸ Although there was no significant difference in global MoCA scores, it is evident that this injury could impact visuospatial functioning, as a significant difference in visuospatial scores emerged between groups. In addition, tissue injury appeared in white matter projections from dorsal and ventral temporal regions, passing through the thalamus and frontal lobule. Injury to this circuitry is known to cause severe cognitive deficits.^49,50 Finally, subcortical structures such as the midline pons, and white matter tracts extending to the medial medulla showed lowered ratio values as well. Damage to these areas and surrounding white matter structures would compromise regulation of blood pressure, respiration, and integration of baroreceptor and chemoreceptor afferents.^51–54

We found widespread altered gray and white matter tissue integrity in regions responsible for cognitive, mood, autonomic, and respiratory functions. Several sites that showed more widespread tissue differences, included the corpus callosum, cingulum bundle, frontal, parietal, temporal, and occipital white matter, cerebellar peduncles, and anterior, mid, and posterior corona radiate over findings described earlier.^7,9,31 These data suggest that aggressive breathing intervention is needed as soon as possible in OSA subjects after initial diagnosis.

Limitations

Several limitations of this study should be acknowledged. The spatial resolution of T2-weighted images (slice thickness of 3.5, 4 mm) are lower, compared with T1-weighted counterparts, which may have affected co-registration, and thus, tissue-related contrast. Due to time and technical constraints during scanning, we were unable to acquire high-resolution T2-weighted data. However, we co-registered T1-weighted to T2-weighted images, and resampled T1-weighted to T2-weighted resolution, which may have minimized variation in tissue contrast. Furthermore, we smoothed the data by 8 mm, which may affect the resolution in the slice direction, but this is an essential step for voxel-based procedures. In addition, the receiver coil bias across the two different scanners differs, and dividing T1-weighted by T2-weighted images may have not completely eradicated these inhomogeneities. Thus, variations in the bias field may have affected our findings, which could be improved if these issues are corrected in future studies. Data collected from two different scanners were pooled in this study without including scanner as a covariate, though histogram obtained from same control subjects from two different scanners showed no differences in shape and consistency. Although the T1-weighted and T2-weighted images were acquired successively, they may have undergone different receiver gain calibrations along with different shimming and auto adjustments that might have affected our results. However, histogram analyses of all the control subjects showed consistent pattern with minimal inter-subject variability across dataset, indicating minimum influence from receive gain changes, if any. Significant differences emerged in BMI between groups. Although BMI was used as a covariate in our statistical analysis, we cannot rule out any impact from BMI on these findings. The microstructural tissue differences found in this study could be enhanced due to multiple factors, including contribution from comorbidities, such as hypertension and type 2 diabetes. Therefore, future studies performing logistic regression to predict these changes by different comorbidities would be helpful. Factors including blood pressure, heart rate, and cholesterol levels have been known to influence cerebral microstructure. However, non-availability of such information from some OSA participants hindered us to account for these contributing factors in the statistical analyses.

Conclusion

Regional brain T1-weighted/T2-weighted ratio values are significantly lowered, indicating differences in gray and white matter microstructural tissue in multiple areas in newly diagnosed, treatment-naive OSA over control subjects. These areas of lowered ratio values include regions of critical autonomic, somatic motor inhibition, respiratory, mood, and cognitive control sites, including the insular cortices and surrounding white matter, hippocampus, fornix fibers, amygdala, basal forebrain fibers, mid-cingulate and cingulum bundle, corpus callosum, cerebellar cortices and peduncles, and midline and caudal pons. Some projections, particularly in thalamic sites, were less affected. Altered brain tissue in OSA is more widespread than determined previously, indicating that aggressive interventions are required in newly diagnosed OSA subjects. A range of pathological processes, including hypoxia and ischemia, repeated arousals, hypercarbia, and inflammatory-induced mechanisms may contribute to these differences. The findings suggest that T1-/T2-weighted ratio procedures could be useful procedures to examine gray and white matter microstructural differences with disease progression and help understand the OSA brain.

Footnotes

Acknowledgements

We thank Ms. Rebecca Harper, Ms. Karan Harada, and Ms. Cristina Cabrera-Mino for their help with data collection.

Declarations

ORCID iD

Rajesh Kumar

References

Sforza

de Saint Hilaire

Pelissolo

, et al. Personality, anxiety and mood traits in patients with sleep-related breathing disorders: effect of reduced daytime alertness. Sleep Med 2002; 3: 139–145.

Rosenzweig

Glasser

Polsek

, et al. Sleep apnoea and the brain: a complex relationship. Lancet Respir Med 2015; 3: 404–414.

Patil

Schneider

Schwartz

, et al. Adult obstructive sleep apnea: pathophysiology and diagnosis. Chest 2007; 132: 325–337.

Lurie

Obstructive sleep apnea in adults: epidemiology, clinical presentation, and treatment options. Adv Cardiol 2011; 46: 1–42.

Kumar

Farahvar

Ogren

, et al. Brain putamen volume changes in newly-diagnosed patients with obstructive sleep apnea. NeuroImage Clin 2014; 4: 383–391.

Tummala

Roy

Park

, et al. Associations between brain white matter integrity and disease severity in obstructive sleep apnea. J Neurosci Res 2016; 94: 915–923.

Kumar

Chavez

Macey

, et al. Altered global and regional brain mean diffusivity in patients with obstructive sleep apnea. J Neurosci Res 2012; 90: 2043–2052.

Iwatani

Ishida

Donishi

, et al. Use of T1-weighted/T2-weighted magnetic resonance ratio images to elucidate changes in the schizophrenic brain. Brain Behav 2015; 5: e00399.

Kumar

Pham

Macey

, et al. Abnormal myelin and axonal integrity in recently diagnosed patients with obstructive sleep apnea. Sleep 2014; 37: 723–732.

10.

Melie-Garcia

Slater

Ruef

, et al. Networks of myelin covariance. Hum Brain Mapp 2018; 39: 1532–1554.

11.

Coremans

Luypaert

Verhelle

, et al. A method for myelin fiber orientation mapping using diffusion-weighted MR images. Magn Reson Imaging 1994; 12: 443–454.

12.

Alonso-Ortiz

Levesque

Pike

GB.

Multi-gradient-echo myelin water fraction imaging: comparison to the multi-echo-spin-echo technique. Magn Reson Med 2018; 79: 1439–1446.

13.

Prasloski

Rauscher

MacKay

, et al. Rapid whole cerebrum myelin water imaging using a 3D GRASE sequence. NeuroImage 2012; 63: 533–539.

14.

Castronovo

Scifo

Castellano

, et al. White matter integrity in obstructive sleep apnea before and after treatment. Sleep 2014; 37: 1465–1475.

15.

Martínez-García

Soler-Cataluña

Ejarque-Martínez

, et al. Continuous positive airway pressure treatment reduces mortality in patients with ischemic stroke and obstructive sleep apnea: a 5-year follow-up study. Am J Respir Crit Care Med 2009; 180: 36–41.

16.

Glasser

Van Essen

DC.

Mapping human cortical areas in vivo based on myelin content as revealed by T1- and T2-weighted MRI. J Neurosci 2011; 31: 11597–11616.

17.

Ganzetti

Wenderoth

Mantini

Whole brain myelin mapping using T1- and T2-weighted MR imaging data. Front Hum Neurosci 2014; 8: 671.

18.

Knutson

Rathouz

Yan

, et al. Stability of the Pittsburgh Sleep Quality Index and the Epworth Sleepiness Questionnaires over 1 year in early middle-aged adults: the CARDIA study. Sleep 2006; 29: 1503–1506.

19.

Beck

Epstein

Brown

, et al. An inventory for measuring clinical anxiety: psychometric properties. J Consult Clin Psychol 1988; 56: 893–897.

20.

Beck

Steer

Ball

, et al. Comparison of Beck Depression Inventories -IA and -II in psychiatric outpatients. J Pers Assess 1996; 67: 588–597.

21.

Rorden

Karnath

Bonilha

Improving lesion-symptom mapping. J Cogn Neurosci 2007; 19: 1081–1088.

22.

Song

Dong

Z-Y

Long

X-Y

, et al. REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One 2011; 6: e25031.

23.

Fields

RD.

White matter in learning, cognition and psychiatric disorders. Trends Neurosci 2008; 31: 361–370.

24.

Fuster

JM.

Network memory. Trend Neurosci 1997; 20: 451–459.

25.

Hunt

Jr Theilmann

Smith

, et al. Cerebral diffusion and T(2): MRI predictors of acute mountain sickness during sustained high-altitude hypoxia. J Cereb Blood Flow Metab 2013; 33: 372–380.

26.

Rukhadze

Fenik

Benincasa

, et al. Chronic intermittent hypoxia alters density of aminergic terminals and receptors in the hypoglossal motor nucleus. Am J Resp Crit Care Med 2010; 182: 1321–1329.

27.

Harper

Macey

Henderson

, et al. fMRI responses to cold pressor challenges in control and obstructive sleep apnea subjects. J Appl Physiol (1985) 2003; 94: 1583–1595.

28.

Yadav

Kumar

Macey

, et al. Insular cortex metabolite changes in obstructive sleep apnea. Sleep 2014; 37: 951–958.

29.

Macey

Henderson

Macey

, et al. Brain morphology associated with obstructive sleep apnea. Am J Respir Crit Care Med 2002; 166: 1382–1387.

30.

Kumar

Birrer

BVX

Macey

, et al. Reduced mammillary body volume in patients with obstructive sleep apnea. Neurosci Lett 2008; 438: 330–334.

31.

Macey

Kumar

Woo

, et al. Brain structural changes in obstructive sleep apnea. Sleep 2008; 31: 967–977.

32.

Kumar

Macey

Cross

, et al. Neural alterations associated with anxiety symptoms in obstructive sleep apnea syndrome. Depress Anxiety 2009; 26: 480–491.

33.

Cross

Kumar

Macey

, et al. Neural alterations and depressive symptoms in obstructive sleep apnea patients. Sleep 2008; 31: 1103–1109.

34.

Marins

Limborço-Filho

Xavier

, et al. Functional topography of cardiovascular regulation along the rostrocaudal axis of the rat posterior insular cortex. Clin Exp Pharmacol Physiol 2016; 43: 484–493.

35.

Oppenheimer

Cechetto

The insular cortex and the regulation of cardiac function. Compr Physiol 2016; 6: 1081–1133.

36.

Park

Palomares

Woo

, et al. Aberrant insular functional network integrity in patients with obstructive sleep apnea. Sleep 2016; 39: 989–1000.

37.

Cao

Huang

Jiang

, et al. Cerebellar fastigial nuclear glutamatergic neurons regulate immune function via hypothalamic and sympathetic pathways. J Neuroimmune Pharmacol 2015; 10: 162–178.

38.

Mendez

Joshi

Daianu

, et al. White matter changes associated with resting sympathetic tone in frontotemporal dementia vs Alzheimer’s disease. PLoS ONE 2015; 10: e0142445.

39.

de Morree

Szabó

Rutten

, et al. Central nervous system involvement in the autonomic responses to psychological distress. Neth Heart J 2013; 21: 64–69.

40.

Macey

Kumar

Woo

, et al. Heart rate responses to autonomic challenges in obstructive sleep apnea. PLoS ONE 2013; 8: e76631.

41.

Henderson

Macey

, et al. Brain responses associated with the Valsalva maneuver revealed by functional magnetic resonance imaging. J Neurophysiol 2002; 88: 3477–3486.

42.

Somers

Dyken

Clary

, et al. Sympathetic neural mechanisms in obstructive sleep apnea. J Clin Invest 1995; 96: 1897–1904.

43.

Drevets

Price

Furey

ML.

Brain structural and functional abnormalities in mood disorders: implications for neurocircuitry models of depression. Brain Struct Funct 2008; 213: 93–118.

44.

Sullivan

Neale

Kendler

KS.

Genetic epidemiology of major depression: review and meta-analysis. Am J Psychiatry 2000; 157: 1552–1562.

45.

Hastings

Parsey

Oquendo

, et al. Volumetric analysis of the prefrontal cortex, amygdala, and hippocampus in major depression. Neuropsychopharmacology 2004; 29: 952–959.

46.

Weber

Giannakopoulos

Delaloye

, et al. Volumetric MRI changes, cognition and personality traits in old age depression. J Affect Disord 2010; 124: 275–282.

47.

Kujala

Portin

Ruutiainen

The progress of cognitive decline in multiple sclerosis. Brain 1997; 120(Pt2): 289–297.

48.

Schulte

Muller-Oehring

EM.

Contribution of callosal connections to the interhemispheric integration of visuomotor and cognitive processes. Neuropsychol Rev 2010; 20: 174–190.

49.

Rodriguez-Cruces

Concha

White matter in temporal lobe epilepsy: clinico-pathological correlates of water diffusion abnormalities. Quant Imaging Med Surg 2015; 5: 264–278.

50.

Bunce

Anstey

Cherbuin

, et al. Cognitive deficits are associated with frontal and temporal lobe white matter lesions in middle-aged adults living in the community. PLoS ONE 2010; 5: e13567.

51.

Gozal

Daniel

Dohanich

GP.

Behavioral and anatomical correlates of chronic episodic hypoxia during sleep in the rat. J Neurosci 2001; 21: 2442–2450.

52.

Pae

Chien

Harper

RM.

Intermittent hypoxia damages cerebellar cortex and deep nuclei. Neurosci Lett 2005; 375: 123–128.

53.

Hatam

Rasoulpanah

Nasimi

GABA modulates baroreflex in the ventral tegmental area in rat. Synapse 2015; 69: 592–599.

54.

van Den Buuse

Catanzariti

. Stimulation of the ventral tegmental area enhances the effect of vasopressin on blood pressure in conscious rats. Br J Pharmacol 2000; 129: 29–36.