Next Article in Journal
Current Status of CRISPR/Cas9 Application in Clinical Cancer Research: Opportunities and Challenges
Next Article in Special Issue
A Novel and Automated Approach to Classify Radiation Induced Lung Tissue Damage on CT Scans
Previous Article in Journal
The Bright and the Dark Side of TGF-β Signaling in Hepatocellular Carcinoma: Mechanisms, Dysregulation, and Therapeutic Implications
Previous Article in Special Issue
Systematic Review of Single-Fraction Stereotactic Body Radiation Therapy for Early Stage Non-Small-Cell Lung Cancer and Lung Oligometastases: How to Stop Worrying and Love One and Done
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Analysis of Radiation-Associated Parenchymal Lung Change

1
Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, UK
2
UCL Cancer Institute, University College London, London WC1E 6BT, UK
3
Sussex Cancer Centre, Royal Sussex County Hospital, Brighton BN2 5BE, UK
4
UCL Respiratory Department, University College London Hospital, London NW1 2PG, UK
5
Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, UK
6
Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool L69 3GE, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Cancers 2022, 14(4), 946; https://doi.org/10.3390/cancers14040946
Submission received: 24 December 2021 / Revised: 7 February 2022 / Accepted: 8 February 2022 / Published: 14 February 2022

Abstract

:

Simple Summary

Radiotherapy is commonly used to treat inoperable locally advanced lung cancer. Despite the use of sophisticated modern planning and imaging techniques to target the tumour and minimise dose to normal lung tissue, patients can suffer from acute and chronic respiratory problems after treatment. Currently, our understanding of the impact that radiotherapy has on patients’ lungs is inadequate. We have, therefore, proposed a novel classification of the damage to the lung tissue, as seen on CT scans after a course of radiotherapy to a lung tumour. We have used deep learning algorithms to allow large numbers of CT scans to be labelled at the level of the individual voxel according to the degree of damage. The dose delivered to the tumour and the change in lung function of the patient after treatment both correlated well to the degree of radiological change measured. Our novel, automated classification combined with a dedicated image registration method has demonstrated an important clinical application that could be used to improve radiotherapy delivery in the future by allowing us to precisely track the changes seen after radiation treatment.

Abstract

We present a novel classification system of the parenchymal features of radiation-induced lung damage (RILD). We developed a deep learning network to automate the delineation of five classes of parenchymal textures. We quantify the volumetric change in classes after radiotherapy in order to allow detailed, quantitative descriptions of the evolution of lung parenchyma up to 24 months after RT, and correlate these with radiotherapy dose and respiratory outcomes. Diagnostic CTs were available pre-RT, and at 3, 6, 12 and 24 months post-RT, for 46 subjects enrolled in a clinical trial of chemoradiotherapy for non-small cell lung cancer. All 230 CT scans were segmented using our network. The five parenchymal classes showed distinct temporal patterns. Moderate correlation was seen between change in tissue class volume and clinical and dosimetric parameters, e.g., the Pearson correlation coefficient was ≤0.49 between V30 and change in Class 2, and was 0.39 between change in Class 1 and decline in FVC. The effect of the local dose on tissue class revealed a strong dose-dependent relationship. Respiratory function measured by spirometry and MRC dyspnoea scores after radiotherapy correlated with the measured radiological RILD. We demonstrate the potential of using our approach to analyse and understand the morphological and functional evolution of RILD in greater detail than previously possible.

1. Introduction

Lung cancer is the leading cause of cancer-related death worldwide [1]. In inoperable locally advanced non-small cell lung cancer (NSCLC), chemo-radiation (CRT) is the standard treatment. Eligible patients are also offered up to 12 months of a Programmed death-ligand 1 inhibitor [2]. The lung is one of the most sensitive tissues to ionising radiation, and this limits the dose that can be delivered to lung tumours [3]. Radiation-induced lung damage (RILD) is divided into an acute, reversible phase, radiation pneumonitis, and a late, permanent, fibrotic phase [4].
Radiation damage is caused by both direct nuclear and mitochondrial DNA damage, and the generation of free radicals. Within minutes of irradiation, altered gene expression can be detected and growth factors, such as TGF-ß, PDGF and interleukin 1, are released [3]. Two distinct pathophysiological mechanisms of RILD have been described [4]. Classic RP occurs in-field and follows radiation to large volumes of lung parenchyma. A three-phase model of the histological changes (an early, intermediate and late phase) is described. The early or latent phase is visible only with electron microscopy, which reveals degenerative changes to type I and type II pneumocytes, thickened secretions of mucus from goblet cells, basement membrane swelling and changes to endothelial cells [4,5]. This early phase is dose dependent and does not generally occur at doses <10 Gy [5]. Damaged cells lead to cytokine release, e.g., TNF alpha, which leads to an acute inflammatory, exudative pneumonitis: the intermediate phase. The final, fibrotic phase is a result of pathological repair of the acute inflammatory insults. Fibroblasts produce collagen after stimulation by the acute phase proteins, particularly TGF-beta. This leads to a lack of lung elasticity, volume shrinkage and scarring. There is accompanying vascular damage, a decrease in type I pneumocytes and a return of type II pneumocytes.
This three-stage model, which presents a dose-dependent response to radiation, though widely cited, is a simplification. RILD is modulated by a host of genetic, environmental and psychological variables that are not yet well understood. An individual patient’s response to thoracic irradiation is somewhat unpredictable and sometimes out of proportion to the dose or volume of lung irradiated. Furthermore, radiological appearances and clinical symptoms are often poorly correlated [6]. It has been suggested that the three-stage model, derived largely from animal studies, or human studies where single dose, whole lung irradiation is administered, may not be clinically applicable to fractionated radiotherapy [5].
The incidence of symptomatic RILD is estimated to be 15–40% with a mortality rate of <2–4% [7,8]. Clinically, radiation pneumonitis manifests within weeks to months after the completion of radiotherapy. It presents with the classic triad of dyspnoea, non-productive cough and hypoxaemia. Pulmonary function tests usually reveal a restrictive defect due to volume loss. Gas exchange is also impaired and a fall in diffusion capacity may be seen.
Patient, tumour and treatment factors affect the likelihood of developing RILD [9] and its severity. Chemotherapy can have a synergistic effect with radiotherapy to exacerbate RILD [10]. Immuno-oncology (IO) therapy can cause pneumonitis, which is a widely recognised complication [11]. Rates of interstitial lung disease with IO agents are reported at rates between 2.7% and 3.5% [10]. It is likely that by priming the immune system, RT exacerbates this [12,13]. A secondary analysis of KEYNOTE-001 demonstrated that previous RT increased the incidence of pulmonary toxicity in patients receiving pembrolizumab [13].
A dose-response relationship can be plotted between Mean Lung Dose (MLD) and probability of radiation pneumonitis (or RP). Various Vx values (% of lung volume receiving ≥ X Gy) are associated with RP risk. There is no sharp dose threshold below which there is no risk of RP. Different dosimetric parameters are closely correlated within individual datasets, suggesting that there is no optimum threshold [14]. Widely accepted clinically-useful normal lung tissue tolerance doses allow the reasonably effective and safe deployment of radiation for lung tumours [14]. However, dose escalation remains an important ambition for radiation oncologists with the hope of improving survival [15]. The improved outcomes achieved for this patient cohort mean that clinicians must now pay more attention to the morbidity of treatment which patients can expect to endure for a greater period of time than their historical counterparts due to prolonged overall survival [16].
Despite a number of published classification systems [17,18,19,20,21], there are not clinically useful RILD classification to allow detailed, quantitative analysis of the radiological changes seen after RT, or that have known correlations with clinical or dosimetric measures. In order to address this, in 2018, our group published a novel method to quantify RILD with twelve automated image-based biomarkers to allow precise, objective and continuous measures of changes on CT scans following RT [22]. The biomarkers track changes to the size and shape of the lungs, the parenchyma and pleural thickness. The biomarkers have been used to track the evolution of RILD over time and were able to differentiate patients into two distinct subtypes according to the temporal pattern of radiological RILD manifested [23].
However, the biomarkers characterised parenchymal change by simply thresholding the lung tissue into ‘normal lung’ and ‘consolidation’ based on the HU value. This approach was an over-simplification and was unable to distinguish the different morphological subtypes of parenchymal change that could be seen in the images. Furthermore, the dosimetric correlations presented in [23] used only global DVH metrics and did not study the relationships between local dose and parenchymal change as, at that time, it was not possible to perform meaningful lung registrations due to the magnitude of the geometrical changes resulting from long-term RILD.
Our group has recently been undertaking work to address the issues raised above. Firstly, we developed a novel deformable registration method that enabled us to register the follow-up CT scans and the baseline scan despite the extreme geometrical and anatomical changes secondary to RILD occurring between the images [24]. This technique, rather than using the image intensities to guide the registrations, instead uses salient features between the images such as the lung boundary, the major airways and the blood vessels, to successfully register CT scans of the same patient up to 24 months apart. Secondly, we have developed a novel automated classification of parenchymal texture at the voxel-wise level, which enables us to study changes to parenchymal tissue due to RILD. To achieve this, we trained a deep-learning network using a two-stage approach akin to active learning. The technical details of this methodology are presented in an accompanying paper in this special issue [25].
The work presented in this paper utilises the recent advances outlined above and applied them to a clinical dataset to quantify parenchymal change due to RILD at a level of detail not previously performed. This includes longitudinal analysis of the proportion of the different tissue classes in the lungs, analysis of the relationships between the tissue classes at each time point and global dose metrics and measures of respiratory function, and investigating the relationship between local dose and changes to the parenchymal tissue over time.

