Diffusion-weighted MRI Signal Research Articles

Assessing tumor microstructure with time-dependent diffusion imaging: Considerations and feasibility on clinical MRI and MRI-Linac.

Quantitative imaging biomarkers (QIBs) can characterize tumor heterogeneity and provide information for biological guidance in radiotherapy (RT). Time-dependent diffusion MRI (TDD-MRI) derived parameters are promising QIBs, as they describe tissue microstructure with more specificity than traditional diffusion-weighted MRI (DW-MRI). Specifically, TDD-MRI can provide information about both restricted diffusion and diffusional exchange, which are the two time-dependent effects affecting diffusion in tissue, and relevant in tumors. However, exhaustive modeling of both effects can require long acquisitions and complex model fitting. Furthermore, several introduced TDD-MRI measurements can require high gradient strengths and/or complex gradient waveforms that are possibly not available in RT settings. In this study, we investigated the feasibility of a simple analysis framework for the detection of restricted diffusion and diffusional exchange effects in the TDD-MRI signal. To promote the clinical applicability, we use standard gradient waveforms on a conventional 1.5 T MRI system with moderate gradient strength (Gmax=45 mT/m), and on a hybrid 1.5 T MRI-Linac system with low gradient strength (Gmax=15 mT/m). Restricted diffusion and diffusional exchange were simulated in geometries mimicking tumor microstructure to investigate the DW-MRI signal behavior and to determine optimal experimental parameters. TDD-MRI was implemented using pulsed field gradient spin echo with the optimized parameters on a conventional MRI system and a MRI-Linac. Experiments in green asparagus and 10 patients with brain lesions were performed to evaluate the time-dependent diffusion (TDD) contrast in the source DW-images. Simulations demonstrated how the TDD contrast was able to differentiate only dominating diffusional exchange in smaller cells from dominating restricted diffusion in larger cells. The maximal TDD contrast in simulations with typical cancer cell sizes and in asparagus measurements exceeded 5% on the conventional MRI but remained below 5% on the MRI-Linac. In particular, the simulated TDD contrast in typical cancer cell sizes (r=5-10µm) remained below or around 2% with the MRI-Linac gradient strength. In patients measured with the conventional MRI, we found sub-regions reflecting either dominating restricted diffusion or dominating diffusional exchange in and around brain lesions compared to the noisy appearing white matter. On the conventional MRI system, the TDD contrast maps showed consistent tumor sub-regions indicating different dominating TDD effects, potentially providing information on the spatial tumor heterogeneity. On the MRI-Linac, the available TDD contrast measured in asparagus showed the same trends as with the conventional MRI but remained close to typical measurement noise levels when simulated in common cancer cell sizes. On conventional MRI systems with moderate gradient strengths, the TDD contrast could potentially be used as a tool to identify which time-dependent effects to include when choosing a biophysical model for more specific tumor characterization.

Case series multilocular cystic hemangioma of the liver: Three cases and literature review.

Multilocular cystic hemangioma is a rare benign tumor classified as an atypical hemangioma. Currently, there are limited imaging reports available, and the imaging characteristics can be challenging to distinguish from other malignant multilocular cystic liver diseases such as cystadenocarcinoma, necessitating confirmation through pathological diagnosis. Here, we discuss the imaging features of 3 cases of multilocular cystic hemangiomas. Case 1 was a 24-year-old young female, and Case 2 involved a 60-year-old elderly male. Both patients were asymptomatic and physical examination revealed hepatic space-occupying lesions. Imaging findings revealed multilocular cystic lesions in the left liver with septa, calcification, a high diffusion-weighted magnetic resonance imaging (DWI) signal at the edge of the lesion, and progressive enhancement of the cyst wall and septa. Case 3 involved a 50-year-old male patient with epigastric distending pain for 1 month and sudden severe abdominal pain for 14 hours. Imaging results revealed a multilocular cystic lesion in the left liver with septa and tumor bleeding, a high DWI signal, and an enhanced cyst wall and septa. The pathological diagnosis confirmed a hepatic hemangioma. All 3 patients underwent liver tumor resection. All 3 patients recovered successfully without any intraoperative or postoperative complications during the follow-up periods of 5 years, 6 months, and 5 months. There were no signs of recurrence. Liver imaging revealed multilocular cystic lesions with features, such as compartmentalization, calcification, or bleeding. Multilocular cystic hemangiomas should be considered in imaging diagnosis. Enhancing our understanding of multilocular cystic hemangiomas can aid in improving the differential diagnosis of other malignant multilocular cystic liver diseases, ultimately reducing unnecessary liver resection.

