PLoS ONE
Home Insular functional organization during handgrip in females and males with obstructive sleep apnea
Insular functional organization during handgrip in females and males with obstructive sleep apnea
Insular functional organization during handgrip in females and males with obstructive sleep apnea

Competing Interests: The authors have declared that no competing interests exist.

Article Type: research-article Article History
Abstract

Study objectives

Brain regulation of autonomic function in obstructive sleep apnea (OSA) is disrupted in a sex-specific manner, including in the insula, which may contribute to several comorbidities. The insular gyri have anatomically distinct functions with respect to autonomic nervous system regulation; yet, OSA exerts little effect on the organization of insular gyral responses to sympathetic components of an autonomic challenge, the Valsalva. We further assessed neural responses of insular gyri in people with OSA to a static handgrip task, which principally involves parasympathetic withdrawal.

Methods

We measured insular function with blood oxygen level dependent functional MRI. We studied 48 newly-diagnosed OSA (age mean±std:46.5±9 years; AHI±std:32.6±21.1 events/hour; 36 male) and 63 healthy (47.2±8.8 years;40 male) participants. Subjects performed four 16s handgrips (1 min intervals, 80% subjective maximum strength) during scanning. fMRI time trends from five insular gyri—anterior short (ASG); mid short (MSG); posterior short (PSG); anterior long (ALG); and posterior long (PLG)—were assessed for within-group responses and between-group differences with repeated measures ANOVA (p<0.05) in combined and separate female-male models; age and resting heart-rate (HR) influences were also assessed.

Results

Females showed greater right anterior dominance at the ASG, but no differences emerged between OSA and controls in relation to functional organization of the insula in response to handgrip. Males showed greater left anterior dominance at the ASG, but there were also no differences between OSA and controls. The males showed a group difference between OSA and controls only in the ALG. OSA males had lower left activation at the ALG compared to control males. Responses were mostly influenced by HR and age; however, age did not impact the response for right anterior dominance in females.

Conclusions

Insular gyri functional responses to handgrip differ in OSA vs controls in a sex-based manner, but only in laterality of one gyrus, suggesting anterior and right-side insular dominance during sympathetic activation but parasympathetic withdrawal is largely intact, despite morphologic injury to the overall structure.

Pal,Ogren,Aysola,Kumar,Henderson,Harper,Macey,and King: Insular functional organization during handgrip in females and males with obstructive sleep apnea

Introduction

Cardiovascular disease (CVD) in obstructive sleep apnea (OSA) is difficult to treat [1, 2], possibly due to long-term changes in the autonomic nervous system (ANS) that underpin high sympathetic tone and disrupt blood pressure regulation [3]. Associations between OSA and hypertension, coronary artery disease, stroke, arrhythmias and death are well established [4]. However, most studies exploring the comorbidities associated with OSA include both male and female participants, and evidence of whether these links are affected by sex remain under investigation [5]. A recent study of almost 2 million OSA subjects and 2 million matched controls reported that hypertension was more prevalent in women than men in the condition [6]. Furthermore, there is evidence of a female predominance of OSA-related stoke as well as more severe detrimental changes in endothelial function, peak blood flow, systemic inflammation, and digital vascular function [7, 8]. Other clinical cardiovascular characteristics, such as morning BP patterns and responses to acute BP challenges in OSA also vary by sex [3, 9]. These differences in physiologic responses between sexes in OSA presumably have an underlying neural basis.

Our earlier work has shown some sex-specific effects of OSA on neural regulation of cardiovascular stimuli. We found central autonomic responses differed in male OSA versus healthy groups [1013]. Furthermore, although autonomic maneuvers such as the cold pressor test, hand grip, and Valsalva maneuver evoked lower amplitude, delayed onset, and slower heart rate changes in combined male and female OSA patients over healthy people [3, 14, 15], these alterations were more pronounced in females with the sleep disorder. These autonomic dysfunctions may reflect sex-related central injury. There is also evidence that hormones play a significant role in the development of OSA itself. For example, the prevalence of OSA increases dramatically following menopause in females [16], and polycystic ovary syndrome which is characterized by menstrual disturbances, excess androgen, and often obesity, is associated with increased prevalence of OSA [17]. A considerable body of literature links alterations in various hormones, including testosterone, with OSA [18], although mechanisms outlining sex-based processes in that regulation have yet to be established.

The hypothalamus is a major regulator of hormone release, and this structure in turn closely interacts with the anterior insula [19, 20]. The insula also has substantial projections to brainstem output nuclei [21]; thus, the structure is a major contributor to autonomic regulation [22, 23]. Multiple studies revealed that the insula as a whole shows significant OSA-related fMRI signal changes [15, 2427], structural injury or adaptations [2833], reduced perfusion [34], and altered metabolic state [3537]. When considered across large scale networks, the resting connectivity of the insula consistently shows OSA-related differences [32, 33, 38], suggesting either an insula-specific deficiency, or network-wide alterations. While precise regional variations in insular changes in OSA are currently being explored, the insula contains various regions with differing developmental beginnings and projection patterns. For example, the posterior insula is granular/dysgranular in nature and projects to other high-order cortical regions, whereas the anterior insular is predominantly agranular in nature and projects to lower brain regions such as the amygdala, ventral striatum, and autonomic regulatory regions, including the hypothalamus and various brainstem nuclei [19, 39]. These projection and morphologic patterns leave a potential for sex to exert differential influences on autonomic patterns in different regions.

To probe autonomic regulation, we used functional MRI (fMRI) measures of neural function during tasks that elicit an autonomic response [40]. Three standard autonomic challenges often used in combination are the cold pressor, handgrip and Valsalva maneuver; other standard tests such as tilt, Mueller maneuver, or pharmacological manipulations are less amenable to fMRI scanning [41]. The cold pressor, handgrip and Valsalva are all pressor challenges, that is they raise blood pressure, but through different mechanisms. Thus, common fMRI patterns across these challenges could be interpreted as reflecting blood pressure regulatory processes, whereas challenge-specific responses could reflect task-related effects. Both cold pressor and Valsalva challenges predominantly increase sympathetic activity during the active phase of the challenge; whereas, the handgrip adds an element of parasympathetic withdrawal in the short term (10 sec) [4245].

In healthy adults, we showed the functional organization of the insular cortex is gyri-specific for the handgrip, a challenge that combines autonomic perturbation and an intentional motor activity [46, 47]. Specifically, the anterior insula is more activated during the early, predominantly parasympathetic withdrawal phase [48, 49]. In contrast, the strain phase of the handgrip challenge, which is associated with a moderate HR and sympathetic increase, elicited the greatest responses in the middle insular gyri. These sex-specific neural patterns co-occur with sex-specific peripheral differences, with females displaying HR smaller increases [3, 47]. The data are consistent with a recent human study showing tachycardic responses elicited by stimulating the posterior insula and bradycardic responses from stimulating the more anterior insular cortex [50]. The insula also displays a laterality effect, with the right side more closely aligned with sympathetic, and the left side parasympathetic, activity changes [48, 51, 52]. Thus, we expect insular responses to OSA during handgrip to differ by sex, hemisphere and gyral subregions [46, 53].

The objective here was to determine the nature of insular functional organization during a handgrip challenge in OSA, both controlling for sex statistically, and considering females and males separately. This study was a secondary analysis of a dataset we collected earlier. Given the anterior autonomic and left-sided parasympathetic role of the insula, and the reduced cardiovascular responses to a handgrip in OSA, we hypothesized an anterior dominance of fMRI responses in that condition, and since the left insula serves more- parasympathetic aspects, we hypothesized a greater left-side response in affected subjects. Since cardiovascular responses differ by sex in healthy people, we further hypothesized that alterations in insular organization contributed to those sex differences in people with and without OSA.

Methods

Participants