2. Materials and Methods

IDEAL-CRT [26] was a multicentre, UK, phase 1/2, single arm, dose escalation trial in non-small cell lung cancer. Patients received isotoxically individualised tumour radiation doses of 63 to 71 Gy in 5 weeks or 63 to 73 Gy in 6 weeks, delivered concurrently with two cycles of cisplatin and vinorelbine. A total of 118 patients (Stage IIB-IV) were recruited from 9 UK centres. All patients gave informed, written consent. Scans of 46 patients who had CT scans at all follow-up points (3, 6, 12 and 24 months post RT) were available. For the baseline scans, a diagnostic CT obtained just prior to RT was used rather than the RT planning scan, as these were generally of superior quality. The RT planning scan, RT structures and dose distribution data were also available for all patients. MRC dyspnoea scores and spirometry data at baseline and follow-up were available for most patients (see Table 1 and Table 2 for full details).
We have proposed a novel classification of lung parenchyma texture at the voxel-wise level and developed a deep learning-based method of its automated annotation. Following a detailed study of the texture changes seen in our patient cohort, and a review of established classifications in the literature [17,18,20,21,27,28,29,30,31,32], we devised a five-class classification based on the morphology and texture of the lung parenchyma. The full details are presented in our accompanying paper in this Special Issue [25]. Table 3 lists the different classes and Figure 1 shows examples of the classes on several CT slices.
To automate the classification, the data were divided into a development set (40 patients) and a hold-out test set (6 patients). An initial manual labelling of the development CT scans was performed by a radiation oncologist (E.C.) with review and refinement by a thoracic radiologist (J.J.). As the manual labelling was challenging, and involved some uncertainty in the labelling process, we decided to adopt a two-stage learning approach, akin to active learning. We used the initial manual labels on the development dataset to train an initial ensemble of Convolutional Neural Networks (CNNs), of which the results were combined to form an initial automatic labelling of the development data. The initial manual and automatic labels were carefully reviewed and edited to produce a single revised set of ground-truth labels for the second stage of training and evaluation. An expanded ensemble of CNNs were then trained on the revised labels, and the results were evaluated on the hold-out test set. The ground truth labels on the hold-out test set were generated entirely manually, but only after the revised ground truth labels on the development set had been generated to ensure they were consistent with the revised labels.
In order to propagate the Planning Tumour Volume (PTV) and local dose distribution onto the baseline and follow-up CT scans, the planning RT scans were registered with each of the five available diagnostic CT scans obtained before and following treatment; see Figure 2 for an example. This was performed using our recently developed novel deformable registration method, which aligns images based on the salient features between the images. Briefly, this method does not use the intensity information in the images to align the images, but instead uses segmentations of the lungs and major airways, together with blood vessels detected using a ‘vesselness’ filter [33]. These features are then used to drive a multi-channel deformable registration. The lung segmentations were performed manually and the airways segmentations were performed automatically using the open-source Pulmonary Toolkit [34] and subsequently reviewed and manually edited if needed using ITK-SNAP [35]. The ‘vesselness’ [33] of each scan was calculated using the Pulmonary Toolkit. The multi-channel registrations were performed using the opensource NiftyReg software [36]. For each patient, the planning scans were registered to all of the diagnostic CT scans, and the registration results were used to propagate the PTV and local dose distribution onto the scans at all time points. All registrations were manually reviewed and judged to be of sufficient quality, except for one patient, who could not be registered due to the severity of lung damage, meaning airways and vessels could not be aligned. Full details have been published elsewhere [24].
To analyse the prevalence of each parenchymal class at the different timepoints, the relative volume (as a percentage of the total lung volume) of each parenchymal class at every time point was calculated for each individual, and the mean values over all 46 patients were calculated and visualised as pie charts.
To analyse the evolution of the parenchymal classes over time, the difference in the relative volumes of the parenchymal classes was calculated between each follow-up CT scan and the pre-RT scan, for each individual. For this analysis, the voxels within the PTV were excluded from all timepoints so that presence of the tumour in the baseline scan did not influence the analysis and only the effects on lung tissue were studied. The registration results are required to transform the PTV on to each of the follow-up CT scans, so the one patient for whom the registrations failed was excluded from this (and all further) analysis. The differences in the relative volumes of the parenchymal classes, and their distribution over all individuals, were visualised using boxplots. A Friedman test was used to test for statistical differences between the distributions at different time points.
The relationship between the global dose and the change in the prevalence of the different tissue classes at each timepoints was analysed. For each follow-up CT scan, the Pearson correlation was calculated between the relative volume change of each tissue class and the following dosimetric parameters: V5, V10, …, V60, Mean Lung Dose (MLD), Max Lung Dose, the GTV volume and the PTV volume.
The relationships between the changes in the tissue classes and the changes in lung function were also investigated. Four different respiratory parameters related to lung function, FVC, FEV1, TLCO and MRC score, were measured at most timepoints for most patients (Table 1b). The changes in the respiratory parameters between each follow-up timepoint and the baseline measurements were calculated, and the Pearson correlation was calculated between the changes in respiratory parameters and the changes in tissue classes.
The relationship between the local dose distribution and the tissue classes distribution, and how this changed over the timepoints, was also explored. Each diagnostic CT scan was divided into subvolumes based on the propagated dose distributions, with a subvolume for each 5 Gy physical dose band up to 65 Gy, and then a final subvolume for all lung voxels receiving more than 65 Gy. For each patient and timepoint, the relative volumes of the tissue classes within each subvolume were calculated, and the average values were calculated for all patients. These results were visualised and qualitatively analysed using a stacked bar graph.
All statistical analysis was performed using Microsoft Excel for Mac Version 16.16.27, RStudio Version 1.2.1335 and IBM SPSS Statistics Version 27.