Exploring white matter dynamics and morphology through interactive numerical phantoms: the White Matter Generator.

Brain white matter is a dynamic environment that continuously adapts and reorganizes in response to stimuli and pathological changes. Glial cells, especially, play a key role in tissue repair, inflammation modulation, and neural recovery. The movements of glial cells and changes in their concentrations can influence the surrounding axon morphology. We introduce the White Matter Generator (WMG) tool to enable the study of how axon morphology is influenced through such dynamical processes, and how this, in turn, influences the diffusion-weighted MRI signal. This is made possible by allowing interactive changes to the configuration of the phantom generation throughout the optimization process. The phantoms can consist of myelinated axons, unmyelinated axons, and cell clusters, separated by extra-cellular space. Due to morphological flexibility and computational advantages during the optimization, the tool uses ellipsoids as building blocks for all structures; chains of ellipsoids for axons, and individual ellipsoids for cell clusters. After optimization, the ellipsoid representation can be converted to a mesh representation which can be employed in Monte-Carlo diffusion simulations. This offers an effective method for evaluating tissue microstructure models for diffusion-weighted MRI in controlled bio-mimicking white matter environments. Hence, the WMG offers valuable insights into white matter's adaptive nature and implications for diffusion-weighted MRI microstructure models, and thereby holds the potential to advance clinical diagnosis, treatment, and rehabilitation strategies for various neurological disorders and injuries.

High-resolution prostate diffusion MRI using eddy current-nulled convex optimized diffusion encoding and random matrix theory-based denoising.

To develop and evaluate a technique combining eddy current-nulled convex optimized diffusion encoding (ENCODE) with random matrix theory (RMT)-based denoising to accelerate and improve the apparent signal-to-noise ratio (aSNR) and apparent diffusion coefficient (ADC) mapping in high-resolution prostate diffusion-weighted MRI (DWI). MATERIALS ANDMETHODS: Eleven subjects with clinical suspicion of prostate cancer were scanned at 3T with high-resolution (HR) (in-plane: 1.0 × 1.0 mm2) ENCODE and standard-resolution (1.6 × 2.2 mm2) bipolar DWI sequences (both had 7 repetitions for averaging, acquisition time [TA] of 5 min 50 s). HR-ENCODE was retrospectively analyzed using three repetitions (accelerated effective TA of 2 min 30 s). The RMT-based denoising pipeline utilized complex DWI signals and Marchenko-Pastur distribution-based principal component analysis to remove additive Gaussian noise in images from multiple coils, b-values, diffusion encoding directions, and repetitions. HR-ENCODE with RMT-based denoising (HR-ENCODE-RMT) was compared with HR-ENCODE in terms of aSNR in prostate peripheral zone (PZ) and transition zone (TZ). Precision and accuracy of ADC were evaluated by the coefficient of variation (CoV) between repeated measurements and mean difference (MD) compared to the bipolar ADC reference, respectively. Differences were compared using two-sided Wilcoxon signed-rank tests (P < 0.05 considered significant). HR-ENCODE-RMT yielded 62% and 56% higher median aSNR than HR-ENCODE (b = 800 s/mm2) in PZ and TZ, respectively (P < 0.001). HR-ENCODE-RMT achieved 63% and 70% lower ADC-CoV than HR-ENCODE in PZ and TZ, respectively (P < 0.001). HR-ENCODE-RMT ADC and bipolar ADC had low MD of 22.7 × 10-6 mm2/s in PZ and low MD of 90.5 × 10-6 mm2/s in TZ. HR-ENCODE-RMT can shorten the acquisition time and improve the aSNR of high-resolution prostate DWI and achieve accurate and precise ADC measurements in the prostate.

