# Altered Relationship Between Heart Rate Variability and fMRI-Based Functional Connectivity in People With Epilepsy

Michalis Kassinopoulos, Ronald Harper, Maxime Guye, Louis Lemieux, Beate Diehl

► To cite this version:

Michalis Kassinopoulos, Ronald Harper, Maxime Guye, Louis Lemieux, Beate Diehl. Altered Relationship Between Heart Rate Variability and fMRI-Based Functional Connectivity in People With Epilepsy. Frontiers in Neurology, 2021, 12, 10.3389/fneur.2021.671890. hal-03501141

HAL Id: hal-03501141

https://amu.hal.science/hal-03501141

Submitted on 25 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Altered Relationship Between Heart Rate Variability and fMRI-Based Functional Connectivity in People With Epilepsy

Michalis Kassinopoulos¹,², Ronald M. Harper³,⁴, Maxime Guye⁵,⁶, Louis Lemieux¹,² and Beate Diehl¹,²*

¹Department of Clinical and Experimental Epilepsy, Institute of Neurology, University College London, London, United Kingdom, ²Epilepsy Society, Buckinghamshire, United Kingdom, ³Brain Research Institute, University of California, Los Angeles, Los Angeles, CA, United States, ⁴Department of Neurobiology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, United States, ⁵Aix Marseille Univ, CNRS, CRMBM, Marseille, France, ⁶APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France

**Background:** Disruptions in central autonomic processes in people with epilepsy have been studied through evaluation of heart rate variability (HRV). Decreased HRV appears in epilepsy compared to healthy controls, suggesting a shift in autonomic balance toward sympathetic dominance; recent studies have associated HRV changes with seizure severity and outcome of interventions. However, the processes underlying these autonomic changes remain unclear. We examined the nature of these changes by assessing alterations in whole-brain functional connectivity, and relating those alterations to HRV.

**Methods:** We examined regional brain activity and functional organization in 28 drug-resistant epilepsy patients and 16 healthy controls using resting-state functional magnetic resonance imaging (fMRI). We employed an HRV state-dependent functional connectivity (FC) framework with low and high HRV states derived from the following four cardiac-related variables: 1. RR interval, 2. root mean square of successive differences (RMSSD), 4. low-frequency HRV (0.04–0.15 Hz; LF-HRV) and high-frequency HRV (0.15–0.40 Hz; HF-HRV). The effect of group (epilepsy vs. controls), HRV state (low vs. high) and the interactions of group and state were assessed using a mixed analysis of variance (ANOVA). We assessed FC within and between 7 large-scale functional networks consisting of cortical regions and 4 subcortical networks, the amygdala, hippocampus, basal ganglia and thalamus networks.

**Results:** Consistent with previous studies, decreased RR interval (increased heart rate) and decreased HF-HRV appeared in people with epilepsy compared to healthy controls. For both groups, fluctuations in heart rate were positively correlated with BOLD activity in bilateral thalamus and regions of the cerebellum, and negatively correlated with BOLD activity in the insula, putamen, superior temporal gyrus and inferior frontal gyrus. Connectivity strength in patients between right thalamus and ventral attention network (mainly insula) increased in the high LF-HRV state compared to low LF-HRV; the opposite trend appeared in healthy controls. A similar pattern emerged for connectivity between the thalamus and basal ganglia.

**OPEN ACCESS**

**Edited by:**
Maria Centeno,
University College London,
United Kingdom

**Reviewed by:**
Piergiorgio Salvan,
University of Oxford, United Kingdom
Friederike Moeller,
Great Ormond Street Hospital for
Children NHS Foundation Trust,
United Kingdom

**\*Correspondence:**
Beate Diehl
b.diehl@ucl.ac.uk

**Specialty section:**
This article was submitted to
Applied Neuroimaging,
a section of the journal
Frontiers in Neurology

**Received:** 24 February 2021
**Accepted:** 18 May 2021
**Published:** 10 June 2021

**Citation:**
Kassinopoulos M, Harper RM, Guye M, Lemieux L and Diehl B (2021) Altered Relationship Between Heart Rate Variability and fMRI-Based Functional Connectivity in People With Epilepsy. Front. Neurol. 12:671890. doi: 10.3389/fneur.2021.671890
**Conclusion:** The findings suggest that resting connectivity patterns between the thalamus and other structures underlying HRV expression are modified in people with drug-resistant epilepsy compared to healthy controls.

**Keywords:** state-dependent functional connectivity, sympathovagal balance, SUDEP, thalamic connectivity, ventral attention network, insula cortex

# INTRODUCTION

Heart rate varies on a moment-to-moment basis in response to changing physiological demands, and is regulated by sympathetic and parasympathetic components of the autonomic nervous system (ANS). Evaluation of the momentary changes in heart rate variability (HRV) can provide insights into the interplay of central mechanisms controlling sympathetic and parasympathetic (vagal) activity (1, 2). A shift toward parasympathetic dominance is typically accompanied by heart rate declines and increased HRV; whereas, increased sympathetic dominance is typically associated with an accelerated heart rate and decreased HRV [although deviations from this generality occur (3)]. Considerable evidence exists that HRV provides an indication of sympathovagal balance and can be useful as a marker for certain cardiovascular diseases (4), mortality, and sudden death (5).

Epilepsy is accompanied by significantly different patterns of HRV (6). A lower interictal HRV is often reported in drug-resistant epilepsy, suggesting a shift toward sympathetic predominance (6–9). In addition, a link between peri-ictal HRV and major motor seizure severity (10) has been outlined, as well as an indication of seizure reduction following vagal nerve stimulation in patients with drug-resistant epilepsy (11). HRV determination of low parasympathetic activity and increased risk of sudden unexpected death in epilepsy SUDEP has been described (12, 13), as well as altered circadian rhythms of HRV in epilepsy (14, 15); the latter finding may explain the larger number of night-time SUDEP cases (16). However, a poor understanding of the mechanisms underlying expression of cardiac functions in epilepsy hampers interpretation of alterations in brain regulatory sites controlling HRV and the potential to gain insights into dysfunctions within those processes.

Functional magnetic resonance imaging (fMRI), a non-invasive tool for probing brain activity and functional connectivity (FC), has been used to study the neural substrates of autonomic regulation (17–20). Initial studies primarily relied on tasks to excite the ANS (21–26), while subsequent studies have used resting-state fMRI (27–29), which has a benefit of not being confounded by task-related changes in local brain activity and FC. Differences in HRV across participants as well as fluctuations in HRV within-individuals have been related to spontaneous regional blood-oxygen-level-dependent (BOLD) fluctuations and connectivity between distinct regions (27–29). Regions found in fMRI studies to be associated with autonomic regulation, such as the anterior cingulate (ACC), medial prefrontal (mPFC) and insular cortices, form part of the central autonomic network (CAN) described in preclinical

studies, a system of brain structures involved in ANS functions (30, 31).

Functional connectivity measures between brain sites are altered in people with epilepsy (32–36); however, it is unclear whether these alterations are linked to impaired cardiac regulation. Here, we investigated alterations in brain functional organization in relation to cardiac rhythms in people with epilepsy. We employed an HRV state-dependent FC framework with two levels of variability states estimated from electrocardiogram (ECG) recorded during resting-state fMRI. Given the association between HRV measures and time-varying FC reported in the literature (27), a state-dependent FC framework informed by concurrent cardiac recordings appeared more suitable for studying cardiac dysfunction than static FC approaches that do not utilize physiological recordings. Moreover, we examined whole-brain FC in a data-driven manner rather than restricting the analysis to interactions between regions of the CAN, as recent studies suggested that the neural correlates of cardiac regulation are more widespread than initially thought (17, 29, 37).

# MATERIALS AND METHODS

## Subjects

