Temporal stability of fMRI in medetomidine-anesthetized rats

Medetomidine has become a popular choice for anesthetizing rats during long-lasting sessions of blood-oxygen-level dependent (BOLD) functional magnetic resonance imaging (fMRI). Despite this, it has not yet been thoroughly established how commonly reported fMRI readouts evolve over several hours of medetomidine anesthesia and how they are affected by the precise timing, dose, and route of administration. We used four different protocols of medetomidine administration to anesthetize rats for up to six hours and repeatedly evaluated somatosensory stimulus-evoked BOLD responses and resting state functional connectivity. We found that the temporal evolution of fMRI readouts strongly depended on the method of administration. Intravenous administration of a medetomidine bolus (0.05 mg/kg), combined with a subsequent continuous infusion (0.1 mg/kg/h), led to temporally stable measures of stimulus-evoked activity and functional connectivity throughout the anesthesia. Deviating from the above protocol—by omitting the bolus, lowering the medetomidine dose, or using the subcutaneous route—compromised the stability of these measures in the initial two-hour period. We conclude that both an appropriate protocol of medetomidine administration and a suitable timing of fMRI experiments are crucial for obtaining consistent results. These factors should be considered for the design and interpretation of future rat fMRI studies.

(1 for active voxels, 0 elsewhere) and averaged across all EFS-fMRI runs to produce an activation probability map. The probability map revealed a consistently active cluster, anatomically corresponding to the left (contralateral to the stimulus) forelimb region of the primary somatosensory cortex-abbreviated as S1FL (Fig. 3a1). This cluster's center was active in 85.16% of all EFS-fMRI runs, whereas no other area was active in more than 7% of runs. To ensure that there was no systematic shift in the location of the active S1FL cluster, the 283 first-level activation maps were split into 12 groups according to the applied medetomidine protocol and the time-window since the start of medetomidine administration (early: 0-2 h; middle: 2-4 h; late: 4-6 h). Examination of activation probability maps from all groups (Fig. 3a2) verified that the active cluster's location remained stable across time and medetomidine protocols.
Shape and strength of stimulus-evoked responses. The portion of the S1FL that was significantly active in at least 30% of all 283 EFS-fMRI runs was taken as a functionally-defined region-of-interest (ROI). For each EFS-fMRI run, this ROI's mean BOLD time course was extracted, normalized to the pre-stimulus baseline,  Table 1), isoflurane was used to induce unconsciousness (5%) and during animal preparation (2-3%). For protocols 1, 2, and 4, a bolus of medetomidine was given after the preparation phase, followed by a gradual reduction of isoflurane and its eventual discontinuation 10 min later; continuous infusion of medetomidine commenced 15 min after the bolus. The bolus was omitted for protocol 3, with continuous infusion starting directly after preparation, and isoflurane being gradually reduced to zero over the course of 20-25 min. Anesthesia was maintained for a maximum of six hours since the start of medetomidine administration (time = 0). In the end, animals were provided with 2% isoflurane and freed from all equipment. Atipamezole was injected SC to antagonize medetomidine effects and to facilitate a smooth recovery. (b) Multiple fMRI runs were acquired per anesthesia session, with consecutive runs being alternated between somatosensory fMRI with electrical forepaw stimulation (EFS-fMRI) and resting state fMRI (RS-fMRI) with no stimulus. The stimulation paradigm applied during EFS-fMRI runs is shown in (c).
and averaged across stimulation blocks to produce an event-related average. This was used to extract the peak % signal change (peak ΔBOLD), as a measure of stimulus-evoked BOLD response strength (Fig. 3b). The examination of event-related averages revealed that the shape of S1FL BOLD responses was very similar across the four medetomidine protocols (Fig. 3c). The signal exhibited a sharp peak about 3 s after stimulus onset, followed by a plateau lasting till the end of the 30 s stimulation. We tested the dependence of peak ΔBOLD on time (early, middle and late two-hour-long time-windows), medetomidine protocol, and their interaction. The medetomidine protocol alone had no effect (ANOVA: F = 1.01, p = 0.38), but we found a strong time dependence (F = 37.65, p < 0.001) and a significant interaction between medetomidine protocol and time-window (F = 4.35, p < 0.001). Post-hoc pair-wise t-tests with p-value adjustment for multiple comparisons (Holm method) revealed varying temporal trends depending on the medetomidine protocol (Fig. 4a,b). Protocol 2 (IV with bolus) led to temporally stable BOLD responses (mean peak ΔBOLD of about 2.3% across time-windows), while the other three protocols led to a significant decrease in response strength between the early and the middle time-windows. Despite their differing temporal trends, all four protocols converged to a mean peak ΔBOLD of 2-3% after the initial two hours of medetomidine anesthesia (Fig. 4c). The results of all statistical comparisons can be found in Supplementary Tables S9-S11.
Resting state functional connectivity. RSFC was probed by examining the pair-wise correlations between the BOLD time courses of 28 anatomically defined ROIs (Fig. 5a). Examination of the pair-wise correlation matrices (Fig. 5b1), and of their network representations (Fig. 5c), showed that the hierarchical structure of the network, i.e. the strength of individual connections relative to each other, was consistent over time for all medetomidine protocols. The network's global RSFC-the mean correlation (Fisher's Z-score) across all unique ROI pairs-was computed for all RS-fMRI runs, and tested for its dependence on time (early, middle and late time-windows), medetomidine protocol, and their interaction. All three factors were found to have a significant effect on global RSFC, according to the results of ANOVA (medetomidine protocol: F = 7.24, p = 0.001; time-window: F = 5.15, p = 0.006; protocol-time interaction: F = 3.41, p = 0.003). Post-hoc pair-wise t-tests followed with p-value adjustment for multiple comparisons (Holm method). Bolus-based medetomidine administration (protocols 1, 2 and 4) was found to sustain stable global RSFC throughout the three time-windows. Omitting the bolus (protocol 3, IV no bolus) led to global RSFC values that were stronger in the early period (p = 0.002 compared to protocol 2; p < 0.001 compared to protocol 4) but decreased significantly between the early and the middle time-windows (p = 0.017). The results of all statistical comparisons can be found in Supplementary Tables S9,S10,S12.