Robust fiber orientation distribution function estimation using deep constrained spherical deconvolution for diffusion-weighted magnetic resonance imaging.

Diffusion-weighted magnetic resonance imaging (DW-MRI) is a critical imaging method for capturing and modeling tissue microarchitecture at a millimeter scale. A common practice to model the measured DW-MRI signal is via fiber orientation distribution function (fODF). This function is the essential first step for the downstream tractography and connectivity analyses. With recent advantages in data sharing, large-scale multisite DW-MRI datasets are being made available for multisite studies. However, measurement variabilities (e.g., inter- and intrasite variability, hardware performance, and sequence design) are inevitable during the acquisition of DW-MRI. Most existing model-based methods [e.g., constrained spherical deconvolution (CSD)] and learning-based methods (e.g., deep learning) do not explicitly consider such variabilities in fODF modeling, which consequently leads to inferior performance on multisite and/or longitudinal diffusion studies. In this paper, we propose a data-driven deep CSD method to explicitly constrain the scan-rescan variabilities for a more reproducible and robust estimation of brain microstructure from repeated DW-MRI scans. Specifically, the proposed method introduces a three-dimensional volumetric scanner-invariant regularization scheme during the fODF estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter- and intrasite scan/rescan data). The Baltimore Longitudinal Study of Aging dataset is employed for external validation. From the experimental results, the proposed data-driven framework outperforms the existing benchmarks in repeated fODF estimation. By introducing the contrastive loss with scan/rescan data, the proposed method achieved a higher consistency while maintaining higher angular correlation coefficients with the CSD modeling. The proposed method is assessing the downstream connectivity analysis and shows increased performance in distinguishing subjects with different biomarkers. We propose a deep CSD method to explicitly reduce the scan-rescan variabilities, so as to model a more reproducible and robust brain microstructure from repeated DW-MRI scans. The plug-and-play design of the proposed approach is potentially applicable to a wider range of data harmonization problems in neuroimaging.

Deep Constrained Spherical Deconvolution for Robust Harmonization.

Diffusion weighted magnetic resonance imaging (DW-MRI) captures tissue microarchitecture at millimeter scale. With recent advantages in data sharing, large-scale multi-site DW-MRI datasets are being made available for multi-site studies. However, DW-MRI suffers from measurement variability (e.g., inter- and intra-site variability, hardware performance, and sequence design), which consequently yields inferior performance on multi-site and/or longitudinal diffusion studies. In this study, we propose a novel, deep learning-based method to harmonize DW-MRI signals for a more reproducible and robust estimation of microstructure. Our method introduces a data-driven scanner-invariant regularization scheme to model a more robust fiber orientation distribution function (FODF) estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter- and intra-site scan/rescan data). The 8th order spherical harmonics coefficients are employed as data representation. The results show that the proposed harmonization approach maintains higher angular correlation coefficients (ACC) with the ground truth signals (0.954 versus 0.942), while achieves higher consistency of FODF signals for intra-scanner data (0.891 versus 0.826), as compared with the baseline supervised deep learning scheme. Furthermore, the proposed data-driven framework is flexible and potentially applicable to a wider range of data harmonization problems in neuroimaging.

Tri-Compartmental Restriction Spectrum Imaging Breast Model Distinguishes Malignant Lesions from Benign Lesions and Healthy Tissue on Diffusion-Weighted Imaging.

