Advances in MRI‐guided precision radiotherapy

Reconstruction of multi-phase parametric maps in 4D-magnetic resonance fingerprinting (4D-MRF) by optimization of local T1 and T2 sensitivities

BACKGROUND

Time-resolved magnetic resonance fingerprinting (MRF), or 4D-MRF, has been demonstrated its feasibility in motion management in radiotherapy (RT). However, the prohibitive long acquisition time is one of challenges of the clinical implementation of 4D-MRF. The shortening of acquisition time causes data insufficiency in each respiratory phase, leading to poor accuracies and consistencies of the predicted tissues' properties of each phase.

PURPOSE

To develop a technique for the reconstruction of multi-phase parametric maps in four-dimensional magnetic resonance fingerprinting (4D-MRF) through the optimization of local T1 and T2 sensitivities.

METHODS

The proposed technique employed an iterative optimization to tailor the data arrangement of each phase by manipulation of inter-phase frames, such that the T1 and T2 sensitivities, which were quantified by the modified Minkowski distance, of the truncated signal evolution curve was maximized. The multi-phase signal evolution curves were modified by sliding window reconstruction and inter-phase frame sharing (SWIFS). Motion correction (MC) and dot product matching were sequentially performed on the modified signal evolution and dictionary to reconstruct the multi-parametric maps. The proposed technique was evaluated by numerical simulations using the extended cardiac-torso (XCAT) phantom with regular and irregular breathing patterns, and by in vivo MRF data of three health volunteers and six liver cancer patients acquired at a 3.0 T scanner.

RESULTS

In simulation study, the proposed SWIFS approach achieved the overall mean absolute percentage error (MAPE) of 8.62% ± 1.59% and 16.2% ± 3.88% for the eight-phases T1 and T2 maps, respectively, in the sagittal view with irregular breathing patterns. In contrast, the overall MAPE of T1 and T2 maps generated by the conventional approach with multiple MRF repetitions were 22.1% ± 11.0% and 30.8% ± 14.9%, respectively. For in-vivo study, the predicted mean T1 and T2 of liver by the proposed SWIFS approach were 795 ms ± 38.9 ms and 58.3 ms ± 11.7 ms, respectively.

CONCLUSIONS

Both simulation and in vivo results showed that the approach empowered by T1 and T2 sensitivities optimization and sliding window under the shortened acquisition of MRF had superior performance in the estimation of multi-phase T1 and T2 maps as compared to the conventional approach with oversampling of MRF data.

Is it necessary to perform measurement-based patient-specific quality assurance for online adaptive radiotherapy with Elekta Unity MR-Linac?

This study aimed to investigate the necessity of measurement-based patient-specific quality assurance (PSQA) for online adaptive radiotherapy by analyzing measurement-based PSQA results and calculation-based 3D independent dose verification results with Elekta Unity MR-Linac. There are two workflows for Elekta Unity enabled in the treatment planning system: adapt to position (ATP) and adapt to shape (ATS). ATP plans are those which have relatively slighter shifts from reference plans by adjusting beam shapes or weights, whereas ATS plans are the new plans optimized from the beginning with probable re-contouring targets and organs-at-risk. PSQA gamma passing rates were measured using an MR-compatible ArcCHECK diode array for 78 reference plans and corresponding 208 adaptive plans (129 ATP plans and 79 ATS plans) of Elekta Unity. Subsequently, the relationships between ATP, or ATS plans and reference plans were evaluated separately. The Pearson's r correlation coefficients between ATP or ATS adaptive plans and corresponding reference plans were also characterized using regression analysis. Moreover, the Bland-Altman plot method was used to describe the agreement of PSQA results between ATP or ATS adaptive plans and reference plans. Additionally, Monte Carlo-based independent dose verification software ArcherQA was used to perform secondary dose check for adaptive plans. For ArcCHECK measurements, the average gamma passing rates (ArcCHECK vs. TPS) of PSQA (3%/2 mm criterion) were 99.51% ± 0.88% and 99.43% ± 0.54% for ATP and ATS plans, respectively, which were higher than the corresponding reference plans 99.34% ± 1.04% (p < 0.05) and 99.20% ± 0.71% (p < 0.05), respectively. The Pearson's r correlation coefficients were 0.720 between ATP and reference plans and 0.300 between ATS and reference plans with ArcCHECK, respectively. Furthermore, >95% of data points of differences between both ATP and ATS plans and reference plans were within ±2σ (standard deviation) of the mean difference between adaptive and reference plans with ArcCHECK measurements. With ArcherQA calculation, the average gamma passing rates (ArcherQA vs. TPS) were 98.23% ± 1.64% and 98.15% ± 1.07% for ATP and ATS adaptive plans, separately. It might be unnecessary to perform measurement-based PSQA for both ATP and ATS adaptive plans for Unity if the gamma passing rates of both measurements of corresponding reference plans and independent dose verification of adaptive plans have high gamma passing rates. Periodic machine QA and verification of adaptive plans were recommended to ensure treatment safety.

Implementation of a retrofit MRI simulator for radiation therapy planning

Magnetic resonance imaging (MRI) is being integrated into routine radiation therapy (RT) planning workflows. To reap the benefits of this imaging modality, patient positioning, image acquisition parameters and a quality assurance programme must be considered for accurate use. This paper will report on the implementation of a retrofit MRI Simulator for RT planning, demonstrating an economical, resource efficient solution to improve the accuracy of MRI in this setting.

A multi-task deep learning model for EGFR genotyping prediction and GTV segmentation of brain metastasis

