On the dependence of quantitative diffusion-weighted imaging on scanner system characteristics and acquisition parameters: A large multicenter and multiparametric phantom study with unsupervised clustering analysis

Unsupervised clustering analysis-based characterization of spatial profiles of inaccuracy in apparent diffusion coefficient values with varying acquisition plan orientation and diffusion weighting gradient direction - a large multicenter phantom study

This large multicenter study of 37 magnetic resonance imaging scanners aimed at characterizing, for the first time, spatial profiles of inaccuracy (namely, Δ-profiles) in apparent diffusion coefficient (ADC) values with varying acquisition plan orientation and diffusion weighting gradient direction, using a statistical approach exploiting unsupervised clustering analysis. A diffusion-weighted imaging (DWI) protocol (b-value: 0-200-400-600-800-1000 s mm-2) with different combinations of acquisition plan orientation (axial/sagittal/coronal) and diffusion weighting gradient direction (anterior-posterior/left-right/feet-head) was acquired on a standard water phantom. For each acquisition setup, Δ-profiles along the 3 main orthogonal directions were characterized by fitting data with a second order polynomial function (ar2+ br + c). Moreover, for each Δ-profile, the maximum minus minimum of the fitting function (δmax) was calculated. The parametersa,b,c, andδmaxshowed some significant variations between scanner systems by different manufacturers or with different static magnetic field strengths, as well as between different acquisition/estimation setups. Unsupervised clustering analysis showed two evident clusters with significantly different values of parametera(p< 0.0001), which can be grouped by acquisition protocol/Δ-profile direction but not scanner system. The results of ∆-profiles confirm an appreciable inter-scanner variability in ADC measurement and corroborate the importance of guarantying the reliability of ADC estimations in clinical or research studies, considering for each scanner system the specific acquisition sequence in terms of acquisition plan orientation and diffusion weighting gradient direction.

Reduction of ADC bias in diffusion MRI with deep learning-based acceleration: A phantom validation study at 3.0 T

PURPOSE

Further acceleration of DWI in diagnostic radiology is desired but challenging mainly due to low SNR in high b-value images and associated bias in quantitative ADC values. Deep learning-based reconstruction and denoising may provide a solution to address this challenge.

METHODS

The effects of SNR reduction on ADC bias and variability were investigated using a commercial diffusion phantom and numerical simulations. In the phantom, performance of different reconstruction methods, including conventional parallel (SENSE) imaging, compressed sensing (C-SENSE), and compressed SENSE acceleration with an artificial intelligence deep learning-based technique (C-SENSE AI), was compared at different acceleration factors and flip angles using ROI-based analysis. ADC bias was assessed by Lin's Concordance correlation coefficient (CCC) followed by bootstrapping to calculate confidence intervals (CI). ADC random measurement error (RME) was assessed by the mean coefficient of variation (CV¯) and non-parametric statistical tests.

RESULTS

The simulations predicted increasingly negative bias and loss of precision towards lower SNR. These effects were confirmed in phantom measurements of increasing acceleration, for which CCC decreased from 0.947 to 0.279 and CV¯ increased from 0.043 to 0.439, and of decreasing flip angle, for which CCC decreased from 0.990 to 0.063 and CV¯ increased from 0.037 to 0.508. At high acceleration and low flip angle, C-SENSE AI reconstruction yielded best denoised ADC maps. For the lowest investigated flip angle, CCC = {0.630, 0.771 and 0.987} and CV¯={0.508, 0.426 and 0.254} were obtained for {SENSE, C-SENSE, C-SENSE AI}, the improvement by C-SENSE AI being significant as compared to the other methods (CV: p = 0.033 for C-SENSE AI vs. C-SENSE and p < 0.001 for C-SENSE AI vs. SENSE; CCC: non-overlapping CI between reconstruction methods). For the highest investigated acceleration factor, CCC = {0.479,0.926,0.960} and CV¯={0.519,0.119,0.118} were found, confirming the reduction of bias and RME by C-SENSE AI as compared to C-SENSE (by trend) and to SENSE (CV: p < 0.001; CCC: non-overlapping CI).

CONCLUSION

ADC bias and random measurement error in DWI at low SNR, typically associated with scan acceleration, can be effectively reduced by deep-learning based C-SENSE AI reconstruction.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

A Review on the Use of Imaging Biomarkers in Oncology Clinical Trials: Quality Assurance Strategies for Technical Validation