Simple SummaryWomen at high risk for breast cancer are regularly screened using contrast-enhanced magnetic resonance imaging (MRI) to identify potential malignancy. Diffusion-weighted MRI (DW-MRI) is a non-contrast technique that measures the differential movement of water molecules in tissues and is sensitive to cancer cells. In this study, we use a multi-exponential advanced DW-MRI model called restriction spectrum imaging (RSI) to characterize the diffusion characteristics of malignant lesions, benign lesions, and healthy breast tissue to help differentiate benign from malignant disease. In a cohort of patients, we show that cancer exhibits more restricted diffusion compared to benign breast lesions and healthy tissue, whereas benign and healthy tissue are not different from each other.Diffusion-weighted MRI (DW-MRI) offers a potential adjunct to dynamic contrast-enhanced MRI to discriminate benign from malignant breast lesions by yielding quantitative information about tissue microstructure. Multi-component modeling of the DW-MRI signal over an extended b-value range (up to 3000 s/mm2) theoretically isolates the slowly diffusing (restricted) water component in tissues. Previously, a three-component restriction spectrum imaging (RSI) model demonstrated the ability to distinguish malignant lesions from healthy breast tissue. We further evaluated the utility of this three-component model to differentiate malignant from benign lesions and healthy tissue in 12 patients with known malignancy and synchronous pathology-proven benign lesions. The signal contributions from three distinct diffusion compartments were measured to generate parametric maps corresponding to diffusivity on a voxel-wise basis. The three-component model discriminated malignant from benign and healthy tissue, particularly using the restricted diffusion C1 compartment and product of the restricted and intermediate diffusion compartments (C1 and C2). However, benign lesions and healthy tissue did not significantly differ in diffusion characteristics. Quantitative discrimination of these three tissue types (malignant, benign, and healthy) in non-pre-defined lesions may enhance the clinical utility of DW-MRI in reducing excessive biopsies and aiding in surveillance and surgical evaluation without repeated exposure to gadolinium contrast.

Neoadjuvant Chemoradiotherapy for Oral Cavity Cancer: Predictive Factors for Response and Interim Analysis of the Prospective INVERT-Trial.

BackgroundTo study neoadjuvant chemoradiotherapy (nCRT) and potential predictive factors for response in locally advanced oral cavity cancer (LA-OCC).MethodsThe INVERT trial is an ongoing single-center, prospective phase 2, proof-of-principle trial. Operable patients with stage III-IVA squamous cell carcinomas of the oral cavity were eligible and received nCRT consisting of 60 Gy with concomitant cisplatin and 5-fluorouracil. Surgery was scheduled 6-8 weeks after completion of nCRT. Explorative, multiplex immunohistochemistry (IHC) was performed on pretreatment tumor specimen, and diffusion-weighted magnetic resonance imaging (DW-MRI) was conducted prior to, during nCRT (day 15), and before surgery to identify potential predictive biomarkers and imaging features. Primary endpoint was the pathological complete response (pCR) rate.ResultsSeventeen patients with stage IVA OCC were included in this interim analysis. All patients completed nCRT. One patient died from pneumonia 10 weeks after nCRT before surgery. Complete tumor resection (R0) was achieved in 16/17 patients, of whom 7 (41%, 95% CI: 18-67%) showed pCR. According to the Clavien-Dindo classification, grade 3a and 3b complications were found in 4 (25%) and 5 (31%) patients, respectively; grade 4-5 complications did not occur. Increased changes in the apparent diffusion coefficient signal intensities between MRI at day 15 of nCRT and before surgery were associated with better response (p=0.022). Higher abundances of programmed cell death protein 1 (PD1) positive cytotoxic T-cells (p=0.012), PD1+ macrophages (p=0.046), and cancer-associated fibroblasts (CAFs, p=0.036) were associated with incomplete response to nCRT.ConclusionnCRT for LA-OCC followed by radical surgery is feasible and shows high response rates. Larger patient cohorts from randomized trials are needed to further investigate nCRT and predictive biomarkers such as changes in DW-MRI signal intensities, tumor infiltrating immune cells, and CAFs.

Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology.

