A predictor-corrector phase unwrapping algorithm for temporally undersampled gradient-echo MRI

PURPOSE

To develop a method for unwrapping temporally undersampled and nonlinear gradient recalled echo (GRE) phase.

THEORY AND METHODS

Temporal unwrapping is performed as a sequential one step prediction of the echo phase, followed by a correction to the nearest integer wrap-count. A spatio-temporal extension of the 1D predictor corrector unwrapping (PCU) algorithm improves the prediction accuracy, and thereby maintains spatial continuity. The proposed method is evaluated using numerical phantom, physical phantom, and in vivo brain data at both 3 T and 9.4 T. The unwrapping performance is compared with the state-of-the-art temporal and spatial unwrapping algorithms, and the spatio-temporal iterative virtual-echo based Nyquist sampled (iVENyS) algorithm.

RESULTS

Simulation results showed significant reduction in unwrapping errors at higher echoes compared with the state-of-the-art algorithms. Similar to the iVENyS algorithm, the PCU algorithm was able to generate spatially smooth phase images for in vivo data acquired at 3 T and 9.4 T, bypassing the use of additional spatial unwrapping step. A key advantage over iVENyS algorithm is the superior performance of PCU algorithm at higher echoes.

CONCLUSION

PCU algorithm serves as a robust phase unwrapping method for temporally undersampled and nonlinear GRE phase, particularly in the presence of high field gradients.

MP2RAGE vs. MPRAGE surface-based morphometry in focal epilepsy

In drug-resistant focal epilepsy, detecting epileptogenic lesions using MRI poses a critical diagnostic challenge. Here, we assessed the utility of MP2RAGE-a T1-weighted sequence with self-bias correcting properties commonly utilized in ultra-high field MRI-for the detection of epileptogenic lesions using a surface-based morphometry pipeline based on FreeSurfer, and compared it to the common approach using T1w MPRAGE, both at 3T. We included data from 32 patients with focal epilepsy (5 MRI-positive, 27 MRI-negative with lobar seizure onset hypotheses) and 94 healthy controls from two epilepsy centres. Surface-based morphological measures and intensities were extracted and evaluated in univariate GLM analyses as well as multivariate unsupervised 'novelty detection' machine learning procedures. The resulting prediction maps were analyzed over a range of possible thresholds using alternative free-response receiver operating characteristic (AFROC) methodology with respect to the concordance with predefined lesion labels or hypotheses on epileptogenic zone location. We found that MP2RAGE performs at least comparable to MPRAGE and that especially analysis of MP2RAGE image intensities may provide additional diagnostic information. Secondly, we demonstrate that unsupervised novelty-detection machine learning approaches may be useful for the detection of epileptogenic lesions (maximum AFROC AUC 0.58) when there is only a limited lesional training set available. Third, we propose a statistical method of assessing lesion localization performance in MRI-negative patients with lobar hypotheses of the epileptogenic zone based on simulation of a random guessing process as null hypothesis. Based on our findings, it appears worthwhile to study similar surface-based morphometry approaches in ultra-high field MRI (≥ 7 T).

Distribution of corpora amylacea in the human midbrain: using synchrotron radiation phase-contrast microtomography, high-field magnetic resonance imaging, and histology

Corpora amylacea (CA) are polyglucosan aggregated granules that accumulate in the human body throughout aging. In the cerebrum, CA have been found in proximity to ventricular walls, pial surfaces, and blood vessels. However, studies showing their three-dimensional spatial distribution are sparse. In this study, volumetric images of four human brain stems were obtained with MRI and phase-contrast X-ray microtomography, followed up by Periodic acid Schiff stain for validation. CA appeared as hyperintense spheroid structures with diameters up to 30 μm. An automatic pipeline was developed to segment the CA, and the spatial distribution of over 200,000 individual corpora amylacea could be investigated. A threefold-or higher-density of CA was detected in the dorsomedial column of the periaqueductal gray (860-4,200 CA count/mm3) than in the superior colliculus (150-340 CA count/mm3). We estimated that about 2% of the CA were located in the immediate vicinity of the vessels or in the peri-vascular space. While CA in the ependymal lining of the cerebral aqueduct was rare, the sub-pial tissue of the anterior and posterior midbrain contained several CA. In the sample with the highest CA density, quantitative maps obtained with MRI revealed high R2* values and a diamagnetic shift in a region which spatially coincided with the CA dense region.

Prediction of motion induced magnetic fields for human brain MRI at 3 T

OBJECTIVE

Maps of B0 field inhomogeneities are often used to improve MRI image quality, even in a retrospective fashion. These field inhomogeneities depend on the exact head position within the static field but acquiring field maps (FM) at every position is time consuming. Here we propose a forward simulation strategy to obtain B0 predictions at different head-positions.

METHODS

FM were predicted by combining (1) a multi-class tissue model for estimation of tissue-induced fields, (2) a linear k-space model for capturing gradient imperfections, (3) a dipole estimation for quantifying lower-body perturbing fields (4) and a position-dependent tissue mask to model FM alterations caused by large motion effects. The performance of the combined simulation strategy was compared with an approach based on a rigid body transformation of the FM measured in the reference position to the new position.

RESULTS

The transformed FM provided inconsistent results for large head movements (> 5° rotation, approximately), while the simulation strategy had a superior prediction accuracy for all positions. The simulated FM was used to optimize B0 shims with up to 22.2% improvement with respect to the transformed FM approach.

CONCLUSION

The proposed simulation strategy is able to predict movement-induced B0 field inhomogeneities yielding more precise estimates of the ground truth field homogeneity than the transformed FM.

Simulated and experimental phantom data for multi-center quality assurance of quantitative susceptibility maps at 3 T, 7 T and 9.4 T

PURPOSE

To develop methods for quality assurance of quantitative susceptibility mapping (QSM) using MRI at different magnetic field strengths, and scanners, using different MR-sequence protocols, and post-processing pipelines.

METHODS

We built a custom phantom based on iron in two forms: homogeneous susceptibility ('free iron') and with fine-scaled variations in susceptibility ('clustered iron') at different iron concentrations. The phantom was measured at 3.0 T (two scanners), 7.0 T and 9.4 T using multi-echo, gradient echo acquisition sequences. A digital phantom analogue to the iron-phantom, tailored to obtain similar results as in experimentation was developed, with similar geometry and susceptibility values. Morphology enabled dipole inversion was applied to the phase images to obtain QSM for experimental and simulated data using the MEDI + 0 approach for background regularization.

RESULTS

Across all scanners, QSM-values showed a linear increase with iron concentrations. The QSM-relaxivity was 0.231 ± 0.047 ppm/mM for free and 0.054 ± 0.013 ppm/mM for clustered iron, with adjusted determination coefficients (DoC) ≥ 0.87. Similarly, the simulations yielded linear increases (DoC ≥ 0.99). In both the experimental and digital phantoms, the estimated molar susceptibility was lower with clustered iron, because clustering led to highly localized field effects.

CONCLUSION

Our iron phantom can be used to evaluate the capability of QSM to detect local variations in susceptibility across different field strengths, when using different MR-sequence protocols. The devised simulation method captures the effect of iron clustering in QSM as seen experimentally and could be used in the future to optimize QSM processing pipelines and achieve higher accuracy for local field effects, as also seen in Alzheimer's beta-amyloid plaques.

