Back to Journals » OncoTargets and Therapy » Volume 12

Bio-physic constraint model using spatial registration of delta 18F-fluorodeoxyglucose positron emission tomography/computed tomography images for predicting radiation pneumonitis in esophageal squamous cell carcinoma patients receiving neoadjuvant chemoradiation

Authors Hou TC , Dai KY, Wu MC , Hua KL , Tai HC, Huang WC, Chen YJ 

Received 18 February 2019

Accepted for publication 12 July 2019

Published 13 August 2019 Volume 2019:12 Pages 6439—6451

DOI https://doi.org/10.2147/OTT.S205803

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 2

Editor who approved publication: Prof. Dr. Nicola Silvestris



Tien-Chi Hou,1 Kun-Yao Dai,1 Ming-Che Wu,2 Kai-Lung Hua,3 Hung-Chi Tai,1 Wen-Chien Huang,4 Yu-Jen Chen1,5

1Department of Radiation Oncology, Mackay Memorial Hospital, Taipei, Taiwan; 2Department of Nuclear Medicine, Mackay Memorial Hospital, Taipei, Taiwan; 3Department of Computer Science and Information Engineering, National Taiwan University of Science and Technology, Taipei, Taiwan; 4Department of Surgery, Division of Thoracic Surgery, Mackay Memorial Hospital, Taipei City 10449, Taiwan; 5Department of Medical Research, China Medical University Hospital, Taichung 40402, Taiwan

Purpose: This study integrated clinical outcomes and radiomics of advanced thoracic esophageal squamous cell carcinoma patients receiving neoadjuvant concurrent chemoradiotherapy (NACCRT) to establish a novel constraint model for predicting radiation pneumonitis (RP).
Patients and methods: We conducted a retrospective review for thoracic advanced esophageal cancer patients who received NACCRT. From 2013 to 2018, 89 patients were eligible for review. Staging workup and response evaluation included positron emission tomography/computed tomography (PET/CT) scans and endoscopic ultrasound. Patients received RT with 48 Gy to gross tumor and 43.2 Gy to elective nodal area in simultaneous integrated boost method divided in 24 fractions. Weekly platinum-based chemotherapy was administered concurrently. Side effects were evaluated using CTCAE v4. Images of 2-fluoro-2-deoxyglucose PET/CT before and after NACCRT were registered to planning CT images to create a region of interest for dosimetry parameters that spatially matched RP-related regions, including V10, V20, V50%, V27, and V30. Correlation between bio-physic parameters and toxicity was used to establish a constraint model for avoiding RP.
Results: Among the investigated cohort, clinical downstaging, complete pathological response, and 5-year overall survival rates were 59.6%, 40%, and 34.4%, respectively. Multivariate logistic regression analysis demonstrated that each individual set standardized uptake value ratios (SUVRs), neither pre- nor post-NACCRT, was not predictive. Interestingly, cutoff increments of 6.2% and 8.9% in SUVRs (delta-SUVR) in registered V20 and V27 regions were powerful predictors for acute and chronic RP, respectively.
Conclusion: Spatial registration of metabolic and planning CT images with delta-radiomics analysis using fore-and-aft image sets can establish a unique bio-physic prediction model for avoiding RP in esophageal cancer patients receiving NACCRT.

Keywords: esophageal cancer, neoadjuvant concurrent chemoradiation, radiation pneumonitis, PET/CT, constraint model

Introduction

Esophageal cancer is the eighth most common cancer, causing about 400,000 deaths worldwide every year. The 5-year overall survival rate for esophageal cancer ranges from 15% to 25%.1 It is difficult to detect early and is usually diagnosed at advanced stages.2 Adenocarcinoma and squamous cell carcinoma are the main histological types of esophageal cancer with significantly heterogenous distribution between occidental and oriental worlds.2,3 Squamous cell carcinoma accounts for more than 90% of all esophageal cancer cases in the pan-Asia region and appears to have poorer prognoses and earlier lymph node metastasis in comparison to adenocarcinoma.46 Nearly 50% of locoregional esophageal cancers are unresectable or borderline resectable disease at diagnosis.7 Neoadjuvant treatment, especially chemoradiation, has been shown to increase the R0 resection rate with acceptable pathological complete response (pCR) ratio and improved survival in locally advanced disease. Therefore, several guidelines recommend neoadjuvant chemoradiation for downstaging before radical surgery.8,9 However, subjects of current reported neoadjuvant chemoradiation therapy (NACCRT) trials were recruited from occidental countries and mainly consisted of adenocarcinoma cases.10,11 Therefore, conclusions from these studies do not provide strong supporting evidence for squamous cell carcinoma.12

