DR-BUDDI (Diffeomorphic Registration for Blip-Up blip-Down Diffusion Imaging) method for correcting echo planar imaging distortions

PyHySCO: GPU-enabled susceptibility artifact distortion correction in seconds

Over the past decade, reversed gradient polarity (RGP) methods have become a popular approach for correcting susceptibility artifacts in echo-planar imaging (EPI). Although several post-processing tools for RGP are available, their implementations do not fully leverage recent hardware, algorithmic, and computational advances, leading to correction times of several minutes per image volume. To enable 3D RGP correction in seconds, we introduce PyTorch Hyperelastic Susceptibility Correction (PyHySCO), a user-friendly EPI distortion correction tool implemented in PyTorch that enables multi-threading and efficient use of graphics processing units (GPUs). PyHySCO uses a time-tested physical distortion model and mathematical formulation and is, therefore, reliable without training. An algorithmic improvement in PyHySCO is its use of the one-dimensional distortion correction method by Chang and Fitzpatrick to initialize the non-linear optimization. PyHySCO is published under the GNU public license and can be used from the command line or its Python interface. Our extensive numerical validation using 3T and 7T data from the Human Connectome Project suggests that PyHySCO can achieve accuracy comparable to that of leading RGP tools at a fraction of the cost. We also validate the new initialization scheme, compare different optimization algorithms, and test the algorithm on different hardware and arithmetic precisions.

Arterial pulsation dependence of perivascular cerebrospinal fluid flow measured by dynamic diffusion tensor imaging in the human brain

Perivascular cerebrospinal fluid (pCSF) flow is a key component of the glymphatic system. Arterial pulsation has been proposed as the main driving force of pCSF influx along the superficial and penetrating arteries; however, evidence of this mechanism in humans is limited. We proposed an experimental framework of dynamic diffusion tensor imaging with low b-values and ultra-long echo time (dynDTIlow-b) to capture pCSF flow properties during the cardiac cycle in human brains. Healthy adult volunteers (aged 17-28 years; seven men, one woman) underwent dynDTIlow-b using a clinical 3T scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) with simultaneously recorded cardiac output. The results showed that diffusion tensors reconstructed from pCSF were mainly oriented in the direction of the neighboring arterial flow. When switching from vasoconstriction to vasodilation, the axial and radial diffusivities of the pCSF increased by 5.7% and 4.94%, respectively, suggesting that arterial pulsation alters the pCSF flow both parallel and perpendicular to the arterial wall. DynDTIlow-b signal intensity at b=0 s/mm2 (i.e., T2-weighted, [S(b=0 s/mm2)]) decreased in systole, but this change was ∼7.5% of a cardiac cycle slower than the changes in apparent diffusivity, suggesting that changes in S(b=0 s/mm2) and apparent diffusivity arise from distinct physiological processes and potential biomarkers associated with perivascular space volume and pCSF flow, respectively. Additionally, the mean diffusivities of white matter showed cardiac-cycle dependencies similar to pCSF, although a delay relative to the peak time of S(b=0 s/mm2) was present, suggesting that dynDTIlow-b could potentially reveal the dynamics of magnetic resonance imaging-invisible pCSF surrounding small arteries and arterioles in white matter; this delay may result from pulse wave propagation along penetrating arteries. In conclusion, the vasodilation-induced increases in axial and radial diffusivities of pCSF and mean diffusivities of white matter are consistent with the notion that arterial pulsation can accelerate pCSF flow in human brain. Furthermore, the proposed dynDTIlow-b technique can capture various pCSF dynamics in artery pulsation.

Influence of preprocessing, distortion correction and cardiac triggering on the quality of diffusion MR images of spinal cord

Diffusion MRI of the spinal cord (SC) is susceptible to geometric distortion caused by field inhomogeneities, and prone to misalignment across time series and signal dropout caused by biological motion. Several modifications of image acquisition and image processing techniques have been introduced to overcome these artifacts, but their specific benefits are largely unproven and warrant further investigations. We aim to evaluate two specific aspects of image acquisition and processing that address image quality in diffusion studies of the spinal cord: susceptibility corrections to reduce geometric distortions, and cardiac triggering to minimize motion artifacts. First, we evaluate 4 distortion preprocessing strategies on 7 datasets of the cervical and lumbar SC and find that while distortion correction techniques increase geometric similarity to structural images, they are largely driven by the high-contrast cerebrospinal fluid, and do not consistently improve the geometry within the cord nor improve white-to-gray matter contrast. We recommend at a minimum to perform bulk-motion correction in preprocessing and posit that improvements/adaptations are needed for spinal cord distortion preprocessing algorithms, which are currently optimized and designed for brain imaging. Second, we design experiments to evaluate the impact of removing cardiac triggering. We show that when triggering is foregone, images are qualitatively similar to triggered sequences, do not have increased prevalence of artifacts, and result in similar diffusion tensor indices with similar reproducibility to triggered acquisitions. When triggering is removed, much shorter acquisitions are possible, which are also qualitatively and quantitatively similar to triggered sequences. We suggest that removing cardiac triggering for cervical SC diffusion can be a reasonable option to save time with minimal sacrifice to image quality.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Distinct alterations in white matter properties and organization related to maternal treatment initiation in neonates exposed to HIV but uninfected

HIV exposed-uninfected (HEU) infants and children are at risk of developmental delays as compared to HIV uninfected unexposed (HUU) populations. The effects of exposure to in utero HIV and ART regimens on the HEU the developing brain are not well understood. In a cohort of 2-week-old newborns, we used diffusion tensor imaging (DTI) tractography and graph theory to examine the influence of HIV and ART exposure in utero on neonate white matter integrity and organisation. The cohort included HEU infants born to mothers who started ART before conception (HEUpre) and after conception (HEUpost), as well as HUU infants from the same community. We investigated HIV exposure and ART duration group differences in DTI metrics (fractional anisotropy (FA) and mean diffusivity (MD)) and graph measures across white matter. We found increased MD in white matter connections involving the thalamus and limbic system in the HEUpre group compared to HUU. We further identified reduced nodal efficiency in the basal ganglia. Within the HEUpost group, we observed reduced FA in cortical-subcortical and cerebellar connections as well as decreased transitivity in the hindbrain area compared to HUU. Overall, our analysis demonstrated distinct alterations in white matter integrity related to the timing of maternal ART initiation that influence regional brain network properties.

A framework for optimizing the acquisition protocol multishell diffusion-weighted imaging for multimodel assessment

Complementary aspects of tissue microstructure can be studied with diffusion-weighted imaging (DWI). However, there is no consensus on how to design a diffusion acquisition protocol for multiple models within a clinically feasible time. The purpose of this study is to provide a flexible framework that is able to optimize the shell acquisition protocol given a set of DWI models. Eleven healthy subjects underwent an extensive DWI acquisition protocol, including 15 candidate shells, ranging from 10 to 3500 s/mm2. The proposed framework aims to determine the optimized acquisition scheme (OAS) with a data-driven procedure minimizing the squared error of model-estimated parameters. We tested the proposed method over five heterogeneous DWI models exploiting both low and high b-values (i.e., diffusion tensor imaging [DTI], free water, intra-voxel incoherent motion [IVIM], diffusion kurtosis imaging [DKI], and neurite orientation dispersion and density imaging [NODDI]). A voxel-level and region of interest (ROI)-level analysis was conducted over the white matter and in 48 fiber bundles, respectively. Results showed that acquiring data for the five abovementioned models via OAS requires 14 min, compared with 35 min for the joint recommended acquisition protocol. The parameters derived from the reference acquisition scheme and the OAS are comparable in terms of estimated values, noise, and tissue contrast. Furthermore, the power analysis showed that the OAS retains the potential sensitivity to group-level differences in the parameters of interest, with the exception of the free water model. Overall, there is a linear correspondence (R2 = 0.91) between OAS and reference-derived parameters. In conclusion, the proposed framework optimizes the shell acquisition scheme for a given set of DWI models (i.e., DTI, free water, IVIM, DKI, and NODDI), combining low and high b-values while saving acquisition time.

Evolutionary continuity and divergence of auditory dorsal and ventral pathways in primates revealed by ultra-high field diffusion MRI

Auditory dorsal and ventral pathways in the human brain play important roles in supporting speech and language processing. However, the evolutionary root of the dual auditory pathways in the primate brain is unclear. By parcellating the auditory cortex of marmosets (a New World monkey species), macaques (an Old World monkey species), and humans using the same individual-based analysis method and tracking the pathways from the auditory cortex based on multi-shell diffusion-weighted MRI (dMRI), homologous auditory dorsal and ventral fiber tracks were identified in these primate species. The ventral pathway was found to be well conserved in all three primate species analyzed but extend to more anterior temporal regions in humans. In contrast, the dorsal pathway showed a divergence between monkey and human brains. First, frontal regions in the human brain have stronger connections to the higher-level auditory regions than to the lower-level auditory regions along the dorsal pathway, while frontal regions in the monkey brain show opposite connection patterns along the dorsal pathway. Second, the left lateralization of the dorsal pathway is only found in humans. Moreover, the connectivity strength of the dorsal pathway in marmosets is more similar to that of humans than macaques. These results demonstrate the continuity and divergence of the dual auditory pathways in the primate brains along the evolutionary path, suggesting that the putative neural networks supporting human speech and language processing might have emerged early in primate evolution.

