3D mapping of somatotopic reorganization with small animal functional MRI

Cortical Reorganization after Limb Loss: Bridging the Gap between Basic Science and Clinical Recovery

Despite the increasing incidence and prevalence of amputation across the globe, individuals with acquired limb loss continue to struggle with functional recovery and chronic pain. A more complete understanding of the motor and sensory remodeling of the peripheral and central nervous system that occurs postamputation may help advance clinical interventions to improve the quality of life for individuals with acquired limb loss. The purpose of this article is to first provide background clinical context on individuals with acquired limb loss and then to provide a comprehensive review of the known motor and sensory neural adaptations from both animal models and human clinical trials. Finally, the article bridges the gap between basic science researchers and clinicians that treat individuals with limb loss by explaining how current clinical treatments may restore function and modulate phantom limb pain using the underlying neural adaptations described above. This review should encourage the further development of novel treatments with known neurological targets to improve the recovery of individuals postamputation.Significance Statement In the United States, 1.6 million people live with limb loss; this number is expected to more than double by 2050. Improved surgical procedures enhance recovery, and new prosthetics and neural interfaces can replace missing limbs with those that communicate bidirectionally with the brain. These advances have been fairly successful, but still most patients experience persistent problems like phantom limb pain, and others discontinue prostheses instead of learning to use them daily. These problematic patient outcomes may be due in part to the lack of consensus among basic and clinical researchers regarding the plasticity mechanisms that occur in the brain after amputation injuries. Here we review results from clinical and animal model studies to bridge this clinical-basic science gap.

Thalamic activations in rat brain by fMRI during tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimulations

The thalamus is a crucial subcortical hub that impacts cortical activity. Tracing experiments in animals and post-mortem humans suggest rich morphological specificity of the thalamus. Very few studies reported rodent thalamic activations by functional MRI (fMRI) as compared to cortical activations for different sensory stimuli. Here, we show different portions of the rat thalamus in response to tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimuli with high field fMRI (11.7T) using a custom-build quadrature surface coil to capture high sensitivity signals from superficial and deep brain regions simultaneously. Results demonstrate reproducible thalamic activations during both tactile and non-tactile stimuli. Forepaw and whisker stimuli activated broader regions within the thalamus: ventral posterior lateral (VPL), ventral posterior medial (VPM), lateral posterior mediorostral (LPMR) and posterior medial (POm) thalamic nuclei. Visual stimuli activated dorsal lateral geniculate nucleus (DLG) of the thalamus but also parts of the superior/inferior colliculus, whereas olfactory stimuli activated specifically the mediodorsal nucleus of the thalamus (MDT). BOLD activations in LGN and MDT were much stronger than in VPL, VPM, LPMR and POm. These fMRI-based thalamic activations suggest that forepaw and whisker (i.e., tactile) stimuli engage VPL, VPM, LPMR and POm whereas visual and olfactory (i.e., non-tactile) stimuli, respectively, recruit DLG and MDT exclusively.

High spatiotemporal vessel-specific hemodynamic mapping with multi-echo single-vessel fMRI

High-resolution fMRI enables noninvasive mapping of the hemodynamic responses from individual penetrating vessels in animal brains. Here, a 2D multi-echo single-vessel fMRI (MESV-fMRI) method has been developed to map the fMRI signal from arterioles and venules with a 100 ms sampling rate at multiple echo times (TE, 3-30 ms) and short acquisition windows (<1 ms). The T₂*-weighted signal shows the increased extravascular effect on venule voxels as a function of TE. In contrast, the arteriole voxels show an increased fMRI signal with earlier onset than venules voxels at the short TE (3 ms) with increased blood inflow and volume effects. MESV-fMRI enables vessel-specific T2* mapping and presents T2*-based fMRI time courses with higher contrast-to-noise ratios (CNRs) than the T2*-weighted fMRI signal at a given TE. The vessel-specific T2* mapping also allows semi-quantitative estimation of the oxygen saturation levels (Y) and their changes (ΔY) at a given blood volume fraction upon neuronal activation. The MESV-fMRI method enables vessel-specific T2* measurements with high spatiotemporal resolution for better modeling of the fMRI signal based on the hemodynamic parameters.

Predicting the fMRI Signal Fluctuation with Recurrent Neural Networks Trained on Vascular Network Dynamics

