Next Article in Journal
Risk Classification for Overall Survival by the Neutrophil–Lymphocyte Ratio and the Number of Metastatic Sites in Patients Treated with Pembrolizumab—A Multicenter Collaborative Study in Japan
Next Article in Special Issue
Radiation-Induced Dyspnea in Lung Cancer Patients Treated with Stereotactic Body Radiation Therapy
Previous Article in Journal
Detection of Circulating Tumor Cells and Microbial DNA Fragments in Stage III Colorectal Cancer Patients under Three versus Six Months of Adjuvant Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis

1
Institute of Biostructures and Bioimaging, National Research Council, 80145 Napoli, Italy
2
Department of Advanced Biomedical Sciences, University of Naples “Federico II”, 80131 Napoli, Italy
3
Department of Radiation Oncology, MD Anderson Cancer Center, Houston, TX 77030, USA
4
Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA
5
Department of Radiation Physics, MD Anderson Cancer Center, Houston, TX 77030, USA
*
Authors to whom correspondence should be addressed.
Cancers 2021, 13(14), 3553; https://doi.org/10.3390/cancers13143553
Submission received: 31 May 2021 / Revised: 7 July 2021 / Accepted: 14 July 2021 / Published: 15 July 2021

Abstract

:

Simple Summary

The pathophysiology of radiation pneumonitis (RP) after thoracic cancer radiation treatments is still not completely understood although the identification of underlying RP mechanisms may improve the therapeutic window of thoracic cancer patients. The aim of our retrospective study was to explore the dose–response patterns associated with RP by a multi-center voxel-based analysis. In a heterogeneously treated population of 382 thoracic cancer patients, we confirmed the previously described heart–lung interaction in the development of RP. The empowerment of VBA with a novel description of dose map spatial properties based on probabilistic independent component analysis (PICA) and connectograms provided valuable additional and independent information on the radiobiology of RP.

Abstract

This study investigates the dose–response patterns associated with radiation pneumonitis (RP) in patients treated for thoracic malignancies with different radiation modalities. To this end, voxel-based analysis (VBA) empowered by a novel strategy for the characterization of spatial properties of dose maps was applied. Data from 382 lung cancer and mediastinal lymphoma patients from three institutions treated with different radiation therapy (RT) techniques were analyzed. Each planning CT and biologically effective dose map (α/β = 3 Gy) was spatially normalized on a common anatomical reference. The VBA of local dose differences between patients with and without RP was performed and the clusters of voxels with dose differences that significantly correlated with RP at a p-level of 0.05 were generated accordingly. The robustness of VBA inference was evaluated by a novel characterization for spatial properties of dose maps based on probabilistic independent component analysis (PICA) and connectograms. This lays robust foundations to the obtained findings that the lower parts of the lungs and the heart play a prominent role in the development of RP. Connectograms showed that the dataset can support a radiobiological differentiation between the main heart and lung substructures.

1. Introduction

Radiation therapy (RT) represents a fundamental treatment strategy in thoracic oncology. Radiation-induced morbidities—including radiation pneumonitis (RP)—may, however, limit the effectiveness of RT [1,2]. The current clinical practice in treatment planning relies on safe parameters such as mean lung dose (MLD) and lung volume receiving at least a certain level of dose (e.g., V20 Gy) to estimate lung toxicity risk by assuming that the entire lung tissue shares the same radiobiological patterns [3]. Moreover, several other variables have been shown to be predictive, including the equivalent uniform dose parameter.
Regional difference in lung radiosensitivity is an old story [4,5,6]. Indeed, while a commonly accepted pathophysiological picture considers the lungs as parallel organs, the higher incidence of lung injury in the lower lobes, together with the differences in regional functionality revealed by nuclear imaging studies, indicates highly heterogeneous structural composition, functional capacity, and sensitivity to radiation. A greater risk of RP after irradiation of caudally located lung tumors has been acknowledged in [7]; however, the underlying local mechanisms still elude a thorough understanding.
The momentum of technological changes, including particle beam therapy [8] or hypofractionated RT [9], of the last decade in radiation therapy poses new challenges in outcome modeling and at the same time emphasizes the limit of the traditional dose-volume histogram (DVH)-based toxicity analysis and the normal tissue complication probability (NTCP) modeling philosophy [10,11].
In this study, we investigate the thoracic dose–response patterns for RP in patients treated for thoracic malignancies taking advantage of the pooled data analysis of different thoracic tumors and different RT modalities. A novel approach based on spatial dose distribution analysis—namely, voxel-based analysis (VBA)—is performed to evaluate the significance of dose differences between groups of patients at a voxel level [12,13]. In addition, we apply an innovative characterization of the spatial properties of the dose distribution, which allows to assess the robustness of regional VBA inferences based on probabilistic independent component analysis and dosimetric connectograms [14].

2. Materials and Methods

2.1. Patients