We studied 111 adults consisting of 48 newly-diagnosed, untreated OSA patients (36 males, 12 females) and 63 healthy control participants (43 males, 23 females); details are in Table 1. Inclusion criteria for all participants included age 21–65 and weight <125kg (MRI constraint). Participants in the control group were in good health and those in the OSA group had a diagnosis of OSA. Recruitment was principally via fliers posted at the UCLA Sleep Disorders sleep clinic and on the campus and nearby communities, with additional fliers emailed or given in person upon request based on word-of-mouth. OSA patients were recruited from the UCLA Sleep Disorders Center from Dec 2005 to Aug 2008, and were recently (< 2 months) diagnosed with OSA according to the 1999 American Academy of Sleep Medicine guidelines based on an in-clinic full polysomnography test [54]. After completing a phone screening (see S1 File), participants were invited to UCLA. At the visit prior to consenting, control participants were screened for OSA using a semi-structured interview to assess daytime sleepiness, snoring, bed partner report of breathing difficulties during sleep, and nighttime gasping episodes, and referred to a full sleep study if those symptoms were present (see S1 File). Exclusion criteria for all participants included other sleep disorders, major illness or head injury, stroke, major cardiovascular disease, diabetes, current use of psychotropic or cardiovascular medications other than statins, and diagnosed mental disorder. Exclusion criteria also included MRI contraindications, including metallic implants not classified as safe at 3 Tesla, pregnancy and claustrophobia. Within the larger OSA population the present OSA sample could be considered representative of the minority of people who are relatively healthy and receiving standard healthcare (no major comorbidities, no mental health diagnoses or medications, no autonomically-active anti-hypertensive medications, limited obesity, UCLA patients).

Table 1
Participant information.
AllMaleFemale
CONTROL Mean ± std [Range] N = 63OSA Mean ± std [Range] N = 48p OSA vs. CONTROLCONTROL Mean ± std [Range] N = 40OSA Mean ± std [Range] N = 36p OSA vs. CONTROLCONTROL Mean ± std [Range] N = 23OSA Mean ± std [Range] N = 12p OSA vs. CONTROL
Age (years)47.5 ± 8.8 [30.9–65.8]46.5 ± 9.0 [30.8–62.7]0.5745.9 ± 9.1 [30.9–64.5]44.9 ± 8.9 [30.8–62.7]0.6550.3 ± 7.8 [40.2–65.8]51.3 ± 7.9 [37.0–62.2]0.73
BMI (m²/kg)24.7 ± 3.7 [16.6–35.5]30.5 ± 5.1 [21.3–43.2]<0.00125.2 ± 2.8 [17.6–29.8]29.6 ± 4.7 [21.4–43.2]<0.00123.94 ± 5.0 [16.6–35.5]32.9 ± 5.7 [21.3–41.4]<0.001
Resting HR (bpm)68.9 ± 11.6 [45.8–102.1]71.7 ± 9.9 [51.6–93.5]0.1966.7 ± 11.1 [45.8–95.7]72.8 ± 10.4 [53.1–93.5]0.01672.8 ± 11.5 [57.7–102.1]68.2 ± 7.9 [51.6–82.3]0.23
Sleep parameters for OSA
AHI (events/hour)n/a32.6 ± 21.1 [5.0–100.7]n/an/a34.6 ± 19.7 [10.0–100.7]n/an/a26.7 ± 24.8 [5.0–89.4]n/a
SaO2 (minimum%)n/a80.3 ± 9.4 [50.0–96.0]n/an/a78.2 ± 9.4 [50.0–96.0]n/an/a86.2 ± 6.4 [73.0–96.0]n/a
SaO2 (mean%)n/a94.9 ± 1.9 [88.0–97.0]n/an/a94.9 ± 2.0 [88.0–97.0]n/an/a94.8 ± 1.5 [92.0–97.0]n/a

Characteristics of OSA and control groups, with separation by sex. Group differences were tested with two-way ANOVA for OSA parameters, p values have been indicated (italicized if ≤ .05). HR was recorded in the scanner over 1 minute immediately prior to the first handgrip task. Sleep parameters were based on the patients’ polysomnographic study.

p for two way ANOVA F-test, group comparison OSA vs. CONTROL

The procedures were approved by the UCLA Institutional Review Board. All participants were provided a description of the procedures priori to visiting UCLA, and upon their visit those procedures were reviewed and participants provided written, informed consent.

Measurements

Brain blood-oxygen level dependent (BOLD) fMRI signals were recorded in a 3.0 Tesla MRI scanner (Siemens, Magnetom, Trio) with an 8-channel head coil. We used a standard echo-planar imaging protocol (repetition time [TR] = 2000 ms; echo time [TE] = 30 ms; flip angle = 90°; matrix size = 64 x 64; field-of-view = 230 mm x 230 mm; slice thickness = 4.5 mm). A pulse oximeter (Nonin 8600FO) with a sensor on the left index finger was used to record O2 saturation and heart rate, and the plethysmographic waveform (SaO2) was recorded at 1 kHz. For spatial localization, two high resolution, T1-weighted anatomical images were acquired with a magnetization prepared rapid acquisition gradient echo sequence (TR = 2200 ms; TE = 2.2 ms; inversion time = 900 ms; flip angle = 9°; matrix size = 256 × 256; field-of-view = 230 × 230 mm; slice thickness = 1.0 mm). These two scans were realigned and averaged for each participant to result in one anatomical reference.

Protocol

Participants were asked to refrain from coffee and other substances with stimulants for 12 hours prior to the study. While lying in a scanner, following a 1 minute baseline, participants performed four 16 s handgrips (80% subjective maximum grip strength) at 1 minute intervals against a squeeze ball without metal components, with one minute baseline after the fourth challenge. An air-filled plastic bag, connected to a pressure transducer, was placed in the participants’ right hand. During the practice period, participants briefly squeezed at 100% subjective maximum strength at least two times and then at 80% for a 10–20 seconds. In the scanner, a light signal was used to indicate the onset of each grip period. Participants were instructed to squeeze to maintain the 80% pressure upon seeing the light signal. Participants practiced the static handgrip exercise maneuver prior to scanning, both outside and supine inside the MRI scanner. At least 30 min of rest (structural scanning) separated the practice from the trial periods. A pressure signal was monitored to verify that all participants performed the four static handgrip exercise tasks at the correct time. Timing was synchronized to the fMRI scans.

Analysis: Physiology and participant characteristics

We measured HR from the SaO2 plethysmographic waveform using peak detection. The median HR was calculated over the 60 second baseline period immediately prior to the first handgrip. Age, BMI, HR, and sleep parameters for the OSA group were described and compared with control groups using an ANOVA model. The full assessment of beat-to-beat HR responses was presented earlier [3].

Analysis: MRI

