- Open Access
18F-FET PET radiomics-based survival prediction in glioblastoma patients receiving radio(chemo)therapy
Radiation Oncology volume 17, Article number: 198 (2022)
Quantitative image analysis based on radiomic feature extraction is an emerging field for survival prediction in oncological patients. 18F-Fluorethyltyrosine positron emission tomography (18F-FET PET) provides important diagnostic and grading information for brain tumors, but data on its use in survival prediction is scarce. In this study, we aim at investigating survival prediction based on multiple radiomic features in glioblastoma patients undergoing radio(chemo)therapy.
A dataset of 37 patients with glioblastoma (WHO grade 4) receiving radio(chemo)therapy was analyzed. Radiomic features were extracted from pre-treatment 18F-FET PET images, following intensity rebinning with a fixed bin width. Principal component analysis (PCA) was applied for variable selection, aiming at the identification of the most relevant features in survival prediction. Random forest classification and prediction algorithms were optimized on an initial set of 25 patients. Testing of the implemented algorithms was carried out in different scenarios, which included additional 12 patients whose images were acquired with a different scanner to check the reproducibility in prediction results.
First order intensity variations and shape features were predominant in the selection of most important radiomic signatures for survival prediction in the available dataset. The major axis length of the 18F-FET-PET volume at tumor to background ratio (TBR) 1.4 and 1.6 correlated significantly with reduced probability of survival. Additional radiomic features were identified as potential survival predictors in the PTV region, showing 76% accuracy in independent testing for both classification and regression.
18F-FET PET prior to radiation provides relevant information for survival prediction in glioblastoma patients. Based on our preliminary analysis, radiomic features in the PTV can be considered a robust dataset for survival prediction.
Glioblastoma is the most aggressive primary brain tumor with a median survival of less than two years [1,2,3,4]. Combining surgery, chemotherapy and radiotherapy is the current standard of care, but prognosis remains poor because of the heterogeneous and complex pathogenesis . Therefore, survival prediction based on images acquired before treatment may represent an important step towards a more personalized approach. In recent years, a lot of research has been done in this area with a strong focus on quantitative medical image analysis .
There is an increasing interest towards the use of 18F-Fluorethyltyrosine (18F-FET) Positron Emission Tomography (PET) imaging for the characterization of aggressive brain tumors such as glioblastoma. 18F-FET PET images were reported to provide valuable information, able to effectively support treatment response assessment in glioblastoma [6,7,8,9,10]. In the framework of the recent randomized multicenter phase II trial ARTE, higher 18F-FET PET TBR (tumor-to-background ratio) on pre-treatment scans and persistent 18F-FET PET signal of longer contrast-enhancing tumor was found to associate with inferior overall survival . However, dedicated studies are needed to assess the predictive power related to 18F-FET PET radiomics for survival prediction in glioblastoma. The effective use of radiomics to predict overall survival as an indicator of treatment response comes with a large variability of different strategies regarding image acquisition, tumor segmentation, image processing and feature extraction, data handling and implementation of machine learning models [5, 11, 12]. For example, reliable tumor segmentation methods should reduce the inherent intra- and inter-observer variability . At the same time, segmented contours should be readily available to facilitate their integration in the treatment planning process. In addition, site-specific image acquisition protocols introduce variations in imaging scanners, image reconstruction and even SUV intensities, which can have an effect on the extracted radiomic features [12, 14, 15]. This calls for specific studies to evaluate the reproducibility of radiomic features, in order to maximize their use in survival prediction models.
The objective of this work was to identify key radiomic features in 18F-FET PET images acquired before treatment for glioblastoma patients undergoing radiochemotherapy, and to assess their reproducibility, aiming at consistent prediction accuracy of overall survival.
Radiomic features were extracted from four different volumes of interest (VOIs) to identify the VOIs with the most reliable information in terms of survival prediction. After features were extracted, principal component analysis (PCA) was employed to reduce dimensionality of the feature space, working either as data reduction or feature selection method. Four different datasets were generated with the PCA and used as input data for the machine learning algorithms. Overall, three different prediction models were implemented: (1) a univariate Cox regression model to predict survival time, providing preliminary investigation of potentially relevant features, (2) a random forest classification model to predict overall survival at 1 year, (3) a random forest regression model to predict survival time.
The retrospective study included 37 patients with glioblastoma (WHO grade 4) diagnosed between 2009 and 2017. Survival time was defined, in days, as the time between the day of treatment planning CT acquisition and the day of decease. Due to the high malignancy of glioblastoma and the related poor survival prognosis, the event of decease occurred in all of the cases. The following inclusion criteria were applied for patient selection:
Evidence of macroscopic tumor before the start of treatment planning
PET image acquired not more than 3 weeks before the treatment planning CT to avoid significant anatomical changes and tumor growth/shrinkage between PET imaging and start of treatment
IDH-wildtype glioblastoma, in cases of available IDH mutational status
Follow-up of at least one year.
18F-FET PET images were acquired at the Department of Nuclear Medicine, University Hospital of Munich prior to primary radiation treatment. 18F-FET PET images were used to delineate the tumor based on standardized uptake values (SUVs) and subsequently extract radiomic features within the selected volumes of interest (VOIs). Planning target volumes (PTVs) were delineated within the treatment planning software Oncentra External Beam (version 4.5, Nucletron, Elekta AB, Sweden). When available (32 out of 37 patients), the reconstructed treatment planning CTs and structures were restored from the archive using the treatment planning software.
18F-FET PET images acquired with a Siemens ECAT Exact HR+ (Siemens Healthineers, Erlangen, Germany) scanner according to the standard protocol of the department, described in detail previously , were used as initial training set (N = 25 cases). A transmission scan was performed during 15 min with a 68Ge rod source. Thereafter, approximately 180 MBq of 18F-FET were injected as an intravenous bolus. Once the patient was in position, a 40 min dynamic emission recording in 3D mode consisting of 16 frames was started. The ordered subset expectation maximization (OSEM) algorithm was used for image reconstruction slice per slice (OSEM-2D). The OSEM-2D reconstruction consisted of 6 iterations and used 16 projection subsets. The voxel size in images reconstructed with the described OSEM-2D algorithm was 2.035 × 2.035 × 2.425 mm3, with a resulting image matrix of size 128 × 128 × 63.
The remaining 12 patients were acquired in list mode with a Biograph 64 PET/CT (Siemens Healthineers, Erlangen, Germany), and OSEM-3D was used to reconstruct images. OSEM-3D images were reconstructed with 4 iterations, 21 subsets, a 5 mm Gaussian post-reconstruction filter and standard corrections (attenuation, random and scattered coincidences, dead time decay). The reconstruction matrix was set to 336 × 336 × 109, with voxel size 1.018 × 1.018 × 2.027 mm3.
The key features of the available clinical dataset are summarized in Table 1.
Four different volumes of interest (VOIs) were used to compare the prognostic power of the radiomics analysis when limited to a specific region. The PMOD software (PMOD Technologies, Zurich, Switzerland) was used to generate PET-based VOIs and to register PTV contours to the PET images. PET-based VOIs were defined using a semi-automatic segmentation method, as described by Unterrainer et al. . A tumor-to-background ratio (TBR) threshold of 1.6 has been considered to be optimal to differentiate tumor tissue and surrounding healthy tissue [16, 17]. Besides 1.6, two additional TBR thresholds were chosen, one smaller value of 1.4 and one larger value of 1.8. This procedure resulted in a set of three PET-based VOIs for every patient, further on referred to as voi14, voi16 and voi18 corresponding to the TBR threshold 1.4, 1.6 and 1.8, respectively.
The fourth VOI was generated using the PTV, which was defined by experienced radiation oncologists during the treatment planning process. The gross tumor volume (GTV) was contoured based on the contrast enhanced T1 weighted MRI, complementing the treatment planning CT. T2 and Fluid Attenuated Inversion Recovery (FLAIR) MRI images were used in addition if available. The clinical target volume (CTV) was then created by adding a 2 cm margin to the GTV . The final planning target volume (PTV) was obtained by adding 0.3–0.5 cm to the CTV. As the goal was to extract radiomics features from the PET image, the PTV structures had to be registered to the PET images. The PET image was loaded into PMOD, serving as reference image, whereas the treatment planning CT image was iteratively registered via rigid matching. The resulting transformation matrix was applied to the PTV contour to be overlaid onto the PET image. This provided PET intensity values within the PTV with no need to interpolate raw image intensities. The workflow of the image registration procedure is depicted in Fig. 1.
The PTV, when defined on the treatment planning CT, may include bone and other structures with high CT numbers that do not provide any information about tracer uptake in tumor tissue and surrounding brain tissue. As those regions would cause noise in the radiomic features and impair identification of intensity related patterns in tumor regions, these structures were excluded from the PTV. As shown schematically in Fig. 1, the PTV was much larger compared to the metabolically active tumor, consistently with previously reported studies . This provided a way to assess the role of areas featuring a low intensity 18F-FET PET signal, included in the PTV but outside the metabolically active area (TBR-based VOIs).
Image intensities were rebinned using a fixed bin width approach [11, 20]. To determine the appropriate bin width, the mean TBR intensity range of the initial 25 patients was calculated from voi14 and subsequently divided by the number of bins. Such a number was selected aiming at optimal feature robustness, as investigated by Tsujikawa et al. for PET imaging , and therefore set to 64. According to this strategy, a bin width of 0.04 was obtained for TBR-based resampling, by dividing the mean TBR intensity range (voi14) in 64 bins.
Following rebinning, the python package pyradiomics was used to extract quantitative image features for the selected VOIs (i.e. voi14, voi16, voi18 and PTV) in the 18F-FET PET image [22, 23]. In total 107 radiomic features were extracted from the 3D PET images, including the following feature classes: First Order (FO), Shape (SH), Gray Level Co-occurrence Matrix (GLCM), Gray Level Size Zone Matrix (GLSZM), Gray Level Run Length Matrix (GLRLM), Neighbouring Gray Tone Difference Matrix (NGTDM), Gray Level Dependence Matrix (GLDM).
Dimensionality reduction and feature selection
Unsupervised methods were applied to reduce the data dimensionality and perform feature selection. For both tasks we applied PCA based on the singular value decomposition (SVD) of the feature set, as provided within the scikit-learn package (version 0.22.2) [24,25,26]. PCA was utilized following the standardization of radiomic features, which resulted in feature values with zero mean and unit variance.
Four different types of feature datasets were considered:
Top: the dataset consisted of 10 features per principal component and was generated using PCA-based feature selection: the features with the highest PCA loadings were retained, focusing on the components that explained 90–95% of the total variance in the data
PC: the dataset consisted of the principal component values of each patient, calculated with the entire feature set of 107 features, aiming at the same total variance applied for the Top features
SH-FO: the dataset consisted of all the shape and first order features (32 features overall)
PC(SH-FO): the dataset consisted of the principal component values of each patient, calculated with shape and first order features (32 overall), with the same goal in terms of total variance explained.
Identification of relevant radiomic features
The extracted radiomic features were analyzed to highlight the relevant ones for the purpose of survival prediction.
The Cox proportional hazard model was applied to identify relevant radiomic features in survival prediction based on the calculated hazard ratios. The Cox regression was implemented using the lifelines package in python (version 0.24.4) . Univariate Cox models were fitted to the four feature datasets and potentially relevant features were identified based on the reported p-value of the Wald statistics. A threshold corresponding to 95% significance (i.e. 5% p-value) was selected to classify whether the corresponding radiomic features had a statistically relevant impact on survival prediction.
Survival prediction based on multiple features
We used random forests for both classification and regression analysis based on multiple radiomic features. Random forest is an ensemble method, which means that it combines a number of base estimators built with decision trees in order to reduce variance and improve generalizability [28, 29]. Random forest algorithms were implemented with the python package Scikit-learn . Survival prediction was applied for classification (categorize survival at one-year follow-up) and regression (prediction of survival time). In Scikit-learn, the random forest algorithm for classification provided the use of different class weights to adjust for the imbalance in class frequencies.
To select the best hyper-parameter settings in the models, a coarse to fine strategy was applied using a randomized search cross-validation on training data. To reduce computation time, the randomized search cross-validation was applied twice: (1) in a very broad range of hyper-parameters to identify the order of magnitude of the best performing hyper-parameters and (2) once the magnitude was assessed, multiple values within this magnitude were included in the final parameter grid.
The performance in survival prediction was quantified on the optimal model (i.e. with the determined hyperparameters): we rated survival classification at 12 months based on the Area Under the Receiver Operating Characteristic Curve (AUC), and survival time prediction based on the Concordance Index (C-index). We considered different testing scenarios, where data from the 25 OSEM-2D cases were mixed with the 12 OSEM-3D cases:
Only OSEM-3D data in the test dataset, only OSEM-2D data in the training dataset
A fivefold train and test cross-validation using the entire dataset (OSEM-2D and OSEM-3D data combined).
Patients and treatment
Patients were treated with conventionally fractionated radiochemotherapy with temozolomide (TMZ) based on the EORTC26981/22981-NCIC CE3 protocol  in 34 cases, in one case with conventionally fractionated radiochemotherapy with TMZ and lomustine according to the NOA-09 protocol , in one case with hypofractionated radiotherapy only and in one case with conventionally fractionated radiotherapy only due to a contraindication for chemotherapy. Radiotherapy dosage was 30 × 2 Gy in 32 cases, 29 × 2 Gy due to discontinuation in one case, 33 × 1.8 Gy in two cases, 12 × 1.8 Gy/2.1 Gy on a simultaneous integrated boost (SIB) due to discontinuation in one case and with 10 × 3 Gy plus 1 × 4 Gy in one case.
Identification of relevant radiomic features
When considering the initial 25 patients, a total number of 5–7 principal components turned out to be sufficient to explain 90–95% of the total variance in the data. The Top dataset included the 10 features with the highest loadings for the first 5–7 principal components: Fig. 2 shows the features that appeared most frequently in the Top dataset across the 4 VOIs (voi14, voi16, voi18 and PTV). Shape and first order features were predominantly selected (87.5% and 55.6% average selection frequency, respectively) compared to all other feature types (37.3% average selection frequency). This indicates that shape and first order features have larger variations across the analyzed glioblastoma patients, which might be relevant for outcome prediction.
When analyzing the statistical relevance of single radiomic features via the Cox proportional hazard model, very few radiomic features turned out to be significant in the Top dataset. Figure 3 depicts graphically which radiomic features were classified as significant in the different regions at 95% confidence in univariate Cox proportional hazard models.
Figure 4 shows the effect on survival of the shape feature Major Axis Length in voi14 and voi16 which resulted in significantly different prediction performance for the univariate Cox proportional hazard model. The log-rank test comparison between the Kaplan–Meier curves, where patients were stratified according to the median value, exhibited a p-value smaller than 0.05. This indicates the ability of the radiomic feature Major Axis Length to stratify patients in different risk classes, connected to a significantly different probability of survival at 95% statistical confidence.
Survival prediction based on multiple features
Table 2 reports the average results for the two testing scenarios. The best performance was obtained when considering the PTV as VOI for independent testing on OSEM-3D images. Conversely, results for the CV case were in general lower than independent testing. This indicates that features selected on OSEM-2D images can adequately predict survival in OSEM-3D, whereas data mixing in training does not allow the implemented models to achieve satisfactory prediction performance.
We analyzed the reproducibility of radiomic features extracted from 18F-FET PET images, and their use for survival analysis in glioblastoma patients. Two different patient cohorts were available which included cases with PET images reconstructed via an OSEM-2D algorithm, as well as more recent cases with OSEM-3D reconstruction. The OSEM-2D cohort served as a homogeneous training dataset, providing a consistent source of data to optimize the performance of survival prediction models based on multiple radiomic features. The random forests survival models, featuring the parameters optimized on the first patient cohort, were then tested in multiple scenarios including newer data. This provided a framework to check the reproducibility of radiomic features extracted from 18F-FET PET images to achieve accurate survival prediction in glioblastoma patients.
When analyzing the role of single radiomic features, shape and first order features extracted from 18F-FET PET images were highlighted as the most relevant in the analyzed dataset. This was confirmed both in unsupervised PCA feature selection (Fig. 2) and in the analysis of hazard ratios in univariate Cox proportional models for PET-based VOIs (Fig. 3). As visible in Fig. 3, different features were deemed significant when analyzing the univariate Cox proportional models in multiple regions (VOIs). Univariate Cox models were used to assess significant variables, instead of multivariate analysis. This choice was justified by the intrinsic limitations in Cox proportional models for reduced number of events per variable, as discussed in .
A closer look at these results highlights that increasing values of the shape feature Major Axis Length (i.e. the maximal extension of the biologically active area in the 18F-FET PET image) in voi14 and voi16 are significantly correlated to a reduced overall survival time. This means that larger regions of biologically active tumor in the 18F-FET PET image, when the image is thresholded at 1.4 or 1.6 TBR, are clear indicators of lower survival probability. When extending the analysis to the PTV, the Elongation was also found to be significantly correlated to survival. The Elongation is defined as the square root of the ratio between the minor and major axis of the PTV, with values ranging between 0 (maximal elongation) and 1 (sphere-like shape) . As increasing values are associated to poorer survival, this means that spherical PTVs seem less favorable in terms of survival compared to elongated PTVs. It should be stressed here that Elongation values in our dataset were relatively high (above 0.57), meaning that the minor axis of the PTV is at least 33% of the major axis. Specific further studies are therefore required to check performance in case of elongated PTVs, i.e. corresponding to small Elongation values.
Higher order features were not linked to significant survival prediction in TBR-based VOIs. Conversely, increasing values of GLCM Joint Entropy in the PTV, which encodes variability in intensity, were shown to decrease survival. The existence of multiple significant features calls for survival prediction models that can better handle multiple radiomics signatures simultaneously, such as the implemented random forest algorithms.
The generalization capabilities were analyzed by investigating the reproducibility of radiomics-based predictions in two testing scenarios, where images reconstructed with the OSEM-3D algorithm were also included. In independent testing, significantly better results were obtained when considering the PTV as VOI (Table 2). In this setting, feature sets that included shape and first order 18F-FET PET features obtained more consistent results, thus confirming the findings of univariate Cox models. The Cox regression models resulted in general lower performance than the implemented Random Forest predictors. When testing the univariate Cox models on OSEM-3D data, the C-index for the Major Axis Length in VOI16 and VOI14 reached 0.58 and 0.56, respectively.
Conversely, the performance in CV when mixing different images (OSEM-2D vs. OSEM-3D) were inferior, with the exception of the PC(SH-FO) dataset for survival classification in the PTV (0.7 AUC). This can be partly explained by the lower degree of reproducibility between VOIs segmented in the OSEM-2D reconstructed images versus the OSEM-3D ones. Results from independent testing on OSEM-3D images indicated that features selected on OSEM-2D images can adequately predict survival on different images. Overall best testing performance was found for the SH-FO dataset in classification for the PTV, thus confirming the importance of shape and first order features.
The optimal performance in the PTV shows that also areas of the 18F-FET PET image at lower image intensities, i.e. outside of the high intensity region at 1.4–1.6 TBR, may play a key role in radiomics-based survival prediction. This also strengthens the relevance of expert-based contouring for PTV delineation, which can potentially highlight relevant areas of low/medium biological activity, due to the intrinsic relevance in survival prediction.
A comparative analysis of our results with other studies is challenging, due to the lack of extensive literature on the topic. 18F-FET PET images have been mostly used for diagnostic purposes, and/or to differentiate progression and treatment-related changes [5, 9]. Carles et al.  have reported an analysis of relevant 18F-FET PET features in recurrent glioblastoma patients, including also overall survival as clinical endpoint. More complex features were highlighted as significant for recurrent glioblastoma patients, whereas our results suggest that shape and first order features are more predictive for overall survival in primary glioblastoma. Conversely, best results in our study were achieved for the PTV, where the GLCM Joint Entropy turned out to be significant in univariate analysis. This confirms that higher order features are indeed relevant for primary glioblastoma patients, especially when considered within the PTV region.
MR images have been used more frequently to predict survival in a radiomics-based approach. The performance reported in our study based on 18F-FET PET images are in line with previously reported outcomes on MRI-based radiomics in glioblastoma [33,34,35,36,37,38]. These studies include also more recent extensions from traditional radiomics to deep learning models, designed to include traditional risk factors as potential predictors . Therefore, 18F-FET PET images are suggested as valuable candidates for the prediction of overall survival in primary glioblastoma patients. The complimentary use of MRI and 18F-FET PET images has been recently shown in recurrent glioblastoma patients to derive independent biomarkers of response to treatment .
The main limitation in our study is the reduced number of patients, also in reason of the unique characteristics of the analyzed patient cohort, which consists solely of glioblastoma patients receiving 18F-FET PET prior to primary radio(chemo)therapy. This motivated the use of stringent criteria for feature selection, but hindered a more extensive analysis of the relevance of single features in survival prediction. For this reason, the univariate analysis was restricted to the Top features dataset, which might have underestimated the relevance of complex textural features, as these latter were under-represented in the PCA feature selection procedure (Fig. 2). Reported results indeed confirm the appropriateness of such a selection, as PCA including only SH-FO features outperformed PCA including more complex features (Table 2). As further limitation, no comparison with MRI-based radiomic features has been carried out in this work. Different studies have investigated the use of quantitative features from MRI and 18F-FET PET for treatment response and survival [6, 10, 38], highlighting the complimentary role of these imaging modalities. A comparative radiomics-based analysis is therefore suggested to highlight the need of dedicated imaging protocols for glioblastoma survival prediction.
The findings of feature analysis and prediction performance based on 18F-FET PET images show that shape features in voi14, voi16, and PTV are significant predictors of overall survival in glioblastoma patients. When radiomics-based prediction is enlarged to the PTV a more robust prediction of survival is achieved, with 76% accuracy in independent testing both for classification and regression analyses. These preliminary results should be extended to a larger patient cohort for further validation of our findings.
Availability of data and materials
Due to privacy and ethical concerns, the source data cannot be made available.
Area under the curve
Fluid attenuated inversion recovery
Gray level co-occurrence matrix
Gray level dependence matrix
Gray level run length matrix
Gray level size zone matrix
Mean absolute error
Magnetic resonance imaging
Mean square error
Neighbouring gray tone difference matrix
Ordered subset expectation maximization
Principal component analysis
Positron emission tomography
Planning target volume
Standardized uptake value
Tumor to BackgRound
Volume of interest
Wirsching HG, Galanis E, Weller M. Glioblastoma. In: Berger MS, Weller M, editors. Handbook of clinical neurology, vol. 134. Amsterdam: Elsevier; 2016. p. 381–97.
Hanif F, Muzaffar K, Perveen K, Malhi SM, Simjee SU. Glioblastoma multiforme: a review of its epidemiology and pathogenesis through clinical presentation and treatment. Asian Pac J Cancer Prev. 2017;18(1):3–9.
Tamimi AF, Juweid M. Epidemiology and outcome of glioblastoma. In: De Vleeschouwer S, editor. Glioblastoma [internet]. Brisbane: Codon Publications; 2017. https://doi.org/10.15586/codon.glioblastoma.2017.ch8.
Fernandes C, Costa A, Osório L, et al. Current standards of care in glioblastoma therapy. In: De Vleeschouwer S, editor., et al., Glioblastoma [internet]. Brisbane: Codon Publications; 2017. https://doi.org/10.15586/codon.glioblastoma.2017.ch11.
Kocher M, Ruge MI, Galldiks N, Lohmann P. Applications of radiomics and machine learning for radiotherapy of malignant brain tumors. Strahlenther Onkol. 2020;196(10):856–67.
Galldiks N, Dunkl V, Ceccon G, Tscherpel C, Stoffels G, Law I, et al. Early treatment response evaluation using FET PET compared to MRI in glioblastoma patients at first progression treated with bevacizumab plus lomustine. Eur J Nucl Med Mol Imaging. 2018;45(13):2377–86.
Pyka T, Gempt J, Hiob D, Ringel F, Schlegel J, Bette S, et al. Textural analysis of pre-therapeutic [18F]-FET-PET and its correlation with tumor grade and patient survival in high-grade gliomas. Eur J Nucl Med Mol Imaging. 2016;43(1):133–41.
Carles M, Popp I, Starke MM, Mix M, Urbach H, Schimek-Jasch T, Eckert F, Niyazi M, Baltas D, Grosu AL. FET-PET radiomics in recurrent glioblastoma: prognostic value for outcome after re-irradiation? Radiat Oncol. 2021;16(1):46.
Suchorska B, Jansen NL, Linn J, Kretzschmar H, Janssen H, Eigenbrod S, Simon M, Pöpperl G, Kreth FW, la Fougere C, Weller M, Tonn JC, Network GG. Biological tumor volume in 18FET-PET before radiochemotherapy correlates with survival in GBM. Neurology. 2015;84(7):710–9.
Wirsching HG, Roelcke U, Weller J, Hundsberger T, Hottinger AF, von Moos R, Caparrotti F, Conen K, Remonda L, Roth P, Ochsenbein A, Tabatabai G, Weller M. MRI and 18FET-PET predict survival benefit from bevacizumab plus radiotherapy in patients with isocitrate dehydrogenase wild-type glioblastoma: results from the randomized ARTE trial. Clin Cancer Res. 2021;27(1):179–88.
Leijenaar RTH, Nalbantov G, Carvalho S, Van Elmpt WJC, Troost EGC, Boellaard R, et al. The effect of SUV discretization in quantitative FDG-PET radiomics: the need for standardized methodology in tumor texture analysis. Sci Rep. 2015;5(5):11075.
Altazi BA, Zhang GG, Fernandez DC, Montejo ME, Hunt D, Werner J, et al. Reproducibility of F18-FDG PET radiomic features for different cervical tumor segmentation methods, gray-level discretization, and reconstruction algorithms. J Appl Clin Med Phys. 2017;18(6):32–48.
Lee M, Woo B, Kuo MD, Jamshidi N, Kim JH. Quality of radiomic features in glioblastoma multiforme: impact of semi-automated tumor segmentation software. Korean J Radiol. 2017;18(3):498–509.
Liu Z, Wang S, Dong D, Wei J, Fang C, Zhou X, et al. The applications of radiomics in precision diagnosis and treatment of oncology: opportunities and challenges. Theranostics. 2019;9(5):1303–22.
Rizzo S, Botta F, Raimondi S, Origgi D, Fanciullo C, Morganti AG, et al. Radiomics: the facts and the challenges of image analysis. Eur Radiol Exp. 2018;2(1):36.
Unterrainer M, Vettermann F, Brendel M, Holzgreve A, Lifschitz M, Zähringer M, et al. Towards standardization of 18F-FET PET imaging: do we need a consistent method of background activity assessment? EJNMMI Res. 2017;7:1–8.
Vomacka L, Unterrainer M, Holzgreve A, Mille E, Gosewisch A, Brosch J, et al. Voxel-wise analysis of dynamic 18F-FET PET: a novel approach for non-invasive glioma characterisation. EJNMMI Res. 2018;8(1):91.
Niyazi M, Brada M, Chalmers AJ, Combs SE, Erridge SC, Fiorentino A, Grosu AL, Lagerwaard FJ, Minniti G, Mirimanoff RO, Ricardi U, Short SC, Weber DC, Belka C. ESTRO-ACROP guideline “target delineation of glioblastomas.” Radiother Oncol. 2016;118(1):35–42.
Fleischmann DF, Unterrainer M, Schön R, Corradini S, Maihöfer C, Bartenstein P, Belka C, Albert NL, Niyazi M. Margin reduction in radiotherapy for glioblastoma through 18F-fluoroethyltyrosine PET? A recurrence pattern analysis. Radiother Oncol. 2020;145:49–55.
Branchini M, Zorz A, Zucchetta P, Bettinelli A, De Monte F, Cecchin D, et al. Impact of acquisition count statistics reduction and SUV discretization on PET radiomic features in pediatric 18F-FDG-PET/MRI examinations. Phys Med. 2019;59:117–26.
Tsujikawa T, Tsuyoshi H, Kanno M, Yamada S, Kobayashi M, Narita N, et al. Selected PET radiomic features remain the same. Oncotarget. 2018;9(29):20734–46.
PyRadiomics Radiomic Features. https://pyradiomics.readthedocs.io/en/latest/features.html.
van Griethuysen J, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, Beets-Tan R, Fillion-Robin JC, Pieper S, Aerts H. Computational radiomics system to decode the radiographic phenotype. Can Res. 2017;77(21):e104–7.
Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12(Oct):2825–30.
Geron A. Dimensionality reduction. In: Géron A, editor. Hands-on machine learning with scikit-learn and TensorFlow. 1st ed. Beijing: O’Reilly Media; 2017. p. 207–27.
sklearn.decomposition.PCA | scikit-learn. https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.PCA.html.
Davidson-Pilon C, Kalderstam J, Jacobson N, Zivich P, Kuhn B, Williamson M, et al. CamDavidsonPilon/lifelines: v0.24.4. 2020. https://lifelines.readthedocs.io/en/latest/
Fratello M, Tagliaferri R. Decision trees and random forests. In: Ranganathan S, Gribskov M, Nakai K, Schoenbach C, editors. Encyclopedia of bioinformatics and computational biology, vol. 1–3. Elsevier: Amsterdam; 2019. p. 374–83.
Ensemble methods | scikit-learn 0.22.2 documentation. https://scikit-learn.org/stable/modules/ensemble.html#forests-of-randomized-trees.
EORTC-NCIC 26981/22981. https://www.eortc.org/research_field/clinical-detail/26981/.
Herrlinger U, Tzaridis T, Mack F, Steinbach JP, Schlegel U, Sabel M, Hau P, Kortmann RD, Krex D, Grauer O, Goldbrunner R, Schnell O, Bähr O, Uhl M, Seidel C, Tabatabai G, Kowalski T, Ringel F, Schmidt-Graf F, Suchorska B, Neurooncology Working Group of the German Cancer Society. Lomustine-temozolomide combination therapy versus standard temozolomide therapy in patients with newly diagnosed glioblastoma with methylated MGMT promoter (CeTeG/NOA-09): a randomised, open-label, phase 3 trial. Lancet (London, England). 2019;393(10172):678–88.
Ramia-de-Cap M, Kaul S. Limitations of multivariate survival analysis. Mod Pathol. 2021. https://doi.org/10.1038/s41379-021-00891-x.
Chaddad A, Sabri S, Niazi T, et al. Prediction of survival with multi-scale radiomic analysis in glioblastoma patients. Med Biol Eng Comput. 2018;56:2287–300.
Vils A, Bogowicz M, Tanadini-Lang S, Vuong D, Saltybaeva N, Kraft J, Wirsching HG, Gramatzki D, Wick W, Rushing E, Reifenberger G, Guckenberger M, Weller M, Andratschke N. Radiomic analysis to predict outcome in recurrent glioblastoma based on multi-center MR imaging from the prospective DIRECTOR trial. Front Oncol. 2021;14(11):636672.
Lao J, Chen Y, Li ZC, Li Q, Zhang J, Liu J, Zhai G. A deep learning-based radiomics model for prediction of survival in glioblastoma multiforme. Sci Rep. 2017;7(1):10353.
Bae S, Choi YS, Ahn SS, Chang JH, Kang SG, Kim EH, Kim SH, Lee SK. Radiomic MRI phenotyping of glioblastoma: improving survival prediction. Radiology. 2018;289(3):797–806.
Wang B, Zhang S, Wu X, Li Y, Yan Y, Liu L, Xiang J, Li D, Yan T. Multiple survival outcome prediction of glioblastoma patients based on multiparametric MRI. Front Oncol. 2021;25(11):778627.
Lombardi G, Spimpolo A, Berti S, Campi C, Anglani MG, Simeone R, Evangelista L, Causin F, Zorzi G, Gorgoni G, Caccese M, Padovan M, Zagonel V, Cecchin D. PET/MR in recurrent glioblastoma patients treated with regorafenib: [18F]FET and DWI-ADC for response assessment and survival prediction. Br J Radiol. 2022;95(1129):20211018.
Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
The retrospective analysis was approved by the local ethics committee of the LMU Munich (Study Number 20-255).
Consent for publication
The Department of Radiation Oncology of the University Hospital of LMU Munich has research agreements with Brainlab, Elekta and Viewray.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wiltgen, T., Fleischmann, D.F., Kaiser, L. et al. 18F-FET PET radiomics-based survival prediction in glioblastoma patients receiving radio(chemo)therapy. Radiat Oncol 17, 198 (2022). https://doi.org/10.1186/s13014-022-02164-6
- Quantitative image analysis
- Survival analysis