Resting-state functional MRI (rs-fMRI) studies have revealed specific low-frequency hemodynamic signal fluctuations (<0.1 Hz) in the brain, which could be related to neuronal oscillations through the neurovascular coupling mechanism. Given the vascular origin of the fMRI signal, it remains challenging to separate the neural correlates of global rs-fMRI signal fluctuations from other confounding sources. However, the slow-oscillation detected from individual vessels by single-vessel fMRI presents strong correlation to neural oscillations. Here, we use recurrent neural networks (RNNs) to predict the future temporal evolution of the rs-fMRI slow oscillation from both rodent and human brains. The RNNs trained with vessel-specific rs-fMRI signals encode the unique brain oscillatory dynamic feature, presenting more effective prediction than the conventional autoregressive model. This RNN-based predictive modeling of rs-fMRI datasets from the Human Connectome Project (HCP) reveals brain state-specific characteristics, demonstrating an inverse relationship between the global rs-fMRI signal fluctuation with the internal default-mode network (DMN) correlation. The RNN prediction method presents a unique data-driven encoding scheme to specify potential brain state differences based on the global fMRI signal fluctuation, but not solely dependent on the global variance.

Effects of Hemodynamic Response Function Selection on Rat fMRI Statistical Analyses

The selection of the appropriate hemodynamic response function (HRF) for signal modeling in functional magnetic resonance imaging (fMRI) is important. Although the use of the boxcar-shaped hemodynamic response function (BHRF) and canonical hemodynamic response (CHRF) has gained increasing popularity in rodent fMRI studies, whether the selected HRF affects the results of rodent fMRI has not been fully elucidated. Here we investigated the signal change and t-statistic sensitivities of BHRF, CHRF, and impulse response function (IRF). The effect of HRF selection on different tasks was analyzed by using data collected from two groups of rats receiving either 3 mA whisker pad or 3 mA forepaw electrical stimulations (n = 10 for each group). Under whisker pad stimulation with large blood-oxygen-level dependent (BOLD) signal change (4.31 ± 0.42%), BHRF significantly underestimated signal changes (P < 0.001) and t-statistics (P < 0.001) compared with CHRF or IRF. CHRF and IRF did not provide significantly different t-statistics (P > 0.05). Under forepaw stimulation with small BOLD signal change (1.71 ± 0.34%), different HRFs provided insignificantly different t-statistics (P > 0.05). Therefore, the selected HRF can influence data analysis in rodent fMRI experiments with large BOLD responses but not in those with small BOLD responses.

Somatotopic Organization and Temporal Characteristics of Cerebrocortical Excitation in Response to Nasal Mucosa Stimulation With and Without an Odor in the Rat: An Optical Imaging Study

Nasal mucosa has roles in warming and humidifying inspired air and is highly sensitive to mechanical stimuli. Moreover, the upper part of the nasal mucosa expresses olfactory receptors processing olfactory information. Although the somatosensory map of the face in the primary (S1) and secondary (S2) somatosensory cortices is clearly documented, the map of the nasal mucosa and the effect of odors on their activities are largely unknown. This study aimed to identify the cortical regions in S1 and their temporal features in response to somatosensory stimulation of the nasal mucosa using an optical imaging technique in urethane-anesthetized rats. An air puff application response to nasal mucosa first occurred in a part of contralateral S1 and subsequently, spread toward the rostrally and ventrally adjacent sites. Upper pharynx stimulation initially activated this rostrally expanded site and the excitatory propagation from the initially activated region toward ventral region likely represented S2. Signal intensity and activated area increased dependent on air pressure. Nasal tip stimulation initially excited S1 region caudally adjacent to that of nasal mucosa. Moreover, the amplitude of S1 excitation was similar between air puff stimulation with and without an odor, amyl acetate. In contrast to contralateral S1, air puff stimulation with the odor showed a faint optical signal increase in the ipsilateral piriform cortex. These results suggest that somatosensory information from the nasal mucosa and skin, and upper pharynx are processed in spatially continuous regions of S1, and interaction between somatosensory and olfactory systems is relatively small in contralateral S1.

Ultra-Slow Single-Vessel BOLD and CBV-Based fMRI Spatiotemporal Dynamics and Their Correlation with Neuronal Intracellular Calcium Signals

Functional MRI has been used to map brain activity and functional connectivity based on the strength and temporal coherence of neurovascular-coupled hemodynamic signals. Here, single-vessel fMRI reveals vessel-specific correlation patterns in both rodents and humans. In anesthetized rats, fluctuations in the vessel-specific fMRI signal are correlated with the intracellular calcium signal measured in neighboring neurons. Further, the blood-oxygen-level-dependent (BOLD) signal from individual venules and the cerebral-blood-volume signal from individual arterioles show correlations at ultra-slow (<0.1 Hz), anesthetic-modulated rhythms. These data support a model that links neuronal activity to intrinsic oscillations in the cerebral vasculature, with a spatial correlation length of ∼2 mm for arterioles. In complementary data from awake human subjects, the BOLD signal is spatially correlated among sulcus veins and specified intracortical veins of the visual cortex at similar ultra-slow rhythms. These data support the use of fMRI to resolve functional connectivity at the level of single vessels.