We preprocessed the fMRI scans using SPM12 (https://www.fil.ion.ucl.ac.uk/spm). Images were realigned for motion correction, and linear detrended over each series. For each participant, scans were spatially normalized in two steps, first coregistering the mean fMRI to the T1 anatomical scan and then warping to the “VBM8” template in Montreal Neurological Institute (MNI) space based on the T1 “DARTEL” spatial normalization algorithm [55]. These steps resulted in all participants’ fMRI images being in the template space.

The five major gyri were parcellated from the average of the high-resolution T1-weighted scans: Three short (anterior) gyri and two long (posterior) gyri: anterior short gyrus (ASG), mid short gyrus (MSG), posterior short gyrus (PSG), anterior long gyrus (ALG), and posterior long gyrus (PLG). We included these regions as mask files in nifti format (S2 File). Two experienced research team members determined the parcellation based on manual tracing with reference to a brain atlas [56]. The regions were outlined in normalized space; although this approach is slightly less accurate than individual tracing, the resolution of the fMRI data (>50 mm3) relative to the anatomical scans (< 1 mm3) is such that any differences in accuracy would not be meaningful. Signal intensity changes over time were extracted from each voxel in each gyrus from the processed images. For each gyrus in each participant, a mean time trend over all voxels was then calculated. Time trends were converted to percent change relative to the mean of the 1-minute baseline period. For each participant, the signals from the four challenges were separated and averaged to create one single handgrip percent change time trend that was passed to the group level analysis. While this averaging could theoretically result in reduced sensitivity, in practice, the statistical approach we chose takes advantage of repeated measures, and could detect small effect sizes.

To assess posterior-anterior effects, signal intensity changes were calculated relative to those in the PLG. As discussed above, the importance of the anterior insula has been described in clinical and animal studies. Our previous work showed this anterior-specific role could be demonstrated by comparing the fMRI signal in the anterior vs poster-most insula (PLG), and we repeated this technique here [46]. At each time point, signal intensity changes within the PLG were subtracted from those in the ASG, MSG, PSG and ALG for each hemisphere so that direct comparisons between these regions and the PLG could be assessed. Lateralization was assessed by subtracting signal changes in each of the five left gyri from the corresponding gyri on the right side; for example, ASG laterality was calculated by subtracting the left ASG time trend from the right ASG time trend.

The resulting fMRI signals were assessed for within- and between-group differences using repeated measures ANOVA (RMANOVA). The analysis was implemented with SAS “proc mixed”, as described earlier [57, 58]. In brief, this approach assesses within-group changes and between-group differences over time, with each 2 sec time-point during and after the challenge assessed relative to baseline time-points. We applied the Tukey-Fisher criterion for multiple comparisons; that is, we assessed the overall model for significance (p ≤ 0.05), and then effects of interest (time, group by time), before considering individual time-points of difference. The latter tests are performed within the “proc mixed” procedure, as the output includes time-point tests of significance (hence no post-hoc tests were needed). We assessed the effects in combined male-female models with sex as a covariate and in sex-specific models.

The RMANOVA mixed model approach allows for continuous variables to be included, so we performed secondary analyses of age and resting HR. We created four models that included different age effects added to the main model (group + time + group x time):

    Main + age: age effects independent of group over the entire protocol, independent of time;

    Main + age + age x group: group-specific effects of age over the entire protocol, independent of time;

    Main + age + age x time: age effects on handgrip responses, independent of group.

    Main + age + age x group + age x time + age x group x time: age effects on between-group differences in handgrip responses.

We repeated these calculations for models with HR in place of age. For the purposes of this study, we only focused on the within- and between-OSA and control group responses in the different models. The age-by-time and HR-by-time measures are not independent of the main effects of interest, but the degree to which these secondary models affect the within and between-group p-values reflects potential associations between the clinical and fMRI measures.

Results

Participants

Table 1 shows participant characteristics. Age was similar between OSA and control groups, and as expected, BMI was high in OSA over control. Unlike females, resting HR was higher in male OSA compared with male control participants.

Resting and evoked HR responses

The HR response during the baseline, handgrip and recovery for the 4 challenges have been shown previously [3]. Control and OSA groups displayed significant increases in HR during each of the 4 handgrip periods (Fig 1). Splitting by sex, male OSA subjects had significantly greater baseline HR compared with male controls, while female OSA subjects displayed reduced baseline HR levels. Furthermore, while male OSA subjects displayed robust increases in HR during each handgrip period, female OSA subjects show few such HR changes.

HR physiological changes moment to moment.
Fig 1

HR physiological changes moment to moment.

The HR (beats per minute) were illustrated with the mean (smooth lines) and stdev (faded lines) for control (grey) and OSA (black) groups in “all” (both sexes), males and females. The four handgrip challenges, along with the baseline and recovery during the experiment have been illustrated here. The first 60 seconds baseline denoted by the box was averaged to get the resting HR reported.

Intrinsic PLG changes in left and right insula

Table 2 and Fig 2 show the left and right insula intrinsic PLG responses. The overall model was not significant for females (model statistics in Table 2). For males, the model was significant, with a small decrease in signal intensity during each handgrip task. However, there were no within- or between-group effects of the handgrip task (Table 2).

PLG intrinsic changes in left and right insula.
Fig 2

PLG intrinsic changes in left and right insula.

Left and right hemisphere fMRI signals of PLG are demonstrated for combined (top panel), males (middle panel) and females (bottom panel). The graphs reflect baseline (group mean ± SE), averaged over challenges, with time-points of significant increase or decrease relative to baseline within-group, and time-points of between-group differences (RMANOVA, p < 0.05).

Table 2
Intrinsic PLG changes in the left and right insula.
p-values & model statistics for RMANOVAModel detailsAllMaleFemale
Left PLGRight PLGLeft PLGRight PLGLeft PLGRight PLG
Main model group, time group, timeChiSq (p value)64.1 (<0.001)52.35 (<0.001)82.84 (<0.001)62.78(<0.001)3.68 (0.30)2.6 (0.45)
Fit (−2 log-likelihood)23525.123750.316054.516214.973307390.3
Group effect p-value (mean over entire series for each group)0.730.710.920.87N/AN/A
Handgrip response: within groupTime (within-group effect of time) p-values
Main0.490.480.310.11N/AN/A
Age0.490.480.310.11N/AN/A
Age x Group0.490.480.310.11N/AN/A
Age x Time0.810.810.490.52N/AN/A
Age x Group x Time0.810.810.490.52N/AN/A
HR0.490.480.310.11N/AN/A
HR x Group0.490.480.310.11N/AN/A
HR x Time0.360.450.530.64N/AN/A
HR x Group x Time0.360.450.530.64N/AN/A
Handgrip response: between-groupTime X Group (between-group effect of time) p-values
Main0.510.540.280.32N/AN/A
Age0.510.540.280.32N/AN/A
Age x Group0.510.540.280.32N/AN/A
Age x Time0.500.520.280.32N/AN/A
Age x Group x Time0.500.520.280.32N/AN/A
HR0.510.540.280.32N/AN/A
HR x Group0.510.540.280.32N/AN/A
HR x Time0.660.690.600.64N/AN/A
HR x Group x Time0.660.690.600.64N/AN/A

Salient statistics and p-values from 9 RMANOVA models for left and right PLG in three sets (mixed, male, female). Full data are available online [59]. The main model (bold) is the interaction of group-by-time (fMRI = group + time + group x time), and statics of significance and fit are in the top rows of the table. The “Group” effect is the mean over the entire series and does not represent responses, and is not discussed. The two effects of interest “Time”, which represents within-group responses over time, and “Time x Group”, which represents between-group differences in responses. The p-values for these effects are shown for the 9 models. All models include the main effects plus additional mean or interaction terms. All interaction models also include means. For example, “Age x Time” is fMRI = group + time + group x time + age + age x time.

Laterality

Fig 3 shows the laterality of the insula responses during the handgrip, with positive values reflecting right dominance. Table 3 shows the model statistics. The right-minus-left percentage change fMRI signal for each gyrus is shown for all participants, and separately for males and females. In the combined participants, and in the females there was no group difference in the response of the right compared with the left during the handgrip for any insula gyri. That is, the overall model statistics in Table 3 show that no between-group differences appeared in the combined groups and in females. However, males showed a significant difference in OSA vs control for the ALG insular gyrus. At the ASG anterior gyrus, there were significant right activation effects by 0.2% in females, whereas males showed a significant left activation by 0.2% at the initial 8s of sympathetic activation phase for handgrip. However, these opposite lateralization effects were not significantly different in OSA vs controls. The influence of age or HR on the main effects was noted where p-values changed substantially. These covariates did modify the time it took for the subject to respond to the handgrip task.

Lateralization by gyri.
Fig 3

Lateralization by gyri.

Right hemisphere fMRI signals relative to left hemisphere for all gyri (group mean ± SE), averaged over challenges, with time-points of significant increase or decrease relative to baseline within-group, and time-points of between-group differences (RMANOVA, p < 0.05) for the right-left laterality insular gyri responses in participants (all in left, male in middle and females in right columns respectively).

Table 3
Laterality of insular fMRI organization.
Model detailsAllMaleFemale
ASGMSGPSGALGPLGASGMSGPSGALGPLGASGMSGPSGALGPLG
Main model: group, timeChiSq (p value)227.1 (<0.001)83.1 (<0.001)180.6 (<0.001)121.1 (<0.001)113.7 (<0.001)193.5 (<0.001)82.2 (<0.001)152.7 (<0.001)95.0 (<0.001)111.6 (<0.001)16.7 (<0.001)10.6 (0.01)22.8 (<0.001)31.5 (<0.001)7.0 (0.07)
Fit (−2 log-likelihood)5461.115181.710712.87682.99458.74077.910590.377045541.26691.81389.74569.93012.622042834.4
Group effect p-value (mean over entire series for each group)0.630.780.470.090.930.820.930.460.0130.870.070.460.370.16N/A
Handgrip response: within-groupTime (within-group effect of time) p-values
Main<0.0010.650.67<0.001<0.01<0.0010.520.500.0030.1<0.010.930.450.22N/A
Age<0.0010.650.67<0.001<0.01<0.0010.520.500.0030.1<0.010.930.450.22N/A
Age x Group<0.0010.650.67<0.001<0.01<0.0010.520.500.0030.1<0.010.930.450.22N/A
Age x Time0.530.310.090.470.810.550.0480.120.150.720.070.390.300.06N/A
Age x Group x Time0.530.310.090.470.810.550.0480.120.150.720.070.390.300.06N/A
HR<0.0010.650.67<0.001<0.01<0.0010.520.500.0030.1<0.010.930.450.22N/A
HR x Group<0.0010.650.67<0.001<0.01<0.0010.520.500.0030.1<0.010.930.450.22N/A
HR x Time0.410.220.210.350.390.340.130.600.370.750.660.890.720.94N/A
HR x Group x Time0.410.220.210.350.390.340.130.600.370.750.660.890.720.94N/A
Handgrip response: between-groupTime X Group (between-group effect of time) p-values
Main0.600.510.640.120.390.210.200.960.0160.850.780.930.660.83N/A
Age0.600.510.640.120.390.210.200.960.0160.850.780.930.660.83N/A
Age x Group0.600.510.640.120.390.210.200.960.0160.850.780.930.660.83N/A
Age x Time0.620.500.580.130.360.230.220.940.0150.830.760.940.640.78N/A
Age x Group x Time0.620.500.580.130.360.230.220.940.0150.830.760.940.640.78N/A
HR0.600.510.640.120.390.210.200.960.0160.850.780.930.660.83N/A
HR x Group0.600.510.640.120.390.210.200.960.0160.850.780.930.660.83N/A
HR x Time0.830.710.710.190.500.640.500.990.100.940.600.870.460.86N/A
HR x Group x Time0.830.710.710.190.500.640.500.990.100.940.600.870.460.86N/A

Salient statistics and p-values from 9 RMANOVA models for right-minus-left insula organization in three sets (mixed, male, female). Full data are available online [59]. The main model (bold) is the interaction of group-by-time (fMRI = group + time + group x time), and statics of significance and fit are in the top rows of the table. The “Group” effect is the mean over the entire series and does not represent responses, and is not discussed. The two effects of interest “Time”, which represents within-group responses over time, and “Time x Group”, which represents between-group differences in responses. The p-values for these effects are shown for the 9 models. All models include the main effects plus additional mean or interaction terms. All interaction models also include means. For example, “Age x Time” is fMRI = group + time + group x time + age + age x time.

Left side anterior-posterior organization

Fig 4 represents the left insula anterior-to-posterior functional organization during handgrip, with positive values representing greater anterior dominance. Table 4 shows the model statistics. The mean±SEM percentage change fMRI signal for each of left ASG, MSG, PSG and ALG with respect to PLG is shown for all participants, and separately for males and females. The signal intensity in the left ASG changes during handgrip were approximately 0.5% higher relative to those changes in the left PLG in both control and OSA groups for both males and females. The overall model statistics in Table 4 show that no between-group differences between controls and OSA appeared in the combined, males or females groups for the left ASG, MSG, PSG or ALG. The males and females showed sympathetic activation effects in the initial sympathetic phase of handgrip at the left insula. The influence of age or HR on the main effects was noted where p-values changed substantially.

Anterior-posterior organization of left insula.
Fig 4

Anterior-posterior organization of left insula.

Left hemisphere fMRI signals relative to PLG such that positive change reflects anterior dominance. Baseline (group mean ± SE), averaged over challenges, with time-points of significant increase or decrease relative to baseline within-group, and time-points of between-group differences (RMANOVA, p < 0.05) for the anterior-posterior responses by the left insular gyri in participants (all in left, male in middle and females in right columns respectively).

Table 4
Left insula anterior-posterior fMRI organization with respect to PLG.
Model detailsAllMaleFemale
ASGMSGPSGALGASGMSGPSGALGASGMSGPSGALG
Main model: group, timeChiSq (p value)72.54 (<0.001)66.82 (<0.001)52 (<0.001)55.07 (<0.001)68.07 (<0.001)56.94 (<0.001)39.3 (<0.001)43.18 (<0.001)12.98 (<0.01)13.46 (<0.01)16.13 (0.001)15 (0.002)
Fit (−2 log-likelihood)13329.121139.526125.226548.59352.314636.218108.818403.83973.16411.77852.47971
Group effect p-value (mean over entire series for each group)0.830.720.970.730.920.920.880.760.850.330.620.89
Handgrip response: within-groupTime (within-group effect of time) p-values
Main<0.0010.350.810.81<0.0010.320.770.680.0540.570.770.78
Age<0.0010.350.810.81<0.0010.320.770.680.0540.570.770.78
Age x Group<0.0010.350.810.81<0.0010.320.770.680.0540.570.770.78
Age x Time0.260.550.550.260.470.360.160.200.500.440.380.30
Age x Group x Time0.260.550.550.260.470.360.160.200.500.440.380.30
HR<0.0010.350.810.81<0.0010.320.770.680.0540.570.770.78
HR x Group<0.0010.350.810.81<0.0010.320.770.680.0540.570.770.78
HR x Time0.680.310.360.400.830.230.380.380.320.810.840.78
HR x Group x Time0.680.310.360.400.830.230.380.380.320.810.840.78
Handgrip response: between-groupTime X Group (between-group effect of time) p-values
Main0.810.400.600.600.970.320.800.840.440.910.680.57
Age0.810.400.600.600.970.320.800.840.440.910.680.57
Age x Group0.810.400.600.600.970.320.800.840.440.910.680.57
Age x Time0.570.380.380.570.970.350.800.840.420.930.690.58
Age x Group x Time0.570.380.380.570.970.350.800.840.420.930.690.58
HR0.810.400.600.600.970.320.800.840.440.910.680.57
HR x Group0.810.400.600.600.970.320.800.840.440.910.680.57
HR x Time0.900.570.730.730.990.590.950.960.210.820.420.27
HR x Group x Time0.900.570.730.730.990.590.950.960.210.820.420.27

Salient statistics and p-values from 9 RMANOVA models for left insula anterior-posterior organization in three sets, namely “all” (both sexes), female and male. Full data are available online [59]. The main model (bold) is the interaction of group-by-time (fMRI = group + time + group x time), and measures of significance and fit are in the top rows of the table. The “Group” effect is the mean over the entire series and does not represent responses, and is not discussed. The two effects of interest “Epoch”, which represents within-group responses over time, and “Epoch x Group”, which represents between-group differences in responses. The p-values for these effects are shown for the 9 models. All models include the main effects plus additional mean or interaction terms. All interaction models also include means. For example, “Age x Time” is fMRI = group + time + group x time + age + age x time.

Right side anterior-posterior organization

Fig 5 represents the right insula anterior-to-posterior functional organization during handgrip, with positive values representing greater anterior dominance. Table 5 shows the model statistics. The percentage change fMRI signal for each of right ASG, MSG, PSG and ALG with respect to right PLG is shown for all participants, and separately for males and females. Similar to the left side, while signal intensity during a handgrip was greater in right ASG (by 0.4% in females and by 0.2% in males), relative to the right PLG, there were no significant differences between the control and OSA groups. The other three gyri did not respond significantly to the handgrip task. Table 5 shows that no between-group differences appeared in the combined, male or female groups. However, unlike the left insula, age did not impact the right insula’s anterior dominance of ASG during the handgrip task in females. Although in males, the interaction of age by time (meaning age influences on handgrip responses) affected all within-group effects at the ASG anterior gyrus, but age had no influence on any between-group effects. The interaction of resting HR and time (meaning HR influence on handgrip responses) influenced within-group effects for females and males in the ASG.

Anterior-posterior organization of right insula.
Fig 5

Anterior-posterior organization of right insula.

Right hemisphere fMRI signals relative to PLG such that positive change reflects anterior dominance. Baseline (group mean ± SE), averaged over challenges, with time-points of significant increase or decrease relative to baseline within-group, and time-points of between-group differences (RMANOVA, p < 0.05) for the anterior-posterior responses by the right insular gyri in participants (all in left, male in middle and females in right columns respectively).

Table 5
Right insula anterior-posterior fMRI organization with respect to PLG.
Model detailsAllMaleFemale
ASGMSGPSGALGASGMSGPSGALGASGMSGPSGALG
Main model: group, timeChiSq (p value)205.5 (<0.001)76.74 (<0.001)58.29 (<0.001)58.18 (<0.001)175.52 (<0.001)59.34 (<0.001)48.71 (<0.001)48.1 (<0.001)29.72 (<0.001)29.33 <0.001)15.75 (0.001)16.9 (<0.001)
Fit (−2 log-likelihood)3633.61566024697.828268.32740.110938.217143.619610.71007.24661.57392.78443.7
Group effect p-value (mean over entire series for each group)0.400.720.760.550.470.720.610.440.250.200.680.67
Handgrip response: within-groupTime (within-group effect of time) p-values
Main<0.001<0.0010.640.84<0.0010.0160.590.59<0.001<0.0010.450.76
Age<0.001<0.0010.640.84<0.0010.0160.590.59<0.001<0.0010.450.76
Age x Group<0.001<0.0010.640.84<0.0010.0160.590.59<0.001<0.0010.450.76
Age x Time0.020.430.200.350.140.560.140.320.0060.430.320.25
Age x Group x Time0.020.430.200.350.140.560.140.320.0060.430.320.25
HR<0.001<0.0010.640.84<0.0010.0160.590.59<0.001<0.0010.450.76
HR x Group<0.001<0.0010.640.84<0.0010.0160.590.59<0.001<0.0010.450.76
HR x Time0.850.620.840.730.840.730.730.720.730.170.540.75
HR x Group x Time0.850.620.840.730.840.730.730.720.730.170.540.75
Handgrip response: between-groupTime X Group (between-group effect of time) p-values
Main0.850.570.690.670.800.510.800.900.920.710.460.36
Age0.850.570.690.670.800.510.800.900.920.710.460.36
Age x Group0.850.570.690.670.800.510.800.900.920.710.460.36
Age x Time0.850.550.660.650.780.540.810.910.860.750.470.35
Age x Group x Time0.850.550.660.650.780.540.810.910.860.750.470.35
HR0.850.570.690.670.800.510.800.900.920.710.460.36
HR x Group0.850.570.690.670.800.510.800.900.920.710.460.36
HR x Time0.920.740.810.780.940.770.950.980.730.520.200.13