Phase-based masking for quantitative susceptibility mapping of the human brain at 9.4T

PURPOSE

To develop improved tissue masks for QSM.

METHODS

Masks including voxels at the brain surface were automatically generated from the magnitude alone (MM) or combined with test functions from the first (PG) or second (PB) derivative of the sign of the wrapped phase. Phase images at 3T and 9.4T were simulated at different TEs and used to generate a mask, PItoh , with between-voxel phase differences less than π. MM, PG, and PB were compared with PItoh . QSM were generated from 3D multi-echo gradient-echo data acquired at 9.4T (21 subjects aged: 20-56y), and from the QSM2016 challenge 3T data using different masks, unwrapping, background removal, and dipole inversion algorithms. QSM contrast was quantified using age-based iron concentrations.

RESULTS

Close to air cavities, phase wraps became denser with increasing field and echo time, yielding increased values of the test functions. Compared with PItoh , PB had the highest Dice coefficient, while PG had the lowest and MM the highest percentage of voxels outside PItoh. Artifacts observed in QSM at 9.4T with MM were mitigated by stronger background filters but yielded a reduced QSM contrast. With PB, QSM contrast was greater and artifacts diminished. Similar results were obtained with challenge data, evidencing larger effects of mask close to air cavities.

CONCLUSION

Automatic, phase-based masking founded on the second derivative of the sign of the wrapped phase, including cortical voxels at the brain surface, was able to mitigate artifacts and restore QSM contrast across cortical and subcortical brain regions.

Characterization of the blood oxygen level dependent hemodynamic response function in human subcortical regions with high spatiotemporal resolution

Subcortical brain regions are absolutely essential for normal human function. These phylogenetically early brain regions play critical roles in human behaviors such as the orientation of attention, arousal, and the modulation of sensory signals to cerebral cortex. Despite the critical health importance of subcortical brain regions, there has been a dearth of research on their neurovascular responses. Blood oxygen level dependent (BOLD) functional MRI (fMRI) experiments can help fill this gap in our understanding. The BOLD hemodynamic response function (HRF) evoked by brief (<4 s) neural activation is crucial for the interpretation of fMRI results because linear analysis between neural activity and the BOLD response relies on the HRF. Moreover, the HRF is a consequence of underlying local blood flow and oxygen metabolism, so characterization of the HRF enables understanding of neurovascular and neurometabolic coupling. We measured the subcortical HRF at 9.4T and 3T with high spatiotemporal resolution using protocols that enabled reliable delineation of HRFs in individual subjects. These results were compared with the HRF in visual cortex. The HRF was faster in subcortical regions than cortical regions at both field strengths. There was no significant undershoot in subcortical areas while there was a significant post-stimulus undershoot that was tightly coupled with its peak amplitude in cortex. The different BOLD temporal dynamics indicate different vascular dynamics and neurometabolic responses between cortex and subcortical nuclei.

Microvascular imaging of the unstained human superior colliculus using synchrotron-radiation phase-contrast microtomography

Characterizing the microvasculature of the human brain is critical to advance understanding of brain vascular function. Most methods rely on tissue staining and microscopy in two-dimensions, which pose several challenges to visualize the three-dimensional structure of microvessels. In this study, we used an edge-based segmentation method to extract the 3D vasculature from synchrotron radiation phase-contrast microtomography (PC-μCT) of two unstained, paraffin-embedded midbrain region of the human brain stem. Vascular structures identified in PC-μCT were validated with histology of the same specimen. Using the Deriche-Canny edge detector that was sensitive to the boundary between tissue and vascular space, we could segment the vessels independent of signal variations in PC-μCT images. From the segmented volumetric vasculature, we calculated vessel diameter, vessel length and volume fraction of the vasculature in the superior colliculi. From high resolution images, we found the most frequent vessel diameter to be between 8.6-10.2 µm. Our findings are consistent with the known anatomy showing two types of vessels with distinctive morphology: peripheral collicular vessels and central collicular vessels. The proposed method opens up new possibilities for vascular research of the central nervous system using synchrotron radiation PC-μCT of unstained human tissue.

Perception is associated with the brain's metabolic response to sensory stimulation

Processing of incoming sensory stimulation triggers an increase of cerebral perfusion and blood oxygenation (neurovascular response) as well as an alteration of the metabolic neurochemical profile (neurometabolic response). Here we show in human primary visual cortex (V1) that perceived and unperceived isoluminant chromatic flickering stimuli designed to have similar neurovascular responses as measured by blood oxygenation level dependent functional MRI (BOLD-fMRI) have markedly different neurometabolic responses as measured by functional MRS. In particular, a significant regional buildup of lactate, an index of aerobic glycolysis, and glutamate, an index of malate-aspartate shuttle, occurred in V1 only when the flickering was perceived, without any relation with behavioral or physiological variables. Whereas the BOLD-fMRI signal in V1, a proxy for input to V1, was insensitive to flickering perception by design, the BOLD-fMRI signal in secondary visual areas was larger during perceived than unperceived flickering, indicating increased output from V1. These results demonstrate that the upregulation of energy metabolism induced by visual stimulation depends on the type of information processing taking place in V1, and that 1H-fMRS provides unique information about local input/output balance that is not measured by BOLD fMRI.

Disturbed Balance of Inhibitory Signaling Links Hearing Loss and Cognition

Neuronal hyperexcitability in the central auditory pathway linked to reduced inhibitory activity is associated with numerous forms of hearing loss, including noise damage, age-dependent hearing loss, and deafness, as well as tinnitus or auditory processing deficits in autism spectrum disorder (ASD). In most cases, the reduced central inhibitory activity and the accompanying hyperexcitability are interpreted as an active compensatory response to the absence of synaptic activity, linked to increased central neural gain control (increased output activity relative to reduced input). We here suggest that hyperexcitability also could be related to an immaturity or impairment of tonic inhibitory strength that typically develops in an activity-dependent process in the ascending auditory pathway with auditory experience. In these cases, high-SR auditory nerve fibers, which are critical for the shortest latencies and lowest sound thresholds, may have either not matured (possibly in congenital deafness or autism) or are dysfunctional (possibly after sudden, stressful auditory trauma or age-dependent hearing loss linked with cognitive decline). Fast auditory processing deficits can occur despite maintained basal hearing. In that case, tonic inhibitory strength is reduced in ascending auditory nuclei, and fast inhibitory parvalbumin positive interneuron (PV-IN) dendrites are diminished in auditory and frontal brain regions. This leads to deficits in central neural gain control linked to hippocampal LTP/LTD deficiencies, cognitive deficits, and unbalanced extra-hypothalamic stress control. Under these conditions, a diminished inhibitory strength may weaken local neuronal coupling to homeostatic vascular responses required for the metabolic support of auditory adjustment processes. We emphasize the need to distinguish these two states of excitatory/inhibitory imbalance in hearing disorders: (i) Under conditions of preserved fast auditory processing and sustained tonic inhibitory strength, an excitatory/inhibitory imbalance following auditory deprivation can maintain precise hearing through a memory linked, transient disinhibition that leads to enhanced spiking fidelity (central neural gain⇑) (ii) Under conditions of critically diminished fast auditory processing and reduced tonic inhibitory strength, hyperexcitability can be part of an increased synchronization over a broader frequency range, linked to reduced spiking reliability (central neural gain⇓). This latter stage mutually reinforces diminished metabolic support for auditory adjustment processes, increasing the risks for canonical dementia syndromes.