Sensory and optogenetically driven single-vessel fMRI

Magnetic resonance imaging (MRI) sensitivity approaches vessel specificity. We developed a single-vessel functional MRI (fMRI) method to image the contribution of vascular components to blood oxygenation level-dependent (BOLD) and cerebral blood volume (CBV) fMRI signal. We mapped individual vessels penetrating the rat somatosensory cortex with 100-ms temporal resolution by MRI with sensory or optogenetic stimulation. The BOLD signal originated primarily from venules, and the CBV signal from arterioles. The single-vessel fMRI method and its combination with optogenetics provide a platform for mapping the hemodynamic signal through the neurovascular network with specificity at the level of individual arterioles and venules.

Optogenetic fMRI reveals distinct, frequency-dependent networks recruited by dorsal and intermediate hippocampus stimulations

Although the connectivity of hippocampal circuits has been extensively studied, the way in which these connections give rise to large-scale dynamic network activity remains unknown. Here, we used optogenetic fMRI to visualize the brain network dynamics evoked by different frequencies of stimulation of two distinct neuronal populations within dorsal and intermediate hippocampus. Stimulation of excitatory cells in intermediate hippocampus caused widespread cortical and subcortical recruitment at high frequencies, whereas stimulation in dorsal hippocampus led to activity primarily restricted to hippocampus across all frequencies tested. Sustained hippocampal responses evoked during high-frequency stimulation of either location predicted seizure-like afterdischarges in video-EEG experiments, while the widespread activation evoked by high-frequency stimulation of intermediate hippocampus predicted behavioral seizures. A negative BOLD signal observed in dentate gyrus during dorsal, but not intermediate, hippocampus stimulation is proposed to underlie the mechanism for these differences. Collectively, our results provide insight into the dynamic function of hippocampal networks and their role in seizures.

Interhemispheric plasticity protects the deafferented somatosensory cortex from functional takeover after nerve injury

Functional changes across brain hemispheres have been reported after unilateral cortical or peripheral nerve injury. Interhemispheric callosal connections usually underlie this cortico-cortical plasticity. However, the effect of the altered callosal inputs on local cortical plasticity in the adult brain is not well studied. Ipsilateral functional magnetic resonance imaging (fMRI) activation has been reliably detected in the deafferented barrel cortex (BC) at 2 weeks after unilateral infraorbital denervation (IO) in adult rats. The ipsilateral fMRI signal relies on callosal-mediated interhemispheric plasticity. This form of interhemispheric plasticity provides a good chronic model to study the interaction between callosal inputs and local cortical plasticity. The receptive field of forepaw in the primary somatosensory cortex (S1), which is adjacent to the BC, was mapped with fMRI. The S1 receptive field expanded to take over a portion of the BC in 2 weeks after both ascending inputs and callosal inputs were removed in IO rats with ablated contralateral BC (IO+ablation). This expansion, estimated specifically by fMRI mapping, is significantly larger than what has been observed in the IO rats with intact callosal connectivity, as well as in the rats with sham surgery. This work indicates that altered callosal inputs prevent the functional takeover of the deafferented BC from adjacent cortices and may help preserve the functional identity of the BC.

Chronic assessment of cerebral hemodynamics during rat forepaw electrical stimulation using functional ultrasound imaging

Functional ultrasound imaging is a method recently developed to assess brain activity via hemodynamics in rodents. Doppler ultrasound signals allow the measurement of cerebral blood volume (CBV) and red blood cells' (RBCs') velocity in small vessels. However, this technique originally requires performing a large craniotomy that limits its use to acute experiments only. Moreover, a detailed description of the hemodynamic changes that underlie functional ultrasound imaging has not been described but is essential for a better interpretation of neuroimaging data. To overcome the limitation of the craniotomy, we developed a dedicated thinned skull surgery for chronic imaging. This procedure did not induce brain inflammation nor neuronal death as confirmed by immunostaining. We successfully acquired both high-resolution images of the microvasculature and functional movies of the brain hemodynamics on the same animal at 0, 2, and 7 days without loss of quality. Then, we investigated the spatiotemporal evolution of the CBV hemodynamic response function (HRF) in response to sensory-evoked electrical stimulus (1 mA) ranging from 1 (200 μs) to 25 pulses (5s). Our results indicate that CBV HRF parameters such as the peak amplitude, the time to peak, the full width at half-maximum and the spatial extent of the activated area increase with stimulus duration. Functional ultrasound imaging was sensitive enough to detect hemodynamic responses evoked by only a single pulse stimulus. We also observed that the RBC velocity during activation could be separated in two distinct speed ranges with the fastest velocities located in the upper part of the cortex and slower velocities in deeper layers. For the first time, functional ultrasound imaging demonstrates its potential to image brain activity chronically in small animals and offers new insights into the spatiotemporal evolution of cerebral hemodynamics.

