Generalized Richardson-Lucy (GRL) for analyzing multi-shell diffusion MRI data

Improving advanced intraoperative MRI methods during pediatric neurosurgery

Advanced intraoperative MR images (ioMRI) acquired during the resection of pediatric brain tumors could offer additional physiological information to preserve healthy tissue. With this work, we aimed to develop a protocol for ioMRI with increased sensitivity for arterial spin labeling (ASL) and diffusion MRI (dMRI), optimized for patient positioning regularly used in the pediatric neurosurgery setting. For ethical reasons, ASL images were acquired in healthy adult subjects that were imaged in the prone and supine position. After this, the ASL cerebral blood flow (CBF) was quantified and compared between both positions. To evaluate the impact of the RF coils setups on image quality, we compared different setups (two vs. four RF coils) by looking at T1-weighted (T1w) signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR), as well as undertaking a qualitative evaluation of T1w, T2w, ASL, and dMR images. Mean ASL CBF did not differ between the surgical prone and supine positions in any of the investigated regions of interest or the whole brain. T1w SNR (gray matter: p = 0.016, 34% increase; white matter: p = 0.016, 32% increase) and CNR were higher (p = 0.016) in the four versus two RF coils setups (18.0 ± 1.8 vs. 13.9 ± 1.8). Qualitative evaluation of T1w, T2w, ASL, and dMR images resulted in acceptable to good image quality and did not differ statistically significantly between setups. Only the nonweighted diffusion image maps and corticospinal tract reconstructions yielded higher image quality and reduced susceptibility artifacts with four RF coils. Advanced ioMRI metrics were more precise with four RF coils as the standard deviation decreased. Taken together, we have investigated the practical use of advanced ioMRI during pediatric neurosurgery. We conclude that ASL CBF quantification in the surgical prone position is valid and that ASL and dMRI acquisition with two RF coils can be performed adequately for clinical use. With four versus two RF coils, the SNR of the images increases, and the sensitivity to artifacts reduces.

An intracellular isotropic diffusion signal is positively associated with pubertal development in white matter

Puberty is a key event in adolescent development that involves significant, hormone-driven changes to many aspects of physiology including the brain. Understanding how the brain responds during this time period is important for evaluating neuronal developments that affect mental health throughout adolescence and the adult lifespan. This study examines diffusion MRI scans from the cross-sectional ABCD Study baseline cohort, a large multi-site study containing thousands of participants, to describe the relationship between pubertal development and brain microstructure. Using advanced, 3-tissue constrained spherical deconvolution methods, this study is able to describe multiple tissue compartments beyond only white matter (WM) axonal qualities. After controlling for age, sex, brain volume, subject handedness, scanning site, and sibling relationships, we observe a positive relationship between an isotropic, intracellular diffusion signal fraction and pubertal development across a majority of regions of interest (ROIs) in the WM skeleton. We also observe regional effects from an intracellular anisotropic signal fraction compartment and extracellular isotropic free water-like compartment in several ROIs. This cross-sectional work suggests that changes in pubertal status are associated with a complex response from brain tissue that cannot be completely described by traditional methods focusing only on WM axonal properties.

A new open, high-resolution, multishell, diffusion-weighted imaging dataset of the living squirrel monkey

Although very well adapted to brain study, Magnetic Resonance Imaging (MRI) remains limited by the facilities and capabilities required to acquire data, especially for non-human primates. Addressing the data gaps resulting from these limitations requires making data more accessible and open. In contempt of the regular use of Saimiri sciureus in neuroscience research, in vivo diffusion has yet to be openly available for this species. Here we built and made openly available a unique new resource consisting of a high-resolution, multishell diffusion-weighted dataset in the anesthetized Saimiri sciureus. The data were acquired on 11 individuals with an 11.7 T MRI scanner (isotropic resolution of 400 µm³). This paper presents an overview of our dataset and illustrates some of its possible use through example analyses. To assess the quality of our data, we analyzed long-range connections (whole-brain tractography), microstructure (Neurite Orientation Dispersion and Density Imaging), and axon diameter in the corpus callosum (ActiveAx). Constituting an essential new resource for primate evolution studies, all data are openly available.

A transcranial magnetic stimulation study on the role of the left intraparietal sulcus in temporal orienting of attention

Adaptive behavior requires the ability to orient attention to the moment in time at which a relevant event is likely to occur. Temporal orienting of attention has been consistently associated with activation of the left intraparietal sulcus (IPS) in prior fMRI studies. However, a direct test of its causal involvement in temporal orienting is still lacking. The present study tackled this issue by transiently perturbing left IPS activity with either online (Experiment 1) or offline (Experiment 2) transcranial magnetic stimulation (TMS). In both experiments, participants performed a temporal orienting task, alternating between blocks in which a temporal cue predicted when a subsequent target would appear and blocks in which a neutral cue provided no information about target timing. In Experiment 1 we used an online TMS protocol, aiming to interfere specifically with cue-related temporal processes, whereas in Experiment 2 we employed an offline protocol whereby participants performed the temporal orienting task before and after receiving TMS. The right IPS and/or the vertex were stimulated as active control regions. While results replicated the canonical pattern of temporal orienting effects on reaction time, with faster responses for temporal than neutral trials, these effects were not modulated by TMS over the left IPS (as compared to the right IPS and/or vertex regions) regardless of the online or offline protocol used. Overall, these findings challenge the causal role of the left IPS in temporal orienting of attention inviting further research on its underlying neural substrates.

