ARTICLE

fMRI responses to cold pressor challenges in control and obstructive sleep apnea subjects

Abstract

Obstructive sleep apnea (OSA) patients exhibit altered sympathetic outflow, which may reveal mechanisms underlying the syndrome. We used functional MRI (fMRI) in 16 control and 10 OSA subjects who were free of cardiovascular or mood-altering drugs to examine neural responses to a forehead cold pressor challenge, which elicits respiratory slowing, bradycardia, and enhanced sympathetic outflow. The magnitude of cold-induced bradycardia was smaller, and respiratory slowing showed greater intersubject variability and reached a nadir later in OSA patients. Both groups showed similar signal changes to cold stimulation in multiple brain sites. However, signal increases emerged in OSA over controls in anterior and posterior cingulate and cerebellar and frontal cortex, whereas signals markedly declined in the ventral thalamus, hippocampus, and insula rather than rising as in controls. Anomalous responses often paralleled changes in breathing and heart rate. Medullary, midbrain areas and lentiform and cerebellar dentate nuclei also showed lower signals in OSA cases. Cold pressor physiological responses are modified in OSA and may result from both diminished and exaggerated responses in multiple brain structures.

obstructive sleep apnea (osa) patients exhibit repetitive periods of upper airway muscle atonia in the presence of enhanced diaphragmatic efforts during sleep. The resultant cessation of airflow from a collapsed airway induces enhanced venous return, large transient elevations in arterial pressure, and a sequence of tachycardia and bradycardia associated with significant alterations in sympathetic nervous system outflow (30). Anatomic characteristics of patients, including restricted airway dimensions resulting from deviated septa, macroglossia, micrognathia, hypertrophied tonsils, or obesity, play a significant role in the genesis of the syndrome (32). However, several features of OSA implicate a role for central neural dysfunction in initiating, maintaining, or worsening the physiological characteristics. The maintenance of upper airway atonia while diaphragmatic musculature receives enhanced drive suggests a centrally mediated loss of integration of airway or cardiovascular afferent activity with efferent flow to the oral musculature. The findings that OSA prevalence increases with age (39), an effect not related to obesity, and that moderate sleep apnea often progresses to severe apnea without significant changes in anatomic features (41) suggest that neural, perhaps including degenerative, changes play a role in OSA. Furthermore, atrial pacing reduces obstructed events (11); the mechanisms presumably involve central integration of cardiovascular afferent signals with upper airway motor control. Patients with OSA show significant gray matter loss in cortical regions mediating sensory and motor control of the upper airway, as well as areas involved in respiratory patterning and blood pressure control, such as limbic sites and the cerebellum (23). It is unclear whether the gray matter loss found in OSA subjects preceded the onset of symptoms, or whether damage was a consequence of repetitive hypoxic episodes associated with apneic events. In either case, the neural morphological changes have the potential to exacerbate characteristics of OSA and to result in further deterioration of symptoms. Examination of functional responses of neural sites mediating integration of sensory, motor, and cardiovascular integration in OSA may reveal the nature of aberrant function.

The reflexive interactions between the cardiovascular and respiratory systems provide a potential means for dysfunctional cardiovascular control to induce upper airway atonia or alter appropriate timing of muscle action. Transient blood pressure elevation normally diminishes respiratory muscle activity and slows breathing (15), an effect that is preferentially enhanced in upper airway musculature (25). The occurrence of enhanced rather than diminished diaphragmatic effort that is out of synchrony with oral muscle activity in OSA suggests disorders of neural reflexive timing or recruitment mechanisms for both blood pressure and somatic breathing actions in afflicted patients. Responses to other autonomic challenges that significantly modify blood pressure, including the Valsalva maneuver, are deficient in patients with OSA (14, 47). The presence of aberrant autonomic reflex expression has the potential to contribute to the disordered action of upper airway musculature in the syndrome.

An autonomic reflex that exerts considerable influence on breathing is the cold pressor “dive reflex,” elicited by application of a cold stimulus to the forehead. Such stimulation exerts profound involuntary effects on blood pressure and breathing, including slowed breathing rate, marked blood pressure elevation, and bradycardia. The challenge has been used to examine the integrity of autonomic nervous system components in a number of disorders, including OSA (45); data from this latter study suggest an altered bradycardic response to the cold application in a subset of patients. Using functional MRI, which allows noninvasive visualization of neural responses, we examined neural responses to a cold stimulus to the forehead in OSA and control subjects. The objective was to outline the normal physiological and neural responses to a cold pressor challenge, to determine whether neural responses to this involuntary autonomic challenge differed in OSA patients, and whether any brain sites showing differences overlapped those previously reported as showing gray matter loss in OSA (23).

METHODS

Subjects.

Twenty-one male patients with a confirmed sleep laboratory diagnosis of OSA and 21 healthy control men participated. Subjects were matched by age. Scanner size limitations or claustrophobia precluded inclusion of some OSA patients, including several women. An earlier analysis of all cases showed that pharmacological agents contaminated both the neural and physiological responses. Thus only data from drug-free cases are described here. Subjects were excluded who were currently treated with cardiovascular-related medications (angiotension conversion enzyme inhibitors, diuretics, β-blockers, α-blockers, angiotension blockers) or mood-altering drugs (e.g., serotonin reuptake inhibitors). Data from subjects in whom significant artifactual signals rising from tissue-cerebrospinal fluid or tissue-bony interfaces that were beyond the maximal 6% change expected by normal blood oxygen level-dependent processes were also excluded.

Of the initial 21 OSA subjects, 10 drug-free cases remained [mean age (range): 46 ± 12 (SD) yr (28–64 yr); and mean body mass index: 31 ± 6 kg/m2 (24–41 kg/m2)], and, of the 21 control subjects, 16 were used for further analysis [mean age: 47 ± 10 yr (30–60 yr); body mass index: 27 ± 3 kg/m2 20–33 kg/m2)].

OSA severity was diagnosed according to standard practices (1) as mild (n = 1), mild-moderate (n = 1), moderate (n = 2), moderate-severe (n = 2), or severe (n = 4). Mean apnea/hypopnea index was 38 ± 27 (range: 8–95), and the mean respiratory disturbance index was 38 ± 24 (range: 19–95). The interval between OSA diagnosis and the study ranged from 0 to 6 yr (mean: 16 ± 21 mo). Three controls and two OSA subjects were lefthanded. The research protocol was approved by the Institutional Review Board. All subjects gave their written, informed consent.