Imaging biomarkers (IBs) have been proposed in medical literature that exploit images in a quantitative way, going beyond the visual assessment by an imaging physician. These IBs can be used in the diagnosis, prognosis, and response assessment of several pathologies and are very often used for patient management pathways. In this respect, IBs to be used in clinical practice and clinical trials have a requirement to be precise, accurate, and reproducible. Due to limitations in imaging technology, an error can be associated with their value when considering the entire imaging chain, from data acquisition to data reconstruction and subsequent analysis. From this point of view, the use of IBs in clinical trials requires a broadening of the concept of quality assurance and this can be a challenge for the responsible medical physics experts (MPEs). Within this manuscript, we describe the concept of an IB, examine some examples of IBs currently employed in clinical practice/clinical trials and analyze the procedure that should be carried out to achieve better accuracy and reproducibility in their use. We anticipate that this narrative review, written by the components of the EFOMP working group on "the role of the MPEs in clinical trials"-imaging sub-group, can represent a valid reference material for MPEs approaching the subject.

Repeatability of quantitative T1rho magnetic resonance imaging in normal brain tissues at 3.0T

PURPOSE

T1rho imaging is a promising MRI technique for imaging of brain disease. This study aimed to assess the repeatability of quantitative T1rho imaging in the normal brain grey and white matter.

METHODS

The study prospectively recruited 30 healthy volunteers without a history of neurological diseases or brain injury, and T1rho was performed and quantified from three imaging sessions. Repeat measures analysis of variance (ANOVA) and within-subject coefficients of variation (wCoV) was used to detect differences in T1rho values between the three scans.

RESULTS

The results showed low wCoVs of less than 4.3% (range 0.92-4.27%) across all the brain structures. No significant differences were observed in T1rho measurement between the three scans (p > 0.05). The amygdala and hippocampus showed the highest T1rho values of 91.79 ± 2.55 msec and 91.07 ± 2.11 msec respectively, and the palladium and putamen had the lowest values of 67.60 ± 1.84 msec and 71.83 ± 1.85 msec respectively.

CONCLUSION

T1rho shows high test-retest repeatability for whole brain imaging in serial imaging sessions, indicating it to be a reliable sequence for quantitative brain imaging.

Machine learning with textural analysis of longitudinal multiparametric MRI and molecular subtypes accurately predicts pathologic complete response in patients with invasive breast cancer

PURPOSE

To predict pathological complete response (pCR) after neoadjuvant chemotherapy using extreme gradient boosting (XGBoost) with MRI and non-imaging data at multiple treatment timepoints.

MATERIAL AND METHODS

This retrospective study included breast cancer patients (n = 117) who underwent neoadjuvant chemotherapy. Data types used included tumor ADC values, diffusion-weighted and dynamic-contrast-enhanced MRI at three treatment timepoints, and patient demographics and tumor data. GLCM textural analysis was performed on MRI data. An extreme gradient boosting machine learning algorithm was used to predict pCR. Prediction performance was evaluated using the area under the curve (AUC) of the receiver operating curve along with precision and recall.

RESULTS

Prediction using texture features of DWI and DCE images at multiple treatment time points (AUC = 0.871; 95% CI: (0.768, 0.974; p<0.001) and (AUC = 0.903 95% CI: 0.854, 0.952; p<0.001) respectively), outperformed that using mean tumor ADC (AUC = 0.850 (95% CI: 0.764, 0.936; p<0.001)). The AUC using all MRI data was 0.933 (95% CI: 0.836, 1.03; p<0.001). The AUC using non-MRI data was 0.919 (95% CI: 0.848, 0.99; p<0.001). The highest AUC of 0.951 (95% CI: 0.909, 0.993; p<0.001) was achieved with all MRI and all non-MRI data at all time points as inputs.

CONCLUSION

Using XGBoost on extracted GLCM features and non-imaging data accurately predicts pCR. This early prediction of response can minimize exposure to toxic chemotherapy, allowing regimen modification mid-treatment and ultimately achieving better outcomes.

The perfect qMR machine: Measurement variance much less than biological variance

Implementing quantitative MR (qMR) methodology can be a time-consuming task, sometimes seemingly without an end. The concept of the Perfect qMR Machine offers the possibility that the implementation is complete and that no further improvements are needed. This is achieved by making the measurement repeatability variance much less than the biological variance. Thus the proposal is: A Perfect Quantitative MR machine is one that, in making a measurement, contributes no significant extra variation to that which already exists from biological variation. A medal system (platinum, gold, silver and bronze) recognises different sources of biological variance, depending on the type of measurement being carried out (whether a serial study or a group comparison), and different kinds of measurement variance (single machine or multi-centre). A perfect machine can in principle be demonstrated for each quantitative measure (T₁, ADC etc).

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

