Front. Oncol.Frontiers in OncologyFront. Oncol.2234-943XFrontiers Media S.A.10.3389/fonc.2024.1401464OncologyOriginal ResearchPrediction of overall survival in patients with locally advanced pancreatic cancer using longitudinal diffusion-weighted MRIBisgaardAnne L. H.12*BrinkCarsten12SchytteTine23BahijRana3Weisz EjlsmarkMathilde23BernchouUffe12BertelsenAnders S.12PfeifferPer23MahmoodFaisal121Laboratory of Radiation Physics, Department of Oncology, Odense University Hospital, Odense, Denmark2Department of Clinical Research, University of Southern Denmark, Odense, Denmark3Department of Oncology, Odense University Hospital, Odense, Denmark
Edited by: Claudio Fiorino, San Raffaele Hospital (IRCCS), Italy
Reviewed by: Garrett Simpson, University of Miami, United States
Luisa Altabella, National Institute of Health (ISS), Italy
Juan Du, Nanjing Drum Tower Hospital, China
*Correspondence: Anne L. H. Bisgaard, anne.bisgaard@rsyd.dk
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background and purpose
Biomarkers for prediction of outcome in patients with pancreatic cancer are wanted in order to personalize the treatment. This study investigated the value of longitudinal diffusion-weighted magnetic resonance imaging (DWI) for prediction of overall survival (OS) in patients with locally advanced pancreatic cancer (LAPC) treated with stereotactic body radiotherapy (SBRT).
Materials and methods
The study included 45 patients with LAPC who received 5 fractions of 10 Gy on a 1.5T MRI-Linac. DWI was acquired prior to irradiation at each fraction. The analysis included baseline values and time-trends of the apparent diffusion coefficient (ADC) and DWI parameters obtained using a decomposition method. A multivariable Cox proportional hazards model for OS was made using best-subset selection, using cross-validation based on Bootstrap.
Results
The median OS from the first day of SBRT was 15.5 months (95% CI: 13.2-20.6), and the median potential follow-up time was 19.8 months. The best-performing multivariable model for OS included two decomposition-based DWI parameters: one baseline and one time-trend parameter. The C-Harrell index describing the model’s discriminating power was 0.754. High baseline ADC values were associated with reduced OS, whereas no association between the ADC time-trend and OS was observed.
Conclusion
Decomposition-based DWI parameters indicated value in the prediction of OS in LAPC. A DWI time-trend parameter was included in the best-performing model, indicating a potential benefit of acquiring longitudinal DWI during the SBRT course. These findings support both baseline and longitudinal DWI as candidate prognostic biomarkers, which may become tools for personalization of the treatment of patients with LAPC.
Pancreatic cancer is the fourth most common cause of cancer-related death in Europe, with a 5-year survival rate of less than 10% and an increasing incidence (1, 2). Twenty percent of the patients are eligible for curative-intent surgery, which increases the survival rate to about 20% (3, 4). Radiotherapy (RT) can be used to downstage tumours, making them eligible for surgery or as a definitive treatment. The development of image-guided RT (IGRT), especially the recent introduction of hybrid MRI linear accelerators (MRI-Linacs), have made hypofractionated stereotactic body RT (SBRT) well-tolerated (5–7). Still, the best RT treatment for patients with pancreatic cancer remains to be settled. A possible problem in interpreting the result from pancreatic cancer studies is the considerable variation within the studied cohorts, possibly related to different responses to RT. Thus, biomarkers may allow a more informed treatment choice (8). The impact of such biomarkers could be either dose escalation or de-escalation to obtain an optimal balance between survival and toxicity.
One promising candidate for such a biomarker is the apparent diffusion coefficient (ADC) derived from diffusion-weighted MRI (DWI). ADC quantifies the motion of water molecules, which indirectly reflects the tissue microstructure (9). Several groups have investigated the value of ADC in detecting pancreatic cancer (10–14) and the relation between ADC and pathological response (15–17). The correlation between overall survival (OS) and ADC has also been investigated (4, 18, 19), with findings indicating that pre-treatment ADC might contain more predictive information than standard clinical parameters such as age and tumour size. However, the predictive power of these studies was limited. A potential improvement can be sought through the MRI-Linac, where DWI can be feasibly integrated into the workflow in each treatment fraction. This setup allows investigation of whether changes in DWI parameters during the treatment course could add independent information to the prediction of the outcome for the patients.
The current study utilized longitudinal DWI data to investigate how baseline DWI parameters and changes in DWI parameters during the treatment course impact survival prediction in patients with locally advanced pancreatic cancer (LAPC). DWI parameters were derived using a data-driven method recently proposed by Rahbek et al. (20, 21). The potential advantage of data-driven approaches is that there are no initial model assumptions which could lead to biased parameters. The standard ADC based on a mono-exponential model was included for comparison.
Materials and methods
A multivariable Cox survival model was made to predict OS in patients with LAPC. The model was based on both clinical parameters and parameters derived from DWI. The DWI parameters were extracted from the GTV volume using both a model-based method (ADC) and a model-free decomposition method. Image analysis was performed by persons without knowledge of the clinical data.
A statistical analysis plan (SAP) was created prior to data analysis of the outcome (Supplementary Materials 3). The SAP states time to local progression as a second endpoint. However, due to very few local progression events (8), stable statistical models could not be derived, and therefore, time to local progression is not included in the analysis.
Patients and endpoints
The study included patients diagnosed with LAPC, treated with SBRT for downstaging or definitive (consolidation) SBRT. The patients had a primary tumour or local recurrence in the pancreas. No tumours were resectable at the time of diagnosis based on evaluations at multi-disciplinary team conferences (three patients underwent surgery after SBRT). No lymph nodes were involved at the time of SBRT for any of the patients. All patients received induction chemotherapy using various chemotherapy regimens for at least two months prior to SBRT and were treated with five fractions of 10 Gy on a 1.5 T MRI-Linac (Unity by Elekta, Stockholm, Sweden). Only patients without visual artefacts in the GTV regions for at least two treatment fractions were included. The follow-up after SBRT, including clinical examination and CT scan, was scheduled for every three months for two years. Only patients that attended at least the first three-month follow-up were included in the study, due to the initial intention of using local progression as a second endpoint. OS was defined as the time from the first SBRT fraction to death for any cause. Time between diagnosis and the first SBRT fraction was included as a clinical parameter to take into account differences due to different chemotherapy regimens. All patients were prospectively included in the MOMENTUM study (clinicaltrials.gov NCT04075305) (22), although this specific investigation was planned after patient inclusion. The research ethics committee at the local institutional board (26/68031) approved inclusion of data in MOMENTUM, and the Region of Southern Denmark (20/35211) approved storage of data. Informed consent was obtained from all patients.
MRI protocols
Before treatment, T2-weighted MRI (T2W) scans were acquired on a 1.5 T MRI simulator (Ingenia, Phillips). Imaging before each treatment fraction at the MRI-Linac included T2W and DWI scans. The majority of the patients (n=39) were DWI scanned using the b-values 0, 30, 80, 150 and 500 s/mm2 (sequence 1), while six patients were scanned using b-values of 0, 20, 60, 100, 300, 800 and 1000 s/mm2 (sequence 2) (acquisition details in Supplementary Tables S1, S2, Supplementary Materials 1).
DWI scans were visually inspected for artefacts. In some patients, an alternation of DWI signal between neighbouring slices was observed, likely due to cross-talk between slices. This was reduced by a pre-processing step in which all slices were convolved with the neighbouring slices with a weight of 0.25 and a weight of 0.5 on the current slice (see Supplementary Figure S6 in Supplementary Materials 1).
Region of interest
GTV delineations from the treatment pre-plan (based on T2W from the MRI simulator) were used to measure the volume of the GTV at baseline. The GTV delineations from the daily adapted treatment plans (based on T2W from the MRI-Linac) were the basis for extracting DWI parameters. The GTVs were transferred to DWIs without translational adjustment. A 5 mm margin was added to the GTVs to account for delineation uncertainty and possible misalignment of the GTV between T2Ws and DWIs which may occur due to motion or geometric distortions present in DWI (23). In the following, all DWI information obtained from the “GTV” refers to the GTV with the 5 mm margin.
ADC calculation
ADC maps were calculated using the standard mono-exponential Stejskal-Tanner model (24), using the b-values 150 and 500 s/mm2 (DWI sequence 1) and 300 and 800 s/mm2 (DWI sequence 2), to adhere as closely as possible to recommendations from the Elekta MRI-Linac Consortium (25). For each patient, a time-trend was extracted using a linear fit of the median ADC within the GTV as a function of fractions. The slope of the fit as well as the ADC value at fraction one were included in the statistical analysis.
Data-driven DWI decomposition
A decomposition method suited for DWI data, the monotonous slope non-negative matrix factorization method (msNMF) recently introduced by Rahbek et al. (20), was used as a data-driven alternative to the model-based ADC analysis. In brief, the DWI signals can be described as a combination of msNMF components, which represent “typical” behaviours of the DWI signals (linear combination with only positive weights). Thus, each DWI voxel (across all b-values) is represented by weights, one for each component.
The components were determined based on data from all patients scanned with DWI sequence 1 (n=39). For this purpose, all GTV voxels from all patients and all fractions were pooled into one data matrix, i.e. only one set of components was derived in total. This set of components was subsequently projected onto each single voxel from all patients (n=45), to obtain the unique weights for each voxel. In the current study, two, three and four components described 86%, 94%, and 96% of the initial data variance, respectively. Based on these values, it was decided to use three components (k=3). Thus, in the following, W1, W2 and W3 refer to the weights associated with the components C1, C2 and C3, respectively. The extracted components are shown in Figure 1 along with example scans. Details regarding image pre-processing and implementation of msNMF in the current study are provided in Supplementary Materials 2.
Example of T2-weighted image (A), DWI images (B–F) and ADC map (G) for a patient. Example of weight maps (W1, W2, W3) corresponding to the three components (C1, C2, C3) (H–J). Components (C1, C2, C3) resulting from decomposition analysis (msNMF) based on a total of 190 DWI scans from 39 patients (5 scans per patient, except for 2 patients with only 2 and 3 scans, respectively) (K).
The spatial distribution of weights associated with the three components were presented as weight maps, (Figures 2A–C). To characterize changes in the weights across fractions, the 10th and 90th percentiles of the weight distributions within the GTV were calculated (Figures 2D–F). Inspired by Rahbek et al. (21), a time-trend of each percentile was extracted using a linear fit (Figure 2G). The slope of the linear fit and the value at fraction one were used as parameters in the statistical analysis. This resulted in 12 decomposition-based parameters [3 components x 2 percentiles x (time-trend + value at fraction one)].
Derivation of decomposition-based DWI parameters for outcome prediction. Examples of weight maps for a patient (A–C). The red contour represents the GTV. The distribution of weights within the GTV (3D volume) are presented using histograms (D–F). The red, vertical lines represent the 10th and 90th percentiles of the distributions. Time trends are extracted from each percentile using a linear fit to the data as a function of fraction number (G).
Statistical analysis
A multivariable Cox proportional hazards survival model was made with OS as the endpoint. Patients who were alive at the cut-off date for data collection (March 30, 2023) were censored. In total, 14 DWI parameters (12 decomposition-based and 2 ADC) and 6 clinical parameters were included in the analysis. The clinical parameters at the time of SBRT were age, GTV volume at baseline, time between diagnosis and SBRT, sex, performance status (PS), and primary tumour/recurrence. The continuous parameters were standardized by subtracting the mean and dividing by the standard deviation.
Parameter selection was based on best-subset selection using bootstrap-based cross-validation. This was chosen as a robust method to limit the risk of overfitting compared to e.g. stepwise selection (26, 27). Thus, for each combination of potential parameters for the multivariable model, a bootstrap of the initial data was made. Patients included in the boot (in-boot patients) were used to calculate the regression constants, and the cross-validated model likelihood was subsequently calculated using the patients not included in the boot (out-of-boot patients). Since the number of out-of-boot patients will vary per boot, the computed likelihood was divided by the number of out-of-boot patients, as suggested by Schemper (28). The entire bootstrap process was iterated 50 times, to obtain the mean cross-validated likelihood. The selected multivariable model was the model which performed best during cross-validation. For the best model, model regression constants (β) and corresponding 2-sided 95% confidence intervals (CI) defined as the range of the 95% most central values were established using 2000 bootstraps.
The multivariable Cox model was validated using a calibration plot, which allows a visual inspection of how well the model fits the data, as well as the model’s ability to discriminate between patients with at short and long survival time (29). Patients were divided into low-, medium- and high-risk groups and a comparison between the Kaplan-Meier survival curve and model estimate was performed within each risk group. The risk groups were defined based on the model’s ranking of the patients’ survival times (i.e. according to the value of the linear predictor for each patient (∑iβixi), such that the high- and low-risk groups contained the 25% shortest and longest living patients, respectively, and the medium-risk group contained the remaining 50%. Besides the calibration plot, the model’s ability to discriminate between patients with high and low risk was reported using the C-Harrell index, which measures the fraction of patient pairs in which the one with the lower risk survives the longest.
The median survival time for the entire cohort was assessed using the Kaplan-Meier estimator, and the median potential follow-up time was calculated using the reverse Kaplan-Meier method (30). Univariable Cox models were made to provide an overview of the entire data set.
Results
In total, 50 patients fulfilled the inclusion criteria (see section “Patients and endpoints”). Out of these, five patients were excluded before statistical analysis, as their DWI sequences differed substantially from those of the remaining patients (sequence 1 and 2). For two patients, DWIs were missing for two and three fractions, respectively, however, they were kept as part of the analysis. The analysis was thus based on 45 patients (demographics shown in Table 1).
Patient demographics for 45 patients receiving SBRT for localized pancreatic cancer.
Overall(N=45)
Age (years)
Mean (SD)
67.5 (10.1)
Median [Min, Max]
69.5 [45.2, 85.3]
Sex
F
26 (57.8%)
M
19 (42.2%)
Performance status
0
17 (37.8%)
1+
28 (62.2%)
Clinical presentation
Primary tumour (LAPC)
37 (82.2%)
Local recurrence
8 (17.8%)
Time between diagnosis and RT (months)
Mean (SD)
7.40 (5.89)
Median [Min, Max]
6.48 [1.15, 35.4]
GTV volume (cm3)
Mean (SD)
21.6 (15.8)
Median [Min, Max]
19.1 [2.68, 77.0]
T-stage
T3-4
45 (100%)
Chemotherapy
Yes
45 (100%)
The reported T-stage is at the time of diagnosis. No re-staging was performed at the time of recurrence. The reported GTV volumes were obtained from the treatment pre-plan, without the 5 mm expansion.
The decomposition analysis (msNMF) resulted in three well-separated components: one slowly decaying component (C1), one fast decaying component (C2), and one component mainly driven by b=0 s/mm2 (C3) (Figure 1K). The weight maps revealed tissue heterogeneity, highlighting different regions within the GTV (Figures 1H–J). The weight maps did not correlate 1:1 with the DWI images or the ADC map and thus provide complementary information (Figures 1B–G).
The median OS was 15.5 months (95% CI: 13.2-20.6) (see Kaplan-Meier plot in Supplementary Figure S1 in Supplementary Materials 1). The median potential follow-up time was 19.8 months. Based on univariable analyses, none of the clinical parameters showed a statistically significant association with OS on a 5% level; however, the association between age and OS was borderline significant (Figure 3). Several of the DWI-derived parameters, including ADC at the first treatment fraction, showed a statistically significant association with OS.
Hazard ratios and 95% confidence intervals for univariable Cox proportional hazards models for overall survival. In the graphical representation, the confidence interval for the DWI parameter “Component3_prc10_frac1” was cut off at 0.10 in order to make a more well-balanced figure. The DWI parameter names refer to the msNMF components, the percentile (10th or 90th percentile), and whether the value represents the time-trend or the value at fraction one (see section “Data-driven DWI decomposition”).
The best-performing cross-validated multivariable Cox model included two decomposition-based parameters (“Component1_prc10_time-trend” and “Component2_prc90_frac1”) (Figure 4). The C-Harrell conformance index for the best model was 0.754. A calibration plot demonstrated good agreement between the model and data for all risk groups (Figure 5). The individual performance of the two predictors in the best-performing model is shown in Supplementary Figure S2 in Supplementary Materials 1. The model’s performance is notably reduced when the parameters are used separately and have C-Harrell indexes of 0.633 and 0.688 for Component2_prc90_frac1 and Component1_prc10_time-trend, respectively.
Hazard ratios and 95% confidence intervals for the multivariable Cox proportional hazards model for overall survival. The parameters included in this model were selected by the cross-validation process.
Comparison of the multivariable Cox model and the Kaplan-Meier estimator for the best-performing model for overall survival. The model included the decomposition-based parameters “Component1_prc10_time-trend” and Component2_prc90_frac1” (see section “Data-driven DWI decomposition”). Patients were split into high, medium, and low-risk groups based on the 25% and 75% percentiles of the calculated linear predictors (-0.66 and 0.72), i.e. the high- and low-risk groups each contained 25% of the patients, and the medium-risk group contained 50% of the patients. The confidence intervals for each of the risk groups were overlapping (not shown). Only patients who attended the first three months follow up were included in the study, hence the graphs show a constant survival probability of one during the first three months.
Discussion
This study investigated the value of longitudinal DWI in the prediction of OS in patients with LAPC utilizing parameters derived using both a standard, model based approach (ADC) and a model-free decomposition approach (msNMF). To our knowledge, this study is the first to utilize longitudinal DWI for this purpose.
The best-performing model for OS prediction included only two parameters; one time-trend DWI parameter (Component1_prc10_time-trend) and one baseline DWI parameter (Component2_prc90_frac1), both of which were based on DWI decomposition. None of the clinical parameters were selected by the cross-validation process. The best model reached a C-Harrell index of 0.754 indicating that the model is good at determining which of two patients will survive the longest (a value of 0.5 corresponds to a random guess, and a value of 1 means that the model can perfectly rank the patients’ survival times). The C-Harrell index was markedly reduced if the model was based on just one of the parameters compared to the best model (Supplementary Figure S2 in Supplementary Materials 1), indicating that both baseline information (Component2_prc90_frac1) and longitudinal changes (Component1_prc10_time-trend) might be important for the prediction of OS.
Previous studies have indicated that image-based parameters derived from pre-treatment CT, PET and MR imaging can provide prognostic information in pancreatic cancer (4, 31, 32). With the ability to perform daily imaging in a non-invasive and radiation-free manner, MRI-Linacs represent an interesting opportunity for exploring changes in image-based parameters over the course of treatment using longitudinal imaging. The feasibility of extracting such changes with an MRI-Linac have been demonstrated, with results indicating an increased performance compared to parameters derived at a single time-point (33–35). The current study made use of this opportunity to investigate image-based parameters derived from longitudinal DWI, exploring different methods for extraction of parameters. The msNMF method has previously been tested in patients with brain metastasis, showing borderline significant differences between responders and non-responders (21), but has not been tested as a predictor of OS in pancreatic cancer. It was included in the current study to investigate whether it could provide information useful for the prediction of OS in pancreatic cancer, as a proof of concept. Interestingly, the cross-validation process selected decomposition-based parameters instead of ADC information. A possible explanation could be that the decomposition-based parameters might be more stable than the standard ADC values. Increased stability may be achieved as the extracted components are based on data from all patients and fractions, which may reduce the impact of noise. The C-Harrell index for a model based on the baseline ADC (ADC_frac1) and the ADC time-trend (ADC_time-trend) was 0.623 (data not shown), indicating an inferior discriminating power. That being said, the Pearson correlation between ADC_frac1 and Component2_prc90_frac1 is 0.78, showing that most of the information in Component2_prc90_frac1 is also present in the ADC signal (Supplementary Figure S4 in Supplementary Materials 1). Thus, based on the current results, the main difference between the standard ADC approach and the msNMF approach is within the longitudinal changes in which msMNF parameters showed a stronger association with OS compared to ADC_time-trend. The potential added benefit of the msNMF approach compared to the ADC approach should be weighed against the increased complexity of the image post-processing, which may make msNMF more labour intensive.
The multivariable Cox model shown in Figure 4 suggests that high values of Component2_prc90_frac1 are associated with reduced expected lifetime (hazard ratio (HR)>1). According to the univariable Cox models presented in Figure 3, the same applies for high values of ADC_frac1 (HR>1), which is in line with the observed correlation between Component2_prc90_frac1 and ADC_frac1 (Supplementary Figure S4 in Supplementary Materials 1). Initially, these results seem to contradict the findings in some other studies, in which an association between low baseline ADC (comparable to ADC_frac1 in the current study) and poor survival was reported (18, 19, 36). The explanation might be a difference in the definition of the region of interest from which the DWI parameters were derived. The current study used the clinically available GTV so that no additional delineation was needed, whereas previous studies focused on the delineation of the “viable” part of the tumour, which differs from the GTV by excluding necrotic and cystic parts. This explanation is supported by a study by Lyng et al. (37), which shows that some necrotic regions are related to increased ADC values. Furthermore, necrosis has been related to less favourable outcomes in patients with LAPC (38), which matches high values of ADC_frac1 in patients with low OS seen in the current study. That ADC’s capacity as a biomarker of response is affected by the choice of regions of interest has also been demonstrated in other anatomical regions (39).
For Component 2 frac1 and Component 1 time-trend, the 10% and 90% percentiles contain quite similar information, with a correlation of 0.70 and 0.61, respectively (Supplementary Figure S4 in Supplementary Materials 1). Since the average of the 10% and 90% percentile values for symmetrical distributions is very close to the median, it is likely that the median values could have been used instead. In fact, a model based on the related median values of the best model (i.e. a model based on “Component1_prc50_time-trend” and Component2_prc50_frac1”) had a C-Harrell index of 0.757 (Supplementary Figure S5 in Supplementary Materials 1), similar to the C-Harrell index of the best model (0.754). For future studies, it is likely as good (and simpler) to focus on the median values instead of the tails of the weight distributions.
To take into account the heterogeneity of the patient cohort, clinical parameters such as primary tumour/local recurrence were included in the analysis. Interestingly, the cross-validation process did not select any of the clinical parameters, indicating a superiority of baseline and longitudinal DWI to predict OS compared to standard clinical parameters.
It should be noted that the two parameters in the best model showed a low correlation with the clinical parameters. The largest Pearson correlation between the model parameters and the continuous clinical parameters was 0.24, and the lowest p-value of a Mann–Whitney U test of best model parameters between the levels of the categorical clinical parameters was 0.53. Therefore, there is no indication that the DWI parameters could be a proxy for any of the clinical parameters.
In general, OS may depend on many factors, both related to the tumour and to the individual patient’s response to treatments, including SBRT. The observed association between DWI parameters and OS could potentially be related to differences in the tumour microstructure between patients, which might be captured by DWI. While ADC is commonly interpreted in terms of cell density, decomposition-based parameters could potentially reflect the composition of the tumour microstructure, as sub-compartments within each image voxel may contain different tissue types. A possible explanation for the observed association between DWI parameters and OS could be that the tumour microstructure is indicative of the tumour aggressiveness and tendency to form distant metastasis, something that is likely linked to OS (18). Moreover, tumours with different microstructure may respond differently to treatments. Based on the results, it might be relevant to investigate the predictive value of DWI, in order to identify patients who responds well/poorly to SBRT.
Initially, time to local progression was included as a second endpoint in the statistical analysis plan. However, due to a small number of local progression events (n=8), it was not possible to derive stable statistical models, and thus, no results are presented for this endpoint. Hence, it might be a relevant topic for a future study to investigate if DWI could provide information of the local response to a local treatment (SBRT), as this information might be relevant in order to personalize the treatment.
A careful parameter selection based on cross-validation was performed to limit the effect of overfitting, instead of e.g stepwise parameter selection (26, 27). However, it is still considered a major limitation of the current study that the statistical model was based on 45 patients only. Although it is considered a relatively large cohort compared to other RT studies within LAPC, it may be a poor representation of patients with LAPC treated worldwide. Validation of the statistical model in patients from other centres is thus warranted.
In conclusion, a statistical model to predict survival after SBRT in patients with LAPC based on only two DWI parameters has been developed. The model contains both baseline information and DWI changes during the SBRT course and showed ability to discriminate between patients with short and long survival times. It is the hope that in the future, prognostic information from DWI can assist in stratifying patients for individual treatment (e.g. dose escalation), to improve treatment outcomes.
Data availability statement
The data analysed in this study is subject to the following licenses/restrictions: The data are not publicly available. Informed written consent was obtained from all patients to include the data in the Momentum study (a database where data is shared between hospitals in the Elekta MRI-linac Consortium). Requests to access these datasets should be directed to The Momentum data management task force (https://mrlconsortium.org/introduction-to-momentum/).
Ethics statement
The studies involving humans were approved by Research Ethics Committee at the University of Southern Denmark. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by the Danish Cancer Society (Grant no. R231-356 A13852), Danish Comprehensive Cancer Center RT (Danish Cancer Society grant) (Grant no. R191-A11526), and by MANTRA (New MAgNetic resonance Technology for Response Adapted radiotherapy), a Frontline research centre based at Odense University Hospital, Denmark. The funding sources had no role in designing or conducting the study.
Acknowledgments
We thank Sofie Rahbek for providing Matlab code for the msNMF decomposition method and for the assistance with implementation of the prediction framework.
Conflict of interest
TS acknowledges payment from ELEKTA for lectures and travel assistance from ELEKTA for meetings, and is a member of the ELEKTA MRI-Linac consortium steering committee. ABe and UB acknowledges travel assistance from ELEKTA for meetings.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be constructed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2024.1401464/full#supplementary-material
ReferencesKirkegårdJBojesenABNielsenMFMortensenFV. Trends in pancreatic cancer incidence, characteristics, and outcomes in Denmark 1980–2019: A nationwide cohort study. Cancer Epidemiol. (2022) 80:102230. doi: 10.1016/j.canep.2022.102230SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalA. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660MizrahiJDSuranaRValleJWShroffRT. Pancreatic cancer. Lancet. (2020) 395:2008–20. doi: 10.1016/S0140-6736(20)30974-0NodaYTomitaHIshiharaTTsuboiYKawaiNKawaguchiM. Prediction of overall survival in patients with pancreatic ductal adenocarcinoma: histogram analysis of ADC value and correlation with pathological intratumoral necrosis. BMC Med Imaging. (2022) 22:1–10. doi: 10.1186/s12880-022-00751-3RaaymakersBWJürgenliemk-SchulzIMBolGHGlitznerMKotteANTJVan AsselenB. First patients treated with a 1.5 T MRI-Linac: Clinical proof of concept of a high-precision, high-field MRI guided radiotherapy treatment. Phys Med Biol. (2017) 62:L41–50. doi: 10.1088/1361-6560/aa9517HalWAStrazaMWChenXMickeviciusNEricksonBSchultzC. Initial clinical experience of Stereotactic Body Radiation Therapy (SBRT) for liver metastases, primary liver Malignancy, and pancreatic cancer with 4D-MRI based online adaptation and real-time MRI monitoring using a 1.5 Tesla MR-Linac. PloS One. (2020) 15:1–10. doi: 10.1371/journal.pone.0236570EjlsmarkMWSchytteTBernchouUBahijRWeberBMortensenMB. Radiotherapy for locally advanced pancreatic adenocarcinoma-A critical review of randomised trials. Curr Oncol. (2023) 30:6820–37. doi: 10.3390/curroncol30070499LeNSundMVinciA. Prognostic and predictive markers in pancreatic adenocarcinoma. Dig Liver Dis. (2016) 48:223–30. doi: 10.1016/j.dld.2015.11.001KohDMCollinsDJ. Diffusion-weighted MRI in the body: Applications and challenges in oncology. Am J Roentgenol. (2007) 188:1622–35. doi: 10.2214/AJR.06.1403IchikawaTErturkSMMotosugiUSouHIinoHArakiT. High-b value diffusion-weighted MRI for detecting pancreatic adenocarcinoma: Preliminary results. Am J Roentgenol. (2007) 188:409–14. doi: 10.2214/AJR.05.1918MatsukiMInadaYNakaiGTatsugamiFTanikakeMNarabayashiI. Diffusion-weighed MR imaging of pancreatic carcinoma. Abdom Imaging. (2007) 32:481–3. doi: 10.1007/s00261-007-9192-6InanNArslanAAkanselGAnikYDemirciA. Diffusion-weighted imaging in the differential diagnosis of cystic lesions of the pancreas. Am J Roentgenol. (2008) 191:1115–21. doi: 10.2214/AJR.07.3754SeungSLJaeHBBeomJPSeongHPKimNParkB. Quantitative analysis of diffusion-weighted magnetic resonance imaging of the pancreas: Usefulness in characterizing solid pancreatic masses. J Magn Reson Imaging. (2008) 28:928–36. doi: 10.1002/jmri.21508KimBLeeSSSungYSCheongHByunJHKimHJ. Intravoxel incoherent motion diffusion-weighted imaging of the pancreas: Characterization of benign and Malignant pancreatic pathologies. J Magn Reson Imaging. (2017) 45:260–9. doi: 10.1002/jmri.25334NiwaTUenoMOhkawaSYoshidaTDoiuchiTItoK. Advanced pancreatic cancer: The use of the apparent diffusion coefficient to predict response to chemotherapy. Br J Radiol. (2009) 82:28–34. doi: 10.1259/bjr/43911400DalahEEricksonBOshimaKSchottDHallWAPaulsonE. Correlation of ADC with pathological treatment response for radiation therapy of pancreatic cancer. Transl Oncol. (2018) 11:391–8. doi: 10.1016/j.tranon.2018.01.018CuneoKCChenevertTLBen-JosefEFengMUGreensonJKHussainHK. A pilot study of diffusion- weighted mri in patients undergoing neoadjuvant chemoradiation for pancreatic cancer. Transl Oncol. (2014) 7:644–9. doi: 10.1016/j.tranon.2014.07.005Garces-DescovichAMorrisonTCBekerKJaramillo-CardosoAMoserAJMorteleKJ. DWI of pancreatic ductal adenocarcinoma: A pilot study to estimate the correlation with metastatic disease potential and overall survival. Am J Roentgenol. (2019) 212:323–31. doi: 10.2214/AJR.18.20017KurosawaJTawadaKMikataRIshiharaTTsuyuguchiTSaitoM. Prognostic relevance of apparent diffusion coefficient obtained by diffusion-weighted MRI in pancreatic cancer. J Magn Reson Imaging. (2015) 42:1532–7. doi: 10.1002/jmri.24939RahbekSMadsenKHLundellHMahmoodFHansonLG. Data-driven separation of MRI signal components for tissue characterization. J Magn Reson. (2021) 333:107103. doi: 10.1016/j.jmr.2021.107103RahbekSMahmoodFTomaszewskiMHansonLMadsenK. Decomposition-based framework for tumor classification and prediction of treatment response from longitudinal MRI. Phys Med Biol. (2023) 68:025006. doi: 10.1088/1361-6560/acaa85de Mol van OtterlooSRChristodouleasJPBlezerELAAkhiatHBrownKChoudhuryA. The MOMENTUM study: an international registry for the evidence-based introduction of MR-guided adaptive therapy. Front Oncol. (2020) 10. doi: 10.3389/fonc.2020.01328HaslerSWKallehaugeJFHansenRHSamsøeEArpDTNissenHD. Geometric distortions in clinical MRI sequences for radiotherapy: insights gained from a multicenter investigation. Acta Oncol (Madr). (2023) 62:1551–60. doi: 10.1080/0284186X.2023.2266560StejskalEOTannerJE. Spin diffusion measurements: Spin echoes in the presence of a time-dependent field gradient. J Chem Phys. (1965) 42:288–92. doi: 10.1063/1.1695690KooremanESvan HoudtPJKeesmanRPosFJvan PeltVWJNoweeME. ADC measurements on the Unity MR-linac – A recommendation on behalf of the Elekta Unity MR-linac consortium. Radiother Oncol. (2020) 153:106–13. doi: 10.1016/j.radonc.2020.09.046HarrellFEJ. Regression modelling strategies. Springer New York, NY: Springer (2015).BrinkCHansenCRFieldMPriceGThwaitesDSarupN. Distributed learning optimisation of Cox models can leak patient data: Risks and solutions. Arxiv (2022). doi: 10.48550/arXiv.2204.05856SchemperM. Further results on the explained variation in proportional hazards regression. Biometrica. (1992) 79:202–4. doi: 10.1093/biomet/79.1.202RoystonPAltmanDGA. External validation of a Cox prognostic model: principles and methods. BMC Res Methodol. (2013) 13. doi: 10.1186/1471-2288-13-33SchemperMSmithTL. A note on quantifying follow-up in studies of failure time. Control Clin Trials. (1996) 17:343–6. doi: 10.1016/0197-2456(96)00075-XMoriMPassoniPIncertiEBettinardiVBroggiSReniM. Training and validation of a robust PET radiomic-based index to predict distant-relapse-free-survival after radio-chemotherapy for locally advanced pancreatic cancer. Radiother Oncol. (2020) 153:258–64. doi: 10.1016/j.radonc.2020.07.003PalumboDMoriMPratoFCrippaSBelfioriGReniM. Prediction of early distant recurrence in upfront resectable pancreatic adenocarcinoma: A multidisciplinary, machine learning-based approach. Cancers (Basel). (2021) 13. doi: 10.3390/cancers13194938BoldriniLCusumanoDChiloiroGCasàCMasciocchiCLenkowiczJ. Delta radiomics for rectal cancer response prediction with hybrid 0.35 T magnetic resonance-guided radiotherapy (MRgRT): a hypothesis-generating study for an innovative personalized medicine approach. Radiol Med. (2019) 124:145–53. doi: 10.1007/s11547-018-0951-ySimpsonGSpielerBDoganNPortelanceLMellonEAKwonD. Predictive value of 0.35 T magnetic resonance imaging radiomic features in stereotactic ablative body radiotherapy of pancreatic cancer: A pilot study. Med Phys. (2020) 47:3682–90. doi: 10.1002/mp.14200CusumanoDBoldriniLYadavPYuGMusurunuBChiloiroG. Delta radiomics for rectal cancer response prediction using low field magnetic resonance guided radiotherapy: an external validation. Phys Medica. (2021) 84:186–91. doi: 10.1016/j.ejmp.2021.03.038De RobertisRMartiniPTDemozziECorsoFDBassiCPederzoliP. Diffusion-weighted imaging of pancreatic cancer. World J Rdiology. (2015) 7:319–29. doi: 10.4329/wjr.v7.i10.319LyngHHaraldsethORofstadEK. Measurement of cell density and necrotic fraction in human melanoma xenografts by diffusion weighted magnetic resonance imaging. Magn Reson Med. (2000) 836:828–36. doi: 10.1002/(ISSN)1522-2594HiraokaNInoYSekineSTsudaHShimadaKKosugeT. Tumour necrosis is a postoperative prognostic marker for pancreatic cancer patients with a high interobserver reproducibility in histological evaluation. Br J Cancer. (2010) 103:1057–65. doi: 10.1038/sj.bjc.6605854MahmoodFJohannesenHHGeertsenPOpheimGFHansenRH. The effect of region of interest strategies on apparent diffusion coefficient assessment in patients treated with palliative radiation therapy to brain metastases. Acta Oncol (Madr). (2015) 54:1529–34. doi: 10.3109/0284186X.2015.1061211