Comparison of distortion correction preprocessing pipelines for DTI in the upper limb

PURPOSE

DTI characterizes tissue microstructure and provides proxy measures of nerve health. Echo-planar imaging is a popular method of acquiring DTI but is susceptible to various artifacts (e.g., susceptibility, motion, and eddy currents), which may be ameliorated via preprocessing. There are many pipelines available but limited data comparing their performance, which provides the rationale for this study.

METHODS

DTI was acquired from the upper limb of heathy volunteers at 3T in blip-up and blip-down directions. Data were independently corrected using (i) FSL's TOPUP & eddy, (ii) FSL's TOPUP, (iii) DSI Studio, and (iv) TORTOISE. DTI metrics were extracted from the median, radial, and ulnar nerves and compared (between pipelines) using mixed-effects linear regression. The geometric similarity of corrected b = 0 images and the slice matched T1-weighted (T1w) images were computed using the Sörenson-Dice coefficient.

RESULTS

Without preprocessing, the similarity coefficient of the blip-up and blip-down datasets to the T1w was 0·80 and 0·79, respectively. Preprocessing improved the geometric similarity by 1% with no difference between pipelines. Compared to TOPUP & eddy, DSI Studio and TORTOISE generated 2% and 6% lower estimates of fractional anisotropy, and 6% and 13% higher estimates of radial diffusivity, respectively. Estimates of anisotropy from TOPUP & eddy versus TOPUP were not different but TOPUP reduced radial diffusivity by 3%. The agreement of DTI metrics between pipelines was poor.

CONCLUSIONS

Preprocessing DTI from the upper limb improves geometric similarity but the choice of the pipeline introduces clinically important variability in diffusion parameter estimates from peripheral nerves.

Fronto-limbic white matter microstructural changes in psychiatrically healthy adults with childhood trauma

Childhood trauma (CT) may influence brain white matter microstructure; however, few studies have examined the differential impact of distinct CT types on white matter microstructure in psychiatrically healthy adults living in a developing country. In adults without significant medical or psychiatric disorders, we investigated the association(s) between CT, including abuse and neglect, and fractional anisotropy (FA) of limbic tracts previously shown to be associated with CT. Participants underwent diffusion tensor imaging and completed the Childhood Trauma Questionnaire. Multivariate analysis of variance models were used to test the effects of total overall CT, as well as CT subtypes, on FA in six fronto-limbic tracts, adjusting for age, sex, and educational level. The final sample included 69 adults (age 47 ± 17 years; 70% female). Overall, CT had a significant main effect on FA for tracts of interest (p < .001). Greater CT severity was associated with lower FA for the bilateral and left stria terminalis (uncorrected) as well as the bilateral, left, and right anterior limb of the internal capsule (ALIC; corrected). Exposure to total non-violent/deprivational trauma specifically was associated with lower FA of the bilateral, left, and right ALIC, suggesting that distinct types of CT are associated with differential white matter changes in apparently healthy adults. The ALIC predominantly carries fibers connecting the thalamus with prefrontal cortical regions. Microstructural alterations in the ALIC may be associated with functional brain changes, which may be adaptive or increase the risk of accelerated age-related cognitive decline, maladaptive behaviors, and subsyndromal psychiatric symptoms.

Role of autonomic, nociceptive, and limbic brainstem nuclei in core autism features

Although multiple theories have speculated about the brainstem reticular formation's involvement in autistic behaviors, the in vivo imaging of brainstem nuclei needed to test these theories has proven technologically challenging. Using methods to improve brainstem imaging in children, this study set out to elucidate the role of the autonomic, nociceptive, and limbic brainstem nuclei in the autism features of 145 children (74 autistic children, 6.0-10.9 years). Participants completed an assessment of core autism features and diffusion- and T1-weighted imaging optimized to improve brainstem images. After data reduction via principal component analysis, correlational analyses examined associations among autism features and the microstructural properties of brainstem clusters. Independent replication was performed in 43 adolescents (24 autistic, 13.0-17.9 years). We found specific nuclei, most robustly the parvicellular reticular formation-alpha (PCRtA) and to a lesser degree the lateral parabrachial nucleus (LPB) and ventral tegmental parabrachial pigmented complex (VTA-PBP), to be associated with autism features. The PCRtA and some of the LPB associations were independently found in the replication sample, but the VTA-PBP associations were not. Consistent with theoretical perspectives, the findings suggest that individual differences in pontine reticular formation nuclei contribute to the prominence of autistic features. Specifically, the PCRtA, a nucleus involved in mastication, digestion, and cardio-respiration in animal models, was associated with social communication in children, while the LPB, a pain-network nucleus, was associated with repetitive behaviors. These findings highlight the contributions of key autonomic brainstem nuclei to the expression of core autism features.

"Pscore": A Novel Percentile-Based Metric to Accurately Assess Individual Deviations in Non-Gaussian Distributions of Quantitative MRI Metrics

BACKGROUND

Quantitative magnetic resonance imaging (MRI) metrics could be used in personalized medicine to assess individuals against normative distributions. Conventional Zscore analysis is inadequate in the presence of non-Gaussian distributions. Therefore, if quantitative MRI metrics deviate from normality, an alternative is needed.

PURPOSE

To confirm non-Gaussianity of diffusion MRI (dMRI) metrics on a publicly available dataset, and to propose a novel percentile-based method, "Pscore" to address this issue.

STUDY TYPE

Retrospective cohort.

POPULATION

Nine hundred and sixty-one healthy young adults (age: 22-35 years, females: 53%) from the Human Connectome Project.

FIELD STRENGTH/SEQUENCE

3-T, spin-echo diffusion echo-planar imaging, T1-weighted: MPRAGE.

ASSESSMENT

The dMRI data were preprocessed using the TORTOISE pipeline. Forty-eight regions of interest (ROIs) from the JHU atlas were redrawn on a study-specific diffusion tensor (DT) template and average values were computed from various DT and mean apparent propagator (MAP) metrics. For each ROI, percentile ranks across participants were computed to generate "Pscores"-which normalized the difference between the median and a participant's value with the corresponding difference between the median and the 5th/95th percentile values.

STATISTICAL TESTS

ROI-wise distributions were assessed using log transformations, Zscore, and the "Pscore" methods. The percentages of extreme values above-95th and below-5th percentile boundaries (PEV>95 (%), PEV<5 (%)) were also assessed in the overall white matter. Bootstrapping was performed to test the reliability of Pscores in small samples (N = 100) using 100 iterations.

RESULTS

The dMRI metric distributions were systematically non-Gaussian, including positively skewed (eg, mean and radial diffusivity) and negatively skewed (eg, fractional and propagator anisotropy) metrics. This resulted in unbalanced tails in Zscore distributions (PEV>95  ≠ 5%, PEV<5  ≠ 5%) whereas "Pscore" distributions were symmetric and balanced (PEV>95  = PEV<5  = 5%); even for small bootstrapped samples (averagePEV > 95 ¯ = PEV < 5 ¯ = 5 ± 0 % $$ \overline{{\mathrm{PEV}}_{>95}}=\overline{{\mathrm{PEV}}_{<5}}=5\pm 0\% $$ [SD]).

DATA CONCLUSION

The inherent skewness observed for dMRI metrics may preclude the use of conventional Zscore analysis. The proposed "Pscore" method may help estimating individual deviations more accurately in skewed normative data, even from small datasets.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

Distinct alterations in white matter properties and organization related to maternal treatment initiation in neonates exposed to HIV but uninfected

HIV exposed-uninfected (HEU) infants and children are at risk of developmental delays as compared to uninfected unexposed (HUU) populations. The effects of exposure to in utero HIV and ART regimens on the HEU the developing brain are not well understood. In a cohort of 2-week-old newborns, we used diffusion tensor imaging (DTI) tractography and graph theory to examine the influence of HIV and ART exposure in utero on neonate white matter integrity and organisation. The cohort included HEU infants born to mothers who started ART before conception (HEUpre) and after conception (HEUpost), as well as HUU infants from the same community. We investigated HIV exposure and ART duration group differences in DTI metrics (fractional anisotropy (FA) and mean diffusivity (MD)) and graph measures across white matter. We found increased MD in white matter connections involving the thalamus and limbic system in the HEUpre group compared to HUU. We further identified reduced nodal efficiency in the basal ganglia. Within the HEUpost group, we observed reduced FA in cortical-subcortical and cerebellar connections as well as decreased transitivity in the hindbrain area compared to HUU. Overall, our analysis demonstrated distinct alterations in white matter integrity related to the timing of maternal ART initiation that influence regional brain network properties.

Neuronal microstructural changes in the human brain are associated with neurocognitive aging

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 unimpaired participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

'Pscore' - A Novel Percentile-Based Metric to Accurately Assess Individual Deviations in Non-Gaussian Distributions of Quantitative MRI Metrics