Successful radical surgery at 5–8 weeks after completion of NACCRT is critical to the clinical outcome of advanced esophageal cancer patients.13 Therefore, development of acute radiation complications, especially radiation pneumonitis (RP), may compromise the prognosis of this multi-modality treatment strategy.14 The radiotherapy (RT) field design has to cover gross tumors and the regional lymphatic basin; therefore, it is inevitable to irradiate central parts of the lung adjacent to para-esophageal tissue. As a result, RT may increase the risk of RP, which may impair the success of radical surgery and postoperative pulmonary rehabilitation.15,16 Therefore, developing a useful tool for predicting RP is clinically important.

Dose-volume histogram (DVH) generated from RT treatment planning system is widely used as a pre-treatment physic constraint model to evaluate dose distribution for tumor coverage and normal tissue avoidance. Several radiation physic studies demonstrated the role of DVH in predicting the risk of RP.1719 However, DVH-based predictors for RP remain unsatisfactory due to concerns such as having markedly individual variations (such as smoking status, tuberculosis, asthma, chronic obstructive pulmonary disease, and other pulmonary comorbidities), confounding by institutional facilities, and no actual engagement of biological effects. The 2-deoxy-2-[fluorine-18] fluoro-D-glucose positron emission tomography integrated with computed tomography (18F-FDG PET/CT) scan is preferred for the evaluation of treatment responses in esophageal cancer patients after NACCRT.8,9 Song et al reported that changes in standardized uptake value (SUV) levels of gross tumors were associated with metabolic responses and pathological complete resection ratios.20 Cunliffe et al quantified radiomic changes in serial CT scans and correlated RT doses to RP incidence.21 Recent studies used functional lung heterogeneity in PET/CT imaging before treatment to predict RP among patients receiving definitive thoracic radiation.22,23 However, precise and easily used posttreatment radiomics-based predictors of RP are still in development. Given that acute radiation reactions in the lung are a dynamic process and that inflammation can be assessed by dynamic metabolic imaging (before treatment and 5–8 weeks after completion of NACCRT for preoperative survey), it is of interest to explore the combination of CT and metabolic imaging modalities with RT dosimetry parameters to analyze their correlation to RP with regards to spatial registration and delta calculations.

In the present study, we quantified the lung tissue reaction at different RT dose levels and volumes using fore-and-aft 18F-FDG PET/CT scan images of patients who received NACCRT for advanced thoracic esophageal cancer. Our goal was to establish a unique and clinically applicable constraint model for RP and lung fibrosis.

Materials and methods

Patient demographic data and workup

We conducted a retrospective review of advanced thoracic esophageal cancer patients who received NACCRT in our institution from September 2013 to September 2018. Clinical stages with lymph node positive or >T2 disease without tissue-proven distant metastasis were eligible for this study. All eligible cohorts were required to have an Eastern Cooperative Oncology Group performance status scale of between 0 and 1. Patients with previous chemotherapy or thoracic RT were excluded. Eighty-eight patients were eligible for assessment of clinical outcome and toxicity. The cancer staging and post-NACCRT surveillances were conducted using CT scan, endoscopic ultrasound with bronchoscope (optional), esophago-gastro-duodenoscopy, and whole body 18F-FDG PET/CT. Posttreatment imaging studies, CT or 18F-FDG PET/CT, were assessed 5–8 weeks after the completion of NACCRT, to avoid treatment-related interference. The clinical and pathological staging was determined according to the TNM classification system announced by the American Joint Committee on Cancer (seventh edition). Hematological and pulmonary morbidities were graded using the National Cancer Institute’s Common Terminology Criteria for Adverse Events version four (CTCAE v4).

RT delivery and chemotherapy administration

All 89 eligible patients underwent CT simulation in the supine position and were immobilized with alpha cradle® kits (Smithers Medical Products). Planning CT images with 3-mm slice thickness through the entire neck, thorax, and upper abdomen were obtained. Target volume delineation was performed on the eclipse platform in which gross tumor volume (GTV) consisted of the primary tumor and lymphadenopathy as diagnosed by CT, endoscopic ultrasonography, and FDG PET/CT. Two dose levels for the clinical target volume, CTVHigh and CTVLow, were defined as follows: CTVHigh included the entire GTV with 1 cm circumferential and 3–4 cm longitudinal expanding margins whereas CTVLow covered the regional lymphatic nodal basin. Elective nodal irradiation was delivered to the supraclavicular fossa if the primary tumor or any enlarged lymphadenopathy were in the upper mediastinum, or to the celiac lymphatics if the primary tumor or any lymphadenopathy were close to the esophagogastric junction. Adding both CTVHigh and CTVLow with a 0.5 cm expanding margin in all directions was defined as PTVHigh and PTVLow, respectively. Treatment plans with associated dose maps were generated by Eclipse treatment planning system version 13 (Varian Medical Systems Inc., Palo Alto, CA, USA) or Pinnacle treatment planning system version 9.2 and 10.1 (Philips Healthcare Inc., Andover, MA, USA). The beam lines were 6-MV or 10-MV photon beams, composited by intensity-modulated radiation therapy (IMRT) or Volumetric Arc Therapy planning. The prescribed doses of PTVHigh and PTVLow were 48and 43.2 Gy divided into 24 fractions within a simultaneous integrated boost (SIB) method in 85 of 89 patients. Two patients had two-step cone down boost planning instead of SIB that delivered 45 Gy to the PTVLow and 48.6 Gy to the PTVHigh with 1.8 Gy per fraction; two patients underwent reduced planning up to 41.4 Gy with 1.8 Gy per fraction to the PTVLow area without PTVHigh boost due to concerns regarding a trachea-esophageal fistula at the previous trachea stenting sites. The coverage of 100% of the prescribed dose volume exceeded 95% of the targeted PTV while meeting the normal organ constraints and included the following: maximum dose of spinal cord was <45 Gy, V20 of lung <30% of the whole lung, and mean heart dose <30 Gy.