Anatomically informed multi-level fiber tractography for targeted virtual dissection

OBJECTIVES

Diffusion-weighted MRI can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways, such as the corticospinal tract (CST). However, accurate reconstruction of the full extent of the CST remains challenging with existing tractography methods. We suggest a novel tractography algorithm exploiting unused fiber orientations to produce more complete and reliable results.

METHODS

Our novel approach, referred to as multi-level fiber tractography (MLFT), reconstructs fiber pathways by progressively considering previously unused fiber orientations at multiple levels of tract propagation. Anatomical priors are used to minimize the number of false-positive pathways. The MLFT method was evaluated on synthetic data and in vivo data by reconstructing the CST while compared to conventional tractography approaches.

RESULTS

The radial extent of MLFT reconstructions is comparable to that of probabilistic reconstruction: [Formula: see text] for the left and [Formula: see text] for the right hemisphere according to Wilcoxon test, while achieving significantly higher topography preservation compared to probabilistic tractography: [Formula: see text].

DISCUSSION

MLFT provides a novel way to reconstruct fiber pathways by adding the capability of including branching pathways in fiber tractography. Thanks to its robustness, feasible reconstruction extent and topography preservation, our approach may assist in clinical practice as well as in virtual dissection studies.

Advanced intraoperative MRI in pediatric brain tumor surgery

Introduction: In the pediatric brain tumor surgery setting, intraoperative MRI (ioMRI) provides "real-time" imaging, allowing for evaluation of the extent of resection and detection of complications. The use of advanced MRI sequences could potentially provide additional physiological information that may aid in the preservation of healthy brain regions. This review aims to determine the added value of advanced imaging in ioMRI for pediatric brain tumor surgery compared to conventional imaging. Methods: Our systematic literature search identified relevant articles on PubMed using keywords associated with pediatrics, ioMRI, and brain tumors. The literature search was extended using the snowball technique to gather more information on advanced MRI techniques, their technical background, their use in adult ioMRI, and their use in routine pediatric brain tumor care. Results: The available literature was sparse and demonstrated that advanced sequences were used to reconstruct fibers to prevent damage to important structures, provide information on relative cerebral blood flow or abnormal metabolites, or to indicate the onset of hemorrhage or ischemic infarcts. The explorative literature search revealed developments within each advanced MRI field, such as multi-shell diffusion MRI, arterial spin labeling, and amide-proton transfer-weighted imaging, that have been studied in adult ioMRI but have not yet been applied in pediatrics. These techniques could have the potential to provide more accurate fiber tractography, information on intraoperative cerebral perfusion, and to match gadolinium-based T1w images without using a contrast agent. Conclusion: The potential added value of advanced MRI in the intraoperative setting for pediatric brain tumors is to prevent damage to important structures, to provide additional physiological or metabolic information, or to indicate the onset of postoperative changes. Current developments within various advanced ioMRI sequences are promising with regard to providing in-depth tissue information.

Multimodal tract-based MRI metrics outperform whole brain markers in determining cognitive impact of small vessel disease-related brain injury

In cerebral small vessel disease (cSVD), whole brain MRI markers of cSVD-related brain injury explain limited variance to support individualized prediction. Here, we investigate whether considering abnormalities in brain tracts by integrating multimodal metrics from diffusion MRI (dMRI) and structural MRI (sMRI), can better capture cognitive performance in cSVD patients than established approaches based on whole brain markers. We selected 102 patients (73.7 ± 10.2 years old, 59 males) with MRI-visible SVD lesions and both sMRI and dMRI. Conventional linear models using demographics and established whole brain markers were used as benchmark of predicting individual cognitive scores. Multi-modal metrics of 73 major brain tracts were derived from dMRI and sMRI, and used together with established markers as input of a feed-forward artificial neural network (ANN) to predict individual cognitive scores. A feature selection strategy was implemented to reduce the risk of overfitting. Prediction was performed with leave-one-out cross-validation and evaluated with the R² of the correlation between measured and predicted cognitive scores. Linear models predicted memory and processing speed with R² = 0.26 and R² = 0.38, respectively. With ANN, feature selection resulted in 13 tract-specific metrics and 5 whole brain markers for predicting processing speed, and 28 tract-specific metrics and 4 whole brain markers for predicting memory. Leave-one-out ANN prediction with the selected features achieved R² = 0.49 and R² = 0.40 for processing speed and memory, respectively. Our results show proof-of-concept that combining tract-specific multimodal MRI metrics can improve the prediction of cognitive performance in cSVD by leveraging tract-specific multi-modal metrics.

Cross-site harmonization of multi-shell diffusion MRI measures based on rotational invariant spherical harmonics (RISH)