We retrospectively analyzed 382 patients from 4 different groups of patients receiving thoracic RT with different treatment modalities at different institutions. A group of patients was treated for Hodgkin lymphoma (HL) at the University “Federico II” of Napoli (Comitato Etico per le Attività Biomediche, IRB 222-10) with 3D conformal RT (3D-CRT) [15,16]. A second group was treated for locally advanced non-small cell lung cancer (NSCLC) at the University of Texas MD Anderson Cancer Center of Houston (IRB 2008-0133) with image-guided intensity-modulated RT (IMRT), and a third one with passive-scattering proton therapy (PSPT) [17,18]. Lastly, a group of patients was treated for NSCLC at the Memorial Sloan Kettering Cancer Center of New York (IRB #16-142) with stereotactic body RT (SBRT) [19]. All patients were scored for RP graded according to the National Cancer Institute’s Common Terminology Criteria for Adverse Events (CTCAE) version 3.

2.1.1. 3D Conformal Radiation Therapy

Ninety-eight patients received supradiaphragmatic involved-site RT with a median total dose of 30.6 Gy (range: [20.8, 45.0] Gy) in daily fractions of 1.5 to 1.8 Gy. Treatment plans were designed in XiO (Elekta AB, Stockholm, Sweden) using the multigrid superposition dose calculation algorithm. Dose maps were obtained with a dose grid size of 3.0 mm × 3.0 mm × 5.0 mm.

2.1.2. Intensity-Modulated Radiation Therapy

A group of 114 patients was treated for locally advanced non-small cell lung cancer (NSCLC) at the University of Texas MD Anderson Cancer Center of Houston according to an institutional review board-approved protocol (2008-0133) with image-guided intensity-modulated RT (IMRT). Patients received a prescribed dose of 66 or 74 Gy, given in 33 or 37 conventional daily fractions (2 Gy). Plans were designed in Pinnacle3 (Philips Medical Systems, Andover, MA, USA) using the convolution/superposition dose calculation algorithm. Dose maps were obtained with a dose grid size of 2.0 mm × 2.0 mm × 2.5 mm.

2.1.3. Stereotactic Body Radiation Therapy

One hundred and six patients were treated at the Memorial Sloan Kettering Cancer Center of New York with stereotactic body RT (SBRT) for NSCLC. Patients received a median total dose of 50 Gy (range: [40, 54] Gy) in median 4 fractions (range: [3, 5]). Plans were designed in Eclipse v.13 (Varian Medical Systems, Palo Alto, CA, USA) using the AAA dose calculation engine. Dose maps were obtained with a dose grid size of 1.0 mm × 1.0 mm × 2.0 mm.

2.1.4. Passive-Scattering Proton Therapy

A group of 64 patients was treated for locally advanced NSCLC at the University of Texas MD Anderson Cancer Center of Houston according to an institutional review board-approved protocol (2008-0133) with image-guided passive-scattering proton therapy (PSPT). All patients received a prescribed dose of 66 or 74 Gy, given in 33 or 37 conventional daily fractions. Relative biological effectiveness (RBE) was set to 1.1. Plans were obtained by using the proton convolution superposition algorithm implemented in the Eclipse proton therapy planning system (Varian Medical Systems, Inc., Palo Alto, CA, USA). Dose maps were obtained with a dose grid size of 2.0 mm × 2.0 mm × 2.5 mm.
More details on the protocols and treatment characteristics for each cohort are published elsewhere [3,13,18,19].
Physical dose maps were corrected for different fractionation schemes by computing the biologically effective dose (BED) according to an α/β ratio of 3 Gy.

2.2. Spatial Normalization

To allow for voxel-wise comparison of the BED distributions given to different patients, the 382 BED maps were spatially normalized on the common anatomical reference provided by the extended cardiac torso (XCAT) synthetic CT phantom [20]. Briefly, for each patient, the planning CT was registered onto the synthetic CT provided with XCAT after tumor masking [21] via the B-spline elastic image registration algorithm implemented in Elastix [22]. Then, the obtained deformation field was applied to the BED map [23].

2.3. Characterization of Spatial Properties of Dosimetric Data

The homogeneity of the VBA statistical power in the analyzed anatomical region (i.e., the thorax) strongly depends on the uniformity of voxel-wise dose moments. In principle, homogeneous mean (μ) and standard deviation (σ) maps rule out the hypothesis that relevant radiosensitive regions were dampened in the VBA results due to the non-uniform variability of the dose maps. Consequently, μ and σ maps of normalized BED distributions were computed voxel by voxel over the patients and their uniformity was quantitatively evaluated by the Michelson contrast, which for a function I can be computed as:
C M = I max I min I max + I min
where I max and I min are the highest and lowest values of the function. For a positive-valued function I (such as μ and σ maps), 0 C M 1 . For a given fraction 0 < f 1 of the volume V of the support (i.e., the analyzed organ or anatomical apparatus) of I , C M I f can be defined as the minimum C M assumed by the restrictions of I over the subsets of V   S i   |   S i = f V . Since C M I f is a monotonically increasing function, a summary description of the uniformity of I is provided by its area under curve (AUC). The AUC value of C M I f is 0 for the constant maps and tends to 1 for highly inhomogeneous maps.
Moreover, the spatial resolution of the VBA results (maps of general linear model—GLM—coefficients and their significance) may be affected by the correlation between the doses delivered to different anatomical substructures (i.e., pericardium, heart chambers and walls, and lung segments and lobes). Indeed, small-scale radiobiological patterns—anyway larger than the nominal resolution of the dose grid—could be captured by a VBA provided that the correlation length in the different dose distributions is accordingly short. The spatial independence of the doses delivered to each substructure was investigated by means of probabilistic independent component analysis (PICA) [24] and connectogram analysis [25]. The PICA blindly infers the model order of the analyzed dataset, which corresponds to the number of statistically significant independent spatial maps that generate the input. When applied to a dose distributions dataset (in the form of n [patient] × M [voxel]), a PICA model order comparable to the cohort size n points out a significantly mutual independence between patients’ dose distributions and allows for a valuable VBA inference. The connectograms, instead, highlight the most relevant associations between each pair of substructures according to the pairwise significance of the Spearman correlation between the mean doses related to the substructures: the VBA will have more chance to discriminate the GLM properties arising from each substructure in correspondence of a weaker correlation. In the connectogram rings, outwards from the center, the average and standard deviation across patients of the mean doses within each substructure and the average across patients of the dose standard deviations within each substructure were presented.

2.4. Voxel-Based Statistical Analysis

The performed VBA involved the following statistical analysis at a voxel level [26]: a non-parametric permutation test based on the threshold free cluster enhancement of a maximum-T statistics for GLM was used to assess the regional dose differences between patients who developed any-grade RP and who did not. The significance map obtained by the VBA was, thus, implicitly corrected for multiple comparisons. In addition, non-dosimetric variables, which correlated with the considered outcome at a backward stepwise multivariable (MVA) logistic regression model [27], were included in the GLM as nuisance variables. Clusters of voxels with significance level p < 0.05 (S0.05) were generated from the VBA significance map and the mean doses in these clusters were computed.

3. Results

The dose μ and σ maps computed voxel by voxel over the patients for each separate group and for the entire cohort are displayed in Figure 1. The analysis of the uniformity of μ and σ (Figure 2) shows that the entire cohort overall exhibits consistently lower C M μ values and a satisfactorily low AUC value of C M σ compared with each separate group of patients.
The PICA detected increasing model orders (Figure 3) for PSPT (23, corresponding to 36% of the number of patients), IMRT (31, 27%), SBRT (35, 33%), and 3D-CRT (36, 37%) groups. The model order for the whole cohort was 66, which corresponds to 17% of the patient number. While we can observe that the PICA components are spread throughout the heart and lungs, they exhibit apparent patterns specific to the irradiation modalities, most clearly highlighted by the patches associated with a unique combination of supra-threshold PICA components (Figure 4). In particular, the 3D-CRT modality for Hodgkin lymphoma patients is visually characterized in the axial view by a strong antero-posterior field direction, with a massive coverage of the mediastinum, which appears indeed to be covered by a few patches layered in the cranio-caudal direction. Notably, the number of patches sensibly increases moving from the single groups to the entire cohort, and the resolution inhomogeneities, specific to the irradiation modalities, are largely washed out throughout the field-of-view in the merged dataset.
On the other hand, the connectograms identified different dosimetric connectivity patterns depending on the RT techniques and the diseases, providing a valuable tool to visually assess the spatial properties of dose distribution. The connectograms for different anatomical districts of the thorax (Figure 5) show that the weakest (possibly negative) spatial correlations are found in PSPT and IMRT groups, with slightly higher correlations in SBRT dose distributions. Conversely, the 3D-CRT group exhibits—as expected—positively correlated dose values throughout the heart and lungs. A similar trend is found in the connectograms dedicated to the smallest cardiac substructures (Figure 6), while the smaller irradiation fields of SBRT plans enable a valuable disentanglement of most small lung segments (Figure 7 and Figure 8). When the groups of dose maps are combined in the analysis of the whole cohort, the connectograms reveal an overall improvement in dose resolution over the considered anatomical substructures (Figure 9).
The incidence of any-grade RP on the entire dataset was 28%. Of the 382 patients, 37 developed grade 1 RP and 50 developed clinically symptomatic RP (i.e., grade ≥ 2). On the whole dataset, age (median: 64 y; range: [13, 93] y) was the only clinical variable significantly correlated with RP (p = 0.004) that survived at the MVA analysis and was consequently included in the VBA. VBA highlighted (Figure 10) two largely overlapping S0.05 clusters (total volume: 551 cc) significantly associated with RP in the lungs (lungs-S0.05 = 346 cc) and in the heart (heart-S0.05 = 205 cc). The most represented substructures in the S0.05 clusters were the right lower and middle lung lobes, the left and right atrium with their walls, and the pericardium. The mean BED in heart-S0.05 for RP patients was 33.3 Gy, and for patients without RP it was 25.2 Gy. The mean BED in lungs-S0.05 for RP patients was 38.5 Gy, and for patients without RP it was 25.2 Gy.

4. Discussion

Voxel-based analyses were originally introduced in functional imaging to discriminate the local morphological or metabolic patterns between groups of subjects.
In the last decade, techniques developed for imaging VBA have been applied in the field of radiation oncology, and in particular, but not exclusively, to the analysis of radiation-induced toxicity [28,29,30,31,32]. In particular, local dose–response patterns have been evaluated via voxel-wise statistical analysis of dose maps, provided that they were spatially normalized on a common anatomical reference.
VBA applied to both functional imaging and radiobiology typically results in a significance map of radiological or dose differences between groups of patients. In principle, VBA techniques can be formally applied, in a similar way, to both functional imaging and dosimetric datasets; however, there is a substantial difference between the two fields that warrants caution during result interpretation. Indeed, in functional imaging, VBA can identify the metabolic patterns that depend only on the investigated pathophysiology, at least as long as the chosen image contrast, which is usually given by the contrast medium in nuclear medicine and by the pulse sequence in MRI, allows for an effective representation of the underlying functional mechanisms [33,34]. In radiation oncology, on the contrary, VBA can identify those dose patterns that not only depend on the regional radiosensitivity, but are also within the reach of the intrinsic heterogeneity of the dose distributions included in the analyzed dataset.
When applying the VBA to the outcome analysis of radiation treatments, it is therefore essential to account for the intrinsic features of the dose map dataset that could, in principle, limit the very validity of the VBA results. In this study, we tried to preliminarily answer the question, how much radiobiological detail can we identify from a cohort of patients treated with a given radiation therapy technique?
Essentially, we identified two main issues: the power homogeneity and the resolution of VBA, respectively related to the homogeneity of mean and standard deviation of the dose maps and to the correlation between doses to different anatomical structures [14,35]. They seem particularly relevant in the context of thoracic radiation oncology, due to the variety of pathologies and the complexity of pathophysiological pathways proposed in the literature to account for the adverse events in the cardiopulmonary system [15,36,37].
In this context, we exploited for the first time, on a heterogeneously treated population of thoracic cancer patients, a collection of tools devised to assess the spatial properties of dose distributions in order to highlight the power and resolution limits related to the VBAs that can be performed on each dataset [14]. As it has long been expected, it turned out from a rigorously defined analysis that merging different types of dose distributions leads to an improved homogeneity of VBA statistical power as well as to dampened spatial autocorrelation functions. In addition to the straightforward boost of statistical power granted by the increasing number of analyzed patients, which permitted us to obtain by the present study an unparalleled accuracy of the dose–response findings associated with RP of any grade. In particular, it is noteworthy that, despite the valuable independence of the dose values observed in the different substructures in the whole cohort (Figure 9), previously claimed interactions between the heart and the lungs in the RP development [38,39,40,41,42,43] are hereby confirmed. Indeed, we found, in the whole cohort, extended lung and heart regions in which the dose appears significantly correlated with the development of any-grade RP. In these regions an increased mean dose was found in association with RP status: the increase, in patients with RP, was measured to be approximately 13 Gy in lungs-S0.05 and 8 Gy in heart-S0.05.
One potential limitation of this work, which may affect both the strictly speaking VBA and our auxiliary analyses of the intrinsic features of the dose map dataset, is the dosimetric uncertainties associated with the use of planning dose maps instead of delivered accumulated dose maps. However, such uncertainties are not group related; consequently no bias is expected in the analysis described here.

5. Conclusions

The Michelson contrast of voxel-wise dose moments, as well as the dosimetric PICA and connectograms, represents a valuable toolbox to provide essential insights into the homogeneity and resolution limits inherent to a given dataset of dose maps in the VBA context. Heterogeneous dose patterns, steeper gradients, and smaller hot spots all contribute to a more suitable dataset for dosomic inference. We suggest that a spatial characterization of the dose datasets should constitute an ancillary analysis for every voxel-based analysis in the field of radiation oncology in order to better elucidate the pathophysiological mechanisms underlying radiation-induced morbidity development such as RP.

Author Contributions

Conceptualization, G.P.; methodology, G.P., S.M., and L.C.; software, G.P. and S.M.; formal analysis, G.P.; data curation, R.P., J.O.D., Z.L. and L.C.; writing—original draft preparation, G.P., S.M. and L.C.; writing—review and editing, R.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of the University of Texas MD Anderson Cancer Center of Houston (IRB 2008-0133), of the Memorial Sloan Kettering Cancer Center of New York (IRB 16-142), and of the University “Federico II” of Napoli (Comitato Etico per le Attività Biomediche, IRB 222-10).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are not publicly available due to restrictions in the Material Transfer Agreements.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jain, V.; Berman, A.T. Radiation Pneumonitis: Old Problem, New Tricks. Cancers 2018, 10, 222. [Google Scholar] [CrossRef] [Green Version]
  2. Cella, L.; Liuzzi, R.; D’Avino, V.; Conson, M.; Di Biase, A.; Picardi, M.; Pugliese, N.; Solla, R.; Salvatore, M.; Pacelli, R. Pulmonary damage in Hodgkin’s lymphoma patients treated with sequential chemo-radiotherapy: Predictors of radiation-induced lung injury. Acta Oncol. 2014, 53, 613–619. [Google Scholar] [CrossRef] [Green Version]
  3. Palma, G.; Monti, S.; Buonanno, A.; Pacelli, R.; Cella, L. PACE: A Probabilistic Atlas for Normal Tissue Complication Estimation in Radiation Oncology. Front. Oncol. 2019, 9, 130. [Google Scholar] [CrossRef] [Green Version]
  4. Wiegman, E.M.; Meertens, H.; Konings, A.W.; Kampinga, H.H.; Coppes, R.P. Loco-regional differences in pulmonary function and density after partial rat lung irradiation. Radiother. Oncol. 2003, 69, 11–19. [Google Scholar] [CrossRef]
  5. 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]
  6. Madani, I.; De Ruyck, K.; Goeminne, H.; De Neve, W.; Thierens, H.; Van Meerbeeck, J. Predicting risk of radiation-induced lung injury. J. Thorac. Oncol. 2007, 2, 864–874. [Google Scholar] [CrossRef] [Green Version]
  7. Seppenwoolde, Y.; De Jaeger, K.; Boersma, L.J.; Belderbos, J.S.; Lebesque, J.V. Regional differences in lung radiosensitivity after radiotherapy for non-small-cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 2004, 60, 748–758. [Google Scholar] [CrossRef] [PubMed]
  8. Cella, L.; Lomax, A.; Miralbell, R. New techniques in hadrontherapy: Intensity modulated proton beams. Phys. Med. 2001, 17 (Suppl. 1), 100–102. [Google Scholar]
  9. Mirri, M.A.; Arcangeli, G.; Benassi, M.; d’Angelo, A.; Pinzi, V.; Caterino, M.; Rinaldi, M.; Ceribelli, A.; Strigari, L. Hypofractionated Conformal Radiotherapy (HCRT) for primary and metastatic lung cancers with small dimension: Efficacy and toxicity. Strahlenther. Onkol. 2009, 185, 27–33. [Google Scholar] [CrossRef] [PubMed]
  10. Tommasino, F.; Durante, M.; D’Avino, V.; Liuzzi, R.; Conson, M.; Farace, P.; Palma, G.; Schwarz, M.; Cella, L.; Pacelli, R. Model-based approach for quantitative estimates of skin, heart, and lung toxicity risk for left-side photon and proton irradiation after breast-conserving surgery. Acta Oncol. 2017, 56, 730–736. [Google Scholar] [CrossRef] [Green Version]
  11. Palma, G.; Monti, S.; Conson, M.; Pacelli, R.; Cella, L. Normal tissue complication probability (NTCP) models for modern radiation therapy. Semin. Oncol. 2019, 46, 210–218. [Google Scholar] [CrossRef]
  12. Monti, S.; Palma, G.; D’Avino, V.; Gerardi, M.; Marvaso, G.; Ciardo, D.; Pacelli, R.; Jereczek-Fossa, B.A.; Alterio, D.; Cella, L. Voxel-based analysis unveils regional dose differences associated with radiation-induced morbidity in head and neck cancer patients. Sci. Rep. 2017, 7, 7220. [Google Scholar] [CrossRef] [PubMed]
  13. Palma, G.; Monti, S.; D’Avino, V.; Conson, M.; Liuzzi, R.; Pressello, M.C.; Donato, V.; Deasy, J.O.; Quarantelli, M.; Pacelli, R.; et al. A Voxel-Based Approach to Explore Local Dose Differences Associated With Radiation-Induced Lung Damage. Int. J. Radiat. Oncol. Biol. Phys. 2016, 96, 127–133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cella, L.; Monti, S.; Xu, T.; Liuzzi, R.; Stanzione, A.; Durante, M.; Mohan, R.; Liao, Z.; Palma, G. Probing thoracic dose patterns associated to pericardial effusion and mortality in patients treated with photons and protons for locally advanced non-small-cell lung cancer. Radiother. Oncol. 2021. [Google Scholar] [CrossRef]
  15. Cella, L.; Oh, J.H.; Deasy, J.O.; Palma, G.; Liuzzi, R.; D’Avino, V.; Conson, M.; Picardi, M.; Salvatore, M.; Pacelli, R. Predicting radiation-induced valvular heart damage. Acta Oncol. 2015, 54, 1796–1804. [Google Scholar] [CrossRef]
  16. Cella, L.; Palma, G.; Deasy, J.O.; Oh, J.H.; Liuzzi, R.; D’Avino, V.; Conson, M.; Pugliese, N.; Picardi, M.; Salvatore, M.; et al. Complication probability models for radiation-induced heart valvular dysfunction: Do heart-lung interactions play a role? PLoS ONE 2014, 9, e111753. [Google Scholar] [CrossRef]
  17. Liao, Z.; Lee, J.J.; Komaki, R.; Gomez, D.R.; O’Reilly, M.S.; Fossella, F.V.; Blumenschein, G.R., Jr.; Heymach, J.V.; Vaporciyan, A.A.; Swisher, S.G.; et al. Bayesian Adaptive Randomization Trial of Passive Scattering Proton Therapy and Intensity-Modulated Photon Radiotherapy for Locally Advanced Non-Small-Cell Lung Cancer. J. Clin. Oncol. 2018, 36, 1813–1822. [Google Scholar] [CrossRef]
  18. Palma, G.; Monti, S.; Conson, M.; Xu, T.; Hahn, S.; Durante, M.; Mohan, R.; Liao, Z.; Cella, L. NTCP Models for Severe Radiation Induced Dermatitis After IMRT or Proton Therapy for Thoracic Cancer Patients. Front. Oncol. 2020, 10, 344. [Google Scholar] [CrossRef] [PubMed]
  19. Palma, G.; Monti, S.; Thor, M.; Rimner, A.; Deasy, J.O.; Cella, L. Spatial signature of dose patterns associated with acute radiation-induced lung damage in lung cancer patients treated with stereotactic body radiation therapy. Phys. Med. Biol. 2019, 64, 155006. [Google Scholar] [CrossRef] [PubMed]
  20. Segars, W.P.; Sturgeon, G.; Mendonca, S.; Grimes, J.; Tsui, B.M. 4D XCAT phantom for multimodality imaging research. Med. Phys. 2010, 37, 4902–4915. [Google Scholar] [CrossRef]
  21. Conson, M.; Cella, L.; Pacelli, R.; Comerci, M.; Liuzzi, R.; Salvatore, M.; Quarantelli, M. Automated delineation of brain structures in patients undergoing radiotherapy for primary brain tumors: From atlas to dose-volume histograms. Radiother. Oncol. 2014, 112, 326–331. [Google Scholar] [CrossRef] [Green Version]
  22. Klein, S.; Staring, M.; Murphy, K.; Viergever, M.A.; Pluim, J.P. elastix: A toolbox for intensity-based medical image registration. IEEE Trans. Med. Imaging 2010, 29, 196–205. [Google Scholar] [CrossRef]
  23. Palma, G.; Monti, S.; Cella, L. Voxel-based analysis in radiation oncology: A methodological cookbook. Phys. Med. 2020, 69, 192–204. [Google Scholar] [CrossRef] [PubMed]
  24. Beckmann, C.F.; Smith, S.M. Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans. Med. Imaging 2004, 23, 137–152. [Google Scholar] [CrossRef]
  25. Irimia, A.; Chambers, M.C.; Torgerson, C.M.; Van Horn, J.D. Circular representation of human cortical networks for subject and population-level connectomic visualization. Neuroimage 2012, 60, 1340–1351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chen, C.; Witte, M.; Heemsbergen, W.; van Herk, M. Multiple comparisons permutation test for image based data mining in radiotherapy. Radiat. Oncol. 2013, 8, 293. [Google Scholar] [CrossRef] [PubMed]
  27. Palma, G.; Taffelli, A.; Fellin, F.; D’Avino, V.; Scartoni, D.; Tommasino, F.; Scifoni, E.; Durante, M.; Amichetti, M.; Schwarz, M.; et al. Modelling the risk of radiation induced alopecia in brain tumor patients treated with scanned proton beams. Radiother. Oncol. 2020, 144, 127–134. [Google Scholar] [CrossRef]
  28. Mylona, E.; Acosta, O.; Lizee, T.; Lafond, C.; Crehange, G.; Magne, N.; Chiavassa, S.; Supiot, S.; Ospina Arango, J.D.; Campillo-Gimenez, B.; et al. Voxel-Based Analysis for Identification of Urethrovesical Subregions Predicting Urinary Toxicity After Prostate Cancer Radiation Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2019, 104, 343–354. [Google Scholar] [CrossRef] [PubMed]
  29. Drean, G.; Acosta, O.; Lafond, C.; Simon, A.; de Crevoisier, R.; Haigron, P. Interindividual registration and dose mapping for voxelwise population analysis of rectal toxicity in prostate cancer radiotherapy. Med. Phys. 2016, 43, 2721–2730. [Google Scholar] [CrossRef] [PubMed]
  30. Abravan, A.; Faivre-Finn, C.; Kennedy, J.; McWilliam, A.; van Herk, M. Radiotherapy-Related Lymphopenia Affects Overall Survival in Patients With Lung Cancer. J. Thorac. Oncol. 2020, 15, 1624–1635. [Google Scholar] [CrossRef]
  31. Beasley, W.; Thor, M.; McWilliam, A.; Green, A.; Mackay, R.; Slevin, N.; Olsson, C.; Pettersson, N.; Finizia, C.; Estilo, C.; et al. Image-based Data Mining to Probe Dosimetric Correlates of Radiation-induced Trismus. Int. J. Radiat. Oncol. Biol. Phys. 2018, 102, 1330–1338. [Google Scholar] [CrossRef]
  32. Yahya, N.; Ebert, M.A.; House, M.J.; Kennedy, A.; Matthews, J.; Joseph, D.J.; Denham, J.W. Modeling Urinary Dysfunction After External Beam Radiation Therapy of the Prostate Using Bladder Dose-Surface Maps: Evidence of Spatially Variable Response of the Bladder Surface. Int. J. Radiat. Oncol. Biol. Phys. 2017, 97, 420–426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Coda, A.R.; Anzilotti, S.; Boscia, F.; Greco, A.; Panico, M.; Gargiulo, S.; Gramanzini, M.; Zannetti, A.; Albanese, S.; Pignataro, G.; et al. In vivo imaging of CNS microglial activation/macrophage infiltration with combined [(18)F] DPA-714-PET and SPIO-MRI in a mouse model of relapsing remitting experimental autoimmune encephalomyelitis. Eur. J. Nucl. Med. Mol. Imaging 2021, 48, 40–52. [Google Scholar] [CrossRef] [PubMed]
  34. Pontillo, G.; Petracca, M.; Monti, S.; Quarantelli, M.; Criscuolo, C.; Lanzillo, R.; Tedeschi, E.; Elefante, A.; Brescia Morra, V.; Brunetti, A.; et al. Unraveling Deep Gray Matter Atrophy and Iron and Myelin Changes in Multiple Sclerosis. AJNR Am. J. Neuroradiol. 2021. [Google Scholar] [CrossRef]
  35. Shortall, J.; Palma, G.; Mistry, H.; Vasquez Osorio, E.; McWilliam, A.; Choudhury, A.; Aznar, M.; van Herk, M.; Green, A. Flogging a Dead Salmon? Reduced Dose Posterior to Prostate Correlates With Increased PSA Progression in Voxel-Based Analysis of 3 Randomized Phase 3 Trials. Int. J. Radiat. Oncol. Biol. Phys. 2021, 110, 696–699. [Google Scholar] [CrossRef] [PubMed]
  36. Voshart, D.C.; Wiedemann, J.; van Luijk, P.; Barazzuol, L. Regional Responses in Radiation-Induced Normal Tissue Damage. Cancers 2021, 13, 367. [Google Scholar] [CrossRef]
  37. Cella, L.; Gagliardi, G.; Hedman, M.; Palma, G. Injuries From Asymptomatic COVID-19 Disease: New Hidden Toxicity Risk Factors in Thoracic Radiation Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2020, 108, 394–396. [Google Scholar] [CrossRef] [PubMed]
  38. Ghobadi, G.; van der Veen, S.; Bartelds, B.; de Boer, R.A.; Dickinson, M.G.; de Jong, J.R.; Faber, H.; Niemantsverdriet, M.; Brandenburg, S.; Berger, R.M.; et al. Physiological interaction of heart and lung in thoracic irradiation. Int. J. Radiat. Oncol. Biol. Phys. 2012, 84, e639–e646. [Google Scholar] [CrossRef] [Green Version]
  39. Monti, S.; Pacelli, R.; Cella, L.; Palma, G. Inter-patient image registration algorithms to disentangle regional dose bioeffects. Sci. Rep. 2018, 8, 4915. [Google Scholar] [CrossRef]
  40. Tomita, N.; Okuda, K.; Ogawa, Y.; Iida, M.; Eguchi, Y.; Kitagawa, Y.; Uchiyama, K.; Takaoka, T.; Nakanishi, R.; Shibamoto, Y. Relationship between radiation doses to heart substructures and radiation pneumonitis in patients with thymic epithelial tumors. Sci. Rep. 2020, 10, 11191. [Google Scholar] [CrossRef]
  41. Palma, G.; Monti, S.; Xu, T.; Scifoni, E.; Yang, P.; Hahn, S.M.; Durante, M.; Mohan, R.; Liao, Z.; Cella, L. Spatial Dose Patterns Associated With Radiation Pneumonitis in a Randomized Trial Comparing Intensity-Modulated Photon Therapy With Passive Scattering Proton Therapy for Locally Advanced Non-Small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2019, 104, 1124–1132. [Google Scholar] [CrossRef] [PubMed]
  42. Huang, E.X.; Hope, A.J.; Lindsay, P.E.; Trovo, M.; El Naqa, I.; Deasy, J.O.; Bradley, J.D. Heart irradiation as a risk factor for radiation pneumonitis. Acta Oncol. 2011, 50, 51–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Cella, L.; D’Avino, V.; Palma, G.; Conson, M.; Liuzzi, R.; Picardi, M.; Pressello, M.C.; Boboc, G.I.; Battistini, R.; Donato, V.; et al. Modeling the risk of radiation-induced lung fibrosis: Irradiated heart tissue is as important as irradiated lung. Radiother. Oncol. 2015, 117, 36–43. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Axial (columns 1 and 3) and coronal (columns 2 and 4) CT views of the reference anatomy fused with mean (columns 1 and 2) and standard deviations (columns 3 and 4) of BED maps computed voxel by voxel over the patients for the Hodgkin lymphoma group (row a), non-small cell lung cancer (NSCLC) patients treated with intensity-modulated radiation therapy (row b), NSCLC patients treated with stereotactic body RT (row c), NSCLC patients treated with passive-scattering proton therapy (row d), and the whole cohort of patients (row e).