Developing formalin-based fixative agents for post mortem brain MRI at 9.4 T

PURPOSE

To develop fixative agents for high-field MRI with suitable dielectric properties and measure MR properties in immersion-fixed brain tissue.

METHODS

Dielectric properties of formalin-based agents were assessed (100 MHz-4.5 GHz), and four candidate fixatives with/without polyvinylpyrrolidone (PVP) and different salt concentrations were formulated. B₁ field and MR properties (T₁ , R 2 ∗ , R₂ , R 2 ' , and magnetic susceptibility [QSM]) were observed in white and gray matter of pig brain samples during 0.5-35 days of immersion fixation. The kinetics were fitted using exponential functions. The immersion time required to reach maximum R 2 ∗ values at different tissue depths was used to estimate the Medawar coefficient for fixative penetration. The effect of replacing the fixatives with Fluoroinert and phosphate-buffered saline as embedding media was also evaluated.

RESULTS

The dielectric properties of formalin were nonlinearly modified by increasing amounts of additives. With 5% PVP and 0.04% NaCl, the dielectric properties and B₁ field reflected in vivo conditions. The highest B₁ values were found in white matter with PVP and varied significantly with tissue depth and embedding media, but not with immersion time. The MR properties depended on PVP yielding lower T₁ , higher R 2 ∗ , more paramagnetic QSM values, and a lower Medawar coefficient (0.9 mm / h ; without PVP: 1.5). Regardless of fixative, switching to phosphate-buffered saline as embedder caused a paramagnetic shift in QSM and decreased R 2 ∗ that progressed during 1 month of storage, whereas no differences were found with Fluorinert.

CONCLUSION

In vivo-like B₁ fields can be achieved in formalin fixatives using PVP and a low salt concentration, yielding lower T₁ , higher R 2 ∗ , and more paramagnetic QSM than without additives. The kinetics of R 2 ∗ allowed estimation of fixative tissue penetration.

T2-Pseudonormalization and Microstructural Characterization in Advanced Stages of Late-infantile Metachromatic Leukodystrophy

PURPOSE

T2-weighted signal hyperintensities in white matter (WM) are a diagnostic finding in brain magnetic resonance imaging (MRI) of patients with metachromatic leukodystrophy (MLD). In our systematic investigation of the evolution of T2-hyperintensities in patients with the late-infantile form, we describe and characterize T2-pseudonormalization in the advanced stage of the natural disease course.

METHODS

The volume of T2-hyperintensities was quantified in 34 MRIs of 27 children with late-infantile MLD (median age 2.25 years, range 0.5-5.2 years). In three children with the most advanced clinical course (age >4 years) and for whom the T2-pseudonormalization was the most pronounced, WM microstructure was investigated using a multimodal MRI protocol, including diffusion-weighted imaging, MR spectroscopy (MRS), myelin water fraction (MWF), magnetization transfer ratio (MTR), T1-mapping and quantitative susceptibility mapping.

RESULTS

T2-hyperintensities in cerebral WM returned to normal in large areas of 3 patients in the advanced disease stage. Multimodal assessment of WM microstructure in areas with T2-pseudonormalization revealed highly decreased values for NAA, neurite density, isotropic water, mean and radial kurtosis, MWF and MTR, as well as increased radial diffusivity.

CONCLUSION

In late-infantile MLD patients, we found T2-pseudonormalization in WM tissue with highly abnormal microstructure characterizing the most advanced disease stage. Pathological hallmarks might be a loss of myelin, but also neuronal loss as well as increased tissue density due to gliosis and accumulated storage material. These results suggest that a multimodal MRI protocol using more specific microstructural parameters than T2-weighted sequences should be used when evaluating the effect of treatment trials in MLD.

Quantitative Susceptibility Mapping of the Basal Ganglia and Thalamus at 9.4 Tesla

The thalamus (Th) and basal ganglia (BG) are central subcortical connectivity hubs of the human brain, whose functional anatomy is still under intense investigation. Nevertheless, both substructures contain a robust and reproducible functional anatomy. The quantitative susceptibility mapping (QSM) at ultra-high field may facilitate an improved characterization of the underlying functional anatomy in vivo. We acquired high-resolution QSM data at 9.4 Tesla in 21 subjects, and analyzed the thalamic and BG by using a prior defined functional parcellation. We found a more substantial contribution of paramagnetic susceptibility sources such as iron in the pallidum in contrast to the caudate, putamen, and Th in descending order. The diamagnetic susceptibility sources such as myelin and calcium revealed significant contributions in the Th parcels compared with the BG. This study presents a detailed nuclei-specific delineation of QSM-provided diamagnetic and paramagnetic susceptibility sources pronounced in the BG and the Th. We also found a reasonable interindividual variability as well as slight hemispheric differences. The results presented here contribute to the microstructural knowledge of the Th and the BG. In specific, the study illustrates QSM values (myelin, calcium, and iron) in functionally similar subregions of the Th and the BG.

Multi-echo gradient-recalled-echo phase unwrapping using a Nyquist sampled virtual echo train in the presence of high-field gradients

PURPOSE

To develop a spatio-temporal approach to accurately unwrap multi-echo gradient-recalled echo phase in the presence of high-field gradients.

THEORY AND METHODS

Using the virtual echo-based Nyquist sampled (VENyS) algorithm, the temporal unwrapping procedure is modified by introduction of one or more virtual echoes between the first lower and the immediate higher echo, so as to reinstate the Nyquist condition at locations with high-field gradients. An iterative extension of the VENyS algorithm maintains spatial continuity by adjusting the phase rotations to make the neighborhood phase differences less than π. The algorithm is evaluated using simulated data, Gadolinium contrast-doped phantom, and in vivo brain, abdomen, and chest data sets acquired at 3 T and 9.4 T. The unwrapping performance is compared with the standard temporal unwrapping algorithm used in the morphology-enabled dipole inversion-QSM pipeline as a benchmark for validation.

RESULTS

Quantitative evaluation using numerical phantom showed significant reduction in unwrapping errors in regions of large field gradients, and the unwrapped phase revealed an exact match with the linear concentration profile of vials in a gadolinium contrast-doped phantom data acquired at 9.4 T. Without the need for additional spatial unwrapping, the iterative VENyS algorithm was able to generate spatially continuous phase images. Application to in vivo data resulted in better unwrapping performance, especially in regions with large susceptibility changes such as the air/tissue interface.

CONCLUSION

The iterative VENyS algorithm serves as a robust unwrapping method for multi-echo gradient-recalled echo phase in the presence of high-field gradients.

Ultra-High Field MRI in Alzheimer's Disease: Effective Transverse Relaxation Rate and Quantitative Susceptibility Mapping of Human Brain In Vivo and Ex Vivo compared to Histology

Alzheimer's disease (AD) is the most common cause of dementia worldwide. So far, diagnosis of AD is only unequivocally defined through postmortem histology. Amyloid plaques are a classical hallmark of AD and amyloid load is currently quantified by Positron Emission tomography (PET) in vivo. Ultra-high field magnetic resonance imaging (UHF-MRI) can potentially provide a non-invasive biomarker for AD by allowing imaging of pathological processes at a very-high spatial resolution. The first aim of this work was to reproduce the characteristic cortical pattern previously observed in vivo in AD patients using weighted-imaging at 7T. We extended these findings using quantitative susceptibility mapping (QSM) and quantification of the effective transverse relaxation rate (R2*) at 9.4T. The second aim was to investigate the origin of the contrast patterns observed in vivo in the cortex of AD patients at 9.4T by comparing quantitative UHF-MRI (9.4T and 14.1T) of postmortem samples with histology. We observed a distinctive cortical pattern in vivo in patients compared to healthy controls (HC), and these findings were confirmed ex vivo. Specifically, we found a close link between the signal changes detected by QSM in the AD sample at 14.1T and the distribution pattern of amyloid plaques in the histological sections of the same specimen. Our findings showed that QSM and R2* maps can distinguish AD from HC at UHF by detecting cortical alterations directly related to amyloid plaques in AD patients. Furthermore, we provided a method to quantify amyloid plaque load in AD patients at UHF non-invasively.

Depth relationships and measures of tissue thickness in dorsal midbrain

Dorsal human midbrain contains two nuclei with clear laminar organization, the superior and inferior colliculi. These nuclei extend in depth between the superficial dorsal surface of midbrain and a deep midbrain nucleus, the periaqueductal gray matter (PAG). The PAG, in turn, surrounds the cerebral aqueduct (CA). This study examined the use of two depth metrics to characterize depth and thickness relationships within dorsal midbrain using the superficial surface of midbrain and CA as references. The first utilized nearest-neighbor Euclidean distance from one reference surface, while the second used a level-set approach that combines signed distances from both reference surfaces. Both depth methods provided similar functional depth profiles generated by saccadic eye movements in a functional MRI task, confirming their efficacy for delineating depth for superficial functional activity. Next, the boundaries of the PAG were estimated using Euclidean distance together with elliptical fitting, indicating that the PAG can be readily characterized by a smooth surface surrounding PAG. Finally, we used the level-set approach to measure tissue depth between the superficial surface and the PAG, thus characterizing the variable thickness of the colliculi. Overall, this study demonstrates depth-mapping schemes for human midbrain that enables accurate segmentation of the PAG and consistent depth and thickness estimates of the superior and inferior colliculi.

"Wrong Way Up": Temporal and Spatial Dynamics of the Networks for Body Motion Processing at 9.4 T

Body motion delivers a wealth of socially relevant information. Yet display inversion severely impedes biological motion (BM) processing. It is largely unknown how the brain circuits for BM are affected by display inversion. As upright and upside-down point-light BM displays are similar, we addressed this issue by using ultrahigh field functional MRI at 9.4 T providing for high sensitivity and spatial resolution. Whole-brain analysis along with exploration of the temporal dynamics of the blood-oxygen-level-dependent response reveals that in the left hemisphere, inverted BM activates anterior networks likely engaged in decision making and cognitive control, whereas readily recognizable upright BM activates posterior areas solely. In the right hemisphere, multiple networks are activated in response to upright BM as compared with scarce activation to inversion. With identical visual input with display inversion, a large-scale network in the right hemisphere is detected in perceivers who do not constantly interpret displays as shown the "wrong way up." For the first time, we uncover (1) (multi)functional involvement of each region in the networks underpinning BM processing and (2) large-scale ensembles of regions playing in unison with distinct temporal dynamics. The outcome sheds light on the neural circuits underlying BM processing as an essential part of the social brain.

In-vivo quantitative structural imaging of the human midbrain and the superior colliculus at 9.4T

We explored anatomical details of the superior colliculus (SC) by in vivo magnetic resonance imaging (MRI) at 9.4T. The high signal-to-noise ratio allowed the acquisition of high resolution, multi-modal images with voxel sizes ranging between 176 × 132 × 600 μm and (800)³μm. Quantitative mapping of the longitudinal relaxation rate R1, the effective transverse relaxation rate R2*, and the magnetic susceptibility QSM was performed in 14 healthy volunteers. The images were analyzed in native space as well as after normalization to a common brain space (MNI). The coefficient-of-variation (CoV) across subjects was evaluated in prominent regions of the midbrain, reaching the best reproducibility (CoV of 5%) in the R2* maps of the SC in MNI space, while the CoV in the QSM maps remained high regardless of brain-space. To investigate whether more complex neurobiological architectural features could be detected, depth profiles through the SC layers towards the red nucleus (RN) were evaluated at different levels of the SC along the rostro-caudal axis. This analysis revealed alterations of the quantitative MRI parameters concordant with previous post mortem histology studies of the cyto- and myeloarchitecture of the SC. In general, the R1 maps were hyperintense in areas characterized by the presence of abundant myelinated fibers, and likely enabled detection of the deep white layer VII of the SC adjacent to the periaqueductal gray. While R1 maps failed to reveal finer details, possibly due to the relatively coarse spatial sampling used for this modality, these could be recovered in R2* maps and in QSM. In the central part of the SC along its rostro-caudal axis, increased R2* values and decreased susceptibility values were observed 2 mm below the SC surface, likely reflecting the myelinated fibers in the superficial optic layer (layer III). Towards the deeper layers, a second increase in R2* was paralleled by a paramagnetic shift in QSM suggesting the presence of an iron-rich layer about 3 mm below the surface of the SC, attributed to the intermediate gray layer (IV) composed of multipolar neurons. These results dovetail observations in histological specimens and animal studies and demonstrate that high-resolution multi-modal MRI at 9.4T can reveal several microstructural features of the SC in vivo.

Assessing White Matter Microstructure in Brain Regions with Different Myelin Architecture Using MRI

OBJECTIVE

We investigate how known differences in myelin architecture between regions along the cortico-spinal tract and frontal white matter (WM) in 19 healthy adolescents are reflected in several quantitative MRI parameters that have been proposed to non-invasively probe WM microstructure. In a clinically feasible scan time, both conventional imaging sequences as well as microstructural MRI parameters were assessed in order to quantitatively characterise WM regions that are known to differ in the thickness of their myelin sheaths, and in the presence of crossing or parallel fibre organisation.

RESULTS

We found that diffusion imaging, MR spectroscopy (MRS), myelin water fraction (MWF), Magnetization Transfer Imaging, and Quantitative Susceptibility Mapping were myelin-sensitive in different ways, giving complementary information for characterising WM microstructure with different underlying fibre architecture. From the diffusion parameters, neurite density (NODDI) was found to be more sensitive than fractional anisotropy (FA), underlining the limitation of FA in WM crossing fibre regions. In terms of sensitivity to different myelin content, we found that MWF, the mean diffusivity and chemical-shift imaging based MRS yielded the best discrimination between areas.

CONCLUSION

Multimodal assessment of WM microstructure was possible within clinically feasible scan times using a broad combination of quantitative microstructural MRI sequences. By assessing new microstructural WM parameters we were able to provide normative data and discuss their interpretation in regions with different myelin architecture, as well as their possible application as biomarker for WM disorders.

Depth-dependence of visual signals in the human superior colliculus at 9.4 T

The superior colliculus (SC) is a layered structure located in the midbrain. We exploited the improved spatial resolution and BOLD signal strength available at 9.4 T to investigate the depth profile of visual BOLD responses in the human SC based on distortion-corrected EPI data with a 1 mm isotropic resolution. We used high resolution (350 µm in-plane) anatomical images to determine regions-of-interest of the SC and applied a semi-automated method to segment it into superficial, intermediate, and deep zones. A greater than linear increase in sensitivity of the functional signal at 9.4 T allowed us to detect a statistically significant depth pattern in a group analysis with a 20 min stimulation paradigm. Descriptive data showed consistent depth profiles also in single individuals. The highest signals were localized to the superficial layers of the right and left SC during contralateral stimulation, which was in good agreement with its functional architecture known from non-human primates. This study thus demonstrates the potential of 9.4 T MRI for functional neuroimaging even in deeply located, particularly challenging brain structures such as the SC. Hum Brain Mapp 38:574-587, 2017. © 2016 Wiley Periodicals, Inc.

MR spectroscopy for in vivo assessment of the oncometabolite 2-hydroxyglutarate and its effects on cellular metabolism in human brain gliomas at 9.4T

PURPOSE

To examine in vivo metabolic alterations in the isocitrate dehydrogenase (IDH) mutated gliomas using magnetic resonance spectroscopy (MRS) at magnetic field 9.4T.

MATERIALS AND METHODS

Spectra were acquired with a 9.4T whole-body scanner with the use of a custom-built head coil (16 channel transmit and 31 channel receive). A modified stimulated echo acquisition mode (STEAM) sequence was used for localization. Eighteen patients with brain tumors of probable glial origin participated in this study. The study was performed in accordance with the guidelines of the local Ethics Committee.

RESULTS

The increased spectral resolution allowed us to directly address metabolic alterations caused by the specific pathophysiology of IDH mutations including the presence of the oncometabolite 2-hydroxglutarate (2HG) and a significant decrease of the pooled glutamate and glutamine (20%, P = 0.024), which probably reflects an attempt to replenish α-ketoglutarate lost by conversion to 2HG. We also observed significantly reduced glutathione (GSH) levels (39%, P = 0.019), which could be similarly caused by depletion of dihydronicotinamide-adenine dinucleotide phosphate (NADPH) during this conversion in IDH mutant gliomas.

CONCLUSION

We demonstrate that MRS at 9.4T provides a noninvasive measure of 2HG in vivo, which may be used for therapy planning and prognostication, and may provide insights into related pathophysiologic metabolic alterations associated with IDH mutations. J. MAGN. RESON. IMAGING 2016;44:823-833.

Diffusion properties of conventional and calcium-sensitive MRI contrast agents in the rat cerebral cortex

Calcium-sensitive MRI contrast agents can only yield quantitative results if the agent concentration in the tissue is known. The agent concentration could be determined by diffusion modeling, if relevant parameters were available. We have established an MRI-based method capable of determining diffusion properties of conventional and calcium-sensitive agents. Simulations and experiments demonstrate that the method is applicable both for conventional contrast agents with a fixed relaxivity value and for calcium-sensitive contrast agents. The full pharmacokinetic time-course of gadolinium concentration estimates was observed by MRI before, during and after intracerebral administration of the agent, and the effective diffusion coefficient D* was determined by voxel-wise fitting of the solution to the diffusion equation. The method yielded whole brain coverage with a high spatial and temporal sampling. The use of two types of MRI sequences for sampling of the diffusion time courses was investigated: Look-Locker-based quantitative T(1) mapping, and T(1) -weighted MRI. The observation times of the proposed MRI method is long (up to 20 h) and consequently the diffusion distances covered are also long (2-4 mm). Despite this difference, the D* values in vivo were in agreement with previous findings using optical measurement techniques, based on observation times of a few minutes. The effective diffusion coefficient determined for the calcium-sensitive contrast agents may be used to determine local tissue concentrations and to design infusion protocols that maintain the agent concentration at a steady state, thereby enabling quantitative sensing of the local calcium concentration.

Investigation of a calcium-responsive contrast agent in cellular model systems: feasibility for use as a smart molecular probe in functional MRI

Responsive or smart contrast agents (SCAs) represent a promising direction for development of novel functional MRI (fMRI) methods for the eventual noninvasive assessment of brain function. In particular, SCAs that respond to Ca(2+) may allow tracking neuronal activity independent of brain vasculature, thus avoiding the characteristic limitations of current fMRI techniques. Here we report an in vitro proof-of-principle study with a Ca(2+)-sensitive, Gd(3+)-based SCA in an attempt to validate its potential use as a functional in vivo marker. First, we quantified its relaxometric response in a complex 3D cell culture model. Subsequently, we examined potential changes in the functionality of primary glial cells following administration of this SCA. Monitoring intracellular Ca(2+) showed that, despite a reduction in the Ca(2+) level, transport of Ca(2+) through the plasma membrane remained unaffected, while stimulation with ATP induced Ca(2+)-transients suggested normal cellular signaling in the presence of low millimolar SCA concentrations. SCAs merely lowered the intracellular Ca(2+) level. Finally, we estimated the longitudinal relaxation times (T1) for an idealized in vivo fMRI experiment with SCA, for extracellular Ca(2+) concentration level changes expected during intense neuronal activity which takes place upon repetitive stimulation. The values we obtained indicate changes in T1 of around 1-6%, sufficient to be robustly detectable using modern MRI methods in high field scanners. Our results encourage further attempts to develop even more potent SCAs and appropriate fMRI protocols. This would result in novel methods that allow monitoring of essential physiological processes at the cellular and molecular level.

Dual-frequency calcium-responsive MRI agents

Responsive or smart magnetic resonance imaging (MRI) contrast agents are molecular sensors that alter the MRI signal upon changes in a particular parameter in their microenvironment. Consequently, they could be exploited for visualization of various biochemical events that take place at molecular and cellular levels. In this study, a set of dual-frequency calcium-responsive MRI agents are reported. These are paramagnetic, fluorine-containing complexes that produce remarkably high MRI signal changes at the (1)H and (19)F frequencies at varying Ca(2+) concentrations. The nature of the processes triggered by Ca(2+) was revealed, allowing a better understanding of these complex systems and their further improvement. The findings indicate that these double-frequency tracers hold great promise for development of novel functional MRI methods.

A smart (19) F and (1) H MRI probe with self-immolative linker as a versatile tool for detection of enzymes

Here we report on a dual-modal (19) F and (1) H MRI paramagnetic probe with a self-immolative linker, Gd-DOMF-Gal. The enzymatic conversion of this probe by β-galactosidase resulted in a simultaneous turning on of the fluorine signal and changed ability of the Gd(3+) complex to modulate the (1) H MR signal intensity of the surrounding water molecules. A versatile imaging platform for monitoring a variety of enzymes by (19) F and (1) H MRI using this molecular design is proposed.

Phase stability in fMRI time series: effect of noise regression, off-resonance correction and spatial filtering techniques

Although the majority of fMRI studies exploit magnitude changes only, there is an increasing interest regarding the potential additive information conveyed by the phase signal. This integrated part of the complex number furnished by the MR scanners can also be used for exploring direct detection of neuronal activity and for thermography. Few studies have explicitly addressed the issue of the available signal stability in the context of phase time-series, and therefore we explored the spatial pattern of frequency specific phase fluctuations, and evaluated the effect of physiological noise components (heart beat and respiration) on the phase signal. Three categories of retrospective noise reduction techniques were explored and the temporal signal stability was evaluated in terms of a physiologic noise model, for seven fMRI measurement protocols in eight healthy subjects at 3T, for segmented CSF, gray and white matter voxels. We confirmed that for most processing methods, an efficient use of the phase information is hampered by the fact that noise from physiological and instrumental sources contributes significantly more to the phase than to the magnitude instability. Noise regression based on the phase evolution of the central k-space point, RETROICOR, or an orthonormalized combination of these were able to reduce their impact, but without bringing phase stability down to levels expected from the magnitude signal. Similar results were obtained after targeted removal of scan-to-scan variations in the bulk magnetic field by the dynamic off-resonance in k-space (DORK) method and by the temporal off-resonance alignment of single-echo time series technique (TOAST). We found that spatial high-pass filtering was necessary, and in vivo a Gaussian filter width of 20mm was sufficient to suppress physiological noise and bring the phase fluctuations to magnitude levels. Stronger filters brought the fluctuations down to levels dictated by thermal noise contributions, and for 62.5mm(3) voxels the phase stability was as low as 5 mrad (0.27°). In conditions of low SNR(o) and high temporal sampling rate (short TR); we achieved an upper bound for the phase instabilities at 0.0017 ppm, which is close to the dHb contribution to the GM/WM phase contrast.

Smoothing that does not blur: effects of the anisotropic approach for evaluating diffusion tensor imaging data in the clinic

PURPOSE

To compare the effects of anisotropic and Gaussian smoothing on the outcomes of diffusion tensor imaging (DTI) voxel-based (VB) analyses in the clinic, in terms of signal-to-noise ratio (SNR) enhancement and directional information and boundary structures preservation.

MATERIALS AND METHODS

DTI data of 30 Alzheimer's disease (AD) patients and 30 matched control subjects were obtained at 3T. Fractional anisotropy (FA) maps with variable degrees and quality (Gaussian and anisotropic) of smoothing were created and compared with an unsmoothed dataset. The two smoothing approaches were evaluated in terms of SNR improvements, capability to separate differential effects between patients and controls by a standard VB analysis, and level of artifacts introduced by the preprocessing.

RESULTS

Gaussian smoothing regionally biased the FA values and introduced a high variability of results in clinical analysis, greatly dependent on the kernel size. On the contrary, anisotropic smoothing proved itself capable of enhancing the SNR of images and maintaining boundary structures, with only moderate dependence of results on smoothing parameters.

CONCLUSION

Our study suggests that anisotropic smoothing is more suitable in DTI studies; however, regardless of technique, a moderate level of smoothing seems to be preferable considering the artifacts introduced by this manipulation.

Metabolic correlatives of brain activity in a FOS epilepsy patient

The correlation and the interactions between neuronal activity and underlying metabolic dynamics are still a matter of debate, especially in pathological conditions. This study reports findings obtained on a subject suffering from fixation-off sensitivity (FOS) epilepsy, exploited as a model system of triggerable anomalous electrical activity. Functional Magnetic Resonance Spectroscopy was used to investigate the metabolic response to visual spike-inducing stimuli in a single voxel placed in the temporo-occipital lobe of a FOS epilepsy patient. MRS measurements were additionally performed on a control group of five healthy volunteers. The FOS patient also underwent an EEG session with the same stimulus paradigm. Uniquely in the FOS patient, glutamate and glutamine concentration increased during the first 10 min of stimulation and then returned to baseline. On the other hand, FOS-induced epileptic activity (spiking) endured throughout all the stimulation epoch. The observed metabolic dynamics may be likely linked to a complex interplay between alterations of the metabolic pathways of glutamate and modulation of the neuronal activity.

Quantification of gray matter changes in the cerebral cortex after isolated cerebellar damage: a voxel-based morphometry study

There is growing evidence based on behavioral and functional imaging studies about the cerebellar involvement in the modulation of cognitive functions. However, it still remains to be clarified how the cerebellum interacts with brain regions sub-serving different cognitive domains. In this study we used magnetic resonance imaging (MRI) and voxel based morphometry (VBM) to investigate changes of cerebral gray matter (GM) density in 15 patients with a focal cerebellar damage (CD) compared to 15 healthy controls. T2-weighted scans and T1-weighted volumes were collected from each subject. With the exception of the cerebellar lesion, none of the patients showed any additional brain MRI abnormality. T1-volumes were analyzed by voxel-based morphometry. Consistent with their neuropsychological abnormalities, patients with right-CD compared to controls showed a reduction of GM density mainly involving the left frontal, parietal and temporal lobes. Conversely, patients with left-CD did not show any significant neuropsychological or cerebral GM abnormality. The present study indicates that specific GM changes may be detected in patients with isolated CD and cognitive dysfunction. We discuss the findings in terms of cerebellar influence on the neuronal networks involved in higher level functions of the association cortex.

A highly sensitive radial diffusion measurement method for white matter tract investigation

This article describes a novel approach for local estimation of the radial diffusion coefficient (D(perpendicular)) in white matter (WM) regions containing well-oriented nervous fibers. The method is based on the assumption that the diffusion process in well-organized WM regions responds to a cylindrical symmetry. The increased precision in the estimation of D(perpendicular) provided by this local approach compared to standard techniques based on diffusion tensor imaging is demonstrated using numerical simulations. An in vivo validation of the technique is also provided, showing its application to the corpus callosum of six healthy volunteers, highlighting the sensitivity of the method. Assuming that D(perpendicular) is sensitive to myelin integrity, our technique has the potential to investigate pathophysiological aspects of several neurological and psychiatric disorders with improved precision in targeted WM tracts.

Essential head tremor is associated with cerebellar vermis atrophy: a volumetric and voxel-based morphometry MR imaging study

BACKGROUND AND PURPOSE

Our aim was to investigate the presence of brain gray matter (GM) abnormalities in patients with different forms of essential tremor (ET).

MATERIALS AND METHODS

We used optimized voxel-based morphometry (VBM) and manually traced single region-of-interest analysis in 50 patients with familial ET and in 32 healthy subjects. Thirty patients with ET had tremor of the arms (a-ET), whereas the remaining 20 patients had both arm and head tremor (h-ET).

RESULTS

VBM showed marked atrophy of the cerebellar vermis in the patients with h-ET with respect to healthy subjects (P(corrected) < .001). Patients with a-ET showed a trend toward a vermal GM volume loss that did not reach a significant difference with respect to healthy controls (P(uncorrected) < .01). The region-of-interest analysis showed a reduction of the cerebellar volume (CV) in the h-ET group (98.2 +/- 13.6 mm(3)) compared with healthy controls (110.5 +/- 15.5 mm(3), P < .012) as well as in the entire vermal area (790.3 +/- 94.5 mm(2), 898.6 +/- 170.6 mm(2), P < .04 in h-ET and control groups, respectively).

CONCLUSIONS

Atrophy of the cerebellar vermis detected in patients with h-ET strongly supports the evidence for the involvement of the cerebellum in the pathophysiology of ET. The lack of a significant CV loss observed in patients with a-ET suggests that a-ET and h-ET might represent distinct subtypes of the same disease.

The effect of physiological noise in phase functional magnetic resonance imaging: from blood oxygen level-dependent effects to direct detection of neuronal currents

Recently, the possibility to use both magnitude and phase image sets for the statistical evaluation of fMRI has been proposed, with the prospective of increasing both statistical power and the spatial specificity. In the present work, several issues that affect the spatial and temporal stability in fMRI phase time series in the presence of physiologic noise processes are reviewed, discussed and illustrated by experiments performed at 3 T. The observed phase value is a fingerprint of the underlying voxel averaged magnetic field variations. Those related to physiological processes can be considered static or dynamic in relation to the temporal scale of a 2D acquisition and will play out on different spatial scales as well: globally across the entire images slice, and locally depending on the constituents and their relative fractions inside the MRI voxel. The 'static' respiration-induced effects lead to magneto-mechanic scan-to-scan variations in the global magnetic field but may also contribute to local BOLD fluctuations due to respiration-related variations in arterial carbon dioxide. Likewise, the 'dynamic' cardiac-related effects will lead to global susceptibility effects caused by pulsatile motion of the brain as well as local blood pressure-related changes in BOLD and changes in blood flow velocity. Finally, subject motion may lead to variations in both local and global tissue susceptibility that will be especially pronounced close to air cavities. Since dissimilar manifestations of physiological processes can be expected in phase and in magnitude images, a direct relationship between phase and magnitude scan-to-scan fluctuations cannot be assumed a priori. Therefore three different models were defined for the phase stability, each dependent on the relation between phase and magnitude variations and the best will depend on the underlying noise processes. By experiments on healthy volunteers at rest, we showed that phase stability depends on the type of post-processing and can be improved by reducing the low-frequency respiration-induced mechano-magnetic effects. Although the manifestations of physiological noise were in general more pronounced in phase than in magnitude images, due to phase wraps and global Bo effects, we suggest that a phase stability similar to that found in magnitude could theoretically be achieved by adequate correction methods. Moreover, as suggested by our experimental data regarding BOLD-related phase effects, phase stability could even supersede magnitude stability in voxels covering dense microvascular networks with BOLD-related fluctuations as the dominant noise contributor. In the interest of the quality of both BOLD-based and nc-MRI methods, future studies are required to find alternative methods that can improve phase stability, designed to match the temporal and spatial scale of the underlying neuronal activity.

Realistic simulations of neuronal activity: a contribution to the debate on direct detection of neuronal currents by MRI

Many efforts have been done in order to preview the properties of the magnetic resonance (MR) signals produced by the neuronal currents using simulations. In this paper, starting with a detailed calculation of the magnetic field produced by the neuronal currents propagating over single hippocampal CA1 pyramidal neurons placed inside a cubic MR voxel of length 1.2 mm, we proceeded on the estimation of the phase and magnitude MR signals. We then extended the results to layers of parallel and synchronous similar neurons and to ensembles of layers, considering different echo times, voxel volumes and neuronal densities. The descriptions of the neurons and of their electrical activity took into account the real neuronal morphologies and the physiology of the neuronal events. Our results concern: (a) the expected time course of the MR signals produced by the neuronal currents in the brain, based on physiological and anatomical properties; (b) the different contributions of post-synaptic potentials and of action potentials to the MR signals; (c) the estimation of the equivalent current dipole and the influence of its orientation with respect to the external magnetic field on the observable MR signal variations; (d) the size of the estimated neuronal current induced phase and magnitude MR signal changes with respect to the echo time, voxel-size and neuronal density. The inclusion of realistic neuronal properties into the simulation introduces new information that can be helpful for the design of MR sequences for the direct detection of neuronal current effects and the testing of bio-electromagnetic models.

Model-free analysis of brain fMRI data by recurrence quantification

We propose a novel model-free univariate strategy for functional magnetic resonance imaging (fMRI) studies based upon recurrence quantification analysis (RQA). RQA is an auto-regressive method, which identifies recurrences in signals without any a priori assumptions. The performance of RQA is compared to that of univariate statistics based on a general linear model (GLM) and probabilistic independent component analysis (P-ICA) for a set of simulated and real fMRI data. RQA provides an appealing alternative to conventional GLM techniques, due to its exclusive feature of being model-free and of detecting potentially both linear and nonlinear dynamic processes, without requiring signal stationarity. The overall performance of the method compares positively also with P-ICA, another well-known model-free algorithm, which requires prior information to discriminate between different spatio-temporal processes. For simulated data, RQA is endowed with excellent accuracy for contrast-to-noise ratios greater than 0.2, and has a performance comparable to that of GLM for t(CNR)>or=0.8. For cerebral fMRI data acquired from a group of healthy subjects performing a finger-tapping task, (i) RQA reveals activations in the primary motor area contra-lateral to the employed hand and in the supplementary motor area, in agreement with the outcome of GLM analysis and (ii) identifies an additional brain region with transient signal changes. Moreover, RQA identifies signal recurrences induced by physiological processes other than BOLD (movement-related or of vascular origin). Finally, RQA is more robust than the GLM with respect to variations in the shape and timing of the underlying neuronal and hemodynamic responses which may vary between brain regions, subjects and tasks.

Multimodal fMRI tractography in normal subjects and in clinically recovered traumatic brain injury patients

In this study, we defined an operator-independent protocol for reconstructing the anatomical connections originating from fMRI activations in order to demonstrate that results obtained with this protocol are affected by alterations of functional activations. Seven healthy volunteers and two patients who sustained traumatic brain injury underwent an fMRI with a finger tapping task and a DTI scan. Cortical fMRI activations were used directly as seed mask for tractography for the reconstruction of individual motor pathways. On patients we observed a different motor network if compared to healthy subjects. However, when the activations of healthy subjects were used as seed masks for the tractography in patients, we observed for the patients a pattern of connectivity more similar to what was observed for healthy subjects. At the same time, when the activations of patients were used for the tractography on healthy subjects, we obtained patterns of connectivity similar to those obtained for patients. These results show the potential of the integration of fMRI and tractography for clarifying the mechanisms of cortical plasticity in the recovery of motor functions.

Implicit learning deficits in dyslexic adults: An fMRI study

It is assumed that several neuropsychological impairments characterize the cognitive profile of individuals with developmental dyslexia (DD). Phonological and visual processing are often impaired as well as auditory processing, attention, and information processing speed. Although reports in the literature on implicit learning abilities are contradictory, recent neurological and physiological data suggest that these abilities are deficient in individuals with DD. To evaluate implicit learning we administered a classical version of the serial reaction time task (SRTT) related to sequence learning. Using functional magnetic resonance imaging we investigated brain activation patterns associated with implicit learning deficits in 14 adults with DD matched with 14 normal readers. SRTT results indicated the absence of implicit learning in the DD group and different activations between groups mainly in SMA, inferior parietal areas and cerebellar lobule 6. These results can be interpreted in the light of the different capacities for the two groups to build an internal model to guide movements. Further, they explain DD individuals' difficulty in domains not directly related to reading ability.

Functional changes in the activity of cerebellum and frontostriatal regions during externally and internally timed movement in Parkinson's disease

We used fMRI to investigate the neurofunctional basis of externally and internally timed movements in Parkinson's disease (PD) patients. Ten PD patients whose medication had been withheld for at least 18h and 11 age- and sex-matched healthy controls were scanned while performing continuation paradigm with a visual metronome. Compared with the controls, PD patients displayed an intact capability to store and reproduce movement frequencies but with a significantly increased movement latencies. No differences in BOLD response were found in both groups when comparing the continuation with the preceding synchronization phase and viceversa, except for activity in visually related regions. Relative to healthy controls during the synchronization phase, PD patients exhibited an overall signal increase in the cerebellum and frontostriatal circuit (putamen, SMA and thalamus) activity together with specific brain areas (right inferior frontal gyrus and insula cortex) that are also implicated in primary timekeeper processes. By contrast, in the continuation phase the only neural network involved to a greater extent by the PD group was the cerebello-thalamic pathway. The lack of neurofunctional differences between the two timing phases suggests that rhythmic externally and internally guided movements engage similar neural networks in PD and matched healthy controls. Moreover, between-group comparison indicates that PD patients OFF medication may compensate for their basal ganglia-cortical loop's dysfunction using different motor pathways involving cerebellum and basal ganglia relays during the two phases of rhythmic movement.

Validation studies on the 5-hydroxy-L-[beta-11C]-tryptophan/PET method for probing the decarboxylase step in serotonin synthesis

The two-tissue compartment model, including irreversible trapping in the second compartment (2TCM) is used to describe the kinetics of 5-Hydroxy-L-[beta-(11)C]-tryptophan ([(11)C]HTP), a radioligand used in positron emission tomography (PET) for probing the second enzymatic step in the biosynthesis of serotonin. In this study, we examined the capacity of the model to track pharmacological changes in this biological process. We also investigated the potential loss of [(11)C]HTP-derived radioactivity during a PET study, since loss should be negligible not to alter quantification. Six rhesus monkeys were investigated using bolus [(11)C]HTP/PET methodology before and after pharmacological intervention. The second enzymatic step in serotonin synthesis was inhibited using the aromatic L-amino acid decarboxylase inhibitor NSD1015 (10 mg/kg). The extent of [(11)C]-derived radioactivity loss from the brain was studied by inhibition of the enzyme responsible for formation of the tissue metabolite, monoamine oxidase A, using clorgyline (2 mg/kg). After NSD1015, the uptake of [(11)C]HTP-derived radioactivity was increased in all the investigated brain regions, while the parameter used to reflect decarboxylase activity, the net accumulation rate constant (K(acc)), was decreased by 37% in the striatum, compared with baseline. Pretreatment with clorgyline did not change the brain uptake of [(11)C]HTP-derived radioactivity or K(acc). This study demonstrates that the 2TCM for [(11)C]HTP/PET is able to detect changes occurring during alteration of the biological process (i.e., the conversion of HTP to serotonin). Elimination of the radiotracer metabolite [(11)C]HIAA from the brain may be considered negligible if the PET study is limited to 60 min.

Imaging nervous pathways with MR tractography

This paper describes the state of the art of tractography, a technique which enables the virtual reconstruction of axon bundles of the central nervous system using diffusion-weighted magnetic resonance images. This technique has raised enthusiasm and expectations among specialists because it is the only non-invasive method for studying the three-dimensional architecture of axonal fibres in vivo. Tractography is a new technique used to assess the anatomy of the central nervous system, and it will be available for routine clinical use in the future. Understanding its potential applications and limitations is therefore important.

Challenges for detection of neuronal currents by MRI

Neuronal current MRI (nc-MRI) is an imaging method that directly maps magnetic field changes caused by neuronal currents with, at the same time, a high spatial and temporal resolution. A viable nc-MRI method would be of great benefit, both for the study of human brain function and for clinical applications in the field of epilepsy, especially for the noninvasive presurgical mapping of epileptogenic foci. A survey of fundamental issues in nc-MRI is reviewed, and challenges for future developments of the method are described within this context. Particularly, an overview of the models for signal generation is given, and the origin and physiology of different sources of neuronal currents are described. Prospects for predicting neuronal currents by electromagnetic field mapping and using this information, both a priori and a posteriori, for nc-MRI are considered. Ways of increasing specificity in nc-MRI by minimizing secondary hemodynamic and metabolic effects are described as well as means of optimizing the nc-MRI method for pushing the detection limit. Previously published works are described within these categories and future directions for nc-MRI are proposed.

Simultaneous EEG-fMRI acquisition: how far is it from being a standardized technique?

Simultaneous EEG-fMRI is a powerful tool to study spontaneous and evoked brain activity because of the complementary advantages of the two techniques in terms of temporal and spatial resolution. In recent years, a significant number of scientific works have been published on this subject. However, many technical problems related to the intrinsic incompatibility of EEG and MRI methods are still not fully solved. Furthermore, simultaneous acquisition of EEG and event-related fMRI requires precise synchronization of all devices involved in the experimental setup. Thus, timing issue must be carefully considered in order to avoid significant methodological errors. The aim of the present work is to highlight and discuss some of technical and methodological open issues associated with the combined use of EEG and fMRI. These issues are presented in the context of preliminary data regarding simultaneous acquisition of event-related evoked potentials and BOLD images during a visual odd-ball paradigm.

Fast detection of diffuse axonal damage in severe traumatic brain injury: comparison of gradient-recalled echo and turbo proton echo-planar spectroscopic imaging MRI sequences

BACKGROUND AND PURPOSE

Diffuse axonal injury (DAI) is a common type of primary neuronal injury in patients with severe traumatic brain injury (TBI), and is frequently accompanied by tissue tear hemorrhage. T2-weighted gradient-recalled echo (GRE) sequences are more sensitive than T2-weighted spin-echo images for detection of hemorrhage. The purpose of this study is to compare turbo Proton Echo Planar Spectroscopic Imaging (t-PEPSI), an extremely fast sequence, with GRE sequence in the detection of DAI.

METHODS

Twenty-one patients (mean age 26.8 years) with severe TBI occurred at least 3 months earlier, underwent a brain MR Imaging study on a 1.5-T scanner. A qualitative evaluation of the t-PEPSI sequences was performed by identifying the optimal echo time and in-plane resolution. The number and size of DAI lesions, as well as the signal intensity contrast ratio (SI CR), were computed for each set of GRE and t-PEPSI images, and divided according to their anatomic location as lobar and/or deep brain.

RESULTS

There was no significant difference between GRE and t-PEPSI sequences in the detection of the total number of DAI lesions (291 vs. 230, respectively). GRE sequence delineated a higher number of DAI in the temporal lobe compared to the t-PEPSI sequence (74 vs. 37, P < .004), while no differences were found for the other regions. The SI CR was significantly lower with the t-PEPSI than the GRE sequence (P < .00001).

CONCLUSION

Owing to its very short scan time and high sensitivity to the hemorrhage foci, the t-PEPSI sequence may be used as an alternative to the GRE to assess brain DAI in severe TBI patients, especially if uncooperative and medically unstable.