Investigation of the BOLD and CBV fMRI responses to somatosensory stimulation in awake marmosets (Callithrix jacchus).

Understanding the spatiotemporal features of the hemodynamic response function (HRF) to brain stimulation is essential for the correct application of neuroimaging methods to study brain function. Here, we investigated the spatiotemporal evolution of the blood oxygen level-dependent (BOLD) and cerebral blood volume (CBV) HRF in conscious, awake marmosets (Callithrix jacchus), a New World non-human primate with a lissencephalic brain and with growing use in biomedical research. The marmosets were acclimatized to head fixation and placed in a 7-T magnetic resonance imaging (MRI) scanner. Somatosensory stimulation (333-µs pulses; amplitude, 2 mA; 64 Hz) was delivered bilaterally via pairs of contact electrodes. A block design paradigm was used in which the stimulus duration increased in pseudo-random order from a single pulse up to 256 electrical pulses (4 s). For CBV measurements, 30 mg/kg of ultrasmall superparamagnetic ironoxide particles (USPIO) injected intravenously, were used. Robust BOLD and CBV HRFs were obtained in the primary somatosensory cortex (S1), secondary somatosensory cortex (S2) and caudate at all stimulus conditions. In particular, BOLD and CBV responses to a single 333-µs-long stimulus were reliably measured, and the CBV HRF presented shorter onset time and time to peak than the BOLD HRF. Both the size of the regions of activation and the peak amplitude of the HRFs grew quickly with increasing stimulus duration, and saturated for stimulus durations greater than 1 s. Onset times in S1 and S2 were faster than in caudate. Finally, the fine spatiotemporal features of the HRF in awake marmosets were similar to those obtained in humans, indicating that the continued refinement of awake non-human primate models is essential to maximize the applicability of animal functional MRI studies to the investigation of human brain function.

Improved mapping of interictal epileptiform discharges with EEG-fMRI and voxel-wise functional connectivity analysis.

OBJECTIVE: We describe a novel method to spatially map interictal epileptiform discharges (IEDs) through voxel-wise functional connectivity analysis of the functional magnetic resonance imaging (fMRI) portion of simultaneous electroencephalography (EEG)-fMRI data. This method measures the local synchronicity of fMRI signals associated with IED and, in contrast to conventional methods, does not require modeling of neural activities or hemodynamic response. METHODS: Simultaneous EEG-fMRI was performed on six patients with focal epilepsy. IED events were detected from the EEG data. The fMRI data was subdivided into time segments of 20 s in length, and then reorganized into one set of concatenated time series containing the IED events and many sets without IEDs. Local degree centrality (LDC), a metric of functional connectivity, was computed for each brain voxel to summarize its signal correlations to brain voxels within 14 mm of physical distance. This computation was repeated for each set of concatenated time series, yielding one whole-brain LDC map for time with the IED events and many maps for time without IED. A statistical score was computed for each voxel to detect the voxels with significant LDC value differences associated with IEDs. The fMRI data were also processed separately by conventional methods for comparison. RESULTS: In all six patients, regions with significant LDC increase during IEDs were concordant in location to both simultaneous EEG and the epileptogenic focus determined from separate clinical studies. In contrast, results from the conventional methods were concordant in only three patients. SIGNIFICANCE: We show that for focal epilepsy, voxel-wise functional connectivity analysis of EEG-fMRI data may improve IED localization and EEG concordance compared to the conventional analysis. This new analytic method may improve the robustness of interictal EEG-fMRI as a technique for mapping the epileptogenic focus, and helps study the local synchronization aspect of the epileptic network.

fMRI in the awake marmoset: somatosensory-evoked responses, functional connectivity, and comparison with propofol anesthesia.

Functional neuroimaging in animal models is essential for understanding the principles of neurovascular coupling and the physiological basis of fMRI signals that are widely used to study sensory and cognitive processing in the human brain. While hemodynamic responses to sensory stimuli have been characterized in humans, animal studies are able to combine very high resolution imaging with invasive measurements and pharmacological manipulation. To date, most high-resolution studies of neurovascular coupling in small animals have been carried out in anesthetized rodents. Here we report fMRI experiments in conscious, awake common marmosets (Callithrix jacchus), and compare responses to animals anesthetized with propofol. In conscious marmosets, robust BOLD fMRI responses to somatosensory stimulation of the forearm were found in contralateral and ipsilateral regions of the thalamus, primary (SI) and secondary (SII) somatosensory cortex, and the caudate nucleus. These responses were markedly stronger than those in anesthetized marmosets and showed a monotonic increase in the amplitude of the BOLD response with stimulus frequency. On the other hand, anesthesia significantly attenuated responses in thalamus, SI and SII, and abolished responses in caudate and ipsilateral SI. Moreover, anesthesia influenced several other aspects of the fMRI responses, including the shape of the hemodynamic response function and the interareal (SI-SII) spontaneous functional connectivity. Together, these findings demonstrate the value of the conscious, awake marmoset model for studying physiological responses in the somatosensory pathway, in the absence of anesthesia, so that the data can be compared most directly to fMRI in conscious humans.