Repeated BOLD-fMRI imaging of deep brain stimulation responses in rats

Functional magnetic resonance imaging (fMRI) provides a picture of the global spatial activation pattern of the brain. Interest is growing regarding the application of fMRI to rodent models to investigate adult brain plasticity. To date, most rodent studies used an electrical forepaw stimulation model to acquire fMRI data, with α-chloralose as the anesthetic. However, α-chloralose is harmful to animals, and not suitable for longitudinal studies. Moreover, peripheral stimulation models enable only a limited number of brain regions to be studied. Processing between peripheral regions and the brain is multisynaptic, and renders interpretation difficult and uncertain. In the present study, we combined the medetomidine-based fMRI protocol (a noninvasive rodent fMRI protocol) with chronic implantation of an MRI-compatible stimulation electrode in the ventroposterior (VP) thalamus to repetitively sample thalamocortical responses in the rat brain. Using this model, we scanned the forebrain responses evoked by the VP stimulation repeatedly of individual rats over 1 week. Cortical BOLD responses were compared between the 2 profiles obtained at day1 and day8. We discovered reproducible frequency- and amplitude-dependent BOLD responses in the ipsilateral somatosensory cortex (S1). The S1 BOLD responses during the 2 sessions were conserved in maximal response amplitude, area size (size ratio from 0.88 to 0.91), and location (overlap ratio from 0.61 to 0.67). The present study provides a long-term chronic brain stimulation protocol for studying the plasticity of specific neural circuits in the rodent brain by BOLD-fMRI.

Optogenetic patterning of whisker-barrel cortical system in transgenic rat expressing channelrhodopsin-2

The rodent whisker-barrel system has been an ideal model for studying somatosensory representations in the cortex. However, it remains a challenge to experimentally stimulate whiskers with a given pattern under spatiotemporal precision. Recently the optogenetic manipulation of neuronal activity has made possible the analysis of the neuronal network with precise spatiotemporal resolution. Here we identified the selective expression of channelrhodopsin-2 (ChR2), an algal light-driven cation channel, in the large mechanoreceptive neurons in the trigeminal ganglion (TG) as well as their peripheral nerve endings innervating the whisker follicles of a transgenic rat. The spatiotemporal pattern of whisker irradiation thus produced a barrel-cortical response with a specific spatiotemporal pattern as evidenced by electrophysiological and functional MRI (fMRI) studies. Our methods of generating an optogenetic tactile pattern (OTP) can be expected to facilitate studies on how the spatiotemporal pattern of touch is represented in the somatosensory cortex, as Hubel and Wiesel did in the visual cortex.

Monitoring the stress-level of rats with different types of anesthesia: a tail-artery cannulation protocol

INTRODUCTION

Functional MRI in rats under anesthesia can largely minimize motion artifacts and attenuate the stress of the animal. However, two issues remain to be clarified and improved. First, fMRI results obtained with different types of anesthesia during surgical preparation and imaging show a large variability, which could be caused by the variable stress level of the rodents. Second, the most common surgical procedure used for anesthesia, blood gas analysis and mean arterial blood-pressure (MABP) monitoring is the femoral vein and artery catheterization that makes longitudinal studies difficult.

METHODS

In order to examine the variability of the stress level with three different anesthesia protocols using isoflurane (Iso), medetomidine-ketamine (MK) or propofol-remifentanil (PR), we measured the plasma corticosterone (CORT) concentration with (125)I-radioimmunoassay in blood samples collected prior to, immediately after and 60min after surgery. Tail-artery and vein catheterization was adapted for long-term monitoring of MABP with periodic blood sampling and is proposed as a less invasive and technically simple alternative to femoral vessel catheterization in fMRI preparation protocols.

RESULTS

We show that the CORT concentration depends on the anesthesia protocol with both alternatives providing more efficient stress reduction than the protocol using Iso. However, only the protocol using PR achieved a significant hormone reduction during surgery. Stress was not reliably manifested in changes in heart-rate and breathing-rate. Anesthesia and strain related changes in these two physiological parameters may be assigned to the pharmacological effects of the premedication and anesthetic agents. The results indicate also that MABP can be monitored over a long period of time (e.g. functional imaging session) through an arterial access point in the rat tail after cannulation with the proposed procedure.

DISCUSSION AND CONCLUSION