BACKGROUND

The precise prediction of epidermal growth factor receptor (EGFR) mutation status and gross tumor volume (GTV) segmentation are crucial goals in computer-aided lung adenocarcinoma brain metastasis diagnosis. However, these two tasks present continuous difficulties due to the nonuniform intensity distributions, ambiguous boundaries, and variable shapes of brain metastasis (BM) in MR images.The existing approaches for tackling these challenges mainly rely on single-task algorithms, which overlook the interdependence between these two tasks.

METHODS

To comprehensively address these challenges, we propose a multi-task deep learning model that simultaneously enables GTV segmentation and EGFR subtype classification. Specifically, a multi-scale self-attention encoder that consists of a convolutional self-attention module is designed to extract the shared spatial and global information for a GTV segmentation decoder and an EGFR genotype classifier. Then, a hybrid CNN-Transformer classifier consisting of a convolutional block and a Transformer block is designed to combine the global and local information. Furthermore, the task correlation and heterogeneity issues are solved with a multi-task loss function, aiming to balance the above two tasks by incorporating segmentation and classification loss functions with learnable weights.

RESULTS

The experimental results demonstrate that our proposed model achieves excellent performance, surpassing that of single-task learning approaches. Our proposed model achieves a mean Dice score of 0.89 for GTV segmentation and an EGFR genotyping accuracy of 0.88 on an internal testing set, and attains an accuracy of 0.81 in the EGFR genotype prediction task and an average Dice score of 0.85 in the GTV segmentation task on the external testing set. This shows that our proposed method has outstanding performance and generalization.

CONCLUSION

With the introduction of an efficient feature extraction module, a hybrid CNN-Transformer classifier, and a multi-task loss function, the proposed multi-task deep learning network significantly enhances the performance achieved in both GTV segmentation and EGFR genotyping tasks. Thus, the model can serve as a noninvasive tool for facilitating clinical treatment.

Towards simulation-free MR-linac treatment: utilizing male pelvis PSMA-PET/CT and population-based electron density assignments

Objective. This study aimed to investigate the dosimetric impact of using population-based relative electron density (RED) overrides in lieu of simulation computerized tomography (CT) in a magnetic resonance linear accelerator (MRL) workflow for male pelvis patients. Additionally, the feasibility of using prostate specific membrane antigen positron emission tomography/CT (PSMA-PET/CT) scans to assess patients' eligibility for this proposed workflow was examined.Approach. In this study, 74 male pelvis patients treated on an Elekta Unity 1.5 T MRL were retrospectively selected. The patients' individual RED values for 8 organs of interest were extracted from their simulation-CT images to establish population-based RED values. These values were used to generate individual (IndD) and population-based (PopD) RED dose plans, representing current and proposed MRL workflows, respectively. Lastly, this study compared RED values obtained from CT and PET-CT scanners in a phantom and a subset of patients.Results. Population-based RED values were mostly within two standard deviations of ICRU Report 46 values. PopD plans were comparable to IndD plans, with the average %difference magnitudes of 0.5%, 0.6%, and 0.6% for mean dose (all organs), D0.1cm3(non-target organs) and D95%/D98% (target organs), respectively. Both phantom and patient PET-CT derived RED values had high agreement with corresponding CT-derived values, with correlation coefficients ≥ 0.9.Significance. Population-based RED values were considered suitable in a simulation-free MRL treatment workflow. Utilizing these RED values resulted in similar dosimetric uncertainties as per the current workflow. Initial findings also suggested that PET-CT scans may be used to assess prospective patients' eligibility for the proposed workflow. Future investigations will evaluate the clinical feasibility of implementing this workflow for prospective patients in the clinical setting. This is aimed to reduce patient burden during radiotherapy and increase department efficiencies.

Radiomic feature repeatability and its impact on prognostic model generalizability: A multi-institutional study on nasopharyngeal carcinoma patients

BACKGROUND AND PURPOSE

To investigate the radiomic feature (RF) repeatability via perturbation and its impact on cross-institutional prognostic model generalizability in Nasopharyngeal Carcinoma (NPC) patients.

MATERIALS AND METHODS

286 and 183 NPC patients from two institutions were included for model training and validation. Perturbations with random translations and rotations were applied to contrast-enhanced T1-weighted (CET1-w) MR images. RFs were extracted from primary tumor volume under a wide range of image filtering and discretization settings. RF repeatability was assessed by intraclass correlation coefficient (ICC), which was used to equally separate the RFs into low- and high-repeatable groups by the median value. After feature selection, multivariate Cox regression and Kaplan-Meier analysis were independently employed to develop and analyze prognostic models. Concordance index (C-index) and P-value from log-rank test were used to assess model performance.

RESULTS

Most textural RFs from high-pass wavelet-filtered images were susceptible to image perturbations. It was more prominent when a smaller discretization bin number was used (e.g., 8, mean ICC = 0.69). Using high-repeatable RFs for model development yielded a significantly higher C-index (0.63) in the validation cohort than when only low-repeatable RFs were used (0.57, P = 0.024), suggesting higher model generalizability. Besides, significant risk stratification in the validation cohort was observed only when high-repeatable RFs were used (P < 0.001).

CONCLUSION

Repeatability of RFs from high-pass wavelet-filtered CET1-w MR images of primary NPC tumor was poor, particularly when a smaller bin number was used. Exclusive use of high-repeatable RFs is suggested to safeguard model generalizability for wide-spreading clinical utilization.