Sex differences in forebrain and cardiovagal responses at the onset of isometric handgrip exercise: a retrospective fMRI study
Abstract
In general, cardiac regulation is dominated by the sympathetic and parasympathetic nervous systems in men and women, respectively. Our recent study had revealed sex differences in the forebrain network associated with sympathoexcitatory response to baroreceptor unloading. The present study further examined the sex differences in forebrain modulation of cardiovagal response at the onset of isometric exercise. Forebrain activity in healthy men (n = 8) and women (n = 9) was measured using functional magnetic resonance imaging during 5 and 35% maximal voluntary contraction handgrip exercise. Heart rate (HR), mean arterial pressure (MAP), and muscle sympathetic nerve activity (MSNA) were collected in a separate recording session. During the exercise, HR and MAP increased progressively, while MSNA was suppressed (P < 0.05). Relative to men, women demonstrated smaller HR (8 ± 2 vs. 18 ± 3 beats/min) and MAP (3 ± 2 vs. 11 ± 2 mmHg) responses to the 35% maximal voluntary contraction trials (P < 0.05). Although a similar forebrain network was activated in both groups, the smaller cardiovascular response in women was reflected in a weaker insular cortex activation. Nevertheless, men did not show a stronger deactivation at the ventral medial prefrontal cortex, which has been associated with modulating cardiovagal activity. In contrast, the smaller cardiovascular response in women related to their stronger suppression of the dorsal anterior cingulate cortex activity, which has been associated with sympathetic control of the heart. Our findings revealed sex differences in both the physiological and forebrain responses to isometric exercise.
studies in humans and rodents indicate important male-female differences in cardiovascular control (13, 26, 45). At least part of these differences has been related to differences in autonomic responses to physiological stress (8, 35, 45). Specifically, heart rate (HR) variability approaches have been interpreted to suggest a dominance of parasympathetic modulation in women and sympathetic dominance in men (18, 31, 32, 46). Also, direct measures of muscle sympathetic nerve activity (MSNA) generally indicate lower overall levels in young women vs. young men (16, 45), although some overlap exists among individuals. Recent efforts have been directed at understanding the brain stem and cortical networks that affect the autonomic nervous system (1, 10, 23, 27, 42) and how these might differ in men and women (29).
A possible mechanism affecting sex differences in cardiovascular control is estrogen. Estrogen attenuates cardiovascular responses to central command (22) and to exercise pressor reflex (43, 44). Also, direct injection of estrogen at the nucleus tractus solitarius enhanced parasympathetic activity and suppressed sympathetic outflow (41). Nevertheless, the impact of sex in the central autonomic network, which has a significant role in modulating autonomic outflow to the cardiovascular system, remains unclear.
Recent developments in neuroimaging have resulted in the understanding that cardiovascular arousal is associated with a network of cortical activation (10, 12, 27, 48, 50). Our laboratory's early work has shown that the cortical activation patterns differ to some extent for baroreflex (27), which emphasizes combined sympathetic and parasympathetic regulations vs. short-term handgrip (50), which emphasizes parasympathetic changes. Specifically, baroreceptor unloading was associated with activation changes within the insular cortex (IC), the anterior cingulate cortex (ACC), the medial prefrontal cortex (MPFC), and the amygdala. In contrast, a dominant feature of nonfatiguing handgrip exercise was reduced activation within the ventral MPFC (vMPFC) and the absence of ACC activation (50). Furthermore, previous studies from our laboratory have shown that, while the cortical autonomic network for baroreflex control was similar between male and female participants, the larger sympathetic response to baroreceptor unloading in men was associated with a greater magnitude of change in this cortical network (29).
Given that men and women differ in their autonomic response and that the cortical autonomic network varies between baroreceptor vs. handgrip exercise protocols, it is reasonable to postulate that sex differences should exist in the cortical organization for autonomic changes during the activation of skeletal muscle. The use of short-term handgrip protocols can be emphasized because of the emphasis on cardiovagal regulation rather than sympathoexcitation. As mentioned above, our laboratory's earlier study (50) used a short-term (e.g., 30 s) moderate-intensity handgrip protocol to elicit rapid changes in HR that were not associated with changes in peroneal nerve sympathetic nerve activity. This exercise model was chosen on the basis of previous studies, which showed that a brief isometric handgrip exercise elicits no change in sympathetic nerve activity (34), but a strong tachycardia that is not altered by β-adrenoceptor blockade but vagal blockade (19, 37, 47).
Presently, there is no information on the sex differences in cortical activation patterns associated with maneuvers that elicit predominant cardiovagal changes. Therefore, the present study aimed to examine sex differences in the forebrain organization associated with cardiovagal control using an isometric handgrip exercise. To do so, we employed a retrospective analysis of our laboratory's previous report (50), in which male-female comparisons were possible.
The specific objectives of this report were to determine 1) sex differences in the HR response to a brief isometric exercise; and 2) whether such sex differences would be associated with different forebrain organization for this maneuver that emphasizes cardiovagal control. The generalized understanding that women tend to exhibit a predominance of vagal influence over cardiovascular control suggests the hypothesis that women would evoke different neural responses in the forebrain cardiovagal network than men. Based on the previous inverse relationship between HR and vMPFC deactivation (50), we tested the additional a priori hypothesis that any sex differences in HR responses would be correlated to vMPFC differences. Measures of HR, mean arterial pressure (MAP), and MSNA were assessed to establish the cardiovascular arousal, whereas functional magnetic resonance imaging (fMRI) examined the cortical neural responses in premenopausal women and similarly aged men.
METHODS
Subjects.
Eight men and nine women were recruited for this study. All subjects had no prior history of cardiovascular disease or neurological disorders. A medical screening and magnetic resonance imaging preliminary questionnaire were given to each subject to ensure safe compatibility within a high magnetic field environment. Subjects were instructed to consume a light meal ∼2 h before the experiment and to abstain from food, nicotine, alcohol, caffeine, and intense physical exertion for 12 h before the experiment. Each subject provided informed consent to the study procedures, which had been approved by the Office of Research Ethics at The University of Western Ontario.
Experimental design.
The subjects participated in two separate experimental sessions: 1) the physiological recording (LAB) session, and 2) the neuroimaging (fMRI) session. The two sessions were separated by a minimum period of 1 wk and performed at the same time of day. Subjects were familiarized with the experimental procedures before their first test session. Each session began with a maximal voluntary contraction (MVC) handgrip calibration, in which the subjects were instructed to squeeze the handgrip with their right hand to their maximum ability for at least two times while in a supine position. The strongest contraction was calibrated as 100% MVC. Then subjects were asked to practice the exercise task by squeezing the handgrip to 5 and 35% of their MVC for several times, so that they could achieve the designated force as rapidly as possible. During the experiment, subjects were presented with visual feedback of their achieved force in real time; hence they could monitor and maintain the level of contraction in the exercise task.
Following development of a stable baseline HR and blood pressure, data were collected during a baseline period followed by three 30-s blocks of a specified handgrip exercise that were separated by 1 min of rest. The subject was cued verbally to start and stop the handgrip exercise. The handgrip protocol was repeated twice within each of two tensions (5 and 35% of maximal voluntary strength). The repeated trials were separated by 2 min of recovery. The level of physical stress of the exercise was monitored by asking the subjects to rate their perceived exertion after each run on a 6–20 scale (2).
Data acquisition and analysis in the LAB session.
HR was determined from successive R-wave-R-wave intervals obtained by standard three-lead ECG. Beat-by-beat measures of arterial blood pressure were obtained using finger photoplethysmographic techniques (Finapres 2300, Ohmeda, Englewood, CA) with the hand held at heart level. These blood pressure measures were corrected against sphygmomanometrically collected systolic and diastolic pressures that were made intermittently throughout the experiment. MAP was calculated as ⅔ diastolic blood pressure + ⅓ systolic blood pressure. The nonmagnetic handgrip device consisted of an inflated rubber bladder, connected in series with plastic tubing to a disposable pressure transducer (PX272, Edwards Lifesciences, Irvine, CA), and a bridge amplifier. Analog signals for the blood pressure and handgrip were sampled at 200 Hz, and ECG was sampled at 400 Hz with an on-line data-acquisition system (PowerLab, ADInstruments, Mountain View, CA).
Postganglionic MSNA was assessed using microneurographic techniques from the common peroneal nerve (15, 21). A tungsten microelectrode that tapered to an uninsulated 1- to 5-μm tip was inserted percutaneously into the right common peroneal nerve with a reference electrode positioned subcutaneously 1–3 cm from the recording site. Neural activity was amplified 1,000 times by a preamplifier and an additional 75,000 times through a variable-gain isolated amplifier. The signal was band-pass filtered (0.7–2.0 kHz), full-wave rectified, and integrated with a resistance-capacitance circuit (0.1-s time constant). Criteria for an acceptable MSNA recording included pulse synchrony with the cardiac cycle and increased activity to a voluntary apnea but not to emotional arousal to a loud noise.
Data analysis.
The HR and MAP data were averaged over 2.5-s bins (the time to repetition interval for functional scans) during 30 s of baseline and over each contraction. These values were time aligned to ensure a mean value that would correspond to each fMRI scan (i.e., 100 data points) obtained in the neuroimaging session. Changes in HR and MAP were calculated by subtracting the baseline value (i.e., the average of the 30 s immediately before the exercise) from each value. For each subject, the HR and MAP responses were averaged over the six repeated blocks in the two separate runs. HR and MAP changes at the 5 and 35% MVC trials were compared with a mixed-effect ANOVA for repeated measures. If significance of P < 0.05 was reached, Fisher's post hoc test was used to identify where the significant differences lay. All data are represented as means ± SE.
A detailed description of the MSNA analytic procedures has been reported previously (28). Briefly, only MSNA bursts with characteristic rising and falling slopes, and amplitudes that were at least twice that of the interburst baseline fluctuations (2:1 signal/noise), were included in the analysis. Adequate sites for measuring MSNA were obtained or maintained in seven of eight men and four of nine women. MSNA burst rate (frequency/minute) and MSNA burst height (amplitude/burst) during the baseline and exercise periods were compared with a mixed-effect ANOVA for repeated measures. Statistical significance was accepted at P < 0.05.
Data acquisition and analysis in the fMRI session.
All imaging data were collected on a 4-T whole body imaging system (Varian, Palo Alto, CA; Siemans, Erlangen, Germany). During the scanning session, HR was calculated from the pulse intervals recorded on a MRI-compatible Oximeter (8600FO MRI, Nonin Medical, Plymouth, MN) placed over the middle finger of the left hand. The handgrip device was the same as that described above.
Before imaging, a global shimming procedure (RASTAMAP) using first- and second-order shims was performed to optimize the magnetic field over the imaging volume of interest (30). Twenty-one interleaved contiguous axial slices (5 mm thick, 3 × 3 mm in-plane voxel resolution) were acquired in each volume. Volume acquisition time was 2.5 s, with a time to repetition of 0.6250 s (4 shots). A total of 100 volumes were collected per session. Five steady-state volumes, acquired before actual data collection to allow for magnetization equilibrium, were discarded before data analysis. Functional data were collected using a segmented T2*-weighted gradient-echo echo planar imaging pulse sequence (TE = 12 ms, flip angle = 45°, field of view = 192 × 192 mm) with navigator echo correction. A corresponding high-resolution T1-weighted structural volume was acquired at the beginning of the same scanning session with a voxel resolution of 0.75 × 0.75 × 2.5 mm. Participants were immobilized within the head cradle with foam padding and were instructed to avoid head movements during the scanning period.
All fMRI data were analyzed with SPM2 software (Wellcome Department of Imaging Neuroscience, London, UK). The functional echo planar images were realigned to correct head motion and normalized to the Montreal Neurological Institute template. The functional images were smoothed with an 8-mm full-width half-maximum Gaussian kernel. A high-pass filter with a cutoff at 128 s was applied to the images to reduce low-frequency noise. Serial correlations in the time series were modeled by an autoregressive model with white noise.
A two-level statistical paradigm was used for all functional imaging analysis. First, individual design matrices were constructed to analyze the subject-session interactions. The blood oxygenation level-dependent (BOLD) signal changes associated with cardiovascular responses during the exercise period were modeled with the time course of the averaged HR response across all of the subjects. General linear model was used to calculate the parameter estimates for all brain voxels (20). The resulting subject-specific contrast images were then entered into a mixed-effects ANOVA model for random-effects analysis. Statistical parametric maps for each group were generated at an uncorrected threshold of P < 0.005, with an extent threshold of 10 voxels and overlayed onto an averaged T1 structural image provided by SPM2. Based on our laboratory's previous report (50), the left motor cortex (MC) and the vMPFC were identified as a priori regions of interest. At the first-level individual analysis, significant activation and deactivation were observed in the MC and vMPFC, respectively, for every subject (P < 0.05, corrected for multiple comparisons). The adjusted BOLD signals of these two regions were extracted from the activation/deactivation peak of each subject and averaged within each group.
RESULTS
Cardiovascular and sympathetic responses to handgrip exercise.
The baseline physiological data are summarized in Table 1. One of the female subjects was excluded from the analysis, as her HR response was more than 4 SDs from the group mean. The male and female subjects were matched for age. Baseline HR values for women and men, respectively, were 59 ± 3 vs. 55 ± 2 beats/min. Similarly, MAP values for women (83 ± 2 mmHg) and men (81 ± 4 mmHg) were not different. Baseline MSNA burst frequency was 15 ± 3 bursts/min for men vs. 8.0 ± 2.6 burst/min in women (not significant). The baseline HR in the LAB and fMRI sessions were similar for both groups.
| Men | Women | P Value | |
|---|---|---|---|
| Age, yr | 25±1 | 24±1 | 0.54 |
| Weight, kg | 79±3 | 57±2 | <0.001 |
| fMRI HR, beats/min | 54±3 | 58±2 | 0.35 |
| LAB HR, beats/min | 55±2 | 59±3 | 0.28 |
| MAP, mmHg | 81±4 | 83±2 | 0.54 |
| MSNA burst frequency, burst/min | 15.4±2.9 | 8.0±2.6 | 0.12 |
| MSNA burst amplitude, mV | 0.51±0.27 | 0.15±0.03 | 0.35 |
For both men and women, HR responses to the handgrip exercise were similar in the LAB and fMRI sessions (Fig. 1). Neither MAP nor MSNA could be measured in the fMRI session. Nonetheless, based on similar HR responses, we expected that both groups had comparable MAP and MSNA responses during the LAB and fMRI sessions.

