Improving in vivo human cerebral cortical surface reconstruction using data-driven super-resolution.

Accurate and automated reconstruction of the in vivo human cerebral cortical surface from anatomical magnetic resonance (MR) images facilitates the quantitative analysis of cortical structure. Anatomical MR images with sub-millimeter isotropic spatial resolution improve the accuracy of cortical surface and thickness estimation compared to the standard 1-millimeter isotropic resolution. Nonetheless, sub-millimeter resolution acquisitions require averaging multiple repetitions to achieve sufficient signal-to-noise ratio and are therefore long and potentially vulnerable to subject motion. We address this challenge by synthesizing sub-millimeter resolution images from standard 1-millimeter isotropic resolution images using a data-driven supervised machine learning-based super-resolution approach achieved via a deep convolutional neural network. We systematically characterize our approach using a large-scale simulated dataset and demonstrate its efficacy in empirical data. The super-resolution data provide improved cortical surfaces similar to those obtained from native sub-millimeter resolution data. The whole-brain mean absolute discrepancy in cortical surface positioning and thickness estimation is below 100 µm at the single-subject level and below 50 µm at the group level for the simulated data, and below 200 µm at the single-subject level and below 100 µm at the group level for the empirical data, making the accuracy of cortical surfaces derived from super-resolution sufficient for most applications.

Resting-state "physiological networks".

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

Axon diameter index estimation independent of fiber orientation distribution using high-gradient diffusion MRI.

Axon diameter mapping using high-gradient diffusion MRI has generated great interest as a noninvasive tool for studying trends in axonal size in the human brain. One of the main barriers to mapping axon diameter across the whole brain is accounting for complex white matter fiber configurations (e.g., crossings and fanning), which are prevalent throughout the brain. Here, we present a framework for generalizing axon diameter index estimation to the whole brain independent of the underlying fiber orientation distribution using the spherical mean technique (SMT). This approach is shown to significantly benefit from the use of real-valued diffusion data with Gaussian noise, which reduces the systematic bias in the estimated parameters resulting from the elevation of the noise floor when using magnitude data with Rician noise. We demonstrate the feasibility of obtaining whole-brain orientationally invariant estimates of axon diameter index and relative volume fractions in six healthy human volunteers using real-valued diffusion data acquired on a dedicated high-gradient 3-Tesla human MRI scanner with 300 mT/m maximum gradient strength. The trends in axon diameter index are consistent with known variations in axon diameter from histology and demonstrate the potential of this generalized framework for revealing coherent patterns in axonal structure throughout the living human brain. The use of real-valued diffusion data provides a viable solution for eliminating the Rician noise floor and should be considered for all spherical mean approaches to microstructural parameter estimation.

Age-related alterations in axonal microstructure in the corpus callosum measured by high-gradient diffusion MRI.

Cerebral white matter exhibits age-related degenerative changes during the course of normal aging, including decreases in axon density and alterations in axonal structure. Noninvasive approaches to measure these microstructural alterations throughout the lifespan would be invaluable for understanding the substrate and regional variability of age-related white matter degeneration. Recent advances in diffusion magnetic resonance imaging (MRI) have leveraged high gradient strengths to increase sensitivity toward axonal size and density in the living human brain. Here, we examined the relationship between age and indices of axon diameter and packing density using high-gradient strength diffusion MRI in 36 healthy adults (aged 22-72) in well-defined central white matter tracts in the brain. A recently validated method for inferring the effective axonal compartment size and packing density from diffusion MRI measurements acquired with 300 mT/m maximum gradient strength was applied to the in vivo human brain to obtain indices of axon diameter and density in the corpus callosum, its sub-regions, and adjacent anterior and posterior fibers in the forceps minor and forceps major. The relationships between the axonal metrics, corpus callosum area and regional gray matter volume were also explored. Results revealed a significant increase in axon diameter index with advancing age in the whole corpus callosum. Similar analyses in sub-regions of the corpus callosum showed that age-related alterations in axon diameter index and axon density were most pronounced in the genu of the corpus callosum and relatively absent in the splenium, in keeping with findings from previous histological studies. The significance of these correlations was mirrored in the forceps minor and forceps major, consistent with previously reported decreases in FA in the forceps minor but not in the forceps major with age. Alterations in the axonal imaging metrics paralleled decreases in corpus callosum area and regional gray matter volume with age. Among older adults, results from cognitive testing suggested an association between larger effective compartment size in the corpus callosum, particularly within the genu of the corpus callosum, and lower scores on the Montreal Cognitive Assessment, largely driven by deficits in short-term memory. The current study suggests that high-gradient diffusion MRI may be sensitive to the axonal substrate of age-related white matter degeneration reflected in traditional DTI metrics and provides further evidence for regionally selective alterations in white matter microstructure with advancing age.

Highly-accelerated volumetric brain examination using optimized wave-CAIPI encoding.