Tuning the optimal diffusion-weighted MRI parameters on a 0.35-T MR-Linac for clinical implementation: A phantom study

PURPOSE

This study aims to assess the quality of a new diffusion-weighted imaging (DWI) sequence implemented on an MR-Linac MRIdian system, evaluating and optimizing the acquisition parameters to explore the possibility of clinically implementing a DWI acquisition protocol in a 0.35-T MR-Linac.

MATERIALS AND METHODS

All the performed analyses have been carried out on two types of phantoms: a homogeneous 24-cm diameter polymethylmethacrylate (PMMA) sphere (SP) and a homemade phantom (HMP) constating in a PMMA cylinder filled with distilled water with empty sockets into which five cylindrical vials filled with five different concentrations of methylcellulose water solutions have been inserted. SP was used to evaluate the dependence of diffusion gradient inhomogeneity artifacts on gantry position. Four diffusion sequences with b-values of 500 s/mm2 and 3 averages have been acquired: three with diffusion gradients in the three main directions (phase direction, read direction, slice direction) and one with the diffusion gradients switched off. The dependence of diffusion image uniformity and SNR on the number of averages in the MR sequences was also investigated to determine the optimal number of averages. Finally, the ADC values of HMP have been computed and then compared between images acquired in the scanners at 0.35 and 1.5 T.

RESULTS

In order to acquire high-quality artifact-free DWI images, the "slice" gradient direction has been identified to be the optimal one and 0° to be the best gradient angle. Both the SNR ratio and the uniformity increase with the number of averages. A threshold value of 80 for SNR and 85% for uniformity was adopted to choose the best number of averages. By making a compromise between time and quality and limiting the number of b-values, it is possible to reduce the acquisition time to 78 s. The Passing-Bablok test showed that the two methods, with 0.35 and 1.5 T scanners, led to similar results.

CONCLUSION

The quality of the DWI has been accurately evaluated in relation to different sequence parameters, and optimal parameters have been identified to select a clinical protocol for the acquisition of ADC maps sustainable in the workflow of a hybrid radiotherapy system with a 0.35-T MRI scanner.

Clinical translation of quantitative magnetic resonance imaging biomarkers - An overview and gap analysis of current practice

PURPOSE

This overview of the current landscape of quantitative magnetic resonance imaging biomarkers (qMR IBs) aims to support the standardisation of academic IBs to assist their translation to clinical practice.

METHODS

We used three complementary approaches to investigate qMR IB use and quality management practices within the UK: 1) a literature search of qMR and quality management terms during 2011-2015 and 2016-2020; 2) a database search for clinical research studies using qMR IBs during 2016-2020; and 3) a survey to ascertain the current availability and quality management practices for clinical MRI scanners and associated equipment at research institutions across the UK.

RESULTS

The analysis showed increased use of all qMR methods between the periods 2011-2015 and 2016-2020 and diffusion-tensor MRI and volumetry to be popular methods. However, the "translation ratio" of journal articles to clinical research studies was higher for qMR methods that have evidence of clinical translation via a commercial route, such as fat fraction and T2 mapping. The number of journal articles citing quality management terms doubled between the periods 2011-2015 and 2016-2020; although, its proportion relative to all journal articles only increased by 3.0%. The survey suggested that quality assurance (QA) and quality control (QC) of data acquisition procedures are under-reported in the literature and that QA/QC of acquired data/data analysis are under-developed and lack consistency between institutions.

CONCLUSIONS

We summarise current attempts to standardise and translate qMR IBs, and conclude by outlining the ideal quality management practices and providing a gap analysis between current practice and a metrological standard.

Perfusion-Diffusion Ratio: A New IVIM Approach in Differentiating Solid Benign and Malignant Primary Lesions of the Liver

Introduction