PurposeRelationships between diffusion‐weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange? MethodsMonte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L/D0, and clinical protocols (single diffusion encoding; maximum b‐value: {1000, 1500, 2000} s/mm2; 5 unique gradient duration/separation pairs; SNR = {∞, 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion (D) and kurtosis (K) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell‐size contrasts, and deployed on 9.4T ex vivo MRI‐histology data of fixed mouse liversResultsRelationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three‐class cell‐size classification task (for SNR = 20, bmax = 1500 s/mm2, δ = 20 ms, and Δ = 75 ms). MRI detects cell‐size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes.ConclusionSalient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.

Characterization of the diffusion signal of breast tissues using multi-exponential models.

Restriction spectrum imaging (RSI) decomposes the diffusion-weighted MRI signal into separate components of known apparent diffusion coefficients (ADCs). The number of diffusion components and optimal ADCs for RSI are organ-specific and determined empirically. The purpose of this work was to determine the RSI model for breast tissues. The diffusion-weighted MRI signal was described using a linear combination of multiple exponential components. A set of ADC values was estimated to fit voxels in cancer and control ROIs. Later, the signal contributions of each diffusion component were estimated using these fixed ADC values. Relative-fitting residuals and Bayesian information criterion were assessed. Contrast-to-noise ratio between cancer and fibroglandular tissue in RSI-derived signal contribution maps was compared to DCE imaging. A total of 74 women with breast cancer were scanned at 3.0 Tesla MRI. The fitting residuals of conventional ADC and Bayesian information criterion suggest that a 3-component model improves the characterization of the diffusion signal over a biexponential model. Estimated ADCs of triexponential model were D1,3 = 0, D2,3 = 1.5 × 10-3 , and D3,3 = 10.8 × 10-3 mm2 /s. The RSI-derived signal contributions of the slower diffusion components were larger in tumors than in fibroglandular tissues. Further, the contrast-to-noise and specificity at 80% sensitivity of DCE and a subset of RSI-derived maps were equivalent. Breast diffusion-weighted MRI signal was best described using a triexponential model. Tumor conspicuity in breast RSI model is comparable to that of DCE without the use of exogenous contrast. These data may be used as differential features between healthy and malignant breast tissues.

Analysis of tumour microstructure estimation from conventional diffusion MRI and application to skull-base chordoma.

Skull-base chordoma (SBC) is a rare tumour whose molecular and radiological characteristics are still being investigated. In neuro-oncology microstructural imaging techniques, like diffusion-weighted MRI (DW-MRI), have been widely investigated, with the apparent diffusion coefficient (ADC) being one of the most used DW-MRI parameters due to its ease of acquisition and computation. ADC is a potential biomarker without a clear link to microstructure. The aim of this work was to derive microstructural information from conventional ADC, showing its potential for the characterisation of skull-base chordomas. Sixteen patients affected by SBC, who underwent conventional DW-MRI were retrospectively selected. From mono-exponential fits of DW-MRI, ADC maps were estimated using different sets of b-values. DW-MRI signals were simulated from synthetic substrates , which mimic the cellular packing of a tumour tissue with well-defined microstructural features. Starting from a published method, an error-driven procedure was evaluated to improve the estimates of microstructural parameters obtained through the simulated signals. A quantitative description of the tumour microstructure was then obtained from the DW-MRI images. This allowed successfully differentiating patients according to histologically-verified cell proliferation information.Clinical Relevance - The impact on cancer management derives from the expected improvement of radiation treatment quality tailored to a patient-specific non-invasive description of tumour microstructure.

Genetic and environmental influences of variation in diffusion MRI measures of white matter microstructure.

Quantitative neuroimaging studies in twin samples can investigate genetic contributions to brain structure and microstructure. Diffusion tensor imaging (DTI) studies with twin samples have shown moderate to high heritability in white matter microstructure. This study investigates the genetic and environmental contributions of another widely used diffusion MRI model not yet applied to twin studies, neurite orientation dispersion and density imaging (NODDI). The NODDI model is a multicompartment model of the diffusion-weighted MRI signal, providing estimates of neurite density (ND) and the orientation dispersion index (ODI). A cohort of monozygotic (MZ) and same-sex dizygotic (DZ) twins (N = 460 individuals) between 13 and 24years of age were scanned with a multi-shell diffusion weighted imaging protocol. Select white matter (WM) regions of interest (ROI) were extracted. Biometric structural equation modeling estimated the relative contributions from additive genetic (A) and common (C) and unique environmental (E) factors. Genetic factors for the NODDI measures accounted for 91% and 65% of the variation of global ND and ODI, respectively, compared with 83% for FA. We observed higher heritability for ND than both FA and ODI in 25 of 30 discrete white matter regions that we examined, suggesting ND may be more sensitive to underlying genetic sources of variation. This study demonstrated that genetic factors play a key role in the development of white matter microstructure using both DTI and NODDI.

Assessment of Magnetic Resonance Imaging Changes and Functional Outcomes Among Adults With Severe Herpes Simplex Encephalitis

Current guidelines recommend brain magnetic resonance imaging (MRI) for clinical management of patients with severe herpes simplex encephalitis (HSE). However, the prognostic value of brain imaging has not been demonstrated in this setting. To investigate the association between early brain MRI data and functional outcomes of patients with HSE at 90 days after intensive care unit (ICU) admission. This multicenter cohort study was conducted in 34 ICUs in France from 2007 to 2019 and recruited all patients who received a clinical diagnosis of encephalitis and exhibited cerebrospinal fluid positivity for herpes simplex virus DNA in the polymerase chain reaction analysis. Data analysis was performed from January to April 2020. All patients underwent a standard brain MRI during the first 30 days after ICU admission. MRI acquisitions were analyzed by radiologists blinded to patients' outcomes, using a predefined score. Multivariable logistic regression and supervised hierarchical classifiers methods were used to identify factors associated with poor outcome at 90 days, defined by a score of 3 to 6 (indicating moderate-to-severe disability or death) on the Modified Rankin Scale. Overall, 138 patients (median [interquartile range {IQR}] age, 62.6 [54.0-72.0] years; 75 men [54.3%]) with an admission median (IQR) Glasgow Coma Scale score of 9 (6-12) were studied. The median (IQR) delay between ICU admission and MRI was 1 (1-7) days. At 90 days, 95 patients (68.8%) had a poor outcome, including 16 deaths (11.6%). The presence of fluid-attenuated inversion recovery MRI signal abnormalities in more than 3 brain lobes (odds ratio [OR], 25.71; 95% CI, 1.21-554.42), age older than 60 years (OR, 7.62; 95% CI, 2.02-28.91), and the presence of diffusion-weighted MRI signal abnormalities in the left thalamus (OR, 6.90; 95% CI, 1.12-43.00) were independently associated with poor outcome. Machine learning models identified bilateral diffusion abnormalities as an additional factor associated with poor outcome (34 of 39 patients [87.2%] with bilateral abnormalities had poor outcomes) and confirmed the functional burden of left thalamic lesions, particularly in older patients (all 11 patients aged >60 years had left thalamic lesions). These findings suggest that in adult patients with HSE requiring ICU admission, extensive MRI changes in the brain are independently associated with poor functional outcome at 90 days. Thalamic diffusion signal changes were frequently observed and were associated with poor prognosis, mainly in older patients.

Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range.

It has been increasingly reported that in biological tissues diffusion-weighted MRI signal attenuation deviates from mono-exponential decay, especially at high b-values. A number of diffusion models have been proposed to characterize this non-Gaussian diffusion behavior. One of these models is the continuous-time random-walk (CTRW) model, which introduces two new parameters: a fractional order time derivative α and a fractional order spatial derivative β. These new parameters have been linked to intravoxel diffusion heterogeneities in time and space, respectively, and are believed to depend on diffusion times. Studies on this time dependency are limited, largely because the diffusion time cannot vary over a board range in a conventional spin-echo echo-planar imaging sequence due to the accompanying T2 decays. In this study, we investigated the time-dependency of the CTRW model in Sephadex gel phantoms across a broad diffusion time range by employing oscillating-gradient spin-echo, pulsed-gradient spin-echo, and pulsed-gradient stimulated echo sequences. We also performed Monte Carlo simulations to help understand our experimental results. It was observed that the diffusion process fell into the Gaussian regime at extremely short diffusion times whereas it exhibited a strong time dependency in the CTRW parameters at longer diffusion times.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called ’b-tensor’ encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Understanding diffusion-weighted MRI analysis: Repeatability and performance of diffusion models in a benign breast lesion cohort.

Diffusion-weighted MRI (DWI) is an important tool for oncology research, with great clinical potential for the classification and monitoring of breast lesions. The utility of parameters derived from DWI, however, is influenced by specific analysis choices. The purpose of this study was to critically evaluate repeatability and curve-fitting performance of common DWI signal representations, for a prospective cohort of patients with benign breast lesions. Twenty informed, consented patients with confirmed benign breast lesions underwent repeated DWI (3 T) using: sagittal single-shot spin-echo echo planar imaging, bipolar encoding, TR/TE: 11,600/86 ms, FOV: 180 x 180 mm, matrix: 90 x 90, slices: 60 x 2.5 mm, iPAT: GRAPPA 2, fat suppression, and 13 b-values: 0-700 s/mm2 . A phase-reversed scan (b = 0 s/mm2 ) was acquired for distortion correction. Voxel-wise repeat-measures coefficients of variation (CoVs) were derived for monoexponential (apparent diffusion coefficient [ADC]), biexponential (intravoxel incoherent motion: f, D, D*) and stretched exponential (α, DDC) across the parameter histograms for lesion regions of interest (ROIs). Goodness-of-fit for each representation was assessed by Bayesian information criterion. The volume of interest (VOI) definition was repeatable (CoV 13.9%). Within lesions, and across both visits and the cohort, there was no dominant best-fit model, with all representations giving the best fit for a fraction of the voxels. Diffusivity measures from the signal representations (ADC, D, DDC) all showed good repeatability (CoV < 10%), whereas parameters associated with pseudodiffusion (f, D*) performed poorly (CoV > 50%). The stretching exponent α was repeatable (CoV < 12%). This pattern of repeatability was consistent over the central part of the parameter percentiles. Assumptions often made in diffusion studies about analysis choices will influence the detectability of changes, potentially obscuring useful information. No single signal representation prevails within or across lesions, or across repeated visits; parameter robustness is therefore a critical consideration. Our results suggest that stretched exponential representation is more repeatable than biexponential, with pseudodiffusion parameters unlikely to provide clinically useful biomarkers.

Improving the characterization of meningioma microstructure in proton therapy from conventional apparent diffusion coefficient measurements using Monte Carlo simulations of diffusion MRI.

Proton therapy could benefit from noninvasively gaining tumor microstructure information, at both planning and monitoring stages. The anatomical location of brain tumors, such as meningiomas, often hinders the recovery of such information from histopathology, and conventional noninvasive imaging biomarkers, like the apparent diffusion coefficient (ADC) from diffusion-weighted MRI (DW-MRI), are nonspecific. The aim of this study was to retrieve discriminative microstructural markers from conventional ADC for meningiomas treated with proton therapy. These markers were employed for tumor grading and tumor response assessment. DW-MRIs from patients affected by meningioma and enrolled in proton therapy were collected before (n=35) and 3 months after (n=25) treatment. For the latter group, the risk of an adverse outcome was inferred by their clinical history. Using Monte Carlo methods, DW-MRI signals were simulated from packings of synthetic cells built with well-defined geometrical and diffusion properties. Patients' ADC was modeled as a weighted sum of selected simulated signals. The weights that best described a patient's ADC were determined through an optimization procedure and used to estimate a set of markers of tumor microstructure: diffusion coefficient (D), volume fraction (vf), and radius (R). Apparent cellularity (ρapp ) was estimated from vf and R for an easier clinical interpretability. Differences between meningothelial and atypical subtypes, and low- and high-grade meningiomas were assessed with nonparametric statistical tests, whereas sensitivity and specificity with ROC analyses. Similar analyses were performed for patients showing low or high risk of an adverse outcome to preliminary evaluate response to treatment. Significant (P<0.05) differences in median ADC, D, vf, R, and ρapp values were found when comparing meningiomas' subtypes and grades. ROC analyses showed that estimated microstructural parameters reached higher specificity than ADC for subtyping (0.93 for D and vf vs 0.80 for ADC) and grading (0.75 for R vs 0.67 for ADC). High- and low-risk patients showed significant differences in ADC and microstructural parameters. The skewness of ρapp was the parameter with highest AUC (0.90) and sensitivity (0.75). Matching measured with simulated ADC yielded a set of potential imaging markers for meningiomas grading and response monitoring in proton therapy, showing higher specificity than conventional ADC. These markers can provide discriminative information about spatial patterns of tumor microstructure implying important advantages for patient-specific proton therapy workflows.

Machine learning based white matter models with permeability: An experimental study in cuprizone treated in-vivo mouse model of axonal demyelination

The intra-axonal water exchange time (τi), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology.We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R2 =0.99, τi: R2 =0.84, intrinsic diffusivity d: R2 =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τi). Finally, we find a statistically significant decrease in τi in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τi>=310ms/330ms/350ms) compared to the WT group (<τi>=370ms/370ms/380ms). This is in line with our expectations that τi is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .

Diffusion MRI Reversibility in Ischemic Stroke Following Thrombolysis: A Meta-Analysis.

Diffusion-weighted magnetic resonance imaging (DWI) detects early infarction in acute stroke. With the substantial progress in stroke therapies, the frequency of posttreatment DWI reversibility in modern stroke cohorts is currently unknown. The purpose of this study was to perform a systematic literature review examining the relationship between characteristics of patients with ischemic stroke and DWI reversibility following treatment with lytic therapy. A systematic review was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines, yielding a total of 422 unique articles. Studies that were nonclinical or did not report data pertaining to DWI reversibility in the context of an acute stroke series were excluded. Characteristics regarding presentation, diagnosis, intervention, and the timing of DWI reversibility were collected for each study. After full-text review, 10 studies were identified as meeting inclusion criteria. The number of patients with DWI reversal ranged from .9% to 50%, whereas the extent of reversal ranged from 1.8% to 72.7%. Studies reporting on younger patients describe greater rates of reversibility following stroke treatment. These data suggest that early DWI signal may not represent the definitive DWI burden in recanalized populations. However, substantial heterogeneity exists regarding the rate of DWI reversal following recanalization. Additional studies are needed to elucidate the relationship among time to treatment, early reversal rates, and clinical outcomes. Physicians should use caution when basing clinical decisions on DWI lesion volumes, as these likely change to some degree with recanalization.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b ​≤ ​3,000 ​s/mm2 (or 3 ​ms/μm2), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm2 or 3 ​ms/μm2). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (bmax ​= ​40,000 ​s/mm2, or 40 ​ms/μm2) and human (bmax ​= ​10,000 ​s/mm2, or 10 ​ms/μm2) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.