Animals can experience stress during fMRI preparation protocols without obvious signs in commonly monitored physiological parameters. Our results challenge the efficiency of surgical protocols using Iso as mono-anesthetic agent, even when extended with topical analgesia. It was demonstrated that the CORT-based stress-level measurement through tail-artery cannulation can be used for developing anesthesia protocols (i.e. the presented PR protocol) when setting up future fMRI studies. The proposed surgical method for the tail is expected to facilitate longitudinal fMRI studies with permanent arterial access.

Rapid functional reorganization of the forelimb cortical representation after thoracic spinal cord injury in adult rats

Thoracic spinal cord injured rats rely largely on forelimbs to walk, as their hindlimbs are dysfunctional. This increased limb use is accompanied by expansion of the cortical forelimb sensory representation. It is unclear how quickly the representational changes occur and whether they are at all related to the behavioral adaptation. Using blood oxygenation level dependent functional mangetic resonance imaging (BOLD-fMRI) we show that major plastic changes of the somato-sensory map can occur as early as one day after injury. The extent of map increase was variable between animals, and some animals showed a reduction in map size. However, at three or seven days after injury a significant enhancement of the forelimb representation was evident in all the animals. In a behavioral test for precise limb control, crossing of a horizontal ladder, the injured rats relied almost entirely on their forelimbs; they initially made more mistakes than at 7 days post injury. Remarkably, in the individual animals the behavioral performance seen at seven days was proportional to the physiological change present at one day after injury. The rapid increase in cortical representation of the injury-spared body part may provide the additional neural substrate necessary for high level behavioral adaptation.

Genetic tools to manipulate MRI contrast

Advances in molecular biology in the early 1970s revolutionized research strategies for the study of complex biological processes, which, in turn, created a high demand for new means to visualize these dynamic biological changes noninvasively and in real time. In this respect, MRI was a perfect fit, because of the versatile possibility to alter the different contrast mechanisms. Genetic manipulations are now being translated to MRI through the development of reporters and sensors, as well as the imaging of transgenic and knockout mice. In the past few years, a new molecular biology toolset, namely optogenetics, has emerged, which allows for the manipulation of cellular behavior using light. This technology provides a few particularly attractive features for combination with newly developed MRI techniques for the probing of in vivo cellular and, in particular, neural processes, specifically the ability to control focal, genetically defined cellular populations with high temporal resolution using equipment that is magnetically inert and does not interact with radiofrequency pulses. Recent studies have demonstrated that the combination of optogenetics and functional MRI (fMRI) can provide an appropriate platform to investigate in vivo, at the cellular and molecular levels, the neuronal basis of fMRI signals. In addition, this novel combination of optogenetics with fMRI has the potential to resolve pre-synaptic versus post-synaptic changes in neuronal activity and changes in the activity of large neuronal networks in the context of plasticity associated with development, learning and pathophysiology.

Large-scale reorganization of the somatosensory cortex of adult macaque monkeys revealed by fMRI

Somatosensory cortex of adult primates undergoes topographic reorganization following spinal cord or peripheral nerve injuries. Electrophysiological studies in monkeys show that after chronic lesions of dorsal columns of the spinal cord at cervical levels, there is an expansion of face representation into the deafferented hand region of area 3b of cortex. However, these techniques can sample only a limited portion of the brain. In order to help understand mechanisms of brain reorganization use of noninvasive tools in non-human primate experimental model is important. Use of blood oxygen level dependent-functional magnetic resonance imaging (BOLD-fMRI) to study brain reorganization in non-human primates has been extremely limited. Here, we show that in monkeys with long-term unilateral lesions of the dorsal columns at cervical levels, tactile stimulation of the chin showed BOLD activation in the deafferented hand region of contralesional area 3b in the post-central gyrus. In a monkey with a partial lesion of the dorsal columns, stimulations of both hand and chin activated the partially deafferented hand region. We also show that the somatotopic organization in the non-deafferented ipsilesional somatosensory cortex remained normal.

Functional ultrasound imaging of the brain: theory and basic principles

Hemodynamic changes in the brain are often used as surrogates of neuronal activity to infer the loci of brain activity. A major limitation of conventional Doppler ultrasound for the imaging of these changes is that it is not sensitive enough to detect the blood flow in small vessels where the major part of the hemodynamic response occurs. Here, we present a μDoppler ultrasound method able to detect and map the cerebral blood volume (CBV) over the entire brain with an important increase in sensitivity. This method is based on imaging the brain at an ultrafast frame rate (1 kHz) using compounded plane wave emissions. A theoretical model demonstrates that the gain in sensitivity of the μDoppler method is due to the combination of 1) the high signal-to-noise ratio of the gray scale images, resulting from the synthetic compounding of backscattered echoes; and 2) the extensive signal averaging enabled by the high temporal sampling of ultrafast frame rates. This μDoppler imaging is performed in vivo on trepanned rats without the use of contrast agents. The resulting images reveal detailed maps of the rat brain vascularization with an acquisition time as short as 320 ms per slice. This new method is the basis for a real-time functional ultrasound (fUS) imaging of the brain.