Discussion
In the present study we evaluated the capacity of medetomidine-administered through four different protocols-to anesthetize rats for up to six hours and sustain temporally stable fMRI measures of stimulus-evoked activity and functional connectivity. We found that anesthesia duration exceeded five hours in most sessions An activation probability map produced by pooling significantly active clusters from first-level analyses of 283 EFS-fMRI runs. The asterisk marks the crossing of the anterior commissure (AC, −0.36 mm relative to the bregma, according to the Paxinos-Watson rat brain atlas). The rest of the slices are taken at 1 mm intervals from the AC slice. In (a2), the EFS-fMRI runs are grouped according to the applied medetomidine protocol and to the time-window since the start of medetomidine administration (early: 0-2 h; middle: 2-4 h; late: 4-6 h), to produce separate activation probability maps for each group. Only the slice containing the peak activation, 2 mm rostral to AC, is shown for each group. All maps in (a1,a2) are thresholded at 3% and overlaid on a T2-weighted structural study template. (b) The only consistently active cluster across all runs corresponds to the forelimb region of the left primary somatosensory cortex (S1FL). The location of this cluster is shown alongside the anatomical delineation of the same area from the Paxinos-Watson rat brain atlas. The functionally defined S1FL (area active in >30% of all EFS-fMRI runs) is set as a region-of-interest (ROI) for the extraction of BOLD signal time courses. Such a time course is shown for one example EFS-fMRI run, with the stimulation blocks marked by horizontal lines. Averaging the three stimulation blocks results in an event-related average, from which the peak % signal change (peak ΔBOLD) can be extracted. Event-related average responses (mean ± s.d. across rats) are plotted for all four medetomidine protocols (c). N: number of rats; n: number of fMRI runs. www.nature.com/scientificreports www.nature.com/scientificreports/ (35/48 or 72.9%; Fig. 2a), heart and respiratory rates were mostly stable after the first hour of anesthesia (Fig. 2b), and the majority of performed EFS-fMRI runs (241/283 or 85.16%) led to a significant and selective activation of the expected cortical area (Fig. 3a). However, both stimulus-evoked activity and functional connectivity exhibited varying temporal trends depending on the protocol of medetomidine administration (Figs 4 and 5).
Limited anesthesia duration has been referred to as a drawback of medetomidine protocols 34 , even though the exact duration is not reported in most papers. According to a study that has addressed this issue, rats anesthetized by constant IV infusion woke up spontaneously 3.5-4 hours into the experiment 25 . However, the rats in that study To calculate pair-wise RSFC, 28 regions-ofinterest (ROIs) were defined based on the Paxinos-Watson rat brain atlas-14 on each hemisphere. Pearson's correlations were calculated between the BOLD time courses of all unique ROI pairs and transformed into Fisher's Z-scores. A pair-wise correlation matrix is shown for one example RS-fMRI run; the mean correlation (Z-score) across all ROI pairs constitutes the global RSFC. The matrix is also represented as a weighted network graph, with the ROIs as nodes and their pairs as edges. For visualization clarity, only the strongest 30% of edges are shown; edge thickness and opacity scale linearly with the relative rank of the correlation value (the highest Z-score corresponds to the thickest edge); node radius scales with the weighted degree (weighted sum of edges passing through the node). (b) The 295 RS-fMRI runs are split into twelve groups according to the applied medetomidine protocol and the time-window since the start of medetomidine administration (early: 0-2 h; middle: 2-4 h; late: 4-6 h). (b1) For each group, the upper triangular matrix represents the mean (across runs) pair-wise correlation, while the lower triangular matrix shows the change in correlation compared to each protocol's early period. N: number of rats; n: number of fMRI runs. In (b2) global RSFC values are plotted across time-windows for each of the medetomidine protocols, while in (b3) they are plotted across medetomidine protocols for each of the time-windows. Dots represent single run values, bars show the estimated Least Squares Means, while error margins correspond to the estimated 95% confidence intervals. Asterisks indicate the significance of pair-wise t-tests, after adjusting for multiple comparisons (p < 0.05: *, p < 0.01: **, p < 0.001: ***). Global RSFC time courses for individual rats are provided in Supplementary  Fig. S7b. A complete list of Least Square Means, confidence intervals and pair-wise comparisons can be found in Supplementary Tables S10 and S12. (c) The mean correlation matrices (upper triangles of (b1)) are also visualized as network graphs, similarly to the example graph in (a). www.nature.com/scientificreports www.nature.com/scientificreports/ were tracheotomized and mechanically ventilated, while most rat fMRI experiments-including ours-are conducted in freely breathing animals. The anesthesia durations achieved in our experiments should satisfy the needs of most researchers. That said, the risk of spontaneous wake-ups must be taken into consideration. In our study, these wake-ups occurred more often in fMRI than in bench sessions (13/24 versus 8/24)-a possible effect of the scanner's acoustic noise-and were typically preceded by an increase in RR and arrhythmic breathing. A timely identification of such an RR increase allows the researcher to promptly supply isoflurane through the mask in order to mitigate an un-controlled wake-up, and the associated animal stress. If early wake-ups must be avoided altogether, other administration practices may be warranted, such as stepping up the infusion rate 25 , or adding a constant low dose of isoflurane throughout the imaging session 35 . The latter strategy has gained popularity in recent years, on the basis that it provides near-normal physiological conditions 35 and protects against epileptic seizures-which have been reported in animals with medetomidine-only anesthesia 36 . However, this might require mechanical ventilation, since even low doses of isoflurane, when combined with medetomidine, have been shown to suppress the amplitude of stimulus-evoked BOLD responses in spontaneously breathing rats 37 .
The recorded cardiorespiratory parameters (HR and RR) underwent rapid changes in the first hour of anesthesia but stabilized thereafter (Fig. 2b). The observed changes agree with medetomidine pharmacology and with previous rat fMRI studies 20,23,38 . The decrease in HR immediately after medetomidine injection can be attributed to its well-described α 2 -adrenergic effects. The activation of central nervous system α 2 -adrenoreceptors causes bradycardia due to a decrease in central sympathetic tone, while α 2 -adrenoreceptors in the peripheral vasculature mediate vasoconstriction. The latter mechanism leads to a transient elevation in blood pressure, which in turn exacerbates the bradycardia through a baroreceptor-mediated reflex 20 . The respiratory effects of medetomidine alone are considered to be minor 38 . However, medetomidine enhances the potency of isoflurane, and can thus amplify isoflurane-induced respiratory depression 20,39 . This is the most likely explanation for the observed sharp drop in RR right after medetomidine injection, as well as for its gradual recovery following isoflurane discontinuation. This effect necessitates the careful monitoring of animal physiology, especially in the transition period between isoflurane and medetomidine. In order to avoid severe respiratory depression, the speed with which isoflurane is faded out needs to be adapted to each individual animal.
The EFS-evoked activation of the contralateral S1FL (Fig. 3) is described by numerous other rat fMRI studies using the same stimulus [23][24][25] . For a subset of trials, these studies have also reported activations in the secondary somatosensory area and in sensory thalamic nuclei. If we increase statistical power by grouping multiple fMRI runs in a second-level analysis, we also find the same areas being responsive to EFS (see Supplementary Fig. S8). The shape of the EFS-evoked BOLD response in the S1FL (sharp peak at about 3 s after the onset of stimulation, followed by a plateau sustained till the end of EFS period; see Fig. 3c) is consistent with previous rat fMRI studies 24,40 , and the time-to-peak is close to the one reported in rats anesthetized with α-chloralose (2.49 ± 0.31 s) 41 . It is worth noting that this is considerably shorter than the 4-8 s delay commonly assumed in modelling the human hemodynamic response, which could be attributed to a species difference in vascular anatomy and dynamics 41 .
The temporal evolution of stimulus-evoked responses exhibited an interesting dependence on the medetomidine administration protocol (Fig. 4). All four protocols converged to steady-state peak ΔBOLD of 2-3% after the first two hours of medetomidine administration, but arrived there through different temporal trajectories. An IV bolus of 0.05 mg/kg followed by a continuous infusion of 0.1 mg/kg/h (protocol 2) had already achieved the steady-state since the early time-window. Deviating from the above dosage, either by omitting the bolus (protocol 3) or by downscaling both bolus and infusion doses (protocol 4), resulted in a time-dependent decrease from the early to the middle time-windows. This implies that early responses under protocols 3 and 4 were stronger than under protocol 2-an effect that appears in the data (Fig. 4c) but does not survive the multiple comparisons correction (see Supplementary Table S11). Omitting the bolus (protocol 3) also led to a time-dependent decrease in global RSFC: the overall pair-wise correlation strength dropped significantly between the early and middle time-windows. This time-dependent attenuation of BOLD readouts in sessions that lacked the bolus could reflect a negative relationship between the strength of these readouts and the concentration of medetomidine. Drug levels in the central nervous system might take hours to stabilize under the continuous infusion regime, given that the pharmacokinetics of medetomidine is characterized by a long terminal half-life of about 57 min, and a hysteresis between plasma and cerebrospinal fluid concentrations 38,42 . Bolus administration likely mediates a faster wash-in of the drug and an earlier establishment of the steady state. For this interpretation, we need to accept that medetomidine dose-dependently suppresses stimulus-evoked BOLD responses and RSFC, at least up to a certain level. Previous studies have in fact found such a dose-dependency for RSFC, but not for stimulus-evoked responses 40,43 . That said, the relevant experiments were restricted to a higher infusion rate range of 0.1-0.3 mg/kg/h (always preceded by a bolus), within which any effects on stimulus-evoked responses could have been saturated.
The case of protocol 1-the SC variant of protocol 2 and the most widely used administration scheme (see Supplementary Table 1)-is a peculiar one. Like protocols 3 and 4, it also exhibited a decrease in peak BOLD between the early and middle time-windows. However, a closer examination of individual rat trajectories (see Supplementary Fig. S7) reveals that this effect was driven by only two rats with exceptionally strong responses within the first 90 minutes; the other four rats produced temporally stable responses. This observation can be interpreted within the framework of dose-dependent response attenuation that was proposed above. Small variations in the injection site and/or the skin temperature may have delayed the drug's absorption into the blood-stream, leading to the observed stronger responses in a sub-set of rats. However, this interpretation cannot be verified without measuring blood concentrations of medetomidine, which we have not performed.
At this point, we would like to emphasize several limitations of the present study. We have no way of dissecting the vascular and neuronal contributions to the observed effects on BOLD readouts: we lack electrophysiological recordings of neural activity and important measures of vascular physiology, such as cerebral blood flow and arterial blood gas concentrations. We also lack fMRI data for very early time points: the earliest fMRI runs Scientific RepoRtS | (2019) 9:16673 | https://doi.org/10.1038/s41598-019-53144-y www.nature.com/scientificreports www.nature.com/scientificreports/ were acquired 30-50 minutes after the onset of medetomidine administration, due to the time spent for animal positioning, structural image acquisition and shimming, and the required overlap between medetomidine administration and isoflurane anesthesia. This delay was most pronounced for protocol 3, during which the said overlap was necessarily prolonged. To put our results into proper context, we would also like to point out the main methodological differences between our study and the existing rat fMRI literature. We used female rats, contrary to almost all other similar studies (Supplementary Table S1). Sex can affect drug metabolism through multiple potential mechanisms 44 , so the optimal dosing scheme may differ between males and females. Moreover, the rats in our study were lying supinely in the custom-built rat bed ( Supplementary Fig. S2), as opposed to the conventional prone position. In our experience, the supine position reduces respiration-related head motion, to a degree that allows us to avoid using ear bars (See Supplementary Fig. S6). The unconventional position and the absence of ear bars, which are a potentially painful stimulus, could have influenced physiological parameters and the duration of anesthesia.
Despite these limitations, our study does provide an empirical answer to the question of how BOLD readouts evolve over long-lasting medetomidine anesthesia sessions, allowing us to recommend best practices for fMRI studies in rats. Protocols 1 (SC with bolus) and 2 (IV with bolus), which follow the dosage of the originally published protocol 23 , are both well-suited for long-lasting fMRI sessions. The choice between them depends on the priorities of the researchers. Whenever maximizing the duration of the steady-state is deemed crucial, protocol 2 is the optimal choice. The main advantage of this protocol lies within the early two hours of anesthesia, since all four protocols converge to a similar steady-state after this time. This is however an important benefit: because the probability of spontaneous wake-ups increases over time, an earlier establishment of the steady-state can significantly prolong the useful experimental time. That said, we recommend always allowing a waiting period of 60 minutes between the bolus injection and the start of functional experiments. The first hour is characterized by rapid changes in cardio-respiratory physiology and is thus best reserved for structural image acquisition. If researchers prefer the practicality of the SC route, they should consider extending this waiting period to 90-120 minutes. We do not claim that these protocols (1 and 2) are necessarily the best choice for rat fMRI, since our study did not include all existing (see Supplementary Table S1) or possible administration protocols. It is however reasonable to consider them as the current 'default' choice, in light of the wealth of available data, and in the interest of promoting comparability among studies.

