Skip to main content

Longitudinal assessment of diffusion-weighted imaging during magnetic resonance-guided radiotherapy in head and neck cancer

Abstract

Background

For radiotherapy of head and neck cancer (HNC) magnetic resonance imaging (MRI) plays a pivotal role due to its high soft tissue contrast. Moreover, it offers the potential to acquire functional information through diffusion weighted imaging (DWI) with the potential to personalize treatment. The aim of this study was to acquire repetitive DWI during the course of online adaptive radiotherapy on an 1.5 T MR-linear accelerator (MR-Linac) for HNC patients and to investigate temporal changes of apparent diffusion coefficient (ADC) values of the tumor and subvolume levels.

Methods

27 patients treated with curative RT on the 1.5 T MR-Linac with at least weekly DWI in treatment position were included into this prospective analysis and divided in four risk groups (HPV-status and localisation). Tumor and lymph node volumes (GTV-P/GTV-N) were delineated on b = 500 s/mm2 images while ADC maps were calculated using b = 150/200 and 500 s/mm2 images. Absolute and relative temporal changes of mean ADC values, tumor volumes and a high-risk subvolume (HRS) defined by low ADC tumor voxels (600 < ADC < 900 × 10−6 mm2/s) were analyzed. Relative changes of mean ADC values, tumor volumes and HRS were statistically tested using Wilcoxon-signed-rank test.

Results

Median pretreatment ADC value for all patients resulted in 1167 × 10−6 mm2/s for GTV-P and 1002 × 10−6 mm2/s for GTV-N while absolute pretreatment tumor volume yielded 9.1 cm3 for GTV-P and 6.0 cm3 for GTV-N, respectively. Pretreatment HRS volumes were 1.5 cm3 for GTV-P and 1.3 cm3 for GTV-P and GTV-N. Median ADC values increase during 35 fractions of RT was 49% for GTV-P and 24% for GTV-N during RT. Median tumor volume decrease was 68% and 52% for GTV-P and GTV-N with a median HRS decrease of 93% and 87%. Significant differences from 0 for mean ADC were observed starting from week 1, for tumor volumes from week 2 for GTV-P and week 1 for GTV-N and for HRS in week 1 for GTV-P and week 2 for GTV-N.

Conclusion

Longitudinal DWI acquisition in HNC is feasible on a MR-Linac during the course of online adaptive MR-guided radiotherapy. Changes in ADC and volumes can be assessed, but future work needs to explore the potential for biologically guided treatment individualization.

Trial registration: NCT04172753, actual study start: 09.05.2018.

Background

Radiotherapy (RT) is one of the cornerstones in the treatment of head and neck cancer (HNC). Different prognostic factors exist for HNC, with the most important being human papilloma virus 16 (HPV) and smoking status in oropharyngeal cancers (OPC) and mainly tumor stage for other tumor sites [1,2,3]. In depth research has been done on imaging biomarkers in HNC, with partially contradictory results. Positron emission tomography either with FDG or hypoxia specific tracers were described to offer prognostic potential and have been used to escalate or de-escalate treatment by dose painting [4,5,6,7,8].

Due to its superior soft tissue contrast magnetic resonance imaging (MRI) plays a pivotal role in target volume definition and might prove beneficial for automation of these time comsuming steps in HNC [9,10,11]. Besides the anatomical information from MRI, it offers the possibility to acquire functional imaging such as dynamic contrast enhanced (DCE) and diffusion-weighted imaging (DWI). Several single center studies looked at different imaging biomarkers with regards to prognostic information and technical validation with contradicting results, especially in HPV associated OPC [12,13,14,15,16,17]. Lambrecht et al. showed prognostic information of DWI acquired before radiotherapy with regards to local control [18]. Contradictory to these results, there are observations in HPV associated OPC, that DWI is not prognostic and may not be compared to non-HPV associated HNC [15, 19].

Due to the complexity and limited availability of MRI, only a few studies were reported focusing on longitudinal imaging during the course of fractionated RT [17, 20,21,22]. In recent years the technique of magnetic resonance linear accelerators (MR-Linac) has been introduced and first data were published for MR-guided RT for HNC [23, 24]. Initial results showed the feasibility of MR-guided RT for HNC, providing the potential to acquire longitudinal anatomical imaging with high resolution and functional imaging during the course of fractionated radiotherapy [23,24,25]. Moreover, in our department a potential quantitative imaging biomarker (QIB) based on a cluster of ADC values has been established in a preclinical model, suggesting a high-risk subvolume (HRS), which could be used for individualized dose prescription [26]. A validation on a clinical cohort of patients has been performed in which the retrained HRS was found to be significantly associated to outcome after primary radiochemotherapy, so this ADC-based HRS might serve as a biomarker for stratification or therapy individualization of HNC patients [27].