BACKGROUND

Quantitative MRI metrics could be used in personalized medicine to assess individuals against normative distributions. Conventional Zscore analysis is inadequate in the presence of non-Gaussian distributions. Therefore, if quantitative MRI metrics deviate from normality, an alternative is needed.

PURPOSE

To confirm non-Gaussianity of diffusion MRI (dMRI) metrics on a publicly available dataset, and to propose a novel percentile-based method, 'Pscore' to address this issue.

STUDY TYPE

Retrospective cohort.

POPULATION

961 healthy young-adults (age:22-35 years, Females:53%) from the Human Connectome Project.

FIELD STRENGTH/SEQUENCE

3-T, spin-echo diffusion echo-planar imaging, T1-weighted: MPRAGE.

ASSESSMENT

The dMRI data were preprocessed using the TORTOISE pipeline. Forty-eight regions of interest (ROIs) from the JHU-atlas were redrawn on a study-specific diffusion tensor (DT) template and average values were computed from various DT and mean apparent propagator (MAP) metrics. For each ROI, percentile ranks across participants were computed to generate 'Pscores'- which normalized the difference between the median and a participant's value with the corresponding difference between the median and the 5th/95th percentile values.

STATISTICAL TESTS

ROI-wise distributions were assessed using Log transformations, Zscore, and the 'Pscore' methods. The percentages of extreme values above-95th and below-5th percentile boundaries (PEV < 5 ( % ) ,PEV < 5 ( % ) ) were also assessed in the overall white matter. Bootstrapping was performed to test the reliability of Pscores in small samples (n=100) using 100 iterations.

RESULTS