All 89 patients underwent concomitant chemoradiotherapy with weekly platinum-based regimens as follows: 1) cisplatin 30 mg/m2, 2) carboplatin (AUC 2), or 3) paclitaxel 50 mg/m2 plus either cisplatin 30 mg/m2 or carboplatin (AUC 2).

Defining dose-level volume on registration of PET/CT fusion image

Radioisotope activity was controlled by a Capintec CRC®-25PET Dose Calibrator. The standard was set as administration of a weight-based injection of 18F-FDG under a 60-min uptake period. PET/CT scan images were acquired on a Discovery ST (GE Healthcare, Waukesha, WI, USA) scanner. PET images within 5 mins per axial field of view were acquired and underwent post-reconstruction with Gaussian filtration.

Image registration and establishment algorithm for candidate predictors of pneumonitis

Thirty-three patients with applicable pre- and post-NACCRT 18F-FDG PET/CT images were eligible for image registration and analysis. All PET/CT scans before and after NACCRT were imported into the Eclipse treatment planning system. Planning CT images were rigidly co-registered to PET/CT images to avoid unwanted deformation of target organs. The original whole lung volume, GTV, PTVHigh, PTVLow, and dose-level volumes of the lung V10, V20, V27, V30, and V50% of the maximal prescribed dose, where dose volume Vx was defined as the total organ volume exceeding a radiation dose of “x”, were generated as regions of interest (ROIs) at planning CT and were co-registered to the corresponding relative PET/CT.

Whole lung contouring was automatically generated by the Eclipse treatment planning system (version 13) according to the CT window width between 400 and 2000 HU at the chest, followed by manual removal of the trachea and bilateral main bronchus contours. The accumulated radioactivity concentration, in Bq/c.c., of five different dose-level volume ROIs as described above were calculated using the planning system program and divided by the radioactivity concentration of the whole lung to define the SUVR of each ROI.

Percentage of delta SUVR (%ΔSUVR) is the relative change of SUVR before and after NACCRT and is described by the following equation:

Percentage of Vx, pre-NACCRT SUVR of Vx, post-NACCRT SUVR of Vx, %ΔSUVR of Vx, with x set to 10%, 20%, 27%, 30%, or 50% of the max, were evaluated as dosimetric candidate predictors. The dose-level volume at 50% of the maximum prescription dose (V50%) served as both a candidate predictor and parameter for the conformality check. The case with V50% line located out of range of 24–27Gy dose-level lines would consider as extreme value.

Follow up and diagnosis of RP

The first follow-up was performed 4–8 weeks after the completion of NACCRT and prior to surgery. After surgery, our follow-up policy was modified from the NCCN guideline (Version 2018.2) as follows: if asymptomatic, medical history and physical examination were performed every 3–6 months for 1–2 years, and then every 6–12 months for 3–5 years, followed by annual follow-up. Imaging studies and upper GI endoscopy were performed as clinically indicated. RP was diagnosed through a combination of clinical symptom and radiographic changes. In our institute, RP was considered when the patient presented with non-productive cough, dyspnea, chest pain, and consistent low-grade fever weeks to months after irradiation. Following these findings, arrangements were made for laboratory data collection, chest X-ray, and CT evaluation in order to exclude other possible causes such as infection or tumor recurrence.

Statistical analysis

Statistical analysis was performed using the Statistical Package for the Social Sciences for Windows, SPSS® software v. 18.0 (IBM Corp., New York, NY, USA; formerly SPSS Inc., Chicago, IL, USA). The comparison of patient characteristics between different grades of RP or chronic radiation lung fibrosis were evaluated using independent t-test and chi-square with omnibus test under Cox and Snell R square model. A two-tailed probability p-value of <0.05 was considered statistically significant. The association between RP and candidate predictors described above was analyzed by simple logistic regression. Multivariate logistic regression with controlling co-variants including age, clinical staging, and position of primary lesion at esophagus were subjected to test the effect of candidate predictors and pneumonitis. Receiver-operating characteristic (ROC) analysis was used to evaluate the discriminative power of pneumonitis predictors to distinguish between patients with and without acute RP ≥grade 2 and chronic radiation pulmonary fibrosis ≥grade 3. Cutoff values for % delta-SUVR pneumonitis predictors were also assessed by the Youden index (Sensitivity+Specificity–1). Kaplan–Meier analysis was used to calculate overall and progression-free survival, which were determined from the date of biopsy to the date of last follow-up.

