Vahid Khalilzad Sharghi
Eric A. Maltbie
Wen-Ju Pan
Shella D. Keilholz
Kaundinya S. Gopinath- 1Department of Biomedical Engineering, Emory University-Georgia Tech, Atlanta, GA, United States
- 2Department of Radiology & Imaging Sciences, Emory University, Atlanta, GA, United States
A number of studies point to slow (0.1–2 Hz) brain rhythms as the basis for the resting-state functional magnetic resonance imaging (rsfMRI) signal. Slow waves exist in the absence of stimulation, propagate across the cortex, and are strongly modulated by vigilance similar to large portions of the rsfMRI signal. However, it is not clear if slow rhythms serve as the basis of all neural activity reflected in rsfMRI signals, or just the vigilance-dependent components. The rsfMRI data exhibit quasi-periodic patterns (QPPs) that appear to increase in strength with decreasing vigilance and propagate across the brain similar to slow rhythms. These QPPs can complicate the estimation of functional connectivity (FC) via rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel-time courses of functionally connected brain regions. In this study, we examined the relationship between cortical slow rhythms and the rsfMRI signal, using a well-established pharmacological model of slow wave suppression. Suppression of cortical slow rhythms led to significant reduction in the amplitude of QPPs but increased rsfMRI measures of intrinsic FC in rats. The results suggest that cortical slow rhythms serve as the basis of only the vigilance-dependent components (e.g., QPPs) of rsfMRI signals. Further attenuation of these non-specific signals enhances delineation of brain functional networks.
Introduction
Blood oxygenation level-dependent (BOLD) resting-state fMRI (rsfMRI) is a powerful tool for mapping the brain connectomes during normal function and illness (Bullmore and Sporns, 2009; Buckner et al., 2013). Numerous studies have shown that the performance in cognitive and sensory tasks is strongly influenced by the resting-state functional connectivity (FC) in relevant networks (Smith et al., 2009; Laird et al., 2011; Tavor et al., 2016). Impairments in brain FC networks have been shown to serve as biomarkers in a number of disease models (Greicius, 2008; Jalilianhasanpour et al., 2019). However, the neurophysiological basis underlying rsfMRI signals are not completely understood, which impedes the interpretation of these studies.
Brain slow rhythms and resting-state functional magnetic resonance imaging (rsfMRI)
Several studies have surmised that rsfMRI signals are related to expression of cortical brain rhythms (He and Raichle, 2009; Scholvinck et al., 2010; Wang et al., 2012; Amemiya et al., 2016; Leong et al., 2016; Lu et al., 2016; Matsui et al., 2016; Wen and Liu, 2016; Chan et al., 2017), especially in gamma rhythms (Scholvinck et al., 2010; Liu et al., 2015a,b) and slow/delta rhythms (Lu et al., 2016; Matsui et al., 2016; Chan et al., 2017) frequency bands. Of particular interest are sub-1 Hz (∼0.1–1 Hz) fluctuations in delta rhythms (1–4 Hz) (Steriade et al., 1993a,b; Steriade, 2000; Crunelli and Hughes, 2010). From here onward, we refer to sub-1Hz waves as slow rhythms, and coherent oscillations in the 1–4-Hz band delta rhythms to avoid confusion. Slow rhythms share many properties with rsfMRI signals. They exist in the absence of stimulation, during wakefulness (Vyazovskiy et al., 2011; McCormick et al., 2015), sleep (Massimini et al., 2004; Crunelli and Hughes, 2010; Nir et al., 2011; David et al., 2013), and under sedation (Steriade et al., 1993b; Volgushev et al., 2006; Chauvette et al., 2011; David et al., 2013) and are influenced by vigilance/arousal (Steriade, 2000; Llinas and Steriade, 2006; Chauvette et al., 2011; McCormick et al., 2015), and propagate across the brain (Massimini et al., 2004; Sheroziya and Timofeev, 2014; Matsui et al., 2016). These properties are seen in rsfMRI signal patterns (Majeed et al., 2011; Mitra et al., 2015; Matsui et al., 2016; Yousefi et al., 2018; Liu and Falahpour, 2020) as well. In a series of studies (Leong et al., 2016; Chan et al., 2017; Wang et al., 2019), one group has shown that that low-frequency (e.g., 1 Hz) optogenetic stimulation of specific hippocampal and thalamic regions in rats evoked increased interhemispheric homotopic FC (IHFC) in various brain regions, whereas high-frequency stimulation of the same hippocampal region at 40 Hz had no effect on FCs (Chan et al., 2017; Wang et al., 2019), indicating that FC of rsfMRI signals may reflect neuronal dynamic related to slow rhythms. Another study in mice (Matsui et al., 2016) found that hemodynamic signals were closely associated with transient neuronal co-activations embedded in slow globally propagating waves of neural activity during rest, strengthening the view that the rsfMRI signals reflect slow wave activity. However, resting-state FC can in principle also be generated through other mechanisms such as gamma rhythms (Scholvinck et al., 2010) and spontaneous fluctuations in action potential (Ponce-Alvarez et al., 2015) between functionally connected brain regions. Importantly, unlike slow/delta rhythms, the strength of FC in canonical brain function networks generally decreases with reductions in vigilance and arousal levels (Larson-Prior et al., 2009; Hutchison et al., 2014; Bettinardi et al., 2015). Relatedly, a recent optical imaging study (Brier et al., 2019) indicated that the slow globally propagating neural activity in resting-state calcium and hemoglobin signals is separable from those that drive correlations between functionally connected brain regions. Another study (Okun et al., 2019) showed that the populations of neurons, which exhibit coherent fluctuations across the brain over slow and fast time-scales, are distinct from each other. Thus, a number of indicators exist that suggest rsfMRI signals encoding FC in brain functional and structural networks are distinct from fMRI signals reflecting slow rhythms. It is also possible that these rhythms serve as the basis of just the vigilance-dependent components of rsfMRI signals. For instance, the rsfMRI data exhibit quasi-periodic patterns (QPPs) wherein sets of cortical brain regions exhibit transient coherent activity (Thompson et al., 2014, 2015; Yousefi et al., 2018; Abbas et al., 2019). QPPs increase in strength with decreasing vigilance/arousal (Thompson et al., 2013a; Billings and Keilholz, 2018) and propagate across the brain (Majeed et al., 2009; Belloy et al., 2018) similar to brain slow rhythms. Another vigilance/arousal-dependent feature of rsfMRI signals is a transient cortex-wide neural signal component called global signal (Scholvinck et al., 2010; Falahpour et al., 2018) that induces correlations between voxel time-series across the brain that are not specific to any brain function networks. Global signal also increases with decreasing vigilance (Falahpour et al., 2018). In fact, some studies have shown that global signal is related to QPPs (Belloy et al., 2018). Both QPP and global signal complicate/confound the estimation of intrinsic FC in brain function networks.
Role of T-type calcium channels in generation of slow rhythms
Given the contradictory hypotheses and findings regarding the relationship of BOLD rsfMRI signal and slow rhythms, a critical need exists to investigate the relationship between them. One way to do this is by examining the effects of systematically manipulating the expression of slow rhythms on rsfMRI signals. This can be achieved through pharmacological interventions. Although the slow rhythms are primarily generated through interaction of excitatory and inhibitory neurons in the cortex (Crunelli et al., 2015; McCormick et al., 2015), one mechanism for expression and maintenance of slow rhythms in the brain is through the interplay between cortical and thalamic slow rhythm generating oscillators (Crunelli et al., 2015). This mechanism is reviewed in detail in Crunelli and Hughes (2010) and Crunelli et al. (2015). In brief, cortical slow rhythms are characterized by long periods (“DOWN” states) of low-neural activity interspersed with “UP” states of increased excitability. The prolonged “UP” states of slow rhythms in cortical layer 5/6 neurons lead to long-lasting excitatory postsynaptic potentials (EPSPs) in thalamocortical neurons. These EPSPs induce slow rhythms in thalamic nuclei, which receive input from the cortical layer 5/6 neurons. The “UP” states in thalamus are preceded by a high-frequency burst generated by low-threshold T-type calcium (Ca2+) channels (TTCCs). These high-frequency bursts generate a new “UP” state of slow rhythms in cortical layer 4 neurons. This cycle of cortex to thalamus to cortex induction of UP and DOWN states helps maintain the rhythms. Slow waves in different thalamocortical circuits are not fully synchronized, which likely provide the basis for propagative properties of slow rhythms (Crunelli et al., 2015). A selective TTCC blocker (Dreyfus et al., 2010) TTA-P2 (3,5-dichloro-N-[1-(2,2-dimethyl-tetrahydro-pyran-4-ylmethyl)-4-fluoro-piperidin-4-ylmethyl]-benzamide), has been shown to inhibit this thalamic bursting activity, resulting in suppression of slow rhythms (Crunelli and Hughes, 2010; Crunelli et al., 2014). Selective blockade of voltage-gated Cav3.1 TTCCs in first-order thalamic nuclei leads to significant reductions in the power of slow waves in corresponding primary cortices in rats (David et al., 2013). Systemic administration of the TTA-P2 induces an even greater (up to 60%) suppression of slow cortical oscillations (David et al., 2013) in rats, both during natural sleep and under anesthesia. TTA-P2 reversibly inhibits TTCCs but does not engage L-type calcium channels, sodium channels, or potassium channels and does not interfere with glutamatergic or GABAergic synaptic currents (Dreyfus et al., 2010). Hence, it provides ideal means for examining the effects of suppression of slow waves without affecting other forms of neural activity not related to TTCCs, e.g., tonic neuronal firing that subserve cognition (Sherman, 2001). Importantly, though TTA-P2 administration suppresses slow rhythms, it paradoxically increases somnolence in rats (Kraus et al., 2009; McCafferty et al., 2012; Crunelli et al., 2014). TTA-P2 does not have any renal or cardiovascular effects (Shipe et al., 2008), and it does not affect arterial blood flow (Masicampo et al., 2018). Hence, it is not likely to induce global changes in cerebral blood flow.
In this preliminary study, we examined the effects of TTA-P2 on rsfMRI signals in rodents. Based on the studies described above, we hypothesized that systemic administration of TTA-P2 will induce significant decrease in the strength of rsfMRI QPPs. And this reduction in the strength of signals, not specific to brain function networks, will increase rsfMRI estimates of intrinsic FC, thereby enhancing the power of rsfMRI technique to probe brain function networks.
Materials and methods
Rodents and preparation
All protocols for animal studies were reviewed and approved by the Emory University Institutional Animal Care and Use Committee (IACUC) and were in compliance with the NIH guidelines. Twelve Sprague-Dawley rats (male, 300–350 g) were employed in this study. Seven rats were administered the TTCC TTA-P2, and five were administered the vehicle. Animal preparation and imaging were conducted with long established techniques (Pan et al., 2013; Thompson et al., 2014, 2015) and more details are provided in Supplementary materials. The fMRI data were collected under dexmedetomidine (Dexmed) anesthesia using long-established techniques (Thompson et al., 2013b,2014, 2015). Dexmed was administered in the form of a 0.025-mg/kg bolus and maintained under sedation with a constant subcutaneous (s.c.) infusion of 0.05 mg/kg/h. Pancuronium bromide was simultaneously administered at a constant rate of 1 mg/kg/h to minimize animal motion. Body temperature, respiration rate, end-tidal CO2, blood oxygen saturation (SpO2), and heart rate were all continuously monitored. After acquiring rsfMRI data under these conditions for 20–90 min, seven rats were injected (s.c.) with 6 mg/kg TTA-P2 dissolved in 2.5 ml of vehicle (4% dimethyl sulfoxide (DMSO) saline) and studied for another 60–90 min. Similarly, data from 5 other rats were acquired with just the vehicle control.
Magnetic resonance (MR) imaging
The MRI data were acquired in a 9.4-T Bruker MRI scanner under Dexmed anesthesia regime as described above. A custom-made oval surface transceiver radiofrequency (RF) coil with an internal diameter of 2 cm × 2.5 cm (Pan et al., 2021) was placed directly over the rat brain to acquire MRI signals. The fMRI scans were acquired with a gradient-echo EPI sequence with the following parameters: TR = 2,000 ms, TE = 20 ms, flip angle = 60°, matrix size 70 × 70, field of view 3.45 cm × 3.45 cm. Twenty-four 0.5-mm thick coronal slices covering the almost the entire brain (see Supplementary Figure 1). A 2D fast spin echo MRI scan (TR/TE = 3,500 ms/44 ms; 24 coronal slices) was also acquired with the same resolution and coverage as the EPI scans.
Data analysis
Functional magnetic resonance imaging (fMRI) data preprocessing
An fMRI data preprocessing was conducted with FSL (Smith et al., 2004) and AFNI (Cox, 1996) tools, and in-house MATLAB scripts. The rsfMRI voxel time-series were corrected for magnetic susceptibility-induced geometric distortions (Andersson et al., 2018), temporally shifted to account for differences in slice acquisition times and 3D volume registered to a base volume to account for global rigid motion. These time-series were co-registered to the T2-weighted high-resolution anatomic scan, and spatially normalized to the Paxinos atlas template (Valdes-Hernandez et al., 2011; Paxinos and Watson, 2014), with the warp computed from alignment of the high-resolution T2-weighted anatomic to the Paxinos atlas template. The spatially normalized fMRI time-series were bandpass (0.01–0.20 Hz) filtered, detrended of similarly filtered estimated motion parameters time-courses through linear regression, and standardized. The resultant time-series were spatially smoothed with an isotropic Gaussian filter (FWHM = 1 mm).
Quasi-periodic pattern (QPP) analysis
Since global signal regression (GSR) enhances the delineation of QPPs (Yousefi et al., 2018; Abbas et al., 2019), the whole-brain average fMRI signal was filtered (0.01–0.20 Hz) and regressed out during the preprocessing fMRI time-series (see Section “Functional magnetic resonance imaging (fMRI) data preprocessing”). Finally, QPPs were extracted from the concatenated pre-injection (Baseline) of all rats’ fMRI data. Supplementary Figure 2 illustrates the process of identifying the expression of a QPP from a representative rsfMRI time-series dataset. Briefly, the QPP algorithm (Majeed et al., 2011) randomly selects an epoch (initial guess for a QPP template) of spatiotemporal data (5- to 10-s duration) and calculates sliding-window correlation between the template and the entire time-course, identifying time points of high correlation (with a set threshold arrived at with bootstrapping; see Majeed et al., 2011). The spatiotemporal blocks centered on these time points are then averaged to create a revised spatiotemporal template of that QPP, and this process is repeated until the template converges. The changes in the strength of the expression of QPPs over time for each fMRI scan (Baseline and TTA-P2) for each rat were estimated through the sliding window spatiotemporal correlation (STC) of the corresponding fMRI time-series with the rats’ QPP template (Majeed et al., 2011). The effects of TTA-P2 on QPPs were assessed with between-session (TTA-P2 vs. Baseline) paired t-tests on the mean of positive excursions of the STC curve above zero. Similar analysis was performed to assess differences between Vehicle and Baseline sessions for the five rats, which were only administered the vehicle.
Interhemispheric homotopic functional connectivity (FC) analysis
The FCs in intrinsic brain function networks were assessed by examining the FC between homotopic brain regions in the two hemispheres, since these are in general known to possess high FC (Stark et al., 2008; Shen et al., 2015; Wang et al., 2019). Further IHFC has been employed as a measure of intrinsic FC in some recent studies examining the neurophysiological basis of rsfMRI, which implicates slow rhythms as a likely basis for intrinsic FC (Matsui et al., 2016; Chan et al., 2017). To examine the effects of slow wave suppression on fMRI IHFC, the preprocessed time-series (see Section “Functional magnetic resonance imaging (fMRI) data preprocessing”) was broken down into 10-min segments (in order to examine any systematic evolution of the FC over time). Subsequently, region of interest (ROI) average reference vectors were obtained from the 48 homotopic brain regions that were labeled in the Paxinos atlas space in a well-established rat MRI template (Valdes-Hernandez et al., 2011). IHFC for each region (for a given 10-min segment) was assessed during Baseline and TTA-P2/Vehicle sessions through the cross-correlation coefficient (CC) between the left and right hemisphere ROIs of the corresponding regions in a general linear regression framework (Cox, 1996). Since the dynamic effects of TTA-P2 on FC were found to be transient and variable across the rats, the 10-min segment for each rat where most of the homotopic ROIs showed maximum FC was selected as the one where the effects of TTA-P2 was strongest across the brain. Since all the 10-min segments in the Baseline condition are equivalent, Baseline IHFC was arrived at by taking the median of the IHFC in all 10-min segments. The effects of TTA-P2 (and Vehicle) were assessed through paired t-test on the IHFC of different brain regions. Multiple-comparison correction was conducted through false discovery rate (FDR) (Benjamini, 2010). Further reproducibility of TTA-P2 effects across the rats for each ROI was assessed via the number of rats, which exhibited increased IHFC after TTA-P2 administration as a fraction of the total number of rats.
Thalamocortical functional connectivity (FC) analysis
Since the mechanism of slow rhythms suppression involves blockade of thalamic TTCCs, the FCs in different thalamocortical networks were estimated through the CC between ipsilateral Paxinos ROIs of primary somatosensory, motor, auditory and visual cortex and their corresponding thalamic centers (ventrobasal (VB) nucleus, ventrolateral (VL) nucleus, medial geniculate nucleus (MGN), and lateral geniculate nucleus (LGN) ROIs). Similarly, IHFC of these four first-order thalamic nuclei and intrathalamic FC between them were also obtained. TTA-P2 effects on FC were assessed as described above (Section “Interhemispheric homotopic functional connectivity (FC) analysis”). Intrathalamic FC were assessed between LGN and the ventrobasal (VB) nucleus complex, and between MGN and ventrolateral (VL) nucleus. This was because VB and VL nuclei, and LGN and MGN are adjacent, which could lead to overlap between some of their fMRI voxels, and vasculature. The thalamic nuclei ROIs were obtained using a publicly available atlas (Papp et al., 2014).
Seed-based cross-correlation maps
Finally, the whole-brain FC maps of the rat barrel cortex (S1BF) and auditory cortex [both of which are unimodal regions, which are known to possess high IHFC (Stark et al., 2008)], were obtained through CC between all voxels in the brain and the right S1BF and right auditory cortex seed ROI average fMRI time-series vectors, for all 10-min fMRI segments of each session. TTA-P2/Vehicle vs. Baseline differences in the FC group paired t-tests are as described above (Section “Interhemispheric homotopic functional connectivity (FC) analysis”). Statistical inference test maps were clustered, and the significance of activations, accounting for multiple comparisons, was derived by means of Monte Carlo simulation of the process of image generation, spatial correlation of voxels, intensity thresholding, masking, and cluster identification (Cox et al., 2017; Gopinath et al., 2018).
Results
None of the rats exhibited significant motion (i.e., frame-to-frame displacement was always less than 0.02 mm). Furthermore, no meaningful changes were observed in the physiological signals monitored: respiration rate, hear rate, oxygen saturation, and body temperature between Baseline and TTA-P2/Vehicle sessions.
TTA-P2 administration reduced strength of quasi-periodic patterns (QPPs)
Figure 1A shows a QPP with anterior-to-posterior propagation. The fluctuation in the strength of the QPP before and after injection of TTA-P2 is shown for a representative rat administered TTA-P2 (Figure 1B), and one administered vehicle (Figure 1C). TTA-P2 administration significantly (paired t-test, p < 0.002) reduced the strength (mean of positive STC values) of QPPs compared to Baseline. All the rats exhibited suppression of QPPs after TTA-P2 administration. The amount of suppression of QPPs induced by TTA-P2 varied from 18 to 58% (mean 48%) across the rats (see Figure 1D). Performing QPP analysis without GSR (Supplementary Figure 3) reduced the significance of TTA-P2 vs. Baseline differences to paired t-test, p < 0.005. This analysis yielded 8–37% (mean 27%) suppression in QPPs after TTA-P2 administration. Vehicle did not appreciably alter the strength of the QPPs, with or without GSR.
Figure 1. (A) Eight frames from the quasi-periodic pattern (QPP) template obtained from both groups’ pre-injection concatenated functional time-series (window length = 10 s). (B,C) The evolution of the QPP strength with time assessed with spatiotemporal correlation of the fMRI time-series with corresponding QPP template on one rat (B) after systematic administration of TTA-P2; and one rat (C) after systematic administration of Vehicle. (D) The QPP strength changes for each rat in the TTA-P2 group. The values are estimated as the mean of positive excursions of the spatiotemporal correlation (STC) curve above zero, normalized by the maximum correlations for each subject.
TTA-P2 administration increased interhemispheric homotopic functional connectivity (IHFC)
Next, we examined IHFC of the Paxinos atlas regions. This Paxinos template (Valdes-Hernandez et al., 2011) has 48 homotopic ROIs (as listed in Supplementary Table 1), but 8 of these ROIs did not yield any fMRI signal in one or both hemispheres. Thirty-one of the remaining forty ROIs exhibited (Figure 2) significantly [paired t-test, p < 0.05 (Table 1)] increased IHFC after injection of TTA-P2. FDR correction rendered TTA-P2 vs. Baseline differences in IHFC in 4 out of these 31 ROIs above significance (FDR q = 0.05) after accounting for multiple comparisons. Reproducibility of the IHFC increases after TTA-P2 administration varied between 86% and 100% (mean 93%) among the IHFCs of the ROIs, which exhibited significant increase in IHFC after TTA-P2 administration (Supplementary Figure 4). Analyzing IHFC after GSR yielded TTA-P2 engendered enhancements in 22 ROIs after FDR correction (Supplementary Table 2). Vehicle administration did not evoke appreciable (significant) changes in IHFC.
Figure 2. Paxinos region of interest (ROIs) exhibiting significantly (p < 0.05) enhanced interhemispheric homotopic functional connectivity (IHFC) relative to Baseline after TTA-P2 injection overlaid on a Paxinos atlas brain. ROIs shown in blue did not achieve significance.
Table 1. TTA-P2 vs. baseline interhemispheric homotopic functional connectivity (IHFC) paired t-test for different Paxinos region of interest (ROIs).
TTA-P2 administration did not affect thalamocortical functional connectivity (FC)
The FC between primary sensory and motor cortices and their corresponding first-order thalamic nuclei were not significant (1-sample t-tests, p > 0.05), both in Baseline and TTA-P2 sessions. TTA-P2 administration did not have appreciable effect on these connections. The IHFC of these thalamic nuclei, and FC between VB and LGN thalamus, and between MGN and VL thalamus increased with TTA-P2 administration, but the results were not significant at paired t-test, p < 0.05. However, these intrathalamic FCs were significantly strong (1-sample t-tests, p < 0.05) in both sessions. GSR did not change the results appreciably.
TTA-P2 administration increased functional connectivity (FC) in canonical brain function networks
The TTA-P2 significantly (paired t-test, FWE α < 0.05) increased the rsfMRI FC between right S1BF and some areas in somatosensory, motor, auditory, visual, and parietal cortices bilaterally (Figure 3). Right auditory cortex exhibited significantly increased FC to contralateral auditory, visual, and somatosensory areas (Figure 4). Both seed ROIs exhibited strong FC with most of the above functionally connected regions in the Baseline session (Supplementary Figures 5, 6). These effects persisted after GSR (Supplementary Figures 7, 8), but the extent and strength of the TTA-P2 vs. Baseline FC differences decreased. Vehicle administration did not evoke appreciable changes in FC.
Figure 3. TTA-P2 vs. Baseline t-statistic maps highlighting regions with enhanced functional connectivity (FC) to right S1BF ROI after TTA-P2 administration. The slice-location coordinates are in Paxinos space. Left hemisphere is on the left-hand side of the maps.
Figure 4. TTA-P2 vs. Baseline t-statistic maps highlighting regions with enhanced functional connectivity (FC) to right auditory cortex region of interest (ROI) after TTA-P2 administration. The slice-location coordinates are in Paxinos space. Left hemisphere is on the left-hand side of the maps.
Discussion
Context and rationale for study
Although 20–40% of the rsfMRI signal can be attributed to vasomotion, and physiological noise, there remains a large portion of the rsfMRI signal that can be attributed to neural signals (Mayhew et al., 1996; Biswal et al., 1997; Mitra et al., 1997; Birn et al., 2006; Bianciardi et al., 2009). The neurophysiological basis of these signals are not clear with the first studies directly relating rsfMRI signal to neural activity published not long ago (Leopold et al., 2003; Nir et al., 2008). Some groups have hypothesized that rsfMRI signal underlying FC in brain function networks could arise from spontaneous neural activity in these networks in the absence of any tasks (Chan et al., 2015; Ponce-Alvarez et al., 2015). On the other hand, many groups have observed correlations between fluctuations in the power of different brain rhythms obtained from electrophysiological recordings and BOLD fMRI signal (Scholvinck et al., 2010; Pan et al., 2013; Liu et al., 2015b; Lu et al., 2016). However, it must be noted that these observations do not clarify whether fMRI signal encode all of the neuronal dynamics that compose brain rhythms (Matsui et al., 2016), or discrete events within or driving these rhythms (Schwalm et al., 2017). It must be noted that different components of the rsfMRI signal could reflect different neural phenomena. For instance, fMRI signals linked to slow rhythms could be distinct from those that serve as the basis for rsfMRI FC in canonical brain function networks. In contrast with slow rhythms, the strength of FC in resting-state brain function networks generally decreases with reductions in vigilance and arousal levels (Larson-Prior et al., 2009; Hutchison et al., 2014). Relatedly, results from a few recent studies (Brier et al., 2019; Okun et al., 2019) indicated that hemodynamic signals related to globally propagating slow waves of neural activity are distinct from those that induce correlations between functionally connected brain regions. Hence, we hypothesized that only the arousal/vigilance-dependent components of rsfMRI signals arise from neural events related to expression of slow rhythms. And removing these rsfMRI signal components will enhance the delineation of brain function networks.
TTA-P2 attenuates functional magnetic resonance imaging (fMRI) quasi-periodic patterns
The fMRI data exhibit QPPs wherein sets of cortical brain regions exhibit transient coherent activity (Majeed et al., 2009, 2011). These patterns propagate across the brain similar to slow rhythms brain (Majeed et al., 2009; Belloy et al., 2018). In addition, just like slow rhythms, QPPs increase in strength with decreasing vigilance/arousal (Thompson et al., 2013a; Billings and Keilholz, 2018), QPPs also complicate the estimation of FC with rsfMRI, either by existing as unmodeled signal or by inducing additional wide-spread correlation between voxel time-courses of functionally disparate brain regions. To examine the dependence of QPPs on slow rhythms, we employed a well-established pharmacological model, which has been shown to reproducibly suppress slow rhythms in rats (Dreyfus et al., 2010; Crunelli et al., 2015). As described in “Role of T-type calcium channels in generation of slow rhythms,” slow waves are maintained through the interplay between cortical and thalamic slow rhythm generating oscillators (Crunelli and Hughes, 2010; Crunelli et al., 2015). Importantly, the initiation of UP states of slow rhythms in thalamus is preceded by a TTCC-mediated high-frequency burst of thalamocortical EPSPs that initiate a new UP state of slow rhythm in the cortex. Administration of the TTCC blocker TTA-P2 suppresses this high-frequency burst of thalamocortical EPSPs, thereby suppressing slow rhythms (Dreyfus et al., 2010; David et al., 2013). In our experiment, systemic administration of TTA-P2 led to significant attenuation in the strength of the QPPs. These results indicate that putative suppression of slow rhythms diminishes the strength of QPPs, thereby reducing the arousal/vigilance-dependent contributions to fMRI signal. This effect was observed in all seven rats, which were administered TTA-P2 in this experiment, indicating that it is very reproducible. The suppression of QPPs occurred almost immediately after TTA-P2 administration. The absence of electrophysiological recordings time-locked to the fMRI signal in this study rendered it impossible to estimate relationship between QPPs and slow rhythms. However, one group (Schwalm et al., 2017) has shown cortex-wide correlations between fMRI signal and onset of cortical slow waves. This indicates that the global increases in signal often seen with transition between vigilance/arousal states (Liu et al., 2015a; Liu and Falahpour, 2020) could be related to the high-frequency burst of EPSPs that thalamus sends to the cortex at onsets of slow wave UP states. Similarly, QPPs could be related to the propagation of these slow wave onsets across the cortex (Crunelli et al., 2015). The fact that the strength of both global signal and QPPs are low during wakefulness could be related to the decreased frequency of thalamic bursting events, and hence onsets of slow wave UP states during wakefulness (Sherman et al., 2006; Vyazovskiy et al., 2011).
TTA-P2 enhances functional connectivity (FC) in brain function networks
In addition to suppression of QPPs, TTA-P2 administration led to increased IHFC between homotopic cortical ROIs, especially those belonging to unimodal regions (e.g., primary sensory and motor cortices), which have been shown to possess highest IHFCs in the brain (Stark et al., 2008). Nine Paxinos cortical ROIs did not exhibit significant increase in IHFC after TTA-P2 administration. Seven of these nine ROIs were in heteromodal regions (frontal and associative cortices), which generally do not exhibit strong IHFC in rsfMRI studies (Stark et al., 2008; Shen et al., 2015). The remaining two were very small ROIs, thereby more susceptible to noise. The reproducibility of these results was not as strong as the suppression of QPPs (Supplementary Figure 4). Reproducibility varied between 86 and 100% (mean 93%) among the IHFCs of the ROIs, which exhibited significant increase in IHFC after TTA-P2 administration. This lack of perfect reproducibility could be due to variations in effects of anesthesia and/or the drug across the rats. The effects of the drug were variable with some rats exhibiting immediate decreases in IHFC after TTA-P2 injection, whereas in some cases, IHFC increases were more delayed.
Similar to IHFC, intrinsic FC in brain areas functionally connected to rat barrel and auditory cortices also increased after TTA-P2 administration. TTA-P2 significantly increased the rsfMRI FC between right S1BF and some areas in somatosensory, motor, auditory, visual, and parietal cortices bilaterally. Anatomical and functional connections between rat barrel cortex and these brain areas have been observed in a number of studies (Frostig et al., 2008; Bedwell et al., 2014; Zakiewicz et al., 2014). Right auditory cortex exhibited significantly increased FC to contralateral auditory, visual, and somatosensory areas, consistent with prior functional and anatomic tracing studies (Smith et al., 2010; Schormans et al., 2016; Thomas et al., 2020). Thus, the suppression of arousal-dependent fMRI signals (QPPs) by TTA-P2 increases apparent intrinsic FC in cortical brain function networks as hypothesized.
Finally, thalamocortical FC between first-order thalamic nuclei and cortical regions to which they project to did not change appreciably after TTA-P2 administration. However, the strengths of thalamocortical FC in both Baseline and TTA-P2 sessions were also low. Studies have shown that the light anesthetic regimen employed in our study could be analogous to slow wave sleep, both in terms of the slow rhythms (Torao-Angosto et al., 2021) and FC patterns (Greicius et al., 2008; Bettinardi et al., 2015). Thalamus can be decoupled from cortex (in terms of tonic thalamic activity) during some general anesthesia conditions (Suzuki and Larkum, 2020), and thalamocortical connectivity is often reduced during anesthesia and deep sleep (Esser et al., 2009; Hudetz, 2012; Picchioni et al., 2014). On the other hand, IHFC of thalamic nuclei and intrathalamic FC, both increased after TTA-P2 administration, though the differences were not significant (paired t-test, p > 0.05). Furthermore, both IHFC and intrathalamic FC were significant and strong in both Baseline and TTA-P2 session. Since there are no interhemispheric homotopic monosynaptic connections between these thalamic nuclei, this indicates that the strong rsfMRI IHFC observed between them is polysynaptic in nature. Interestingly, thalamic fMRI signal has been observed to decrease with arousal/vigilance in contrast to cortical global signal and QPPs in many studies (Chang et al., 2016; Liu and Falahpour, 2020). Blockade of TTCCs suppress slow rhythms in thalamus along with the cortex. However, thalamic signal increased in our study after TTA-P2 administration (i.e., slow wave suppression). Hence, decreases observed in thalamic signal with vigilance/arousal (Chang et al., 2016; Liu and Falahpour, 2020) are likely not related to slow rhythms.
Limitations and pitfalls
Since this was a preliminary study, there are still a lot of unresolved issues. The major drawback of this study is that electrophysiological data were not acquired during the fMRI scans. Hence, there is no direct evidence that cortical slow rhythms were attenuated by TTA-P2. Although a well-established pharmacological model (David et al., 2013) for slow wave suppression was employed, lack of electrophysiological recordings time-locked to the fMRI data hindered both confirmation of slow wave suppression and ability to obtain quantitative relationships between slow rhythms and fMRI QPPs. Furthermore, although strength of QPPs decreased after TTA-P2 administration in all rats, the magnitude of suppression of QPPs varied significantly across the rats. In addition, the dynamic changes in IHFC with time exhibited large variations across the rats in both Baseline and TTA-P2 sessions. Obtaining simultaneous electrophysiological data could have helped to separate the variabilities that could be attributed to drug and anesthetic effects.
Second, although physiological parameters: respiration rate, end-tidal CO2, hear rate, oxygen saturation, and body temperature were monitored, these data could not be digitally recorded and time-locked to the fMRI data. Observationally, no substantial changes were observed in these physiological responses between Baseline and TTA-P2 sessions, which is consistent with previous studies (Shipe et al., 2008; David et al., 2013). However, the contributions of these physiological responses, which are known to affect BOLD fMRI data (Birn et al., 2006; Chang et al., 2009), could not be regressed out. Physiological noise was minimized before estimation of QPPs with GSR. Performing QPP analysis without GSR did not change the conclusion that QPP strengths were significantly attenuated upon TTA-P2 administration. The FC results presented in Figures 2–4 were obtained without GSR due to well-known confounds (Murphy et al., 2009) during estimation of correlations. Analyzing seed ROI based FC after GSR yielded TTA-P2 engendered IHFC enhancements in only 22 ROIs after FDR correction (Supplementary Table 2). However, GSR did not affect the maps of TTAP2 vs. Baseline paired t-tests in FC to right S1BF and right auditory cortices ROIs significantly (see Supplementary Figures 7, 8).
Ideally, the rats should be scanned when they are awake. However, this was not feasible within the constraints of this study. Slow rhythms are suppressed by TTA-P2 under even deep levels of anesthesia in the studies (David et al., 2013) validating the pharmacological model employed in this project. These studies employed xylazine and ketamine, whereas Dexmed and pancuronium bromide were employed in this project. At the low levels of anesthesia employed in this study, the two drug combinations should have similar effects. FC networks have been reliably obtained under both these anesthesia regimens in rats (Hutchison et al., 2010; Magnuson et al., 2014). However, Dexmed anesthesia best approximates natural sleep, in terms of its electrophysiological signatures and sensitivity to slow rhythms among different classes of anesthesia including isoflurane and ketamine (Musizza et al., 2007; Akeju and Brown, 2017). Furthermore, we could not employ isoflurane in this study as it can affect the action of the drug TTA-P2 (Joksovic et al., 2005).
The fMRI cortico-cortical FC decreases but is by and large intact in general going from wakefulness to the low levels of anesthesia employed in this project (Hutchison et al., 2014; Magnuson et al., 2014; Bettinardi et al., 2015). This is consistent with the results of this project. Similarly, arousal-dependent fMRI signal and slow rhythms, while reduced, are also present during light anesthesia and wakefulness. Hence, results of our project are translatable to common fMRI studies conducted while subjects are awake.
Finally, the sample size of this preliminary study was small, which reduced the effect size. Hence, the experiments will need to be repeated with larger sample sizes to confirm these results. All of these above limitations will be eliminated in future studies planned for this project.
Conclusion
The results indicate that the vigilance-dependent components of the rsfMRI signal (e.g., QPPs) reflect the dynamics of cortical slow rhythms. Suppression of slow rhythms reduces the strength of vigilance-dependent rsfMRI signals and enhances intrinsic FC derived via rsfMRI in canonical brain function networks. These results have profound implications to our understanding of neurophysiological basis of rsfMRI signals. Future work would include large sample sizes, and simultaneous EEG/optical imaging recordings to directly examine cortical slow rhythms, and intrathalamic administration of TTA-P2 to specifically target only those TTCCs that are part of thalamocortical slow wave generating unit.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by Emory University Institutional Animal Care and Use Committee.
Author contributions
VK conducted all the experiments in the project including imaging and data analysis. EM conducted some of the data analysis. W-JP provided critical advice in setting up the experiment. SK conceived the idea, developed the experiment and supervised VK, EM, and W-JP. KG conceived the idea, developed the experiment conducted data analysis and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
Research reported in this publication was supported by National Institute of Neurological Disorders and Stroke of the National Institutes of Health under award number 1R21NS122013-01A1. Additional support was provided by Departments of Radiology and Imaging Sciences, and Biomedical Engineering, Emory University. The authors would also like to thank Dr. Ronald Voll, Centers of Systems Imaging Core, Emory University for constructive input.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.909999/full#supplementary-material
References
Abbas, A., Belloy, M., Kashyap, A., Billings, J., Nezafati, M., Schumacher, E. H., et al. (2019). Quasi-periodic patterns contribute to functional connectivity in the brain. Neuroimage 191, 193–204. doi: 10.1016/j.neuroimage.2019.01.076
Akeju, O., and Brown, E. N. (2017). Neural oscillations demonstrate that general anesthesia and sedative states are neurophysiologically distinct from sleep. Curr. Opin. Neurobiol. 44, 178–185. doi: 10.1016/j.conb.2017.04.011
Amemiya, S., Takao, H., Hanaoka, S., and Ohtomo, K. (2016). Global and structured waves of rs-fMRI signal identified as putative propagation of spontaneous neural activity. Neuroimage 133, 331–340. doi: 10.1016/j.neuroimage.2016.03.033
Andersson, J. L. R., Graham, M. S., Drobnjak, I., Zhang, H., and Campbell, J. (2018). Susceptibility-induced distortion that varies due to motion: Correction in diffusion MR without acquiring additional data. Neuroimage 171, 277–295. doi: 10.1016/j.neuroimage.2017.12.040
Bedwell, S. A., Billett, E. E., Crofts, J. J., and Tinsley, C. J. (2014). The topology of connections between rat prefrontal, motor and sensory cortices. Front. Syst. Neurosci. 8:177. doi: 10.3389/fnsys.2014.00177
Belloy, M. E., Naeyaert, M., Abbas, A., Shah, D., Vanreusel, V., van Audekerke, J., et al. (2018). Dynamic resting state fMRI analysis in mice reveals a set of Quasi-Periodic Patterns and illustrates their relationship with the global signal. Neuroimage 180(Pt B), 463–484. doi: 10.1016/j.neuroimage.2018.01.075
Benjamini, Y. (2010). Simultaneous and selective inference: Current successes and future challenges. Biom. J. 52, 708–721. doi: 10.1002/bimj.200900299
Bettinardi, R. G., Tort-Colet, N., Ruiz-Mejias, M., Sanchez-Vives, M. V., and Deco, G. (2015). Gradual emergence of spontaneous correlated brain activity during fading of general anesthesia in rats: Evidences from fMRI and local field potentials. Neuroimage 114, 185–198. doi: 10.1016/j.neuroimage.2015.03.037
Bianciardi, M., Fukunaga, M., van Gelderen, P., Horovitz, S. G., Shmueli, K., de Zwart, J. A., et al. (2009). Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: A 7 T study. Magn. Reson. Imaging 27, 1019–1029. doi: 10.1016/j.mri.2009.02.004
Billings, J. C. W., and Keilholz, S. D. (2018). The not-so-global blood oxygen level-dependent signal. Brain Connect. 8, 121–128. doi: 10.1089/brain.2017.0517
Birn, R. M., Diamond, J. B., Smith, M. A., and Bandettini, P. A. (2006). Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31, 1536–1548. doi: 10.1016/j.neuroimage.2006.02.048
Biswal, B. B., Van Kylen, J., and Hyde, J. S. (1997). Simultaneous assessment of flow and BOLD signals in resting-state functional connectivity maps. NMR Biomed. 10, 165–170. doi: 10.1002/(sici)1099-1492(199706/08)10:4/5<165::aid-nbm454>3.0.co;2-7
Brier, L. M., Landsness, E. C., Snyder, A. Z., Wright, P. W., Baxter, G. A., Bauer, A. Q., et al. (2019). Separability of calcium slow waves and functional connectivity during wake, sleep, and anesthesia. Neurophotonics 6:035002. doi: 10.1117/1.NPh.6.3.035002
Buckner, R. L., Krienen, F. M., and Yeo, B. T. (2013). Opportunities and limitations of intrinsic functional connectivity MRI. Nat. Neurosci. 16, 832–837. doi: 10.1038/nn.3423
Bullmore, E., and Sporns, O. (2009). Complex brain networks: Graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Chan, A. W., Mohajerani, M. H., LeDue, J. M., Wang, Y. T., and Murphy, T. H. (2015). Mesoscale infraslow spontaneous membrane potential fluctuations recapitulate high-frequency activity cortical motifs. Nat. Commun. 6:7738. doi: 10.1038/ncomms8738
Chan, R. W., Leong, A. T. L., Ho, L. C., Gao, P. P., Wong, E. C., Dong, C. M., et al. (2017). Low-frequency hippocampal-cortical activity drives brain-wide resting-state functional MRI connectivity. Proc. Natl. Acad. Sci. U.S.A. 114, E6972–E6981. doi: 10.1073/pnas.1703309114
Chang, C., Cunningham, J. P., and Glover, G. H. (2009). Influence of heart rate on the BOLD signal: The cardiac response function. Neuroimage 44, 857–869. doi: 10.1016/j.neuroimage.2008.09.029
Chang, C., Leopold, D. A., Scholvinck, M. L., Mandelkow, H., Picchioni, D., Liu, X., et al. (2016). Tracking brain arousal fluctuations with fMRI. Proc. Natl. Acad. Sci. U.S.A. 113, 4518–4523. doi: 10.1073/pnas.1520613113
Chauvette, S., Crochet, S., Volgushev, M., and Timofeev, I. (2011). Properties of slow oscillation during slow-wave sleep and anesthesia in cats. J. Neurosci. 31, 14998–15008. doi: 10.1523/JNEUROSCI.2339-11.2011
Cox, R. W. (1996). AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29, 162–173.
Cox, R. W., Chen, G., Glen, D. R., Reynolds, R. C., and Taylor, P. A. (2017). FMRI clustering in AFNI: False-positive rates redux. Brain Connect. 7, 152–171. doi: 10.1089/brain.2016.0475
Crunelli, V., David, F., Leresche, N., and Lambert, R. C. (2014). Role for T-type Ca2+ channels in sleep waves. Pflugers Arch. 466, 735–745. doi: 10.1007/s00424-014-1477-3
Crunelli, V., David, F., Lorincz, M. L., and Hughes, S. W. (2015). The thalamocortical network as a single slow wave-generating unit. Curr. Opin. Neurobiol. 31, 72–80. doi: 10.1016/j.conb.2014.09.001
Crunelli, V., and Hughes, S. W. (2010). The slow (<1 Hz) rhythm of non-REM sleep: A dialogue between three cardinal oscillators. Nat. Neurosci. 13, 9–17. doi: 10.1038/nn.2445
David, F., Schmiedt, J. T., Taylor, H. L., Orban, G., Di Giovanni, G., Uebele, V. N., et al. (2013). Essential thalamic contribution to slow waves of natural sleep. J. Neurosci. 33, 19599–19610. doi: 10.1523/jneurosci.3169-13.2013
Dreyfus, F. M., Tscherter, A., Errington, A. C., Renger, J. J., Shin, H. S., Uebele, V. N., et al. (2010). Selective T-type calcium channel block in thalamic neurons reveals channel redundancy and physiological impact of I(T)window. J. Neurosci. 30, 99–109. doi: 10.1523/jneurosci.4305-09.2010
Esser, S. K., Hill, S., and Tononi, G. (2009). Breakdown of effective connectivity during slow wave sleep: Investigating the mechanism underlying a cortical gate using large-scale modeling. J. Neurophysiol. 102, 2096–2111. doi: 10.1152/jn.00059.2009
Falahpour, M., Nalci, A., and Liu, T. T. (2018). The effects of global signal regression on estimates of resting-state blood oxygen-level-dependent functional magnetic resonance imaging and electroencephalogram vigilance correlations. Brain Connect. 8, 618–627. doi: 10.1089/brain.2018.0645
Frostig, R. D., Xiong, Y., Chen-Bee, C. H., Kvasnak, E., and Stehberg, J. (2008). Large-scale organization of rat sensorimotor cortex based on a motif of large activation spreads. J. Neurosci. 28, 13274–13284. doi: 10.1523/JNEUROSCI.4074-08.2008
Gopinath, K., Krishnamurthy, V., and Sathian, K. (2018). Accounting for non-Gaussian sources of spatial correlation in parametric functional magnetic resonance imaging paradigms I: Revisiting cluster-based inferences. Brain Connect. 8, 1–9. doi: 10.1089/brain.2017.0521
Greicius, M. (2008). Resting-state functional connectivity in neuropsychiatric disorders. Curr. Opin. Neurol. 21, 424–430. doi: 10.1097/WCO.0b013e328306f2c5
Greicius, M. D., Kiviniemi, V., Tervonen, O., Vainionpaa, V., Alahuhta, S., Reiss, A. L., et al. (2008). Persistent default-mode network connectivity during light sedation. Hum. Brain Mapp. 29, 839–847. doi: 10.1002/hbm.20537
He, B. J., and Raichle, M. E. (2009). The fMRI signal, slow cortical potential and consciousness. Trends Cogn. Sci. 13, 302–309. doi: 10.1016/j.tics.2009.04.004
Hudetz, A. G. (2012). General anesthesia and human brain connectivity. Brain Connect. 2, 291–302. doi: 10.1089/brain.2012.0107
Hutchison, R. M., Hutchison, M., Manning, K. Y., Menon, R. S., and Everling, S. (2014). Isoflurane induces dose-dependent alterations in the cortical connectivity profiles and dynamic properties of the brain’s functional architecture. Hum. Brain Mapp. 35, 5754–5775. doi: 10.1002/hbm.22583
Hutchison, R. M., Mirsattari, S. M., Jones, C. K., Gati, J. S., and Leung, L. S. (2010). Functional networks in the anesthetized rat brain revealed by independent component analysis of resting-state FMRI. J. Neurophysiol. 103, 3398–3406. doi: 10.1152/jn.00141.2010
Jalilianhasanpour, R., Beheshtian, E., Sherbaf, G., Sahraian, S., and Sair, H. I. (2019). Functional connectivity in neurodegenerative disorders: Alzheimer’s disease and frontotemporal dementia. Top. Magn. Reson. Imaging 28, 317–324. doi: 10.1097/RMR.0000000000000223
Joksovic, P. M., Brimelow, B. C., Murbartian, J., Perez-Reyes, E., and Todorovic, S. M. (2005). Contrasting anesthetic sensitivities of T-type Ca2+ channels of reticular thalamic neurons and recombinant Ca(v)3.3 channels. Br. J. Pharmacol. 144, 59–70. doi: 10.1038/sj.bjp.0706020
Kraus, A., Geuze, E., Schmahl, C., Greffrath, W., Treede, R. D., Bohus, M., et al. (2009). Differentiation of pain ratings in combat-related posttraumatic stress disorder. Pain 143, 179–185. doi: 10.1016/j.pain.2008.12.018
Laird, A. R., Fox, P. M., Eickhoff, S. B., Turner, J. A., Ray, K. L., McKay, D. R., et al. (2011). Behavioral interpretations of intrinsic connectivity networks. J. Cogn. Neurosci. 23, 4022–4037. doi: 10.1162/jocn_a_00077
Larson-Prior, L. J., Zempel, J. M., Nolan, T. S., Prior, F. W., Snyder, A. Z., and Raichle, M. E. (2009). Cortical network functional connectivity in the descent to sleep. Proc. Natl. Acad. Sci. U.S.A. 106, 4489–4494. doi: 10.1073/pnas.0900924106
Leong, A. T., Chan, R. W., Gao, P. P., Chan, Y. S., Tsia, K. K., Yung, W. H., et al. (2016). Long-range projections coordinate distributed brain-wide neural activity with a specific spatiotemporal profile. Proc. Natl. Acad. Sci. U.S.A. 113, E8306–E8315. doi: 10.1073/pnas.1616361113
Leopold, D. A., Murayama, Y., and Logothetis, N. K. (2003). Very slow activity fluctuations in monkey visual cortex: Implications for functional brain imaging. Cereb. Cortex 13, 422–433. doi: 10.1093/cercor/13.4.422
Liu, T. T., and Falahpour, M. (2020). Vigilance Effects in Resting-State fMRI. Front. Neurosci. 14:321. doi: 10.3389/fnins.2020.00321
Liu, X., Yanagawa, T., Leopold, D. A., Chang, C., Ishida, H., Fujii, N., et al. (2015a). Arousal transitions in sleep, wakefulness, and anesthesia are characterized by an orderly sequence of cortical events. Neuroimage 116, 222–231. doi: 10.1016/j.neuroimage.2015.04.003
Liu, X., Yanagawa, T., Leopold, D. A., Fujii, N., and Duyn, J. H. (2015b). Robust Long-Range Coordination of Spontaneous Neural Activity in Waking, Sleep and Anesthesia. Cereb. Cortex 25, 2929–2938. doi: 10.1093/cercor/bhu089
Llinas, R. R., and Steriade, M. (2006). Bursting of thalamic neurons and states of vigilance. J. Neurophysiol. 95, 3297–3308. doi: 10.1152/jn.00166.2006
Lu, H., Wang, L., Rea, W. W., Brynildsen, J. K., Jaime, S., Zuo, Y., et al. (2016). Low- but not high-frequency LFP correlates with spontaneous BOLD fluctuations in rat whisker barrel cortex. Cereb. Cortex 26, 683–694. doi: 10.1093/cercor/bhu248
Magnuson, M. E., Thompson, G. J., Pan, W. J., and Keilholz, S. D. (2014). Time-dependent effects of isoflurane and dexmedetomidine on functional connectivity, spectral characteristics, and spatial distribution of spontaneous BOLD fluctuations. NMR Biomed. 27, 291–303. doi: 10.1002/nbm.3062
Majeed, W., Magnuson, M., Hasenkamp, W., Schwarb, H., Schumacher, E. H., Barsalou, L., et al. (2011). Spatiotemporal dynamics of low frequency BOLD fluctuations in rats and humans. Neuroimage 54, 1140–1150. doi: 10.1016/j.neuroimage.2010.08.030
Majeed, W., Magnuson, M., and Keilholz, S. D. (2009). Spatiotemporal dynamics of low frequency fluctuations in BOLD fMRI of the rat. J. Magn. Reson. Imaging 30, 384–393. doi: 10.1002/jmri.21848
Masicampo, M. L., Shan, H. Q., Xu, V., Speagle, M., and Godwin, D. W. (2018). Selective blockade of T-Type Ca2+ channels is protective against alcohol-withdrawal induced seizure and mortality. Alcohol Alcohol. 53, 526–531. doi: 10.1093/alcalc/agy042
Massimini, M., Huber, R., Ferrarelli, F., Hill, S., and Tononi, G. (2004). The sleep slow oscillation as a traveling wave. J. Neurosci. 24, 6862–6870. doi: 10.1523/JNEUROSCI.1318-04.2004
Matsui, T., Murakami, T., and Ohki, K. (2016). Transient neuronal coactivations embedded in globally propagating waves underlie resting-state functional connectivity. Proc. Natl. Acad. Sci. U.S.A. 113, 6556–6561. doi: 10.1073/pnas.1521299113
Mayhew, J. E., Askew, S., Zheng, Y., Porrill, J., Westby, G. W., Redgrave, P., et al. (1996). Cerebral vasomotion: A 0.1-Hz oscillation in reflected light imaging of neural activity. Neuroimage 4(3 Pt 1), 183–193. doi: 10.1006/nimg.1996.0069
McCafferty, C., David, F., Venzi, M., Di Giovanni, G., Orban, G., Ubele, V. N., et al. (2012). T-type calcium channels of cortical and thalamocortical neurons are necessary for absence seizures. Program No. 159.04, (New Orleans, LA: Society for Neuroscience).
McCormick, D. A., McGinley, M. J., and Salkoff, D. B. (2015). Brain state dependent activity in the cortex and thalamus. Curr. Opin. Neurobiol. 31, 133–140. doi: 10.1016/j.conb.2014.10.003
Mitra, A., Snyder, A. Z., Tagliazucchi, E., Laufs, H., and Raichle, M. E. (2015). Propagated infra-slow intrinsic brain activity reorganizes across wake and slow wave sleep. eLife 4:e10781. doi: 10.7554/eLife.10781
Mitra, P. P., Ogawa, S., Hu, X., and Ugurbil, K. (1997). The nature of spatiotemporal changes in cerebral hemodynamics as manifested in functional magnetic resonance imaging. Magn. Reson. Med. 37, 511–518. doi: 10.1002/mrm.1910370407
Murphy, K., Birn, R. M., Handwerker, D. A., Jones, T. B., and Bandettini, P. A. (2009). The impact of global signal regression on resting state correlations: Are anti-correlated networks introduced? Neuroimage 44, 893–905. doi: 10.1016/j.neuroimage.2008.09.036
Musizza, B., Stefanovska, A., McClintock, P. V., Palus, M., Petrovcic, J., Ribaric, S., et al. (2007). Interactions between cardiac, respiratory and EEG-delta oscillations in rats during anaesthesia. J. Physiol. 580(Pt 1), 315–326. doi: 10.1113/jphysiol.2006.126748
Nir, Y., Mukamel, R., Dinstein, I., Privman, E., Harel, M., Fisch, L., et al. (2008). Interhemispheric correlations of slow spontaneous neuronal fluctuations revealed in human sensory cortex. Nat. Neurosci. 11, 1100–1108. doi: 10.1038/nn.2177
Nir, Y., Staba, R. J., Andrillon, T., Vyazovskiy, V. V., Cirelli, C., Fried, I., et al. (2011). Regional slow waves and spindles in human sleep. Neuron 70, 153–169. doi: 10.1016/j.neuron.2011.02.043
Okun, M., Steinmetz, N. A., Lak, A., Dervinis, M., and Harris, K. D. (2019). Distinct structure of cortical population activity on fast and infraslow timescales. Cereb. Cortex 29, 2196–2210. doi: 10.1093/cercor/bhz023
Pan, W. J., Thompson, G. J., Magnuson, M. E., Jaeger, D., and Keilholz, S. (2013). Infraslow LFP correlates to resting-state fMRI BOLD signals. Neuroimage 74, 288–297. doi: 10.1016/j.neuroimage.2013.02.035
Pan, W. J., Zhou, L., Clavijo, G. P., Khalilzad-Sharghi, V., and Keilholz, S. D. (2021). Optimized single-loop coil with 3D-shaped design for simultaneous fMRI and optical imaging in rodent. Proc. ISMRM 29.
Papp, E. A., Leergaard, T. B., Calabrese, E., Johnson, G. A., and Bjaalie, J. G. (2014). Waxholm Space atlas of the Sprague Dawley rat brain. Neuroimage 97, 374–386. doi: 10.1016/j.neuroimage.2014.04.001
Paxinos, G., and Watson, C. (2014). Paxino’s and Watson’s the rat brain in stereotaxic coordinates. Amsterdam: Elsevier.
Picchioni, D., Pixa, M. L., Fukunaga, M., Carr, W. S., Horovitz, S. G., Braun, A. R., et al. (2014). Decreased connectivity between the thalamus and the neocortex during human nonrapid eye movement sleep. Sleep 37, 387–397. doi: 10.5665/sleep.3422
Ponce-Alvarez, A., Deco, G., Hagmann, P., Romani, G. L., Mantini, D., and Corbetta, M. (2015). Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity. PLoS Comput. Biol. 11:e1004100. doi: 10.1371/journal.pcbi.1004100
Scholvinck, M. L., Maier, A., Ye, F. Q., Duyn, J. H., and Leopold, D. A. (2010). Neural basis of global resting-state fMRI activity. Proc. Natl. Acad. Sci. U.S.A. 107, 10238–10243. doi: 10.1073/pnas.0913110107
Schormans, A. L., Scott, K. E., Vo, A. M., Tyker, A., Typlt, M., Stolzberg, D., et al. (2016). Audiovisual temporal processing and synchrony perception in the rat. Front. Behav. Neurosci. 10:246. doi: 10.3389/fnbeh.2016.00246
Schwalm, M., Schmid, F., Wachsmuth, L., Backhaus, H., Kronfeld, A., Aedo Jury, F., et al. (2017). Cortex-wide BOLD fMRI activity reflects locally-recorded slow oscillation-associated calcium waves. eLife 6:e27602. doi: 10.7554/eLife.27602
Shen, K., Misic, B., Cipollini, B. N., Bezgin, G., Buschkuehl, M., Hutchison, R. M., et al. (2015). Stable long-range interhemispheric coordination is supported by direct anatomical projections. Proc. Natl. Acad. Sci. U.S.A. 112, 6473–6478. doi: 10.1073/pnas.1503436112
Sherman, S. M. (2001). Tonic and burst firing: Dual modes of thalamocortical relay. Trends Neurosci. 24, 122–126.
Sherman, S. M., Guillery, R. W., and Sherman, S. M. (2006). Exploring the thalamus and its role in cortical function. Cambridge, MA: MIT Press.
Sheroziya, M., and Timofeev, I. (2014). Global intracellular slow-wave dynamics of the thalamocortical system. J. Neurosci. 34, 8875–8893. doi: 10.1523/JNEUROSCI.4460-13.2014
Shipe, W. D., Barrow, J. C., Yang, Z. Q., Lindsley, C. W., Yang, F. V., Schlegel, K. A., et al. (2008). Design, synthesis, and evaluation of a novel 4-aminomethyl-4-fluoropiperidine as a T-type Ca2+ channel antagonist. J. Med. Chem. 51, 3692–3695. doi: 10.1021/jm800419w
Smith, P. H., Manning, K. A., and Uhlrich, D. J. (2010). Evaluation of inputs to rat primary auditory cortex from the suprageniculate nucleus and extrastriate visual cortex. J. Comp. Neurol. 518, 3679–3700. doi: 10.1002/cne.22411
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, (Suppl. 1), S208–S219. doi: 10.1016/j.neuroimage.2004.07.051
Stark, D. E., Margulies, D. S., Shehzad, Z. E., Reiss, P., Kelly, A. M., Uddin, L. Q., et al. (2008). Regional variation in interhemispheric coordination of intrinsic hemodynamic fluctuations. J. Neurosci. 28, 13754–13764. doi: 10.1523/JNEUROSCI.4544-08.2008
Steriade, M. (2000). Corticothalamic resonance, states of vigilance and mentation. Neuroscience 101, 243–276. doi: 10.1016/s0306-4522(00)00353-5
Steriade, M., McCormick, D. A., and Sejnowski, T. J. (1993a). Thalamocortical oscillations in the sleeping and aroused brain. Science 262, 679–685.
Steriade, M., Nunez, A., and Amzica, F. (1993b). A novel slow (< 1 Hz) oscillation of neocortical neurons in vivo: Depolarizing and hyperpolarizing components. J. Neurosci. 13, 3252–3265. doi: 10.1523/JNEUROSCI.13-08-03252.1993
Suzuki, M., and Larkum, M. E. (2020). General Anesthesia Decouples Cortical Pyramidal Neurons. Cell 180, 666–676.e13. doi: 10.1016/j.cell.2020.01.024
Tavor, I., Parker Jones, O., Mars, R. B., Smith, S. M., Behrens, T. E., and Jbabdi, S. (2016). Task-free MRI predicts individual differences in brain activity during task performance. Science 352, 216–220. doi: 10.1126/science.aad8127
Thomas, M. E., Lane, C. P., Chaudron, Y. M. J., Cisneros-Franco, J. M., and de Villers-Sidani, E. (2020). Modifying the adult rat tonotopic map with sound exposure produces frequency discrimination deficits that are recovered with training. J. Neurosci. 40, 2259–2268. doi: 10.1523/JNEUROSCI.1445-19.2019
Thompson, G. J., Magnuson, M. E., Merritt, M. D., Schwarb, H., Pan, W. J., McKinley, A., et al. (2013a). Short-time windows of correlation between large-scale functional brain networks predict vigilance intraindividually and interindividually. Hum. Brain Mapp. 34, 3280–3298. doi: 10.1002/hbm.22140
Thompson, G. J., Merritt, M. D., Pan, W. J., Magnuson, M. E., Grooms, J. K., Jaeger, D., et al. (2013b). Neural correlates of time-varying functional connectivity in the rat. Neuroimage 83, 826–836. doi: 10.1016/j.neuroimage.2013.07.036
Thompson, G. J., Pan, W. J., and Keilholz, S. D. (2015). Different dynamic resting state fMRI patterns are linked to different frequencies of neural activity. J. Neurophysiol. 114, 114–124. doi: 10.1152/jn.00235.2015
Thompson, G. J., Pan, W. J., Magnuson, M. E., Jaeger, D., and Keilholz, S. D. (2014). Quasi-periodic patterns (QPP): Large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity. Neuroimage 84, 1018–1031. doi: 10.1016/j.neuroimage.2013.09.029
Torao-Angosto, M., Manasanch, A., Mattia, M., and Sanchez-Vives, M. V. (2021). Up and down states during slow oscillations in slow-wave sleep and different levels of anesthesia. Front. Syst. Neurosci. 15:609645. doi: 10.3389/fnsys.2021.609645
Valdes-Hernandez, P. A., Sumiyoshi, A., Nonaka, H., Haga, R., Aubert-Vasquez, E., Ogawa, T., et al. (2011). An in vivo MRI template set for morphometry, tissue segmentation, and fMRI localization in rats. Front. Neuroinform. 5:26. doi: 10.3389/fninf.2011.00026
Volgushev, M., Chauvette, S., Mukovski, M., and Timofeev, I. (2006). Precise long-range synchronization of activity and silence in neocortical neurons during slow-wave oscillations [corrected]. J. Neurosci. 26, 5665–5672. doi: 10.1523/JNEUROSCI.0279-06.2006
Vyazovskiy, V. V., Olcese, U., Hanlon, E. C., Nir, Y., Cirelli, C., and Tononi, G. (2011). Local sleep in awake rats. Nature 472, 443–447. doi: 10.1038/nature10009
Wang, L., Saalmann, Y. B., Pinsk, M. A., Arcaro, M. J., and Kastner, S. (2012). Electrophysiological low-frequency coherence and cross-frequency coupling contribute to BOLD connectivity. Neuron 76, 1010–1020. doi: 10.1016/j.neuron.2012.09.033
Wang, X., Leong, A. T. L., Chan, R. W., Liu, Y., and Wu, E. X. (2019). Thalamic low frequency activity facilitates resting-state cortical interhemispheric MRI functional connectivity. Neuroimage 201:115985. doi: 10.1016/j.neuroimage.2019.06.063
Wen, H., and Liu, Z. (2016). Broadband electrophysiological dynamics contribute to global resting-state fMRI signal. J. Neurosci. 36, 6030–6040. doi: 10.1523/JNEUROSCI.0187-16.2016
Yousefi, B., Shin, J., Schumacher, E. H., and Keilholz, S. D. (2018). Quasi-periodic patterns of intrinsic brain activity in individuals and their relationship to global signal. Neuroimage 167, 297–308. doi: 10.1016/j.neuroimage.2017.11.043
Keywords: cortical slow wave activity, T-type calcium channel blocker, vigilance, arousal, animal studies, resting state fMRI, quasiperiodic patterns
Citation: Khalilzad Sharghi V, Maltbie EA, Pan W-J, Keilholz SD and Gopinath KS (2022) Selective blockade of rat brain T-type calcium channels provides insights on neurophysiological basis of arousal dependent resting state functional magnetic resonance imaging signals. Front. Neurosci. 16:909999. doi: 10.3389/fnins.2022.909999
Received: 31 March 2022; Accepted: 19 July 2022;
Published: 08 August 2022.
Edited by:
Yen-Yu Ian Shih, University of North Carolina at Chapel Hill, United StatesReviewed by:
Marcelo Febo, University of Florida, United StatesBasavaraju G. Sanganahalli, Yale University, United States
Copyright © 2022 Khalilzad Sharghi, Maltbie, Pan, Keilholz and Gopinath. 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.
*Correspondence: Kaundinya S. Gopinath, kaundinya.s.gopinath@emory.edu