Rapid high-resolution three-dimensional mapping of T1 and age-dependent variations in the non-human primate brain using magnetization-prepared rapid gradient-echo (MPRAGE) sequence.

The use of quantitative T(1) mapping in neuroscience and neurology has raised strong interest in the development of T(1)-mapping techniques that can measure T(1) in the whole brain, with high accuracy and precision and within short imaging and computation times. Here, we present a new inversion-recovery (IR) based T(1)-mapping method using a standard 3D magnetization-prepared rapid gradient-echo (MPRAGE) sequence. By varying only the inversion time (TI), but keeping other parameters constant, MPRAGE image signals become linear to exp(-TI/T(1)), allowing for accurate T(1) estimation without flip angle correction. We also show that acquiring data at just 3 TIs, with the three different TI values optimized, gives maximum T(1) precision per unit time, allowing for new efficient approaches to measure and compute T(1). We demonstrate the use of our method at 7 T to obtain 3D T(1) maps of the whole brain in common marmosets at 0.60mm resolution and within 11 min. T(1) maps from the same individuals were highly reproducible across different days. Across subjects, the peak of cerebral gray matter T(1) distribution was 1735±52 ms, and the lower edge of cerebral white matter T(1) distribution was 1270±43 ms. We found a significant decrease of T(1) in both gray and white matter of the marmoset brain with age over a span of 14 years, in agreement with previous human studies. This application illustrates that MPRAGE-based 3D T(1) mapping is rapid, accurate and precise, and can facilitate high-resolution anatomical studies in neuroscience and neurological diseases.

Longitudinal functional magnetic resonance imaging in animal models.

Functional magnetic resonance imaging (fMRI) has had an essential role in furthering our understanding of brain physiology and function. fMRI techniques are nowadays widely applied in neuroscience research, as well as in translational and clinical studies. The use of animal models in fMRI studies has been fundamental in helping elucidate the mechanisms of cerebral blood-flow regulation, and in the exploration of basic neuroscience questions, such as the mechanisms of perception, behavior, and cognition. Because animals are inherently non-compliant, most fMRI performed to date have required the use of anesthesia, which interferes with brain function and compromises interpretability and applicability of results to our understanding of human brain function. An alternative approach that eliminates the need for anesthesia involves training the animal to tolerate physical restraint during the data acquisition. In the present chapter, we review these two different approaches to obtaining fMRI data from animal models, with a specific focus on the acquisition of longitudinal data from the same subjects.

Visualizing the entire cortical myelination pattern in marmosets with magnetic resonance imaging.

Myeloarchitecture, the pattern of myelin density across the cerebral cortex, has long been visualized in histological sections to identify distinct anatomical areas of the cortex. In humans, two-dimensional (2D) magnetic resonance imaging (MRI) has been used to visualize myeloarchitecture in select areas of the cortex, such as the stripe of Gennari in the primary visual cortex and Heschl's gyrus in the primary auditory cortex. Here, we investigated the use of MRI contrast based on longitudinal relaxation time (T(1)) to visualize myeloarchitecture in vivo over the entire cortex of the common marmoset (Callithrix jacchus), a small non-human primate that is becoming increasingly important in neuroscience and neurobiology research. Using quantitative T(1) mapping, we found that T(1) at 7T in a cortical region with a high myelin content was 15% shorter than T(1) in a region with a low myelin content. To maximize this T(1) contrast for imaging cortical myelination patterns, we optimized a magnetization-prepared rapidly acquired gradient echo (MP-RAGE) sequence. In whole-brain, 3D T(1)-weighted images made in vivo with the sequence, we identified six major cortical areas with high myelination and confirmed the results with histological sections stained for myelin. We also identified several subtle features of myeloarchitecture, showing the sensitivity of our technique. The ability to image myeloarchitecture over the entire cortex may prove useful in studies of longitudinal changes of the topography of the cortex associated with development and neuronal plasticity, as well as for guiding and confirming the location of functional measurements.

Assessment of stimulus-induced changes in human V1 visual field maps.

Visual cortex contains a set of field maps in which nearby scene points are represented in the responses of nearby neurons. We tested a recent hypothesis that the visual field map in primary visual cortex (V1) is dynamic, changing in response to stimulus motion direction. The original experimental report replicates, but further experimental and analytical investigations do not support, the interpretation of the results. The V1 map remains invariant when measured using stimuli moving in different directions. The measurements can be explained by small and systematic response amplitude differences that arise when probing with stimuli moving in different directions.