Thalamocortical inputs show post-critical-period plasticity

Experience-dependent plasticity in the adult brain has clinical potential for functional rehabilitation following central and peripheral nerve injuries. Here, plasticity induced by unilateral infraorbital (IO) nerve resection in 4-week-old rats was mapped using MRI and synaptic mechanisms were elucidated by slice electrophysiology. Functional MRI demonstrates a cortical potentiation compared to thalamus 2 weeks after IO nerve resection. Tracing thalamocortical (TC) projections with manganese-enhanced MRI revealed circuit changes in the spared layer 4 (L4) barrel cortex. Brain slice electrophysiology revealed TC input strengthening onto L4 stellate cells due to an increase in postsynaptic strength and the number of functional synapses. This work shows that the TC input is a site for robust plasticity after the end of the previously defined critical period for this input. Thus, TC inputs may represent a major site for adult plasticity in contrast to the consensus that adult plasticity mainly occurs at cortico-cortical connections.

Imaging the spatio-temporal dynamics of supragranular activity in the rat somatosensory cortex in response to stimulation of the paws

We employed voltage-sensitive dye (VSD) imaging to investigate the spatio-temporal dynamics of the responses of the supragranular somatosensory cortex to stimulation of the four paws in urethane-anesthetized rats. We obtained the following main results. (1) Stimulation of the contralateral forepaw evoked VSD responses with greater amplitude and smaller latency than stimulation of the contralateral hindpaw, and ipsilateral VSD responses had a lower amplitude and greater latency than contralateral responses. (2) While the contralateral stimulation initially activated only one focus, the ipsilateral stimulation initially activated two foci: one focus was typically medial to the focus activated by contralateral stimulation and was stereotaxically localized in the motor cortex; the other focus was typically posterior to the focus activated by contralateral stimulation and was stereotaxically localized in the somatosensory cortex. (3) Forepaw and hindpaw somatosensory stimuli activated large areas of the sensorimotor cortex, well beyond the forepaw and hindpaw somatosensory areas of classical somatotopic maps, and forepaw stimuli activated larger cortical areas with greater activation velocity than hindpaw stimuli. (4) Stimulation of the forepaw and hindpaw evoked different cortical activation dynamics: forepaw responses displayed a clear medial directionality, whereas hindpaw responses were much more uniform in all directions. In conclusion, this work offers a complete spatio-temporal map of the supragranular VSD cortical activation in response to stimulation of the paws, showing important somatotopic differences between contralateral and ipsilateral maps as well as differences in the spatio-temporal activation dynamics in response to forepaw and hindpaw stimuli.

Structure-function relationship of cerebral networks in experimental neuroscience: contribution of magnetic resonance imaging

The analysis of neuronal networks, their interactions in resting condition as well as during brain activation have become of great interest for a better understanding of the signal processing of the brain during sensory stimulus or cognitive tasks. Parallel to the study of the functional networks and their dynamics, the underlying network structure is highly important as it provides the basis of the functional interaction. Moreover, under pathological conditions, some nodes in such a net may be impaired and the function of the whole network affected. Mechanisms such as functional deficit and improvement, and plastic reorganization are increasingly discussed in the context of existing structural and functional networks. While many of these aspects have been followed in human and clinical studies, the experimental range is limited for obvious reasons. Here, animal experimental studies are needed as they permit longer scan times and, moreover, comparison with invasive histology. Experimental non-invasive imaging modalities are now able to perform impressive contributions. In this review we try to highlight most recent new cutting-edge developments and applications in experimental neuroscience of functional and structural networks of the brain, relying on non-invasive imaging. We focus primarily on the potential of experimental Magnetic Resonance Imaging (MRI), but also touch upon micro positron emission tomography (μPET) and optical imaging developments where they are applicable to the topic of the present review.

BOLD fMRI investigation of the rat auditory pathway and tonotopic organization