Fig. 1.Average heart rate (HR) and mean arterial pressure (MAP) responses during baseline and 30 s of isometric handgrip exercise. HR and MAP were elevated during the exercise, whereas the 35% maximal voluntary contraction (MVC) trials elicited greater HR and MAP responses relative to the 5% MVC trials (P < 0.001). During the 35% MVC trials, men demonstrated greater HR and MAP responses than women (P < 0.001). Comparable HR responses were elicited in the neuroimaging (fMRI) and physiological recording (LAB) sessions. The HR and MAP data were averaged over the six repeated blocks, and the shaded area represents the exercise period. M35, men 35% MVC trials; F35, women 35% MVC trials; M5, men 5% MVC trials; F5, women 5% MVC trials; Δ, change.
In both groups, HR and MAP increased progressively during the 35% trials and peaked at the end of the exercise period (Fig. 1). A main effect of exercise tension was observed for both HR and MAP, where the 35% trials elicited greater responses than the 5% trials (P < 0.001; Fig. 1). A significant group × intensities interaction was also observed for HR (P < 0.023) and MAP (P < 0.002). Post hoc analysis indicated that the HR and MAP responses were greater in men vs. women during the 35% trials (Fig. 1). For the 35% trials, the peak HR changes in men and women were 18 ± 3 and 8 ± 2 beats/min, respectively, whereas the peak MAP changes were 11 ± 2 and 3 ± 2 mmHg in men and women, respectively.
MSNA data were obtained in seven men and four women. Although MSNA burst amplitude was not different between the baseline and exercise period for both groups (Fig. 2), there was an overall decrease of MSNA burst frequency during the 35% MVC exercise (P < 0.01). Nevertheless, the sex × exercise interaction for MSNA burst frequency was P < 0.079.