Figure 1. Axial (columns 1 and 3) and coronal (columns 2 and 4) CT views of the reference anatomy fused with mean (columns 1 and 2) and standard deviations (columns 3 and 4) of BED maps computed voxel by voxel over the patients for the Hodgkin lymphoma group (row a), non-small cell lung cancer (NSCLC) patients treated with intensity-modulated radiation therapy (row b), NSCLC patients treated with stereotactic body RT (row c), NSCLC patients treated with passive-scattering proton therapy (row d), and the whole cohort of patients (row e).
Cancers 13 03553 g001
Figure 2. Plots of C M μ (a) and C M σ (b) as function of the fractional volume V for the different groups of patients and the entire cohort.
Figure 2. Plots of C M μ (a) and C M σ (b) as function of the fractional volume V for the different groups of patients and the entire cohort.
Cancers 13 03553 g002
Figure 3. Probabilistic independent component analysis (PICA) order estimation for the analyzed groups of patients. (a,b) show the explained variance and the model evidence as a function of the included number of components; (c) shows the eigenvalue of each component.
Figure 3. Probabilistic independent component analysis (PICA) order estimation for the analyzed groups of patients. (a,b) show the explained variance and the model evidence as a function of the included number of components; (c) shows the eigenvalue of each component.
Cancers 13 03553 g003
Figure 4. Probabilistic independent component analysis (PICA) of the dose distributions of Hodgkin lymphoma patients (first row), non-small cell lung cancer (NSCLC) patients treated with intensity-modulated radiation therapy (second row), NSCLC patients treated with stereotactic body RT (third row), NSCLC patients treated with passive-scattering proton therapy (fourth row), and the whole cohort of patients (fifth row). For each row, the first four PICA components are displayed anchored to the underlying patient anatomy, as well as the axial (fifth column) and coronal (sixth column) views of the PICA component patches exceeding the 95% confidence interval of the normal distribution (sets of patch intersections marked with the same color are connected in the 3D domain).
Figure 4. Probabilistic independent component analysis (PICA) of the dose distributions of Hodgkin lymphoma patients (first row), non-small cell lung cancer (NSCLC) patients treated with intensity-modulated radiation therapy (second row), NSCLC patients treated with stereotactic body RT (third row), NSCLC patients treated with passive-scattering proton therapy (fourth row), and the whole cohort of patients (fifth row). For each row, the first four PICA components are displayed anchored to the underlying patient anatomy, as well as the axial (fifth column) and coronal (sixth column) views of the PICA component patches exceeding the 95% confidence interval of the normal distribution (sets of patch intersections marked with the same color are connected in the 3D domain).
Cancers 13 03553 g004
Figure 5. Connectograms in heart and lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LLL: left lung lobe; LUL: left upper lobe; RLL: right lower lobe; RML: right middle lobe; RUL: right upper lobe; LA: left atrium; LV: left ventricle; RA: right atrium; RV: right ventricle; PC: pericardium.
Figure 5. Connectograms in heart and lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LLL: left lung lobe; LUL: left upper lobe; RLL: right lower lobe; RML: right middle lobe; RUL: right upper lobe; LA: left atrium; LV: left ventricle; RA: right atrium; RV: right ventricle; PC: pericardium.
Cancers 13 03553 g005
Figure 6. Connectograms in heart substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LA: left atrium; LAW: left atrium wall; LV: left ventricle; LVW: left ventricle wall; RA: right atrium; RAW: right atrium wall; RV: right ventricle; RVW: right ventricle wall; PC: pericardium.
Figure 6. Connectograms in heart substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LA: left atrium; LAW: left atrium wall; LV: left ventricle; LVW: left ventricle wall; RA: right atrium; RAW: right atrium wall; RV: right ventricle; RVW: right ventricle wall; PC: pericardium.
Cancers 13 03553 g006
Figure 7. Connectograms in left lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LLAM: left lung lower lobe anteromedial; LLL: left lung lower lobe lateral; LLP: left lung lower lobe posterior; LLS: left lung lower lobe superior; IL: inferior lingula; SL: superior lingula; LUA: left lung upper lobe anterior; LUAP: left lung upper lobe apicoposterior.
Figure 7. Connectograms in left lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: LLAM: left lung lower lobe anteromedial; LLL: left lung lower lobe lateral; LLP: left lung lower lobe posterior; LLS: left lung lower lobe superior; IL: inferior lingula; SL: superior lingula; LUA: left lung upper lobe anterior; LUAP: left lung upper lobe apicoposterior.
Cancers 13 03553 g007
Figure 8. Connectograms in right lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: RLA: right lung lower lobe anterior; RLL: right lung lower lobe lateral; RLM: right lung lower lobe medial; RLP: right lung lower lobe posterior; RLS: right lung lower lobe superior; RML: right lung middle lobe lateral; RMM: right lung middle lobe medial; RUA: right lung upper lobe anterior; RUAP: right lung upper lobe apical; RUP: right lung upper lobe posterior.
Figure 8. Connectograms in right lung substructures for the four groups of patients: Hodgkin lymphoma group—3DCRT (a), non-small cell lung cancer patients treated with intensity-modulated radiation therapy—IMRT (b), NSCLC patients treated with stereotactic body RT—SBRT (c), NSCLC patients treated with passive-scattering proton therapy—PSPT (d). The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: RLA: right lung lower lobe anterior; RLL: right lung lower lobe lateral; RLM: right lung lower lobe medial; RLP: right lung lower lobe posterior; RLS: right lung lower lobe superior; RML: right lung middle lobe lateral; RMM: right lung middle lobe medial; RUA: right lung upper lobe anterior; RUAP: right lung upper lobe apical; RUP: right lung upper lobe posterior.
Cancers 13 03553 g008
Figure 9. Connectograms for the whole cohort of patients in different sets of substructures: (a) heart and lung substructures, (b) left lung substructures, (c) heart substructures, (d) right lung substructures. The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: (a) LLL: left lung lobe; LUL: left upper lobe; RLL: right lower lobe; RML: right middle lobe; RUL: right upper lobe; LA: left atrium; LV: left ventricle; RA: right atrium; RV: right ventricle; PC: pericardium. (c) LA: left atrium; LAW: left atrium wall; LV: left ventricle; LVW: left ventricle wall; RA: right atrium; RAW: right atrium wall; RV: right ventricle; RVW: right ventricle wall; PC: pericardium. (b) LLAM: left lung lower lobe anteromedial; LLL: left lung lower lobe lateral; LLP: left lung lower lobe posterior; LLS: left lung lower lobe superior; IL: inferior lingula; SL: superior lingula; LUA: left lung upper lobe anterior; LUAP: left lung upper lobe apicoposterior. (d) RLA: right lung lower lobe anterior; RLL: right lung lower lobe lateral; RLM: right lung lower lobe medial; RLP: right lung lower lobe posterior; RLS: right lung lower lobe superior; RML: right lung middle lobe lateral; RMM: right lung middle lobe medial; RUA: right lung upper lobe anterior; RUAP: right lung upper lobe apical; RUP: right lung upper lobe posterior.
Figure 9. Connectograms for the whole cohort of patients in different sets of substructures: (a) heart and lung substructures, (b) left lung substructures, (c) heart substructures, (d) right lung substructures. The pairwise Spearman correlation coefficients between the mean biological effective dose (BED) values to the reported substructures are represented by the lines within the rings. From inside to outside, the rings represent the average of the substructure mean BEDs, standard deviation of the substructure mean BEDs, and average of the dose standard deviations within the substructure. BED is expressed in Gy. Abbreviations: (a) LLL: left lung lobe; LUL: left upper lobe; RLL: right lower lobe; RML: right middle lobe; RUL: right upper lobe; LA: left atrium; LV: left ventricle; RA: right atrium; RV: right ventricle; PC: pericardium. (c) LA: left atrium; LAW: left atrium wall; LV: left ventricle; LVW: left ventricle wall; RA: right atrium; RAW: right atrium wall; RV: right ventricle; RVW: right ventricle wall; PC: pericardium. (b) LLAM: left lung lower lobe anteromedial; LLL: left lung lower lobe lateral; LLP: left lung lower lobe posterior; LLS: left lung lower lobe superior; IL: inferior lingula; SL: superior lingula; LUA: left lung upper lobe anterior; LUAP: left lung upper lobe apicoposterior. (d) RLA: right lung lower lobe anterior; RLL: right lung lower lobe lateral; RLM: right lung lower lobe medial; RLP: right lung lower lobe posterior; RLS: right lung lower lobe superior; RML: right lung middle lobe lateral; RMM: right lung middle lobe medial; RUA: right lung upper lobe anterior; RUAP: right lung upper lobe apical; RUP: right lung upper lobe posterior.
Cancers 13 03553 g009
Figure 10. Coronal view of the XCAT computed tomography fused with: (a) the voxel-wise mean (μ) of the biologically effective dose (BED, in Gy) maps; (b) generalized linear model coefficient (in kGy−1) associated with BED for the development of radiation pneumonitis; (c) the voxel-wise standard deviation (σ) of BED maps (in Gy); (d) significance of the BED coefficient, expressed as −log p.
Figure 10. Coronal view of the XCAT computed tomography fused with: (a) the voxel-wise mean (μ) of the biologically effective dose (BED, in Gy) maps; (b) generalized linear model coefficient (in kGy−1) associated with BED for the development of radiation pneumonitis; (c) the voxel-wise standard deviation (σ) of BED maps (in Gy); (d) significance of the BED coefficient, expressed as −log p.
Cancers 13 03553 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Palma, G.; Monti, S.; Pacelli, R.; Liao, Z.; Deasy, J.O.; Mohan, R.; Cella, L. Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis. Cancers 2021, 13, 3553. https://doi.org/10.3390/cancers13143553

AMA Style

Palma G, Monti S, Pacelli R, Liao Z, Deasy JO, Mohan R, Cella L. Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis. Cancers. 2021; 13(14):3553. https://doi.org/10.3390/cancers13143553

Chicago/Turabian Style

Palma, Giuseppe, Serena Monti, Roberto Pacelli, Zhongxing Liao, Joseph O. Deasy, Radhe Mohan, and Laura Cella. 2021. "Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis" Cancers 13, no. 14: 3553. https://doi.org/10.3390/cancers13143553

APA Style

Palma, G., Monti, S., Pacelli, R., Liao, Z., Deasy, J. O., Mohan, R., & Cella, L. (2021). Radiation Pneumonitis in Thoracic Cancer Patients: Multi-Center Voxel-Based Analysis. Cancers, 13(14), 3553. https://doi.org/10.3390/cancers13143553

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