3. Results

The characteristics of the 46 patients included in this study are summarised in Table 1. The parenchymal tissue classes developed in our accompanying paper [25] are described in Table 3 and examples are given in Figure 1. The classes (1–5) were designed to represent texture with increasing density.
Figure 3 shows the mean proportion of each parenchymal class over all 46 patients at each timepoint. It can be seen that the prevalence of the different classes is broadly similar between the different timepoints, although the exact proportion of each class does vary between the timepoints. Class 1 is most prevalent for all timepoints and ranges from 90.6% at 6 months to 93.5% at the pre-RT scan. The next most prevalent class is Class 2, which ranges from 4.0–6.5%, followed by Class 5, which ranges from 1.4% to 2.9%. Note, as the PTV voxels were not excluded from this analysis, most of the Class 5 tissue in the pre-RT scan corresponds to the tumour. Classes 3 and 4 both represent very small proportions of the lungs, ranging from 0.25–0.74% and 0.29–1.25%, respectively.
Figure 4 shows a boxplot for each parenchymal class, presenting the difference in the relative volume of that class between each of the follow-up CT scans and the pre-RT scan for the 45 patients that were successfully registered. Each of the classes displays its own distinctive temporal pattern. Prevalence of Class 1, which best represents undamaged tissue, shows a consistent decrease compared to pre-RT. This is most pronounced at 6 months and then gradually returns towards baseline at later timepoints. Classes 2 and 3, which represent ground glass and are, therefore, radiological markers of pneumonitis, are most prevalent at 3–6 months and then return towards baseline values at 24 months. Classes 4 and 5, which represent more solid textured parenchyma, behave differently from each other. Class 4 peaks at 6 months, while Class 5 continues to increase in prevalence up to 24 months. The latter correlates with the increased incidence of late lung fibrosis and collapse. The former may be an intermediary stage between the more acute Classes (2 and 3) and Class 5, representing evolving fibrosis. Friedman tests showed that there were statistically significant differences between the volumes at the time points connected by a horizontal line.
Figure 5 shows a Pearson correlation matrix for a number of global dose metrics and the relative volume changes of the different tissue classes at 3 months. Mild to moderate correlations are seen (up to r = 0.5). As expected, volume of Class 1 is negatively correlated with dose, while the other tissue classes show a positive correlation. The strongest correlation between the tissue classes is a negative correlation (r = −0.91) between the volumes of Class 1 and Class 2, suggesting that the decrease in default lung tissue at 3 months is mostly a result of an increase in ground glass texture (a classical radiological hallmark of radiation pneumonitis). The dosimetric variables that most strongly correlated with change in tissue class volume are lung V20–V50 Gy and mean lung dose. Lung V30 Gy was correlated with change in Classes 1, 2, 3 and 4 with r values of −0.45, 0.5, 0.31 and 0.42, respectively. Change in Class 5 shows the least strong correlations with global dose. Similar patterns are seen at other follow-up time points (see Figures S1–S3).
Figure 6 shows a Pearson correlation matrix presenting the relationship between change in respiratory outcome metrics between baseline and 24 months and the changes in volume of tissue classes at 24 months. A moderate correlation was seen between change in Classes 1 and 2, and changes in the FVC, DCLO and MRC dyspnoea score. A weak correlation was observed between changes in volume of Classes 4 and 5 and changes in the MRC dyspnoea score. However, at earlier time points, the correlations were weaker. See Figures S4–S6 for the other time points.
Figure 7 represents the relationship between the local dose and the distribution of the tissue classes, and how this changed over the different timepoints. It can be seen that the prevalence of tissue Classes 2–5 increases more in regions of higher dose, and remains fairly constant in regions of low dose, indicating that the changes to the tissue classes are indeed caused by radiation. It can also be seen that the dose response and temporal evolution is different for the different tissue classes. There is an initial increase in Class 2 at 3 months, which is approximately linearly related to local dose, and then a gradual decrease over later timepoints. There is also an increase in Class 3 at 3 months, but this occurs more evenly over the mid- and high-dose regions, and then rapidly decreases at later time points. Likewise, Class 4 increases at 3 months in the mid- and high-dose regions, but interestingly, decreases in the mid-dose regions but continues to increase in the high-dose regions at 6 months, and then gradually decreases in all regions at 12 and 24 months. Class 5 increases over all time points, but at 3 months, the increase is approximately the same in the mid- and high-dose regions, whereas at later time points, there is an increasingly linear relationship between the increase in Class 5 and the local dose.

4. Discussion