Fig. 2.Muscle sympathetic nerve activity (MSNA) burst frequency and amplitude during the baseline and 35% MVC handgrip exercise. Although MSNA burst amplitude was not different between the baseline and exercise period for both groups, there was an overall decrease of MSNA burst frequency during the 35% MVC trials (P < 0.01). Nevertheless, the sex × exercise interaction for MSNA burst frequency was P < 0.079. au, Arbitrary units.
The Borg scale responses reflected the level of exertion to the exercise and were obtained from five women and eight men in the LAB session. Consistent ratings were obtained from the men and women. Men and women rated, respectively, the 5% trials as 6 ± 0.1 and 6 ± 0.2 and the 35% trials as 11 ± 0.3 and 11 ± 0.7. These ratings indicated that both the 5 and 35% handgrip contractions were nonfatiguing or “very light” to “light” exercise to the participants (2).
Forebrain BOLD signal response to handgrip exercise.
The statistical parametric maps reflecting the forebrain activity during the handgrip exercise are overlayed onto a structural MRI template image and are presented in Fig. 3. During the 35% trials, a similar forebrain network was activated in both men and women. This included the left MC, somatosensory cortex, supplementary motor area (SMA), the cingulate motor area, bilateral IC, thalamus, left putamen, left caudate, brain stem, right cerebellum, and vermis (Table 2). Significant activation in the right hippocampus was observed in men but not in women. On the contrary, women exhibited activation in the left inferior frontal gyrus, while men did not.