In order to improve the efficacy of intravoxel incoherent motion (IVIM) parameters in characterising specific tissues, a new concept is introduced: the perfusion–diffusion ratio (PDR), which expresses the relationship between the signal S(b) decline rate as a result of IVIM and the rate of signal S(b) decline due to diffusion. The aim of this study was to investigate this novel approach in the differentiation of solid primary liver lesions. Material and Methods. Eighty-three patients referred for liver MRI between August 2017 and January 2020 with a suspected liver tumour were prospectively examined with the standard liver MRI protocol extended by DWI-IVIM sequence. Patients with no liver lesions, haemangiomas, or metastases were excluded. The final study population consisted of 34 patients with primary solid liver masses, 9 with FNH, 4 with regenerative nodules, 10 with HCC, and 11 with CCC. The PDR coefficient was introduced, defined as the ratio of the rate of signal S(b) decrease due to the IVIM effect to the rate of signal S(b) decrease due to the diffusion process, for b = 0.

Results

No significant differences were found between benign and malignant lesions in the case of IVIM parameters (f, D, or D^∗) and ADC. Significant differences were observed only for PDR, with lower values for malignant lesions (p = 0.03). The ROC analysis yielded an AUC value for PDR equal to 0.74, with a cut-off value of 5.06, sensitivity of 81%, specificity of 77%, and accuracy of 79%.

Conclusion

PDR proved to be more effective than IVIM parameters and ADC in the differentiation of solid benign and malignant primary liver lesions.

Factors affecting the value of diffusion-weighted imaging for identifying breast cancer patients with pathological complete response on neoadjuvant systemic therapy: a systematic review

This review aims to identify factors causing heterogeneity in breast DWI-MRI and their impact on its value for identifying breast cancer patients with pathological complete response (pCR) on neoadjuvant systemic therapy (NST). A search was performed on PubMed until April 2020 for studies analyzing DWI for identifying breast cancer patients with pCR on NST. Technical and clinical study aspects were extracted and assessed for variability. Twenty studies representing 1455 patients/lesions were included. The studies differed with respect to study population, treatment type, DWI acquisition technique, post-processing (e.g., mono-exponential/intravoxel incoherent motion/stretched exponential modeling), and timing of follow-up studies. For the acquisition and generation of ADC-maps, various b-value combinations were used. Approaches for drawing regions of interest on longitudinal MRIs were highly variable. Biological variability due to various molecular subtypes was usually not taken into account. Moreover, definitions of pCR varied. The individual areas under the curve for the studies range from 0.50 to 0.92. However, overlapping ranges of mean/median ADC-values at pre- and/or during and/or post-NST were found for the pCR and non-pCR groups between studies. The technical, clinical, and epidemiological heterogeneity may be causal for the observed variability in the ability of DWI to predict pCR accurately. This makes implementation of DWI for pCR prediction and evaluation based on one absolute ADC threshold for all breast cancer types undesirable. Multidisciplinary consensus and appropriate clinical study design, taking biological and therapeutic variation into account, is required for obtaining standardized, reliable, and reproducible DWI measurements for pCR/non-pCR identification.

Same Brain, Different Look?-The Impact of Scanner, Sequence and Preprocessing on Diffusion Imaging Outcome Parameters

In clinical diagnostics and longitudinal studies, the reproducibility of MRI assessments is of high importance in order to detect pathological changes, but developments in MRI hard- and software often outrun extended periods of data acquisition and analysis. This could potentially introduce artefactual changes or mask pathological alterations. However, if and how changes of MRI hardware, scanning protocols or preprocessing software affect complex neuroimaging outcomes from, e.g., diffusion weighted imaging (DWI) remains largely understudied. We therefore compared DWI outcomes and artefact severity of 121 healthy participants (age range 19–54 years) who underwent two matched DWI protocols (Siemens product and Center for Magnetic Resonance Research sequence) at two sites (Siemens 3T Magnetom Verio and Skyra^fit). After different preprocessing steps, fractional anisotropy (FA) and mean diffusivity (MD) maps, obtained by tensor fitting, were processed with tract-based spatial statistics (TBSS). Inter-scanner and inter-sequence variability of skeletonised FA values reached up to 5% and differed largely in magnitude and direction across the brain. Skeletonised MD values differed up to 14% between scanners. We here demonstrate that DTI outcome measures strongly depend on imaging site and software, and that these biases vary between brain regions. These regionally inhomogeneous biases may exceed and considerably confound physiological effects such as ageing, highlighting the need to harmonise data acquisition and analysis. Future studies thus need to implement novel strategies to augment neuroimaging data reliability and replicability.

New developments in MRI: System characterization, technical advances and radiotherapy applications

A Special Issue of Physica Medica - European Journal of Medical Physics, focused on some important points of contact between the world of magnetic resonance and that of medical physics, was published during 2021. This Editorial describes and comments on the content of this Focus Issue, which contains articles from leading groups invited by the Guest Editors.