Results

Patient characteristics

Table 1 shows the characteristics of the 89 patients with advanced thoracic esophageal squamous carcinoma treated with NACCRT. The median age at diagnosis was 62 years old (range, 40–87). The major population of the cohort was clinical stage III (n=72, 80.1%), cT3 (n=64, 71.9%), lymph node positive (n=84, 95.5%), and histological grade 2 differentiation (n=59, 66.3%). Primary tumor located at middle esophagus accounted for 53.9% of the patients (n=48). The average time interval between NACCRT completion and first posttreatment image surveillance was 6.1±1.2 weeks (mean ±STD). The median follow-up time was 32.0 months (range, 1.7–60.8 months).

Table 1 Baseline characteristics of patients with advanced thoracic esophageal cancer

Treatment compliance and outcomes

Descriptive statistics for treatment compliance are listed in Table 2. Seventy-eight of 89 (87.6%) patients completed the whole course of chemotherapy, and 82 of 89 (92.1%) patients completed the RT. Interruption of RT occurred in four patients who developed distant metastasis during treatment, two cases progressed to death during NACCRT. Treatment was discontinued in two elderly patients with intolerable fatigue, and in one patient due to severe nausea and weakness after chemotherapy. The 85.4% (76/89) overall compliance rate indicated favorable tolerance. Table 3 shows the clinical treatment outcomes. The clinical downstaging rate, defined as decreases in either T or N stage at the first follow-up after NACCRT and prior to surgery, was 59.6%. Subsequently, 65 patients (73%) were subjected to receive esophagectomy with gastrointestinal tube reconstruction. The pCR rate was 40.0% (26/65). The 1-, 3-, and 5-year overall survival rates were 62.2%, 41.5%, and 34.4%, respectively. The progression-free survival rate after 1, 3, and 5 years were 46.3%, 23.5%, and 13.9%, respectively. The median overall survival and median progression-free survival were 25.0 months (95% CI range: 9.7–40.3 months) and 9.4 months (95% CI range: 3.4–15.5 months), respectively. The local recurrence rate was 20.0% (13/65) and the distant recurrence rate was 49.2% (32/65). Kaplan–Meier survival curves of overall survival and progression-free survival are shown in Figure 1.

Table 2 Treatment and compliance of patients with advanced thoracic esophageal cancer

Table 3 Clinical outcomes

Figure 1 Kaplan–Meier survival curve of overall survival and progression-free survival. (A) Overall survival (OS) (B) progression-free survival (PFS) of all 89 reviewed patients in presented study at 1,3, 5 years relatively.

Acute and chronic pulmonary toxicity

Based on favorable compliance, acceptable pCR rate, and survival rates, toxicity profiles were analyzed and are presented in Table 4. Among these patients, serial CT scans from 33 patients having both applicable pre- and post-NACCRT images of 18F-FDG PET/CT scan were subjected to review for evaluation of acute and chronic RP. In total, 11 of 33 patients developed grade 2 acute RP. Four of 33 patients (12.2%) developed moderate to severe acute RP (≥grade 3), two of whom died from acute respiratory failure. The above two cases that resulted in mortalities developed severe pulmonary symptoms at 3–4 months after the end of NACCRT, with one case having undergone esophagectomy during this time. For the chronic pulmonary toxicity, five patients had grade 2 pulmonary fibrosis, and required monthly steroid and antitussive treatment for symptom control. Four patients with grade 3 RP developed organic consolidative patches and were hospitalized for repeated occurrence of pneumonia.

Table 4 Frequency of hematological and pulmonary toxicities

Association between clinical and physic parameters and risks of acute and chronic pulmonary toxicity

Among all clinical parameters, only age was significantly correlated with the incidence of chronic pulmonary fibrosis in univariate logistic regression (hazard ratio=1.116, 95% CI 1.008–1.236, p=0.035). Other parameters showed no significant association with lung toxicity. Both simple and multivariate logistic regression analyses demonstrated that all physic parameters, including the DVH, were not associated with statistically significant differences in the incidence of acute RP (≥grade 2) and chronic pulmonary fibrosis (≥grade 3). GTV showed a trend toward chronic radiation pulmonary fibrosis (p=0.08). The results of simple and multivariate logistic regression of physic parameters associated with lung toxicity were described in Tables 5 and 6.

Table 5 Simple logistic regression for physic and metabolic image parameters