Fig. 3.Brain regions activated and deactivated during the 35% MVC handgrip exercise. A similar forebrain network was activated in men and women during the exercise. It included the left motor cortex (MC), supplementary motor area (SMA), and bilateral insular cortices (IC). In addition, the ventral medial prefrontal cortex (vMPFC) was commonly deactivated in both men and women (P < 0.005 uncorrected). The red blobs denote the activation clusters, whereas the blue blobs denote the deactivation clusters. RH 35%, right hand 35% MVC handgrip exercise.
| Activation | Men | Women | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Coordinates | Z-Score | Coordinates | Z-Score | |||||||||||||
| Location | Side | x | y | z | x | y | z | |||||||||
| Precentral gyrus (MC) | L | −40 | −22 | 52 | 5.47 | −38 | −20 | 50 | 5.31 | |||||||
| Postcentral gyrus (SC) | L | −22 | −30 | 66 | 5.33 | −44 | −24 | 58 | 5.02 | |||||||
| R | 48 | −28 | 40 | 5.02 | 52 | −34 | 44 | 4.73 | ||||||||
| SMA | L | −8 | 6 | 44 | 4.12 | −6 | 8 | 44 | 4.14 | |||||||
| R | 12 | 0 | 62 | 4.38 | 4 | 12 | 40 | 4.41 | ||||||||
| Midcingulate | R | 6 | −16 | 34 | 4.29 | 8 | −4 | 34 | 3.58 | |||||||
| Inferior frontal gyrus | L | −56 | 4 | 12 | 4.64 | |||||||||||
| Insula | L | −36 | 4 | 4 | 4.85 | −32 | 10 | 8 | 3.56 | |||||||
| R | 36 | 6 | 4 | 5.46 | 40 | 12 | 10 | 3.89 | ||||||||
| Thalamus | L | −4 | −14 | 14 | 6.03 | −14 | −22 | 10 | 4.91 | |||||||
| R | 10 | −6 | 12 | 5.78 | 14 | −8 | 12 | 4.78 | ||||||||
| Putamen | L | −30 | −4 | 0 | 4.88 | −20 | 0 | 6 | 2.95 | |||||||
| Hippocampus | R | 28 | −20 | −6 | 4.44 | |||||||||||
| Caudate | L | −18 | −22 | 20 | 6.6 | −18 | −4 | 24 | 3.85 | |||||||
| Brainstem | L | −4 | −20 | −18 | 5.06 | 0 | −20 | −20 | 3.72 | |||||||
| R | 6 | −26 | −12 | 3.48 | ||||||||||||
| Cerebelum | R | 16 | −48 | −28 | 5.65 | 18 | −48 | −28 | 5.61 | |||||||
| Vermis | 6 | −72 | −22 | 5.58 | 6 | −62 | −14 | 4.79 | ||||||||
In addition, men and women showed a comparable deactivation pattern in the forebrain (Fig. 3). As summarized in Table 3, the vMPFC, posterior cingulate cortex, right inferior insula, bilateral lingual gyri, and left calcarine were commonly deactivated in both men and women. Nevertheless, deactivation in the bilateral temporal pole was observed in men only, whereas deactivations in the left angular gyrus, left superior medial frontal gyrus, right fusiform gyrus, bilateral superior temporal gyri, bilateral superior frontal gyri, parahippocampus, and dorsal ACC (dACC) were observed in women only.
| Deactivation | Men | Women | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Coordinates | Z-Score | Coordinates | Z-Score | |||||||||||||
| Location | Side | x | y | z | x | y | z | |||||||||
| vMPFC | 4 | 42 | −6 | 4.18 | 0 | 44 | −6 | 4.61 | ||||||||
| PCC | −2 | −50 | 30 | 3.16 | −4 | −46 | 30 | 4.6 | ||||||||
| dACC | −2 | 24 | 18 | 2.89 | ||||||||||||
| Superior posterior insula | R | 40 | −18 | 10 | 2.85 | |||||||||||
| Inferior insula | R | 32 | 12 | −18 | 3.22 | 40 | 10 | −18 | 3.75 | |||||||
| Temporal pole | L | −38 | 8 | −26 | 2.97 | |||||||||||
| R | 48 | 8 | −28 | 3.03 | ||||||||||||
| Superior temporal gyrus | L | −46 | −10 | −8 | 3.11 | |||||||||||
| R | 50 | −6 | 12 | 3.3 | ||||||||||||
| Angular gyrus | L | −44 | −58 | 28 | 4.79 | |||||||||||
| Superior medial frontal gyrus | L | −4 | 62 | 18 | 3.25 | |||||||||||
| Superior frontal gyrus | L | −22 | 30 | 44 | 3.44 | |||||||||||
| R | 22 | 24 | 44 | 4.31 | ||||||||||||
| Parahippocampus | L | −30 | −42 | −10 | 3.81 | |||||||||||
| R | 34 | −20 | −20 | 3.17 | ||||||||||||
| Fusiform gyrus | R | 26 | −40 | −12 | 4.11 | |||||||||||
| Lingual gyrus | L | −12 | −48 | −6 | 3.9 | −10 | −60 | 2 | 4.52 | |||||||
| R | 14 | −50 | 0 | 4.29 | 14 | −54 | 4 | 4.8 | ||||||||
| Calcarine | L | −12 | −62 | 14 | 4.01 | −10 | −62 | 18 | 5.06 | |||||||
The averaged time course of the BOLD responses in two a priori regions, namely, the left MC and the vMPFC, are plotted in Fig. 4. The 35% trials elicited stronger activation than the 5% trials in the left MC. However, the BOLD responses at MC were not different between men and women. Moreover, stronger deactivation was observed in the vMPFC during the 35% trails than the 5% trials. In both men and women, vMPFC activity decreased progressively during the exercise period and reached minimum at the end of the exercise. These were inversely correlated with the HR response during the exercise period (Fig. 1). In addition, women elicited slightly stronger vMPFC deactivation than men during the 35% trials. However, the differences between the two groups had a P value >0.005 (uncorrected) and, therefore, did not reach our predefined criterion for statistical significance.

Fig. 4.Time course of the blood oxygenation level-dependent (BOLD) response in the left primary MC and vMPFC. Top: average BOLD response of the MC for both groups during the 5 and 35% exercise trials. The 35% trials elicited stronger activation in the left MC than the 5% trials. However, the BOLD response of the MC was not different between men and women. Bottom: average BOLD response of the vMPFC. Stronger deactivation was elicited in the vMPFC during the 35% trails than the 5% trials. In both men and women, vMPFC activity decreased progressively during the exercise period and reached minimum at the end of the exercise. However, women elicited slightly stronger vMPFC deactivation than men during the 35% trials. The red and blue regions represent the activation and deactivation clusters, respectively. The BOLD time course data were averaged over the six repeated blocks, and scaled to a zero baseline. The shaded area represents the exercise period.
To further examine the forebrain activation associated with sex differences in HR, direct subtraction analysis was performed to expose the extent of the differences in the BOLD response between men and women in the 35% trials. Men demonstrated stronger activation in the right IC, left thalamus, left putamen, right hippocampus, and brain stem, whereas women exhibited stronger deactivation in the right superior frontal gyrus, left angular gyrus, dACC, and left parahippocampus (P < 0.005 uncorrected; Fig. 5, Table 4).

