This protocol describes the quantification of volumetric cerebral hemodynamic variations in the mouse brain using functional ultrasound (fUS). Procedures for 3D functional activation map following sensory stimulation as well as resting-state functional connectivity are provided as illustrative examples, in anesthetized and awake mice.
Functional ultrasound (fUS) imaging is a novel brain imaging modality that relies on the high-sensitivity measure of the cerebral blood volume achieved by ultrafast doppler angiography. As brain perfusion is strongly linked to local neuronal activity, this technique allows the whole-brain 3D mapping of task-induced regional activation as well as resting-state functional connectivity, non-invasively, with unmatched spatio-temporal resolution and operational simplicity. In comparison with fMRI (functional magnetic resonance imaging), a main advantage of fUS imaging consists in enabling a complete compatibility with awake and behaving animal experiments. Moreover, fMRI brain mapping in mice, the most used preclinical model in Neuroscience, remains technically challenging due to the small size of the brain and the difficulty to maintain stable physiological conditions. Here we present a simple, reliable and robust protocol for whole-brain fUS imaging in anesthetized and awake mice using an off-the-shelf commercial fUS system with a motorized linear transducer, yielding significant cortical activation following sensory stimulation as well as reproducible 3D functional connectivity pattern for network identification.
Over the last two decades, neuroimaging has become an important tool for studying brain function and organization, enabling researchers to make important discoveries in the field of neuroscience. Today, functional magnetic resonance imaging (fMRI) has become the gold standard clinical neuroimaging technique to assess task or drug-evoked brain activation and to map functional connectivity at rest. While human fMRI has high reliability and sensibility, mouse fMRI remains technically challenging for numerous reasons1. First, fMRI has a poor spatial and temporal resolution. The small size of the mouse brain necessitates the use of strong magnetic fields using expensive scanners to achieve reasonable spatial resolution. Second, maintaining stable physiological parameters within the narrow range allowing efficient neuro-vascular coupling is very difficult in anesthetized mice. Finally, the blood oxygen level dependent (BOLD) signal on which fMRI studies rely has relatively poor sensitivity, leading to low signal-to-noise ratio when applied to mice and often requires repeated stimulus presentation over long acquisition to detect small variations. The mouse being the most widely used animal model in biomedical preclinical research, these limitations are partly responsible for the translational gap in neuropsychiatry, hindering new promising therapeutic targets on the bench to be transposed into effective treatments on bedside.
Functional ultrasound (fUS) is a recently developed neuroimaging technique based on ultrafast doppler2. By directly sampling cerebral blood volume, this technique allows probing brain activity in real-time through the neurovascular coupling. Compared to other neuroimaging techniques, fUS yields a spatial resolution of 100 µm and a temporal resolution in the tens of milliseconds. This technique allows whole-brain imaging of complete coronal sections of the mouse brain, completely non-invasively. Furthermore, it is fully compatible with conscious and behaving animals3,4,5. One of the main current limitations of fUS is its 2D feature, allowing to record a single coronal plane at the same time. While volumetric 3D fUS using 2D matrix array transducers has already been successfully demonstrated in rats6 and confirmed in mice7, its current lack of sensitivity requires a full craniotomy as well as averaging an important number of trials to detect a slight change of activity. Alternatively, linear transducers can be stepped across multiple positions and perform functional imaging plane by plane to cover the whole brain. However, this technique requires numerous experimental paradigm repetitions and as such long acquisition times (3-4 hours for the mouse brain)8,9.
In the present work, we describe a robust experimental platform including a commercially available functional ultrasound scanner and a fast plane-switching linear transducer with procedures to acquire 3D fUS data in anesthetized and awake mice, allowing volumetric and transcranial functional mapping of the mouse brain, non-invasively, without contrast-agent and within short acquisition times. We illustrate this feature by mapping somatosensory cortex activation following whisker stimulation as well as resting-state functional connectivity. Aside from animal preparation and data collection, we also describe the procedure for visualization, atlas registration and analysis of real-time fUS signals.
All the procedures presented here have been performed in agreement with the European Community Council Directive of 22 September 2010 (010/63/UE) and our local ethics committee (Comité d'éthique en matière d'expérimentation animale number 59, 'Paris Centre et Sud', project #2017-23). Adult mice (male C57BL/6 Rj, age 2-3 months, 20-30 g, from Janvier Labs, France) were housed 4 per cage with a 12h light/dark cycle, constant temperature at 22 °C and food and water ad libitum. Before the beginning of the experiments, animals are given a one-week minimum acclimatization period to housing conditions.
1. Animal preparation for anesthetized fUS imaging
Figure 1: Experimental setup for anesthetized fUS experiments. Description of the experimental setup showing all the scientific equipment needed during an anesthetized experiment. 1. Physiological monitoring : live display of both respiratory and cardiac frequencies. 2. Four-axis motor module (three translations and one rotation) monitored by Iconeus One system (9) and allowing to perform transcranial 3D tomographic scans or 4D acquisitions. 3a. Servo-Motor driving the whisker stimulator (3b.) The servo-motor is controlled by an arduino uno card which is interfaced with the Iconeus One system (9) in order to synchronize stimulation patterns with imaging sequences. 4.a. Syringe pump controller. 4.b. Syringe holder. 5.a. Temperature plate monitor controlling the heating plate. 5.b. Heating plate and rectal thermometer interfaced with the temperature plate monitor (5.a.). 6. Ultrasound gel placed between the animal's head and the ultrasound probe, providing acoustic coupling between them. 7. 15 MHz ultrasound probe. 8. Probe holder linking the probe (7) to the motor module (2). 9. Iconeus One equipment and software, allowing programing different imaging sequences and controlling the motors module (2) driving the probe (7). Please click here to view a larger version of this figure.
2. Animal preparation for awake head-fixed mice experiments
Figure 2: Experimental setup for awake fUS experiments. A. Schematic illustration of the headplate magnetic cover protecting the imaging window (created with BioRender.com). During imaging sessions (Left), the cover is removed to scan the brain in the large aperture offered by the head plate. B. Photograph of the experimental setup for transcranial awake imaging in head-fixed freely-behaving mice. 1. Iconeus One system and software, allowing to set up different imaging sequences and controlling the motors module. 2. Four-axis motors module (three translations and one rotation) monitored by Iconeus One system (1) and allowing 3D tomographic scans or 4D acquisitions. 3. Air dispensing table. 4. Mobile Home Cage (MHC). 5a,5b. Photographs showing closer views of the animal's environment inside the MHC. 6. Head fixation system clamping the head plate. 7. Probe holder linking the probe to the motor module (2). 8. 15 MHz ultrasonic probe. 9. Ultrasound gel placed between the mouse head and the ultrasound probe, providing acoustic coupling between them. 10. Servo-Motor driving the whisker stimulator. The Servo-Motor is controlled by an Arduino Uno card which is interfaced with the Iconeus One system through TTL signal (1) in order to synchronize stimulation patterns with imaging sequences. C. Illustration of the different spatial sampling possibilities (created with BioRender.com): in each case, the probe is stepped from the first position to the last one and a Doppler image is recorded at each position to reconstruct the stacked volume. This process is continuously repeated during the whole acquisition time. Dense Scan (left): the step between slices must be small enough (typically 400 µm, which corresponds to the elevation resolution) to allow volumetric imaging. Sparse Scan (right): if distant functional regions are targeted (at different positions), it is also possible to decrease spatial sampling to image different slices that intersect these regions while not compromising the temporal sampling. Please click here to view a larger version of this figure.
3. Probe positioning
4. Angiographic scan and atlas registration
5. Brain Positioning System (BPS)
Figure 3: Fast transcranial angiographic Scan and Brain registration for precise probe positioning. A. Schematic representation of the mouse brain being scanned transcranially by the ultrasonic probe from the first coronal slice (green) to the last coronal slice (blue) during a fast angiographic scan. The current imaged slice (represented in red) moves step by step from the back (green) to the front (blue) of the brain. Created with BioRender.com B. Screenshot of IcoScan acquisition software in the Angio 3D panel. The preset parameters on the right configure the fast scan. The positions in mm of the first slice, the last slice and the step size must be well chosen to scan linearly the whole brain. C. Screenshot of the IcoStudio processing software. The fast Angio 3D scan is automatically registered to a reference template of the mouse brain. The three-views (left) shows the superposition of the vasculature and the mouse brain Allen atlas in the coronal, sagittal and axial views. D. Linear lay-out (montage) of 16 slices (out of 31) from the 3D angio scan, with the registered Allen reference atlas superimposed onto the vasculature. E. Screenshot of the Brain Navigation panel showing the predicted imaging plane corresponding to the motor coordinates computed by the software thanks to the two markers placed at the center of the left and right primary somatosensory cortex, barrel fields region. Please click here to view a larger version of this figure.
6. Task-evoked experiment: whisker stimulation
7. 4D functional connectivity
This protocol describes the 3D quantification of cerebral hemodynamic variations transcranially in the mouse brain, at rest or in response to sensory stimulation. Whisker stimulation, a standard paradigm to map brain functional activation in rodents, has been selected as an example of sensory stimulation-evoked response. Figure 4 shows a representative activation map in response to mechanical whisker stimulation in an anesthetized mouse obtained using transcranial fUS imaging. The total trial time was 760 s, with a 60 s baseline (before and after the stimulation), an 80 s stimulation and a 60 s recovery time, repeated 5x. Significant activation was determined with the resolution of a general linear model (GLM) using a default mouse hemodynamic response function (HRF). The activated regions (Z scores with p-value >0.0000006 after stringent Bonferroni correction for multiple comparison) are displayed as color-coded values overlaid onto the Allen common coordinate framework template. Voxel-wise time course of the contralateral primary somatosensory cortex, barrel field region (S1BF) revealed a 15-20% increase of the CBV compared to baseline.
Figure 4: Transcranial activation Maps and rCBV time course following whiskers stimulation in ketamine/xylazine anesthetized mouse. A. Activation map showing significantly activated voxels following mechanical stimulation of the right whiskers (80 s ON, 60 s OFF ,5x) under ketamine/xylazine anesthesia. Maps were obtained by computing Z-scores based on general linear model analysis (GLM) with Bonferroni correction for multiple comparison. Z-scores (color coded) are overlaid on the Allen brain 3D template (after registration with the brain positioning system) and displayed in three-views: coronal (left), sagittal (middle) and axial (right). Anatomical regions from the Allen mouse brain common coordinate framework are displayed for reference. Activated voxels are well located inside the left S1BF cortex. Scale bar: 1 mm. Each sample volume was scanned over 2.8 mm (corresponding to 7 slices in the elevation direction) in 3.85 s allowing to record 20 volumic samples during each functional response. B. 3D rendering of whisker stimulation-evoked relative cerebral blood volume (rCBV) increase compared to baseline level. The anatomical delineation of the S1BF is indicated in blue. C. Time course of CBV variations in the left S1BF (blue) and the corresponding stimulus applied (red). Please click here to view a larger version of this figure.
The same paradigm has been applied in a head-fixed behaving mouse in the mobile homecage using the awake preset of IcoScan. Figure 5 presents the activation map after a multiple whisker stimulation experiment using the experimental setup described in Figure 2. A few posterior and caudal whiskers were stimulated with the following pattern: 30 s baseline followed by five consecutive trials of 30 s ON (4 Hz) and 30 s OFF (Figure 5C). Stimulation was delivered using a servo-motor driven by an Arduino UNO card triggering the image acquisition sequence for synchronization. Significant activation was determined with the resolution of a general linear model (GLM) using a default mouse hemodynamic response function (HRF). Multiple comparison correction was performed with the Bonferroni method. Conventional alpha level of 0.05 was normalized by the total number of voxels in the acquisition volume, resulting in a final stringent threshold of 0.000003.
Figure 5: Activation Maps and rCBV time course following whiskers stimulation in awake behaving mouse. A. Activation map showing significantly activated voxels following mechanical stimulation of the right whiskers (30 s ON, 30 s OFF, 5x) in an awake mouse in the mobile homecage. Maps were obtained by computing Z-scores based on general linear model analysis (GLM) with Bonferroni correction for multiple comparison (normalization by the total number of voxels). Z-scores (color coded) are overlaid on the Allen brain 3D template (after registration with the Brain Positioning System) and displayed in three-views: coronal (left), sagittal (middle) and axial (right). Anatomical regions from the Allen Mouse Brain Common Coordinate Framework are displayed for reference. Activated voxels are well located inside the left S1BF cortex. Scales bars, 1 mm. Each sample volume was scanned over 1.6 mm (corresponding to 3 slices in the elevation direction) in 3.85 s allowing to record 17 volumic samples during each functional response. B. 3D rendering of whisker stimulation-evoked relative Cerebral Blood Volume (rCBV) increase compared to baseline level. The anatomical delineation of the S1BF is indicated in blue. C. Illustration of the mouse in the mobile homecage during the right whisker stimulation experiment, during which five 30 s trials were performed for a total acquisition time of 330 s. D. Instantaneous relative CBV time course extracted inside the activated area (blue), with the corresponding stimulus superimposed (red). Please click here to view a larger version of this figure.
Figure 6 shows the temporal correlations of normalized low-frequency (<0.2 Hz) spontaneous CBV fluctuations between 3D brain regions (identified from registration to the Allen common coordinate framework) in a ketamine-xylazine anesthetized mouse. Total acquisition time was 20 min (1200 s). Atlas-supervised analysis revealed strong interhemispheric connectivity patterns, with resulting correlation coefficient values up to 0.8. Seed-based analysis in the dorsal hippocampus revealed a significant interhemispheric connectivity between right and left hippocampus as well as deep retro-hippocampal regions and piriform cortices. A seed region selected in the S1BF also resulted in a symmetrical (cortico-cortical) correlation pattern, as previously described.
Figure 6: Transcranial volumetric resting-state functional connectivity of the mouse brain under ketamine/xylazine anesthesia assessed on a 20 min 3D fUS acquisition. A. Correlation matrix based on 3D regions of the Allen common coordinate framework registered on the transcranial functional acquisition. The matrix is obtained by computing the normalized Pearson's correlation of spontaneous low-frequency fluctuations (<0.1 Hz) of the average time signals from all the voxels included in each identified ROI after slice timing correction. Each sampled volume was scanned over 1.6 mm in the elevation direction (corresponding to 4 slices) acquired over 2.2 s. B. Seed-based analysis projected onto a 3D template. The seed was selected within the right dorsal hippocampus at β – 2.1 mm. Correlation map is obtained by computing the Pearson Correlation coefficient between the temporal signals of the seed and each voxel of the whole acquisition after slice timing correction. C. 3D correlation map based on seed-based analysis with seed region selected within the S1BF at β – 2.1 mm. Scale bars: 1 mm. Please click here to view a larger version of this figure.
Whole brain imaging methods are crucial tools to better understand brain physiology and pathology. The method described here allows the precise quantification of hemodynamic signals in the living brain directly at the bench. The unmatched sensitivity and spatio-temporal resolution of functional ultrasound is particularly well suited for the mouse physiology. Functional responses and resting-state networks can be mapped within short acquisition times, longitudinally and without having to average trials or subjects to obtain a reliable measure. The relevant combination of high sensitivity ultrasonic linear probes and fast motorized setups enables one to perform transcranial volumetric fUS imaging in mice within reasonable acquisition times. This protocol can be performed either on anesthetized or awake mice using a mobile home cage.
Whisker stimulation, the sensory stimulus used as an illustrating example in this manuscript, is a standard functional activation paradigm in rodents and a reliable read-out to study sensory processing, neurovascular coupling and their alterations5,6,10, 11. While coarse manual brushing of the whiskers may be preferred for its ease of use, this method lacks spatial and temporal precision. The use of an automatic stimulator, such as the one described here triggered with the fUS imaging scanner, allows for a better control of several parameters including the time of onset, the amplitude displacement, the frequency as well as the angle of the Q-tip/comb, resulting in a better inter-animal reproducibility. Additionally, a more precise timing of stimulation enables the modeling of the Hemodynamic Response Function (HRF) by determining the time to onset and time to peak parameters12,13. To ensure better precision on the number of whiskers deflected during the stimulation (and thus the area of the activated region), more sophisticated stimulators can be adapted to this protocol. Many other stimuli such as light8, sound14 or odor presentation15 can be implemented using the same protocol.
The compatibility of functional ultrasound with awake and behaving animals is an important advantage compared to other neuroimaging techniques, enabling functional activation mapping without the anesthesia bias. Using an air-lifted mobile homecage is a good alternative to other existing head-fixed apparatus such as linear or spherical treadmills. While being firmly head-fixed, the motion of the homecage gives the mouse the illusion to navigate the environment, allowing a wide-range of behavioral testings to be coupled to fUS imaging16. However, the habituation procedure to head-fixing constitutes an important step to decrease stress, especially for experiments where it can be considered a confounding factor. The procedure detailed here (6-days of handling and habituation to head fixation) gives robust results for sensory stimulation and resting-state functional connectivity. However, it might be necessary to extend the habituation period for more refined behavioral tests17.
The authors have nothing to disclose.
This work was supported by the European Research Council (ERC) Advanced Grant N° 339244-FUSIMAGINE, the National Agency for Research funding 'Pinch' (ANR-18-CE37-005), the Inserm Research Technology Accelerator in Biomedical Ultrasound, the ElfUS technical core of the IPNP, Inserm U1266, the European research program FUSIMICE of the Human Brain Project, and EMBO Short-Term Fellowship 8439 to Andrea Kliewer.
BD Plastipak 1 mL syringes | Dutscher, France | 303172 | |
BD Microlance 26 Gauge needles | Dutscher, France | 303800 | |
Animal Temperature Controller (heating Plate coupled with a rectal probe) | Physitemp | TCAT-2DF | |
Arduino | Arduino | Arduino Uno-Rev3 | |
Atipamezole | Orion Pharma, France | Antisedan® | 5 mg/ml injectable solution |
Dental Ciment | Sun Médical, Shiga, japan | Superbond C&B | |
Depilatory cream | Klorane | N/A | |
Eye Ointment | TVM, UK | Ocry-gel | |
Hair trimmer | Wella Profesionnals | N/A | |
Head plates | Neurotar, Finland | Model 14 | |
Iconeus One standard package for fUS | Iconeus, France | Iconeus One | |
IcoScan acquisition software (v1.0) | Iconeus, France | IcoScan | |
IcoStudio analysis software (v1.0) | Iconeus, France | IcoStudio | |
Isoflurane Anesthesia station | Minerve, Esternay, France | ||
Ketamine | Virbac, France | Ketamine1000 | 100 mg/ml injectable solution |
Lidocaine | Vetoquinol | Lurocaine® | 20 mg/ml injectable solution |
Medetomidine | Orion Pharma, France | Domitor® | 1 mg/ml injectable solution |
Meloxicam | Boehringer lingelheim | Metacam® | 0.5 mg/ml injectable solution |
Mobile HomeCage Large with tracking capability | Neurotar, Finland | MHC-L-T-V4 | |
Monitoring of ECG and breathing rate | AD Systems, (USA) and LabChart software | ||
Servomotor | Feetech | FT90B | |
Stereotaxic frame | David Kopf (Tujunga, USA) | 900-WA | Using Mouse Adaptor (Ref: 922) and Non-Rupture Ear Bars (ref: 922) |
Surgical glue | 3M, USA | Vetbond | |
Syringe Pump | KD Scientific, USA | Legato® 130, Cat# 788130 | |
Ultrasound gel | DREXCO medical, France | Medi'Gel | |
Xylazine 2% | Bayer, France | Rompun® | 20 mg/ml injectable solution |