The aim of the present study was to explore the feasibility of online acquired, longitudinal DWI with an exploratory analysis of changes in different DWI parameters during the course of MR-guided RT in HNC patients.

Methods

Patient and treatment characteristics

A total of n = 28 patients with HNC, which were treated on the 1.5 T MR-Linac (Unity, Elekta AB, Sweden) between October 2018 and December 2021 have been prospectively included into this analysis within a enclosed project funded by the German Research Council. All patients gave written informed consent to be treated within this prospective trial of MR-guided adaptive RT which was approved by the local ethics committee (no. 659/2017BO1, NCT04172753). All patients were immobilized with a 5-point radiotherapy mask on a dedicated head step for the use on the 1.5 T MR-Linac in a neutral neck position [13]. The planning-CT of 2 to 3 mm slice thickness and all MRIs have been acquired in this setup as shown in detail elsewhere [13, 14, 23, 24]. Patients were instructed to limit swallowing whenever possible. Treatment consisted of fractionated RT with 60 Gy and 54 Gy in 30 fractions to the macroscopic tumor, the high risk area and the elective volume with a sequential boost to the macroscopic tumor of 10 Gy in 5 fractions, consistent with international guidelines. Concomitant weekly Cisplatin was administered. The treatment was performed on the MR-Linac with an online “adapt-to-position” workflow and offline adaptation in case of large anatomical deviations based on the treating physicians discretion [24]. Detailed patient characteristics are shown in Table 1. Follow-up was according to clinical routine with a FDG-PET/CT after 3 months in case of nodal involvement and further clinical and radiological examinations every 3 months. For OPC, risk groups have been defined according to Ang et al. with stratification based on p16-, smoking-status and TNM stage [1]. The remaining six non-OPC patients have been included into a fourth patient group.

Table 1 Patient characteristics

Imaging protocol

Imaging of all patients was performed on the 1.5 T MR-Linac except for pretreatment images of patient 25 which were acquired on a 1.5 T diagnostic scanner (Ingenia, Philips). The protocol consisted of a single-shot echo-planar imaging (SS-EPI) DWI as well as T1- and/or T2-weighted anatomical imaging. For the first eleven patients the acquired b-values of the DWI sequence were 0, 200, 500 and 800 s/mm2 while starting with patient 12, the DWI sequence was adapted according to the recommendations of the MR-Linac consortium [28] and b-values 0, 150 and 500 s/mm2 were acquired. For patient 2 the b-value 500 s/mm2 was not part of the DWI sequence pretreatment and in fraction four of RT, while patient one had to be excluded from analysis due to missing DWI data. Details about the sequence parameters of both DWI as well as the two anatomical MRI protocols are provided in Table 2. Patients were imaged in RT position either directly before the start, during beam on or after completion of RT on the MR-Linac. The imaging protocol was applied during simulation before the start of radiotherapy and sequentially during RT, approximately once per week.

Table 2 Details of the sequence parameters

Image processing

Volumes of interest (VOIs) defined as primary tumors (GTV-P) and all conspicuous lymph nodes (GTV-N) were delineated by a board-certified radiation oncologist (SB) on the b-500 images with the open-source software 3D Slicer (Version 4.10) taking the anatomical images as visual input to ensure the capture of the whole tumor volume. Then, ADC maps were calculated using an in-house written python script (version 3.8.10) and b-values 150/200 and 500 s/mm.2 with the mono-exponential model [29]

$${\text{SI}} = {\text{SI}}_{0} \cdot {\text{e}}^{{ - {\text{b}} \cdot {\text{ADC}}}} .$$

Here SI is defined as the signal intensity and SI0 as the calculated signal intensity at b-value 0 s/mm2.

Furthermore, a high-risk subvolume (HRS) was defined, representing potentially radioresistant tumor subregions [26] inside both GTV-P and GTV-N, separately. This HRS was defined as all voxels inside the GTV-P or GTV-N with ADC values in the range of 600–900 × 10−6 mm2/s. ADC thresholds were used from a prior investigation [27] and adapted for the use on the MR-Linac to the known underestimation of ADC on the MR-Linac [13, 30, 31].

Statistics

GTV-P and GTV-N volumes as well as mean ADC values (ADCmean) and HRS volume were calculated using Matlab 2020a (MathWorks, Natick, MA, USA). Target volumes smaller than 1 cm3 and HRS smaller than 0.2 cm3 were excluded from further analysis. Absolute values of ADC and volumes were evaluated for all patients for baseline and every week of radiotherapy. Furthermore, linear regression as a function of treatment time regarding mean ADC values as well as volumes of GTV-P, GVT-N and HRS was fitted for all risk groups with Matlab. Additionally, relative changes of ADC and volume from baseline to every week of RT were evaluated for every risk group. One-sample Wilcoxon-signed-rank test with a significance level of 5% was used to evaluate relative changes of ADCmean, absolute tumor volumes and HRS volumes in every week of radiotherapy for a significant difference from 0 using SPSS statistical package 28.0.0.0 (SPSS Inc, Chicago, Illinois).

