Implementing Radiation Dose-Volume Liver Response in Biomechanical Deformable Image Registration

A vessel bifurcation liver CT landmark pair dataset for evaluating deformable image registration algorithms

PURPOSE

Evaluating deformable image registration (DIR) algorithms is vital for enhancing algorithm performance and gaining clinical acceptance. However, there is a notable lack of dependable DIR benchmark datasets for assessing DIR performance except for lung images. To address this gap, we aim to introduce our comprehensive liver computed tomography (CT) DIR landmark dataset library. This dataset is designed for efficient and quantitative evaluation of various DIR methods for liver CTs, paving the way for more accurate and reliable image registration techniques.

ACQUISITION AND VALIDATION METHODS

Forty CT liver image pairs were acquired from several publicly available image archives and authors' institutions under institutional review board (IRB) approval. The images were processed with a semi-automatic procedure to generate landmark pairs: (1) for each case, liver vessels were automatically segmented on one image; (2) landmarks were automatically detected at vessel bifurcations; (3) corresponding landmarks in the second image were placed using two deformable image registration methods to avoid algorithm-specific biases; (4) a comprehensive validation process based on quantitative evaluation and manual assessment was applied to reject outliers and ensure the landmarks' positional accuracy. This workflow resulted in an average of ∼56 landmark pairs per image pair, comprising a total of 2220 landmarks for 40 cases. The general landmarking accuracy of this procedure was evaluated using digital phantoms and manual landmark placement. The landmark pair target registration errors (TRE) on digital phantoms were 0.37 ± 0.26 and 0.55 ± 0.34 mm respectively for the two selected DIR algorithms used in our workflow, with 97% of landmark pairs having TREs below 1.5 mm. The distances from the calculated landmarks to the averaged manual placement were 1.27 ± 0.79 mm.

DATA FORMAT AND USAGE NOTES