Table 6 Multivariate logistic regression for physic and metabolic image parameters

Association between changes of metabolic response imaging and risks of acute and chronic pulmonary toxicity

ROIs generated according to dose-level volume of described physic parameters on planning CT images were rigidly co-registered to PET/CT images (representative data presented in Figure 2). Percentage of change between pre- and post-NACCRT SUVR was defined and described with %ΔSUVR. All %ΔSUVR of Vx parameters had significant correlation to acute RP in both simple and multivariate regression. For chronic radiation pulmonary fibrosis, both simple and multivariate regression analyses demonstrated a significant correlation between %ΔSUVR of V50%, V27, and V30. By contrast, separate analyses using the pre- or post-NACCRT SUVR parameters independently showed no statistically significant differences between the incidence of acute RP and chronic radiation pulmonary fibrosis in simple and multivariate logistic regression. Tables 5 and 6 display results from the correlation analyses between PET/CT dosimetric parameters and incidence of RP.

Figure 2 Representative image of patient with grade 3 acute radiation pneumonitis predicted by presented bio-physic parameters. Upper panel showed the PET-dose-volume fusion images. Lower panel showed the CT-dose-volume fusion images; blue, yellow, and pink line circled the volume of prescribed dose exceeding 20, 27, and 30 Gy relatively.

Establishment of the bio-physic constraint model for acute and chronic RP

From the statistically significant multivariant predictors shown in Table 6, the cutoff value maximizing the Youden index to predict patients with high risk of acute and chronic radiation pulmonary toxicity was calculated and is summarized in Table 7. The Youden thresholds corresponded to a cutoff value point with >90% sensitivity and >75% specificity for prediction of acute RP. Nevertheless, cutoff values for predictors of chronic radiation pulmonary fibrosis corresponded to about 90% for both sensitivity and specificity. These predictors were above AUC >0.9, which indicates good discriminative power. The cutoff value at the 6.2% increment of V20 area ΔSUVR and the 8.9% increment of V27 area ΔSUVR were selected as the new bio-physic constraint. ROC curves for the above two predictors are presented in Figure 3.

Table 7 Receiver-operator characteristic operating point of acute radiation pneumonitis and chronic pulmonary fibrosis predictors that maximized the Youden index

Figure 3 Receiver-operating characteristic (ROC) curves for predictors of the presented bio-physic constraints for acute and chronic radiation pneumonitis (A) acute radiation pneumonitis (%ΔSUVR of V20); (B) chronic pulmonary fibrosis (%ΔSUVR of V27). Abbreviations: AUC, area under the ROC curve; %ΔSUVR of Vx, percentage of SUV change between pre- and post-NACCRT at dose volume Vx.

Discussion

This study associated the changes in serial metabolic response imaging to the NACCRT-induced acute and chronic pulmonary toxicity in esophageal cancer. We demonstrated that the increment of SUVR (delta-SUVR) of PET/CT in specific dose-level volumes positively correlated to the incidence of radiation-induced acute RP and chronic lung fibrosis, whereas single time point SUVR or traditional physic parameters, used separately, did not. The cutoff values at the 6.2% increment of V20 area ΔSUVR and the 8.9% increment of V27 area ΔSUVR are powerful predictors of acute RP and chronic lung fibrosis, respectively, and could potentially represent a new bio-physic constraint model in treatment planning for NACCRT in patients with advanced thoracic esophageal squamous cell carcinoma.

In the present study, all physic parameters based on the lung DVH failed to show statistically significant correlation with pulmonary toxicity. Among the study cohorts, only four patients had V20 exceeding 30%. However, the incidence of RP developed much more frequently than expected, implying that DVH parameters alone may not be the key for predicting RP. During the past decade, radiation oncologists have evaluated thoracic treatment planning according to physic parameters based on DVH.24 The QUANTEC board serially reported advise of diametric constraints based on conclusions that the mean lung dose and V20 were best supported for routine clinical practice. The latest statement released in 2010 recommended maintaining the V20 at <30–35% and mean lung dose between 20 and 23 Gy in order to lower the risk of pneumonitis to <20%.25 However, the association between calculated dosimetry and normal organ outcome was estimated mainly using results of definitive RT and 3D-conformal RT planning, which may present difficulty in interpreting toxicity of multimodality treatment when using CCRT and IMRT.26,27 Several chemotherapy agents are known to be radiation sensitizers and therefore increase the risk of radiation reactions, including lung injury. A meta-analysis demonstrated a 1.6-fold increase in hazard ratio of RP in patients receiving CCRT compared to those receiving sequential regimen.28 When multimodality treatments consisted of RT and concurrent administration of chemotherapeutics, targeted agents, or immunotherapy, the establishment of predictors using a practical constraint model to avoid RP remains an urgent task.

As a widely used metabolic imaging method, PET/CT may provide a better tool in interpreting biological effects.29 Nonetheless, previous studies based on single time point PET scans, either during pretreatment or posttreatment, reported moderate or unsatisfactory predictive power.22,23,30 Zhang et al showed a correlation between mean lung SUV and RP (≥grade 2) in a group of patients whose radiation dose was more than 60 Gy; however, in a patient group that received <60 Gy, statistically significant correlations were not achieved.30 Also, their study had a limited power of establishing clear cutoff values for mean lung SUV in predicting RP due to a limited number of cases and limited discriminative capability of the selected parameters. Two studies based on pre-radiation PET scan images established their own predictors with significant correlation to the RP incidence with moderate sensitivity and lower specificity.22,23 Lee et al most recently successfully established strong predictors of RP using functional lung heterogenicity and dosimetry parameters using pretreatment perfusion SPECT/CT and FDG PET/CT.31 The functional lung heterogenicity describing the preexisting area of emphysematous lung indicated a potentially hypersensitive microenvironment that may be at increased risk of radiation injury. In agreement with our study, Lee et al found dose-function histogram parameters to be competent predictors of RP that were better than other physic parameters. It is important to realize that predictors solely based on pretreatment functional images may not thoroughly represent the dynamic changes in 3D anatomical features, clinical lung conditions, image textures, and adaptive RT parameters before application in clinical RT. Few studies have mentioned developing serial PET scan imaging for use as RP predictors. A recent study found a relationship between maximal SUV in irradiated lung parenchyma and RP.32 Still, the authors only focused on the discrimination of observed phenomena without quantifying the correlation using statistic tools. Also, maximal SUV may not be a precise parameter since lung parenchyma, consisting of bronchial trees and pulmonary vasculatures, has distinct inhomogeneity, which interferes with FDG uptake values of PET/CT scans. Considering that acute RP is a dynamic process with inflammatory responses in radiation fields, our bio-physic constraint model using spatial registration and delta calculation algorithms, rather than pre- or post-CCRT image sets alone, could be a practical way to establish specific constraint models for RP across institutions. This model may play a role in developing a model for validation of previous and future RP constraints. In real-world clinical practice, the prediction model indicated that medical physics should eliminate 20and 27 Gy distribution area at lung in treatment planning in esophageal cancer patients going to receive NACCRT. Our algorithm may potentially assist surgeons in evaluating pulmonary morbidity prior to surgery. Additionally, for patients who are predicted to be at high risk of pulmonary morbidity, intensive lung function monitoring and prophylactic steroid administration may be helpful. For proton therapy, the radiobiology of the lung remains unclear,33 and as such our model may be used as a reference for establishing new RP constraints for particle therapy. Moreover, predictors based on the concept of this study may better predict RP risks of combination modalities consisting of RT and immune therapy or targeted therapy, compared to physic and dosimetric parameters. The delta model could serve as a platform basis for further incorporation of data from radiomics and genomics.

Although the delta algorithm requires fewer cases, the case numbers used in our study remains insufficient and our results still require further validation. The 18F-FDG PET/CT has been routinely applied in staging workup and response evaluation for esophageal cancer patients in our institution over the last 3 years but has not yet fully covered all esophageal cancer patients. These retrospective cohorts came from the time points determined for routine use of fore-and-aft PET/CT images, raising the concern of potential selection bias.

Conclusion

By using this bio-physic constraint model, the specified %ΔSUVR of pre- and post-NACCRT PET/CT can potentially be a predictor for acute RP and chronic pulmonary fibrosis in patients with advanced esophageal cancer.

Ethics approval and informed consent

The study protocol was granted institutional review board approval (serial number: 18MMHIS194e, MacKay Memorial Hospital, Taipei, Taiwan). Written informed consent was obtained from all study participants and research was conducted in accordance with the 1975 Declaration of Helsinki, as revised in 1983.

Acknowledgment

This study was funded by MacKay Memorial Hospital (grant number: MMH-E-107-13) and the Ministry of Science and Technology, Taiwan (grant number: MOST 106-2314-B-195-002-MY3 and MOST-106-2623-E-195-001-NU). The preliminary data of this manuscript has been presented in 2017 ASTRO Annual Meeting, which was held on September 24–27, 2017, at the San Diego Convention Center, San Diego, CA. It has not been published or presented elsewhere in part or in entirety and is not under consideration by another journal. All study participants provided informed consent, and the study design was approved by the appropriate ethics review board. We have read and understood the journals policies, and we believe that neither the manuscript nor the study violates any of these.

Author contributions

Tien-Chi Hou and Yu-Jen Chen conceived, designed, discussed the study and wrote the manuscript. Tien-Chi Hou and Kun-Yao Dai did the patient follow-up and data analysis. Ming-Che Wu prepared all the PET/CT fusion and registration images and afford all the nuclear medicine technique among the study. Wen-Chien Huang did the endoscopy and biopsy during patient follow-up. Hung-Chi Tai did the calibration of radiation physic parameters and reviewed the treatment planning of radiotherapy. Kai-Long Hua afforded the algorithm of image processing. Yu-Jen Chen and Wen-Chien Huang revised and supervised the study. All authors contributed to data analysis, drafting or revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Disclosure

The authors report no conflicts of interests in this work.

References

1. Gupta B, Kumar N. Worldwide incidence, mortality and time trends for cancer of the oesophagus. Eur J Cancer Prev. 2017;26(2):107–118. doi:10.1097/CEJ.0000000000000249

2. Domper Arnal MJ, Ferrandez Arenas A, Lanas Arbeloa A. Esophageal cancer: risk factors, screening and endoscopic treatment in Western and Eastern countries. World J Gastroenterol. 2015;21(26):7933–7943. doi:10.3748/wjg.v21.i10.2937

3. Trivers KF, Sabatino SA, Stewart SL. Trends in esophageal cancer incidence by histology, United States, 1998-2003. Int J Cancer. 2008;123(6):1422–1428. doi:10.1002/ijc.v123:6

4. Szumilo J. [Epidemiology and risk factors of the esophageal squamous cell carcinoma]. Pol Merkur Lekarski. 2009;26(151):82–85.

5. Xi M, Xu C, Liao Z, et al. The impact of histology on recurrence patterns in esophageal cancer treated with definitive chemoradiotherapy. Radiother Oncol. 2017;124(2):318–324. doi:10.1016/j.radonc.2017.06.019

6. Bollschweiler E, Holscher AH. Prognosis of early esophageal cancer: differences between squamous cell carcinoma and adenocarcinoma. Ann Surg. 2007;245(2):334. doi:10.1097/01.sla.0000253072.33351.a4

7. Enzinger PC, Mayer RJ. Esophageal cancer. N Engl J Med. 2003;349(23):2241–2252. doi:10.1056/NEJMoa030969

8. Network TNCC. NCCN clinical practice guidelines in oncology (NCCN guidelines) - esophageal and esophagogastric junction cancers. 2018; Version 2. Available from: https://www.nccn.org/professionals/physician_gls/pdf/esophageal.pdf.. Accessed July 12, 2019.

9. Lordick F, Mariette C, Haustermans K, Obermannova R, Arnold D, Committee EG. Oesophageal cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2016;27(suppl 5):v50–v57.

10. Rath P, Fahrig A. [Active surveillance and postponed operation after neoadjuvant radiochemotherapy for oesophageal cancer-findings of the preSANO study]. Strahlenther Onkol. 2018. doi:10.1007/s00066-018-1371-6.

11. Shapiro J, van Lanschot JJB, Hulshof M, et al. Neoadjuvant chemoradiotherapy plus surgery versus surgery alone for oesophageal or junctional cancer (CROSS): long-term results of a randomised controlled trial. Lancet Oncol. 2015;16(9):1090–1098. doi:10.1016/S1470-2045(15)00040-6

12. Tsukada Y, Higashi T, Shimada H, Kikuchi Y, Terahara A. The use of neoadjuvant therapy for resectable locally advanced thoracic esophageal squamous cell carcinoma in an analysis of 5016 patients from 305 designated cancer care hospitals in Japan. Int J Clin Oncol. 2018;23(1):81–91. doi:10.1007/s10147-017-1178-y

13. Rizvi FH, Syed AA, Khattak S, Rizvi SS, Kazmi SA, Khan MQ. Complete pathological response after neoadjuvant treatment in locally advanced esophageal cancer predicts long term survival: a retrospective cohort study. Int J Surg. 2014;12(6):621–625. doi:10.1016/j.ijsu.2014.04.014

14. Hamai Y, Hihara J, Taomoto J, Yamakita I, Ibuki Y, Okada M. Effects of neoadjuvant chemoradiotherapy on postoperative morbidity and mortality associated with esophageal cancer. Dis Esophagus. 2015;28(4):358–364. doi:10.1111/dote.12207

15. Graves PR, Siddiqui F, Anscher MS, Movsas B. Radiation pulmonary toxicity: from mechanisms to management. Semin Radiat Oncol. 2010;20(3):201–207. doi:10.1016/j.semradonc.2010.01.010

16. Morota M, Gomi K, Kozuka T, et al. Late toxicity after definitive concurrent chemoradiotherapy for thoracic esophageal carcinoma. Int J Radiat Oncol Biol Phys. 2009;75(1):122–128. doi:10.1016/j.ijrobp.2008.10.075

17. Zhao Y, Chen L, Zhang S, et al. Predictive factors for acute radiation pneumonitis in postoperative intensity modulated radiation therapy and volumetric modulated arc therapy of esophageal cancer. Thorac Cancer. 2015;6(1):49–57. doi:10.1111/1759-7714.12142