Results

A total of 28 patients have been treated on the MR-Linac out of which 27 primary lesions as well as 39 lymph nodes from 27 patients were available for analysis. Median (interquartile range) number of scans per patient during treatment was 7 (1). Exemplary pretreatment anatomical as well as DW images for patient 17 including delineations of the primary tumor and calculated HRS are shown in Fig. 1. At the time of analysis median follow up was 40.4 months, with five patients having died due to metastatic disease and one likely due to pneumonia after aspiration. One patient with a combined hypopharyngeal and proximal esophageal carcinoma (in initial endoscopy distant to each other) showed a recurrence in the proximal esophagus 18 months after completion of the treatment with no recurrence in the hypopharynx. No isolated local or regional recurrence have been observed in the cohort.

Fig. 1
figure 1

Exemplary pretreatment images of patient #17 including T1w (A), T2w (B), b0 image (C), b150 image (D), b500 image (E) and calculated ADC map (F) with GTV depicted in red and high-risk subvolume (600 < ADC < 900 × 10−6 mm2/s) in blue

Median (interquartile range) pretreatment ADC value for all patients amounted to 1167 (282) × 10−6 mm2/s for GTV-P and 1002 (398) × 10−6 mm2/s for GTV-N, respectively. Over the course of radiotherapy, the ADCmean of primary lesions increased to 1687 (381) × 10−6 mm2/s and 1420 (424) × 10−6 mm2/s of conspicuous lymph nodes in the last week of radiotherapy. The temporal changes including linear regression of ADCmean for GTV-P and GTV-N over the course of radiotherapy depending on the patient’s risk classification are displayed in Figs. 2 and 3 for all patients.

Fig. 2
figure 2