Rodents share general anatomical, physiological and behavioral features in the central auditory system with humans. In this study, monaural broadband noise and pure tone sounds are presented to normal rats and the resulting hemodynamic responses are measured with blood oxygenation level-dependent (BOLD) fMRI using a standard spin-echo echo planar imaging sequence (without sparse temporal sampling). The cochlear nucleus (CN), superior olivary complex, lateral lemniscus, inferior colliculus (IC), medial geniculate body and primary auditory cortex, all major auditory structures, are activated by broadband stimulation. The CN and IC BOLD signal changes increase monotonically with sound pressure level. Pure tone stimulation with three distinct frequencies (7, 20 and 40 kHz) reveals the tonotopic organization of the IC. The activated regions shift from dorsolateral to ventromedial IC with increasing frequency. These results agree with electrophysiology and immunohistochemistry findings, indicating the feasibility of auditory fMRI in rats. This is the first fMRI study of the rodent ascending auditory pathway.

Repeatedly pairing vagus nerve stimulation with a movement reorganizes primary motor cortex

Although sensory and motor systems support different functions, both systems exhibit experience-dependent cortical plasticity under similar conditions. If mechanisms regulating cortical plasticity are common to sensory and motor cortices, then methods generating plasticity in sensory cortex should be effective in motor cortex. Repeatedly pairing a tone with a brief period of vagus nerve stimulation (VNS) increases the proportion of primary auditory cortex responding to the paired tone (Engineer ND, Riley JR, Seale JD, Vrana WA, Shetake J, Sudanagunta SP, Borland MS, Kilgard MP. 2011. Reversing pathological neural activity using targeted plasticity. Nature. 470:101-104). In this study, we predicted that repeatedly pairing VNS with a specific movement would result in an increased representation of that movement in primary motor cortex. To test this hypothesis, we paired VNS with movements of the distal or proximal forelimb in 2 groups of rats. After 5 days of VNS movement pairing, intracranial microstimulation was used to quantify the organization of primary motor cortex. Larger cortical areas were associated with movements paired with VNS. Rats receiving identical motor training without VNS pairing did not exhibit motor cortex map plasticity. These results suggest that pairing VNS with specific events may act as a general method for increasing cortical representations of those events. VNS movement pairing could provide a new approach for treating disorders associated with abnormal movement representations.

Direct imaging of macrovascular and microvascular contributions to BOLD fMRI in layers IV-V of the rat whisker-barrel cortex

The spatiotemporal characteristics of the hemodynamic response to increased neural activity were investigated at the level of individual intracortical vessels using BOLD-fMRI in a well-established rodent model of somatosensory stimulation at 11.7 T. Functional maps of the rat barrel cortex were obtained at 150 × 150 × 500 μm spatial resolution every 200 ms. The high spatial resolution allowed separation of active voxels into those containing intracortical macro vessels, mainly vein/venules (referred to as macrovasculature), and those enriched with arteries/capillaries and small venules (referred to as microvasculature) since the macro vessel can be readily mapped due to the fast T2 decay of blood at 11.7 T. The earliest BOLD response was observed within layers IV-V by 0.8s following stimulation and encompassed mainly the voxels containing the microvasculature and some confined macrovasculature voxels. By 1.2s, the BOLD signal propagated to the macrovasculature voxels where the peak BOLD signal was 2-3 times higher than that of the microvasculature voxels. The BOLD response propagated in individual venules/veins far from neuronal sources at later times. This was also observed in layers IV-V of the barrel cortex after specific stimulation of separated whisker rows. These results directly visualized that the earliest hemodynamic changes to increased neural activity occur mainly in the microvasculature and spread toward the macrovasculature. However, at peak response, the BOLD signal is dominated by penetrating venules even at layers IV-V of the cortex.

MRI of neuronal plasticity in rodent models

Modifications in the behavior and architecture of neuronal networks are well documented to occur in association with learning and memory, as well as following injury. These plasticity mechanisms are crucial to ensure adequate processing of stimuli, and they also dictate the degree of recovery following peripheral or central nervous system injury. Nevertheless, the underlying neuronal mechanisms that determine the degree of plasticity of neuronal pathways are not fully understood. Recent developments in animal-dedicated magnetic resonance imaging (MRI) scanners and related hardware afford a high spatial and temporal resolution, making functional MRI and manganese-enhanced MRI emerging tools for studying reorganization of neuronal pathways in rodent models. Many of the observed changes in neuronal functions in rodent's brains following injury discussed here agree with clinical human fMRI findings. This demonstrates that animal model imaging can have a significant clinical impact in the neuronal plasticity and rehabilitation arenas.

Mapping cortical representations of the rodent forepaw and hindpaw with BOLD fMRI reveals two spatial boundaries