Salient statistics and p-values from 9 RMANOVA models for right insula anterior-posterior organization in three sets (mixed, male, female). Full data are available online [59]. The main model (bold) is the interaction of group-by-time (fMRI = group + time + group x time), and statics of significance and fit are in the top rows of the table. The “Group” effect is the mean over the entire series and does not represent responses, and is not discussed. The two effects of interest “Time”, which represents within-group responses over time, and “Time x Group”, which represents between-group differences in responses. The p-values for these effects are shown for the 9 models. All models include the main effects plus additional mean or interaction terms. All interaction models also include means. For example, “Age x Time” is fMRI = group + time + group x time + age + age x time.

The full model results, including averaged timetrends as presented in the Figs, are available in a data repository [59].

Consistent with prior studies in healthy people [46], we showed functional anterior-posterior and right-left organization of functional responses of gyri within the insula to a handgrip autonomic challenge. We observed minimal differences between OSA and control participants, suggesting that previously-described insular dysfunction occurs across the whole structure [15, 2427]. Sex differences were apparent as previously shown, with anterior-most gyri exhibiting enhanced right-sided activation in females, and greater left-side activation in males [47]. We did not see substantial differences between OSA and controls in either combined-sex or separate male and female analyses. Age and resting HR showed associations with the magnitude of fMRI response to handgrip, but these associations did not change the patterns of gyral organization in OSA over control.

