Comparison of diffusion MRI and CLARITY fiber orientation estimates in both gray and white matter regions of human and primate brain.

Diffusion MRI (dMRI) represents one of the few methods for mapping brain fiber orientations non-invasively. Unfortunately, dMRI fiber mapping is an indirect method that relies on inference from measured diffusion patterns. Comparing dMRI results with other modalities is a way to improve the interpretation of dMRI data and help advance dMRI technologies. Here, we present methods for comparing dMRI fiber orientation estimates with optical imaging of fluorescently labeled neurofilaments and vasculature in 3D human and primate brain tissue cuboids cleared using CLARITY. The recent advancements in tissue clearing provide a new opportunity to histologically map fibers projecting in 3D, which represents a captivating complement to dMRI measurements. In this work, we demonstrate the capability to directly compare dMRI and CLARITY in the same human brain tissue and assess multiple approaches for extracting fiber orientation estimates from CLARITY data. We estimate the three-dimensional neuronal fiber and vasculature orientations from neurofilament and vasculature stained CLARITY images by calculating the tertiary eigenvector of structure tensors. We then extend CLARITY orientation estimates to an orientation distribution function (ODF) formalism by summing multiple sub-voxel structure tensor orientation estimates. In a sample containing part of the human thalamus, there is a mean angular difference of 19o±15o between the primary eigenvectors of the dMRI tensors and the tertiary eigenvectors from the CLARITY neurofilament stain. We also demonstrate evidence that vascular compartments do not affect the dMRI orientation estimates by showing an apparent lack of correspondence (mean angular difference = 49o±23o) between the orientation of the dMRI tensors and the structure tensors in the vasculature stained CLARITY images. In a macaque brain dataset, we examine how the CLARITY feature extraction depends on the chosen feature extraction parameters. By varying the volume of tissue over which the structure tensor estimates are derived, we show that orientation estimates are noisier with more spurious ODF peaks for sub-voxels below 30 µm3 and that, for our data, the optimal gray matter sub-voxel size is between 62.5 µm3 and 125 µm3. The example experiments presented here represent an important advancement towards robust multi-modal MRI-CLARITY comparisons.

Pushing the limits of in vivo diffusion MRI for the Human Connectome Project.

Perhaps more than any other "-omics" endeavor, the accuracy and level of detail obtained from mapping the major connection pathways in the living human brain with diffusion MRI depend on the capabilities of the imaging technology used. The current tools are remarkable; allowing the formation of an "image" of the water diffusion probability distribution in regions of complex crossing fibers at each of half a million voxels in the brain. Nonetheless our ability to map the connection pathways is limited by the image sensitivity and resolution, and also the contrast and resolution in encoding of the diffusion probability distribution. The goal of our Human Connectome Project (HCP) is to address these limiting factors by re-engineering the scanner from the ground up to optimize the high b-value, high angular resolution diffusion imaging needed for sensitive and accurate mapping of the brain's structural connections. Our efforts were directed based on the relative contributions of each scanner component. The gradient subsection was a major focus since gradient amplitude is central to determining the diffusion contrast, the amount of T2 signal loss, and the blurring of the water PDF over the course of the diffusion time. By implementing a novel 4-port drive geometry and optimizing size and linearity for the brain, we demonstrate a whole-body sized scanner with G(max) = 300 mT/m on each axis capable of the sustained duty cycle needed for diffusion imaging. The system is capable of slewing the gradient at a rate of 200 T/m/s as needed for the EPI image encoding. In order to enhance the efficiency of the diffusion sequence we implemented a FOV shifting approach to Simultaneous MultiSlice (SMS) EPI capable of unaliasing 3 slices excited simultaneously with a modest g-factor penalty allowing us to diffusion encode whole brain volumes with low TR and TE. Finally we combine the multi-slice approach with a compressive sampling reconstruction to sufficiently undersample q-space to achieve a DSI scan in less than 5 min. To augment this accelerated imaging approach we developed a 64-channel, tight-fitting brain array coil and show its performance benefit compared to a commercial 32-channel coil at all locations in the brain for these accelerated acquisitions. The technical challenges of developing the over-all system are discussed as well as results from SNR comparisons, ODF metrics and fiber tracking comparisons. The ultra-high gradients yielded substantial and immediate gains in the sensitivity through reduction of TE and improved signal detection and increased efficiency of the DSI or HARDI acquisition, accuracy and resolution of diffusion tractography, as defined by identification of known structure and fiber crossing.