Mean ADC values for GTV-P over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients (according to Ang et al. [1]. Solid lines represent linear regression for each group of patients

Fig. 3
figure 3

Mean ADC values for GTV-N over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients (according to Ang et al. [1]. Solid lines represent linear regression for each group of patients

Regarding tumor volumes, median pretreatment GTV-P volume was 9.1 (15.4) cm3 while GTV-N displayed a median size of 6.0 (12.2) cm3. Median primary tumor volumes decreased over the course of radiotherapy to 3.6 (3.1) cm3 in the last week of RT. In contrast, the conspicuous lymph nodes had a median volume of 2.8 (4.4) cm3 in the last week of RT. For all patients, trends of absolute tumor volumes and linear regression are visualized in Figs. 4 and 5 for GTV-P and GTV-N, respectively.

Fig. 4
figure 4

Volumes of GTV-P over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients. Solid lines represent linear regression for each group of patients

Fig. 5
figure 5

Volumes of GTV-N over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients. Solid lines represent linear regression for each group of patients

The ADC-based HRS of the GTV-P and GTV-N showed a median pretreatment size of 1.5 (3.1) cm3 and 1.3 (2.4) cm3, respectively. In the last week of RT, the median HRS had shrunk to 0.4 (0.0) cm3 for GTV-P and 0.6 (0.5) cm3 for GTV-N. Figures 6 and 7 show the trend of HRS for all analyzed patients and the linear regression. Detailed information about ADCmean, absolute tumor volume and HRS volumes for every risk group pretreatment and in the last week of RT are presented in Table 3. Additionally, absolute values for ADCmean, tumor volume and HRS are summarized on a weekly basis for all patients in Fig. 8 and separated by risk groups in Fig. 9.

Fig. 6
figure 6

Trend for HRSs inside GTV-P over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients. Solid lines represent linear regression for each group of patients

Fig. 7
figure 7

Trend for HRSs inside each individual GTV-N over the course of radiotherapy for high (red), intermediate (gold), low risk (blue) and non-oropharyngeal (green) patients. Solid lines represent linear regression for each group of patients

Table 3 ADCmean, absolute tumor volumes and HRS volumes before the start of radiotherapy (RT), in week 7 of RT and the relative change over course of RT separated by risk groups
Fig. 8
figure 8

Boxplots of mean ADC, absolute tumor volumes and absolute HRS volumes for all patients before and in every week of radiotherapy

Fig. 9
figure 9

Boxplots of mean ADC values (top) and absolute tumor volumes (bottom) pretreatment (w0) and during radiotherapy (w1–w7) for high risk (red), intermediate risk (gold), low risk (blue) and non-oropharyngeal (green) patients

Median ADCmean increase to the last week of radiotherapy was 49 (52)% for GTV-P and 24 (48)% and GTV-N. GTV-P volumes showed a median decrease of 68 (24)% while GTV-N volumes decreased by 52 (62)% to the last week of RT. Furthermore, median decrease of HRS volumes yielded 93 (9)% for GTV-P and 87 (18)% for GTV-N. Significant differences (p < 0.05) for relative change of ADCmean of GTV-P and GTV-N were found starting in week 1 of RT while tumor volumes showed significant differences (p < 0.05) beginning in week 2 for GTV-P and week 1 for GTV-N, respectively. HRS size displayed significant differences (p < 0.05) between baseline and week 1 of RT for GTV-P and week 2 for GTV-N. Furthermore, we evaluated the relative changes of ADCmean, absolute tumor volumes and HRS size depending on patient risk stratification from baseline to every week of radiotherapy (cf. Figure 10). All clinical risk categories showed an increase of ADC parameters during treatment and decrease in HRS during treatment.

Fig. 10
figure 10

Boxplots of relative change of mean ADC values (top) and tumor volumes (bottom) from baseline to every week of radiotherapy (w1–w7) for high risk (red), intermediate risk (gold), low risk (blue) and non-oropharyngeal (green) patients. Outliers are shown with a plus sign. Significant different changes (p < 0.05) from 0 are indicated by a star

Discussion

As new QIBs need to undergo technical and also clinical validation with respect to different features such as reliability and repeatability, the aim of the present analysis was to show clinical feasibility of online acquired, longitudinal DWI on a MR-Linac during fractionated RT for HNC and to characterize ADC dynamics. Previously, we reported relative repeatability coefficients (relRCs) of DWI in HNC on the MR-Linac to be 31% for GTV-P and 23% for GTV-N [14] and performed a comparison with a diagnostic scanner as a benchmark [13]. In this study, we observed ADCmean changes larger than these relRCs in weeks 3–4 for GTV-P and weeks 2–3 for GTV-N which correspond to the timepoints where significant ADC changes from 0 could be measured. Hence, week 3 seems to be a reasonable timepoint for ADC based interventions on a MR-Linac.

The current patient cohort is one of the largest to date reporting longitudinal ADC during RT and one of the first to report on online acquired serial DWI. As anticipated and in congruence with previously published data, GTV-P and GTV-N volume decreased during treatment [22, 32, 33]. El-Habashy et al. showed a decrease of GTV-P and GTV-N volumes during fractionated RT for 30 patients treated on an 1.5 T MR-Linac [34]. These volume changes built the rationale for prospective adaptive RT trials such as the MR-ADAPTOR [35, 36]. Our results also show a consistent increase in mean ADC in GTV-P and GTV-N, as reported by others in an offline setting, too [12, 15, 17, 22].

One major limitation is the inherent high technical uncertainty related to the acquisition and analysis of DWI. While echo-planar-imaging (EPI) for DWI has the advantage of fast image acquisition, it severely suffers from geometrical distortions caused by inhomogeneities of the magnetic field [37,38,39]. EPI applications in HNC patients are especially prone to geometric distortions because of many air-tissue boundaries in the respective area which cause magnetic field inhomogeneity and therefore geometric distortions. Improvements could be achieved by using image acquisition techniques like turbo-spin-echo (TSE) or split acquisition of fast spin echo signal (SPLICE), but with the downside of lower signal-to-noise ratio (SNR) and longer acquisition times [30, 38, 40]. Furthermore, comparability of ADC values and the translation of results from other studies is complicated due to inconsistencies in selected b-values for the DWI sequences and different algorithms used for ADC calculation, different magnetic field strengths and imaging techniques which all might affect the reported ADC values [30, 41, 42]. Another limitation is the small number of patients receiving their treatment within a planned period of time, limiting the potential prognostic value of the analysis as discussed above. Therefore, the MR-Linac consortium is working on publishing consensus DWI sequences or guidelines for the 1.5 T MR-Linac to form a prospective large group of patients, scanned with a dedicated DWI sequence and to make the translation and interpretation of ADC values from different studies easier and to shed further light on the role of DWI and its prognostic value for HNC.

In a preclinical trial, tumour subvolumes based on a cluster of ADC values were found to be significantly associated with radiation sensitivity and local tumor control [26]. In this clinical cohort this association was not reproduced, as the number of events was very small and moreover the reproducibility for serial DWI measurements might be worse than in a classical diagnostic setup and especially compared to a preclinical setting. In contrast to other reported data [15], we could not find a difference between OPC and non OPC or between HPV associated OPC and non-HPV associated OPC, possibly due to the small subgroups.

With no isolated local or regional recurrence, likely due to sampling bias, no correlation of DWI with outcome in this small cohort could be generated. Moreover, the established clinical risk groups for OPC could not be differentiated by the means of mean ADC or ADC-based HRS, potentially due to the small number of patients in each subgroup. Although no obvious differences could be seen between the subgroups in OPC or non-OPC, the high risk subgroup had larger tumor volumes and showed initially a higher HRS, consistent with the potential as a QIB, whereas intermediate and low risk subgroups had already a median HRS below the proposed threshold of 5.8 cm3 [27]. The missing correlation with outcome seen in the present analysis, emphasizes the need of larger cohorts of patients and a group of patients suffering from less selection bias. In addition, the HRS even in the high risk and the non-OPC subgroup decreased very rapidly and were below the established threshold of 5.8 cm3 already in week 2 for all but one of the patients (Fig. 6). This, together with the significant decrease in various parameters already early during treatment (Fig. 10) emphazises the need of repeated imaging to access the dynamic change of imaging based biomarkers and needs to be taken into account for future QIB-based clinical trials. Moreover, tumor volume in some cases was very small and the resolution of the DWI sequence may affect the analysis especially in the ADC-based HRS at late time points during RT. The seen GTV volume and ADC changes early during treatment offer a further explanation for the so far inconclusive and conflicting body of literature on the potential role of DWI as a prognostic or even predictive marker in HNC [15, 16, 18, 27, 43,44,45]. Nevertheless, the potential of repeated online quantitative imaging and anlysis of ADC-based subvolumes may facilitate the possibility for an ADC-based dose individualization and emphasizes the consideration of adaptive approaches during QIB-based trials.

Conclusions

In conclusion, this study shows the feasibility of online acquired, longitudinal DWI on a tumor level as well as in ADC-based subvolumes during the course of MR-guided RT for HNC, which might build the basis of ADC-based biological individualized online adaptive RT trials. The validation of these results in a prospective multicenter study seems to be an important next step.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due to inclusion in an ongoing prospective registry for secondary multi-center analysis but are available from the corresponding author on reasonable request.

Abbreviations

ADC:

Apparent diffusion coefficient

CT:

Computed tomography

DCE:

Dynamic contrast enhanced

DWI:

Diffusion-weighted imaging

EPI:

Echo-planar imaging

FDG:

Fluorodeoxyglucose

GTV-N:

Conspicious lymph nodes

GTV-P:

Primary tumor

HNC:

Head and neck cancer

HPV:

Human papilloma virus 16

HRS:

High-risk subvolume

MRI:

Magnetic resonance imaging

MR-Linac:

Magnetic resonance linear accelerator

OPC:

Oropharyngeal cancer

PET:

Positron emission tomography

QIB:

Quantitative imaging biomarker

relRC:

Relative repeatability coefficient

RT:

Radiotherapy

SNR:

Signal-to-noise ratio

SPLICE:

Split acquisition of fast spin echo signal

SS-EPI:

Single-shot echo-planar imaging

TSE:

Turbo-spin-echo

VOI:

Volume of interest

References

  1. Ang KK, Harris J, Wheeler R, Weber R, Rosenthal DI, Nguyen-Tân PF, et al. Human papillomavirus and survival of patients with oropharyngeal cancer. N Engl J Med. 2010;363:24–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Lassen P, Lacas B, Pignon JP, Trotti A, Zackrisson B, Zhang Q, et al. Prognostic impact of HPV-associated p16-expression and smoking status on outcomes following radiotherapy for oropharyngeal cancer: the MARCH-HPV project. Radiother Oncol. 2018;126:107–15.

    Article  CAS  PubMed  Google Scholar 

  3. Linge A, Lohaus F, Löck S, Nowak A, Gudziol V, Valentini C, et al. HPV status, cancer stem cell marker expression, hypoxia gene signatures and tumour volume identify good prognosis subgroups in patients with HNSCC after primary radiochemotherapy: a multicentre retrospective study of the German Cancer Consortium Radiation Oncology Group (DKTK-ROG). Radiother Oncol. 2016;121:364–73.

    Article  PubMed  Google Scholar 

  4. Berwouts D, Madani I, Duprez F, Olteanu AL, Vercauteren T, Boterberg T, et al. Long-term outcome of (18) F-fluorodeoxyglucose-positron emission tomography-guided dose painting for head and neck cancer: Matched case-control study. Head Neck. 2017;39:2264–75.

    Article  PubMed  Google Scholar 

  5. Lee NY, Sherman EJ, Schöder H, Wray R, Boyle JO, Singh B, et al. Hypoxia-directed treatment of human papillomavirus-related oropharyngeal carcinoma. J Clin Oncol. 2024;42:940–50.

    Article  CAS  PubMed  Google Scholar 

  6. Michaelidou A, Adjogatse D, Suh Y, Pike L, Thomas C, Woodley O, et al. 18F-FDG-PET in guided dose-painting with intensity modulated radiotherapy in oropharyngeal tumours: a phase I study (FiGaRO). Radiother Oncol. 2021;155:261–8.

    Article  CAS  PubMed  Google Scholar 

  7. Welz S, Paulsen F, Pfannenberg C, Reimold M, Reischl G, Nikolaou K, et al. Dose escalation to hypoxic subvolumes in head and neck cancer: a randomized phase II study using dynamic [(18)F]FMISO PET/CT. Radiother Oncol. 2022;171:30–6.

    Article  CAS  PubMed  Google Scholar 

  8. Zschaeck S, Löck S, Hofheinz F, Zips D, Saksø Mortensen L, Zöphel K, et al. Individual patient data meta-analysis of FMISO and FAZA hypoxia PET scans from head and neck cancer patients undergoing definitive radio-chemotherapy. Radiother Oncol. 2020;149:189–96.

    Article  CAS  PubMed  Google Scholar 

  9. Lavigne D, Ng SP, O’Sullivan B, Nguyen-Tan PF, Filion E, Létourneau-Guillon L, et al. Magnetic resonance-guided radiation therapy for head and neck cancers. Curr Oncol. 2022;29:8302–15.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wahid KA, Ahmed S, He R, van Dijk LV, Teuwen J, McDonald BA, et al. Evaluation of deep learning-based multiparametric MRI oropharyngeal primary tumor auto-segmentation and investigation of input channel effects: results from a prospective imaging registry. Clin Transl Radiat Oncol. 2022;32:6–14.

    CAS  PubMed  Google Scholar 

  11. Zukauskaite R, Rumley CN, Hansen CR, Jameson MG, Trada Y, Johansen J, et al. Delineation uncertainties of tumour volumes on MRI of head and neck cancer patients. Clin Transl Radiat Oncol. 2022;36:121–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Guerreiro F, van Houdt PJ, Navest RJM, Hoekstra N, de Jong M, Heijnen BJ, et al. Validation of quantitative magnetic resonance imaging techniques in head and neck healthy structures involved in the salivary and swallowing function: accuracy and repeatability. Phys Imaging Radiat Oncol. 2024;31:100608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Habrich J, Boeke S, Fritz V, Koerner E, Nikolaou K, Schick F, et al. Reproducibility of diffusion-weighted magnetic resonance imaging in head and neck cancer assessed on a 1.5 T MR-Linac and comparison to parallel measurements on a 3 T diagnostic scanner. Radiother Oncol. 2024;191:110046.

    Article  PubMed  Google Scholar 

  14. Habrich J, Boeke S, Nachbar M, Nikolaou K, Schick F, Gani C, et al. Repeatability of diffusion-weighted magnetic resonance imaging in head and neck cancer at a 1.5 T MR-Linac. Radiother Oncol. 2022;174:141–8.

    Article  PubMed  Google Scholar 

  15. Sijtsema ND, Lauwers I, Verduijn GM, Hoogeman MS, Poot DHJ, Hernandez-Tamames JA, et al. Relating pre-treatment non-Gaussian intravoxel incoherent motion diffusion-weighted imaging to human papillomavirus status and response in oropharyngeal carcinoma. Phys Imaging Radiat Oncol. 2024;30:100574.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Smits HJG, Vink SJ, de Ridder M, Philippens MEP, Dankbaar JW. Prognostic value of pretreatment radiological MRI variables and dynamic contrast-enhanced MRI on radiotherapy treatment outcome in laryngeal and hypopharyngeal tumors. Clin Transl Radiat Oncol. 2024;49:100857.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Trada Y, Keall P, Jameson M, Moses D, Lin P, Chlap P, et al. Changes in serial multiparametric MRI and FDG-PET/CT functional imaging during radiation therapy can predict treatment response in patients with head and neck cancer. Eur Radiol. 2023;33:8788–99.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Lambrecht M, Van Calster B, Vandecaveye V, De Keyzer F, Roebben I, Hermans R, et al. Integrating pretreatment diffusion weighted MRI into a multivariable prognostic model for head and neck squamous cell carcinoma. Radiother Oncol. 2014;110:429–34.

    Article  PubMed  Google Scholar 

  19. Connor S, Anjari M, Burd C, Guha A, Lei M, Guerrero-Urbano T, et al. The impact of Human Papilloma Virus status on the prediction of head and neck cancer chemoradiotherapy outcomes using the pre-treatment apparent diffusion coefficient. Br J Radiol. 2022;95:20210333.

    Article  PubMed  Google Scholar 

  20. Aliotta E, Paudyal R, Diplas B, Han J, Hu YC, Hun OhJ, et al. Multi-modality imaging parameters that predict rapid tumor regression in head and neck radiotherapy. Phys Imaging Radiat Oncol. 2024;31:100603.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Mohamed ASR, Abusaif A, He R, Wahid KA, Salama V, Youssef S, et al. Prospective validation of diffusion-weighted MRI as a biomarker of tumor response and oncologic outcomes in head and neck cancer: results from an observational biomarker pre-qualification study. Radiother Oncol. 2023;183:109641.

    Article  CAS  PubMed  Google Scholar 

  22. Ng SP, Cardenas CE, Bahig H, Elgohari B, Wang J, Johnson JM, et al. Changes in apparent diffusion coefficient (ADC) in serial weekly MRI during radiotherapy in patients with head and neck cancer: results from the PREDICT-HN study. Curr Oncol. 2022;29:6303–13.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Boeke S, Mönnich D, van Timmeren JE, Balermpas P. MR-guided radiotherapy for head and neck cancer: current developments, perspectives, and challenges. Front Oncol. 2021;11:616156.

    Article  PubMed  PubMed Central  Google Scholar 

  24. McDonald BA, Vedam S, Yang J, Wang J, Castillo P, Lee B, et al. Initial feasibility and clinical implementation of daily MR-guided adaptive head and neck cancer radiation therapy on a 1.5T MR-Linac system: prospective R-IDEAL 2a/2b systematic clinical evaluation of technical innovation. Int J Radiat Oncol Biol Phys. 2021;109:1606–18.

    Article  PubMed  Google Scholar 

  25. Thorwarth D, Ege M, Nachbar M, Mönnich D, Gani C, Zips D, et al. Quantitative magnetic resonance imaging on hybrid magnetic resonance linear accelerators: perspective on technical and clinical validation. Phys Imaging Radiat Oncol. 2020;16:69–73.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Boeke S, Winter RM, Leibfarth S, Krueger MA, Bowden G, Cotton J, et al. Machine learning identifies multi-parametric functional PET/MR imaging cluster to predict radiation resistance in preclinical head and neck cancer models. Eur J Nucl Med Mol Imaging. 2023;50:3084–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Winter RM, Boeke S, Leibfarth S, Habrich J, Clasen K, Nikolaou K, et al. Clinical validation of a prognostic preclinical magnetic resonance imaging biomarker for radiotherapy outcome in head-and-neck cancer. Radiother Oncol. 2025;204:110702.

    Article  CAS  Google Scholar 

  28. Kooreman ES, van Houdt PJ, Keesman R, Pos FJ, van Pelt VWJ, Nowee ME, et al. ADC measurements on the Unity MR-linac—a recommendation on behalf of the Elekta Unity MR-linac consortium. Radiother Oncol. 2020;153:106–13.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Stejskal EO. Use of spin echoes in a pulsed magnetic-field gradient to study anisotropic, restricted diffusion and flow. J Chem Phys. 1965;43:3597–603.

    Article  Google Scholar 

  30. McDonald BA, Salzillo T, Mulder S, Ahmed S, Dresner A, Preston K, et al. Prospective evaluation of in vivo and phantom repeatability and reproducibility of diffusion-weighted MRI sequences on 1.5 T MRI-linear accelerator (MR-Linac) and MR simulator devices for head and neck cancers. Radiother Oncol. 2023;185:109717.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Wong OL, Yuan J, Poon DMC, Chiu ST, Yang B, Chiu G, et al. Prostate diffusion-weighted imaging (DWI) in MR-guided radiotherapy: reproducibility assessment on 1.5 T MR-Linac and 1.5 T MR-simulator. Magn Reson Imaging. 2024;111:47–56.

    Article  PubMed  Google Scholar 

  32. Hamming-Vrieze O, van Kranen SR, Heemsbergen WD, Lange CAH, van den Brekel MWM, Verheij M, et al. Analysis of GTV reduction during radiotherapy for oropharyngeal cancer: implications for adaptive radiotherapy. Radiother Oncol. 2017;122:224–8.

    Article  PubMed  Google Scholar 

  33. Subesinghe M, Scarsbrook AF, Sourbron S, Wilson DJ, McDermott G, Speight R, et al. Alterations in anatomic and functional imaging parameters with repeated FDG PET-CT and MRI during radiotherapy for head and neck cancer: a pilot study. BMC Cancer. 2015;15:137.

    Article  PubMed  PubMed Central  Google Scholar 

  34. El-Habashy DM, Wahid KA, He R, McDonald B, Rigert J, Mulder SJ, et al. Longitudinal diffusion and volumetric kinetics of head and neck cancer magnetic resonance on a 1.5 T MR-Linear accelerator hybrid system: a prospective R-IDEAL Stage 2a imaging biomarker characterization/pre-qualification study. Clin Transl Radiat Oncol. 2023;42:100666.

    Google Scholar 

  35. Bahig H, Yuan Y, Mohamed ASR, Brock KK, Ng SP, Wang J, et al. Magnetic resonance-based response assessment and dose adaptation in human papilloma virus positive tumors of the oropharynx treated with radiotherapy (MR-ADAPTOR): an R-IDEAL stage 2a–2b/Bayesian phase II trial. Clin Transl Radiat Oncol. 2018;13:19–23.

    PubMed  PubMed Central  Google Scholar 

  36. Mohamed ASR, Bahig H, Aristophanous M, Blanchard P, Kamal M, Ding Y, et al. Prospective in silico study of the feasibility and dosimetric advantages of MRI-guided dose adaptation for human papillomavirus positive oropharyngeal cancer patients compared with standard IMRT. Clin Transl Radiat Oncol. 2018;11:11–8.

    PubMed  PubMed Central  Google Scholar 

  37. Jezzard P, Balaban RS. Correction for geometric distortion in echo planar images from B0 field variations. Magn Reson Med. 1995;34:65–73.

    Article  CAS  PubMed  Google Scholar 

  38. Sakamoto J, Yoshino N, Okochi K, Imaizumi A, Tetsumura A, Kurohara K, et al. Tissue characterization of head and neck lesions using diffusion-weighted MR imaging with SPLICE. Eur J Radiol. 2009;69:260–8.

    Article  PubMed  Google Scholar 

  39. Treiber JM, White NS, Steed TC, Bartsch H, Holland D, Farid N, et al. Characterization and correction of geometric distortions in 814 diffusion weighted images. PLoS ONE. 2016;11:e0152472.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Schakel T, Hoogduin JM, Terhaard CHJ, Philippens MEP. Technical note: diffusion-weighted MRI with minimal distortion in head-and-neck radiotherapy using a turbo spin echo acquisition method. Med Phys. 2017;44:4188–93.

    Article  CAS  PubMed  Google Scholar 

  41. Bisgaard ALH, Keesman R, van Lier A, Coolens C, van Houdt PJ, Tree A, et al. Recommendations for improved reproducibility of ADC derivation on behalf of the Elekta MRI-linac consortium image analysis working group. Radiother Oncol. 2023;186:109803.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Ermongkonchai T, Khor R, Wada M, Lau E, Xing DT, Ng SP. A review of diffusion-weighted magnetic resonance imaging in head and neck cancer patients for treatment evaluation and prediction of radiation-induced xerostomia. Radiat Oncol. 2023;18:20.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Hatakenaka M, Nakamura K, Yabuuchi H, Shioyama Y, Matsuo Y, Ohnishi K, et al. Pretreatment apparent diffusion coefficient of the primary lesion correlates with local failure in head-and-neck cancer treated with chemoradiotherapy or radiotherapy. Int J Radiat Oncol Biol Phys. 2011;81:339–45.

    Article  PubMed  Google Scholar 

  44. Hauser T, Essig M, Jensen A, Laun FB, Münter M, Maier-Hein KH, et al. Prediction of treatment response in head and neck carcinomas using IVIM-DWI: evaluation of lymph node metastasis. Eur J Radiol. 2014;83:783–7.

    Article  PubMed  Google Scholar 

  45. Peltenburg B, Driessen JP, Vasmel JE, Pameijer FA, Janssen LM, Terhaard CHJ, et al. Pretreatment ADC is not a prognostic factor for local recurrences in head and neck squamous cell carcinoma when clinical T-stage is known. Eur Radiol. 2020;30:1228–31.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We acknowledge support from the Open Access Publication Fund of the University of Tübingen.

Funding

Open Access funding enabled and organized by Projekt DEAL. This study received funding by the German Research Council (DFG Grants PAK 997/1: TH1528/6-1, SCHI498/14-1, NI707/7-1 and ZI 736/2-1).

Author information

Authors and Affiliations

Authors

Contributions

SB did a majority of the patients recruitement and data acquisition, contouring, selected analysis methods, interpreted the results and was a major contributor to wirting of the manuscript. JH selected the analysis methods, did data acquisition, solely analysed the data, interpretated the results and was a major contributor to writing the manuscript. SK helped acquiring the data. JB supported in the data acquisition. FS contributed substantially to the research concept, selection of methods and interpretation of the results. KN made major contributions to the research concept, selection of methods and interpretation of the results. JK helped with the data acquisition. CG contributed majorly to the patient recruitement and data acquisition. MN added to selection of methods, critical interpretation of results and writing of the manuscript. DZ substantially accounted for the research concept, selection of methods and interpretation of the results. DT added majorly to the research concept, selection of methods, interpretation of the results and writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jonas Habrich.

Ethics declarations

Ethics approval and consent to participate

The trial was approved by the local ethical committee (no. 659/2017BO1, NCT04172753) and all patients gave their written informant consent.

Consent for publication

Approval for publication of patient images was included in all patients’ written informant consent.

Competing interests

The Department of Radiation Oncology Tübingen receives financial and technical support by Elekta, Philips, Dr. Sennewald Medizintechnik, Kaiku Health, TheraPanacea, PTW Freiburg, ITV and Brainlab in the context of research cooperations. The Department of Diagnostic and Interventional Radiology, University of Tübingen reports institutional collaborations, consultation and grants supported Siemens Healthineers, Erlangen, Germany and Bayer Healthcare, Leverkusen, Germany outside this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Boeke, S., Habrich, J., Kübler, S. et al. Longitudinal assessment of diffusion-weighted imaging during magnetic resonance-guided radiotherapy in head and neck cancer. Radiat Oncol 20, 15 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13014-025-02589-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13014-025-02589-9

Keywords