Discussion

Increased anterior insular activity during the handgrip occurs in OSA and healthy groups, consistent with the established autonomic function of this sub-region found in both human neuroimaging and animal lesion studies [46, 48]. The absence of substantial OSA-related differences within the insula occurred despite other studies showing functional variations in networks containing the structure [32, 33, 38]. Together with existing evidence, the present findings therefore suggest that within-structure organization of the insula is intact in OSA, but function of the structure as a whole is impaired.

The insula plays a key role in cardiovascular regulation, and hence is relevant to this main comorbidity of OSA. Primate studies have revealed that the anterior, agranular insula receives inputs from brainstem regions such as the nucleus of the solitary tract (NTS) and sends direct projections to caudal regions including the midbrain periaqueductal gray, parabrachial nucleus and NTS [60]. A motor command alone can evoke increased heart rate and blood pressure [61], even if the target muscle is paralyzed [62, 63], and animal electrophysiological studies have revealed that unmyelinated and small-diameter myelinated nerve fibers are responsible for muscle contraction-evoked cardiovascular and respiratory responses [64, 65]. Furthermore, these muscle afferents project to the premotor sympathetic neurons in the rostral ventrolateral medulla (RVLM) via the NTS [66], and healthy individuals show increased signal intensity changes in the parabrachial nucleus, NTS and RVLM during static handgrip [67]. As discussed previously [68], the insula modulates autonomic function via projections to autonomic outflow regions in the brainstem, often indirectly via the hypothalamus [21]. The hypothalamus, in particular, is tightly linked with sympathetic activity and has direct projections to and from the insula [19, 6972]. This close connection pattern is consistent with resting state neuroimaging studies that revealed altered baseline function in insular cortices and altered functional connectivity from the insula to other autonomic-related brain regions in OSA [25, 26], so a next area of investigation might be studies of functional interactions between insular, hypothalamic and brainstem activity during autonomic responses.

The magnitude of the fMRI responses to the task and of group differences was in the order of 0.3%, as seen in the time trends. For neuroimaging studies generally, a change of 1% in the BOLD signal is considered close to maximal, and the effects seen here are consistent with moderate effect sizes [73, 74]. Additionally, some effect sizes were minimal as the hypothesized differences were not observed; these negative findings were included in the results.

While a role for the brainstem in generating cardiovascular changes to muscle contractions is clear, cardiovascular changes can also occur during imagination of muscle contraction, that is without peripheral input [75]. The insular cortex and anterior cingulate cortex (ACC) have been implicated in integrating motor command and cardiovascular changes during exercise, which is supported by numerous brain imaging investigations [67, 7577]. The anterior insular also receives inputs from and projects to the ACC, and cardiovascular changes during imagined movement are apparently associated with activity changes in the ACC and insular cortex [75]. The findings from previous studies are consistent with the data presented here; that is, during a handgrip challenge, the increase in HR is associated with signal intensity changes within the anterior agranular insula and not in other more-posterior regions of the structure.

While only limited evidence exists on neural function-related sex differences in OSA, here, we found that both OSA and control females showed higher anterior signal dominance in both left and right insulae during an autonomic challenge. However, relative to control participants, OSA males showed left sided anterior dominance. This influence did not appear to be modulated by resting HR. However, resting HR was high in male, but not female in OSA relative to healthy groups, which could reflect different resting sympathetic tone, and potentially contribute to a sex-specific ceiling effect in fMRI responses. Given previous resting-state and neurotransmitter findings of the insula in OSA [25, 26, 3537], the question arises for future studies whether the baseline neural state is altered in a different manner in OSA females and males. Estrogen exerts sympathetic influences on the insula [78]. Furthermore, these effects appear to be mediated by GABA, and given the lower GABA in OSA, complex interactions between estrogen and OSA-related GABA reductions affecting autonomic regulation by the insula could arise. This left-sided parasympathetic and right-sided sympathetic insular laterality is more a bias than an absolute distinction between sympathetic and parasympathetic activity, since there are common cardiovascular responses to stimulation across multiple insular regions [50, 79]. Oppenheimer and colleagues [51] showed lateralization in insular cortex stimulation-elicited differential cardiovascular rhythm changes in epileptic patients, with right insula stimulation triggering sympathetic and left insula, parasympathetic effects. Removal of the right insula in rats leads to increased parasympathetic activity [80]. Under these assumptions, the findings suggest that males showed more parasympathetic withdrawal and females showed more sympathetic activation during the handgrip task, consistent with previous findings [47]. Consistent with our study in healthy individuals [46], here we report that the direction of left-right organization is similar in both OSA and control groups, with higher activity on the right side only for the females and not males.