Electrical stimulation of the rat forepaw and hindpaw was employed to study the spatial distribution of BOLD fMRI. Averaging of multiple fMRI sessions significantly improved the spatial stability of the BOLD signal and enabled quantitative determination of the boundaries of the BOLD fMRI maps. The averaged BOLD fMRI signal was distributed unevenly over the extent of the map and the data at the boundaries could be modeled with major and minor spatial components. Comparison of three-dimensional echo-planar imaging (EPI) fMRI at isotropic 300 μm resolution demonstrated that the border locations of the major spatial component of BOLD signal did not overlap between the forepaw and hindpaw maps. Interestingly, the border positions of the minor BOLD fMRI spatial components extended significantly into neighboring representations. Similar results were found for cerebral blood volume (CBV) weighted fMRI obtained using iron oxide particles, suggesting that the minor spatial components may not be due to vascular mislocalization typically associated with BOLD fMRI. Comparison of the BOLD fMRI maps of the forepaw and hindpaw to histological determination of these representations using cytochrome oxidase (CO) staining demonstrated that the major spatial component of the BOLD fMRI activation maps accurately localizes the borders. Finally, 2-3 weeks following peripheral nerve denervation, cortical reorganization/plasticity at the boundaries of somatosensory limb representations in adult rat brain was studied. Denervation of the hindpaw caused a growth in the major component of forepaw representation into the adjacent border of hindpaw representation, such that fitting to two components no longer led to a better fit as compared to using one major component. The border of the representation after plasticity was the same as the border of its minor component in the absence of any plasticity. It is possible that the minor components represent either vascular effects that extend from the real neuronal representations or the neuronal communication between neighboring regions. Either way the results will be useful for studying mechanisms of plasticity that cause alterations in the boundaries of neuronal representations.

Mapping plasticity in the forepaw digit barrel subfield of rat brains using functional MRI

The topographic organization of the forepaw barrel subfield in layer IV of rat primary somatosensory cortex (S1) is a good model for studying neural function and plasticity. The goal of this study was to test the feasibility of functional MRI (fMRI) to map the forepaw digit representations in the S1 of the rat and its plasticity after digit amputation. Three dimensional echo-planar imaging with 300 micron isotropic resolution at 11.7 T was used to achieve high signal-to-noise ratios and laminar layer resolution. By alternating electrical stimulation of the 2nd (D2) and 4th (D4) digits, functional activation in layer IV of the barrel subfields could be distinguished using a differential analysis. Furthermore, 2 and a half months after the amputation of the 3rd digit in baby rats, the overlapping area between D2 and D4 representations was increased. This indicates that the forepaw barrel subfield previously associated with the ablated digit is now associated with the representation of nearby digits, which is consistent with studies using electrophysiology and cytochrome oxidase staining.

Mislocalization of near-threshold tactile stimuli in humans: a central or peripheral phenomenon?

Principles of brain function can be disclosed by studying their limits during performance. Tactile stimuli with near-threshold intensities have been used to assess features of somatosensory processing. When stimulating fingers of one hand using near-threshold intensities, localization errors are observed that deviate significantly from responses obtained by guessing - incorrectly located stimuli are attributed more often to fingers neighbouring the stimulated one than to more distant fingers. Two hypotheses to explain the findings are proposed. The 'central hypothesis' posits that the degree of overlap of cortical tactile representations depends on stimulus intensity, with representations less separated for near-threshold stimuli than for suprathreshold stimuli. The 'peripheral hypothesis' assumes that systematic mislocalizations are due to activation of different sets of skin receptors with specific thresholds. The present experiments were designed to decide between the two hypotheses. Taking advantage of the frequency tuning of somatosensory receptors, their contribution to systematic misclocalizations was studied. In the first experiment, mislocalization profiles were investigated using vibratory stimuli with frequencies of 10, 20 and 100 Hz. Unambiguous mislocalization effects were only obtained for the 10-Hz stimulation, precluding the involvement of Pacinian corpuscles in systematic mislocalization. In the second experiment, Pacinian corpuscles were functionally eliminated by applying a constant 100-Hz vibratory masking stimulus together with near-threshold pulses. Despite masking, systematic mislocation patterns were observed rendering the involvement of Pacinian corpuscles unlikely. The results of both experiments are in favor of the 'central hypothesis' assuming that the extent of overlap in somatosensory representations is modulated by stimulus intensity.

CNS animal fMRI in pain and analgesia

Animal imaging of brain systems offers exciting opportunities to better understand the neurobiology of pain and analgesia. Overall functional studies have lagged behind human studies as a result of technical issues including the use of anesthesia. Now that many of these issues have been overcome including the possibility of imaging awake animals, there are new opportunities to study whole brain systems neurobiology of acute and chronic pain as well as analgesic effects on brain systems de novo (using pharmacological MRI) or testing in animal models of pain. Understanding brain networks in these areas may provide new insights into translational science, and use neural networks as a "language of translation" between preclinical to clinical models. In this review we evaluate the role of functional and anatomical imaging in furthering our understanding in pain and analgesia.