Thirty-two (32) patients with drug-resistant epilepsy were selected from an ongoing investigation into the localization of epileptic activity in the brain using simultaneous EEG-fMRI with ECG (38), with a case ascertainment period between 2005 and 2014. The inclusion criteria were: (1) the availability of a resting-state EEG-fMRI scan; and (2) a high-resolution T1-weighted scan. The exclusion criteria were: (1) large brain lesion or previous neurosurgery [we considered large to be anything greater than a small area of focal cortical dysplasia (FCD) or sclerosis – e.g., tumors, cavernomas] (2) incomplete clinical or imaging data (e.g., abandoned scans). Sixteen (16) healthy controls were also considered with comparable age and sex characteristics; the group demographics and clinical details are shown in **Supplementary Tables 1, 2**. The study was approved by the National Research Ethics Committee (United Kingdom; 04/Q0502/89) and all patients gave written informed consent.

## Simultaneous EEG-fMRI Acquisition

Scanning was performed at the Epilepsy Society (Chalfont St Peter, Buckinghamshire, UK) on a 3.0 Tesla GE (Signa excite HDX) scanner. A 20-min (400 vol) T₂-weighted fMRI scan was collected from each subject except for two patients that were scanned for 10-min instead. The fMRI scan was done using a gradient-echo echo-planar-imaging with the following
characteristics: repetition time (TR) = 3,000 ms, echo time (TE) = 30 ms, flip angle = 90, matrix size = 64 × 64, field of view (FOV) = 24 × 24 cm², slice thickness = 2.4 mm with 0.6 mm gap, 44 slices, and voxel size = 3.75 × 3.75 × 3 mm³. Subjects were instructed to keep their eyes closed, avoid falling asleep, and not think about anything in particular. A T₁-weighted image was also acquired using an FSPGR (fast spoiled gradient recalled echo) sequence, with the following parameters: matrix size = 256 × 256, FOV = 24 × 24 cm², slice thickness = 1.5 mm, 150 slices, and voxel size = 0.94 × 0.94 × 1.5 mm³.

Scalp EEG signals and an ECG signal were simultaneously acquired during fMRI scanning at 5 kHz using a 64 channel MR-compatible EEG system with ring Ag/AgCl electrodes (BrainAmp MR+; Brain Products GmbH, Munich, Germany). The electrodes were placed according to the 10/20 system and referenced to electrode FCz.

## Preprocessing of fMRI Data