All data, including image files and landmark information, are publicly available at Zenodo (https://zenodo.org/records/13738577). Instructions for using our data can be found on our GitHub page at https://github.com/deshanyang/Liver-DIR-QA.

POTENTIAL APPLICATIONS

The landmark dataset generated in this work is the first collection of large-scale liver CT DIR landmarks prepared on real patient images. This dataset can provide researchers with a dense set of ground truth benchmarks for the quantitative evaluation of DIR algorithms within the liver.

PET/CT-Based Absorbed Dose Maps in 90Y Selective Internal Radiation Therapy Correlate with Spatial Changes in Liver Function Derived from Dynamic MRI

Functional liver parenchyma can be damaged from treatment of liver malignancies with 90Y selective internal radiation therapy (SIRT). Evaluating functional parenchymal changes and developing an absorbed dose (AD)-toxicity model can assist the clinical management of patients receiving SIRT. We aimed to determine whether there is a correlation between 90Y PET AD voxel maps and spatial changes in the nontumoral liver (NTL) function derived from dynamic gadoxetic acid-enhanced MRI before and after SIRT. Methods: Dynamic gadoxetic acid-enhanced MRI scans were acquired before and after treatment for 11 patients undergoing 90Y SIRT. Gadoxetic acid uptake rate (k1) maps that directly quantify spatial liver parenchymal function were generated from MRI data. Voxel-based AD maps, derived from the 90Y PET/CT scans, were binned according to AD. Pre- and post-SIRT k1 maps were coregistered to the AD map. Absolute and percentage k1 loss in each bin was calculated as a measure of loss of liver function, and Spearman correlation coefficients between k1 loss and AD were evaluated for each patient. Average k1 loss over the patients was fit to a 3-parameter logistic function based on AD. Patients were further stratified into subgroups based on lesion type, baseline albumin-bilirubin scores and alanine transaminase levels, dose-volume effect, and number of SIRT treatments. Results: Significant positive correlations (ρ = 0.53-0.99, P < 0.001) between both absolute and percentage k1 loss and AD were observed in most patients (8/11). The average k1 loss over 9 patients also exhibited a significant strong correlation with AD (ρ ≥ 0.92, P < 0.001). The average percentage k1 loss of patients across AD bins was 28%, with a logistic function model demonstrating about a 25% k1 loss at about 100 Gy. Analysis between patient subgroups demonstrated that k1 loss was greater among patients with hepatocellular carcinoma, higher alanine transaminase levels, larger fractional volumes of NTL receiving an AD of 70 Gy or more, and sequential SIRT treatments. Conclusion: Novel application of multimodality imaging demonstrated a correlation between 90Y SIRT AD and spatial functional liver parenchymal degradation, indicating that a higher AD is associated with a larger loss of local hepatocyte function. With the developed response models, PET-derived AD maps can potentially be used prospectively to identify localized damage in liver and to enhance treatment strategies.

Fully automated deep learning based auto-contouring of liver segments and spleen on contrast-enhanced CT images

Manual delineation of liver segments on computed tomography (CT) images for primary/secondary liver cancer (LC) patients is time-intensive and prone to inter/intra-observer variability. Therefore, we developed a deep-learning-based model to auto-contour liver segments and spleen on contrast-enhanced CT (CECT) images. We trained two models using 3d patch-based attention U-Net ([Formula: see text] and 3d full resolution of nnU-Net ([Formula: see text] to determine the best architecture ([Formula: see text]. BA was used with vessels ([Formula: see text] and spleen ([Formula: see text] to assess the impact on segment contouring. Models were trained, validated, and tested on 160 ([Formula: see text]), 40 ([Formula: see text]), 33 ([Formula: see text]), 25 (CCH) and 20 (CPVE) CECT of LC patients. [Formula: see text] outperformed [Formula: see text] across all segments with median differences in Dice similarity coefficients (DSC) ranging 0.03-0.05 (p < 0.05). [Formula: see text], and [Formula: see text] were not statistically different (p > 0.05), however, both were slightly better than [Formula: see text] by DSC up to 0.02. The final model, [Formula: see text], showed a mean DSC of 0.89, 0.82, 0.88, 0.87, 0.96, and 0.95 for segments 1, 2, 3, 4, 5-8, and spleen, respectively on entire test sets. Qualitatively, more than 85% of cases showed a Likert score [Formula: see text] 3 on test sets. Our final model provides clinically acceptable contours of liver segments and spleen which are usable in treatment planning.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Analysis and prediction of liver volume change maps derived from computational tomography scans acquired pre- and post-radiation therapy

External beam radiation therapy (EBRT) of liver cancers can cause local liver atrophy as a result of tissue damage or hypertrophy as a result of liver regeneration. Predicting those volumetric changes would enable new strategies for liver function preservation during treatment planning. However, understanding of the spatial dose/volume relationship is still limited. This study leverages the use of deep learning-based segmentation and biomechanical deformable image registration (DIR) to analyze and predict this relationship. Pre- and Post-EBRT imaging data were collected for 100 patients treated for hepatocellular carcinomas, cholangiocarcinoma or CRC with intensity-modulated radiotherapy (IMRT) with prescription doses ranging from 50 to 100 Gy delivered in 10-28 fractions. For each patient, DIR between the portal and venous (PV) phase of a diagnostic computed tomography (CT) scan acquired before radiation therapy (RT) planning, and a PV phase of a diagnostic CT scan acquired after the end of RT (on average 147 ± 36 d) was performed to calculate Jacobian maps representing volume changes in the liver. These volume change maps were used: (i): to analyze the dose/volume relationship in the whole liver and individual Couinaud's segments; and (ii): to investigate the use of deep-learning to predict a Jacobian map solely based on the pre-RT diagnostic CT and planned dose distribution. Moderate correlations between mean equivalent dose in 2 Gy fractions (EQD2) and volume change was observed for all liver sub-regions analyzed individually with Pearson correlationrranging from -0.36 to -067. The predicted volume change maps showed a significantly stronger voxel-wise correlation with the DIR-based volume change maps than when considering the original EQD2 distribution (0.63 ± 0.24 versus 0.55 ± 23, respectively), demonstrating the ability of the proposed approach to establish complex relationships between planned dose and liver volume response months after treatment, which represents a promising prediction tool for the development of future adaptive and personalized liver radiation therapy strategies.

Ablative Margins of Colorectal Liver Metastases Using Deformable CT Image Registration and Autosegmentation

Background Confirming ablation completeness with sufficient ablative margin is critical for local tumor control following colorectal liver metastasis (CLM) ablation. An image-based confirmation method considering patient- and ablation-related biomechanical deformation is an unmet need. Purpose To evaluate a biomechanical deformable image registration (DIR) method for three-dimensional (3D) minimal ablative margin (MAM) quantification and the association with local disease progression following CT-guided CLM ablation. Materials and Methods This single-institution retrospective study included patients with CLM treated with CT-guided microwave or radiofrequency ablation from October 2015 to March 2020. A biomechanical DIR method with AI-based autosegmentation of liver, tumors, and ablation zones on CT images was applied for MAM quantification retrospectively. The per-tumor incidence of local disease progression was defined as residual tumor or local tumor progression. Factors associated with local disease progression were evaluated using the multivariable Fine-Gray subdistribution hazard model. Local disease progression sites were spatially localized with the tissue at risk for tumor progression (<5 mm) using a 3D ray-tracing method. Results Overall, 213 ablated CLMs (mean diameter, 1.4 cm) in 124 consecutive patients (mean age, 57 years ± 12 [SD]; 69 women) were evaluated, with a median follow-up interval of 25.8 months. In ablated CLMs, an MAM of 0 mm was depicted in 14.6% (31 of 213), from greater than 0 to less than 5 mm in 40.4% (86 of 213), and greater than or equal to 5 mm in 45.1% (96 of 213). The 2-year cumulative incidence of local disease progression was 72% for 0 mm and 12% for greater than 0 to less than 5 mm. No local disease progression was observed for an MAM greater than or equal to 5 mm. Among 117 tumors with an MAM less than 5 mm, 36 had local disease progression and 30 were spatially localized within the tissue at risk for tumor progression. On multivariable analysis, an MAM of 0 mm (subdistribution hazard ratio, 23.3; 95% CI: 10.8, 50.5; P < .001) was independently associated with local disease progression. Conclusion Biomechanical deformable image registration and autosegmentation on CT images enabled identification and spatial localization of colorectal liver metastases at risk for local disease progression following ablation, with a minimal ablative margin greater than or equal to 5 mm as the optimal end point. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Sofocleous in this issue.

Incorporating the synthetic CT image for improving the performance of deformable image registration between planning CT and cone-beam CT

Objective

To develop a contrast learning-based generative (CLG) model for the generation of high-quality synthetic computed tomography (sCT) from low-quality cone-beam CT (CBCT). The CLG model improves the performance of deformable image registration (DIR).

Methods

This study included 100 post-breast-conserving patients with the pCT images, CBCT images, and the target contours, which the physicians delineated. The CT images were generated from CBCT images via the proposed CLG model. We used the Sct images as the fixed images instead of the CBCT images to achieve the multi-modality image registration accurately. The deformation vector field is applied to propagate the target contour from the pCT to CBCT to realize the automatic target segmentation on CBCT images. We calculate the Dice similarity coefficient (DSC), 95 % Hausdorff distance (HD95), and average surface distance (ASD) between the prediction and reference segmentation to evaluate the proposed method.

Results

The DSC, HD95, and ASD of the target contours with the proposed method were 0.87 ± 0.04, 4.55 ± 2.18, and 1.41 ± 0.56, respectively. Compared with the traditional method without the synthetic CT assisted (0.86 ± 0.05, 5.17 ± 2.60, and 1.55 ± 0.72), the proposed method was outperformed, especially in the soft tissue target, such as the tumor bed region.

Conclusion

The CLG model proposed in this study can create the high-quality sCT from low-quality CBCT and improve the performance of DIR between the CBCT and the pCT. The target segmentation accuracy is better than using the traditional DIR.

Relationship between absorbed dose and changes in liver volume after chemoradiotherapy for esophageal cancer

PURPOSE

The liver is the largest organ in the abdomen and is often irradiated in radiotherapy for non-hepatic malignancies. As most of the studies on changes in liver volume are on hepatocellular carcinoma based on liver dysfunction, there are few studies on healthy liver. In this study, we investigated the relationship between absorbed dose and changes in liver volume after chemoradiotherapy for esophageal cancer in patients without apparent pre-treatment liver dysfunction.

MATERIALS AND METHODS

Liver volume was compared between pre-treatment, acute (< 4 months) and late post-treatment (≥ 4 and < 13 months) phases in 12 patients using abdominal plain CT images. Volume changes were evaluated separately for the right and left lobes. We investigated the relationship between the volume change and V_xGy (percentage of volume received x Gy or more dose). In addition, volume change for each absorbed dose was investigated using deformable image registration.

RESULTS

The volume of the left lobe showed a significant decrease between pre-treatment and acute post-treatment phases (p < 0.001), while the volume of right lobe and between acute and late post-treatment phase of left lobe did not. The mean value of the volume reduction rate of the left lobe was 51.1% and equivalent to the mean value of V_30Gy. As a result of the volume change for each absorbed dose, the volume reduction rate increased as the absorbed dose increased, and a significant volume loss was observed at doses above 11 Gy.

CONCLUSION

Volume of the liver significantly decreased only in the acute phase after chemoradiotherapy for esophageal cancer. The tolerable dose for a healthy liver is generally considered to be 30 Gy, but attention should be paid to lower doses to avoid radiation-induced liver injury.

Leveraging deep learning-based segmentation and contours-driven deformable registration for dose accumulation in abdominal structures

Purpose

Discrepancies between planned and delivered dose to GI structures during radiation therapy (RT) of liver cancer may hamper the prediction of treatment outcomes. The purpose of this study is to develop a streamlined workflow for dose accumulation in a treatment planning system (TPS) during liver image-guided RT and to assess its accuracy when using different deformable image registration (DIR) algorithms.

Materials and Methods

Fifty-six patients with primary and metastatic liver cancer treated with external beam radiotherapy guided by daily CT-on-rails (CTOR) were retrospectively analyzed. The liver, stomach and duodenum contours were auto-segmented on all planning CTs and daily CTORs using deep-learning methods. Dose accumulation was performed for each patient using scripting functionalities of the TPS and considering three available DIR algorithms based on: (i) image intensities only; (ii) intensities + contours; (iii) a biomechanical model (contours only). Planned and accumulated doses were converted to equivalent dose in 2Gy (EQD2) and normal tissue complication probabilities (NTCP) were calculated for the stomach and duodenum. Dosimetric indexes for the normal liver, GTV, stomach and duodenum and the NTCP values were exported from the TPS for analysis of the discrepancies between planned and the different accumulated doses.

Results

Deep learning segmentation of the stomach and duodenum enabled considerable acceleration of the dose accumulation process for the 56 patients. Differences between accumulated and planned doses were analyzed considering the 3 DIR methods. For the normal liver, stomach and duodenum, the distribution of the 56 differences in maximum doses (D2%) presented a significantly higher variance when a contour-driven DIR method was used instead of the intensity only-based method. Comparing the two contour-driven DIR methods, differences in accumulated minimum doses (D98%) in the GTV were >2Gy for 15 (27%) of the patients. Considering accumulated dose instead of planned dose in standard NTCP models of the duodenum demonstrated a high sensitivity of the duodenum toxicity risk to these dose discrepancies, whereas smaller variations were observed for the stomach.

Conclusion

This study demonstrated a successful implementation of an automatic workflow for dose accumulation during liver cancer RT in a commercial TPS. The use of contour-driven DIR methods led to larger discrepancies between planned and accumulated doses in comparison to using an intensity only based DIR method, suggesting a better capability of these approaches in estimating complex deformations of the GI organs.

Effect of the normal liver mean dose on intrahepatic recurrence in patients with hepatocellular carcinoma after receiving liver stereotactic body radiation therapy

•

The normal liver mean dosage is rarely studied in its association with metastasis related survival.

•

Higher normal liver mean dosage was related to a better intra-hepatic progression free survival.

•

No associations were found between normal liver mean dose with Local control rate, Overall Survival, Distant Metastasis Free Survival, Progression Free Survival.

Background and purpose

This study aims to evaluate whether dosimetric parameters affect the intrahepatic out-field recurrence or distant metastasis-free survival following the stereotactic body radiation therapy (SBRT) in patients with hepatocellular carcinoma (HCC).

Materials and methods

A total of 76 patients with HCC who were treated with SBRT from January 2015 to May 2020 were included in this retrospective study. The main clinical endpoints considered were intrahepatic out-field free survival (OutFFS) and distant metastasis-free survival (DMFS). The target parameters and the liver were documented including tumor diameters, gross tumor volume (GTV), Liver minus GTV volume (LGV), and Liver minus GTV mean dose (LGD). Multivariable Cox regression with forward stepwise selection was performed to identify independent risk factors for OutFFS and DMFS. Maximally selected rank statistics were used to determine the most informative cut-off value for age and LGD.

Results

The median follow-up was 28.2 months (range, 7.7–74.5 months). LGD higher than 12.54 Gy [HR, 0.861(0.747–0.993); p = 0.040] and age greater than 67-year-old [HR, 0.966(0.937–0.997); p = 0.030] are two independent predictors of OutFFS, previous TACE treatment [HR, 0.117(0.015–0.891); p = 0.038] was an independent predictor of DMFS.

Conclusions

The results of this study suggested that the higher the dose received by the normal liver (greater than 12.54 Gy) the better the intrahepatic out-field recurrence-free survival (RFS) rate. Further study is warranted to confirm and to better understand this phenomenon.

Detection of vessel bifurcations in CT scans for automatic objective assessment of deformable image registration accuracy

PURPOSE

Objective assessment of deformable image registration (DIR) accuracy often relies on the identification of anatomical landmarks in image pairs, a manual process known to be extremely time-expensive. The goal of this study is to propose a method to automatically detect vessel bifurcations in images and assess their use for the computation of target registration errors (TREs).

MATERIALS AND METHODS

Three image datasets were retrospectively analyzed. The first dataset included 10 pairs of inhale/exhale phases from lung 4DCTs and full inhale and exhale breath-hold CT scans from 10 patients presenting with chronic obstructive pulmonary disease, with 300 corresponding landmarks available for each case (DIR-Lab). The second dataset included six pairs of inhale/exhale phases from lung 4DCTs (POPI dataset), with 100 pairs of landmarks for each case. The third dataset included 28 pairs of pre/post-radiotherapy liver contrast-enhanced CT scans, each with five manually picked vessel bifurcation correspondences. For all images, the vasculature was autosegmented by computing and thresholding a vesselness image. Images of the vasculature centerline were computed, and bifurcations were detected based on centerline voxel neighbors' count. The vasculature segmentations were independently registered using a Demons algorithm between representations of their surface with distance maps. Detected bifurcations were considered as corresponding when distant by less than 5 mm after vasculature DIR. The selected pairs of bifurcations were used to calculate TRE after registration of the images considering three algorithms: rigid registration, Anaconda, and a Demons algorithm. For comparison with the ground truth, TRE values calculated using the automatically detected correspondences were interpolated in the whole organs to generate TRE maps. The performance of the method in automatically calculating TRE after image registration was quantified by measuring the correlation with the TRE obtained when using the ground truth landmarks.

RESULTS

The median Pearson correlation coefficients between ground truth TRE and corresponding values in the generated TRE maps were r = 0.81 and r = 0.67 for the lung and liver cases, respectively. The correlation coefficients between mean TRE for each case were r = 0.99 and r = 0.64 for the lung and liver cases, respectively.

CONCLUSION

For lungs or liver CT scans DIR, a strong correlation was obtained between TRE calculated using manually picked or landmarks automatically detected with the proposed method. This tool should be particularly useful in studies requiring assessing the reliability of a high number of DIRs.

Rigid and Deformable Image Registration for Radiation Therapy: A Self-Study Evaluation Guide for NRG Oncology Clinical Trial Participation

PURPOSE

The registration of multiple imaging studies to radiation therapy computed tomography simulation, including magnetic resonance imaging, positron emission tomography-computed tomography, etc. is a widely used strategy in radiation oncology treatment planning, and these registrations have valuable roles in image guidance, dose composition/accumulation, and treatment delivery adaptation. The NRG Oncology Medical Physics subcommittee formed a working group to investigate feasible workflows for a self-study credentialing process of image registration commissioning.

METHODS AND MATERIALS

The American Association of Physicists in Medicine (AAPM) Task Group 132 (TG132) report on the use of image registration and fusion algorithms in radiation therapy provides basic guidelines for quality assurance and quality control of the image registration algorithms and the overall clinical process. The report recommends a series of tests and the corresponding metrics that should be evaluated and reported during commissioning and routine quality assurance, as well as a set of recommendations for vendors. The NRG Oncology medical physics subcommittee working group found incompatibility of some digital phantoms with commercial systems. Thus, there is still a need to provide further recommendations in terms of compatible digital phantoms, clinical feasible workflow, and achievable thresholds, especially for future clinical trials involving deformable image registration algorithms. Nine institutions participated and evaluated 4 commonly used commercial imaging registration software and various versions in the field of radiation oncology.

RESULTS AND CONCLUSIONS

The NRG Oncology Working Group on image registration commissioning herein provides recommendations on the use of digital phantom/data sets and analytical software access for institutions and clinics to perform their own self-study evaluation of commercial imaging systems that might be employed for coregistration in radiation therapy treatment planning and image guidance procedures. Evaluation metrics and their corresponding values were given as guidelines to establish practical tolerances. Vendor compliance for image registration commissioning was evaluated, and recommendations were given for future development.

Biomechanical modeling of radiation dose-induced volumetric changes of the parotid glands for deformable image registration

PURPOSE

Early animal studies suggest that parotid gland (PG) toxicity prediction could be improved by an accurate estimation of the radiation dose to sub-regions of the PG. Translation to clinical investigation requires voxel-level dose accumulation in this organ that responds volumetrically throughout treatment. To date, deformable image registration (DIR) has been evaluated for the PG using only surface alignment. We sought to develop and evaluate an advanced DIR technique capable of modeling these complex PG volume changes over the course of radiation therapy.

MATERIALS AND METHODS

Planning and mid-treatment magnetic resonance images from 19 patients and computed tomography images from nine patients who underwent radiation therapy for head and neck cancer were retrospectively evaluated. A finite element model (FEM)-based DIR algorithm was applied between the corresponding pairs of images, based on boundary conditions on the PG surfaces only (Morfeus-spatial). To investigate an anticipated improvement in accuracy, we added a population model-based thermal expansion coefficient to simulate the dose distribution effect on the volume change inside the glands (Morfeus-spatialDose). The model accuracy was quantified using target registration error for magnetic resonance images, where corresponding anatomical landmarks could be identified. The potential clinical impact was evaluated using differences in mean dose, median dose, D98, and D50 of the PGs.

RESULTS

In the magnetic resonance images, the mean (±standard deviation) target registration error significantly reduced by 0.25 ± 0.38 mm (p = 0.01) when using Morfeus-spatialDose instead of Morfeus-spatial. In the computed tomography images, differences in the mean dose, median dose, D98, and D50 of the PGs reached 2.9 ± 0.8, 3.8, 4.1, and 3.8 Gy, respectively, between Morfeus-spatial and Morfeus-spatialDose.

CONCLUSION

Differences between Morfeus-spatial and Morfeus-spatialDose may be impactful when considering high-dose gradients of radiation in the PGs. The proposed DIR model can allow more accurate PG alignment than the standard model and improve dose estimation and toxicity prediction modeling.

Challenges in Reirradiation of Intrahepatic Tumors

Definitive reirradiation using a stereotactic technique is an effective local treatment option for both recurrent liver metastases and recurrent primary liver cancers. The tolerance of the liver, bile ducts, and surrounding gastrointestinal luminal organs must be respected to ensure safe retreatment. The risks associated with retreatment to these organs must be carefully balanced with the probability of clinical benefit. We present 2 cases for consideration of repeat irradiation along with the opinions of 4 experts, along with conclusions about recommendations.

Dosimetric Analysis and Normal-Tissue Complication Probability Modeling of Child-Pugh Score and Albumin-Bilirubin Grade Increase After Hepatic Irradiation

PURPOSE

This study aimed to develop robust normal-tissue complication probability (NTCP) models for patients with hepatocellular carcinoma treated with radiation therapy (RT) using Child-Pugh (CP) score and albumin-bilirubin (ALBI) grade increase as endpoints for hepatic toxicity.

METHODS AND MATERIALS

Data from 108 patients with hepatocellular carcinoma treated with RT between 2008 and 2017 were evaluated, of which 47 patients (44%) were treated with proton RT. Of these patients, 29 received stereotactic body RT and 79 moderately hypofractionated RT to median physical tumor doses of 43 Gy in 5 fractions and 59 Gy in 15 fractions, respectively. A generalized Lyman-Kutcher-Berman (LKB) model was used to model the NTCP using 2 clinical endpoints, both evaluated at 3 months after RT: CP score increase of ≥2 and ALBI grade increase of ≥1 from the pre-RT baseline. Confidence intervals on LKB fit parameters were determined using bootstrap resampling.

RESULTS

Compared with previous NTCP models, this study found a stronger correlation between normal liver volume receiving low doses of radiation (5-10 Gy) and a CP score or ALBI grade increase. A CP score increase exhibited a stronger correlation to normal liver volumes irradiated than an ALBI grade increase. LKB models for CP increase found values for the volume-effect parameter of a = 0.06 for all patients, and a = 0.02/0.09 when fit to photon/proton patients separately. Subset analyses for patients with superior initial liver functions showed consistent dose-volume effects (a = 0.1) and consistent dose-response relationships.

CONCLUSIONS

This study presents an update of liver NTCP models in the era of modern RT techniques using relevant endpoints of hepatic toxicity, CP score and ALBI grade increase. The results show a stronger influence of low-dose bath on hepatic toxicity than those found in previous studies, indicating that RT techniques that minimize the low-dose bath may be beneficial for patients.

Vasculature-Driven Biomechanical Deformable Image Registration of Longitudinal Liver Cholangiocarcinoma Computed Tomographic Scans

PURPOSE

Deformable image registration (DIR) of longitudinal liver cancer computed tomographic (CT) images can be challenging owing to anatomic changes caused by radiation therapy (RT) or disease progression. We propose a workflow for the DIR of longitudinal contrast-enhanced CT scans of liver cancer based on a biomechanical model of the liver driven by boundary conditions on the liver surface and centerline of an autosegmentation of the vasculature.

METHODS AND MATERIALS

Pre- and post-RT CT scans acquired with a median gap of 112 (32-217) days for 28 patients who underwent RT for intrahepatic cholangiocarcinoma were retrospectively analyzed. For each patient, 5 corresponding anatomic landmarks in pre- and post-RT scans were identified in the liver by a clinical expert for evaluation of the accuracy of different DIR strategies. The first strategy corresponded to the use of a biomechanical model-based DIR method with boundary conditions specified on the liver surface (BM_DIR). The second strategy corresponded to the use of an expansion of BM_DIR consisting of the auto-segmentation of the liver vasculature to determine additional boundary conditions in the biomechanical model (BM_DIR_VBC). The 2 strategies were also compared with an intensity-based DIR strategy using a Demons algorithms.

RESULTS

The group mean target registration errors were 12.4 ± 7.5, 7.7 ± 3.7 and 4.4 ± 2.5 mm, for the Demons, BM_DIR and BM_DIR_VBC, respectively.

CONCLUSIONS

In regard to the large and complex deformation observed in this study and the achieved accuracy of 4.4 mm, the proposed BM_DIR_VBC method might reveal itself as a valuable tool in future studies on the relationship between delivered dose and treatment outcome.

Accuracy of deformable image registration techniques for alignment of longitudinal cholangiocarcinoma CT images

PURPOSE

Response assessment of radiotherapy for the treatment of intrahepatic cholangiocarcinoma (IHCC) across longitudinal images is challenging due to anatomical changes. Advanced deformable image registration (DIR) techniques are required to correlate corresponding tissues across time. In this study, the accuracy of five commercially available DIR algorithms in four treatment planning systems (TPS) was investigated for the registration of planning images with posttreatment follow-up images for response assessment or re-treatment purposes.

METHODS

Twenty-nine IHCC patients treated with hypofractionated radiotherapy and with pretreatment and posttreatment contrast-enhanced computed tomography (CT) images were analyzed. Liver segmentations were semiautomatically generated on all CTs and the posttreatment CT was then registered to the pretreatment CT using five commercially available algorithms (Demons, B-splines, salient feature-based, anatomically constrained and finite element-based) in four TPSs. This was followed by an in-depth analysis of 10 DIR strategies (plus global and liver-focused rigid registration) in one of the TPSs. Eight of the strategies were variants of the anatomically constrained DIR while the two were based on a finite element-based biomechanical registration. The anatomically constrained techniques were combinations of: (a) initializations with the two rigid registrations; (b) two similarity metrics - correlation coefficient (CC) and mutual information (MI); and (c) with and without a controlling region of interest (ROI) - the liver. The finite element-based techniques were initialized by the two rigid registrations. The accuracy of each registration was evaluated using target registration error (TRE) based on identified vessel bifurcations. The results were statistically analyzed with a one-way analysis of variance (ANOVA) and pairwise comparison tests. Stratified analysis was conducted on the inter-TPS data (plus the liver-focused rigid registration) using treatment volume changes, slice thickness, time between scans, and abnormal lab values as stratifying factors.

RESULTS

The complex deformation observed following treatment resulted in average TRE exceeding the image voxel size for all techniques. For the inter-TPS comparison, the Demons algorithm had the lowest TRE, which was significantly superior to all the other algorithms. The respective mean (standard deviation) TRE (in mm) for the Demons, B-splines, salient feature-based, anatomically constrained, and finite element-based algorithms were 4.6 (2.0), 7.4 (2.7), 7.2 (2.6), 6.3 (2.3), and 7.5 (4.0). In the follow-up comparison of the anatomically constrained DIR, the strategy with liver-focused rigid registration initialization, CC as similarity metric and liver as a controlling ROI had the lowest mean TRE - 6.0 (2.0). The maximum TRE for all techniques exceeded 10 mm. Selection of DIR strategy was found to be a statistically significant factor for registration accuracy. Tumor volume change had a significant effect on TRE for finite element-based registration and B-splines DIR. Time between scans had a substantial effect on TRE for all registrations but was only significant for liver-focused rigid, finite element-based and salient feature-based DIRs.

CONCLUSIONS

This study demonstrates the limitations of commercially available DIR techniques in TPSs for alignment of longitudinal images of liver cancer presenting complex anatomical changes including local hypertrophy and fibrosis/necrosis. DIR in this setting should be used with caution and careful evaluation.

Comparison of regional lung perfusion response on longitudinal MAA SPECT/CT in lung cancer patients treated with and without functional tissue-avoidance radiation therapy

OBJECTIVE

The effect of functional lung avoidance planning on radiation dose-dependent changes in regional lung perfusion is unknown. We characterized dose-perfusion response on longitudinal perfusion single photon emission computed tomography (SPECT)/CT in two cohorts of lung cancer patients treated with and without functional lung avoidance techniques.

METHODS

The study included 28 primary lung cancer patients: 20 from interventional (NCT02773238) (FLARE-RT) and eight from observational (NCT01982123) (LUNG-RT) clinical trials. FLARE-RT treatment plans included perfused lung dose constraints while LUNG-RT plans adhered to clinical standards. Pre- and 3 month post-treatment macro-aggregated albumin (MAA) SPECT/CT scans were rigidly co-registered to planning four-dimensional CT scans. Tumour-subtracted lung dose was converted to EQD2 and sorted into 5 Gy bins. Mean dose and percent change between pre/post-RT MAA-SPECT uptake (%ΔPERF), normalized to total tumour-subtracted lung uptake, were calculated in each binned dose region. Perfusion frequency histograms of pre/post-RT MAA-SPECT were analyzed. Dose-response data were parameterized by sigmoid logistic functions to estimate maximum perfusion increase (%ΔPERFmaxincrease), maximum perfusion decrease (%ΔPERFmaxdecrease), dose midpoint (Dmid), and dose-response slope (k).

RESULTS

Differences in MAA perfusion frequency distribution shape between time points were observed in 11/20 (55%) FLARE-RT and 2/8 (25%) LUNG-RT patients (p < 0.05). FLARE-RT dose response was characterized by >10% perfusion increase in the 0-5 Gy dose bin for 8/20 patients (%ΔPERFmaxincrease = 10-40%), which was not observed in any LUNG-RT patients (p = 0.03). The dose midpoint Dmid at which relative perfusion declined by 50% trended higher in FLARE-RT compared to LUNG-RT cohorts (35 GyEQD2 vs 21 GyEQD2, p = 0.09), while the dose-response slope k was similar between FLARE-RT and LUNG-RT cohorts (3.1-3.2, p = 0.86).

CONCLUSION

Functional lung avoidance planning may promote increased post-treatment perfusion in low dose regions for select patients, though inter-patient variability remains high in unbalanced cohorts. These preliminary findings form testable hypotheses that warrant subsequent validation in larger cohorts within randomized or case-matched control investigations.

ADVANCES IN KNOWLEDGE

This novel preliminary study reports differences in dose-response relationships between patients receiving functional lung avoidance radiation therapy (FLARE-RT) and those receiving conventionally planned radiation therapy (LUNG-RT). Following further validation and testing of these effects in larger patient populations, individualized estimation of regional lung perfusion dose-response may help refine future risk-adaptive strategies to minimize lung function deficits and toxicity incidence.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.