An implanted 8-channel array coil for high-resolution macaque MRI at 3T.

An 8-channel receive coil array was constructed and implanted adjacent to the skull in a male rhesus monkey in order to improve the sensitivity of (functional) brain imaging. The permanent implant was part of an acrylic headpost assembly and only the coil element loop wires were implanted. The tuning, matching, and preamplifier circuitry was connected via a removable external assembly. Signal-to-noise ratio (SNR) and noise amplification for parallel imaging were compared to single-, 4-, and 8-channel external receive-only coils routinely used for macaque fMRI. In vivo measurements showed significantly improved SNR within the brain for the implanted versus the external coils. Within a region-of-interest covering the cerebral cortex, we observed a 5.4-, 3.6-fold, and 3.4-fold increase in SNR compared to the external single-, 4-, and 8-channel coils, respectively. In the center of the brain, the implanted array maintained a 2.4×, 2.5×, and 2.1× higher SNR, respectively compared to the external coils. The array performance was evaluated for anatomical, diffusion tensor and functional brain imaging. This study suggests that a stable implanted phased-array coil can be used in macaque MRI to substantially increase the spatial resolution for anatomical, diffusion tensor, and functional imaging.

T2* mapping and B0 orientation-dependence at 7 T reveal cyto- and myeloarchitecture organization of the human cortex.

Ultra-high field MRI (≥ 7 T) has recently shown great sensitivity to depict patterns of tissue microarchitecture. Moreover, recent studies have demonstrated a dependency between T2* and orientation of white matter fibers with respect to the main magnetic field B0. In this study we probed the potential of T2* mapping at 7 T to provide new markers of cortical architecture. We acquired multi-echo measurements at 7 T and mapped T2* over the entire cortex of eight healthy individuals using surface-based analysis. B0 dependence was tested by computing the angle θ(z) between the normal of the surface and the direction of B0, then fitting T2*(θ(z)) using model from the literature. Average T2* in the cortex was 32.20 +/- 1.35 ms. Patterns of lower T2* were detected in the sensorimotor, visual and auditory cortices, likely reflecting higher myelin content. Significantly lower T2* was detected in the left hemisphere of the auditory region (p<0.005), suggesting higher myelin content, in accordance with previous investigations. B0 orientation dependence was detected in some areas of the cortex, the strongest being in the primary motor cortex (∆R2*=4.10 Hz). This study demonstrates that quantitative T2* measures at 7 T MRI can reveal patterns of cytoarchitectural organization of the human cortex in vivo and that B0 orientation dependence can probe the coherency and orientation of gray matter fibers in the cortex, shedding light into the potential use of this type of contrast to characterize cyto-/myeloarchitecture and to understand the pathophysiology of diseases associated with changes in iron and/or myelin concentration.

Tissue oxygen tension measurements in the Shionogi model of prostate cancer using 19F MRS and MRI.

OBJECTIVES: To investigate changes in tumour tissue oxygenation throughout the tumour growth-regression-relapse cycle in an androgen-dependent animal tumour model. MATERIALS AND METHODS: 19F T1 relaxometry of Perfluoro-15-Crown-5-Ether was used to measure in vivo partial oxygen pressure (pO2) of Shionogi tumours on a 2.35-T MR scanner. Perfluoro-15-Crown-5-Ether was administered as an emulsion injected intravenously or as a neat compound injected directly into the tumour. Non-localized, tumour 19F T1 measurements, made at multiple time points throughout the tumour cycle, were translated into pO2 levels. RESULTS: No correlation between tumour size and pO2 values was found. Values of pO2 for growing tumours (50 +/- 30 torr) were significantly lower than for regressing and relapsing tumours after 9 days post-castration (70 +/- 10 torr, p<0.05). Maximum pO2 values (90 +/- 30 torr) were reached between fifth and eighth day post-castration, when tumour pO2 was significantly higher than both pre-castration (p<0.001) and after 9 days post-castration (p<0.05). CONCLUSION: We demonstrate that longitudinal pO2 measurements in vivo are feasible. Values of pO2 for growing androgen-dependent tumours were significantly lower than for regressing and relapsing androgen-independent tumours. These results have potential clinical importance in optimizing the timing of chemotherapy and/or radiotherapy of hormone dependent tumours.