BACKGROUND: Rapid volumetric imaging protocols could better utilize limited scanner resources. PURPOSE: To develop and validate an optimized 6-minute high-resolution volumetric brain MRI examination using Wave-CAIPI encoding. STUDY TYPE: Prospective. POPULATION/SUBJECTS: Ten healthy subjects and 20 patients with a variety of intracranial pathologies. FIELD STRENGTH/SEQUENCE: At 3 T, MPRAGE, T2 -weighted SPACE, SPACE FLAIR, and SWI were acquired at 9-fold acceleration using Wave-CAIPI and for comparison at 2-4-fold acceleration using conventional GRAPPA. ASSESSMENT: Extensive simulations were performed to optimize the Wave-CAIPI protocol and minimize both g-factor noise amplification and potential T1 /T2 blurring artifacts. Moreover, refinements in the autocalibrated reconstruction of Wave-CAIPI were developed to ensure high-quality reconstructions in the presence of gradient imperfections. In a randomized and blinded fashion, three neuroradiologists assessed the diagnostic quality of the optimized 6-minute Wave-CAIPI exam and compared it to the roughly 3× slower GRAPPA accelerated protocol using both an individual and head-to-head analysis. STATISTICAL TEST: A noninferiority test was used to test whether the diagnostic quality of Wave-CAIPI was noninferior to the GRAPPA acquisition, with a 15% noninferiority margin. RESULTS: Among all sequences, Wave-CAIPI achieved negligible g-factor noise amplification (gavg ≤ 1.04) and burring artifacts from T1 /T2 relaxation. Improvements of our autocalibration approach for gradient imperfections enabled increased robustness to gradient mixing imperfections in tilted-field of view (FOV) prescriptions as well as variations in gradient and analog-to-digital converter (ADC) sampling rates. In the clinical evaluation, Wave-CAIPI achieved similar mean scores when compared with GRAPPA (MPRAGE: ØW = 4.03, ØG = 3.97; T2 w SPACE: ØW = 4.00, ØG = 4.00; SPACE FLAIR: ØW = 3.97, ØG = 3.97; SWI: ØW = 3.93, ØG = 3.83) and was statistically noninferior (N = 30, P < 0.05 for all sequences). DATA CONCLUSION: The proposed volumetric brain exam retained comparable image quality when compared with the much longer conventional protocol. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:961-974.

Motion-robust sub-millimeter isotropic diffusion imaging through motion corrected generalized slice dithered enhanced resolution (MC-gSlider) acquisition.

PURPOSE: To develop an efficient MR technique for ultra-high resolution diffusion MRI (dMRI) in the presence of motion. METHODS: gSlider is an SNR-efficient high-resolution dMRI acquisition technique. However, subject motion is inevitable during a prolonged scan for high spatial resolution, leading to potential image artifacts and blurring. In this study, an integrated technique termed Motion Corrected gSlider (MC-gSlider) is proposed to obtain high-quality, high-resolution dMRI in the presence of large in-plane and through-plane motion. A motion-aware reconstruction with spatially adaptive regularization is developed to optimize the conditioning of the image reconstruction under difficult through-plane motion cases. In addition, an approach for intra-volume motion estimation and correction is proposed to achieve motion correction at high temporal resolution. RESULTS: Theoretical SNR and resolution analysis validated the efficiency of MC-gSlider with regularization, and aided in selection of reconstruction parameters. Simulations and in vivo experiments further demonstrated the ability of MC-gSlider to mitigate motion artifacts and recover detailed brain structures for dMRI at 860 µm isotropic resolution in the presence of motion with various ranges. CONCLUSION: MC-gSlider provides motion-robust, high-resolution dMRI with a temporal motion correction sensitivity of 2 s, allowing for the recovery of fine detailed brain structures in the presence of large subject movements.

Comprehensive diffusion MRI dataset for in vivo human brain microstructure mapping using 300 mT/m gradients.

Strong gradient systems can improve the signal-to-noise ratio of diffusion MRI measurements and enable a wider range of acquisition parameters that are beneficial for microstructural imaging. We present a comprehensive diffusion MRI dataset of 26 healthy participants acquired on the MGH-USC 3 T Connectome scanner equipped with 300 mT/m maximum gradient strength and a custom-built 64-channel head coil. For each participant, the one-hour long acquisition systematically sampled the accessible diffusion measurement space, including two diffusion times (19 and 49 ms), eight gradient strengths linearly spaced between 30 mT/m and 290 mT/m for each diffusion time, and 32 or 64 uniformly distributed directions. The diffusion MRI data were preprocessed to correct for gradient nonlinearity, eddy currents, and susceptibility induced distortions. In addition, scan/rescan data from a subset of seven individuals were also acquired and provided. The MGH Connectome Diffusion Microstructure Dataset (CDMD) may serve as a test bed for the development of new data analysis methods, such as fiber orientation estimation, tractography and microstructural modelling.