The relationships between radiotherapy dose, radiological changes and clinical outcomes are complicated. There are a large number of important clinical variables that make its analysis complex. Nonetheless, in this work, we have demonstrated that our novel classification system of lung parenchyma is a useful tool for studying RILD and uncovering patterns in its evolution and relationship to dose. Our classification system tracks longitudinal change that are related to both global and local prescribed RT dose. The global metrics with the strongest correlations were seen at 3 months and include r values of up to 0.5 between Class 2 and lung V30 Gy. Classes 1–4 all showed consistent mild to moderate correlation with a number of global dose metrics. Clinical outcomes (spirometry and MRC dyspnoea score) were available for most patients but only showed weak correlations with parenchymal changes. This is likely due to the number of confounders that can affect these measures. Spirometry is a summary measure of lung physiology and is determined by co-morbidities and patient compliance, amongst other things. Signals from spirometry that are directly affected by RILD can easily be lost amongst the noise of other variables [37]. The MRC dyspnoea score has even greater problems due to its subjective nature, and is significantly influenced by cardiac, orthopaedic and psychological morbidity as well as RILD and other respiratory conditions.
Given the complexity discussed above, this work was motivated by the recognition that the relationship between radiotherapy dose, radiological changes and clinical outcome is not yet adequately understood. There are several widely used clinical scoring systems for RILD based on a combination of symptoms and radiological severity [38]. The morphological changes are described imprecisely and terms such as ‘patchy’ and ‘dense’ are not formally defined. The clinical outcomes are categorised largely on their required therapeutic interventions and there is little emphasis on the functional impact on the patient [6].
Faria et al. have highlighted that two of the most commonly used scoring systems (RTOG/EORTC and NCI-CTC) are very poorly correlated with each other, and that radiological toxicity was rarely associated with symptoms [6]. Tucker et al. [38] have similarly shown that using NCI-CTC, CTCAE and RTOG as scoring systems for RILD give rise to markedly different normal tissue complication probabilities for a particular MLD, suggesting that the existing radiobiological models of RILD are not reliable.
A number of attempts have been made to improve upon these manual classification systems. Computational approaches have the advantages of allowing quantitative, continuous, objective and automated classification of radiological or functional lung damage. One of the most common approaches has been to use Hounsfield Unit (HU) density [31,32,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55]. Palma et al. [52] performed deformable image registration on the phase of the planning scan with the lung volume most similar to the 3-month follow-up CT images of patients treated with SABR. Rigid registration was applied manually and then a modified B-spline Free Form Deformation algorithm was used to warp structures to achieve the required 3D displacement. Voxel HU density histograms were created and mean lung densities were derived. There was very poor correlation between HU density changes in the whole lung and the severity of physician graded radiological pneumonitis; however, local density changes around the PTV correlated strongly with increased radiological pneumonitis (Spearman’s r = 0.75).
A more recent example of the same technique [54] studied 31 patients receiving SABR to investigate the relationship between normal lung CT density changes with dose accuracy and outcome. Each patient was assigned a CTCAE RP grade. HU changes in 5 Gy dose bins from 5–45 Gy were assessed in the peri-tumoural region (ITV+ 3 cm margin). The 0–5 Gy lung volume was used as a baseline correction of the density changes. The average lung density changes in the peritumoural region for each of the 5–45 Gy dose bins were tabulated and compared across different dose algorithms. There was a strong positive relationship between peritumoural lung density changes and RP grade (Spearman’s r = 0.76). Positive correlation was also observed between RP and HU changes in the region covered by V20 for all algorithms (Spearman’s r ≥ 0.738). Additionally, V20, MLD and gEUD (generalised equivalent uniform dose) were significantly correlated with RP grade (p < 0.01).
Bernchou [32,56,57] et al. used a similar approach to investigate longitudinal change of HU density after IMRT for NSCLC. They found that normal tissue showed a significant increase in HU density after RT within even low dose areas. They noted that the evolution of changes differed between the low- and high-dose regions. Lung parenchyma receiving doses <45 Gy underwent a decrease in HU density after 3 months, while areas receiving 50–60 Gy became denser between 3–9 months before then decreasing again. Beyond 12 months, the density changes stabilise across all dose intervals. The bimodal time distribution supports a model of RILD characterised by early (RP) and late (fibrotic) changes.
The advantage of using HU density includes its quantitative, objective and automated characterisation of RILD. There are similarities between RILD graded by traditional, physician-assessed measures and using HU density. However, the technique has important limitations. Firstly, by simply recording HU density, it is not possible to distinguish different causes for HU changes, for example tumour recurrence, infection, intravascular contrast injection or lobar collapse. Parenchymal changes are heterogenous in aetiology and morphology and to reduce them to mean HU discounts information contained in the CT scan. Simply measuring different average HU in an area of lung cannot distinguish between a small change in density at all voxels or a large increase in density across a small number of voxels. Furthermore, HU density can be confounded by a number of scan artefacts, such as the use of contrast, the presence of vessels, differing scan acquisition protocols, respiratory motion, etc. The advantage of our classification system over the ones above is that by developing a novel texture-based analysis, we are able to describe density changes with greater richness than simply using HU density.
Other groups have developed sophisticated techniques for analysing the lung parenchyma after radiotherapy [50,51] by using radiomics. While these techniques can uncover relationships between texture information in the images and the development of radiation pneumonitis and/or the radiotherapy dose, the radiomic features they use to do this do not provide a direct classification of the radiological changes, as seen in the scans by a human observer, as the tissue classes used in this work aim to do. The long-term aim of this work is to better understand the radiological manifestations of RILD, their evolution over time, and the complex relationships they have with dose, clinical outcomes, and other factors. Therefore, we have developed novel tools, such as the automated parenchymal tissue classification system presented in our accompanying paper [25], and our method for registering heavily damaged lungs [24], which will enable us to achieve a better understanding of RILD, rather than employing radiomics or modern learning-based approaches that may produce strong predictive models but do not provide any insights that can help further our understanding of RILD.
One of the potential weaknesses of our classification system is that the individual classes binned several distinct radiological entities. Class 5, for example, which represents opacity on a CT scan, could be the result of tumour, lung collapse or pleural effusion. Pleural effusions, in particular, often take up a significant proportion of the intrathoracic cavity. Being idiosyncratic and temporally unstable, they may have weakened the dosimetric correlations seen with this class. It may be helpful to introduce additional classes to represent for example, tumour, atelectasis and pleural effusion in order to further increase the utility of the classification schema. Nonetheless, we have developed a far richer classification of texture changes than that employed by using HU density alone.
The study presented in this paper is limited due to relatively low patient numbers, and therefore, did not aim to provide definitive clinical conclusions about RILD. Rather, it is a proof of principal study that aimed to demonstrate the potential of our tools to provide valuable insights into RILD, and the relationships between parenchymal tissue damage and the global and local dose and respiratory function. Another limitation of this study is that our dataset is from a trial with an isotoxic dose design, and there may be less heterogeneity between patients than in a non-trial data set, which may have masked correlations between dose and parenchymal texture. The IDEAL trial was published in 2016, and accordingly, most of the patients received 3D conformal RT. In future work, we will fully automate our methods so that they can be applied to large datasets with many patients. We will then employ these to study both large retrospective cohorts of non-trial patients, including those treated with VMAT radiotherapy [58], as well as prospective datasets that will include richer clinical data, such as patients’ reported outcome measures.

5. Conclusions

We have demonstrated that the recent tools we have developed for studying the radiological manifestations of RILD can provide novel insights into the temporal evolution of RILD and its relationship to global and local dose and respiratory outcomes. Our registration method for heavily damaged lungs was able to successfully register the pre-RT scan to the follow-up CT scans in 45 of 46 patients. The parenchymal tissue classes we developed demonstrated statistical correlation to both global and local dose metrics in our study, and have a distinct evolution over time. The dosimetric variables most strongly correlated with change in tissue class volume are lung V20, V30 Gy and mean lung dose. Although less strong, there is a relationship between the tissue class changes and respiratory outcomes, particularly the MRC dyspnoea score, which directly represents a patient’s functional status. We have demonstrated the potential of using our tools to analyse and understand the evolution of the radiological manifestation of RILD in greater detail than previously possible, and we hope this can ultimately be used to inform the refinement of RT treatment in order to reduce the burden of morbidity on lung cancer patients as they begin to live longer with their disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14040946/s1, Figure S1: Pearson Correlation Matrix of Dosimetric metrics against Parenchymal Classes at 6 months, Figure S2: Pearson Correlation Matrix of Dosimetric metrics against Parenchymal Classes at 12 months, Figure S3: Pearson Correlation Matrix of Dosimetric metrics against Parenchymal Classes at 24 months, Figure S4: Pearson Correlation Matrix of Respiratory metrics against Parenchymal Classes at 3 months, Figure S5: Pearson Correlation Matrix of Respiratory metrics against Parenchymal Classes at 6 months, Figure S6: Pearson Correlation Matrix of Respiratory metrics against Parenchymal Classes at 12 months.