Interpretation of these findings is limited by several factors. The handgrip task was short, and likely did not elicit metabolic effects or the increase in sympathetic activity found in healthy people. The grip strength was a subjective rating, but a percentage of maximum is the standard approach for handgrip tasks; this issue may have added variability to the data. The sample was originally powered for OSA-control differences, but not sex-specific effects; nevertheless, since we know there are sex differences in OSA and healthy cardiovascular and neural function, we decided to provide sex-specific results. Another limitation due to the original study intent is the lack of information in females on menstrual cycle, menopausal status, or use of hormones; these factors are all associated with autonomic influences [81, 82]. Sleep study parameters were based on the sleep study reports provided, and did not typically include quantification of awakenings, sleep onset latency, sleep efficiency or other parameters. Sleep quality in particular has been associated with changes in insular function [83]. With respect to neuroimaging, insula cortex folding is not constant across people, and the gyri are likely accurate only to within a few mm; fMRI signals are also limited in spatial accuracy due to the diffuse nature of the BOLD effect. Finally, representation of the OSA group to the broader patient population is limited.

In conclusion, functional response organization of gyri within the insular cortex is not substantially altered during a handgrip challenge in OSA subjects. This finding is similar to the lack of substantial OSA differences in those gyri in response to a Valsalva maneuver, and suggests that the anterior and right dominance of responses within the insula to autonomic stimuli may remain largely intact in the condition. It appears that in response to the motor driven handgrip task the autonomic functional organization of insular gyri appear muted and only the central command motor task aspect of the insula appears to differ in male OSA subjects compared to controls. Females showed higher anterior and right fMRI signal dominance in insula gyri compared with males, but the sample was insufficiently large to generalize with confidence. Since central autonomic regulation is impaired in OSA, given the peripheral weakened responses, questions remain regarding the resting state functional activity and connectivity with other autonomic regions such as the hypothalamus and brainstem.

References

JRTietjens, DClaman, EJKezirian, TDe Marco, AMirzayan, BSadroonri, et al Obstructive Sleep Apnea in Cardiovascular Disease: A Review of the Literature and Proposed Multidisciplinary Clinical Management Strategy. J Am Heart Assoc. 2019;8(1):e010440 Epub 2018/12/29. 10.1161/JAHA.118.010440

GLabarca, JDreyse, LDrake, JJorquera, FBarbe. Efficacy of continuous positive airway pressure (CPAP) in the prevention of cardiovascular events in patients with obstructive sleep apnea: Systematic review and meta-analysis. Sleep Med Rev. 2020;52:101312 Epub 2020/04/06. 10.1016/j.smrv.2020.101312 .

PMMacey, RKumar, MAWoo, FLYan-Go, RMHarper. Heart Rate Responses to Autonomic Challenges in Obstructive Sleep Apnea. Plos One. 2013;8(10). 10.1371/journal.pone.0076631

JMParish, VKSomers. Obstructive sleep apnea and cardiovascular disease. Mayo Clin Proc. 2004;79(8):103646. Epub 2004/08/11. 10.4065/79.8.1036 .

ZNestor, SSiddharth, RYBaba, NShah. The Influence of Sex on the Sleep-Cardiovascular Disease Relationship: a Review. Current Epidemiology Reports. 2015;2(4):25762. 10.1007/s40471-015-0058-0

BMokhlesi, SAHam, DGozal. The effect of sex and age on the comorbidity burden of OSA: an observational analysis from a large nationwide US health claims database. Eur Respir J. 2016;47(4):11629. Epub 2016/01/23. 10.1183/13993003.01618-2015 .

MDFaux, EKLarkin, BDHoit, JEAylor, ATWright, SRedline. Sex influences endothelial function in sleep-disordered breathing. Sleep. 2004;27(6):111320. 10.1093/sleep/27.6.1113

ARandby, SKNamtvedt, HHrubos-Strom, GEinvik, VKSomers, TOmland. Sex-Dependent Impact of OSA on Digital Vascular Function. Chest. 2013;144(3):91522. 10.1378/chest.12-2283

SHHan, HJKim, SALee. The effect of high evening blood pressure on obstructive sleep apnea-related morning blood pressure elevation: does sex modify this interaction effect? Sleep Breath. 2019;23(4):125563. Epub 2019/06/12. 10.1007/s11325-019-01869-5 .

10 

RMHarper, PMMacey, LAHenderson, MAWoo, KEMacey, RCFrysinger, et al fMRI responses to cold pressor challenges in control and obstructive sleep apnea subjects. Journal of Applied Physiology. 2003;94(4):158395. 10.1152/japplphysiol.00881.2002

11 

PMMacey, KEMacey, LAHenderson, JRAlger, RCFrysinger, MAWoo, et al Functional magnetic resonance imaging responses to expiratory loading in obstructive sleep apnea. Respiratory Physiology & Neurobiology. 2003;138(2):27590. 10.1016/j.resp.2003.09.002.

12 

LAHenderson, MAWoo, PMMacey, KEMacey, RCFrysinger, JRAlger, et al Neural responses during Valsalva maneuvers in obstructive sleep apnea syndrome. Journal of Applied Physiology. 2003;94(3):106374. 10.1152/japplphysiol.00702.2002

13 

KEMacey, PMMacey, MAWoo, LAHenderson, RCFrysinger, RKHarper, et al Inspiratory loading elicits aberrant fMRI signal changes in obstructive sleep apnea. Respiratory Physiology & Neurobiology. 2006;151(1):4460. 10.1016/j.resp.2005.05.024.

14 

PMMacey, RKumar, JAOgren, MAWoo, RMHarper. Global Brain Blood-Oxygen Level Responses to Autonomic Challenges in Obstructive Sleep Apnea. Plos One. 2014;9(8). 10.1371/journal.pone.0105261

15 

LAHenderson, MAWoo, PMMacey, KEMacey, RCFrysinger, JRAlger, et al Neural responses during Valsalva maneuvers in obstructive sleep apnea syndrome. J Appl Physiol (1985). 2003;94(3):106374. Epub 2002/11/16. 10.1152/japplphysiol.00702.2002 .

16 

DRDancey, PJHanly, CSoong, BLee, VHoffstein. Impact of menopause on the prevalence and severity of sleep apnea. Chest. 2001;120(1):1515. 10.1378/chest.120.1.151 .

17 

RBFogel, AMalhotra, GPillar, SDPittman, ADunaif, DPWhite. Increased prevalence of obstructive sleep apnea syndrome in obese women with polycystic ovary syndrome. J Clin Endocrinol Metab. 2001;86(3):117580. Epub 2001/03/10. 10.1210/jcem.86.3.7316 .

18 

TSaaresranta, OPolo. Sleep-disordered breathing and hormones. Eur Respir J. 2003;22(1):16172. Epub 2003/07/29. 10.1183/09031936.03.00062403 .

19 

DOngur, XAn, JLPrice. Prefrontal cortical projections to the hypothalamus in macaque monkeys. J Comp Neurol. 1998;401(4):480505. .

20 

PYRisold, RHThompson, LWSwanson. The structural organization of connections between hypothalamus and cerebral cortex. Brain Res Brain Res Rev. 1997;24(2–3):197254. Epub 1997/12/31. 10.1016/s0165-0173(97)00007-6 .

21 

CBSaper. Convergence of autonomic and limbic connections in the insular cortex of the rat. J Comp Neurol. 1982;210(2):16373. 10.1002/cne.902100207 .