The dMRI metric distributions were systematically non-Gaussian, including positively skewed (e.g., mean and radial distributions P E V > 95 ≠ 5 % , P E V < 5 ≠ 5 % whereas 'Pscore' distributions were symmetric and balanced P E V > 95 = P E V < 5 = 5 % ; even for small bootstrapped samples (average P E V > 95 ¯ = P E V < 5 ¯ = 5 ± 0 % S D ).

DATA CONCLUSION

The inherent skewness observed for dMRI metrics may preclude the use of conventional Zscore analysis. The proposed 'Pscore' method may help estimating individual deviations more accurately in skewed normative data, even from small datasets.

FD-Net: An unsupervised deep forward-distortion model for susceptibility artifact correction in EPI

PURPOSE

To introduce an unsupervised deep-learning method for fast and effective correction of susceptibility artifacts in reversed phase-encode (PE) image pairs acquired with echo planar imaging (EPI).

METHODS

Recent learning-based correction approaches in EPI estimate a displacement field, unwarp the reversed-PE image pair with the estimated field, and average the unwarped pair to yield a corrected image. Unsupervised learning in these unwarping-based methods is commonly attained via a similarity constraint between the unwarped images in reversed-PE directions, neglecting consistency to the acquired EPI images. This work introduces a novel unsupervised deep Forward-Distortion Network (FD-Net) that predicts both the susceptibility-induced displacement field and the underlying anatomically correct image. Unlike previous methods, FD-Net enforces the forward-distortions of the correct image in both PE directions to be consistent with the acquired reversed-PE image pair. FD-Net further leverages a multiresolution architecture to maintain high local and global performance.

RESULTS

FD-Net performs competitively with a gold-standard reference method (TOPUP) in image quality, while enabling a leap in computational efficiency. Furthermore, FD-Net outperforms recent unwarping-based methods for unsupervised correction in terms of both image and field quality.

CONCLUSION

The unsupervised FD-Net method introduces a deep forward-distortion approach to enable fast, high-fidelity correction of susceptibility artifacts in EPI by maintaining consistency to measured data. Therefore, it holds great promise for improving the anatomical accuracy of EPI imaging.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

Evaluating Reproducibility of the ADC and Distortion in Diffusion-weighted Imaging (DWI) with Reverse Encoding Distortion Correction (RDC)

PURPOSE

To compare image distortion and reproducibility of quantitative values between reverse encoding distortion correction (RDC) diffusion-weighted imaging (DWI) and conventional DWI techniques in a phantom study and in healthy volunteers.

METHODS

This prospective study was conducted with the approval of our institutional review board. Written informed consent was obtained from each participant. RDC-DWIs were created from images obtained at 3T in three orthogonal directions in a phantom and in 10 participants (mean age, 70.9 years; age range, 63-83 years). Images without distortion correction (noDC-DWI) and those corrected with B0 (B0c-DWI) were also created. To evaluate distortion, coefficients of variation were calculated for each voxel and ROIs were placed at four levels of the brain. To evaluate the reproducibility of apparent diffusion coefficient (ADC) measurements, intra- and inter-scan variability (%CVADC) were calculated from repeated scans of the phantom. Analysis was performed using Wilcoxon signed-rank test with Bonferroni correction, and P < 0.05 was considered statistically significant.

RESULTS

In the phantom, distortion was less in RDC-DWI than in B0c-DWI (P < 0.006), and was less in B0c-DWI than in noDC-DWI (P < 0.006). Intra-scan %CVADC was within 1.30%, and inter-scan %CVADC was within 2.99%. In the volunteers, distortion was less in RDC-DWI than in B0c-DWI in three of four locations (P < 0.006), and was less in B0c-DWI than in noDC-DWI (P < 0.006). At the middle cerebellar peduncle, distortion was less in RDC-DWI than in noDC-DWI (P < 0.006), and was less in noDC-DWI than in B0c-DWI (P < 0.0177).

CONCLUSION

In both the phantom and in volunteers, distortion was the least in RDC-DWI than in B0c-DWI and noDC-DWI.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Double pulsed field gradient diffusion MRI to assess skeletal muscle microstructure

PURPOSE

Preliminary study to determine whether double pulsed field gradient (PFG) diffusion MRI is sensitive to key features of muscle microstructure related to function.

METHODS

The restricted diffusion profile of molecules in models of muscle microstructure derived from histology were systematically simulated using a numerical simulation approach. Diffusion tensor subspace imaging analysis of the diffusion signal was performed, and spherical anisotropy (SA) was calculated for each model. Linear regression was used to determine the predictive capacity of SA on the fiber area, fiber diameter, and surface area to volume ratio of the models. Additionally, a rat model of muscle hypertrophy was scanned using a single PFG and a double PFG pulse sequence, and the restricted diffusion measurements were compared with histological measurements of microstructure.

RESULTS

Excellent agreement between SA and muscle fiber area (r2  = 0.71; p < 0.0001), fiber diameter (r2  = 0.83; p < 0.0001), and surface area to volume ratio (r2  = 0.97; p < 0.0001) in simulated models was found. In a scanned rat leg, the distribution of these microstructural features measured from histology was broad and demonstrated that there is a wide variance in the microstructural features observed, similar to the SA distributions. However, the distribution of fractional anisotropy measurements in the same tissue was narrow.

CONCLUSIONS

This study demonstrates that SA-a scalar value from diffusion tensor subspace imaging analysis-is highly sensitive to muscle microstructural features predictive of function. Furthermore, these techniques and analysis tools can be translated to real experiments in skeletal muscle. The increased dynamic range of SA compared with fractional anisotropy in the same tissue suggests increased sensitivity to detecting changes in tissue microstructure.

Sampling strategies and integrated reconstruction for reducing distortion and boundary slice aliasing in high-resolution 3D diffusion MRI

PURPOSE

To develop a new method for high-fidelity, high-resolution 3D multi-slab diffusion MRI with minimal distortion and boundary slice aliasing.

METHODS

Our method modifies 3D multi-slab imaging to integrate blip-reversed acquisitions for distortion correction and oversampling in the slice direction (kz ) for reducing boundary slice aliasing. Our aim is to achieve robust acceleration to keep the scan time the same as conventional 3D multi-slab acquisitions, in which data are acquired with a single direction of blip traversal and without kz -oversampling. We employ a two-stage reconstruction. In the first stage, the blip-up/down images are respectively reconstructed and analyzed to produce a field map for each diffusion direction. In the second stage, the blip-reversed data and the field map are incorporated into a joint reconstruction to produce images that are corrected for distortion and boundary slice aliasing.

RESULTS

We conducted experiments at 7T in six healthy subjects. Stage 1 reconstruction produces images from highly under-sampled data (R = 7.2) with sufficient quality to provide accurate field map estimation. Stage 2 joint reconstruction substantially reduces distortion artifacts with comparable quality to fully-sampled blip-reversed results (2.4× scan time). Whole-brain in-vivo results acquired at 1.22 mm and 1.05 mm isotropic resolutions demonstrate improved anatomical fidelity compared to conventional 3D multi-slab imaging. Data demonstrate good reliability and reproducibility of the proposed method over multiple subjects.

CONCLUSION

The proposed acquisition and reconstruction framework provide major reductions in distortion and boundary slice aliasing for 3D multi-slab diffusion MRI without increasing the scan time, which can potentially produce high-quality, high-resolution diffusion MRI.

Mechanical stiffness and anisotropy measured by MRE during brain development in the minipig

The relationship between brain development and mechanical properties of brain tissue is important, but remains incompletely understood, in part due to the challenges in measuring these properties longitudinally over time. In addition, white matter, which is composed of aligned, myelinated, axonal fibers, may be mechanically anisotropic. Here we use data from magnetic resonance elastography (MRE) and diffusion tensor imaging (DTI) to estimate anisotropic mechanical properties in six female Yucatan minipigs at ages from 3 to 6 months. Fiber direction was estimated from the principal axis of the diffusion tensor in each voxel. Harmonic shear waves in the brain were excited by three different configurations of a jaw actuator and measured using a motion-sensitive MR imaging sequence. Anisotropic mechanical properties are estimated from displacement field and fiber direction data with a finite element- based, transversely-isotropic nonlinear inversion (TI-NLI) algorithm. TI-NLI finds spatially resolved TI material properties that minimize the error between measured and simulated displacement fields. Maps of anisotropic mechanical properties in the minipig brain were generated for each animal at all four ages. These maps show that white matter is more dissipative and anisotropic than gray matter, and reveal significant effects of brain development on brain stiffness and structural anisotropy. Changes in brain mechanical properties may be a fundamental biophysical signature of brain development.

Properties of white matter tract diffusivity in children with developmental dyslexia and comorbid attention deficit/hyperactivity disorder

BACKGROUND

Developmental dyslexia (DD) and attention deficit/hyperactivity disorder (ADHD) are highly comorbid neurodevelopmental disorders. Individuals with DD or ADHD have both been shown to have deficits in white matter tracts associated with reading and attentional control networks. However, white matter diffusivity in individuals comorbid with both DD and ADHD (DD + ADHD) has not been specifically explored.

METHODS

Participants were 3rd and 4th graders (age range = 7 to 11 years; SD = 0.69) from three diagnostic groups ((DD (n = 40), DD + ADHD (n = 22), and typical developing (TD) (n = 20)). Behavioral measures of reading and attention alongside measures of white matter diffusivity were collected for all participants.

RESULTS

DD + ADHD and TD groups differed in mean fractional anisotropy (FA) for the left and right Superior Longitudinal Fasciculus (SLF)-Parietal Terminations and SLF-Temporal Terminations. Mean FA for the DD group across these SLF tracts fell between the lower DD + ADHD and higher TD averages. No differences in mean diffusivity nor significant brain-behavior relations were found.

CONCLUSIONS

Findings suggest that WM diffusivity in the SLF increases along a continuum across DD + ADHD, DD, and TD.

White matter microstructure and serum biomarkers of inflammation in psychogenic non-epileptic seizures

BACKGROUND

Neuroinflammation may contribute to the pathophysiology of psychogenic non-epileptic seizures (PNES). However, it is unclear whether and to what degree comorbid psychiatric symptoms explain this association. In this study, we investigated the neuroinflammatory signature of PNES and how it compares to that of people with psychiatric conditions (PwPCs).

METHODS

We prospectively assessed differences in neurite density (NDI), orientation dispersion (ODI), and isotropic diffusion (F-ISO) in 23 participants with PNES and 27 PwPCs, and their relationships to serum levels of tumor necrosis factor (TNF)-α, TNF receptor 1 (TNF-R1), TNF-related apoptosis-inducing ligand (TRAIL), interleukin (IL)-6, intercellular adhesion molecule (ICAM)-1, and monocyte chemoattractant protein (MCP)-1 using voxelwise multiple linear regressions. Pearson correlations between serum biomarkers and clinical symptoms were also obtained.

RESULTS

There were no white matter (WM) microstructural differences between groups. In PNES, TNF-R1 was negatively associated with NDI in the right uncinate fasciculus (UF) and positively associated with F-ISO in the left UF. IL-6 was positively associated with NDI and negatively with F-ISO in the left UF. ICAM-1 was positively associated with ODI in the left UF. TNF-α was negatively associated with ODI in the left cingulum bundle. The opposite relationships were observed in PwPCs. Higher TNF-R1 was associated with higher depression, anxiety, lower emotional quality of life, and higher levels of disability in PNES.

CONCLUSIONS

For the first time, we report relationships between peripheral inflammatory biomarkers and WM integrity in PNES, including abnormalities in the UF and cingulum bundle. Our results suggest that serum biomarkers of inflammation may, with additional studies, become a useful aid to PNES diagnosis, especially in settings where video-EEG is not available. The lack of group differences in WM microstructure suggests that previously identified WM abnormalities in PNES versus healthy controls may be related to psychological comorbidities of PNES.

Impact of asthma on the brain: evidence from diffusion MRI, CSF biomarkers and cognitive decline

Chronic systemic inflammation increases the risk of neurodegeneration, but the mechanisms remain unclear. Part of the challenge in reaching a nuanced understanding is the presence of multiple risk factors that interact to potentiate adverse consequences. To address modifiable risk factors and mitigate downstream effects, it is necessary, although difficult, to tease apart the contribution of an individual risk factor by accounting for concurrent factors such as advanced age, cardiovascular risk, and genetic predisposition. Using a case-control design, we investigated the influence of asthma, a highly prevalent chronic inflammatory disease of the airways, on brain health in participants recruited to the Wisconsin Alzheimer's Disease Research Center (31 asthma patients, 186 non-asthma controls, aged 45-90 years, 62.2% female, 92.2% cognitively unimpaired), a sample enriched for parental history of Alzheimer's disease. Asthma status was determined using detailed prescription information. We employed multi-shell diffusion weighted imaging scans and the three-compartment neurite orientation dispersion and density imaging model to assess white and gray matter microstructure. We used cerebrospinal fluid biomarkers to examine evidence of Alzheimer's disease pathology, glial activation, neuroinflammation and neurodegeneration. We evaluated cognitive changes over time using a preclinical Alzheimer cognitive composite. Using permutation analysis of linear models, we examined the moderating influence of asthma on relationships between diffusion imaging metrics, CSF biomarkers, and cognitive decline, controlling for age, sex, and cognitive status. We ran additional models controlling for cardiovascular risk and genetic risk of Alzheimer's disease, defined as a carrier of at least one apolipoprotein E (APOE) ε4 allele. Relative to controls, greater Alzheimer's disease pathology (lower amyloid-β42/amyloid-β40, higher phosphorylated-tau-181) and synaptic degeneration (neurogranin) biomarker concentrations were associated with more adverse white matter metrics (e.g. lower neurite density, higher mean diffusivity) in patients with asthma. Higher concentrations of the pleiotropic cytokine IL-6 and the glial marker S100B were associated with more salubrious white matter metrics in asthma, but not in controls. The adverse effects of age on white matter integrity were accelerated in asthma. Finally, we found evidence that in asthma, relative to controls, deterioration in white and gray matter microstructure was associated with accelerated cognitive decline. Taken together, our findings suggest that asthma accelerates white and gray matter microstructural changes associated with aging and increasing neuropathology, that in turn, are associated with more rapid cognitive decline. Effective asthma control, on the other hand, may be protective and slow progression of cognitive symptoms.

Age- and time-of-day dependence of glymphatic function in the human brain measured via two diffusion MRI methods

Advanced age, accompanied by impaired glymphatic function, is a key risk factor for many neurodegenerative diseases. To study age-related differences in the human glymphatic system, we measured the influx and efflux activities of the glymphatic system via two non-invasive diffusion magnetic resonance imaging (MRI) methods, ultra-long echo time and low-b diffusion tensor imaging (DTI_low-b) measuring the subarachnoid space (SAS) flow along the middle cerebral artery and DTI analysis along the perivascular space (DTI-ALPS) along medullary veins in 22 healthy volunteers (aged 21-75 years). We first evaluated the circadian rhythm dependence of the glymphatic activity by repeating the MRI measurements at five time points from 8:00 to 23:00 and found no time-of-day dependence in the awake state under the current sensitivity of MRI measurements. Further test-retest analysis demonstrated high repeatability of both diffusion MRI measurements, suggesting their reliability. Additionally, the influx rate of the glymphatic system was significantly higher in participants aged >45 years than in participants aged 21-38, while the efflux rate was significantly lower in those aged >45 years. The mismatched influx and efflux activities in the glymphatic system might be due to age-related changes in arterial pulsation and aquaporin-4 polarization.

Regional microstructural differences in ADHD youth with and without a family history of bipolar I disorder

BACKGROUND

Having a first-degree relative with bipolar I disorder (BD) in conjunction with prodromal attention deficit/hyperactivity disorder (ADHD) may represent a unique phenotype that confers greater risk for developing BD than ADHD alone. However, underlying neuropathoetiological mechanisms remain poorly understood. This cross-sectional study compared regional microstructure in psychostimulant-free ADHD youth with ('high-risk', HR) and without ('low-risk', LR) a first-degree relative with BD, and healthy controls (HC).

METHODS

A total of 140 (high-risk, n = 44; low-risk, n = 49; and HC, n = 47) youth (mean age: 14.1 ± 2.5 years, 65 % male) were included in the analysis. Diffusion tensor images were collected and fractional anisotropy (FA) and mean diffusivity (MD) maps calculated. Both tract-based and voxel-based analyses were performed. Correlations between clinical ratings and microstructural metrics that differed among groups were examined.

RESULTS

No significant group differences in major long-distance fiber tracts were observed. The high-risk ADHD group exhibited predominantly higher FA and lower MD in frontal, limbic, and striatal subregions compared with the low-risk ADHD group. Both low-risk and high-risk ADHD groups exhibited higher FA in unique and overlapping regions compared with HC subjects. Significant correlations between regional microstructural metrics and clinical ratings were observed in ADHD groups.

LIMITATIONS

Prospective longitudinal studies will be required to determine the relevance of these findings to BD risk progression.

CONCLUSIONS

Psychostimulant-free ADHD youth with a BD family history exhibit different microstructure alterations in frontal, limbic, and striatal regions compared with ADHD youth without a BD family history, and may therefore represent unique phenotypes relevant to BD risk progression.

Simultaneous Acquisition of Diffusion Tensor and Dynamic Diffusion MRI

BACKGROUND

Dynamic diffusion magnetic resonance imaging (ddMRI) metrics can assess transient microstructural alterations in tissue diffusivity but requires additional scan time hindering its clinical application.

PURPOSE

To determine whether a diffusion gradient table can simultaneously acquire data to estimate dynamic and diffusion tensor imaging (DTI) metrics.

STUDY TYPE

Prospective.

SUBJECTS

Seven healthy subjects, 39 epilepsy patients (15 female, 31 male, age ± 15).

FIELD STRENGTH/SEQUENCE

Two-dimensional diffusion MRI (b = 1000 s/mm2 ) at a field strength of 3 T. Sessions in healthy subjects-standard ddMRI (30 directions), standard DTI (15 and 30 directions), and nested cubes scans (15 and 30 directions). Sessions in epilepsy patients-two 30 direction (standard ddMRI, 10 nested cubes) or two 15 direction scans (standard DTI, 5 nested cubes).

ASSESSMENT

Fifteen direction DTI was repeated twice for within-session test-retest measurements in healthy subjects. Bland-Altman analysis computed bias and limits of agreement for DTI metrics using test-retest scans and standard 15 direction vs. 5 nested cubes scans. Intraclass correlation (ICC) analysis compared tensor metrics between 15 direction DTI scans (standard vs. 5 nested cubes) and the coefficients of variation (CoV) of trace and apparent diffusion coefficient (ADC) between 30 direction ddMRI scans (standard vs. 10 nested cubes).

STATISTICAL TESTS

Bland-Altman and ICC analysis using a P-value of 0.05 for statistical significance.

RESULTS

Correlations of mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were strong and significant in gray (ICC > 0.95) and white matter (ICC > 0.95) between standard vs. nested cubes DTI acquisitions. Correlation of white matter fractional anisotropy was also strong (ICC > 0.95) and significant. ICCs of the CoV of dynamic ADC measured using repeated cubes and nested cubes acquisitions were modest (ICC >0.60), but significant in gray matter.

CONCLUSION

A nested cubes diffusion gradient table produces tensor-based and dynamic diffusion measurements in a single acquisition.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

The entorhinal-DG/CA3 pathway in the medial temporal lobe retains visual working memory of a simple surface feature

Classic models consider working memory (WM) and long-term memory as distinct mental faculties that are supported by different neural mechanisms. Yet, there are significant parallels in the computation that both types of memory require. For instance, the representation of precise item-specific memory requires the separation of overlapping neural representations of similar information. This computation has been referred to as pattern separation, which can be mediated by the entorhinal-DG/CA3 pathway of the medial temporal lobe (MTL) in service of long-term episodic memory. However, although recent evidence has suggested that the MTL is involved in WM, the extent to which the entorhinal-DG/CA3 pathway supports precise item-specific WM has remained elusive. Here, we combine an established orientation WM task with high-resolution fMRI to test the hypothesis that the entorhinal-DG/CA3 pathway retains visual WM of a simple surface feature. Participants were retrospectively cued to retain one of the two studied orientation gratings during a brief delay period and then tried to reproduce the cued orientation as precisely as possible. By modeling the delay-period activity to reconstruct the retained WM content, we found that the anterior-lateral entorhinal cortex (aLEC) and the hippocampal DG/CA3 subfield both contain item-specific WM information that is associated with subsequent recall fidelity. Together, these results highlight the contribution of MTL circuitry to item-specific WM representation.

DACO: Distortion/artefact correction for diffusion MRI data

In this paper, we propose a registration-based algorithm to correct various distortions or artefacts (DACO) commonly observed in diffusion-weighted (DW) magnetic resonance images (MRI). The registration in DACO is accomplished by means of a pseudo b₀ image, which is synthesized from the anatomical images such as T1-weighted image or T2-weighted image, and a pseudo diffusion MRI (dMRI) data, which is derived from the Gaussian model of diffusion tensor imaging (DTI) or the Hermite model of mean apparent propagator (MAP)-MRI. DACO corrects (1) the susceptibility-induced distortions and (2) the misalignment between the dMRI data and anatomical images by registering the real b₀ image to the pseudo b₀ image, and corrects (3) the eddy current-induced distortions and (4) the head motions by registering each image in the real dMRI data to the corresponding image in the pseudo dMRI data. DACO estimates the models of artefacts simultaneously in an iterative and interleaved manner. The mathematical formulation of the models and the estimation procedures are detailed in this paper. Using the human connectome project (HCP) data the evaluation shows that DACO could estimate the model parameters accurately. Furthermore, the evaluation conducted on the real human data acquired from clinical MRI scanners reveals that the method could reduce the artefacts effectively. The DACO method leverages the anatomical image, which is routinely acquired in clinical practice, to correct the artefacts, omitting the additional acquisitions needed to conduct the algorithm. Therefore, our method should be beneficial to most dMRI data, particularly to those acquired without field maps or reverse phase-encoding images.

Assessment of intraoperative diffusion EPI distortion and its impact on estimation of supratentorial white matter tract positions in pediatric epilepsy surgery

The effectiveness of correcting diffusion Echo Planar Imaging (EPI) distortion and its impact on tractography reconstruction have not been adequately investigated in the intraoperative MRI setting, particularly for High Angular Resolution Diffusion Imaging (HARDI) acquisition. In this study, we evaluated the effectiveness of EPI distortion correction using 27 legacy intraoperative HARDI datasets over two consecutive surgical time points, acquired without reverse phase-encoded data, from 17 children who underwent epilepsy surgery at our institution. The data was processed with EPI distortion correction using the Synb0-Disco technique (Schilling et al., 2019) and without distortion correction. The corrected and uncorrected b0 diffusion-weighted images (DWI) were first compared visually. The mutual information indices between the original T1-weighted images and the fractional anisotropy images derived from corrected and uncorrected DWI were used to quantify the effect of distortion correction. Sixty-four white matter tracts were segmented from each dataset, using a deep-learning based automated tractography algorithm for the purpose of a standardized and unbiased evaluation. Displacement was calculated between tracts generated before and after distortion correction. The tracts were grouped based on their principal morphological orientations to investigate whether the effects of EPI distortion vary with tract orientation. Group differences in tract distortion were investigated both globally, and regionally with respect to proximity to the resecting lesion in the operative hemisphere. Qualitatively, we observed notable improvement in the corrected diffusion images, over the typically affected brain regions near skull-base air sinuses, and correction of additional distortion unique to intraoperative open cranium images, particularly over the resection site. This improvement was supported quantitatively, as mutual information indices between the FA and T1-weighted images were significantly greater after the correction, compared to before the correction. Maximum tract displacement between the corrected and uncorrected data, was in the range of 7.5 to 10.0 mm, a magnitude that would challenge the safety resection margin typically tolerated for tractography-informed surgical guidance. This was particularly relevant for tracts oriented partially or fully in-line with the acquired diffusion phase-encoded direction. Portions of these tracts passing close to the resection site demonstrated significantly greater magnitude of displacement, compared to portions of tracts remote from the resection site in the operative hemisphere. Our findings have direct clinical implication on the accuracy of intraoperative tractography-informed image guidance and emphasize the need to develop a distortion correction technique with feasible intraoperative processing time.

Performance Comparison of Different Neuroimaging Methods for Predicting Upper Limb Motor Outcomes in Patients after Stroke

Several neuroimaging methods have been proposed to assess the integrity of the corticospinal tract (CST) for predicting recovery of motor function after stroke, including conventional structural magnetic resonance imaging (sMRI) and diffusion tensor imaging (DTI). In this study, we aimed to compare the predicative performance of these methods using different neuroimaging modalities and optimize the prediction protocol for upper limb motor function after stroke in a clinical environment. We assessed 28 first-ever stroke patients with upper limb motor impairment. We used the upper extremity module of the Fugl-Meyer assessment (UE-FM) within 1 month of onset (baseline) and again 3 months poststroke. sMRI (T1- and T2-based) was used to measure CST-weighted lesion load (CST-wLL), and DTI was used to measure the fractional anisotropy asymmetry index (FAAI) and the ratio of fractional anisotropy (rFA). The CST-wLL within 1 month poststroke was closely correlated with upper limb motor outcomes and recovery potential. CST‐wLL ≥ 2.068 cc indicated serious CST damage and a poor outcome (100%). CST‐wLL < 1.799 cc was correlated with a considerable rate (>70%) of upper limb motor function recovery. CST-wLL showed a comparable area under the curve (AUC) to that of the CST-FAAI (p = 0.71). Inclusion of extra-CST-FAAI did not significantly increase the AUC (p = 0.58). Our findings suggest that sMRI-derived CST-wLL is a precise predictor of upper limb motor outcomes 3 months poststroke. We recommend this parameter as a predictive imaging biomarker for classifying patients' recovery prognosis in clinical practice. Conversely, including DTI appeared to induce no significant benefits.

Combining MRI and cognitive evaluation to classify concussion in university athletes

Current methods of concussion assessment lack the objectivity and reliability to detect neurological injury. This multi-site study uses combinations of neuroimaging (diffusion tensor imaging and resting state functional MRI) and cognitive measures to train algorithms to detect the presence of concussion in university athletes. Athletes (29 concussed, 48 controls) completed symptom reports, brief cognitive evaluation, and MRI within 72 h of injury. Hierarchical linear regression compared groups on cognitive and neuroimaging measures while controlling for sex and data collection site. Logistic regression and support vector machine models were trained using cognitive and neuroimaging measures and evaluated for overall accuracy, sensitivity, and specificity. Concussed athletes reported greater symptoms than controls (∆R² = 0.32, p < .001), and performed worse on tests of concentration (∆R² = 0.07, p < .05) and delayed memory (∆R² = 0.17, p < .001). Concussed athletes showed lower functional connectivity within the frontoparietal and primary visual networks (p < .05), but did not differ on mean diffusivity and fractional anisotropy. Of the cognitive measures, classifiers trained using delayed memory yielded the best performance with overall accuracy of 71%, though sensitivity was poor at 46%. Of the neuroimaging measures, classifiers trained using mean diffusivity yielded similar accuracy. Combining cognitive measures with mean diffusivity increased overall accuracy to 74% and sensitivity to 64%, comparable to the sensitivity of symptom report. Trained algorithms incorporating both MRI and cognitive performance variables can reliably detect common neurobiological sequelae of acute concussion. The integration of multi-modal data can serve as an objective, reliable tool in the assessment and diagnosis of concussion.

Unsupervised Deep Learning for FOD-Based Susceptibility Distortion Correction in Diffusion MRI

Susceptibility induced distortion is a major artifact that affects the diffusion MRI (dMRI) data analysis. In the Human Connectome Project (HCP), the state-of-the-art method adopted to correct this kind of distortion is to exploit the displacement field from the B0 image in the reversed phase encoding images. However, both the traditional and learning-based approaches have limitations in achieving high correction accuracy in certain brain regions, such as brainstem. By utilizing the fiber orientation distribution (FOD) computed from the dMRI, we propose a novel deep learning framework named DistoRtion Correction Net (DrC-Net), which consists of the U-Net to capture the latent information from the 4D FOD images and the spatial transformer network to propagate the displacement field and back propagate the losses between the deformed FOD images. The experiments are performed on two datasets acquired with different phase encoding (PE) directions including the HCP and the Human Connectome Low Vision (HCLV) dataset. Compared to two traditional methods topup and FODReg and two deep learning methods S-Net and flow-net, the proposed method achieves significant improvements in terms of the mean squared difference (MSD) of fractional anisotropy (FA) images and minimum angular difference between two PEs in white matter and also brainstem regions. In the meantime, the proposed DrC-Net takes only several seconds to predict a displacement field, which is much faster than the FODReg method.

Critical factors in achieving fine-scale functional MRI: Removing sources of inadvertent spatial smoothing

Ultra‐high Field (≥7T) functional magnetic resonance imaging (UHF‐fMRI) provides opportunities to resolve fine‐scale features of functional architecture such as cerebral cortical columns and layers, in vivo. While the nominal resolution of modern fMRI acquisitions may appear to be sufficient to resolve these features, several common data preprocessing steps can introduce unwanted spatial blurring, especially those that require interpolation of the data. These resolution losses can impede the detection of the fine‐scale features of interest. To examine quantitatively and systematically the sources of spatial resolution losses occurring during preprocessing, we used synthetic fMRI data and real fMRI data from the human visual cortex—the spatially interdigitated human V2 “thin” and “thick” stripes. The pattern of these cortical columns lies along the cortical surface and thus can be best appreciated using surface‐based fMRI analysis. We used this as a testbed for evaluating strategies that can reduce spatial blurring of fMRI data. Our results show that resolution losses can be mitigated at multiple points in preprocessing pathway. We show that unwanted blur is introduced at each step of volume transformation and surface projection, and can be ameliorated by replacing multi‐step transformations with equivalent single‐step transformations. Surprisingly, the simple approaches of volume upsampling and of cortical mesh refinement also helped to reduce resolution losses caused by interpolation. Volume upsampling also serves to improve motion estimation accuracy, which helps to reduce blur. Moreover, we demonstrate that the level of spatial blurring is nonuniform over the brain—knowledge which is critical for interpreting data in high‐resolution fMRI studies. Importantly, our study provides recommendations for reducing unwanted blurring during preprocessing as well as methods that enable quantitative comparisons between preprocessing strategies. These findings highlight several underappreciated sources of a spatial blur. Individually, the factors that contribute to spatial blur may appear to be minor, but in combination, the cumulative effects can hinder the interpretation of fine‐scale fMRI and the detectability of these fine‐scale features of functional architecture.

Improving Imaging of the Brainstem and Cerebellum in Autistic Children: Transformation-Based High-Resolution Diffusion MRI (TiDi-Fused) in the Human Brainstem

Diffusion-weighted magnetic resonance imaging (dMRI) of the brainstem is technically challenging, especially in young autistic children as nearby tissue-air interfaces and motion (voluntary and physiological) can lead to artifacts. This limits the availability of high-resolution images, which are desirable for improving the ability to study brainstem structures. Furthermore, inherently low signal-to-noise ratios, geometric distortions, and sensitivity to motion not related to molecular diffusion have resulted in limited techniques for high-resolution data acquisition compared to other modalities such as T1-weighted imaging. Here, we implement a method for achieving increased apparent spatial resolution in pediatric dMRI that hinges on accurate geometric distortion correction and on high fidelity within subject image registration between dMRI and magnetization prepared rapid acquisition gradient echo (MPnRAGE) images. We call this post-processing pipeline T1 weighted-diffusion fused, or "TiDi-Fused". Data used in this work consists of dMRI data (2.4 mm resolution, corrected using FSL's Topup) and T1-weighted (T1w) MPnRAGE anatomical data (1 mm resolution) acquired from 128 autistic and non-autistic children (ages 6-10 years old). Accurate correction of geometric distortion permitted for a further increase in apparent resolution of the dMRI scan via boundary-based registration to the MPnRAGE T1w. Estimation of fiber orientation distributions and further analyses were carried out in the T1w space. Data processed with the TiDi-Fused method were qualitatively and quantitatively compared to data processed with conventional dMRI processing methods. Results show the advantages of the TiDi-Fused pipeline including sharper brainstem gray-white matter tissue contrast, improved inter-subject spatial alignment for group analyses of dMRI based measures, accurate spatial alignment with histology-based imaging of the brainstem, reduced variability in brainstem-cerebellar white matter tracts, and more robust biologically plausible relationships between age and brainstem-cerebellar white matter tracts. Overall, this work identifies a promising pipeline for achieving high-resolution imaging of brainstem structures in pediatric and clinical populations who may not be able to endure long scan times. This pipeline may serve as a gateway for feasibly elucidating brainstem contributions to autism and other conditions.

Short-Range Structural Connections Are More Severely Damaged in Early-Stage MS

BACKGROUND AND PURPOSE

Long-range connections are more severely damaged and relevant for cognition in long-standing MS. However, the evolution of such coordinated network damage in patients with MS is unclear. We investigated whether short- and long-range structural connections sustained equal damage in early-stage MS.

MATERIALS AND METHODS

Sixteen patients with early-stage MS and 17 healthy controls were scanned by high-resolution, multishell diffusion imaging on 7T MR imaging and assessed cognitively. We investigated macrostructural properties in short- and long-range fibers and of microstructural metrics derived from 2 quantitative diffusion MR imaging models: DTI and neurite orientation dispersion and density imaging.

RESULTS

Patients had significant WM integrity damage-that is, higher radial diffusivity and a lower intracellular volume fraction in the focal WM lesions. Compared with the healthy controls, the patients had noticeable microstructure changes in both short- and long-range fibers, including increased radial diffusivity, mean diffusivity, and axial diffusivity. Z scores further indicated greater damage in the short-range fibers than in the long-range fibers.

CONCLUSIONS

Our findings demonstrate that more severe demyelination preceding axonal degeneration occurs in short-range connections but not in long-range connections in early-stage MS, suggesting the possibility that there are cortical lesions that are undetectable by current MR imaging.

Translationally Relevant Magnetic Resonance Imaging Markers in a Ferret Model of Closed Head Injury

Pre-clinical models of traumatic brain injury (TBI) have been the primary experimental tool for understanding the potential mechanisms and cellular alterations that follow brain injury, but the human relevance and translational value of these models are often called into question. Efforts to better recapitulate injury biomechanics and the use of non-rodent species with neuroanatomical similarities to humans may address these concerns and promise to advance experimental studies toward clinical impact. In addition to improving translational aspects of animal models, it is also advantageous to establish pre-clinical outcomes that can be directly compared with the same outcomes in humans. Non-invasive imaging and particularly MRI is promising for this purpose given that MRI is a primary tool for clinical diagnosis and at the same time increasingly available at the pre-clinical level. The objective of this study was to identify which commonly used radiologic markers of TBI outcomes can be found also in a translationally relevant pre-clinical model of TBI. The ferret was selected as a human relevant species for this study with folded cortical geometry and relatively high white matter content and the closed head injury model of engineered rotation and acceleration (CHIMERA) TBI model was selected for biomechanical similarities to human injury. A comprehensive battery of MRI protocols based on common data elements (CDEs) for human TBI was collected longitudinally for the identification of MRI markers and voxelwise analysis of T2, contrast enhancement and diffusion tensor MRI values. The most prominent MRI findings were consistent with focal hemorrhage and edema in the brain stem region following high severity injury as well as vascular and meningeal injury evident by contrast enhancement. While conventional MRI outcomes were not highly conspicuous in less severe cases, quantitative voxelwise analysis indicated diffusivity and anisotropy alterations in the acute and chronic periods after TBI. The main conclusions of this study support the translational relevance of closed head TBI models in intermediate species and identify brain stem and meningeal vulnerability. Additionally, the MRI findings highlight a subset of CDEs with promise to bridge pre-clinical studies with human TBI outcomes.

Prenatal depression exposure alters white matter integrity and neurodevelopment in early childhood

Prenatal exposure to maternal depression increases the risk for onset of emotional and behavioral disorders in children. We investigated the effects of exposure to prenatal depression on white matter microstructural integrity at birth and at 2-3 years, and associated neurodevelopment. Diffusion-weighted images were acquired for children of the Drakenstein Child Health Study at 2-4 weeks postpartum (n=70, 47% boys) and at 2-3 years of age (n=60, 58% boys). Tract-Based Spatial Statistics was used to compare, using an ROI based approach, diffusion tensor metrics across groups defined by presence (>19 on Beck's Depression Inventory and/or >12 on the Edinburgh Postnatal Depression Scale) or absence (below depression thresholds) of depression, and associations with neurodevelopmental measures at age 2-3 years were determined. We did not detect group differences in white matter integrity at neonatal age, but at 2-3 years, children in the exposed group demonstrated higher fractional anisotropy, and lower mean and radial diffusivity in association tracts compared to controls. This was notable in the sagittal stratum (radial diffusivity: p<0.01). Altered white matter integrity metrics were also observed in projection tracts, including the corona radiata, which associated with cognitive and motor outcomes in exposed 2-3-year-olds (p<0.05). Our findings of widespread white matter alterations in 2-3-year-old children with prenatal exposure to depression are consistent with previous findings, as well as with neuroimaging findings in adults with major depression. Further, we identified novel associations of altered white matter integrity with cognitive development in depression-exposed children, suggesting that these neuroimaging findings may have early functional impact.

Manifold-based respiratory phase estimation enables motion and distortion correction of free-breathing cardiac diffusion tensor MRI

PURPOSE

For in vivo cardiac DTI, breathing motion and B₀ field inhomogeneities produce misalignment and geometric distortion in diffusion-weighted (DW) images acquired with conventional single-shot EPI. We propose using a dimensionality reduction method to retrospectively estimate the respiratory phase of DW images and facilitate both distortion correction (DisCo) and motion compensation.

METHODS

Free-breathing electrocardiogram-triggered whole left-ventricular cardiac DTI using a second-order motion-compensated spin echo EPI sequence and alternating directionality of phase encoding blips was performed on 11 healthy volunteers. The respiratory phase of each DW image was estimated after projecting the DW images into a 2D space with Laplacian eigenmaps. DisCo and motion compensation were applied to the respiratory sorted DW images. The results were compared against conventional breath-held T₂ half-Fourier single shot turbo spin echo. Cardiac DTI parameters including fractional anisotropy, mean diffusivity, and helix angle transmurality were compared with and without DisCo.

RESULTS

The left-ventricular geometries after DisCo and motion compensation resulted in significantly improved alignment of DW images with T₂ reference. DisCo reduced the distance between the left-ventricular contours by 13.2% ± 19.2%, P < .05 (2.0 ± 0.4 for DisCo and 2.4 ± 0.5 mm for uncorrected). DisCo DTI parameter maps yielded no significant differences (mean diffusivity: 1.55 ± 0.13 × 10^-3 mm² /s and 1.53 ± 0.13 × 10^-3 mm² /s, P = .09; fractional anisotropy: 0.375 ± 0.041 and 0.379 ± 0.045, P = .11; helix angle transmurality: 1.00% ± 0.10°/% and 0.99% ± 0.12°/%, P = .44), although the orientation of individual tensors differed.

CONCLUSION

Retrospective respiratory phase estimation with LE-based DisCo and motion compensation in free-breathing cardiac DTI resulting in significantly reduced geometric distortion and improved alignment within and across slices.

What's New and What's Next in Diffusion MRI Preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B₁ bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing.

MR Imaging of Human Brain Mechanics In Vivo: New Measurements to Facilitate the Development of Computational Models of Brain Injury

Computational models of the brain and its biomechanical response to skull accelerations are important tools for understanding and predicting traumatic brain injuries (TBIs). However, most models have been developed using experimental data collected on animal models and cadaveric specimens, both of which differ from the living human brain. Here we describe efforts to noninvasively measure the biomechanical response of the human brain with MRI-at non-injurious strain levels-and generate data that can be used to develop, calibrate, and evaluate computational brain biomechanics models. Specifically, this paper reports on a project supported by the National Institute of Neurological Disorders and Stroke to comprehensively image brain anatomy and geometry, mechanical properties, and brain deformations that arise from impulsive and harmonic skull loadings. The outcome of this work will be a publicly available dataset ( http://www.nitrc.org/projects/bbir ) that includes measurements on both males and females across an age range from adolescence to older adulthood. This article describes the rationale and approach for this study, the data available, and how these data may be used to develop new computational models and augment existing approaches; it will serve as a reference to researchers interested in using these data.

Disruption of the white matter structural network and its correlation with baseline progression rate in patients with sporadic amyotrophic lateral sclerosis

OBJECTIVE

There is increasing evidence that amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease impacting large-scale brain networks. However, it is still unclear which structural networks are associated with the disease and whether the network connectomics are associated with disease progression. This study was aimed to characterize the network abnormalities in ALS and to identify the network-based biomarkers that predict the ALS baseline progression rate.

METHODS

Magnetic resonance imaging was performed on 73 patients with sporadic ALS and 100 healthy participants to acquire diffusion-weighted magnetic resonance images and construct white matter (WM) networks using tractography methods. The global and regional network properties were compared between ALS and healthy subjects. The single-subject WM network matrices of patients were used to predict the ALS baseline progression rate using machine learning algorithms.

RESULTS

Compared with the healthy participants, the patients with ALS showed significantly decreased clustering coefficient Cp (P = 0.0034, t = 2.98), normalized clustering coefficient γ (P = 0.039, t = 2.08), and small-worldness σ (P = 0.038, t = 2.10) at the global network level. The patients also showed decreased regional centralities in motor and non-motor systems including the frontal, temporal and subcortical regions. Using the single-subject structural connection matrix, our classification model could distinguish patients with fast versus slow progression rate with an average accuracy of 85%.

CONCLUSION

Disruption of the WM structural networks in ALS is indicated by weaker small-worldness and disturbances in regions outside of the motor systems, extending the classical pathophysiological understanding of ALS as a motor disorder. The individual WM structural network matrices of ALS patients are potential neuroimaging biomarkers for the baseline disease progression in clinical practice.

White matter integrity, duration of untreated psychosis, and antipsychotic treatment response in medication-naïve first-episode psychosis patients

It is becoming increasingly clear that longer duration of untreated psychosis (DUP) is associated with adverse clinical outcomes in patients with psychosis spectrum disorders. Because this association is often cited when justifying early intervention efforts, it is imperative to better understand underlying biological mechanisms. We enrolled 66 antipsychotic-naïve first-episode psychosis (FEP) patients and 45 matched healthy controls in this trial. At baseline, we used a human connectome style diffusion-weighted imaging (DWI) sequence to quantify white matter integrity in both groups. Patients then received 16 weeks of treatment with risperidone, 51 FEP completed the trial. We compared whole-brain fractional anisotropy (FA), mean diffusivity, axial diffusivity (AD), and radial diffusivity between groups. To test if structural white matter integrity mediates the relationship between longer DUP and poorer treatment response, we fit a mediator model and estimated indirect effects. We found decreased whole-brain FA and AD in medication-naive FEP compared with controls. In patients, lower FA was correlated with longer DUP (r = -0.32; p = 0.03) and poorer subsequent response to antipsychotic treatment (r = 0.40; p = 0.01). Importantly, we found a significant mediation effect for FA (indirect effect: -2.70; p = 0.03), indicating that DUP exerts its effects on treatment response through affecting white matter integrity. Our data provide empirical support to the idea the DUP may have fundamental pathogenic effects on the natural history of psychosis, suggest a biological mechanism underlying this phenomenon, and underscore the importance of early intervention efforts in this disabling neuropsychiatric syndrome.

Reducing Susceptibility Distortion Related Image Blurring in Diffusion MRI EPI Data

Diffusion magnetic resonance imaging (MRI) is an increasingly popular technique in basic and clinical neuroscience. One promising application is to combine diffusion MRI with myelin maps from complementary MRI techniques such as multi-parameter mapping (MPM) to produce g-ratio maps that represent the relative myelination of axons and predict their conduction velocity. Statistical Parametric Mapping (SPM) can process both diffusion data and MPMs, making SPM the only widely accessible software that contains all the processing steps required to perform group analyses of g-ratio data in a common space. However, limitations have been identified in its method for reducing susceptibility-related distortion in diffusion data. More generally, susceptibility-related image distortion is often corrected by combining reverse phase-encoded images (blip-up and blip-down) using the arithmetic mean (AM), however, this can lead to blurred images. In this study we sought to (1) improve the susceptibility-related distortion correction for diffusion MRI data in SPM; (2) deploy an alternative approach to the AM to reduce image blurring in diffusion MRI data when combining blip-up and blip-down EPI data after susceptibility-related distortion correction; and (3) assess the benefits of these changes for g-ratio mapping. We found that the new processing pipeline, called consecutive Hyperelastic Susceptibility Artefact Correction (HySCO) improved distortion correction when compared to the standard approach in the ACID toolbox for SPM. Moreover, using a weighted average (WA) method to combine the distortion corrected data from each phase-encoding polarity achieved greater overlap of diffusion and more anatomically faithful structural white matter probability maps derived from minimally distorted multi-parameter maps as compared to the AM. Third, we showed that the consecutive HySCO WA performed better than the AM method when combined with multi-parameter maps to perform g-ratio mapping. These improvements mean that researchers can conveniently access a wide range of diffusion-related analysis methods within one framework because they are now available within the open-source ACID toolbox as part of SPM, which can be easily combined with other SPM toolboxes, such as the hMRI toolbox, to facilitate computation of myelin biomarkers that are necessary for g-ratio mapping.

Central white matter integrity alterations in 2-3-year-old children following prenatal alcohol exposure

BACKGROUND

Prenatal alcohol exposure (PAE) remains a potentially preventable, but pervasive risk factor to neurodevelopment. Yet, evidence is lacking on the impact of alcohol on brain development in toddlers. This study aimed to investigate the impact of PAE on brain white matter integrity in 2-3-year-old children.

METHODS

Children (n = 83, 30-37 months old) of the Drakenstein Child Health Study birth cohort, underwent diffusion MRI on a 3 T Siemens scanner during natural sleep. Parameters were extracted in children with PAE (n = 25, 56 % boys) and unexposed controls (n = 58, 62 % boys) using Tract-based Spatial Statistics, and compared by group. The contribution of maternal tobacco smoking to white matter differences was also explored.

RESULTS

Children with PAE had altered fractional anisotropy, radial diffusivity and axial diffusivity in brain stem, limbic and association tracts compared to unexposed controls. Notably lower fractional anisotropy was found in the uncinate fasciculus, and lower mean and radial diffusivity were found in the fornix stria terminalis and corticospinal tract (FDR corrected p < 0.05). There was a significant interaction effect of PAE and prenatal tobacco exposure which lowered mean, radial and axial diffusivity in the corticospinal tract significantly in the PAE group but not controls.

CONCLUSION

Widespread altered white matter microstructural integrity at 2-3 years of age is consistent with findings in neonates in the same and other cohorts, indicating persistence of effects of PAE through early life. Findings also highlight that prenatal tobacco exposure impacts the association of PAE on white matter alterations, amplifying effects in tracts underlying motor function.

PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images

PURPOSE

Diffusion weighted MRI imaging (DWI) is often subject to low signal-to-noise ratios (SNRs) and artifacts. Recent work has produced software tools that can correct individual problems, but these tools have not been combined with each other and with quality assurance (QA). A single integrated pipeline is proposed to perform DWI preprocessing with a spectrum of tools and produce an intuitive QA document.

METHODS

The proposed pipeline, built around the FSL, MRTrix3, and ANTs software packages, performs DWI denoising; inter-scan intensity normalization; susceptibility-, eddy current-, and motion-induced artifact correction; and slice-wise signal drop-out imputation. To perform QA on the raw and preprocessed data and each preprocessing operation, the pipeline documents qualitative visualizations, quantitative plots, gradient verifications, and tensor goodness-of-fit and fractional anisotropy analyses.

RESULTS

Raw DWI data were preprocessed and quality checked with the proposed pipeline and demonstrated improved SNRs; physiologic intensity ratios; corrected susceptibility-, eddy current-, and motion-induced artifacts; imputed signal-lost slices; and improved tensor fits. The pipeline identified incorrect gradient configurations and file-type conversion errors and was shown to be effective on externally available datasets.

CONCLUSIONS

The proposed pipeline is a single integrated pipeline that combines established diffusion preprocessing tools from major MRI-focused software packages with intuitive QA.

Multiple parietal pathways are associated with rTMS-induced hippocampal network enhancement and episodic memory changes

Repetitive transcranial magnetic stimulation (rTMS) of the inferior parietal cortex (IPC) increases resting-state functional connectivity (rsFC) of the hippocampus with the precuneus and other posterior cortical areas and causes proportional improvement of episodic memory. The anatomical pathway(s) responsible for the propagation of these effects from the IPC is unknown and may not be direct. In order to assess the relative contributions of candidate pathways from the IPC to the MTL via the parahippocampal cortex and precuneus, to the effects of rTMS on rsFC and memory improvement, we used diffusion tensor imaging to measure the extent to which individual differences in fractional anisotropy (FA) in these pathways accounted for individual differences in response. FA in the IPC-parahippocampal pathway and several MTL pathways predicted changes in rsFC. FA in both parahippocampal and hippocampal pathways was related to changes in episodic, but not procedural, memory. These results implicate pathways to the MTL in the enhancing effect of parietal rTMS on hippocampal rsFC and memory.

Improved reproducibility of diffusion MRI of the human brain with a four-way blip-up and down phase-encoding acquisition approach

PURPOSE

To assess the effects of blip-up and -down echo planar imaging (EPI) acquisition designs, with different choices of phase-encoding directions (PEDs) on the reproducibility of diffusion MRI (dMRI)-derived metrics in the human brain.

METHODS

Diffusion MRI data in seven subjects were acquired five times, each with five different protocols. The base design included 64 diffusion directions acquired with anterior-posterior (AP) PED, the first and second protocols added reverse phase-encoded b = 0 s / mm 2 posterior-anterior (PA) PED images. The third one included 32 directions all with PED acquisitions with opposite polarity (AP and PA). The fourth protocol, also with 32 unique directions used four PEDs (AP, PA, right-left (RL), and left-right (LR)). The scan time was virtually identical for all protocols. The variability of diffusion MRI metrics for each subject and each protocol was computed across the different sessions.

RESULTS

The highest reproducibility for all dMRI metrics was obtained with protocol four (AP/PA-RL/LR, ie, four-way PED). Protocols that used only b = 0 s / mm 2 for distortion correction, which are the most widely used designs, had the lowest reproducibility.

CONCLUSIONS

An acquisition design with four PEDs, including all DWIs in addition to b = 0 s / mm 2 images should be used to achieve high reproducibility in diffusion MRI studies.

Correcting Susceptibility Artifacts of MRI Sensors in Brain Scanning: A 3D Anatomy-Guided Deep Learning Approach

Echo planar imaging (EPI), a fast magnetic resonance imaging technique, is a powerful tool in functional neuroimaging studies. However, susceptibility artifacts, which cause misinterpretations of brain functions, are unavoidable distortions in EPI. This paper proposes an end-to-end deep learning framework, named TS-Net, for susceptibility artifact correction (SAC) in a pair of 3D EPI images with reversed phase-encoding directions. The proposed TS-Net comprises a deep convolutional network to predict a displacement field in three dimensions to overcome the limitation of existing methods, which only estimate the displacement field along the dominant-distortion direction. In the training phase, anatomical T1-weighted images are leveraged to regularize the correction, but they are not required during the inference phase to make TS-Net more flexible for general use. The experimental results show that TS-Net achieves favorable accuracy and speed trade-off when compared with the state-of-the-art SAC methods, i.e., TOPUP, TISAC, and S-Net. The fast inference speed (less than a second) of TS-Net makes real-time SAC during EPI image acquisition feasible and accelerates the medical image-processing pipelines.