Physiological recording.

Electrophysiological recordings were continuously collected during scanning. An ECG was obtained by using magnetic resonance (MR)-compatible surface electrodes, and oxygen saturation was collected by using a Nonin oximeter with a finger transducer (Nonin Medical, Plymouth, MN). Low-noise amplifiers, constructed to operate in high-magnetic field environments, were used to provide signal gain and initial filtering with optical transfer of signals outside the scanner room (38). A measure of respiration was provided by pressure changes from thoracic wall movement by using an air-filled bag placed over the thoracic wall. These pressure signals were conveyed by low-compliance tubing to a transducer outside the scanner room. Arterial pressure readings were obtained immediately before and immediately after each scanning period by using an automated sphygmomanometer. All physiological signals were acquired at 1 kHz on a portable computer with a Quatech analog-to-digital converter (Quatech, Akron, OH). Breathing and heart rates were determined from breath-to-breath and cardiac interbeat intervals, which were calculated from output of a peak-searching algorithm applied to thoracic wall movement and ECG signals, respectively. The heart and breathing rates were filtered by using a median filter with window lengths of 30 and 60 s, respectively. Thresholds for breathing rate changes were determined from the first significant change in rate, as determined by a repeated-measures analysis of variance, and values for the nadir were similarly calculated.

Subsequent to the imaging studies, beat-by-beat blood pressure was collected from a subset of control (n = 6) and OSA (n = 8) subjects by using an identical challenge outside the scanner with a noninvasive blood pressure monitor (Colin 7000; Colin Medical Instruments). Subjects lay in a comparable position, with scanner sounds played through earphones. Data were recorded in a similar fashion as in imaging trials.

MR imaging.

Images were collected with a 1.5-T MR scanner (General Electric Signa EchoSpeed System, Milwaukee, WI). Foam was placed on either side of the head to minimize head movement during scans, and masking tape was applied to the forehead and to the headrest to further minimize motion. To explore rapid signal changes, the gradient-echo, echo-planar procedure, with blood oxygen level-dependent contrast (31), was used for image acquisition. Two series of 25 echo-planar image volumes, comprising 20 oblique sections, were acquired during 150-s time periods. The first series was a baseline with no challenge. The second series consisted of an initial 60-s baseline (10 volumes) followed without pause by a 90-s challenge period (15 volumes). Image parameters were as follows: repetition time (TR) = 6 s per image set, echo time = 60 ms, flip angle = 90°, field of view = 30 × 30 cm, no interslice gap, voxel size = 2.3 × 2.3 × 5.0 mm thick. A series of T1-weighted spin echo images (TR = 500 ms, echo time = 9 ms, field of view = 30 × 30 cm, no interslice gap, voxel size = 1.2 × 1.2 × 5.0 mm thick) was then collected at the same slice locations as the echo-planar images to assist anatomic identification.

The challenge consisted of application of a cold (4°C) bag of deuterium oxide (D2O, heavy water) to the forehead. D2O, rather than water (H2O), was used as a cold medium to reduce the MRI signal artifact produced by liquid water in the bag.

Early in the study, in a subset of subjects, the bag was inadvertently dropped on a forehead area covered with a bandage, which insulated the subject from the cold stimulus. These subjects showed no cardiovascular or respiratory responses and none of the neural signal changes to the challenge. The data were excluded from analysis of cold effects, but served as a partial control for stimulus effects of bag application.

Data analysis.

The first scan (6 s) of each series was excluded from analysis due to signal saturation. The scan representing the transition from baseline to application of the cold challenge was also removed to avoid artifactual contamination with the rapid introduction of the D2O bag near the brain. The remaining 23 image challenge and 24 image baseline sets were evaluated, in part by using SPM99 software (8). SPM99 was used for all of the preprocessing steps, except for intensity normalization (see below). Image sets were initially corrected for slice timing and motion (52) and were then spatially normalized so that each image volume in the time series was in the standardized three-dimensional Montreal Neurological Institute (MNI) reference space (26).

Spatial normalization was performed in two steps. Before normalization, the echo-planar images and T1 images were manually coregistered. In the first step, echo-planar images were normalized to an echo-planar imaging template provided by the MNI. The resulting normalization parameters were applied to the registered T1 image. In a second step, the T1 images from the first normalization step were normalized to the MNI template, and the resulting parameters were applied to the echo-planar images obtained in the first step. This process was necessary due to signal drop-out in the echo-planar images. Gray matter was segmented from white matter components based on image intensity and a priori information on cellular distribution; this procedure was undertaken to remove potential global differences between axons vs. gray matter during the pressor challenges. Images were then Gaussian smoothed (full-width at half-maximum = 8 mm) and intensity normalized; the mean signal intensity from the entire gray matter of the brain at each time point was calculated, and the images scaled so that the mean image intensity was set to an arbitrary level of 500.

Two analysis procedures were employed to determine group and time course patterns in the neural signals. First, regions of interest from anatomic sites known to play significant roles in breathing and cardiovascular control were defined manually on a subject-by-subject basis, and signal patterns from these areas were evaluated over the entire baseline and challenge periods and assessed with repeated-measures ANOVA by using custom software. Averaged trends were plotted for each area. The selection of areas for scrutiny included regions classically associated with autonomic patterning [amygdala, insula, dorsal, medial, and ventral medulla (5a, 13, 18, 19, 33, 42)], as well as regions only recently implicated in exerting a direct or indirect influence on autonomic outflow or motoric-breathing characteristics that could participate in the cold pressor challenge [hippocampus, lentiform nuclei, head of caudate, ventral and dorsal pons, dorsal and ventral midbrain, dentate and fastigial nuclei (13, 40, 53)]. A second “cluster” analysis compared baseline and cold pressor conditions on a voxel-by-voxel basis, with changes in signal intensity modeled with the average heart rate trends. SPM was used for the cluster analysis. A second-level population analysis was used, as this procedure was most effective in finding regions of consistent responses across all subjects compared with a single-level analysis. Extraction of the time trends was performed with custom software. The baseline included the baseline series and the first 60 s of the challenge series. Averaged signal intensity was determined from clusters showing significant differences on the statistical maps at each time point. Time courses of these averaged values were then plotted by using custom software. This second analysis was designed to evaluate differences in the group responses to the challenge while considering a principal physiological accompaniment of the challenge, the bradycardia associated with cold application. After review for potential warping of images relative to anatomic slices, significant voxel changes were mapped onto spatially normalized mean T1-weighted anatomic images. Particular attention was paid to distortion characteristics with the potential to introduce artifacts from inhomogeneous fields or other error sources, such as the vasculature, bone-tissue interfaces, or tissue-cerebrospinal fluid junctures. Maximum temporal resolution was limited by the TR of 6 s.

RESULTS

Physiology.

Application of the cold pressor challenge resulted in heart rate slowing in both control (nadir = −4.8 ± 0.9%) and OSA subjects (nadir = −4.0 ± 1.1%); heart rate returned toward baseline earlier in the OSA group (Fig.1A). Breathing rates declined in controls (nadir = −8.3 ± 3.3%); in OSA patients, breathing rates also declined (nadir = −4.9 ± 12.4%), but with significantly increased intersubject variability (Fig.1B). The onset time for the decline and nadir of the breathing rate was later in the OSA group (P < 0.05). Baseline breathing rates for controls were 14 ± 5 (SD) breaths/min and for OSA were 16 ± 5 breaths/min; baseline heart rates for controls were 64 ± 12 beats/min and for OSA were 73 ± 8 beats/min. Beat-by-beat blood pressure measures, acquired from a subset of subjects outside the scanner, showed higher baseline blood pressures in OSA cases (mean baseline blood pressure: controls, 87 ± 5 mmHg; OSA, 95 ± 12 mmHg) and increased elevation to the challenge in the last 30 s of the challenge (mean rise: controls, 19.3 ± 4.4%; OSA, 26.9 ± 6.6%;P < 0.05, repeated-measures ANOVA).

Fig. 1.

Fig. 1.Mean (± SD) percent change of median heart rate (A) and median breathing rate (B) for 16 control (shaded trace) and 10 obstructive sleep apnea (OSA; solid trace) subjects, plotted over time. Time 0 represents onset of the cold pressor application. * Significant group heart rate differences (P < 0.05); † significant group variability differences (P < 0.05).


Imaging: normalization for overall signal changes.

Overall perfusion changes induced by the challenge or by scanner drift caused global-intensity trends averaging 0.2% difference in gray matter relative to baseline over the course of the scan, and corrections for these trends were incorporated by intensity, normalizing all images before analysis or extraction of time trends.

Control subject responses.

Eight of fourteen anatomically defined regions of interest showed altered signal responses in control subjects (Fig. 2; on these and following images, the neurological presentation convention is adopted, with the right side of images representing the right side of the brain). These areas included the dentate nucleus of the cerebellum; dorsal, medial, and ventral medulla; insula; lentiform nuclei; head of caudate; and ventral midbrain (head of caudate and medial medulla not shown). The hippocampus showed significant signal decrease. Clusters of voxels that increased in a pattern inverse to that of heart rate for control subjects are shown in Fig. 3 and Table 1. These regions included portions of the cerebellar cortex; thalamus; insula; hippocampus; amygdala and superior temporal cortex; and ventral, frontal, and cingulate cortices.

Fig. 2.

Fig. 2.Mean (± SE) signal intensity changes, plotted over time from relevant brain structures in control (n = 16; green) and OSA (n = 10; red) subjects. Selected areas from which these averaged signal intensities were obtained are shown in red and drawn on sagittal or axial views of anatomic slices of the brain from 1 subject and are shown to the left of each time trend. The onset of the cold pressor challenge is at time 0. Time periods of significant (P < 0.05) differences from baseline are indicated by green and red bars for control and OSA subjects, respectively. Group differences are indicated by the yellow bars. * Signal drop-out limited the number of OSA cases in some areas to 9 for the dorsal and ventral medulla and ventral pons. † Similarly, 15 controls were analyzed for the dorsal midbrain.


Fig. 3.

Fig. 3.Top: regions showing signal increases that correlated significantly (P < 0.01) with inverse heart rate in control subjects evoked by the cold pressor challenge. Regions showing significant increases included the cerebellum (1–2), thalamus (5–7; A, C), insular cortex (3–8; A–F), hippocampus, amygdala, and superior temporal cortex (2–4;C–D), frontal cortex (3–4;G–H), and cingulate regions (8–10;G–H). The regions are overlaid onto the mean of the 16 control subjects' anatomic images. The color bar codes tstatistic. Bottom: slice levels for 1–10 in the axial plane and A–H in the coronal plane are shown in the outline drawings.


Table 1. MNI coordinates of significant clusters for 6 response comparisons

tValueVoxelsXYZtValueVoxelsXYZ
Control increase (Fig.3)Control and OSA increase (Fig. 4)
5.272−368125.347640014
4.7199562−85.3393−40−4−14
4.622942−2165.178−83418
4.42936−2−144.7127−30−72−32
4.441−2826−164.5111−18−60−26
4.187−40−4−124.033−58−44
4.0443016−143.76028−72−30
3.9216−6363.618−546−10
3.7256−4−23.5233016−14
3.737−834163.5344−642
3.622−548−183.312400
3.557−18−52−263.214−4−2036
3.51636−24−143.013−36−2614
3.314−12−643.01658−28−2
3.249−34−22143.012−14−200
3.08−30−72−323.014−12−66−6
3.08−58−622.91418−228
2.8156−4642
Control decreaseControl and OSA decrease
4.1418−68144.41478−7012
3.8128−42−303.3156−68
3.66630−64−203.1922−18−18
3.418−14−48−563.129−14−604
3.2758−6023.0124622−2
3.29−60−52162.91532−58−22
3.129−8−5242.87−58−5214
OSA < control (Fig.5)OSA > control (Fig. 6)
4.41478−70124.714−40−74−30
3.3156−683.363626−8
3.1922−18−183.1106584
3.162−38−142.92236−62−24
3.129−14−6042.8520−76−28
3.0124622−22.76−56−1410
2.91532−58−222.738−50−26
2.93180−122.650−5420
2.87−58−52142.6445412
2.66−8−7014

MNI, Montreal Neurological Institute; OSA, obstructive sleep apnea.

Similar responses in OSA and control subjects: cluster analysis.

Significant increases in signal, modeled to the inverse of heart rate, occurred in both groups in the superior and anterior insula, superior temporal cortex, quadrangular lobe of the cerebellum, and the anterior and left posterior cingulate (Fig. 4, Table1). The OSA group showed significantly higher signals than controls in the cerebellar quadrangular cortex and lower signals in the superior (dorsal) insula late in the challenge.

Fig. 4.

Fig. 4.Regions showing significant (P < 0.05) signal increases that correlated with inverse heart rate in major brain structures. Significance levels (t statistic) are conveyed in the color bar and overlaid onto a mean of a set of anatomic images (16 control, 10 OSA) in sagittal, coronal, and axial views. The signal intensity changes, plotted over time, for each cluster at the center of the cross hairs are shown on the right. Time 0represents onset of challenge. * Significant between-group differences are shown at the relevant time points (P < 0.05).


OSA signals different from control subject responses: anatomically defined areas.

Eleven of fourteen anatomically defined regions showed response pattern differences between the OSA and control subjects; 10 of those patterns in which signal intensities were greater in control subjects over OSA cases are shown in Fig. 2. The cerebellar fastigial nucleus (not shown in Fig. 2) was the only one of anatomically defined structures in which OSA signal changes were larger than in control subjects; the remaining three areas (dorsal pons, head of the caudate, medial medulla) showed no group differences. Although the insula signals increased significantly in both groups, the increase in the OSA group was transient, whereas the increase in the control group was sustained. Similarly, the hippocampus in both groups showed decreased signal intensity; however, the change was more sustained in the OSA group.

Cluster analysis: OSA signal values less than control responses.

With modeling to the heart rate, signals in the ventral thalamus, ventral anterior insula, and regions within the hippocampus modestly increased in controls, but decreased in OSA cases (Fig. 5; Table1).

Fig. 5.

Fig. 5.Clusters of voxels in brain areas where OSA group values were significantly lower than controls, correlated with inverse heart rate, are overlaid on anatomic images. Averaged trends for such clusters are also shown. Modest increases emerged in control subjects.


Cluster analysis: OSA greater than control subject responses.

The cluster analysis, modeling to inverse heart rate, demonstrated signal intensity increases in OSA cases within the quadrangular lobe of the cerebellar cortex, anterior insula and medial frontal cortex, and anterior and posterior cingulate gyri, with declines or no change in the control subjects (Fig. 6; Table 1). The responses were, on occasion, expressed unilaterally or were greater on one side over the other, e.g., left posterior cingulate, right insula.

Fig. 6.

Fig. 6.Statistically significant voxels, overlaid on anatomic images from various brain areas, together with averaged trends of signal changes in OSA and control subjects during the course of the cold pressor challenge, where signals increased in OSA and remained the same or declined in control subjects. Time 0represents onset of the challenge. Clusters are color-coded for significance according to the color bar for tstatistic.


DISCUSSION

Several generalizations emerged from this study of a trigeminally mediated cold pressor challenge. The classic physiological response of bradycardia and respiratory rate slowing appeared in control subjects, whereas, in OSA cases, onset of respiratory slowing was delayed, the decreased breathing rate was not as sustained and was substantially more variable, and heart rate did not decline as precipitously. Blood pressure rises to the challenge in a subset of OSA cases outside the scanner were higher. The control neural responses to the challenge included increased signals in thalamic, medullary, and cerebellar cortex; dentate nucleus; and insular, cingulate, and frontal cortex areas. Cerebellar, insular, cingulate, and frontal cortex areas have previously been implicated in cardiovascular activation triggered by cognitive tasks and isometric exercises in normal subjects by using positron emission tomography procedures (5a). The principal response pattern differences between groups emerged in the cerebellar cortex, the insular and posterior cingulate cortex, thalamus, and the hippocampus. The response patterns differed inversely in certain areas. Cerebellar and posterior cingulate structures responded to the challenge with decreased signal in controls, but increased signal in OSA; signals in insular, thalamic, and hippocampal areas decreased in OSA cases, but remained largely unchanged in controls. Some of these regions have earlier been shown to exhibit diminished gray matter volume in OSA (23) or in heart failure patients with sleep-disordered breathing (51). The findings suggest multiple impaired neural sites in OSA patients, with cerebellar and limbic structures particularly implicated as deficient in the pressor challenge.

Cold pressor challenge.

The cardiovascular and breathing effects elicited by cold facial stimulation (the dive reflex) originate from primitive processes that, under unusual circumstances, serve to override homeostatic control mechanisms. The processes include induction and maintenance of apnea necessary for underwater motility in lower animals, despite increasing CO2 and reduced O2 availability, and increases in sympathetic outflow, despite an expected decrease after blood pressure elevation accompanying dives (27). The sudden drop in forehead temperature activates trigeminal sensory afferents, which synapse in the spinal trigeminal nucleus of the medulla (19, 27). Neurons within the spinal trigeminal nucleus project to a number of brain stem cardiorespiratory regions, including the nucleus of the solitary tract, parabrachial regions of the dorsal pons, and, indirectly, to A1, A5, and C1 cells of the medulla (7,28, 37). Neurons within the spinal trigeminal nucleus project directly to the medial ventroposterior region of the thalamus, which in turn projects to the lateral aspect of the precentral gyrus for perception of temperature. The challenge elicited less of a rate slowing for both breathing and heart rate and a delayed onset for the respiratory response in OSA. An enhanced blood pressure response over controls was found in a subset of OSA cases, suggesting that reactions to blood pressure in the cold pressor challenge were also altered. However, those values were obtained subsequent to the imaging task (necessary because a MR-compatible noninvasive blood pressure monitor was not available), did not include all of the subjects, and were obtained in an environment similar, but not identical, to the scanning circumstance. Thus correlations of signal changes were taken with heart rate measures.

Thalamic structures.

The ventroposterior thalamic region showed increased signal intensities in controls, but not in OSA subjects, when signal changes were modeled to heart rate. Thus the differences in responses between the two groups may depend partially on a loss of primary sensory input to thalamic areas, a possibility that is supported by findings of diminished responses to oral airway stimulation in OSA subjects (24).

Medullary and midbrain structures.

The absence, in OSA cases, of normal responses within the dorsal medulla may also reflect a diminished recruitment of the nucleus of the solitary tract by reduced trigeminal nucleus input or diminished participation of the solitary tract nucleus. The ventral medullary areas showed an early decline in the OSA cases; this region receives projections from the solitary tract nucleus and plays a significant role in respiratory pattern generation. The early decline in signal in the ventral medulla may contribute to the late onset of respiratory slowing in the OSA group. Similarly, diminished signal changes, relative to controls, in the dorsal and ventral midbrain likely represent reduced participation of periaqueductal gray regions. These areas project to the nucleus of the solitary tract and to the ventral medulla, modify respiratory patterning and cardiovascular control (13, 17), and have the potential to influence the pattern differences found here.

Cerebellar structures.

Localized regions within the quadrangular lobe of the cerebellar cortex showed a signal intensity rise to the challenge in OSA subjects, as opposed to little change in signal intensity in the control group. Other localized regions within the cerebellar cortex showed comparable responses in both groups. The cerebellar cortical areas of aberrant responses overlapped regions of gray matter loss in OSA patients (23), as well as patients with heart failure showing disordered breathing during sleep (51). Cerebellar structures classically have been associated with roles in coordinating somatic afferent with efferent information in the performance of “error correction” on a moment-to-moment basis. A role in breathing control is frequently overlooked, despite recognition that respiratory muscles represent a substantial proportion of all skeletal musculature and that respiratory muscle modulation is heavily dependent on afferent information from thoracic wall stretch receptors, chemoreceptors, and arterial pressure receptors. Experimental evidence from animal models emphasizes the role of the deep cerebellar nuclei in regulating phrenic outflow and regulation of inspiratory onset after apnea (22). Moreover, neurons in the fastigial nucleus respond to CO2 and mechanical manipulation (54), and interpositus nuclei are involved in cough responses, which are associated with transient elevations of blood pressure (54). Damage to cerebellar areas is accompanied by a high incidence of sleep-disordered breathing, including OSA, both in human adults and infants (4, 50).

The breathing response in OSA subjects during the cold pressure challenge was characterized by extreme variability, with some subjects exhibiting increased, and others decreased, breathing rates. These varied responses appear to reflect a “loss of modulatory control” by relevant neural structures to the cold pressure challenge. The increase in signal intensity in the cerebellar cortex of OSA subjects may reflect efforts accounting for this loss of normal modulatory control.

Insula.

The insula signal intensity changes differed significantly in OSA cases compared with controls. The insular cortex plays substantial roles in autonomic control in humans (34). The pattern of responses in the insular cortex paralleled the heart rate responses in both groups. In OSA subjects, both heart rate and insula signal intensity returned to baseline levels significantly earlier than in controls. Discrete portions of the insula also markedly differed in responses to the challenge in OSA over controls, with significantly greater declines in signal in OSA compared with controls in the right ventral insula. The enhanced laterality of the response difference to the right side is of particular interest, because the nature of autonomic influence exerted unilaterally by the insula differs, with parasympathetic outflow modified by the left insula and sympathetic outflow modulated by the right (33, 35, 36). The larger neural signal response difference on the right (sympathetic-related) side may underlie the higher basal heart rates found in the OSA group. Regions within the insula show substantial gray matter loss in heart failure patients with sleep-disordered breathing, a group that also shows high sympathetic tone (51); a similar gray matter decline was not detected in OSA patients (23), possibly because of limited resolution of the anatomic technique. It remains to be determined whether enhanced-resolution MR procedures will reveal gray matter loss in OSA; however, the extent of alterations in neural response indicates significant functional impairments in a brain structure shown by others to contribute substantially to autonomic control.

Hippocampus.

The tissue loss found earlier in OSA also included regions of the hippocampus that are damaged by intermittent hypoxia in animals (12) and an area preferentially susceptible to injury as a consequence of seizure activity triggered by excitation of the principal afferent pathway, the perforant pathway, or other seizure-related hypoxia-ischemia (46). In addition to cognitive roles apparently affected in both OSA and epilepsy (2, 55), hippocampal structures participate in cardiovascular and respiratory regulation. Direct chemical stimulation of the hippocampus evokes marked decreases in arterial pressure and heart rate (44). Furthermore, hippocampal structures contain neurons that discharge with the respiratory and cardiac cycle (10), show significant patterns of activity change before sigh-apnea sequences in animals (42), and, on stimulation, alter breathing patterns (44). The hippocampus may play a modulatory role, adjusting the gain of the cardiorespiratory reflex in the cold pressor challenge, because resting arterial pressure and heart rate are unaffected by hippocampal blockade (49). The diminished hippocampal response in the OSA group may represent a reduced modulatory role, resulting in the significant differences in heart rate and respiratory responses in OSA subjects in response to the cold application.

Cingulate gyrus.

The posterior cingulate also showed significant response differences from those of controls. A motor role for both anterior and posterior portions of the cingulate has been recognized by several stimulation studies (e.g., Ref. 6). The anterior cingulate contains neurons that discharge with both the respiratory and cardiovascular cycles (9) and is involved in vocalization, a task requiring intricate control of the fine musculature of the upper airway (48). The role in vocalization has been especially related to affective components of expressive sounds. Dyspnea activates the right posterior cingulate (40), and this activation may derive from affective contributions to the breathing task. However, the posterior cingulate functions may be even more complex. A growing body of evidence suggests that the structure integrates information of internal movement cues for spatial navigation (16). It is possible that the “internal movement” cues extend to appropriate motor behaviors of a wider range, including respiratory muscle motions in these OSA subjects.

More rostral portions of the cingulate gyrus show gray matter loss in OSA subjects (23). The rat posterior cingulate shows a region-specific sensitivity to excitotoxicity apparently mediated through N-methyl-d-aspartate-related disinhibition (20), a finding that suggests that the area may be at particular risk in humans during the intermittent hypoxia of OSA. Cingulate areas near sites showing significantly different response patterns from controls project to the anterior insula and parahippocampal regions in primates (29). Thus deficient functional and anatomic characteristics appear in regions with anatomic projections to other areas showing gray matter loss. These areas include the insula of heart failure patients with sleep-disordered breathing and the cingulate and hippocampus of OSA (23) and heart failure cases (51). It is unclear whether the commonality of axonal connections and response deficiencies relate to aspects of hypoxic or other damage in sleep-related breathing pathology.

Portions of the cingulate gyrus play a role in mediating responses to blood pressure challenges, as shown by functional MRI evidence from animal and human studies (13, 18, 56). Stimulation results in large changes in blood pressure and heart rate (3). It may be the case that the posterior cingulate is involved in the heart rate slowing found here; however, the history of its role in motor control suggests that participation in the respiratory pattern changes is more likely.

Lateralization of responses.

The cold pressor challenge applied to the entire forehead should elicit bilateral recruitment of primary sensory neurons and, therefore, activity of bilateral central neural structures. The accompanying blood pressure elevation and slowing of heart and respiratory rates would presumably be mediated by areas on both sides of the brain. Instead, pronounced lateralization of signal changes was observed in multiple regions, including the right insula, right cingulate, and cerebellar cortex. The laterality finding was similar to that observed previously with multiple blood pressure elevation stimuli (13) and also noted in response to pronounced inspiratory loading challenges (40). As in the earlier studies, no relationship between lateralization of the signal and handedness of the subjects was observed.

The lateralization of neural responses in some structures in response to autonomic challenges, a finding also shown elsewhere (5a, 13), raises significant theoretical and practical issues. The finding places aspects of breathing and cardiovascular control within the framework of a variety of other behavioral acts, including language expression and voluntary motor behaviors, both of which are highly lateralized. Unilateral neural damage from OSA may interact with ipsilateral breathing control mechanisms and may be especially important during sleep, in which affected areas may be less capable of an adequate response, lacking the parallel redundancy of respiratory control typically available in waking states.

Overall perfusion changes.

Any manipulation that exerts significant blood pressure effects, including cold pressor challenges, can alter overall neural signal properties through gross perfusion changes. Autoregulatory processes normally operate to protect against major or sustained perfusion changes, including those from hypotension (57). However, autoregulatory processes may not operate normally in the OSA population, by reason of repetitive hypoxic, hypercarbic, and transient hypertensive damage associated with obstructed breathing and momentary global perfusion changes that can occur with autonomic manipulations (14). Overall perfusion changes induced by the dive reflex are of particular concern, because the primitive reflex response partially serves to conserve metabolic resources for underwater survival. Manipulations such as increased inspired CO2 can elicit overall signal changes from altered cerebral perfusion. However, determination of regional neural changes in primary visual cortex to photic stimulation can be readily evaluated in such conditions, because the changes are superimposed on the hypercapnia-induced global signal intensity alterations (5, 21), and the global changes can readily be partitioned by intensity normalization procedures. We found an average 0.2% global change in gray matter in both groups and used intensity normalization of gray matter signals to control for the contribution of global changes from regional patterns. Separation of global trends in gray matter from those in fibers precluded differential global effect contributions in the normalization process. Comparable procedures minimized effects of other autonomic manipulations in normal and OSA subjects (14). The use of normalization procedures, together with findings that areas immediately adjacent to recruited sites showed different patterns of signal intensity response, argue that signal changes reported in this study result from differential recruitment of brain sites and not from overall changes in cerebral blood flow.

In summary, cerebellar, insular, and cingulate cortex areas, as well as hippocampal, medullary, midbrain, and thalamic sites, showed significant response pattern differences to the cold pressor challenge in OSA cases. Relative to control subjects, cerebellar and limbic areas appeared to be particularly affected when signal changes modeled to heart rate were considered. The affected limbic structures have been associated with fine control of the upper airway musculature and autonomic outflow. The cerebellar dysfunction may have contributed to a mismatch of afferent input with coordinated motor output. The inappropriate neural responses to the cold pressor challenge suggest that specific neural structures involved in affect and motor control contribute to the respiratory dysfunctions of OSA.

This research was supported by National Heart, Lung, and Blood Institute Grant HL-60296.

FOOTNOTES

  • Address for reprint requests and other correspondence: R. M. Harper, Dept. of Neurobiology, Univ. of California at Los Angeles, Los Angeles, CA 90095-1763 (E-mail:).

  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

  • First published January 3, 2003;10.1152/japplphysiol.00881.2002