22 

MNagai, SHoshide, KKario. The insular cortex and cardiovascular system: a new insight into the brain-heart axis. Journal of the American Society of Hypertension: JASH. 2010;4(4):17482. Epub 2010/07/27. 10.1016/j.jash.2010.05.001 .

23 

DARuggiero, SMraovitch, ARGranata, MAnwar, DJReis. A role of insular cortex in cardiovascular function. J Comp Neurol. 1987;257(2):189207. Epub 1987/03/08. 10.1002/cne.902570206 .

24 

KSTaylor, PJMillar, HMurai, NHaruki, DSKimmerly, TDBradley, et al Cortical autonomic network gray matter and sympathetic nerve activity in obstructive sleep apnea. Sleep. 2018;41(2). Epub 2018/01/09. 10.1093/sleep/zsx208 .

25 

BPark, JAPalomares, MAWoo, DWKang, PMMacey, FLYan-Go, et al Aberrant Insular Functional Network Integrity in Patients with Obstructive Sleep Apnea. Sleep. 2016;39(5):9891000. 10.5665/sleep.5738

26 

QZhang, WQin, XHe, QLi, BChen, YZhang, et al Functional disconnection of the right anterior insula in obstructive sleep apnea. Sleep Med. 2015;16(9):106270. 10.1016/j.sleep.2015.04.018 .

27 

DCPeng, XJDai, HHGong, HJLi, XNie, WZhang. Altered intrinsic regional brain activity in male patients with severe obstructive sleep apnea: a resting-state functional magnetic resonance imaging study. Neuropsychiatr Dis Treat. 2014;10:181926. 10.2147/NDT.S67805

28 

RKumar, ASChavez, PMMacey, MAWoo, FLYan-Go, RMHarper. Altered global and regional brain mean diffusivity in patients with obstructive sleep apnea. Journal of Neuroscience Research. 2012;90(10):204352. 10.1002/jnr.23083

29 

STummala, BRoy, RVig, BPark, DWKang, MAWoo, et al Non-Gaussian Diffusion Imaging Shows Brain Myelin and Axonal Changes in Obstructive Sleep Apnea. J Comput Assist Tomogr. 2017;41(2):1819. 10.1097/RCT.0000000000000537

30 

PMMacey, RKumar, MAWoo, EMValladares, FLYan-Go, RMHarper. Brain structural changes in obstructive sleep apnea. Sleep. 2008;31(7):96777. Epub 2008/07/26. 10.5665/sleep/31.7.967

31 

PMMacey, ARomigi, NHaris, RKumar, MAThomas, MAWoo, et al Obstructive sleep apnea and cortical thickness in females and males. PLoS One. 2018 10.1371/journal.pone.0193854

32 

AWKYeung. Morphometric and functional connectivity changes in the brain of patients with obstructive sleep apnea: A meta-analysis. J Sleep Res. 2019;28(6):e12857 Epub 2019/04/16. 10.1111/jsr.12857 .

33 

MTahmasian, IRosenzweig, SBEickhoff, AASepehry, ARLaird, PTFox, et al Structural and functional neural adaptations in obstructive sleep apnea: An activation likelihood estimation meta-analysis. Neurosci Biobehav Rev. 2016;65:14256. Epub 2016/04/04. 10.1016/j.neubiorev.2016.03.026

34 

CRInnes, PTKelly, MHlavac, TRMelzer, RDJones. Decreased Regional Cerebral Perfusion in Moderate-Severe Obstructive Sleep Apnoea during Wakefulness. Sleep. 2015;38(5):699706. 10.5665/sleep.4658

35 

PMMacey, MKSarma, RNagarajan, RAysola, JMSiegel, RMHarper, et al Obstructive sleep apnea is associated with low GABA and high glutamate in the insular cortex. J Sleep Res. 2016;25(4):3904. Epub 2016/02/05. 10.1111/jsr.12392

36 

SKYadav, RKumar, PMMacey, MAWoo, FLYan-Go, RMHarper. Insular Cortex Metabolite Changes in Obstructive Sleep Apnea. Sleep. 2014;37(5):9518. 10.5665/sleep.3668

37 

JKang, ZTian, MLi. Changes in insular cortex metabolites in patients with obstructive sleep apnea syndrome. Neuroreport. 2018;29(12):9816. Epub 2018/06/19. 10.1097/WNR.0000000000001065 .

38 

LTChen, XLFan, HJLi, CLYe, HHYu, HZXin, et al Aberrant brain functional connectome in patients with obstructive sleep apnea. Neuropsychiatr Dis Treat. 2018;14:105970. Epub 2018/05/02. 10.2147/NDT.S161085

39 

MGuenot, JIsnard, MSindou. Surgical anatomy of the insula. Adv Tech Stand Neurosurg. 2004;29:26588. 10.1007/978-3-7091-0558-0_7 .

40 

PMMacey, JAOgren, RKumar, RMHarper. Functional Imaging of Autonomic Regulation: Methods and Key Findings. Front Neurosci. 2015;9:513 Epub 2016/02/10. 10.3389/fnins.2015.00513

41 

TZiemssen, TSiepmann. The Investigation of the Cardiovascular and Sudomotor Autonomic Nervous System-A Review. Front Neurol. 2019;10:53 Epub 2019/02/28. 10.3389/fneur.2019.00053

42 

HAKluess, RHWood. Heart rate variability and the exercise pressor reflex during dynamic handgrip exercise and postexercise arterial occlusion. Am J Med Sci. 2005;329(3):11723. Epub 2005/03/16. 10.1097/00000441-200503000-00002 .

43 

BCMaciel, LGalloJr., JAMarin-Neto, JTerra-Filho, JCManco. Efficacy of pharmacological blockade of the cardiac parasympathetic system with atropine in normal men. Braz J Med Biol Res. 1985;18(3):3038. Epub 1985/01/01. .

44 

LGalloJr., BCMaciel, JAMarin-Neto, LEMartins, ECLima-Filho, JCManco. The use of isometric exercise as a means of evaluating the parasympathetic contribution to the tachycardia induced by dynamic exercise in normal man. Pflugers Arch. 1988;412(1–2):12832. Epub 1988/07/01. 10.1007/BF00583741 .

45 

JSSanders, ALMark, DWFerguson. Evidence for cholinergically mediated vasodilation at the beginning of isometric exercise in humans. Circulation. 1989;79(4):81524. Epub 1989/04/01. 10.1161/01.cir.79.4.815 .

46 

PMMacey, PWu, RKumar, JAOgren, HLRichardson, MAWoo, et al Differential responses of the insular cortex gyri to autonomic challenges. Auton Neurosci. 2012;168(1–2):7281. Epub 2012/02/22. 10.1016/j.autneu.2012.01.009

47 

PMMacey, NSRieken, JAOgren, KEMacey, RKumar, RMHarper. Sex differences in insular cortex gyri responses to a brief static handgrip challenge. Biol Sex Differ. 2017;8:13 Epub 2017/04/25. 10.1186/s13293-017-0135-9

48 

SOppenheimer, DCechetto. The Insular Cortex and the Regulation of Cardiac Function. Compr Physiol. 2016;6(2):1081133. 10.1002/cphy.c140076 .

49 

APHollander, LNBouman. Cardiac acceleration in man elicited by a muscle-heart reflex. Journal of Applied Physiology. 1975;38(2):2728. 10.1152/jappl.1975.38.2.272

50 

FChouchou, FMauguiere, OVallayer, HCatenoix, JIsnard, AMontavont, et al How the insula speaks to the heart: Cardiac responses to insular stimulation in humans. Hum Brain Mapp. 2019;40(9):261122. Epub 2019/03/01. 10.1002/hbm.24548

51 

SMOppenheimer, AGelb, JPGirvin, VCHachinski. Cardiovascular effects of human insular cortex stimulation. Neurology. 1992;42(9):172732. Epub 1992/09/01. 10.1212/wnl.42.9.1727 .

52 