Author Contributions

Conceptualisation, A.S. (Adam Szmul), E.C., C.V., D.L., M.A.H., J.D.F., S.G., J.W., C.H. and J.R.M.; Data curation, E.C. and J.J.; Funding acquisition, D.L.; Investigation, A.S., E.C., J.J., A.S. (Alkisti Stavropoulou) and J.R.M.; Methodology, A.S. (Adam Szmul), E.C., C.V., J.J., D.L., C.H. and J.R.M.; Project administration, A.S. (Adam Szmul), E.C. and J.R.M.; Software, A.S. (Adam Szmul). All authors have read and agreed to the published version of the manuscript.

Funding

A.S. was supported by the EPSRC-funded UCL Centre for Doctoral Training in Intelligent, Integrated Imaging in Healthcare (i4health) (EP/S021930/1). J.J. was supported by the Wellcome Trust Clinical Research Career Development Fellowship 209,553/Z/17/Z and the NIHR UCLH Biomedical Research Centre, UK. This research was funded in whole or in part by the Wellcome Trust [209553/Z/17/Z]. For the purpose of open access, J.J. has applied a CC-BY public copyright licence to any author-accepted manuscript version arising from this submission. M.A.H. is supported by funding from the NIHR Biomedical Research Centre at University College London Hospitals NHS Foundation Trust. C.V. is supported by the Royal Academy of Engineering under the Research Fellowship scheme (RF\201718\17140). J.R.M. is supported by a CRUK Centres Network Accelerator Award Grant (A21993) to the ART-NET consortium. The IDEAL CRT trial was funded by Cancer Research UK, grant no. C13530/A10424 and C13530/A17007. The authors would like to thank the Cancer Research UK and UCL Trials Unit.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by a data sharing agreement between the University College London Cancer Trials Centre and the Centre for Medical Image Computing, University College London (07/01/2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

J.J. reports fees from Boehringer Ingelheim, Roche, NHSX, Takeda, and GlaxoSmithKline, unrelated to the submitted work. The authors declare no other conflict of interests.

References

  1. Barta, J.A.; Powell, C.A.; Wisnivesky, J.P. Global Epidemiology of Lung Cancer. Ann. Glob. Health 2019, 85, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Faivre-Finn, C.; Vicente, D.; Kurata, T.; Planchard, D.; Paz-Ares, L.; Vansteenkiste, J.F.; Spigel, D.R.; Garassino, M.C.; Reck, M.; Senan, S.; et al. Four-Year Survival with Durvalumab after Chemoradiotherapy in Stage III NSCLC—An Update from the PACIFIC Trial. J. Thorac. Oncol. 2021, 16, 860–867. [Google Scholar] [CrossRef] [PubMed]
  3. Giuranno, L.; Ient, J.; De Ruysscher, D.; Vooijs, M.A. Radiation-Induced Lung Injury (RILI). Front. Oncol. 2019, 9, 877. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hanania, A.N.; Mainwaring, W.; Ghebre, Y.T.; Hanania, N.A.; Ludwig, M. Radiation-Induced Lung Injury: Assessment and Management. Chest 2019, 156, 150–162. [Google Scholar] [CrossRef]
  5. Morgan, G.W.; Pharm, B.; Breit, S.N. Radiation and the lung: A reevaluation of the mechanisms mediating pulmonary injury. Int. J. Radiat. Oncol. Biol. Phys. 1995, 31, 361–369. [Google Scholar] [CrossRef]
  6. Faria, S.L.; Aslani, M.; Tafazoli, F.S.; Souhami, L.; Freeman, C.R. The challenge of scoring radiation-induced lung toxicity. Clin. Oncol. (R. Coll. Radiol.) 2009, 21, 371–375. [Google Scholar] [CrossRef]
  7. Palma, D.A.; Senan, S.; Tsujino, K.; Barriger, R.B.; Rengan, R.; Moreno, M.; Bradley, J.D.; Kim, T.H.; Ramella, S.; Marks, L.B.; et al. Predicting radiation pneumonitis after chemoradiation therapy for lung cancer: An international individual patient data meta-analysis. Int. J. Radiat. Oncol. Biol. Phys. 2013, 85, 444–450. [Google Scholar] [CrossRef] [Green Version]
  8. Minami-Shimmyo, Y.; Ohe, Y.; Yamamoto, S.; Sumi, M.; Nokihara, H.; Horinouchi, H.; Yamamoto, N.; Sekine, I.; Kubota, K.; Tamura, T. Risk factors for treatment-related death associated with chemotherapy and thoracic radiotherapy for lung cancer. J. Thorac. Oncol. 2012, 7, 177–182. [Google Scholar] [CrossRef] [Green Version]
  9. Kong, F.-M.; Wang, S. Nondosimetric Risk Factors for Radiation-Induced Lung Toxicity. Semin. Radiat. Oncol. 2015, 25, 100–109. [Google Scholar] [CrossRef] [Green Version]
  10. Jain, V.; Berman, A.T. Radiation Pneumonitis: Old Problem, New Tricks. Cancers 2018, 10, 222. [Google Scholar] [CrossRef] [Green Version]
  11. Delaunay, M.; Prévot, G.; Collot, S.; Guilleminault, L.; Didier, A.; Mazières, J. Management of pulmonary toxicity associated with immune checkpoint inhibitors. Eur. Respir. Rev. 2019, 28, 190012. [Google Scholar] [CrossRef] [PubMed]
  12. Twyman-Saint Victor, C.; Rech, A.J.; Maity, A.; Rengan, R.; Pauken, K.E.; Stelekati, E.; Benci, J.L.; Xu, B.; Dada, H.; Odorizzi, P.M.; et al. Radiation and dual checkpoint blockade activate non-redundant immune mechanisms in cancer. Nature 2015, 520, 373–377. [Google Scholar] [CrossRef] [Green Version]
  13. Shaverdian, N.; Lisberg, A.E.; Bornazyan, K.; Veruttipong, D.; Goldman, J.W.; Formenti, S.C.; Garon, E.B.; Lee, P. Previous radiotherapy and the clinical activity and toxicity of pembrolizumab in the treatment of non-small-cell lung cancer: A secondary analysis of the KEYNOTE-001 phase 1 trial. Lancet Oncol. 2017, 18, 895–903. [Google Scholar] [CrossRef]
  14. Marks, L.B.; Bentzen, S.M.; Deasy, J.O.; Kong, F.-M.S.; Bradley, J.D.; Vogelius, I.S.; El Naqa, I.; Hubbs, J.L.; Lebesque, J.V.; Timmerman, R.D.; et al. Radiation dose-volume effects in the lung. Int. J. Radiat. Oncol. Biol. Phys. 2010, 76 Suppl. S3, S70–S76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bradley, J.D.; Hu, C.; Komaki, R.U.; Masters, G.; Blumenschein, G.R.; Schild, S.E.; Bogart, J.A.; Forster, K.M.; Magliocco, A.; Kavadi, V.S.; et al. Long-Term Results of RTOG 0617: A Randomized Phase 3 Comparison of Standard Dose Versus High Dose Conformal Chemoradiation Therapy +/- Cetuximab for Stage III NSCLC. Int. J. Radiat. Oncol. Biol. Phys. 2017, 99, S105. [Google Scholar] [CrossRef]
  16. O’Rourke, N.; Roque, I.F.M.; Farre Bernado, N.; Macbeth, F. Concurrent chemoradiotherapy in non-small cell lung cancer. Cochrane Database Syst Rev 2010, 6, Cd002140. [Google Scholar] [CrossRef]
  17. Libshitz, H.I.; Shuman, L.S. Radiation-Induced Pulmonary Change: CT Findings. J. Comput. Assist. Tomogr. 1984, 8, 15–19. [Google Scholar] [CrossRef]
  18. Ikezoe, J.; Takashima, S.; Morimoto, S.; Kadowaki, K.; Takeuchi, N.; Yamamoto, T.; Nakanishi, K.; Isaza, M.; Arisawa, J.; Ikeda, H.; et al. CT appearance of acute radiation-induced injury in the lung. AJR Am. J. Roentgenol. 1988, 150, 765–770. [Google Scholar] [CrossRef] [Green Version]
  19. Trovo, M.; Linda, A.; El Naqa, I.; Javidan-Nejad, C.; Bradley, J. Early and late lung radiographic injury following stereotactic body radiation therapy (SBRT). Lung Cancer 2010, 69, 77–85. [Google Scholar] [CrossRef]
  20. Koenig, T.R.; Munden, R.F.; Erasmus, J.J.; Sabloff, B.S.; Gladish, G.W.; Komaki, R.; Stevens, C.W. Radiation injury of the lung after three-dimensional conformal radiation therapy. AJR Am. J. Roentgenol. 2002, 178, 1383–1388. [Google Scholar] [CrossRef] [PubMed]
  21. Linda, A.; Trovo, M.; Bradley, J.D. Radiation injury of the lung after stereotactic body radiation therapy (SBRT) for lung cancer: A timeline and pattern of CT changes. Eur. J. Radiol. 2011, 79, 147–154. [Google Scholar] [CrossRef]
  22. Veiga, C.; Landau, D.; Devaraj, A.; Doel, T.; White, J.; Ngai, Y.; Hawkes, D.J.; McClelland, J.R. Novel CT-Based Objective Imaging Biomarkers of Long-Term Radiation-Induced Lung Damage. Int. J. Radiat. Oncol. Biol. Phys. 2018, 102, 1287–1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Veiga, C.; Chandy, E.; Jacob, J.; Yip, N.; Szmul, A.; Landau, D.; McClelland, J.R. Investigation of the evolution of radiation-induced lung damage using serial CT imaging and pulmonary function tests. Radiother. Oncol. 2020, 148, 89–96. [Google Scholar] [CrossRef] [PubMed]
  24. Stavropoulou, A.; Szmul, A.; Chandy, E.; Veiga, C.; Landau, D.; McClelland, J.R. A multichannel feature-based approach for longitudinal lung CT registration in the presence of radiation induced lung damage. Phys. Med. Biol. 2021, 66, 175020. [Google Scholar] [CrossRef] [PubMed]
  25. Szmul, A.; Chandy, E.; Veiga, C.; Jacob, J.; Stavropoulou, A.; Landau, D.; Hiley, C.; McClelland, J. A novel and automated approach to classify Radiation Induced Lung Tissue Damage on CT scans. Cancers 2021. [Google Scholar]
  26. Fenwick, J.D.; Landau, D.B.; Baker, A.T.; Bates, A.T.; Eswar, C.; Garcia-Alonso, A.; Harden, S.V.; Illsley, M.C.; Laurence, V.; Malik, Z.; et al. Long-Term Results from the IDEAL-CRT Phase 1/2 Trial of Isotoxically Dose-Escalated Radiation Therapy and Concurrent Chemotherapy for Stage II/III Non-small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2020, 106, 733–742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Libshitz, H.I.; Southard, M.E. Complications of radiation therapy: The thorax. Semin. Roentgenol. 1974, 9, 41–49. [Google Scholar] [CrossRef]
  28. Pagani, J.J.; Libshitz, H.I. CT Manifestations of Radiation-Induced Change in Chest Tissue. J. Comput. Assist. Tomogr. 1982, 6, 243–248. [Google Scholar] [CrossRef]
  29. Libshitz, H.I. Radiation changes in the lung. Semin. Roentgenol. 1993, 28, 303–320. [Google Scholar] [CrossRef]
  30. Bernchou, U.; Christiansen, R.L.; Asmussen, J.T.; Schytte, T.; Hansen, O.; Brink, C. Extent and computed tomography appearance of early radiation induced lung injury for non-small cell lung cancer. Radiother. Oncol. 2017, 123, 93–98. [Google Scholar] [CrossRef]
  31. Bernchou, U.; Hansen, O.; Schytte, T.; Bertelsen, A.; Hope, A.; Moseley, D.; Brink, C. Prediction of lung density changes after radiotherapy by cone beam computed tomography response markers and pre-treatment factors for non-small cell lung cancer patients. Radiother. Oncol. 2015, 117, 17–22. [Google Scholar] [CrossRef] [PubMed]
  32. Bernchou, U.; Schytte, T.; Bertelsen, A.; Bentzen, S.M.; Hansen, O.; Brink, C. Time evolution of regional CT density changes in normal lung after IMRT for NSCLC. Radiother. Oncol. 2013, 109, 89–94. [Google Scholar] [CrossRef] [PubMed]
  33. Frangi, A.F.; Niessen, W.J.; Vincken, K.L.; Viergever, M.A. Multiscale vessel enhancement filtering. In Medical Image Computing and Computer-Assisted Intervention—MICCAI’98; Wells, W.M., Colchester, A., Delp, S., Eds.; Springer: Berlin/Heidelberg, Germany, 1998; pp. 130–137. [Google Scholar]
  34. Doel, T. 2012. Available online: https://github.com/tomdoel/pulmonarytoolkit (accessed on 12 February 2021).
  35. Yushkevich, P.A.; Piven, J.; Hazlett, H.C.; Smith, R.G.; Ho, S.; Gee, J.C.; Gerig, G. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage 2006, 31, 1116–1128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Modat, M.; Ridgway, G.R.; Taylor, Z.A.; Lehmann, M.; Barnes, J.; Hawkes, D.J.; Fox, N.C.; Ourselin, S. Fast free-form deformation using graphics processing units. Comput. Methods Programs Biomed. 2010, 98, 278–284. [Google Scholar] [CrossRef] [Green Version]
  37. Agnew, M. Spirometry in clinical use: Practical issues. Breathe 2010, 6, 196. [Google Scholar] [CrossRef]
  38. Tucker, S.L.; Jin, H.; Wei, X.; Wang, S.; Martel, M.K.; Komaki, R.; Liu, H.H.; Mohan, R.; Chen, Y.; Cox, J.D.; et al. Impact of toxicity grade and scoring system on the relationship between mean lung dose and risk of radiation pneumonitis in a large cohort of patients with non-small cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 2010, 77, 691–698. [Google Scholar] [CrossRef]
  39. Schroder, C.; Engenhart-Cabillic, R.; Kirschner, S.; Blank, E.; Buchali, A. Changes of lung parenchyma density following high dose radiation therapy for thoracic carcinomas—An automated analysis of follow up CT scans. Radiat. Oncol. 2019, 14, 72. [Google Scholar] [CrossRef]
  40. Alharbi, M.; Janssen, S.; Golpon, H.; Bremer, M.; Henkenberens, C. Temporal and spatial dose distribution of radiation pneumonitis after concurrent radiochemotherapy in stage III non-small cell cancer patients. Radiat. Oncol. 2017, 12, 165. [Google Scholar] [CrossRef] [Green Version]
  41. Defraene, G.; van Elmpt, W.; Crijns, W.; De Ruysscher, D. Regional variability in radiation-induced lung damage can be predicted by baseline CT numbers. Radiother. Oncol. 2017, 122, 300–306. [Google Scholar] [CrossRef]
  42. Sharifi, H.; van Elmpt, W.; Oberije, C.; Nalbantov, G.; Das, M.; Ollers, M.; Lambin, P.; Dingmans, A.M.; De Ruysscher, D. Quantification of CT-assessed radiation-induced lung damage in lung cancer patients treated with or without chemotherapy and cetuximab. Acta Oncol. 2016, 55, 156–162. [Google Scholar] [CrossRef] [Green Version]
  43. Defraene, G.; van Elmpt, W.; Crijns, W.; Slagmolen, P.; De Ruysscher, D. CT characteristics allow identification of patient-specific susceptibility for radiation-induced lung damage. Radiother. Oncol. 2015, 117, 29–35. [Google Scholar] [CrossRef] [PubMed]
  44. Ghobadi, G.; Wiegman, E.M.; Langendijk, J.A.; Widder, J.; Coppes, R.P.; van Luijk, P. A new CT-based method to quantify radiation-induced lung damage in patients. Radiother. Oncol. 2015, 117, 4–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Avanzo, M.; Trovo, M.; Furlan, C.; Barresi, L.; Linda, A.; Stancanello, J.; Andreon, L.; Minatel, E.; Bazzocchi, M.; Trovo, M.G.; et al. Normal tissue complication probability models for severe acute radiological lung injury after radiotherapy for lung cancer. Phys. Med. 2015, 31, 1–8. [Google Scholar] [CrossRef] [PubMed]
  46. Diot, Q.; Marks, L.B.; Bentzen, S.M.; Senan, S.; Kavanagh, B.D.; Lawrence, M.V.; Miften, M.; Palma, D.A. Comparison of radiation-induced normal lung tissue density changes for patients from multiple institutions receiving conventional or hypofractionated treatments. Int. J. Radiat. Oncol. Biol. Phys. 2014, 89, 626–632. [Google Scholar] [CrossRef]
  47. De Ruysscher, D.; Sharifi, H.; Defraene, G.; Kerns, S.L.; Christiaens, M.; De Ruyck, K.; Peeters, S.; Vansteenkiste, J.; Jeraj, R.; Van Den Heuvel, F.; et al. Quantification of radiation-induced lung damage with CT scans: The possible benefit for radiogenomics. Acta Oncol. 2013, 52, 1405–1410. [Google Scholar] [CrossRef] [Green Version]
  48. Lawrence, M.V.; Saynak, M.; Fried, D.V.; Bateman, T.A.; Green, R.L.; Hubbs, J.L.; Jaszczak, R.J.; Wong, T.Z.; Zhou, S.; Das, S.K.; et al. Assessing the impact of radiation-induced changes in soft tissue density thickness on the study of radiation-induced perfusion changes in the lung and heart. Med. Phys. 2012, 39, 7644–7649. [Google Scholar] [CrossRef] [Green Version]
  49. Mattonen, S.A.; Palma, D.A.; Haasbeek, C.J.; Senan, S.; Ward, A.D. Distinguishing radiation fibrosis from tumour recurrence after stereotactic ablative radiotherapy (SABR) for lung cancer: A quantitative analysis of CT density changes. Acta Oncol. 2013, 52, 910–918. [Google Scholar] [CrossRef]
  50. Diot, Q.; Kavanagh, B.; Schefter, T.; Gaspar, L.; Stuhr, K.; Miften, M. Regional normal lung tissue density changes in patients treated with stereotactic body radiation therapy for lung tumors. Int. J. Radiat. Oncol. Biol. Phys. 2012, 84, 1024–1030. [Google Scholar] [CrossRef]
  51. Lee, S.; Stroian, G.; Kopek, N.; AlBahhar, M.; Seuntjens, J.; El Naqa, I. Analytical modelling of regional radiotherapy dose response of lung. Phys. Med. Biol. 2012, 57, 3309–3321. [Google Scholar] [CrossRef]
  52. Palma, D.A.; van Sornsen de Koste, J.R.; Verbakel, W.F.; Senan, S. A new approach to quantifying lung damage after stereotactic body radiation therapy. Acta Oncol. 2011, 50, 509–517. [Google Scholar] [CrossRef] [PubMed]
  53. Ghobadi, G.; Hogeweg, L.E.; Faber, H.; Tukker, W.G.; Schippers, J.M.; Brandenburg, S.; Langendijk, J.A.; Coppes, R.P.; van Luijk, P. Quantifying local radiation-induced lung damage from computed tomography. Int. J. Radiat. Oncol. Biol. Phys. 2010, 76, 548–556. [Google Scholar] [CrossRef]
  54. Al Feghali, K.A.; Wu, Q.C.; Devpura, S.; Liu, C.; Ghanem, A.I.; Wen, N.W.; Ajlouni, M.; Simoff, M.J.; Movsas, B.; Chetty, I.J. Correlation of normal lung density changes with dose after stereotactic body radiotherapy (SBRT) for early stage lung cancer. Clin. Transl. Radiat. Oncol. 2020, 22, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Wuschner, A.E.; Wallat, E.M.; Flakus, M.J.; Shanmuganayagam, D.; Meudt, J.; Christensen, G.E.; Reinhardt, J.M.; Miller, J.R.; Lawless, M.J.; Baschnagel, A.M.; et al. Radiation-induced Hounsfield unit change correlates with dynamic CT perfusion better than 4DCT-based ventilation measures in a novel-swine model. Sci. Rep. 2021, 11, 13156. [Google Scholar] [CrossRef] [PubMed]
  56. Hirose, T.-A.; Arimura, H.; Ninomiya, K.; Yoshitake, T.; Fukunaga, J.-I.; Shioyama, Y. Radiomic prediction of radiation pneumonitis on pretreatment planning computed tomography images prior to lung cancer stereotactic body radiation therapy. Sci. Rep. 2020, 10, 20424. [Google Scholar] [CrossRef]
  57. Cunliffe, A.; Armato, S.G., 3rd; Castillo, R.; Pham, N.; Guerrero, T.; Al-Hallaq, H.A. 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, 1048–1056. [Google Scholar] [CrossRef] [Green Version]
  58. Zhang, T.; Zhou, Z.; Bi, N.; Wang, J.; Wang, L.; Deng, L.; Feng, Q.; Liang, J.; Xiao, Z.; Chen, D.; et al. VMAT for Unresectable Locally Advanced NSCLC Does Not Increase the Risk of Radiation Pneumonitis Compared with IMRT. Int. J. Radiat. Oncol. Biol. Phys. 2019, 105, E543. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Examples of morphological classes used in parenchymal classification system. Left plain CT. Right automated labels. Class 1; Red, Class 2; Green, Class 3; Blue, Class 4; Magenta; Class 5; Yellow.
Figure 1. Examples of morphological classes used in parenchymal classification system. Left plain CT. Right automated labels. Class 1; Red, Class 2; Green, Class 3; Blue, Class 4; Magenta; Class 5; Yellow.
Cancers 14 00946 g001
Figure 2. An example of how the dose was propagated from the planning scan onto the 6- and 12-month time points.
Figure 2. An example of how the dose was propagated from the planning scan onto the 6- and 12-month time points.
Cancers 14 00946 g002
Figure 3. Mean % Volume of Tissue Classes for all patients in the entire lung volume including PTV, at each time point. Red = Class 1, Green = Class 2, Blue = Class 3, Purple = Class 4, Yellow = Class 5.
Figure 3. Mean % Volume of Tissue Classes for all patients in the entire lung volume including PTV, at each time point. Red = Class 1, Green = Class 2, Blue = Class 3, Purple = Class 4, Yellow = Class 5.
Cancers 14 00946 g003
Figure 4. Boxplot showing difference in volumes of each class at follow up and baseline. Solid lines signify statistically significant (p < 0.05) differences between time points. The area within the PTV at each time point is excluded.
Figure 4. Boxplot showing difference in volumes of each class at follow up and baseline. Solid lines signify statistically significant (p < 0.05) differences between time points. The area within the PTV at each time point is excluded.
Cancers 14 00946 g004
Figure 5. Pearson correlation matrix of global dose metrics and change in relative volume between pre-RT and 3-months post RT for each tissue class. VX = Lung VX Gray. MLD = Mean Lung Dose. PTV = Planning Target Volume, GTV = Gross Tumour Volume. Max dose = maximum prescribed dose (Gy). White boxes, p ≥ 0.1.
Figure 5. Pearson correlation matrix of global dose metrics and change in relative volume between pre-RT and 3-months post RT for each tissue class. VX = Lung VX Gray. MLD = Mean Lung Dose. PTV = Planning Target Volume, GTV = Gross Tumour Volume. Max dose = maximum prescribed dose (Gy). White boxes, p ≥ 0.1.
Cancers 14 00946 g005
Figure 6. Pearson correlation matrix of respiratory metrics and change in relative volume between pre-RT and 24-months post RT for each tissue class. FVC = Forced Vital Capacity, FEV = Forced Expiratory Volume, DLCO = diffusing capacity for carbon monoxide, MRC = MRC dyspnoea score. White boxes, p ≥ 0.1.
Figure 6. Pearson correlation matrix of respiratory metrics and change in relative volume between pre-RT and 24-months post RT for each tissue class. FVC = Forced Vital Capacity, FEV = Forced Expiratory Volume, DLCO = diffusing capacity for carbon monoxide, MRC = MRC dyspnoea score. White boxes, p ≥ 0.1.
Cancers 14 00946 g006
Figure 7. Stacked bar graph visualising the relationship between the local dose and the distribution of the tissue classes, and how this changed over the different timepoints. The different colours correspond to the different tissue classes and each column represents a 5 Gy dose bin from 0–5 Gy–60–65 Gy (with a final column for >65 Gy) from each of the different time points.
Figure 7. Stacked bar graph visualising the relationship between the local dose and the distribution of the tissue classes, and how this changed over the different timepoints. The different colours correspond to the different tissue classes and each column represents a 5 Gy dose bin from 0–5 Gy–60–65 Gy (with a final column for >65 Gy) from each of the different time points.
Cancers 14 00946 g007
Table 1. Patient characteristics and accompanying respiratory data.
Table 1. Patient characteristics and accompanying respiratory data.
Patient Characteristics Number (%)Median (Range)
Age 64 (42–83)
GenderMale31 (67)
Female15 (33)
AJCC StageIIB3 (7)
IIIA30 (65)
IIIB13 (28)
Fractionation Schedule6 weeks36 (78)
5 weeks10 (22)
Radiotherapy TechniqueConformal 3D45 (98)
VMAT1 (2)
Prescribed dose (Gy) 66.75 (63–73)
PTV (Planning Tumour Volume) (cm3) 360.00 (139–821)
MLD (Mean Lung Dose) 14.56 (8.75–19.96)
Lung V20 Gy 22.58 (13.86–43.61)
ProgressionAll20 (43)
Loco-regional15 (33)
Table 2. IDEAL-CRT Respiratory data available for patients included in this study.
Table 2. IDEAL-CRT Respiratory data available for patients included in this study.
MetricTime PointAvailable (out of 46)
FVCBaseline46
3 m43
6 m40
12 m40
24 m35
FEV1Baseline46
3 m43
6 m40
12 m40
24 m36
TLCOBaseline46
3 m40
6 m37
12 m38
24 m35
MRC ScoreBaseline43
3 m43
6 m43
12 m41
24 m40
Table 3. Lung Parenchyma classification.
Table 3. Lung Parenchyma classification.
ClassDescription
1Normal, healthy or emphysematous lung, without any high-density abnormality and representing most of the lung parenchymal tissue prior to radiation, as well as areas not affected during the radiotherapy.
2Areas mostly characterised by changes similar to ground-glass opacity.
3Areas with mixed ground-glass opacity and overlaid reticulation.
4Mostly solid lung tissue, either aerated opaque tissue or tissue with a density just below dense opacity.
5Homogeneous, dense lung tissue, which could represent a number of pathological entities, including tumour, pleural effusion or collapse.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chandy, E.; Szmul, A.; Stavropoulou, A.; Jacob, J.; Veiga, C.; Landau, D.; Wilson, J.; Gulliford, S.; Fenwick, J.D.; Hawkins, M.A.; et al. Quantitative Analysis of Radiation-Associated Parenchymal Lung Change. Cancers 2022, 14, 946. https://doi.org/10.3390/cancers14040946

AMA Style

Chandy E, Szmul A, Stavropoulou A, Jacob J, Veiga C, Landau D, Wilson J, Gulliford S, Fenwick JD, Hawkins MA, et al. Quantitative Analysis of Radiation-Associated Parenchymal Lung Change. Cancers. 2022; 14(4):946. https://doi.org/10.3390/cancers14040946

Chicago/Turabian Style

Chandy, Edward, Adam Szmul, Alkisti Stavropoulou, Joseph Jacob, Catarina Veiga, David Landau, James Wilson, Sarah Gulliford, John D. Fenwick, Maria A. Hawkins, and et al. 2022. "Quantitative Analysis of Radiation-Associated Parenchymal Lung Change" Cancers 14, no. 4: 946. https://doi.org/10.3390/cancers14040946

APA Style

Chandy, E., Szmul, A., Stavropoulou, A., Jacob, J., Veiga, C., Landau, D., Wilson, J., Gulliford, S., Fenwick, J. D., Hawkins, M. A., Hiley, C., & McClelland, J. R. (2022). Quantitative Analysis of Radiation-Associated Parenchymal Lung Change. Cancers, 14(4), 946. https://doi.org/10.3390/cancers14040946

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop