LASSO-based NTCP model for radiation-induced temporal lobe injury developing after intensity-modulated radiotherapy of nasopharyngeal carcinoma.

We investigated the incidence of temporal lobe injury (TLI) in 132 nasopharyngeal carcinoma (NPC) patients who had undergone intensity-modulated radiotherapy (IMRT) in our hospital between March 2005 and November 2009; and identified significant dosimetric predictors of TLI development. Contrast-enhanced lesions or cysts in the temporal lobes, as detected by magnetic resonance imaging (MRI), were regarded as radiation-induced TLIs. We used the least absolute shrinkage and selection operator (LASSO) method to select Dmax (the maximum point dose) and the D1cc (the top dose delivered to a 1-mL volume) from 15 dose-volume-histogram-associated and four clinically relevant candidate factors; the Dmax and the D1cc were the most significant predictors of TLI development. We drew dose-response curves for Dmax and D1cc. The tolerance dose (TD) for the 5% and 50% probabilities of TLI development were 69.0 ± 1.6 and 82.1 ± 2.4 Gy for Dmax and 62.8 ± 2.2 and 80.9 ± 3.4 Gy for D1cc, respectively. The incidence of TLI in NPC patients after IMRT was higher than expected because the therapeutic window is narrow. High-quality longitudinal studies are needed to gain further insight into the complex spatiotemporal effects of non-uniform irradiation on TLI development in NPC patients.

Temporal lobe necrosis is one of the most serious late complications developing after definitive radiotherapy of nasopharyngeal carcinoma (NPC) patients. Such injury has a tremendous impact on both quality-of-life (QoL) and survival in the era of conventional two-dimensional radiotherapy (2D-RT)123. Currently, the use of intensity-modulated radiotherapy (IMRT) (in preference to 2D-RT or forward-planned 3D-RT) to treat NPC has greatly reduced the occurrence of late toxicities, as IMRT improves sparing of adjacent normal tissues. However, in patients with locally advanced NPC, tumor infiltration into the skull base and intracranial tissue is common. The surrounding normal organs (anteriorly: the optic nerve and chiasm; posteriorly: the brain stem) have relatively low tolerances to radiation, and strict dose constraints must be put in place to preserve their critical physiological functions. High doses to the temporal lobes (TLs), a common undesirable side-effect of treatment, may trigger RT-induced TL injury (TLI). Calculation of appropriate dose thresholds improves treatment planning, thereby mitigating TLI development and preventing the negative impacts thereof on the QoL of locally advanced NPC patients456.

In recent years, the clinical characteristics of, risk factors for, and dose-volume tolerances appropriate to prevent, TLI after IMRT have been preliminarily explored in several retrospective studies4578910. However, significant methodological heterogeneities are evident among these studies, rendering it impossible to compare the results, and the actual tolerance dose (TD) used to avoid the development of TLI after IMRT remains ambiguous. The development of TLI probably depends heavily on dose-volume-histogram (DVH)-associated risk factors. Selection of the most predictive DVH-related factor for incorporation in a risk-estimating model of TDs is technically difficult, because of the general multicollinearity problem of multivariable analysis. Recently, a sophisticated tool, the least absolute shrinkage and selection operator (LASSO), has been used to address the multicollinearity problem effectively and thus to identify important predictors of injury for logistic regression ‘normal-tissue complication probability’ (NTCP) modelling of the lung, parotid gland, small bowel, and oesophagus1112, but not the TL yet. Therefore, in the present study, we sought to upgrade the methodology to allow making valid predictions of the tolerance threshold of the TL. We used actuarial-rate-based NTCP modelling and a ‘custom’ LASSO tool to identify critical predictors of TLI development. The incidence and clinical characteristics of TLI in NPC patients treated via IMRT were also investigated. Finally, the methodological quality of the relevant literature was assessed, and we offer recommendations for future studies.

Materials and Methods

Patient characteristics

This retrospective study was approved by the ethics committee of the Nanjing Medical University Cancer Center. All participants provided written informed consent, and all experiments were performed in accordance with all relevant guidelines and regulations.

From March 2005 to November 2009, 179 newly diagnosed, biopsy-proven, consecutive NPC patients were treated in our hospital. Of these, MRI data allowing TLI evaluation after completion of IMRT were available for 132. MRI of the other 47 cases had been performed only before the earliest time of TLI occurrence (Supplementary Information, part 1). Thus, we used data from these 132 cases in the analysis. The male/female ratio was 2.7:1, and patient age ranged from 12 to 77 years (median, 47 years). Using the 2010 7th AJCC staging system13, 12 patients had stage I disease, 23 stage II, 36 stage III, and 61 stage IVa-b. By T-stage classification, 41 patients were T1, 11 T2, 33 T3, and 47 T4. All patients underwent a series of pre-treatment evaluations and examinations (including history-taking, physical examination, haematological and biochemical profiling, nasopharynx and neck MRI, thoracic-abdominal computed tomography [CT], and whole-body single-photon emission CT bone scanning) to exclude those with contra-indications to treatment and distant metastases.

Treatment Methods

IMRT

Inverse IMRT treatment planning was performed on a Varian Inspiration Platform (version 10.0) using the simultaneous integrated boost technique. The nasopharyngeal primary lesions and their direct extensions (GTVnx) and positive neck lymph nodes (GTVnd) were delineated according to the recommendations of the International Commission on Radiation Units and Measurements Reports nos. 50 and 62. The clinical target volume 1 (CTV1) was defined as the GTVnx with 5–10-mm margins to encompass areas at high risk of microscopic extension and the entire nasopharyngeal mucosa plus a 5-mm depth of sub-mucosal tissue. The CTV2 was defined by addition of 3–10-mm margins to the CTV1 to include areas at low risk of microscopic extension, the level of the identified positive lymph node, and the elective cervical region. The corresponding planning target volumes (PTVs) were generated from the GTVs or CTVs plus 3-mm margins to allow for setup uncertainties. The prescribed doses were 68–75 Gy to the PTV of the GTVnx, in 32–34 fractions; 64–75 Gy to the PTV of the GTVnd, in 32–34 fractions; 60 Gy to the PTV of CTV1, in 32 fractions; and 50 Gy to the PTV of CTV2, in 28 fractions. All patients were given one fraction daily 5 days a week. The DVHs of the organs at risk were evaluated as described in the RTOG0225 protocol, to prevent violation of the tolerance limits14.

Chemotherapy

IMRT alone was recommended for stage I patients and IMRT combined with concurrent cisplatin/nedaplatin-based chemotherapy for stage II-IVb patients15. Neoadjuvant chemotherapy was prescribed for patients with bulky lesions (at the primary site or in the neck); those with residual disease after IMRT received adjuvant chemotherapy.

Image assessment and the diagnostic criteria for TLI

The endpoint of analysis was the development of TLI as identified by MRI after irradiation. Only contrast-enhanced lesions evident on post-contrast T1-weighted images and cysts were regarded as radiation-induced temporal injuries16. All MR images were reviewed independently by two examiners (L.F.Z. and S.F.H.) who specialized in head-and-neck cancer. Consensus was reached by discussion if any initial disagreement was apparent. As both QoL and brain function were of concern, close attention was paid to clinical features relevant to TLI after imaging diagnosis.

Temporal lobe re-contouring and DVH data collection

As the temporal lobes had been delineated arbitrarily during original IMRT planning, we used a recommended method17 to re-contour the temporal lobes. This allowed us to collect accurate data on the following dosimetric parameters: the volume of the TL (TLV), the mean dose, the median dose, the Dmax (the maximum point dose), the V40 (the volume receiving at least 40 Gy; the following five parameters are similar), V50, V60, V65, V70, and V75, and the D0.1cc (the maximum dose delivered to a volume of 0.1 mL; the following four parameters are similar), D1cc, D5cc, D10cc, and D20cc. In addition, the clinical variables age, sex, T stage, and chemotherapy use were included in the analysis.

Follow-up evaluation

Follow-up included clinical assessment and MRI of the head and neck. The follow-up duration was calculated from the end of IMRT to the day of the final scan. All patients were routinely followed-up at least every 3 months during the first year, every 3–6 months during the next 2 years, and annually thereafter. The usage of MRI examination on the head and neck region during follow up was well implemented. A total of 828 MR images were collected during follow-up; an average of approximately six scans was available for each patient. The TLI latency period was measured from the time of IMRT completion to the first appearance of TLI.

Statistical Analysis

Regularization and variable selection using the elastic net (LASSO)

Two dose-response models (Cox’s regression model and a logistic model) were created to calculate the TLI TDs. The explanatory variables were selected from 15 DVH-based and 4 clinical variables. Strong correlations were evident among candidate variables, as revealed non-parametrically by calculation of Spearman’s rank coefficients, and LASSO was thus performed in a setting of 10-fold cross-validation. This regularization technique identified important predictors of TLI development. Empirical studies have suggested that the elastic net technique can outperform LASSO when the data contain highly correlated predictors. Thus, we set the alpha name-value pair to 0.5 to force LASSO to engage in elastic net regularization. The lambda value of the minimal mean squared error plus one standard deviation was the enforced setting. A free online Matlab package was used to fit the entire elastic-net path for Cox’s regression model18.

Dose response model

The dosimetric variable most predictive of TLI development, determined using the elastic net, was analysed in the context of a dose response model. Both Cox’s regression and logistic models, using only the most predictive dosimetric variable as an independent variable, were adopted to construct dose-response relationships. In Cox’s dose-effect model, curves adjusted in terms of various values of the dosimetric covariate were plotted to obtain the exact values (TDs) that rendered the failure probabilities at 60 months to be 0.05, 0.1, 0.2, and 0.5, respectively. The logistic dose-effect model was also used to calculate the TDs and to validate the results yielded by Cox’s regression model. The standard errors of the TDs at the different response probabilities were generated using maximum-likelihood estimation.

In the elastic net-fitting and dose-response models, it was assumed that each temporal lobe of a patient was independently affected. Patients who suffered local relapses and underwent repeat RT that delivered an additional dose to a temporal lobe were censored at the time of the repeat RT; this allowed us to evaluate toxicities exclusively attributed to the primary RT. Death from any cause was also regarded as a censoring event if it occurred prior to TLI diagnosis.

Statistical analysis was performed using Stata SE software, version 13.1 (College Station, TX, USA). Group patient characteristics were compared using the χ2-test and the t-test. All statistical tests were two-sided, and a P value < 0.05 was considered to reflect statistical significance.

More methodological details can be found in the Supplementary Information.

Results

Survival and patterns of treatment failure

The median follow-up duration was 63.5 months (range, 11–106 months). The 5-year overall survival, disease-free survival, and disease-specific survival rates were 84.7%, 74.0%, and 86.9%, respectively. Fifteen (11.4%) patients developed recurrence and 24 (18.2%) developed distant metastases. The median time to recurrence was 25 (range, 11–72) months and that to development of distant metastasis 15.5 (range, 3–90) months. Twenty-three (17.4%) patients died during the follow-up period.

Incidence and characteristics of TLI

A total of 17 cases (17/132, 12.9%, crude rate) developed MRI-indicated TLI. The baseline characteristics of these patients are summarized and compared with those of patients without TLI in Table 1. Of the 17 cases, 13 (76.5%) and 4 (23.5%) had radiation-induced injuries to the unilateral and bilateral TLs, respectively. Two cases simultaneously developed radiation-induced brain stem damage (contrast-enhanced lesions were evident on post-contrast T1-weighted images). Of the 17 patients who developed TLIs, 8 exhibited varying degrees of clinical symptoms, such as vertigo, headache, memory deterioration, muscle weakness, and personality changes; the remaining 9 patients were asymptomatic. Overall, the MRI findings in the TLI patients revealed a continuous spectrum of RT-associated damage. Small solid enhanced nodules were evident in five cases, while moderate and large lesions were apparent in seven and five cases, respectively (lesions were classified as described in ref. 16). The outcomes of TLI evolution were documented; the MRI features of TLI (including signal abnormality on T1 and T2 images) resolved completely in four patients following gradual improvement over a few years. In contrast, one patient developed a large contrast-enhanced lesion, with a central necrotic core, in the right TL, accompanied by severe white matter lesions with a large confluent area that extended superiorly and caused a mass effect. This was the only patient to die from TLI. Apart from one T2 case who had not received chemotherapy, all other TLI patients were of stage T3 (6 cases) or T4 (10 cases) and underwent chemotherapy during their treatment periods. Data from a representative TLI case, and the corresponding dose distributions, are shown in Fig. 1. The median TLI latency period was 43 months (range, 19–68 months). The overall TLI-free survival rate was calculated using the Kaplan-Meier procedure (Fig. 2); the actuarial TLI-free survival rates were 95.5% and 85.8% at 3 and 5 years, respectively.

Table 1

The base-line characteristics of the patients with and without temporal lobe injury (TLI).

Characteristic	Group without TLI n = 115 (%)	Group with TLI n = 17 (%)	P-value
Age (years)			0.15
Median	46	50
Range	12–77	35–65
Gender			1.0
Men	83 (72.2)	13 (76.5)
Women	32 (27.8)	4 (23.5)
Stage			0.04
I-II	34 (29.6)	1 (5.9)
III- IVa-b	81 (70.4)	16 (94.1)
T category			0.002
T1-2	51 (44.4)	1 (5.9)
T3-4	64 (55.7)	16 (94.1)
Chemotherapy			1.0
Yes	98 (85.2)	15 (88.2)
No	17 (14.8)	2 (11.8)

Figure 1

An example of TLI lesion showing (a) homogeneous high signal area in the white matter on T2-weighted images (white arrow) and the medial necrosis (black arrow) (b) contrast-enhanced lesion on postcontrast T1-weighted MR images with internal necrosis (black arrow) (c) corresponding isodose lines.

Figure 2

The probability of the TLI-free survival (n = 132) is shown.

Each censored case is indicated by a spike.

Selection of predictors and estimation of tolerance doses

The maximum, minimum, and mean values of 15 DVH-associated variables are listed in Table 2 by subgroup (with or without TLI). Elastic net fitting of these 19 variables revealed that only 4 (Dmax, D0.1cc, D1cc, and T stage) were selected (Figs s1 and s2). Based on prior knowledge451019 and the law of parsimony, we preferred to use Dmax and D1cc to estimate TLI risks employing single DVH point thresholds in univariate logistic/Cox’s dose-response models.

Table 2

Characteristics of the dose-volume variables for the two subgroups (n = 264).

Variable	Variable value for patients without TLI			Variable value for patients with TLI
	Mean	Minimum	Maximum	Mean	Minimum	Maximum
Dmax	70.1	56.1	80.8	74.3	65.3	83.1
Mean dose	22.7	6.9	47.2	27.3	19.1	37.7
D0.1cc	68.3	54.5	79.1	72.7	63.9	82.5
D1cc	64.5	48.1	77.4	69.9	58.8	81.1
D5cc	53.3	28.1	75.9	62.0	42.9	77.2
D10cc	43.4	14.6	74.8	52.9	32	73.3
D20cc	32.6	8.2	71.9	39.8	24.2	55.3
TLV	100.2	67	137.8	99.5	70.3	122.2
V75	0.2	0	9.0	1.0	0	8.4
V70	1.0	0	23.6	3.0	0	12.2
V65	2.3	0	30.6	5.1	0.003	14.9
V60	4.0	0	36.5	7.3	0.8	17.5
V50	8.0	0.8	48.4	12.2	3.1	23.2
V40	13.6	2.1	60.1	19.9	6.0	33.5
Median dose	18.9	3.3	45.9	22.4	13.6	31.3

Cox’s model

A global test using the Schoenfeld residuals indicated that the PH assumption was not violated. Univariate Cox’s regression showed that both Dmax and D1cc were significant risk factors (P < 0.001). Thus, we evaluated Dmax and D1cc individually in Cox’s dose-response model. The 5-year TLI incidences were predicted to be 5, 10, 20, and 50% for Dmax values of 69.6, 72.7, 76.0, and 81.0 Gy, respectively. For those D1cc values, the doses were 63.4, 67.9, 72.5, and 79.6 Gy, respectively (Fig. 3).

Figure 3

Failure functions adjusted to Dmax using a univariate Cox model (n = 264) are shown.

The failure functions adjusted to D1cc were similar and omitted.

Logistic modelling

Sun et al.8 found that 19/20 (95%) TLI cases had latencies of less than 50 months, and in the present study, 16/17 (94%) of the TLI cases had latencies of less than 50 months. Thus, injury-free TLs evident on MRI for more than 50 months were regarded as normal. Cases censored because of death or repeat RT within 50 months were excluded because of the need for adequate follow-up to determine whether the TLs would develop radiation injuries.

Consequently, the 21 cases censored within 50 months because of death or repeat RT were excluded, 42 TLs were excluded, and 222 TLs were subjected to analysis. Using the two input variables selected from the elastic net, the results of univariate logistic regression were consistent with the findings of Cox’s regression, namely, that both Dmax and D1cc were significant predictors of injury. Thus, either Dmax or D1cc was adopted as the best parsimonious predictor in logistic dose-response modelling. The TDs at different response probabilities are shown in Table 3 and are consistent with the dose constraints acquired from Cox’s model. To allow visualization, the fitted dose-response curves and the superior and inferior standard error bounds are plotted together with the crude incidence rates calculated at discretized dose intervals of 2.5 Gy (Fig. 4).

Table 3

Tolerance dosages for the specified effect probability levels*.

Tolerance doses	TD for D1cc	TD for Dmax
TD1	52.7 ± 4.5	61.7 ± 3.2
TD5	62.8 ± 2.2	69.0 ± 1.6
TD10	67.4 ± 1.5	72.4 ± 1.1
TD20	72.4 ± 1.7	76.0 ± 1.2
TD30	75.7 ± 2.3	78.4 ± 1.6
TD50	80.9 ± 3.4	82.1 ± 2.4

*Note that the tolerance doses were referred to with an underlying background of average of 33 fractions in the model.

Figure 4

Dose effect curves using a logistic model (n = 222) (a) for variable Dmax and (b) variable D1cc. Dotted lines indicate the single standard error of the tolerance dose at discrete probabilities of 0.01 to 0.99.

Discussion

Although IMRT has been widely used to treat NPC over the past decade, neither the incidence of TLI after IMRT nor the TD have been well-defined, and accurate data are required for clinical planning. Recently, in two early studies of IMRT-induced TLI in NPC patients treated at a single centre, the authors discussed their experiences with dosimetric features associated with the occurrence of TLI47. However, non-ignorable selection biases were evident in both studies; the inclusion criteria for injured and normal TLs were not balanced (normal TL, no injury on follow-up for ≥60 months; injured TL, injury present for as long as injury was noted). In essence, the use of crude rates can underestimate the risks, causing deviations in TD calculations, attributable to inaccurate TLI rates20. Therefore, we recommend the use of actuarial-rate-based approaches for TD estimation. This indicates that Cox’s model must be adjusted in terms of the specified independent variable used to calculate tolerance. Such methods may potentially be useful additions to typical NTCP modelling. Our work may represent an initial step in the development of a statistically rigorous method for estimation of TDs.

Radiation dose-response curves are typically sigmoid in shape, and NTCP models such as the logistic and Lyman models are frequently used to construct such curves21. Empirically, the sigmoid-shaped NTCP model applied in our present study and in ref. 10 is superior to the linear model used in refs 4 and 7 for data-fitting and may therefore provide more precise estimations. It is also important to consider multicollinearity effects, which were not addressed in most studies featuring multivariable analyses of dosimetric risk factors for TLI45789. As most dosimetric predictors are highly correlated, direct insertion of DVH-based parameters into a multivariate model may not achieve meaningful results. The most recent study used principal component analysis (PCA) to select variables (in a pre-processing step) and employed a logistic NTCP model similar to ours to calculate TDs10. However, as we saw in ref. 10, the extensive multicollinearity evident among dosimetric parameters was avoided by manually selecting D1cc/VD>80 to represent all other correlated factors, rather than the efficiency of PCA per se. However, it is questionable whether the D1cc/VD>80 is representative of all information in the cluster. Thus, as PCA is intrinsically defective, we applied elastic net fitting to Cox’s regression model to identify the crucial factors. The methodological limitations of recent reports on TL tolerance doses in NPC patients treated via IMRT are listed in Table 4. Considerable methodological heterogeneities are evident among the studies, and therefore, it is very difficult to compare the results. Harmonization of the methodology using actuarial-rate-based NTCP modelling and employing effective predictor-selection tools may be warranted in future studies of TD calculation.

Table 4

The methodological quality of the published literature regarding the tolerance dose of temporal lobe in NPC patients after IMRT.

	Consecutive cohort	Multicolinearity effect	Actuarial rate-based	Selection bias	NTCP model
Ref. 4	Yes	Direct selection (Dmax/D1cc) no MA	No	Yes	Linear
Ref. 7	Yes	Not considered in MA	No	Yes	Linear
Ref. 5	Yes	Not considered in MA	No	Unknown	Not used
Ref. 8	No (only TLI cases)	Not considered in MA	No	—	Not used
Ref. 9	No (1:1 match)	Perform MA for each DVH parameter	No	Yes, serious, maybe	Not used
Ref. 10	Yes	using principal component analysis	No	Yes	Logistic
This study	Yes	using Elastic Net(LASSO)	Yes	No	Logistic/cox

Abbreviations: MA = multivariable analysis; NTCP = normal tissue complication probability; TLI = temporal lobe injury; DVH = dose volume histogram; LASSO = the least absolute shrinkage and selection operator.

In our current series, the crude TLI rate after IMRT was 12.9% (17/132), higher than the values in recent reports (~7.5%)56. The principal reason is that the entire consecutive cohort served as the denominator; this included some patients for whom MRI data were inadequate for use. This also partially explains the relatively higher biologically effective dose to the TL (e.g., Dmax: 70.5 Gy in 33 fractions, on average) noted in the present study. On the other hand, the rate increase may be partially attributed to the fact that the MRI data is relative sufficient during follow-up in our present study. The TLI of the four cases in whom MRI abnormalities disappeared over a few years may have been overlooked if the intervals between MRI scans had been excessively long. In addition, heterogeneities in terms of patient characteristics, event definitions, prescribed total doses and fractionation doses for tumor treatment22, target delineation strategies, TL constraints in the original plans, and chemotherapy all affect the crude and actuarial incidences to some extent. A focus on the dosimetric parameters of the TL per se avoids the need to consider some of these factors, such as the prescribed total and fractionation doses. The dosimetric parameters of the TL per se are better predictors of TLI than are tumor parameters.

Furthermore, we have shown that the NPC therapeutic window is extremely narrow. The dose constraints for Dmax in Table 3 clearly indicate that the risk of TLI development increases from 1 to 5% with a TD increment of 7.3 Gy and from 5 to 10% with a dose increment of only 3.4 Gy. Thus, an additional boost of 5 Gy in two fractions towards the end of a basic RT course of 70 Gy with 2-Gy daily fractions may dramatically exacerbate the TLI risk23. The “narrow therapeutic window” is well-explained by the TL NTCP curve (Fig. 4a). At approximately 70 Gy, the slope increases significantly within a relatively small interval. As a result, a planner should be prudent to provide any additional boost to the primary tumor after the TL Dmax has reached approximately 70 Gy. If it is difficult to suppress the maximum dose to less than 70 Gy in patients with locally advanced disease, then the TL D1cc or the volume receiving the high dose9 should be maximally reduced to minimize the extent of TLI, if it would occur.

Our study had several limitations, including its retrospective study design, heterogeneities in chemotherapy regimens, and a relatively small sample size and event number; these may have limited the statistical power. Another problem was the ad-hoc nature of follow-up MRI, which was performed primarily to detect possible locoregional recurrences rather than TLI. Notably, 9/17 patients with TLIs evident on MRI were asymptomatic. Unless all included patients undergo (at least) annual MRI, the TLI rate will remain under-reported, and dosimetric calculations may thus be in error. Such deficiencies are evident in all relevant retrospective articles, although we made an effort to increase the MRI follow-up frequency in the present study. Ideally, a prospective multicentre study should confirm our data before any definite conclusions are drawn. Given the long latency and relative rarity of overt TLI, the best case scenario, at present, would be retrospective analysis of data from large cohorts of NPC patients treated both prospectively and uniformly in clinical trials. However, our report is useful, because it identifies a valuable direction for future research.

Additional Information

How to cite this article: Kong, C. et al. LASSO-based NTCP model for radiation-induced temporal lobe injury developing after intensity-modulated radiotherapy of nasopharyngeal carcinoma. Sci. Rep. 6, 26378; doi: 10.1038/srep26378 (2016).