As described previously (38), preprocessing of fMRI data was conducted using the Statistical Parametric Mapping software (SPM12, Welcome Trust Centre for Neuroimaging, London, UK, http://www.fil.ion.ucl.ac.uk/spm) (39) in a Matlab environment (R2020a; Mathworks, Natick, Massachusetts, USA). The first five functional volumes were discarded to allow steady-state magnetization to be established, and the remaining volumes were realigned to correct for head movements. The structural image of each subject was co-registered to the mean realigned functional volume and, subsequently, underwent tissue segmentation into gray matter, white matter and cerebrospinal fluid tissue compartments. The functional images as well as the coregistered structural images and tissue compartment masks were spatially normalized to the Montreal Neurological Institute (MNI) reference space using non-linear transformation.

To account for anatomical variability across participants and reduce thermal noise, all individual functional volumes were smoothed using a 5 mm full-width half-maximum (FWHM) Gaussian kernel. Subsequently, the Brainnetome atlas was used to extract mean fMRI time-series from 210 cortical and 36 subcortical parcels (40). The parcel time-series were high-pass filtered at 0.008 Hz to avoid spurious correlations that arise from low-frequency fluctuations (41).

We used the framewise displacement (FD) as defined in Power et al. (42) to identify and exclude subjects with high levels of motion, as motion can obscure neural-related BOLD activity (43, 44) and lead to systematic biases in FC studies (45–47). FD is calculated from the six motion realignment parameters and reflects the extent of motion at each timepoint. Subjects that were characterized by mean FD larger than 0.25 mm were excluded. In addition, for the remaining of the subjects that were considered in the study, timepoints with FD larger than 0.2 mm were disregarded.

Finally, to further mitigate the effects of motion as well as reduce the effects of physiological processes and scanner artifacts, we regressed out the following nuisance regressors from all parcel time-series: the first ten principal components from voxel time-series within the white matter (48), six regressors related to cardiac pulsatility artifacts obtained with the convolution model

proposed in Kassinopoulos and Mitsis (45), and the mean fMRI time-series averaged across all voxels within the gray matter.

## Preprocessing of ECG and Calculation of HRV Measures

The ECG was corrected for gradient artifacts using adaptive template subtraction (49) implemented in BrainVision Analyzer 2 software (Brain Products GmbH, Munich, Germany), and band-pass filtered from 0.5 to 40 Hz. The R-waves were detected using Matlab's function findpeaks with a minimum peak distance varying between 0.5 and 0.9 s depending on the subject's average RR interval (time between successive R-waves).

The RR intervals were used to obtain time-series of the root mean square of successive differences in RR intervals (RMSSD), and the normalized low (0.04–0.15 Hz) and high (0.15–0.40 Hz) frequency components of HRV. The aforementioned three HRV measures were computed in adjacent time windows of 100 s each, and a timestep of 1 s, to probe changes in sympathetic and parasympathetic activity during the 20-min resting-state scan. RMSSD is a time-domain HRV measure that is believed to reflect parasympathetic activity (50), the low-frequency HRV (LF-HRV) is a frequency-domain measure presumably sensitive to both branches of the ANS, and the high-frequency HRV (HF-HRV) is a frequency-domain measure that, similar to RMSSD, reflects parasympathetic activity. To derive the normalized LF- and HF-HRV measures, the time-series of successive differences in RR intervals was uniformly resampled at 10 Hz before estimating the Welch power spectral density. Subsequently, the power within the frequency ranges 0.04–0.15 Hz (LF-HRV) and 0.15–0.40 Hz (HF-HRV) was divided by the power within the range 0.04–0.50 Hz and multiplied by 100%. Apart from the three HRV measures, the moving average of RR intervals was also computed across the time windows. Before the calculation of an HRV measure or mean RR interval within a time window, outliers of RR values, defined as three median absolute deviations (MAD) away from the median, were linearly interpolated. To disentangle fluctuations in HRV from fluctuations in RR interval, the three HRV measures were orthogonalized with respect to fluctuations in RR interval. Furthermore, heart rate was estimated as the inverse of instantaneous RR interval multiplied by 60. The heart rate was preferred over the instantaneous RR interval to facilitate comparison with the activation maps shown in Valenza et al. (29) that link regional BOLD fluctuations to heart rate.

## Relationship Between BOLD Fluctuations and Cardiac Dynamics

The three HRV measures (RMSSD, LF-HRV and HF-HRV), moving-average RR interval and instantaneous heart rate were convolved with the canonical hemodynamic response function (HRF) from SPM12 (39) prior to resampling at the fMRI acquisition rate. The association between the obtained time-series and voxel-wise fMRI time-series within the gray matter was quantified using one-sample t-test on the associated beta parameters derived from the general linear model. Differences in the activation maps between epilepsy patients and healthy controls were assessed with a second-level, mixed-effects analysis
with subjects as the random-effects factor, using a two-sample t-
test on the associated beta parameters. To control for potential
effects, sex, age and levels of head motion (i.e., mean FD) were
treated as covariates. Moreover, to account for false positives, the
statistical maps were thresholded with a voxel-wise threshold of
family-wise error (FWE) rate *p* < 0.05 corrected for multiple
comparisons using the random-field theory (51) and an extent
threshold ≥ 10 voxels.

Whole-Brain HRV State-Dependent
Functional Connectivity

To investigate the effects of HRV state on brain functional
organization, a whole-brain state-dependent FC analysis was
performed on the parcellation of the Brainnetome atlas that
includes 246 parcels covering the neocortex and sub-cortical
regions (40). The connectivity strength between parcel pairs was
determined based on the pairwise Pearson correlation coefficient
of the parcel time-series. Cortical parcels were grouped into
the following seven large-scale functional networks described in
Yeo et al. (52): visual (34 parcels), somatomotor (33 parcels),
dorsal attention (30 parcels), ventral attention (22 parcels), limbic
(26 parcels), frontoparietal (26 parcels) and default mode (36
parcels), using the mapping provided on the Brainnetome atlas'
website (https://atlas.brainnetome.org/), and the subcortical
parcels were grouped into the following four networks: amygdala
(four parcels), hippocampus (four parcels), basal ganglia (12
parcels) and thalamus (16 parcels) (three cortical parcels were
excluded from the analysis as they were not assigned to any of
the networks). To better understand how FC depends on the
state of autonomic activity, we estimated FC in each individual
considering low or high HRV states separately. The low and high
HRV states were defined as the timepoints in a scan at which an
HRV measure (e.g., RMSSD) had values in the lowest and highest
quartile range for that given scan, respectively.

A mixed analysis of variance (ANOVA) was conducted with
the group (epilepsy patients / healthy controls) as a between-
subject factor and HRV state (low/high) as a within-subject
factor, which allowed us to examine the effect of the group
and HRV state on FC within and between networks as well as
their interactions. Potential effects of sex, age and levels of head
motion on FC were regressed out through linear regression at the
group level before conducting the mixed ANOVA. The levels of
head motion for the low and high HRV state were determined
separately considering only the timepoints corresponding to
each state. The connectivity strength between pairs of networks
that was used in the mixed ANOVA was defined as the mean
correlation averaged across all pairs of parcels that belonged to
the two corresponding networks. The HRV state-dependent FC
analysis was performed for the three HRV measures and the
moving-average RR interval. Statistical significance was set at *p* <
0.0125 (i.e., 0.05/4) adjusted for multiple comparison with respect
to pairs of networks using false discovery rate (FDR).

RESULTS

Data from four patients were excluded due to excessive motion
(mean FD > 0.25 mm). The sex and age distributions were

similar between the epilepsy patients (n = 28, mean age of 28.7,
14 women) and healthy controls (n = 16, mean age of 30.6,
seven women); (*p* > 0.48; two-sample permutation test; number
of permutations *q* = 10,000; **Supplementary Table 2**). Motion-
contaminated fMRI volumes (FD > 0.2 mm) were also excluded,
resulting in an average of 353 ± 65 volumes per subject. Based
on a two-sample permutation test (*q* = 10,000) there were no
significant differences in the number of volumes between healthy
controls and epilepsy patients (*p* > 0.10), and the two groups
exhibited similar levels of motion during the fMRI scan (*p* >
0.10; **Supplementary Figure 1**). The consideration of a subset of
volumes for the low and high HRV states (average 88 volumes
per state) did not have any apparent effects on the estimates of
whole-brain FC as compared to the FC matrices obtained from
the entire scan (**Supplementary Figure 2**).

Heart Rate and HRV Measures

Comparisons of cardiac dynamic metrics between patients
and healthy controls were performed using a two-sample
permutation test (q = 10,000) after regressing out potential
effects of sex and age. The mean RR interval during the
20-min scan was significantly lower in patients compared to
healthy controls (950 ± 100 ms vs. 1,100 ± 200 ms; *p* <
0.003; equivalently, the mean heart rate was significantly higher)
(Figure 1). HF-HRV was also lower in epilepsy compared
to controls (71 ± 11% vs. 78 ± 8%; *p* < 0.05) whereas
RMSSD and LF-HRV were similar in the two groups. The
LF-HRV (0.04–0.15 Hz) and HF-HRV (0.15–0.40 Hz) measures
demonstrated a strong negative inter-correlation (r = −0.56;
**Supplementary Figure 3**) whereas the correlations between the
remaining pairs of RR interval and HRV measures were relatively
low (<0.26; **Supplementary Figure 3**).

Relationship Between BOLD Fluctuations
and Cardiac Dynamics

Across all subjects, BOLD signal fluctuations were associated (*p* < 0.05; FWE-corrected) with only one of the cardiac dynamic metrics, namely instantaneous heart rate, in the following regions: positive correlations in the thalamus and several regions of the cerebellum (culmen, declive, uvula, nodulus and inferior semilunar lobule); negative correlations in the bilateral inferior frontal gyrus, orbitofrontal cortex, middle temporal gyrus, precentral gyrus and claustrum, as well as right insula and putamen (Figure 2). A more liberal threshold of *p* < 0.001 uncorrected indicated positive correlations of heart rate in bilateral caudate, and negative correlations in left insula and putamen, as well as bilateral superior temporal gyrus (Supplementary Figure 4). We did not find significant differences in correlations between the two groups.

HRV State-Dependent Functional
Connectivity

When comparing whole-brain connectivity between patients and
healthy controls, for all cardiac dynamic metrics, the strongest
differences were observed in the connectivity strength between
the frontoparietal, limbic and default mode networks, albeit these
did not reach statistical significance (Figure 3 left column). The
connectivity strength of the thalamus with the ventral attention
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Comparison of mean RR interval and heart rate variability (HRV) measures between epilepsy patients and healthy controls. The bottom and top of each box correspond to the 25th and 75th percentiles of the sample distribution, the line in the box corresponds to the median and the crosses indicate outliers, defined as values that are more than 1.5 times the interquartile range away from the edges of the box. The epilepsy patients showed significantly lower values of mean RR interval and HR-HRV than healthy controls. n.s., not significant.

Association of regional BOLD fluctuations with changes in heart rate at the group level

<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Association of regional BOLD fluctuations with changes in heart rate at the group level. Statistical map of one-sample *t*-test considering patients and controls (n = 44), thresholded with a voxel-wise threshold of family-wise error (FWE) rate *p* < 0.05 corrected for multiple comparisons using the random-field theory (51). The unthresholded statistical map is available at https://neurovault.org/collections/9452/.

network and basal ganglia had strong interactions of group and LF-HRV state (*p* < 0.0125, FDR corrected; **Figure 3** right column). The connectivity between the thalamus and ventral attention network demonstrated also strong interactions of group and RMSSD state.

To shed further light on the interactions of group and HRV state in the connectivity of thalamus with the ventral

attention network and basal ganglia, we performed a *post-hoc* analysis on the connectivity between the pairs of parcels of the associated networks with the strongest interactions. Specifically, we investigated the connectivity between the right caudal temporal thalamus and left dorsal granular insula (parcel of the ventral attention network; *t*-test for interactions = 17.5) and the connectivity between the left posterior parietal thalamus
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Mixed analysis of variance (ANOVA) for evaluating the effects of the group (epilepsy vs. controls; 1st column) and HRV state (low vs. high, patients and healthy controls; 2nd column) on whole-brain connectivity, and their interactions (3rd column). The lower triangular matrix corresponds to the *f*-test for pairs of parcels whereas the upper triangular matrix corresponds to the *t*-test for pairs of the eleven networks. The connectivity strength between pairs of networks that was used in the mixed ANOVA was defined as the mean correlation averaged across all pairs of parcels that belonged to the two corresponding networks. Pairs of networks with *p* < 0.0125 after FDR correction are indicated with an asterisk (*). Note that the lower triangular matrices are only shown to provide a qualitative description of the interactions of networks at the parcel level, and that their significance levels are not assessed. We observe that the connectivity strength of the thalamus with the ventral attention network and basal ganglia has strong interactions between the group and the LF-HRV state, which led us to examine thalamic connectivity with the ventral attention network and basal ganglia more carefully.</figcaption></figure>
and left ventromedial putamen (parcel of the basal ganglia; f-test for interactions = 17.7). For both connections of the thalamus, healthy controls exhibited a lower (toward negative values) connectivity strength in the high LF-HRV state compared to the low LF-HRV state, whereas epilepsy patients exhibited the opposite trend (Figure 4).

**Figure 5** shows the degree to which the voxel-level connectivity profile of the right caudal temporal and left posterior parietal thalamus differs between low and high LF-HRV state (red color corresponds to higher correlations in high LF-HRV compared to low LF-HRV, and vice versa for blue color), for both controls and patients as well as the differences between the two groups. When comparing high with low HRV-state, in controls we observe a decrease in the connectivity of the caudal temporal thalamus with the bilateral anterior insula cortex (AIC), the anterior cingulate cortex (ACC), middle frontal gyrus (MFG) and supramarginal gyrus (SMG); whereas, in patients, we observe a small decrease in the connectivity with MFG and an increase with ACC. When examining the effect of LF-HRV state in the connectivity of the posterior parietal thalamus, in controls, a decrease in connectivity with the bilateral AIC, putamen and caudate appears, and in patients, an increase in connectivity with left putamen and AIC as well as bilateral caudate.

DISCUSSION

We examined the association of autonomic cardiac regulation with spontaneous fluctuations in fMRI and whole-brain FC in people with drug-resistant epilepsy, compared to healthy controls. In both groups, heart rate was positively correlated with fMRI signal intensity in bilateral thalamus and regions of the cerebellum, and negatively correlated with lateral regions, including bilateral inferior frontal gyrus, orbitofrontal cortex, middle temporal gyrus and right insula and putamen (Figure 2; Supplementary Figure 4). In addition, fluctuations in RMSSD and LF-HRV exhibited strong associations with changes in FC (Figure 3), despite the absence of correlation with brain activity in individual regions. Importantly, these relations differed between healthy controls and epilepsy patients. In controls, increased levels of RMSSD and LF-HRV were associated with declines in connectivity between thalamus and ventral attention network, whereas in patients, similar HRV changes accompanied increases in connectivity. The different patterns between the two groups were more pronounced for the connectivity between right caudal temporal thalamus and left dorsal granular insula (Figure 4). Note, however, that the interactions between ipsilateral regions [left (right) caudal temporal thalamus with left (right) dorsal granular insula] were also statistically significant, albeit with slightly higher *p*-values (*p* < 0.01). Therefore, it is unclear whether the stronger interactions observed between the right caudal temporal thalamus and left dorsal granular insula compared to ipsilateral interactions have some biological significance. Similar altered interactions emerged for changes in LF-HRV levels and connectivity between the thalamus and basal ganglia, with more pronounced effects for the connectivity of left posterior parietal thalamus and ventromedial putamen.

These findings support the role of thalamus, insula and putamen in autonomic control, as shown in previous studies (17, 20, 30, 53), and add roles for the temporal gyrus whose role in cardiac regulation has been recently suggested (29, 37). Despite the well-established association of activity in amygdala with heart rate fluctuations in task-based experiments (17, 54), no association was observed here. Valenza et al. (29), who also investigated the neural substrates of heart rate in task-free fMRI, found no association of amygdala activity with heart rate either, which may indicate that recruitment of amygdala activity with heart rate occurs mainly during emotional processing tasks rather than the neutral conditions studied; the amygdala traditionally serves affective roles. The neural correlates of heart rate found in our work and in Valenza et al. (29) were not entirely consistent, which may result from the more aggressive physiological correction applied in our work. Artifacts due to cardiac pulsatility were removed using the newly proposed cardiac pulsatility model (55) and systemic low-frequency oscillations were removed through gray matter signal regression (47, 56).

In healthy controls, a seed-to-voxel connectivity analysis revealed that thalamic activity was anticorrelated with core regions of the ventral attention network such as the insula, anterior cingulate cortex (ACC) and supramarginal gyrus, and this anti-correlation was enhanced during elevated levels of HRV (Figure 5). However, the HRV-dependent interplay between the thalamus and ventral attention network was absent in epilepsy. Burianová et al. (57) previously demonstrated a disturbed (static) connectivity between thalamus and the ventral attention network (also referred to as the salience network) in patients with mesial temporal lobe epilepsy which is consistent with our findings (i.e., the insula exhibited increased connectivity with the thalamus and decreased connectivity with the dorsal ACC). However, the present study also shows a strong relationship within healthy individuals between autonomic cardiac regulation and thalamus – ventral attention network connectivity, in line with findings of Chang et al. (27), which is altered in epilepsy.

The thalamus consists of a series of nuclei which are responsible, among others, for the relay of information from cardiovascular receptors to the insular cortex (31). In turn, the insular cortex integrates this information with inputs from ACC, amygdala and high-order polysensory cortex, providing interoceptive awareness. Stimulation of the insula (58, 59), basal ganglia or thalamus (60) lead to marked changes in heart rate and blood pressure. Any impairment in connectivity between these regions, such as found here, may be involved in cardiac rate and variability dysfunction.

A growing body of evidence from functional and structural studies suggests thalamic dysfunction in epilepsy which may underlie the abnormal connectivity of the thalamus with the ventral attention network and the basal ganglia observed in our study (Figures 4, 5). Allen et al. (32), using resting-state fMRI, showed that the nodal participation of thalamus, a measure that reflects the connectivity strength of a region with regions from separate large-scale networks, was increased in epilepsy patients compared to healthy controls, and particularly in patients that succumbed to SUDEP or were at high-risk. Similarly, two recent studies reported altered thalamocortical
Interactions between group (epilepsy vs controls) and LF-HRV states

<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Functional connectivity (FC) for pairs of parcels with strong group (epilepsy vs. controls)-LF-HRV interactions. In the healthy controls (n = 16), the connectivity strength between the right caudal temporal thalamus and left dorsal granular insula [region of ventral attention (VA) network] was lower in times with high levels of LF-HRV (i.e., levels of LF-HRV in the highest quartile of a scan) compared to times with low levels of LF-HRV (i.e., levels of LF-HRV in the lowest quartile of a scan), whereas in people with epilepsy (n = 28) the connectivity strength was higher in times with high levels of LF-HRV. Similar results were observed for the connectivity strength between the left posterior parietal thalamus and the left ventromedial putamen [region of the basal ganglia (BG)].</figcaption></figure>

connectivity (61) as well as impaired connectivity between
thalamus and basal ganglia (62) in individuals with focal to
bilateral tonic-clonic seizures (FBTCS), a group of epilepsy
patients associated with increased risk of seizure-related injuries
and sudden unexpected death. Structural studies have revealed
association of thalamic volume loss with SUDEP and high-
risk patients (63) as well as with patients that present severe
hypoxia during generalized tonic-clonic seizures (64). Moreover,
electrical stimulations of the anterior nucleus of the thalamus
has been shown in clinical trials to reduce seizure frequency
even when seizures are remote from the stimulation site (65–
67). The body of thalamic evidence on mediating seizure
processes, and especially the altered FC between the thalamus
and ventral attention network in epilepsy suggest a target for
intervention. Specialized regions within the thalamus can be
modified by peripheral somatosensory stimulation; activation
of those thalamic sites by active stimulation has the potential

to modify these FC networks, and thus alter the dysfunction
patterns we found here.

RR intervals and, to a less extent, the high frequency
component of HRV (i.e., HF-HRV), were on average lower in
patients compared to controls (Figure 1) which is consistent with
the increased interictal heart rate and decreased HRV reported
in several studies (6–9). A major component of HF-HRV (0.15–
0.40 Hz) is respiratory sinus arrhythmia, a phenomenon where
heart rate fluctuates in synchrony with the breathing cycle at
around 0.3 Hz, and is often considered to reflect parasympathetic
influences on heart rate (68). Therefore, the findings may
indicate reduced parasympathetic influences on the heart in
the patients of our cohort. Interestingly, although RMSSD and
LF-HRV had similar levels in the two groups, when state-
dependent FC was assessed based on these two metrics it revealed
different connectivity patterns between the groups, suggesting
that HRV-state dependent FC has the potential to lead to more
# Differences in FC (seed-based analysis) between low and high LF-HRV state

## Seed: caudal temporal thalamus (R)

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

## Seed: posterior parietal thalamus (L)

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**FIGURE 5** | Differences in seed-based FC between low and high LF-HRV state with seeds placed in the (top) right caudal temporal thalamus and (bottom) left posterior parietal thalamus. The first and second rows of each panel show the fisher-transformed correlations average across all healthy controls and epilepsy</figcaption></figure>

(Continued)
**FIGURE 5** | patients, respectively. Red color indicates higher (toward positive values) connectivity with the seed region in high LF-HRV state whereas blue color indicates lower connectivity. The last row of each panel shows the two-sample *t*-test thresholded at *p* < 0.001 (uncorrected). Even though no significant differences were found between the two groups after correcting for multiple comparison (FWE; *p* < 0.05), the differences observed with *p* < 0.001 (uncorrected) are consistent with the results obtained with the analysis in the atlas space (Figures 3, 4) where the spatial autocorrelation between voxels of the same parcel are implicitly taken into account. The unthresholded statistical maps are available at https://neurovault.org/collections/9452/.

sensitive biomarkers for cardiac dysfunction processes than HRV quantification alone.

This study has limitations that should be considered. While cardiovascular mechanisms are likely impaired in epilepsy and contribute to SUDEP (69), breathing disturbances also appear (70), and may contribute to alterations in FC. Cardiorespiratory arrests monitored via video-electroencephalogram (VEEG) suggest that terminal cardiac arrest was preceded by central apnea in the majority of the cases (71), indicating a potential mediator role for disturbed breathing in cardiac dysfunction. To obtain a more holistic understanding of the neural processes underlying autonomic dysregulation in epilepsy, recognition of the close coupling of respiratory and autonomic control mechanisms should be incorporated in the analysis which was not possible in the present study, as breathing was not monitored during the fMRI scans. Respiratory measures would also be helpful in distinguishing parasympathetic from sympathetic activity in frequency-based HRV measures. HRV parasympathetic activity, lying within the high-frequency range (0.15–0.40 Hz) and associated with respiratory sinus arrhythmia, can decline below 0.15 Hz during periods with low breathing rate, and apnea can completely disrupt respiratory sinus arrhythmia measures. As a consequence, HRV-based measures of parasympathetic and sympathetic activity may be blurred when considering solely cardiac recordings (72).

Several studies describe an inverse relationship between heart rate and HRV measures (73–76). However, as this relationship is not well-understood and heart rate (or RR interval) is already a good measure of ANS activity that can be easily measured, further research is needed to understand the additional information provided with HRV compared to heart rate (73, 77). To this end, in this study, the power spectral density of the HRV that the LF-HRV and HF-HRV measures were derived from, was estimated using the successive difference in RR intervals rather than the RR intervals as this was found to yield HRV measures less correlated with fluctuations in RR interval. In addition, to further disentangle fluctuations in HRV from fluctuations in RR interval, the three HRV measures were orthogonalized with respect to fluctuations in RR interval.

An important caveat of this study in the use of fMRI as a means to study the neural correlates of ANS activity is that there are not well-established methods for disentangling neuronal from physiological effects of autonomic activity (18, 19). While the BOLD (T<sub>2</sub>*) contrast mechanism used in fMRI is a measure sensitive to changes in blood oxygenation induced by local neuronal activity (78), it is also prone to sources of noise that can be categorized to scanner artifacts, motion artifacts, high-frequency physiological artifacts and systemic low-frequency oscillations (79–81). Sources from the first three categories, including fast effects of cardiac pulsatility (~1.0 Hz)

and breathing motion (~0.3 Hz), can be mitigated to a large degree using advanced pulse sequences (e.g., multi-echo fMRI) and noise correction techniques (43–45, 47, 79, 82, 83). However, systemic low-frequency (<0.1 Hz) oscillations which typically refer to BOLD fluctuations driven by changes in heart rate, breathing patterns and blood pressure can be difficult to be separated from neuronal fluctuations as they share the same mechanism; i.e., both neuronal (in an attempt to satisfy increased demands in oxygen) and physiological processes (e.g., heart rate) influence the levels of blood oxygenation (56, 84–87). When studying the neural substrates of the ANS, this is particularly problematic as brain regions not involved in autonomic regulation may share similar BOLD activity with core regions of the CAN due to effects of heart rate in global cerebral blood flow, and therefore the physiological effects of autonomic activity (e.g., fluctuations in heart rate) may lead to artificial connectivity. In this study, to mitigate the effects of systemic low-frequency oscillations, we employed gray matter signal regression which outperforms alternative preprocessing strategies (45, 47, 88, 89). Note that the effects of systemic oscillations are more prominent in visual and sensorimotor areas (85, 89), regions that did not appear to be correlated with heart rate variations in the present study (Figure 2), suggesting that the preprocessing strategy employed here successfully removed the effects of systemic oscillations. However, we cannot exclude the possibility that gray matter regression removed signal of interest as well.

HRV impairment in epilepsy is more pronounced during nocturnal periods (14, 15) and risk for SUDEP is increased during night hours (16). These observations raise the question whether alterations in FC are also enhanced by sleep or during particular phases of the HRV circadian cycles. Note that even though participants often fall asleep during resting-state fMRI and, thus, there might be segments of fMRI data corresponding to sleep in our dataset, the low sample size (N = 44) impedes investigations in relation to sleep effects. The potential for more exaggerated changes in FC during sleep mandates further studies on this issue.

This study represents an exploratory, data-driven approach to investigate whether large-scale networks are involved in cardiac regulation, and is hampered by potential sleep, breathing, and circadian interactions that could interfere with understanding of important brain interactions. A hypothesis-driven analysis, controlling for these interactions may elucidate more precisely the key disruptions in autonomic processes found in epilepsy.

## CONCLUSION

In healthy controls and people with drug-resistant epilepsy, fluctuations in heart rate covaried with brain activity in key
regions of the central autonomic network and in regions associated with cardiac regulation. Functional connectivity of the thalamus with the basal ganglia, a major autonomic regulatory site, and the ventral attention network was strongly linked to levels of LF-HRV, and that relationship differed between healthy controls and epilepsy patients. These findings support a significant role for thalamic contributions to cardiovascular impairments in epilepsy which may lead to cardiac rhythm and blood pressure failings implicated in SUDEP. Because activity in regional thalamic structures can be so readily modified by somatosensory peripheral stimulation, we speculate that the findings suggest a means to interfere with the deficient functional connectivity patterns in epilepsy.

## DATA AVAILABILITY STATEMENT

The datasets presented in this article are not readily available because there is not written informed consent from the patients/participants to share their data with external investigators. Requests to access the datasets should be directed to LL, louis.lemieux@ucl.ac.uk.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by the National Research Ethics Committee (United Kingdom; 04/Q0502/89). The patients/participants provided their written informed consent to participate in this study.

## REFERENCES

1. Acharya UR, Joseph KP, Kannathal N, Lim CM, Suri JS. Heart rate variability: a review. *Med Biol Eng Comput.* (2006) 44:1031-51. doi: 10.1007/s11517-006-0119-0

2. Katona PG, Jih F. Respiratory sinus arrhythmia: noninvasive measure of parasympathetic cardiac control. *J Appl Physiol.* (1975) 39:801-5. doi: 10.1152/jappl.1975.39.5.801

3. Goldberger JJ, Challapalli S, Tung R, Parker MA, Kadish AH. Relationship of heart rate variability to parasympathetic effect. *Circulation.* (2001) 103:1977-83. doi: 10.1161/01.CIR.103.15.1977

4. Haensel A, Mills PJ, Nelesen RA, Ziegler MG, Dimsdale JE. The relationship between heart rate variability and inflammatory markers in cardiovascular diseases. *Psychoneuroendocrinology.* (2008) 33:1305-12. doi: 10.1016/j.psyneuen.2008.08.007

5. Malik M, Camm AJ, Bigger JT, Breithardt G, Cerutti S, Cohen RJ, et al. Heart rate variability. Standards of measurement, physiological interpretation, clinical use. *Eur Heart J.* (1996) 17:354-81. doi: 10.1111/j.1542-474X.1996.tb00275.x

6. Myers KA, Sivathamboo S, Perucca P. Heart rate variability measurement in epilepsy: how can we move from research to clinical practice? *Epilepsia.* (2018) 59:2169-78. doi: 10.1111/epi.14587

7. Ericson M, Ihrman C, Lindblad LE. Heart rate variability in patients with epilepsy. *Epilepsy Res.* (1998) 30:77-83. doi: 10.1016/S0920-1211(97)00094-6

8. Massetani R, Strata G, Galli R, Gori S, Gneri C, Limbruno U, et al. Alteration of cardiac function in patients with temporal lobe epilepsy: different roles of

## AUTHOR CONTRIBUTIONS

MK preprocessed and analyzed data, interpreted the findings, and drafted the initial manuscript. RH interpreted the findings, revised the manuscript, and refined the clinical and physiological interpretation of findings. MG advised on methodological aspects of the study and revised the manuscript. LL and BD designed and conceptualized study, interpreted the findings, and revised the manuscript. All authors contributed to the article and approved the submitted version.

## FUNDING

This work was funded by Epilepsy Research UK (Grant No. P1905; MK, LL, and BD). Funding by the NIH — National Institute of Neurological Disorders and Stroke U01-NS090407 (The Center for SUDEP Research) is acknowledged (BD and RH). This work was supported by the National Institute for Health Research University College London Hospitals Biomedical Research Centre.

## ACKNOWLEDGMENTS

Thanks to Dr. Luke Allen for help. We are grateful to the Epilepsy Society for supporting the Epilepsy Society MRI scanner.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.671890/full#supplementary-material

EEG-ECG monitoring and spectral analysis of RR variability. *Epilepsia*. (1997) 38:363-9. doi: 10.1111/j.1528-1157.1997.tb01129.x

9. Mativo P, Anjum J, Pradhan C, Sathyaprabha TN, Raju TR, Satishchandra P. Study of cardiac autonomic function in drug-naïve, newly diagnosed epilepsy patients. *Epileptic Disord.* (2010) 12:212-6. doi: 10.1684/epd.2010.0325

10. Arbune AA, Jeppesen J, Conradsen I, Ryvlin P, Beniczky S. Peri-ictal heart rate variability parameters as surrogate markers of seizure severity. *Epilepsia*. (2020) 61:S55-60. doi: 10.1111/epi.16491

11. Liu H, Yang Z, Meng F, Guan Y, Ma Y, Liang S. Preoperative heart rate variability as predictors of vagus nerve stimulation outcome in patients with drug-resistant epilepsy. *Sci Rep.* (2018) 8:3856. doi: 10.1038/s41598-018-21669-3

12. DeGiorgio CM, Miller P, Meymandi S, Chin A, Epps J, Gordon S, et al. RMSSD, a measure of vagus-mediated heart rate variability, is associated with risk factors for SUDEP: The SUDEP-7 Inventory. *Epilepsy Behav.* (2010) 19:78-81. doi: 10.1016/j.yebeh.2010.06.011

13. Novak JL, Miller PR, Markovic D, Meymandi SK, DeGiorgio CM. Risk assessment for sudden death in epilepsy: the SUDEP-7 inventory. *Front Neurol.* (2015) 6:1-5. doi: 10.3389/fneur.2015.00252

14. Persson H, Kumlien E, Ericson M, Tomson T. Circadian variation in heart-rate variability in localization-related epilepsy. *Epilepsia.* (2007) 48:917-22. doi: 10.1111/j.1528-1167.2006.00961.x

15. Ronkainen E, Ansakorpi H, Huikuri HV, Myllylä VV, Isojärvi JIT, Korpelainen JT. Suppressed circadian heart rate dynamics in temporal lobe epilepsy. *J Neurol Neurosurg Psychiatry.* (2005) 76:1382-6. doi: 10.1136/jnnp.2004.053777
16. Purnell BS, Thijs RD, Buchanan GF. Dead in the night: sleep-wake and time-of-day influences on sudden unexpected death in epilepsy. *Front Neurol.* (2018) 9:5–8. doi: 10.3389/fneur.2018.01079

17. Beissner F, Meissner K, Bar K-J, Napadow V. The autonomic brain: an activation likelihood estimation meta-analysis for central processing of autonomic function. *J Neurosci.* (2013) 33:10503–11. doi: 10.1523/JNEUROSCI.1103-13.2013

18. Chang C, Raven EP, Duyn JH. Brain-heart interactions: challenges and opportunities with functional magnetic resonance imaging at ultra-high field. *Philos Trans R Soc A Math Phys Eng Sci.* (2016) 374:20150188. doi: 10.1098/rsta.2015.0188

19. Mulcahy JS, Larsson DEO, Garfinkel SN, Critchley HD. Heart rate variability as a biomarker in health and affective disorders: a perspective on neuroimaging studies. *Neuroimage.* (2019) 202:116072. doi: 10.1016/j.neuroimage.2019.116072

20. Thayer JF, Åhs F, Fredrikson M, Sollers JJ, Wager TD. A meta-analysis of heart rate variability and neuroimaging studies: implications for heart rate variability as a marker of stress and health. *Neurosci Biobehav Rev.* (2012) 36:747–56. doi: 10.1016/j.neubiorev.2011.11.009

21. Critchley HD, Mathias CJ, Josephs O, O'Doherty J, Zanini S, Dewar BK, et al. Human cingulate cortex and autonomic control: converging neuroimaging and clinical evidence. *Brain.* (2003) 126:2139–52. doi: 10.1093/brain/awg216

22. Gianaros PJ, Van der Veen FM, Jennings JR. Regional cerebral blood flow correlates with heart period and high-frequency heart period variability during working-memory tasks: implications for the cortical and subcortical regulation of cardiac autonomic activity. *Psychophysiology.* (2004) 41:521–30. doi: 10.1111/1469-8986.2004.00179.x

23. Harper RM, Bandler R, Spriggs D, Alger JR. Lateralized and widespread brain activation during transient blood pressure elevation revealed by magnetic resonance imaging. *J Comp Neurol.* (2000) 417:195–204. doi: 10.1002/(SICI)1096-9861(20000207)417:2<195::AID-CNE5>3.0. CO2;V

24. Harper RM, Gozal D, Bandler R, Spriggs D, Lee J, Alger J. Regional brain activation in humans during respiratory and blood pressure challenges. *Clin Exp Pharmacol Physiol.* (1998) 25:483–6. doi: 10.1111/j.1440-1681.1998.tb02240.x

25. Macey PM, Ogren JA, Kumar R, Harper RM. Functional imaging of autonomic regulation: methods and key findings. *Front Neurosci.* (2016) 9:1–23. doi: 10.3389/fnins.2015.00513

26. Napadow V, Dhond R, Conti G, Makris N, Brown EN, Barbieri R. Brain correlates of autonomic modulation: combining heart rate variability with fMRI. *Neuroimage.* (2008) 42:169–77. doi: 10.1016/j.neuroimage.2008.04.238

27. Chang C, Metzger CD, Glover GH, Duyn JH, Heinze H-JJ, et al. Association between heart rate variability and fluctuations in resting-state functional connectivity. *Neuroimage.* (2013) 68:93–104. doi: 10.1016/j.neuroimage.2012.11.038

28. Jennings JR, Sheu LK, Kuan DCH, Manuck SB, Gianaros PJ. Resting state connectivity of the medial prefrontal cortex covaries with individual differences in high-frequency heart rate variability. *Psychophysiology.* (2016) 53:444–54. doi: 10.1111/psyp.12586

29. Valenza G, Sclocco R, Duggento A, Passamonti L, Napadow V, Barbieri R, et al. The central autonomic network at rest: uncovering functional MRI correlates of time-varying autonomic outflow. *Neuroimage.* (2019) 197:383–90. doi: 10.1016/j.neuroimage.2019.04.075

30. Benarroch EE. The central autonomic network: functional organization, dysfunction, and perspective. *Mayo Clin Proc.* (1993) 68:988–1001. doi: 10.1016/S0025-6196(12)62272-1

31. Palma J-A, Benarroch EE. Neural control of the heart: recent concepts and clinical correlations. *Neurology.* (2014) 83:261–71. doi: 10.1212/WNL.0000000000000605

32. Allen LA, Harper RM, Guye M, Kumar R, Ogren J, Vos S, et al. Altered brain connectivity in sudden unexpected death in epilepsy (SUDEP) revealed using resting-state fMRI. *NeuroImage Clin.* (2019) 24:102060. doi: 10.1016/j.nicl.2019.102060

33. Liao W, Zhang Z, Pan Z, Mantini D, Ding J, Duan X, et al. Altered functional connectivity and small-world in mesial temporal lobe epilepsy. *PLoS ONE.* (2010) 5:27–9. doi: 10.1371/journal.pone.008525

34. Pereira FRS, Alessio A, Sercheli MS, Pedro T, Bilevicius E, Rondina JM, et al. temporal lobe epilepsy: evidence from resting state fMRI. *BMC Neurosci.* (2010) 11:66. doi: 10.1186/1471-2202-11-66

35. Pittau F, Grova C, Moeller F, Dubeau F, Gotman J. Patterns of altered functional connectivity in mesial temporal lobe epilepsy. *Epilepsia.* (2012) 53:1013–23. doi: 10.1111/j.1528-1167.2012.03464.x

36. Waites AB, Briellmann RS, Saling MM, Abbott DF, Jackson GD. Functional connectivity networks are disrupted in left temporal lobe epilepsy. *Ann Neurol.* (2006) 59:335–43. doi: 10.1002/ana.20733

37. Valenza G, Passamonti L, Duggento A, Toschi N, Barbieri R. Uncovering complex central autonomic networks at rest: a functional magnetic resonance imaging study on complex cardiovascular oscillations. *J R Soc Interface.* (2020) 17:20190878. doi: 10.1098/rsif.2019.0878

38. Coan AC, Chaudhary UJ, Grouiller F, Campos BM, Perani S, De Ciantis A, et al. EEG-fMRI in the presurgical evaluation of temporal lobe epilepsy. *J Neurol Neurosurg Psychiatry.* (2016) 87:642–9. doi: 10.1136/jnnp-2015-310401

39. Friston KJ, Ashburner J, Kiebel SJ, Nichols TE, Penny WD, (editors). *Statistical Parametric Mapping: The Analysis of Functional Brain Images*. Academic Press (2007). Available online at: https://www.elsevier.com/books/statistical-parametric-mapping-the-analysis-of-functional-brain-images/penny/978-0-12-372560-8

40. Fan L, Li H, Zhuo J, Zhang Y, Wang J, Chen L, et al. The human brainnetome atlas: a new brain atlas based on connectional architecture. *Cereb Cortex.* (2016) 26:3508–26. doi: 10.1093/cercor/bhw157

41. Leonardi N, Van De Ville D. On spurious and real fluctuations of dynamic functional connectivity during rest. *Neuroimage.* (2015) 104:430–6. doi: 10.1016/j.neuroimage.2014.09.007

42. Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. *Neuroimage.* (2012) 59:2142–54. doi: 10.1016/j.neuroimage.2011.10.018

43. Lemieux L, Salek-Haddadi A, Lund TE, Laufs H, Carmichael D. Modelling large motion events in fMRI studies of patients with epilepsy. *Magn Reson Imaging.* (2007) 25:894–901. doi: 10.1016/j.mri.2007.03.009

44. Salek-Haddadi A, Diehl B, Hamandi K, Merschhemke M, Liston A, Friston K, et al. Hemodynamic correlates of epileptiform discharges: an EEG-fMRI study of 63 patients with focal epilepsy. *Brain Res.* (2006) 1088:148–66. doi: 10.1016/j.brainres.2006.02.098

45. Kassinopoulos M, Mitsis GD. A multi-measure approach for assessing the performance of fMRI preprocessing strategies in resting-state functional connectivity. *bioRxiv.* (2021) 1–61. doi: 10.1101/837609

46. Power JD, Schlaggar BL, Petersen SE. Recent progress and outstanding issues in motion correction in resting state fMRI. *Neuroimage.* (2015) 105:536–51. doi: 10.1016/j.neuroimage.2014.10.044

47. Xifra-Porxas A, Kassinopoulos M, Mitsis GD. Physiological and head motion signatures in static and time-varying functional connectivity and their subject discriminability. *bioRxiv.* (2020). doi: 10.1101/2020.02.04.934554

48. Behzadi Y, Restom K, Liau J, Liu TT. A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. *Neuroimage.* (2007) 37:90–101. doi: 10.1016/j.neuroimage.2007.04.042

49. Allen PJ, Josephs O, Turner R. A method for removing imaging artifact from continuous EEG recorded during functional MRI. *Neuroimage.* (2000) 12:230–9. doi: 10.1006/nimg.2000.0599

50. Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. *Front Public Health.* (2017) 5:1–17. doi: 10.3389/fpubh.2017.00258

51. Worsley KJ, Marrett S, Neelin P, Vandal AC, Friston KJ, Evans AC. A unified statistical approach for determining significant signals in images of cerebral activation. *Hum Brain Mapp.* (1996) 4:58–73. doi: 10.1002/(SICI)1097-0193(1996)4:1<58::AID-HBM4>3.0.CO;2-O

52. Yeo TBT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. *J Neurophysiol.* (2011) 106:1125–65. doi: 10.1152/jn.00338.2011

53. Chapman WP, Livingston RB, E LK. Frontal lobotomy and electrical stimulation of orbital surface of frontal lobes; effect on respiration and on blood pressure in man. *Arch Neurol Psychiatry.* (1949) 62:701. doi: 10.1001/archneurpsyc.1949.02310180002001
54. Yang TT, Simmons AN, Matthews SC, Tapert SF, Bischoff-Grethe A, Frank GKW, et al. Increased amygdala activation is related to heart rate during emotion processing in adolescent subjects. *Neurosci Lett.* (2007) 428:109–14. doi: 10.1016/j.neulet.2007.09.039

55. Kassinopoulos M, Mitsis GD. Physiological noise modeling in fMRI based on the pulsatile component of photoplethysmograph. *bioRxiv.* (2020). doi: 10.1101/2020.06.01.128306

56. Tong Y, Hocke LM, Frederick BB. Low frequency systemic hemodynamic "noise" in resting state BOLD fMRI: characteristics, causes, implications, mitigation strategies, and applications. *Front Neurosci.* (2019) 13:787. doi: 10.3389/fnins.2019.00787

57. Burianová H, Faizo NL, Gray M, Hocking J, Galloway G, Reutens D. Altered functional connectivity in mesial temporal lobe epilepsy. *Epilepsy Res.* (2017) 137:45–52. doi: 10.1016/j.epilepsres.2017.09.001

58. Oppenheimer S. Cerebrogenic cardiac arrhythmias. *Clin Auton Res.* (2006) 16:6–11. doi: 10.1007/s10286-006-0276-0

59. Oppenheimer SM, Gelb A, Girvin JP, Hachinski VC, Med D. Cardiovascular effects of human insular cortex stimulation. *Neurology.* (1992) 42:1727–32. doi: 10.1212/WNL.42.9.1727

60. Thornton JM, Aziz T, Schlugman D, Paterson DJ. Electrical stimulation of the midbrain increases heart rate and arterial blood pressure in awake humans. *J Physiol.* (2002) 539:615–21. doi: 10.1113/jphysiol.2001.014621

61. Caciagli L, Allen LA, He X, Trimmel K, Vos SB, Centeno M, et al. Thalamus and focal to bilateral seizures: a multiscale cognitive imaging study. *Neurology.* (2020) 95:e2427–41. doi: 10.1212/WNL.0000000000010645

62. He X, Chaitanya G, Asma B, Caciagli L, Bassett DS, Tracy JI, et al. Disrupted basal ganglia-thalamocortical loops in focal to bilateral tonic-clonic seizures. *Brain.* (2020) 143:175–90. doi: 10.1093/brain/awz361

63. Allen LA, Vos SB, Kumar R, Ogren JA, Harper RK, Winston GP, et al. Cerebellar, limbic, and midbrain volume alterations in sudden unexpected death in epilepsy. *Epilepsia.* (2019) 60:718–29. doi: 10.1111/epi.14689

64. Allen LA, Harper RM, Vos SB, Scott CA, Lacuey N, Vilella L, et al. Peri-ictal hypoxia is related to extent of regional brain volume loss accompanying generalized tonic-clonic seizures. *Epilepsia.* (2020) 61:1570–80. doi: 10.1111/epi.16615

65. Cooper IS, Upton ARM, Amin I. Reversibility of chronic neurologic deficits. Some effects of electrical stimulation of the thalamus and internal capsule in man. *Stereotact. Funct Neurosurg.* (1980) 43:244–58. doi: 10.1159/000102263

66. Fisher R, Salanova V, Witt T, Worth R, Henry T, Gross R, et al. Electrical stimulation of the anterior nucleus of thalamus for treatment of refractory epilepsy. *Epilepsia.* (2010) 51:899–908. doi: 10.1111/j.1528-1167.2010.02536.x

67. Upton AR, Amin I, Garnett S, Springman M, Nahmias C, Cooper IS. Evoked metabolic responses in the limbicstriate system produced by stimulation of anterior thalamic nucleus in man. *Pacing Clin Electrophysiol.* (1987) 10:217–25. doi: 10.1111/j.1540-8159.1987.tb05952.x

68. Angelone A, Coulter NA. Respiratory sinus arrhythmia: a frequency dependent phenomenon. *J Appl Physiol.* (1964) 19:479–82. doi: 10.1152/jappl.1964.19.3.479

69. Devinsky O, Hesdorffer DC, Thurman DJ, Lhatoo S, Richerson G. Sudden unexpected death in epilepsy: epidemiology, mechanisms, and prevention. *Lancet Neurol.* (2016) 15:1075–88. doi: 10.1016/S1474-4422(16)30158-2

70. Manni R, Terzaghi M, Arbasino C, Sartori I, Galimberti CA, Tartara A. Obstructive sleep apnea in a clinical series of adult epilepsy patients: frequency and features of the comorbidity. *Epilepsia.* (2003) 44:836–40. doi: 10.1046/j.1528-1157.2003.55702.x

71. Ryvlin P, Nashef L, Lhatoo SD, Bateman LM, Bird J, Bleasel A, et al. Incidence and mechanisms of cardiorespiratory arrests in epilepsy monitoring units (MORTEMUS): a retrospective study. *Lancet Neurol.* (2013) 12:966–77. doi: 10.1016/S1474-4422(13)70214-X

72. García-González MA, Vázquez-Seisdedos C, Pallàs-Areny R. Variations in breathing patterns increase low frequency contents in HRV spectra. *Physiol Meas.* (2000) 21:417–23. doi: 10.1088/0967-3334/21/3/307

73. Kazmi SZH, Zhang H, Aziz W, Monfredi O, Abbas SA, Shah SA, et al. Inverse correlation between heart rate variability and heart rate demonstrated by linear and nonlinear analysis. *PLoS ONE.* (2016) 11:1–10. doi: 10.1371/journal.pone.0157557

74. Mangin L, Swynghedauw B, Benis A, Thibault N, Lerebours G, Carré F. Relationships between heart rate and heart rate variability: study in conscious rats. *J Cardiovasc Pharmacol.* (1998) 32:601–7. doi: 10.1097/00005344-199810000-00012

75. Monfredi O, Lyashkov AE, Johnsen AB, Inada S, Schneider H, Wang R, et al. Biophysical characterization of the underappreciated and important relationship between heart rate variability and heart rate. *Hypertension.* (2014) 64:1334–43. doi: 10.1161/HYPERTENSIONAHA.114.03782

76. Yu L, De Mazancourt M, Hess A, Ashadi FR, Klein I, Mal H, et al. Functional connectivity and information flow of the respiratory neural network in chronic obstructive pulmonary disease. *Hum Brain Mapp.* (2016) 37:2736–54. doi: 10.1002/hbm.23205

77. Coumel P, Maison-Blanche P, Catuli D. Heart rate and heart rate variability in normal young adults. *J Cardiovasc Electrophysiol.* (1994) 5:899–911. doi: 10.1111/j.1540-8167.1994.tb01130.x

78. Attwell D, Iadecola C. The neural basis of functional brain imaging signals. *Trends Neurosci.* (2002) 25:621–5. doi: 10.1016/S0166-2236(02)02264-6

79. Caballero-Gaudes C, Reynolds RC. Methods for cleaning the BOLD fMRI signal. *Neuroimage.* (2017) 154:128–49. doi: 10.1016/j.neuroimage.2016.12.018

80. Liu TT. Noise contributions to the fMRI signal: an overview. *Neuroimage.* (2016) 143:141–51. doi: 10.1016/j.neuroimage.2016.09.008

81. Murphy K, Birn RM, Bandettini PA. Resting-state fMRI confounds and cleanup. *Neuroimage.* (2013) 80:349–59. doi: 10.1016/j.neuroimage.2013.04.001

82. Kundu P, Brenowitz ND, Voon V, Worbe Y, Vértes PE, Inati SJ, et al. Integrated strategy for improving functional connectivity mapping using multiecho fMRI. *Proc Natl Acad Sci USA.* (2013) 110:16187–92. doi: 10.1073/pnas.1301725110

83. Liston AD, Lund TE, Salek-Haddadi A, Hamandi K, Friston KJ, Lemieux L. Modelling cardiac signal as a confound in EEG-fMRI and its application in focal epilepsy studies. *Neuroimage.* (2006) 30:827–34. doi: 10.1016/j.neuroimage.2005.10.025

84. Birn RM, Diamond JB, Smith MA, Bandettini PA. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. *Neuroimage.* (2006) 31:1536–48. doi: 10.1016/j.neuroimage.2006.02.048

85. Kassinopoulos M, Mitsis GD. Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration. *Neuroimage.* (2019) 202:116150. doi: 10.1016/j.neuroimage.2019.116150

86. Shmueli K, van Gelderen P, de Zwart J, a., Horovitz SG, Fukunaga M, et al. Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. *Neuroimage.* (2007) 38:306–20. doi: 10.1016/j.neuroimage.2007.07.037

87. Whittaker JR, Driver ID, Venzi M, Bright MG, Murphy K, Chen J, et al. Cerebral autoregulation evidenced by synchronized low frequency oscillations in blood pressure and resting-state fMRI. *Front Neurosci.* (2019) 13:433. doi: 10.3389/fnins.2019.00433

88. Power JD, Plitt M, Gotts SJ, Kundu P, Voon V, Bandettini PA, et al. Ridding fMRI data of motion-related influences: removal of signals with distinct spatial and physical bases in multiecho data. *Proc Natl Acad Sci USA.* (2018) 115:E2105–14. doi: 10.1073/pnas.1720985115

89. Power JD, Plitt M, Laumann TO, Martin A. Sources and implications of whole-brain fMRI signals in humans. *Neuroimage.* (2017) 146:609–25. doi: 10.1016/j.neuroimage.2016.09.038

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Kassinopoulos, Harper, Guye, Lemieux and Diehl. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.