Quantification methods based on the acquisition of diffusion magnetic resonance imaging (dMRI) with multiple diffusion weightings (e.g., multi-shell) are becoming increasingly applied to study the in-vivo brain. Compared to single-shell data for diffusion tensor imaging (DTI), multi-shell data allows to apply more complex models such as diffusion kurtosis imaging (DKI), which attempts to capture both diffusion hindrance and restriction effects, or biophysical models such as NODDI, which attempt to increase specificity by separating biophysical components. Because of the strong dependence of the dMRI signal on the measurement hardware, DKI and NODDI metrics show scanner and site differences, much like other dMRI metrics. These effects limit the implementation of multi-shell approaches in multicenter studies, which are needed to collect large sample sizes for robust analyses. Recently, a post-processing technique based on rotation invariant spherical harmonics (RISH) features was introduced to mitigate cross-scanner differences in DTI metrics. Unlike statistical harmonization methods, which require repeated application to every dMRI metric of choice, RISH harmonization is applied once on the raw data, and can be followed by any analysis. RISH features harmonization has been tested on DTI features but not its generalizability to harmonize multi-shell dMRI. In this work, we investigated whether performing the RISH features harmonization of multi-shell dMRI data removes cross-site differences in DKI and NODDI metrics while retaining longitudinal effects. To this end, 46 subjects underwent a longitudinal (up to 3 time points) two-shell dMRI protocol at 3 imaging sites. DKI and NODDI metrics were derived before and after harmonization and compared both at the whole brain level and at the voxel level. Then, the harmonization effects on cross-sectional and on longitudinal group differences were evaluated. RISH features averaged for each of the 3 sites exhibited prominent between-site differences in the frontal and posterior part of the brain. Statistically significant differences in fractional anisotropy, mean diffusivity and mean kurtosis were observed both at the whole brain and voxel level between all the acquisition sites before harmonization, but not after. The RISH method also proved effective to harmonize NODDI metrics, particularly in white matter. The RISH based harmonization maintained the magnitude and variance of longitudinal changes as compared to the non-harmonized data of all considered metrics. In conclusion, the application of RISH feature based harmonization to multi-shell dMRI data can be used to remove cross-site differences in DKI metrics and NODDI analyses, while retaining inherent relations between longitudinal acquisitions.

Fiber orientation distribution from diffusion MRI: Effects of inaccurate response function calibration

Background and Purpose

Diffusion MRI of the brain enables to quantify white matter fiber orientations noninvasively. Several approaches have been proposed to estimate such characteristics from diffusion MRI data with spherical deconvolution being one of the most widely used methods. Spherical deconvolution requires to define––or derive from the data––a response function, which is used to compute the fiber orientation distribution (FOD). Different characteristics of the response function are expected to affect the FOD computation and the subsequent fiber tracking.

Methods

In this work, we explored the effects of inaccuracies in the shape factors of the response function on the FOD characteristics.

Results

With simulations, we show that the apparent fiber density could be doubled in the presence of underestimated shape factors in the response functions, whereas the overestimation of the shape factor will cause more spurious peaks in the FOD, especially when the signal‐to‐noise ratio is below 15. Moreover, crossing fiber populations with a separation angle smaller than 60° were more sensitive to inaccuracies in the response function than fiber populations with more orthogonal separation angles. Results with in vivo data demonstrate angular deviations in the FODs and spurious peaks as a result of modified shape factors of the response function, while the reconstruction of the main parts of fiber bundles is well preserved.

Conclusions

This work sheds light on how specific aspects of the response function shape can affect the estimated FODs, and highlights the importance of a proper calibration/definition of the response function.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

Constrained spherical deconvolution of nonspherically sampled diffusion MRI data

Constrained spherical deconvolution (CSD) of diffusion-weighted MRI (DW-MRI) is a popular analysis method that extracts the full white matter (WM) fiber orientation density function (fODF) in the living human brain, noninvasively. It assumes that the DW-MRI signal on the sphere can be represented as the spherical convolution of a single-fiber response function (RF) and the fODF, and recovers the fODF through the inverse operation. CSD approaches typically require that the DW-MRI data is sampled shell-wise, and estimate the RF in a purely spherical manner using spherical basis functions, such as spherical harmonics (SH), disregarding any radial dependencies. This precludes analysis of data acquired with nonspherical sampling schemes, for example, Cartesian sampling. Additionally, nonspherical sampling can also arise due to technical issues, for example, gradient nonlinearities, resulting in a spatially dependent bias of the apparent tissue densities and connectivity information. Here, we adopt a compact model for the RFs that also describes their radial dependency. We demonstrate that the proposed model can accurately predict the tissue response for a wide range of b-values. On shell-wise data, our approach provides fODFs and tissue densities indistinguishable from those estimated using SH. On Cartesian data, fODF estimates and apparent tissue densities are on par with those obtained from shell-wise data, significantly broadening the range of data sets that can be analyzed using CSD. In addition, gradient nonlinearities can be accounted for using the proposed model, resulting in much more accurate apparent tissue densities and connectivity metrics.