REFERENCES

  • 1 [Anonymous].Practice parameters for the indications for polysomnography and related procedures. Polysomnography Task Force, American Sleep Disorders Association Standards of Practice Committee.Sleep201997406422
    Crossref | ISI | Google Scholar
  • 2 Beebe DW, Gozal D.Obstructive sleep apnea and the prefrontal cortex: towards a comprehensive model linking nocturnal upper airway obstruction to daytime cognitive and behavioral deficits.J Sleep Res112002116
    Crossref | PubMed | ISI | Google Scholar
  • 3 Burns SM, Wyss JM.The involvement of the anterior cingulate cortex in blood pressure control.Brain Res34019857177
    Crossref | PubMed | ISI | Google Scholar
  • 4 Chokroverty S, Sachdeo R, Masdeu J.Autonomic dysfunction and sleep apnea in olivopontocerebellar degeneration.Arch Neurol411984926931
    Crossref | Google Scholar
  • 5 Corfield DR, Murphy K, Josephs O, Adams L, Turner R.Does hypercapnia-induced cerebral vasodilation modulate the hemodynamic response to neural activation?Neuroimage13200112071211
    Crossref | PubMed | ISI | Google Scholar
  • 5a Critchley HD, Corfield DR, Chandler MP, Mathias CJ, Dolan RJ.Cerebral correlates of autonomic cardiovascular arousal: a functional neuroimaging investigation in humans.J Physiol5232000259270
    Crossref | PubMed | ISI | Google Scholar
  • 6 Diehl B, Dinner DS, Mohamed A, Najm I, Klem G, LaPresto E, Bingaman W, Luders HO.Evidence of cingulate motor representation in humans.Neurology552000725728
    Crossref | ISI | Google Scholar
  • 7 Esser MJ, Pronych SP, Allen GV.Trigeminal-reticular connections: possible pathways for nociception-induced cardiovascular reflex responses in the rat.J Comp Neurol3911998526544
    Crossref | PubMed | ISI | Google Scholar
  • 8 Friston KJ, Holmes AP, Worsley K, Poline JB, Frith CD, Frackowiak RS.Statistic parametric maps in functional imaging: a general linear approach.Hum Brain Mapp21995189210
    Crossref | Google Scholar
  • 9 Frysinger RC, Harper RM.Cardiac and respiratory relationships with neural discharge in the anterior cingulate cortex during sleep-waking states.Exp Neurol941986247263
    Crossref | ISI | Google Scholar
  • 10 Frysinger RC, Harper RM.Cardiac and respiratory correlations with unit discharge in human amygdala and hippocampus.Electroencephalogr Clin Neurophysiol721989463470
    Crossref | PubMed | Google Scholar
  • 11 Garrigue S, Bordier P, Jais P, Shah DC, Hocini M, Raherison C, Tunon De Lara M, Haissaguerre M, Clementy J.Benefit of atrial pacing in sleep apnea syndrome.N Engl J Med3462002404412
    Crossref | ISI | Google Scholar
  • 12 Gozal D, Daniel JM, Dohanich GP.Behavioral and anatomical correlates of chronic episodic hypoxia during sleep in the rat.J Neurosci21200124422450
    Crossref | PubMed | ISI | Google Scholar
  • 13 Harper RM, Bandler R, Spriggs D, Alger JR.Lateralized and widespread brain activation during transient blood pressure elevation revealed by magnetic resonance imaging.J Comp Neurol4172000195204
    Crossref | PubMed | ISI | Google Scholar
  • 14 Henderson LA, Woo MA, Macey PM, Macey KE, Frysinger RC, Alger JR, Yan-Go F, Harper RM.Neural responses during Valsalva maneuvers in obstructive sleep apnea syndrome.J Appl Physiol94200310631074
    Link | ISI | Google Scholar
  • 15 Heymans C, Neil E.Reflexogenic Areas of the Cardiovascular System.1958J & A ChurchillLondon
    Google Scholar
  • 16 Katayama K, Takahashi N, Ogawara K, Hattori T.Pure topographical disorientation due to right posterior cingulate lesion.Cortex351999279282
    Crossref | ISI | Google Scholar
  • 17 Keay KA, Crowfoot LJ, Floyd NS, Henderson LA, Christie MJ, Bandler R.Cardiovascular effects of microinjections of opioid agonists into the “Depressor Region” of the ventrolateral periaqueductal gray region.Brain Res76219976171
    Crossref | PubMed | ISI | Google Scholar
  • 18 King AB, Menon RS, Hachinski V, Cechetto DF.Human forebrain activation by visceral stimuli.J Comp Neurol4131999572582
    Crossref | PubMed | ISI | Google Scholar
  • 19 Kumada M, Dampney RA, Reis D.The trigeminal depressor response: a cardiovascular reflex originating from the trigeminal system.Brain Res921975485489
    Crossref | PubMed | ISI | Google Scholar
  • 20 Li Q, Clark S, Lewis DV, Wilson WA.NMDA receptor antagonists disinhibit rat posterior cingulate and retrosplenial cortices: a potential mechanism of neurotoxicity.J Neurosci22200230703080
    Crossref | PubMed | ISI | Google Scholar
  • 21 Li TQ, Kastrup A, Moseley ME, Glover GH.Changes in baseline cerebral blood flow in humans do not influence regional cerebral blood flow response to photic stimulation.J Magn Reson Imaging122000757762
    Crossref | PubMed | ISI | Google Scholar
  • 22 Lutherer LO, Williams JL.Stimulating fastigial nucleus pressor region elicits patterned respiratory responses.Am J Physiol Regul Integr Comp Physiol2501986R418R426
    Link | ISI | Google Scholar
  • 23 Macey PM, Henderson LA, Macey KE, Alger JR, Frysinger RC, Woo MA, Harper RK, Yan-Go FL, Harper RM.Brain morphology associated with obstructive sleep apnea.Am J Respir Crit Care Med166200213821387
    Crossref | PubMed | ISI | Google Scholar
  • 24 Malhotra A, Fogel RB, Edwards JK, Shea SA, White DP.Local mechanisms drive genioglossus activation in obstructive sleep apnea.Am J Respir Crit Care Med161200017461749
    Crossref | PubMed | ISI | Google Scholar
  • 25 Marks JD, Harper RM.Differential inhibition of the diaphragm and posterior cricoarytenoid muscles induced by transient hypertension across sleep states in intact cats.Exp Neurol951987730742
    Crossref | ISI | Google Scholar
  • 26 Mazziotta JC, Toga AW, Evans A, Fox P, Lancaster J.A probabilistic atlas of the human brain: theory and rationale for its development. The International Consortium for Brain Mapping (ICBM).Neuroimage2199589101
    Crossref | PubMed | ISI | Google Scholar
  • 27 McCulloch PF, Panneton WM, Guyenet PG.The rostral ventrolateral medulla mediates the sympathoactivation produced by chemical stimulation of the rat nasal mucosa.J Physiol5161999471484
    Crossref | PubMed | ISI | Google Scholar
  • 28 Menetrey D, Basbaum AI.Spinal and trigeminal projections to the nucleus of the solitary tract: a possible substrate for somatovisceral and viscerovisceral reflex activation.J Comp Neurol2551987439450
    Crossref | PubMed | ISI | Google Scholar
  • 29 Morecraft RJ, Van Hoesen GW.Convergence of limbic input to the cingulate motor cortex in the rhesus monkey.Brain Res Bull451998209232
    Crossref | PubMed | ISI | Google Scholar
  • 30 Narkiewicz K, Somers VK.Cardiovascular variability characteristics in obstructive sleep apnea.Auton Neurosci9020018994
    Crossref | ISI | Google Scholar
  • 31 Ogawa S, Lee TM, Nayak AS, Glynn P.Oxygenation-sensitive contrast in magnetic resonance image of rodent brain at high magnetic fields.Magn Reson Med1419906878
    Crossref | PubMed | ISI | Google Scholar
  • 32 Olsen KD, Kern EB.Nasal influences on snoring and obstructive sleep apnea.Mayo Clin Proc65199010951105
    Crossref | ISI | Google Scholar
  • 33 Oppenheimer SM, Cechetto DF.Cardiac chronotropic organization of the rat insular cortex.Brain Res53319906672
    Crossref | PubMed | ISI | Google Scholar
  • 34 Oppenheimer SM, Cechetto DF, Hachinski VC.Cerebrogenic cardiac arrhythmias. Cerebral electrocardiographic influences and their role in sudden death.Arch Neurol471990513519
    Crossref | PubMed | Google Scholar
  • 35 Oppenheimer SM, Gelb A, Girvin JP, Hachinski VC.Cardiovascular effects of human insular cortex stimulation.Neurology42199217271732
    Crossref | PubMed | ISI | Google Scholar
  • 36 Oppenheimer SM, Kedem G, Martin WM.Left-insular cortex lesions perturb cardiac autonomic tone in humans.Clin Auton Res61996131140
    Crossref | PubMed | ISI | Google Scholar
  • 37 Panneton WM, Johnson SN, Christensen ND.Trigeminal projections to the peribrachial region in the muskrat.Neuroscience581994605625
    Crossref | PubMed | ISI | Google Scholar
  • 38 Parker JM, Alger JR, Woo MA, Spriggs D, Harper RM.Acquisition of electrophysiologic signals during magnetic resonance imaging.Sleep22199911251126
    Crossref | PubMed | ISI | Google Scholar
  • 39 Partinen M, Hublin C.Epidemiology of sleep disorders.Principles and Practice of Sleep Medicine3rd ed., Kryger MH, Roth T, Dement WC.2000558579SaundersPhiladelphia, PA
    Google Scholar
  • 40 Peiffer C, Poline JB, Thivard Aubier ML, Samson Y.Neural substrates for the perception of acutely induced dyspnea.Am J Respir Crit Care Med1632001951957
    Crossref | PubMed | ISI | Google Scholar
  • 41 Pendlebury ST, Pepin JL, Veale D, Levy P.Natural evolution of moderate sleep apnoea syndrome: significant progression over a mean of 17 months.Thorax521997872878
    Crossref | ISI | Google Scholar
  • 42 Poe GR, Kristensen MP, Rector DM, Harper RM.Hippocampal activity during transient respiratory events in the freely behaving cat.Neuroscience7219963948
    Crossref | PubMed | ISI | Google Scholar
  • 44 Ruit KG, Neafsey EJ.Cardiovascular and respiratory responses to electrical and chemical stimulation of the hippocampus in anesthetized and awake rats.Brain Res4571988310321
    Crossref | PubMed | ISI | Google Scholar
  • 45 Sachs C, Kaijser L.Autonomic regulation of cardiopulmonary functions in sleep apnea syndrome and narcolepsy.Sleep51982227238
    Crossref | PubMed | ISI | Google Scholar
  • 46 Schmidt-Kastner R, Freund TF.Selective vulnerability of the hippocampus in brain ischemia.Neuroscience401991599636
    Crossref | PubMed | ISI | Google Scholar
  • 47 Veale D, Pépin JL, Lévy PA.Autonomic stress tests in obstructive sleep apnea syndrome and snoring.Sleep151992505513
    ISI | Google Scholar
  • 48 Von Cramon D, Jurgens U.The anterior cingulate cortex and the phonatory control in monkey and man.Neurosci Biobehav Rev71983433425
    Crossref | ISI | Google Scholar
  • 49 Wang JQ, Ingenito AJ.Centrally mediated cardiovascular actions of dynorphin A(1–8) on rat hippocampal formation.J Pharmacol Exp Ther2611992678685
    ISI | Google Scholar
  • 50 Waters KA, Forbes P, Morielli A, Hum C, O'Gorman AM, Vernet O, Davis GM, Tewfik TL, Ducharme FM, Brouillette RT.Sleep-disordered breathing in children with myelomeningocele.J Pediatr1321998672681
    Crossref | ISI | Google Scholar
  • 51 Woo MA, Macey PM, Fonarow GC, Hamilton MA, Harper RM.Regional brain cell loss in heart failure.Circulation1062002369
    ISI | Google Scholar
  • 52 Woods RP, Mazziotta JC, Cherry SR.MRI-PET registration with automated algorithm.J Comput Assist Tomogr171993536546
    Crossref | PubMed | ISI | Google Scholar
  • 53 Xu F, Frazier DT.Respiratory-related neurons of the fastigial nucleus in response to chemical and mechanical challenges.J Appl Physiol82199711771184
    Link | ISI | Google Scholar
  • 54 Xu F, Frazier DT, Zhang Z, Baekey DM, Shannon R.Cerebellar modulation of cough motor pattern in cats.J Appl Physiol831997391397
    Link | ISI | Google Scholar
  • 55 Ylinen A, Lahtinen H, Sirvio J, Partanen J, Asikainen A, Gulyas A, Freund TF, Riekkinen P.Behavioural, electrophysiological and histopathological changes following sustained stimulation of the perforant pathway input to the hippocampus: effect of the NMDA receptor antagonist, CGP 39551.Brain Res5531991195200
    Crossref | ISI | Google Scholar
  • 56 Yu PL, Henderson LA, Frysinger RC, Galons JP, Bandler R, and Harper RM. Functional magnetic resonance imaging (fMRI) of feline brain during depressor challenges (Abstract). Soc Neurosc AbstrA690.8, 2001.
    Google Scholar
  • 57 Zaharchuk G, Mandeville JB, Bogdanov AA, Weissleder R, Rosen BR, Marota JJ.Cerebrovascular dynamics of autoregulation and hypoperfusion. An MRI study of CBF and changes in total and microvascular cerebral blood volume during hemorrhagic hypotension.Stroke30199921972205
    Crossref | PubMed | ISI | Google Scholar