Experimental animals. All experiments followed the standards of the German Federal Law on Care and
Use of Laboratory Animals and were approved by the local government authorities (Lower Saxony State Office for Consumer Protection and Food Safety, approval number 33.19-42502-04-15/2042). A total number of 24 female adult Wistar rats (Charles Rivers Laboratories, Sulzfeld Germany) with a median body weight of 308 g (interquartile range 285-350 g) were used for this study. Rats were group-housed in cages with environmental enrichment, at a 12/12-hour light/dark cycle, with 20-24 °C temperature and 45-55% humidity. Water and standard chow were provided ad libitum. The 24 rats were split into four equally sized groups, each assigned to a different protocol of medetomidine administration (see Table 1). No animal was excluded from the experiments or from the analysis. The investigators were not blind to the group allocation.
Anesthesia and monitoring. Each animal was anesthetized on two sessions separated by a minimum of two weeks. The first session took place on a laboratory bench-top to accommodate unrestricted access to the animal and close monitoring of anesthesia duration and cardio-respiratory physiology (bench session). The functional imaging took place during the second session, performed inside a dedicated small animal MR system (fMRI session). All applied anesthetic protocols followed the same general outline, with isoflurane being used during preparation, and medetomidine during data acquisition (Fig. 1a).
Unconsciousness was induced in a chamber filled with 5% isoflurane and maintained throughout preparation with 2-3% isoflurane in medical air, supplied through a nose cone. The animal was placed in a supine position for both bench and fMRI sessions. The eyes were covered with ophthalmic ointment to prevent them from drying. A cannula was inserted in the SC tissue of the left flank (protocol 1) or in a tail vein (protocols 2-4). Two subdermal needle electrodes were placed in the right forepaw, between the 2 nd and the 4 th digit. Monitoring equipment was attached, consisting of a rectal temperature probe, a pneumatic pressure sensor placed on the chest, and three SC needle electrodes for electrocardiogram (ECG). After fixing the equipment with adhesive tape, the animal was transferred to a custom-built MRI-compatible rat bed.
Four different protocols of medetomidine administration were used: 1) SC with bolus; 2) IV with bolus; 3) IV no bolus; 4) IV lower dose. The detailed dosing for all protocols is given in Table 1. For protocols 1, 2, and 4, medetomidine (Dorbene vet, Zoetis Deutschland GmbH, Germany) was initially given as a bolus loading dose, followed by a gradual reduction of isoflurane and its eventual discontinuation 10 min later. To mitigate the risk of respiratory depression during this transition period, we tailored the speed of the gradual isoflurane reduction to each individual animal-aiming for a RR range of 40-70 bpm. Continuous infusion of medetomidine commenced 15 min after the bolus. The bolus was omitted for protocol 3, with continuous infusion starting directly after the preparation phase, and isoflurane being gradually reduced to zero over the course of 20-25 min. Attempts to shut off isoflurane earlier led to fast breathing, indicative of an imminent wake-up (bench sessions). All doses were delivered using an MRI-compatible infusion pump (PHD 2000 Infuse/Withdraw; Harvard apparatus, Holliston, Massachusetts, USA), which was loaded with the medetomidine solution diluted 1/20 in saline (final concentration of 0.05 mg/ml). From the start of medetomidine administration, heart rate (HR), respiratory rate (RR), and rectal temperature were monitored using the MR-compatible Model 1030 monitoring and gating system (Small Animal Instruments Inc., Stony Brook, NY 11790, USA). Rectal temperature was kept at 36.5 ± 1 °C using a pad heated by circulating water. Anesthesia was maintained for six hours since the start of medetomidine administration, except for sessions in which the rat spontaneously woke up earlier. For bench sessions, such wake-ups were identified as spontaneous movements of the rat and were found to be always preceded by the RR getting progressively faster and irregular. Since the animal is not visible during fMRI, the endpoint for fMRI sessions was set based on respiration (RR > 90/min, irregular, and continuously rising for at least 2 min). Anesthesia duration was defined as the time from the start of medetomidine administration (time = 0) till one of the aforementioned endpoints: spontaneous movement, rapid rise in RR, or passage of six hours. Upon reaching an endpoint, the animal was provided with 2% isoflurane through the nose cone and was disconnected from all electrodes, cannulas, and monitoring equipment. Finally, isoflurane was shut off and atipamezole (Atipazole, Prodivet pharmaceuticals, Belgium) was injected SC (0.25 mg/kg for protocols 1-3; 0.175 mg/kg for protocol 4) to facilitate a smooth wake-up.
The monitoring data (HR and RR traces) were recorded at a temporal resolution of 1 s and further processed with in-house python scripts as follows. Firstly, physiologically implausible values-corresponding to data acquisition errors-were dropped. Secondly, fMRI acquisition periods were removed from HR traces, since the rapidly switching magnetic gradients had introduced electrical noise in the ECG recording. Lastly, the traces were smoothed with an exponentially weighted moving average filter (smoothing factor α = 0.02). The processed HR and RR traces from bench and fMRI sessions closely resembled each other and were therefore pooled together. For each of the four medetomidine protocols the mean ( ± s.d.) HR/RR trace was calculated across all anesthesia sessions.
MRI acquisition. The fMRI sessions were performed inside a 9.4 Tesla Bruker BioSpec MR system, equipped with the BGA12 gradient, and operated via ParaVision 6.0.1 software. Signal was transmitted via a volume resonator (inner diameter 86 mm) and received by a rat brain 4-channel coil array (all equipment and software from Bruker BioSpin MRI GmbH, Ettlingen, Germany). Approximately 10 min after the discontinuation of isoflurane, the rat was positioned in the isocenter of the MR system, with its body lying supinely and its head fixed with the help of a bite bar. After obtaining low-resolution images for animal localization, a T2-weighted structural image was acquired using a TurboRARE sequence (repetition time 5.225 s, effective echo time 33 ms; 2 averages; RARE factor 8; 30-50 axial slices with a thickness of 0.5 mm; in-plane resolution 0.137 × 0.137 mm 2 ; matrix size 256 × 256). A field map was measured, and shims were adjusted to ensure homogeneity in an ellipsoidal volume encompassing the rat brain (MAPSHIM). This was followed by multiple 330-second-long BOLD fMRI runs, repeated approximately every 10 min, until the experimental endpoint was reached. All fMRI runs were acquired with a single-shot gradient-echo echo planar imaging sequence (220 repetitions; repetition time 1.5 s; echo time 15 ms; flip angle 90°; 30 axial slices with a thickness of 0.5 mm, ascending interleaved slice order with no slice gap; in-plane resolution 0.2 × 0.2 mm 2 ; matrix size 128 × 96; 4 dummy scans; bandwidth 375 kHz). The fMRI slices covered the entire rat brain, excluding the olfactory bulbs and the caudal 3/4 of the cerebellum. Consecutive fMRI runs were alternated between somatosensory fMRI with electrical forepaw stimulation (EFS-fMRI) and resting state fMRI (RS-fMRI) with no stimulus (Fig. 1b). This resulted in a total number of 283 EFS-fMRI and 295 RS-fMRI runs across all rats, spanning from 0.5 up to 6 hours since the start of medetomidine administration. EFS-fMRI runs included a baseline period of 60 s, followed by stimulation of the right forepaw in three 30 s blocks-each paired with 60 s of rest (Fig. 1c). Each stimulus block comprised square unipolar pulses with 3 mA amplitude and 0.3 ms pulse width, delivered at 9 Hz (Stimulus Generator 4002, Multi Channel Systems MCS GmbH, Reutlingen, Germany). The above parameters were chosen based on previous EFS studies in medetomidine-anesthetized rats 24,25,45 . MRI preprocessing. All MR images were first exported from ParaVision to DICOM format and then converted to NIfTI (Neuroimaging Informatics Technology Initiative; http://nifti.nimh.nih.gov) using the dcm2nii (https://www.nitrc.org/projects/dcm2nii/) tool. The structural T2-weighted images were used to construct a study template with the help of the Advanced Normalization Tools software-ANTs (http://stnava.github.io/ANTs/). The structural image of one of the rats was chosen as the target reference space. Every other structural image was registered to the target in two steps. First a linear rigid (3 translations and 3 rotations) registration was used to bring each structural image to the same origin and orientation as the reference image. This was followed by a non-linear symmetric diffeomorphic (SyN) registration to account for differences in brain size and shape 46 . All 24 structural images were averaged in the reference space to produce a mean anatomical image. A down-sampled (0.2 × 0.2 × 0.5 mm 3 ) version of this image was used to create a brain mask and served as the study template to which all functional datasets were eventually registered. The fMRI image series were preprocessed using functions from multiple neuroimaging toolkits, combined into a pipeline with python's Nipype library 47 . Images were corrected for slice timing with FSL (FMRIB Software Library, https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/) and temporally filtered with AFNI (https://afni.nimh.nih.gov/). A high-pass filter of 0.01 Hz was used to remove slow temporal drifts, with an additional low-pass filter of 0.15 Hz being applied to RS-fMRI datasets only. Spatial smoothing (FSL) was performed using a 0.5 mm 3D Gaussian kernel. We chose to skip motion correction and regression of nuisance variables, as these preprocessing steps were shown to have little effect in anesthetized and head-fixed rodents 48 . Besides, a quantification of motion parameters showed that motion was minimal across all fMRI runs (see Supplementary Fig. S6). A rigid transformation matrix was calculated between the mean image of each functional run and the native structural image. This matrix was combined with the previously calculated linear and non-linear transforms into a composite warp file, which was used to transform the preprocessed fMRI datasets into the study template space (ANTs). fMRI analysis and statistics. Areas activated during each EFS-fMRI run were identified via a first-level general linear model analysis carried out using FEAT (FMRI Expert Analysis Tool) Version 6.00, part of FSL 49 .
The stimulus time course (1 during stimulation, 0 elsewhere) was convolved with a standard double-gamma hemodynamic response function to generate the model predictor. The resulting statistical maps were masked for brain, thresholded non-parametrically using clusters determined by z > 3.1 and a corrected cluster significance threshold of p = 0.05 50,51 , binarized, and averaged across all 283 EFS-fMRI runs to construct an overall activation probability map. The region active in at least 30% of all EFS-fMRI runs (contralateral S1FL; see Fig. 3a) was taken as a region-of-interest (ROI). This ROI's mean BOLD time course was extracted from each EFS-fMRI run, normalized to the pre-stimulus baseline, and averaged across the three stimulation blocks to produce an event-related average. The peak % signal change (peak ΔBOLD) was extracted from the event-related average, as a measure of BOLD response strength (Fig. 3b). We also performed a second-level fixed-effects analysis (FEAT) for each medetomidine protocol, by pooling the corresponding EFS-fMRI runs and computing the mean group effect. The resulting statistical maps were masked for brain and thresholded non-parametrically using maximum height thresholding based on Gaussian random field theory, with a (corrected) significance threshold of P = 0.05 50,51 .
Resting state functional connectivity (RSFC) analysis was performed using in-house python scripts. All 295 preprocessed RS-fMRI datasets were normalized by subtracting the temporal mean and dividing by the standard deviation. Normalized BOLD signal time courses were extracted from 28 ROIs (14 on each hemisphere), manually delineated based on the Paxinos-Watson rat brain atlas 52 (see Fig. 5a, left). Pearson's correlation coefficients were calculated for all ROI pairs, transformed into Fisher's Z-scores, and stored as a pair-wise correlation matrix. The mean Z-score across all 378 unique ROI pairs served as a measure of global RSFC. To visualize the hierarchical structure of the functional connectome, each correlation matrix was also represented as a weighted graph, with the ROIs as nodes and the ROI pairs as edges. The edges were ranked by ascending correlation strength, and the rank was assigned as the numerical weight of the edge (strongest correlation to highest weight). Figure 5a (right) depicts a pair-wise correlation matrix and its graph representation for an example RS-fMRI run.
Statistical analysis was performed in R, version 3.6.1 (https://www.r-project.org/). We used the lmer function of the lme4 package 53 to design a linear mixed effects model, separately for each of two response variables: the EFS-evoked response strength (peak ΔBOLD) and the global RSFC (mean Z-score). This model can account for the repeated-measures design of the study and for missing data points (e.g. due to spontaneous earlier wake-ups). The fMRI runs were grouped into three two-hour-long time-windows (early: 0-2 h; middle: 2-4 h; late: 4-6 h), based on the onset time of their acquisition. The time-window, the medetomidine administration protocol, and their interaction were modeled as fixed effects, while individual rat intercepts were included as random effects. This approach treats the time-window as a within-subject repeated-measures factor and the medetomidine protocol as a between-subject factor. Model inference and statistical tests were carried out using the lmerTest package 54 . Specifically, Satterthwaite's method was used to create Type III Analysis of Variance (ANOVA) tables, and the ls_ means function was used to estimate Least Squares Means for all sub-classes (levels of protocol, time-window, and their combinations). Finally, the difflsmeans function was applied to estimate pair-wise differences between the sub-class means; the resulting p-values were adjusted for multiple comparisons according to the Holm method 55 . The ANOVA tables and the outputs of ls_means and difflsmeans functions can be found in the Supplementary Material (Supplementary Tables S9-S12).

Data availability
The unprocessed anatomical and functional MRI data are publicly available through the OpenNeuro repository (https://openneuro.org/datasets/ds001981/). They are provided in NIfTI format and organized according to the BIDS specification (Brain Imaging Data Structure, https://bids.neuroimaging.io/). The preprocessing scripts (python) can be shared upon request.