Fig. 5.Sex differences in BOLD signal during 35% MVC handgrip exercise. A: brain stem; B: right (R) insula; C: dorsal anterior cingulate cortex (dACC). Men demonstrated stronger activation in the right IC, left thalamus, left putamen, right hippocampus, and brain stem, whereas women exhibited stronger deactivation in the right superior frontal gyrus, left angular gyrus, dACC, and left parahippocampus (P < 0.005 uncorrected). Right: signal intensity change at the brain stem, right insula, and dACC. The BOLD signal time course was extracted from the voxels of peak activation in the subtraction analysis and was averaged over the six repeated trials for all of the subjects. The percentage of signal change is relative to the average of the 30-s preexercise period (white area).
| Men > Women | Coordinates | Z-Score | ||||||
|---|---|---|---|---|---|---|---|---|
| Location | Side | x | y | z | ||||
| Superior frontal gyrus | R | 20 | 26 | 44 | 4.57 | |||
| Insula | R | 40 | 10 | −8 | 3.6 | |||
| Angular gyrus | L | −38 | −58 | 32 | 3.05 | |||
| Dorsal ACC | L | −4 | 24 | 18 | 3.02 | |||
| Thalamus | L | −4 | −8 | 14 | 3.92 | |||
| Putamen | L | −32 | −4 | −2 | 3.19 | |||
| Hippocampus | R | 28 | −20 | −6 | 3.4 | |||
| Parahippocampus | L | −34 | −40 | −6 | 3.27 | |||
| Caudate | L | −22 | −18 | 20 | 4.21 | |||
| Brain stem | −8 | −20 | −16 | 2.94 | ||||
DISCUSSION
The present study examined the sex differences in cardiovascular responses and the associated forebrain activity during a brief isometric handgrip contraction, which elevated cardiovascular arousal but not peripheral sympathetic activity. A robust tachycardia and a small pressor response occurred during the 35% MVC handgrip exercise in both men and women. However, greater HR and blood pressure responses developed in men than in women. These sex differences in cardiovascular responses were associated with stronger neural activity in the right IC, left thalamus, left putamen, right hippocampus, and brain stem in men. In addition, stronger deactivation in the right superior frontal gyrus, left angular gyrus, dACC, and left parahippocampus occurred in women. In contrast to the hypothesis, subtraction analysis indicated that the deactivation of vMPFC region, although visually greater in women than men in Fig. 4, was not different in the two groups, despite notably different HR responses. An unexpected observation was the deactivation of the dACC in women but not men. Our findings suggest that the early cardiovagal response to brief isometric exercise is greater in men than women. Although the hypothesis that a greater HR change in men would be linked with stronger vMPFC deactivation was not supported, there were distinct sex differences in other regions of the cortical autonomic network, particularly the dACC and insular cortices.
Sex differences in the physiological response to brief isometric exercise.
The HR response to a brief isometric handgrip task is largely due to cardiovagal withdrawal. The observations that MSNA did not increase in either group in this task supports this conclusion. From this perspective, the sex-dependent HR responses suggest that men either display greater vagal withdrawal, or have a greater end-organ response to a given level of vagal withdrawal to this type of physiological stress, than women. The finding of greater physiological responses in men during handgrip exercise is consistent with earlier observations during orthostatic stress (45). Nonetheless, the current findings may contrast with the earlier conclusions that cardiovagal control of HR is more pronounced in women (31, 32, 46). However, the differences may relate to task-specific responses and/or to the fact that we measured the HR response rather than spectral power at particular frequencies. Such measures were not possible during this study due to the short duration of handgrip task and the non-steady-state HR responses. Furthermore, changes in sympathetic activation patterns at the heart may have occurred in this handgrip protocol that were not detected in the peroneal recordings.
As with previous studies (16, 45), baseline MSNA levels are lower in women than men at baseline. Ettinger et al. (16) reported that steady-state levels of MSNA after 1–2 min of fatiguing handgrip were different between male and female participants; however, sex differences at the onset of exercise have not been studied. Thus an intriguing finding of this study, although based on a small number of subjects, is that MSNA may be reduced in this early period in women but not men. This observation was not hypothesized and requires further experimentation. Nonetheless, this MSNA response raises the intriguing possibility that the greater cardiovascular responses in men could be driven by variations in the suppression of cardiac vagal or sympathetic control between the sexes.
Sex differences in the forebrain neural responses to a brief isometric exercise.
In this study, activation of cortical areas that were associated with motor execution, including the left MC and SMA, occurred commonly in men and women (Fig. 3). The stereotaxic coordinates of the activation peaks were highly overlapped between the two groups (Table 2), indicating that men and women activated a similar neural network in executing the motor task. In the MC, the magnitude of BOLD responses was identical between men and women (Fig. 4), suggesting that the two groups exerted comparable motor effort in completing the handgrip exercise. The consistent Borg scale exertion rating between men and women supports this proposition.
In addition to the motor control network, BOLD responses were observed in a number of additional cortical areas that are associated with autonomic modulation. These included the vMPFC, dACC, and IC, which are parts of the central autonomic network. Importantly, significant sex differences in the neural responses were present in these areas.
The vMPFC is of particular interest in the present study, as it has been associated with cardiovagal control (9, 36). As noted in the previous study (50), the time course of vMPFC response closely resembled the HR response during the exercise period (Fig. 1). The hypothesis that vMPFC deactivation is related directly to cardiovagal control, and the knowledge that the cardiovascular responses to the current isometric handgrip exercise are predominantly mediated by vagal withdrawal, suggest that the greater HR response in men should be associated with stronger vMPFC deactivation. However, the greater HR response, despite similar (or even reduced) reduction in vMPFC activation in men, stands in contrast to this hypothesis. Therefore, while the vMPFC may be important in overall cardiovagal regulation, this paradoxical observation suggested that it is not the sole factor in regulating cardiovascular responses. Rather, subcortical modulation of the cortical signal, and/or other peripheral factors that affect the autonomic expression at the heart, may also contribute.
For instance, previous studies showed that cardiovascular and ventilatory responses to central command were attenuated by estrogen (22). In premenopausal women, the MSNA response to an extended period of static exercise was affected by the level of estrogen; sympathoexcitation was attenuated during follicular phase in which the level of estrogen was higher (17). In male rats, injection of estrogen at the nucleus tractus solitarius attenuated sympathetic outflow and enhanced parasympathetic tone (41). Also, Del Rio et al. (14) reported that cardiovascular responses to mental stress were attenuated in men with percutaneous estrogen administration. Accordingly, in the present study, the higher level of estrogen in women may have dampened the reduction of cardiovagal tone to the heart and resulted in relatively smaller cardiovascular responses in women.
In addition to the ventral portion of the MPFC, previous neuroimaging studies investigating cardiovascular control reported that, in humans, tachycardia was often accompanied by increased activity at the dorsal part of the MPFC, in particular the dACC (11, 12, 27, 33, 49). Most of these studies used maneuvers that elevated HR and sympathetic outflow. Thus the activity of dACC was associated with modulating sympathetic activity (9). In the present study, dACC activation did not occur during the handgrip task (Table 2, Fig. 3); it was indeed deactivated in women but not in men during the 35% MVC trials (Table 3, Fig. 5). Such sex differences in the dACC response matched well with the MSNA data, where, during the handgrip exercise, women demonstrated a trend toward suppression of peripheral sympathetic activity, whereas men did not (Fig. 2). It has been suggested that the suppression of MSNA is induced by central command during mild isometric muscle contraction (34). According to Mark et al. (34), MSNA burst rate falls during the first minute of 30% MVC exercise, whereas an increase of MSNA occurs during involuntary contraction induced by muscle stimulation. These observations supported the preposition that suppression of MSNA during the first minute of isometric exercise is induced by active top-down inhibitory signals associated with central command (34). Nonetheless, Mark et al. included both men and women in their study, and they did not distinguish whether there were sex differences in the MSNA response. Thus our findings not only support the previous conclusion that the dACC is involved in modulating sympathetic activity, they also provide new evidence that women appear to actively suppress their sympathetic outflow at the onset of isometric exercise. Taken together, the smaller tachycardia and pressor response in women may be affected by both their higher level of estrogen that restrained the reduction of vagal tone to the heart and their stronger dACC deactivation that suppressed the sympathetic activity during the exercise. Whether the reduction in MSNA reflects concurrent reductions in baseline sympathetic outflow to the heart could not be assessed in this retrospective analysis.
In addition to the MPFC, the IC have been associated with autonomic regulation in animals (3, 25, 39, 40, 51) and humans (11, 12, 27, 39, 49). Anatomically, the IC, which are part of the orbital prefrontal network (38), receive visceral afferent input from internal organs and provide a topographic representation of internal states (4). Our fMRI data showed bilateral IC activation during the handgrip exercise in both men and women (Fig. 3). Relative to women, men elicited stronger IC activation (Fig. 5) that corresponded to their greater cardiovascular responses. This observation is in line with the findings of our laboratory's previous study examining sex differences in baroreflex response in which men elicited greater cardiovascular responses and stronger IC activation relative to women (29).
Currently, it is uncertain whether the increase in IC activation represents somatosensory stimulation within the heart and/or blood vessels as HR and MAP changed. Previous experiments showed that the activities of IC neurons correlated with cardiac rhythm and were responsive to muscular contraction (25). The work of Cechetto and Saper (4) suggests that the IC represents viscerotopic sensory aspects of cardiovascular arousal. However, additional clinical and animal studies have revealed that cardiovascular responses can be elicited by electrical stimulation at the IC (3, 39, 40, 51). Also, ischemic stroke in the IC increases the risk of cardiac arrhythmias (5-7). These studies indicate that IC is not merely reflecting the internal states, but also has an active role in feedback modulation of cardiovascular responses through its connection with the medial prefrontal network. Therefore, this retrospective analysis has exposed the possibility that the stronger vMPFC deactivation in women may be associated with a lower level of IC activation.
Limitations.
An assumption in this study was that the rapid HR response to the exercise challenge was related to cardiovagal control. Although we were not able to collect a direct measure of vagal activity during the fMRI experiment, the MSNA data and evidence from previous pharmacological and microneurographic studies (19, 24, 34, 37, 47, 49) provide confidence that the brief isometric handgrip exercise elevated HR mainly by attenuating vagal activity to the heart.
In this study, the fMRI experiment adopted a blocked design paradigm to examine the sex differences in cortical responses over a 30-s period. A possible limitation of the blocked design is that it may have limited capacity in detecting subtle differences over a very brief period. For example, men and women might elicit distinct BOLD responses at the onset of exercise but similar responses at the end of the exercise. This does not appear to be the case, as the BOLD signal time course plots in Fig. 5 provide a rough illustration that the men demonstrated an overall stronger response in the brain stem, right insula, and dACC. It is cautioned, however, that these time course plots may be biased to the subtraction results, because the data were extracted from the peak voxel in the subtraction results.
In addition, the scheduling of tests in the women could not accommodate variations in the menstrual cycle due to challenges in scheduling the multiple testing sessions with fMRI accessibility. Menstrual phase did not appear to affect the general pattern of responses during tilt or cold pressor experiments in a previous study (45). However, menstrual phase did appear to affect the magnitude of cardiovascular responses to fatiguing handgrip exercise, where attenuated responses were observed during the follicular phase (17). Whether this limitation applies to the exercise onset phase before fatigue will require additional study. Nevertheless, the major influence of not testing women at the same phase of menstrual cycle would be an increase of variability in the female groups and hence a reduction of statistical significance. Thus it is not expected that scheduling the tests based on menstrual phase would affect the overall conclusion that one's sex, and possibly the levels of sex hormones, may affect not only the magnitude of the peripheral response, but that of the cortical activation patterns as well.
In summary, our physiological and neuroimaging data revealed that men and women had different cardiovascular and forebrain neural responses to a brief isometric handgrip contraction. The differences in cortical responses were not necessarily related to different anatomical regions in the two groups, but rather to the variations of activation/deactivation patterns within the same and expected regions of the cortical autonomic network. Compared with men, women had smaller HR and blood pressure responses as well as some suppression of MSNA that, in turn, was related to dACC deactivation in women. The greater cardiovascular response in men was coincident with stronger IC activation. However, paradoxical associations between HR and vMPFC deactivation may reflect the influence of estrogen on cardiac behavior and the active suppression of sympathetic outflow in women. Thus these data provide new evidence for sex-based differences in the cortical and peripheral autonomic regulation of cardiovascular function during isometric exercise.
GRANTS
This research was supported by research grants from the Heart and Stroke Foundation of Ontario (NA no. 5020 and T no. 5342), The Ontario March of Dimes, and the Natural Sciences and Engineering Research Council of Canada.
FOOTNOTES
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.
REFERENCES
- 1 Barnabi F, Cechetto DF. Neurotransmitters in the thalamus relaying visceral input to the insular cortex in the rat. Am J Physiol Regul Integr Comp Physiol 281: R1665–R1674, 2001.
Link | ISI | Google Scholar - 2 Borg GA. Psychophysical bases of perceived exertion. Med Sci Sports Exerc 14: 377–381, 1982.
Crossref | PubMed | ISI | Google Scholar - 3 Cechetto DF, Chen SJ. Subcortical sites mediating sympathetic responses from insular cortex in rats. Am J Physiol Regul Integr Comp Physiol 258: R245–R255, 1990.
Link | ISI | Google Scholar - 4 Cechetto DF, Saper CB. Evidence for a viscerotopic sensory representation in the cortex and thalamus in the rat. J Comp Neurol 262: 27–45, 1987.
Crossref | PubMed | ISI | Google Scholar - 5 Cheung RT, Hachinski V. The insula and cerebrogenic sudden death. Arch Neurol 57: 1685–1688, 2000.
Crossref | Google Scholar - 6 Cheung RT, Hachinski V. Cardiac effects of stroke. Curr Treat Options Cardiovasc Med 6: 199–207, 2004.
Crossref | Google Scholar - 7 Colivicchi F, Bassi A, Santini M, Caltagirone C. Cardiac autonomic derangement and arrhythmias in right-sided stroke with insular involvement. Stroke 35: 2094–2098, 2004.
Crossref | PubMed | ISI | Google Scholar - 8 Convertino VA. Gender differences in autonomic functions associated with blood pressure regulation. Am J Physiol Regul Integr Comp Physiol 275: R1909–R1920, 1998.
Link | ISI | Google Scholar - 9 Critchley HD. The human cortex responds to an interoceptive challenge. Proc Natl Acad Sci USA 101: 6333–6334, 2004.
Crossref | PubMed | ISI | Google Scholar - 10 Critchley HD. Neural mechanisms of autonomic, affective, and cognitive integration. J Comp Neurol 493: 154–166, 2005.
Crossref | PubMed | ISI | Google Scholar - 11 Critchley HD, Corfield DR, Chandler MP, Mathias CJ, Dolan RJ. Cerebral correlates of autonomic cardiovascular arousal: a functional neuroimaging investigation in humans. J Physiol 523: 259–270, 2000.
Crossref | PubMed | ISI | Google Scholar - 12 Critchley HD, Mathias CJ, Josephs O, O'Doherty J, Zanini S, Dewar BK, Cipolotti L, Shallice T, Dolan RJ. Human cingulate cortex and autonomic control: converging neuroimaging and clinical evidence. Brain 126: 2139–2152, 2003.
Crossref | PubMed | ISI | Google Scholar - 13 de Oliveira CV, Rosas-Arellano MP, Solano-Flores LP, Babic T, Li Z, Ciriello J. Estrogen alters the bradycardia response to hypocretin-1 in the nucleus tractus solitarius of the ovariectomized female. Brain Res 978: 14–23, 2003.
Crossref | ISI | Google Scholar - 14 Del Rio G, Velardo A, Zizzo G, Avogaro A, Cipolli C, Della CL, Marrama P, MacDonald IA. Effect of estradiol on the sympathoadrenal response to mental stress in normal men. J Clin Endocrinol Metab 79: 836–840, 1994.
PubMed | ISI | Google Scholar - 15 Delius W, Hagbarth KE, Hongell A, Wallin BG. General characteristics of sympathetic activity in human muscle nerves. Acta Physiol Scand 84: 65–81, 1972.
Crossref | PubMed | Google Scholar - 16 Ettinger SM, Silber DH, Collins BG, Gray KS, Sutliff G, Whisler SK, McClain JM, Smith MB, Yang QX, Sinoway LI. Influences of gender on sympathetic nerve responses to static exercise. J Appl Physiol 80: 245–251, 1996.
Link | ISI | Google Scholar - 17 Ettinger SM, Silber DH, Gray KS, Smith MB, Yang QX, Kunselman AR, Sinoway LI. Effects of the ovarian cycle on sympathetic neural outflow during static exercise. J Appl Physiol 85: 2075–2081, 1998.
Link | ISI | Google Scholar - 18 Evans JM, Ziegler MG, Patwardhan AR, Ott JB, Kim CS, Leonelli FM, Knapp CF. Gender differences in autonomic cardiovascular regulation: spectral, hormonal, and hemodynamic indexes. J Appl Physiol 91: 2611–2618, 2001.
Link | ISI | Google Scholar - 19 Fagraeus L, Linnarsson D. Autonomic origin of heart rate fluctuations at the onset of muscular exercise. J Appl Physiol 40: 679–682, 1976.
Link | ISI | Google Scholar - 20 Friston KJ, Holmes AP, Worsley KJ, Poline JB, Frith CD, Frackowiak RS. Statistical parametric maps in functional imaging: a general linear approach. Hum Brain Mapp 2: 189–210, 1995.
Crossref | Google Scholar - 21 Hagbarth KE, Vallbo AB. Pulse and respiratory grouping of sympathetic impulses in human muscle-nerves. Acta Physiol Scand 74: 96–108, 1968.
Crossref | PubMed | Google Scholar - 22 Hayes SG, Moya Del Pino NB, Kaufman MP. Estrogen attenuates the cardiovascular and ventilatory responses to central command in cats. J Appl Physiol 92: 1635–1641, 2002.
Link | ISI | Google Scholar - 23 Henderson LA, Richard CA, Macey PM, Runquist ML, Yu PL, Galons JP, Harper RM. Functional magnetic resonance signal changes in neural structures to baroreceptor reflex activation. J Appl Physiol 96: 693–703, 2004.
Link | ISI | Google Scholar - 24 Hollander AP, Bouman LN. Cardiac acceleration in man elicited by a muscle-heart reflex. J Appl Physiol 38: 272–278, 1975.
Link | ISI | Google Scholar - 25 Ichiyama RM, Waldrop TG, Iwamoto GA. Neurons in and near insular cortex are responsive to muscular contraction and have sympathetic and/or cardiac-related discharge. Brain Res 1008: 273–277, 2004.
Crossref | ISI | Google Scholar - 26 Jackson DN, Milne KJ, Noble EG, Shoemaker JK. Gender-modulated endogenous baseline neuropeptide Y Y1-receptor activation in the hindlimb of Sprague-Dawley rats. J Physiol 562: 285–294, 2005.
Crossref | PubMed | ISI | Google Scholar - 27 Kimmerly DS, O'Leary DD, Menon RS, Gati JS, Shoemaker JK. Cortical regions associated with autonomic cardiovascular regulation during lower body negative pressure in humans. J Physiol 569: 331–345, 2005.
Crossref | PubMed | ISI | Google Scholar - 28 Kimmerly DS, O'Leary DD, Shoemaker JK. Test-retest repeatability of muscle sympathetic nerve activity: influence of data analysis and head-up tilt. Auton Neurosci 114: 61–71, 2004.
Crossref | PubMed | ISI | Google Scholar - 29 Kimmerly DS, Wong S, Menon R, Shoemaker JK. Forebrain neural patterns associated with sex differences in autonomic cardiovascular function during baroreceptor unloading. Am J Physiol Regul Integr Comp Physiol 292: R715–R722, 2007.
Link | ISI | Google Scholar - 30 Klassen LM, Menon RS. Robust automated shimming technique using arbitrary mapping acquisition parameters (RASTAMAP). Magn Reson Med 51: 881–887, 2004.
Crossref | PubMed | ISI | Google Scholar - 31 Kuch B, Hense HW, Sinnreich R, Kark JD, von EA, Sapoznikov D, Bolte HD. Determinants of short-period heart rate variability in the general population. Cardiology 95: 131–138, 2001.
Crossref | PubMed | ISI | Google Scholar - 32 Kuo TB, Lin T, Yang CC, Li CL, Chen CF, Chou P. Effect of aging on gender differences in neural control of heart rate. Am J Physiol Heart Circ Physiol 277: H2233–H2239, 1999.
Link | ISI | Google Scholar - 33 Macefield VG, Gandevia SC, Henderson LA. Neural sites involved in the sustained increase in muscle sympathetic nerve activity induced by inspiratory capacity apnea: a fMRI study. J Appl Physiol 100: 266–273, 2006.
Link | ISI | Google Scholar - 34 Mark AL, Victor RG, Nerhed C, Wallin BG. Microneurographic studies of the mechanisms of sympathetic nerve responses to static exercise in humans. Circ Res 57: 461–469, 1985.
Crossref | PubMed | ISI | Google Scholar - 35 Matsukawa T, Sugiyama Y, Watanabe T, Kobayashi F, Mano T. Gender difference in age-related changes in muscle sympathetic nerve activity in healthy subjects. Am J Physiol Regul Integr Comp Physiol 275: R1600–R1604, 1998.
Link | ISI | Google Scholar - 36 Matthews SC, Paulus MP, Simmons AN, Nelesen RA, Dimsdale JE. Functional subdivisions within anterior cingulate cortex and their relationship to autonomic nervous system function. Neuroimage 22: 1151–1156, 2004.
Crossref | PubMed | ISI | Google Scholar - 37 Mitchell JH, Reeves DR Jr, Rogers HB, Secher NH, Victor RG. Autonomic blockade and cardiovascular responses to static exercise in partially curarized man. J Physiol 413: 433–445, 1989.
Crossref | PubMed | ISI | Google Scholar - 38 Ongur D, Price JL. The organization of networks within the orbital and medial prefrontal cortex of rats, monkeys and humans. Cereb Cortex 10: 206–219, 2000.
Crossref | PubMed | ISI | Google Scholar - 39 Oppenheimer SM, Gelb A, Girvin JP, Hachinski VC. Cardiovascular effects of human insular cortex stimulation. Neurology 42: 1727–1732, 1992.
Crossref | PubMed | ISI | Google Scholar - 40 Ruggiero DA, Mraovitch S, Granata AR, Anwar M, Reis DJ. A role of insular cortex in cardiovascular function. J Comp Neurol 257: 189–207, 1987.
Crossref | PubMed | ISI | Google Scholar - 41 Saleh MC, Connell BJ, Saleh TM. Medullary and intrathecal injections of 17beta-estradiol in male rats. Brain Res 867: 200–209, 2000.
Crossref | PubMed | ISI | Google Scholar - 42 Saper CB. The central autonomic nervous system: conscious visceral perception and autonomic pattern generation. Annu Rev Neurosci 25: 433–469, 2002.
Crossref | PubMed | ISI | Google Scholar - 43 Schmitt PM, Kaufman MP. Estrogen attenuates the exercise pressor reflex in female cats. J Appl Physiol 95: 1418–1424, 2003.
Link | ISI | Google Scholar - 44 Schmitt PM, Kaufman MP. High concentrations of 17beta-estradiol attenuate the exercise pressor reflex in male cats. J Appl Physiol 94: 1431–1436, 2003.
Link | ISI | Google Scholar - 45 Shoemaker JK, Hogeman CS, Khan M, Kimmerly DS, Sinoway LI. Gender affects sympathetic and hemodynamic response to postural stress. Am J Physiol Heart Circ Physiol 281: H2028–H2035, 2001.
Link | ISI | Google Scholar - 46 Sinnreich R, Kark JD, Friedlander Y, Sapoznikov D, Luria MH. Five minute recordings of heart rate variability for population studies: repeatability and age-sex characteristics. Heart 80: 156–162, 1998.
Crossref | PubMed | ISI | Google Scholar - 47 Victor RG, Seals DR, Mark AL. Differential control of heart rate and sympathetic nerve activity during dynamic exercise. Insight from intraneural recordings in humans. J Clin Invest 79: 508–516, 1987.
Crossref | PubMed | ISI | Google Scholar - 48 Williamson JW, Fadel PJ, Mitchell JH. New insights into central cardiovascular control during exercise in humans: a central command update. Exp Physiol 91: 51–58, 2006.
Crossref | PubMed | ISI | Google Scholar - 49 Williamson JW, McColl R, Mathews D. Evidence for central command activation of the human insular cortex during exercise. J Appl Physiol 94: 1726–1734, 2003.
Link | ISI | Google Scholar - 50 Wong SW, Massé N, Kimmerly DS, Menon RS, Shoemaker JK. Ventral medial prefrontal cortex and cardiovagal control in conscious humans. Neuroimage 35: 698–708, 2007.
Crossref | ISI | Google Scholar - 51 Yasui Y, Breder CD, Saper CB, Cechetto DF. Autonomic responses and efferent pathways from the insular cortex in the rat. J Comp Neurol 303: 355–374, 1991.
Crossref | PubMed | ISI | Google Scholar