FRMarins, MLimborco-Filho, BFD’Abreu, PWMachado de Almeida, MGavioli, CHXavier, et al Autonomic and cardiovascular consequences resulting from experimental hemorrhagic stroke in the left or right intermediate insular cortex in rats. Auton Neurosci. 2020;227:102695 Epub 2020/07/07. 10.1016/j.autneu.2020.102695 .

53 

RHFatouleh, EHammam, LCLundblad, PMMacey, DKMcKenzie, LAHenderson, et al Functional and structural changes in the brain associated with the increase in muscle sympathetic nerve activity in obstructive sleep apnoea. Neuroimage Clin. 2014;6:27583. Epub 2014/11/08. 10.1016/j.nicl.2014.08.021

54 

Sleep-related breathing disorders in adults: recommendations for syndrome definition and measurement techniques in clinical research. The Report of an American Academy of Sleep Medicine Task Force. Sleep. 1999;22(5):66789. Epub 1999/08/18. .

55 

JAshburner. A fast diffeomorphic image registration algorithm. Neuroimage. 2007;38(1):95113. 10.1016/j.neuroimage.2007.07.007 .

56 

Mai JK, Majtanik M, Paxinos G. Atlas of the human brain: Academic Press; 2015.

57 

Macey PM, Schluter PJ, Macey KE, Harper RM. Detecting variable responses within fMRI time-series of volumes-of-interest using repeated measures ANOVA [version 1; referees: awaiting peer review]. 2016.

58 

RCLittell, GAMilliken, WWStroup, RDWolfinger. SAS System for Mixed Models. Cary, NC: SAS Institute Inc.; 1996.

59 

Macey PM. Summary fMRI data from handgrip in OSA and control groups with separation by sex. DRAFT VERSION ed: Harvard Dataverse; 2020.

60 

KNakano, TKayahara, TChiba. Afferent connections to the ventral striatum from the medial prefrontal cortex (area 25) and the thalamic nuclei in the macaque monkey. Ann N Y Acad Sci. 1999;877:66770. Epub 1999/07/23. 10.1111/j.1749-6632.1999.tb09297.x .

61 

GMGoodwin, DIMcCloskey, JHMitchell. Cardiovascular and respiratory responses to changes in central command during isometric exercise at constant muscle tension. J Physiol. 1972;226(1):17390. 10.1113/jphysiol.1972.sp009979

62 

SCGandevia, SFHobbs. Cardiovascular responses to static exercise in man: central and reflex contributions. J Physiol. 1990;430:10517. 10.1113/jphysiol.1990.sp018284

63 

SCGandevia, KKillian, DKMcKenzie, MCrawford, GMAllen, RBGorman, et al Respiratory sensations, cardiovascular control, kinaesthesia and transcranial stimulation during paralysis in humans. J Physiol. 1993;470:85107. Epub 1993/10/01. 10.1113/jphysiol.1993.sp019849

64 

McCord JL, Kaufman MP. Reflex Autonomic Responses Evoked by Group III and IV Muscle Afferents. In: Kruger L, Light AR, editors. Translational Pain Research: From Mouse to Man. Frontiers in Neuroscience. Boca Raton (FL)2010.

65 

MNMurphy, MMizuno, JHMitchell, SASmith. Cardiovascular regulation by skeletal muscle reflexes in health and disease. Am J Physiol Heart Circ Physiol. 2011;301(4):H1191204. Epub 2011/08/16. 10.1152/ajpheart.00208.2011

66 

GAIwamoto, MPKaufman. Caudal ventrolateral medullary cells responsive to muscular contraction. J Appl Physiol (1985). 1987;62(1):14957. Epub 1987/01/01. 10.1152/jappl.1987.62.1.149 .

67 

MSander, VGMacefield, LAHenderson. Cortical and brain stem changes in neural activity during static handgrip and postexercise ischemia in humans. J Appl Physiol (1985). 2010;108(6):1691700. 10.1152/japplphysiol.91539.2008 .

68 

APal, JAOgren, APAguila, RAysola, RKumar, LAHenderson, et al Functional Organization of the Insula in Men and Women with Obstructive Sleep Apnea during Valsalva. Sleep. 2020 Epub 2020/06/28. 10.1093/sleep/zsaa124 .

69 

DFCechetto, SJChen. Subcortical sites mediating sympathetic responses from insular cortex in rats. Am J Physiol. 1990;258(1 Pt 2):R24555. Epub 1990/01/01. 10.1152/ajpregu.1990.258.1.R245 .

70 

VGMacefield, LAHenderson. Identification of the human sympathetic connectome involved in blood pressure regulation. Neuroimage. 2019;202:116119 Epub 2019/08/26. 10.1016/j.neuroimage.2019.116119 .

71 

VGMacefield, LAHenderson. "Real-time" imaging of cortical and subcortical sites of cardiovascular control: concurrent recordings of sympathetic nerve activity and fMRI in awake subjects. J Neurophysiol. 2016;116(3):1199207. Epub 2016/06/24. 10.1152/jn.00783.2015

72 

LAHenderson, CJames, VGMacefield. Identification of sites of sympathetic outflow during concurrent recordings of sympathetic nerve activity and fMRI. Anat Rec (Hoboken). 2012;295(9):1396403. 10.1002/ar.22513 .

73 

CMagri, USchridde, YMurayama, SPanzeri, NKLogothetis. The amplitude and timing of the BOLD signal reflects the relationship between local field potential power at different frequencies. J Neurosci. 2012;32(4):1395407. Epub 2012/01/27. 10.1523/JNEUROSCI.3985-11.2012

74 

NKLogothetis, JPauls, MAugath, TTrinath, AOeltermann. Neurophysiological investigation of the basis of the fMRI signal. Nature. 2001;412(6843):1507. 10.1038/35084005 .

75 

JWWilliamson, RMcColl, DMathews, JHMitchell, PBRaven, WPMorgan. Brain activation by central command during actual and imagined handgrip under hypnosis. J Appl Physiol (1985). 2002;92(3):131724. Epub 2002/02/14. 10.1152/japplphysiol.00939.2001 .

76 

JWWilliamson, ACNobrega, RMcColl, DMathews, PWinchester, LFriberg, et al Activation of the insular cortex during dynamic exercise in humans. J Physiol. 1997;503 (Pt 2):27783. Epub 1997/11/05. 10.1111/j.1469-7793.1997.277bh.x

77 

JWWilliamson, RMcColl, DMathews. Evidence for central command activation of the human insular cortex during exercise. J Appl Physiol (1985). 2003;94(5):172634. Epub 2003/01/21. 10.1152/japplphysiol.01152.2002 .

78 

TMSaleh, BJConnell, AECribb. Sympathoexcitatory effects of estrogen in the insular cortex are mediated by GABA. Brain Res. 2005;1037(1–2):11422. 10.1016/j.brainres.2005.01.010 .

79 

LMazzola, FMauguiere, JIsnard. Functional mapping of the human insula: Data from electrical stimulations. Revue Neurologique. 2019;175(3):1506. 10.1016/j.neurol.2018.12.003

80 

HMde Morree, GJRutten, BMSzabo, MMSitskoorn, WJKop. Effects of Insula Resection on Autonomic Nervous System Activity. J Neurosurg Anesthesiol. 2016;28(2):1538. 10.1097/ANA.0000000000000207 .

81 

HRMiddlekauff, JPark, JAGornbein. Lack of effect of ovarian cycle and oral contraceptives on baroreceptor and nonbaroreceptor control of sympathetic nerve activity in healthy women. Am J Physiol Heart Circ Physiol. 2012;302(12):H25606. Epub 2012/05/01. 10.1152/ajpheart.00579.2011

82 

JRSmith, KEKoepp, JDBerg, JGAkinsanya, TPOlson. Influence of Sex, Menstrual Cycle, and Menopause Status on the Exercise Pressor Reflex. Med Sci Sports Exerc. 2019;51(5):87481. Epub 2019/04/16. 10.1249/MSS.0000000000001877

83 

VGuadagni, FBurles, MFerrara, GIaria. Sleep quality and its association with the insular cortex in emotional empathy. Eur J Neurosci. 2018;48(6):2288300. Epub 2018/08/18. 10.1111/ejn.14124 .