18. Shen WB, Zhu SC, Gao HM, et al. [Low dose volume histogram analysis of the lungs in prediction of acute radiation pneumonitis in patients with esophageal cancer treated with three-dimensional conformal radiotherapy]. Zhonghua Zhong Liu Za Zhi. 2013;35(1):45–49.

19. Kumar G, Rawat S, Puri A, et al. Analysis of dose-volume parameters predicting radiation pneumonitis in patients with esophageal cancer treated with 3D-conformal radiation therapy or IMRT. Jpn J Radiol. 2012;30(1):18–24. doi:10.1007/s11604-011-0002-2

20. Song SY, Kim JH, Ryu JS, et al. FDG-PET in the prediction of pathologic response after neoadjuvant chemoradiotherapy in locally advanced, resectable esophageal cancer. Int J Radiat Oncol Biol Phys. 2005;63(4):1053–1059. doi:10.1016/j.ijrobp.2005.03.033

21. Cunliffe A, Armato SG 3rd, Castillo R, Pham N, Guerrero T, Al-Hallaq HA. Lung texture in serial thoracic computed tomography scans: correlation of radiomics-based features with radiation therapy dose and radiation pneumonitis development. Int J Radiat Oncol Biol Phys. 2015;91(5):1048–1056. doi:10.1016/j.ijrobp.2014.11.030

22. Chaudhuri AA, Binkley MS, Rigdon J, et al. Pre-treatment non-target lung FDG-PET uptake predicts symptomatic radiation pneumonitis following stereotactic ablative radiotherapy (SABR). Radiother Oncol. 2016;119(3):454–460. doi:10.1016/j.radonc.2016.05.007

23. Castillo R, Pham N, Castillo E, et al. Pre-radiation therapy fluorine 18 fluorodeoxyglucose PET helps identify patients with esophageal cancer at high risk for radiation pneumonitis. Radiology. 2015;275(3):822–831. doi:10.1148/radiol.14140457

24. Martel MK, Ten Haken RK, Hazuka MB, Turrisi AT, Fraass BA, Lichter AS. Dose-volume histogram and 3-D treatment planning evaluation of patients with pneumonitis. Int J Radiat Oncol Biol Phys. 1994;28(3):575–581. doi:10.1016/0360-3016(94)90181-3

25. Marks LB, Yorke ED, Jackson A, et al. Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys. 2010;76(3 Suppl):S10–19. doi:10.1016/j.ijrobp.2009.07.1754

26. Seppenwoolde Y, De Jaeger K, Lebesque JV. In regard to Tsujino et al.: predictive value of dose-volume histogram parameters for predicting radiation pneumonitis after concurrent chemoradiation for lung cancer. IJROBP 2003;55:110–115. Int J Radiat Oncol Biol Phys. 2003;56(4):1208–1209; author reply 1209. doi:10.1016/S0360-3016(03)00345-6

27. Tsujino K, Hirota S, Endo M, et al. Predictive value of dose-volume histogram parameters for predicting radiation pneumonitis after concurrent chemoradiation for lung cancer. Int J Radiat Oncol Biol Phys. 2003;55(1):110–115. doi:10.1016/S0360-3016(02)03807-5

28. Vogelius IR, Bentzen SM. A literature-based meta-analysis of clinical risk factors for development of radiation induced pneumonitis. Acta Oncol. 2012;51(8):975–983. doi:10.3109/0284186X.2012.718093

29. Cremonesi M, Garibaldi C, Timmerman R, et al. Interim (18)F-FDG-PET/CT during chemo-radiotherapy in the management of oesophageal cancer patients. A systematic review. Radiother Oncol. 2017;125(2):200–212. doi:10.1016/j.radonc.2017.09.022

30. Zhang Y, Yu Y, Yu J, Fu Z, Liu T, Guo S. (18)FDG PET-CT standardized uptake value for the prediction of radiation pneumonitis in patients with lung cancer receiving radiotherapy. Oncol Lett. 2015;10(5):2909–2914. doi:10.3892/ol.2015.3207

31. Lee HJ Jr., Zeng J, Vesselle HJ, Patel SA, Rengan R, Bowen SR. Correlation of functional lung heterogeneity and dosimetry to radiation pneumonitis using perfusion SPECT/CT and FDG PET/CT imaging. Int J Radiat Oncol Biol Phys. 2018;102(4):1255–1264. doi:10.1016/j.ijrobp.2018.05.051

32. Demirev AK, Kostadinova ID, Gabrovski IR. (18)F-FDG PET/CT in patients with parenchymal changes attributed to radiation pneumonitis. Mol Imaging Radionucl Ther. 2018;27(3):107–112. doi:10.4274/mirt.55706

33. Guan F, Geng C, Ma D, et al. RBE model-based biological dose optimization for proton radiobiology studies. Int J Part Ther. 2018;5(1):160–171. doi:10.14338/IJPT-18-00007.1

Creative Commons License © 2019 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.