Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome

Radiomics from Primary Tumor on Dual-Energy CT Derived Iodine Maps can Predict Cervical Lymph Node Metastasis in Papillary Thyroid Cancer

RATIONALE AND OBJECTIVES

To develop and validate 2 iodine maps based radiomics nomograms for preoperatively predicting cervical lymph node metastasis (LNM) and central lymph node metastasis (CLNM) in papillary thyroid cancer (PTC).

MATERIALS AND METHODS

A total of 346 patients with PTC were enrolled and allocated to training (242) and validation (104) sets. Radiomics features were extracted from arterial and venous phase iodine maps, respectively. Aggregated machine-learning strategy was applied for features selection and construction of 2 radiomics scores (LN rad-score; CLN rad-score). Logistic regression model was employed to establish two radiomics nomograms (nomogram 1: predicting LNM; nomogram 2: predicting CLNM) after incorporating LN or CLN rad-score with clinical predictors. Nomograms performance was determined by discrimination, calibration and clinical usefulness.

RESULTS

Nomogram 1 incorporated LN rad-score, age (categorized by 55) and CT reported LN status; Nomogram 2 incorporated CLN rad-score, capsule contact >25% and CT reported CLN status. 2 nomograms both showed good discrimination and calibration in the training (AUC = 0.847; AUC = 0.837) and validation cohorts (AUC = 0.807; AUC = 0.795). Significant improved AUC, net reclassification index (NRI) and integrated discriminatory improvement (IDI) confirmed additional great predictive value of 2 rad-scores, compared with clinical models without radiomics. Decision curve analysis indicated clinical utility of nomograms. 2 nomograms both demonstrated favorable predictive efficacy in CT reported LN or CLN negative subgroup (AUC = 0.766; AUC = 0.744).

CONCLUSION

The presented 2 radiomics nomograms are useful tools for preoperative prediction of LNM and CLNM in PTC.

Clinical Perspectives for 18F-FDG PET Imaging in Pediatric Oncology: Μetabolic Tumor Volume and Radiomics

Pediatric cancer, although rare, requires the most optimized treatment approach to obtain high survival rates and minimize serious long-term side effects in early adulthood. 18F-FDG PET/CT is most helpful and widely used in staging, recurrence detection, and response assessment in pediatric oncology. The well-known 18F-FDG PET metabolic indices of metabolic tumor volume (MTV) and tumor lesion glycolysis (TLG) have already revealed an independent significant prognostic value for survival in oncologic patients, although the corresponding cut-off values remain study-dependent and not validated for use in clinical practice. Advanced tumor “radiomic” analysis sheds new light into these indices. Numerous patterns of texture 18F-FDG uptake features can be extracted from segmented PET tumor images due to new powerful computational systems supporting complex “deep learning” algorithms. This high number of “quantitative” tumor imaging data, although not decrypted in their majority and once standardized for the different imaging systems and segmentation methods, could be used for the development of new “clinical” models for specific cancer types and, more interestingly, for specific age groups. In addition, data from novel techniques of tumor genome analysis could reveal new genes as biomarkers for prognosis and/or targeted therapies in childhood malignancies. Therefore, this ever-growing information of “radiogenomics”, in which the underlying tumor “genetic profile” could be expressed in the tumor-imaging signature of “radiomics”, possibly represents the next model for precision medicine in pediatric cancer management. This paper reviews 18F-FDG PET image segmentation methods as applied to pediatric sarcomas and lymphomas and summarizes reported findings on the values of metabolic and radiomic features in the assessment of these pediatric tumors.

Assessing Tumour Haemodynamic Heterogeneity and Response to Choline Kinase Inhibition Using Clustered Dynamic Contrast Enhanced MRI Parameters in Rodent Models of Glioblastoma

To investigate the utility of DCE-MRI derived pharmacokinetic parameters in evaluating tumour haemodynamic heterogeneity and treatment response in rodent models of glioblastoma, imaging was performed on intracranial F98 and GL261 glioblastoma bearing rodents. Clustering of the DCE-MRI-based parametric maps (using Tofts, extended Tofts, shutter speed, two-compartment, and the second generation shutter speed models) was performed using a hierarchical clustering algorithm, resulting in areas with poor fit (reflecting necrosis), low, medium, and high valued pixels representing parameters Ktrans, ve, Kep, vp, τi and Fp. There was a significant increase in the number of necrotic pixels with increasing tumour volume and a significant correlation between ve and tumour volume suggesting increased extracellular volume in larger tumours. In terms of therapeutic response in F98 rat GBMs, a sustained decrease in permeability and perfusion and a reduced cell density was observed during treatment with JAS239 based on Ktrans, Fp and ve as compared to control animals. No significant differences in these parameters were found for the GL261 tumour, indicating that this model may be less sensitive to JAS239 treatment regarding changes in vascular parameters. This study demonstrates that region-based clustered pharmacokinetic parameters derived from DCE-MRI may be useful in assessing tumour haemodynamic heterogeneity with the potential for assessing therapeutic response.

Correlation of histopathology and multi-modal magnetic resonance imaging in childhood osteosarcoma: Predicting tumor response to chemotherapy

Background

Osteosarcoma, which is the most common malignant pediatric bone cancer, remains dependent on an imprecise systemic treatment largely unchanged in 30 years. In this study, we correlated histopathology with magnetic resonance imaging (MRI), used the correlation to extract MRI-specific features representative of tumor necrosis, and subsequently developed a novel classification model for predicting tumor response to neoadjuvant chemotherapy in pediatric patients with osteosarcoma using multi-modal MRI. The model could ultimately serve as a testable biomarker for a high-risk malignancy without successful precision treatments.

Methods

Patients with newly diagnosed high-grade appendicular osteosarcoma were enrolled in a single-center observational study, wherein patients underwent pre-surgical evaluation using both conventional MRI (post-contrast T1-weighted with fat saturation, pre-contrast T1-weighted, and short inversion-time inversion recovery (STIR)) and advanced MRI (diffusion weighted (DW) and dynamic contrast enhanced (DCE)). A classification model was established based on a direct correlation between histopathology and MRI, which was achieved through histologic-MR image co-registration and subsequent extraction of MR image features for identifying histologic tumor necrosis. By operating on the MR image features, tumor necrosis was estimated from different combinations of MR images using a multi-feature fuzzy clustering technique together with a weighted majority ruling. Tumor necrosis calculated from MR images, for either an MRI plane of interest or whole tumor volume, was compared to pathologist-estimated necrosis and necrosis quantified from digitized histologic section images using a previously described deep learning classification method.

Results

15 patients were enrolled, of whom two withdrew, one became ineligible, and two were subjected to inadequate pre-surgical imaging. MRI sequences of n = 10 patients were subsequently used for classification model development. Different MR image features, depending on the modality of MRI, were shown to be significant in distinguishing necrosis from viable tumor. The scales at which MR image features optimally signified tumor necrosis were different as well depending on the MR image type. Conventional MRI was shown capable of differentiating necrosis from viable tumor with an accuracy averaging above 90%. Conventional MRI was equally effective as DWI in distinguishing necrotic from viable tumor regions. The accuracy of tumor necrosis prediction by conventional MRI improved to above 95% when DCE-MRI was added into consideration. Volume-based tumor necrosis estimations tended to be lower than those evaluated on an MRI plane of interest.

Conclusions

The study has shown a proof-of-principle model for interpreting chemotherapeutic response using multi-modal MRI for patients with high-grade osteosarcoma. The model will continue to be evaluated as MR image features indicative of tumor response are now computable for the disease prior to surgery.

Tumor Prognostic Prediction of Nasopharyngeal Carcinoma Using CT-Based Radiomics in Non-Chinese Patients

Purpose

We aimed to construct predictive models for the overall survival (OS), progression-free survival (PFS), and distant metastasis-free survival (DMFS) for nasopharyngeal carcinoma (NPC) patients by using CT-based radiomics.

Materials and Methods

We collected data from 197 NPC patients. For each patient, radiomic features were extracted from the CT image acquired at pretreatment via PyRadiomics. Feature selection was performed in two steps. First, features with high inter-observer variability based on multiple tumor delineations were excluded. Then, stratified bootstrappings were performed to identify feature combinations that most frequently achieved the highest (i) area under the receiver operating characteristic curve (AUC) for predicting 3-year OS, PFS, and DMFS or (ii) Harrell’s C-index for predicting time to event. Finally, regularized logistic regression and Cox proportional hazard models with the most frequently selected feature combinations as input were tuned using cross-validation. Additionally, we examined the robustness of the constructed model to variation in tumor delineation by simulating 100 realizations of radiomic feature values to mimic features extracted from different tumor boundaries.

Results

The combined model that used both radiomics and clinical features yielded significantly higher AUC and Harrell’s C-index than models using either feature set alone for all outcomes (p < 0.05). The AUCs and Harrell’s C-indices of the clinical-only and radiomics-only models ranged from 0.758 ± 0.091 to 0.789 ± 0.082 and from 0.747 ± 0.062 to 0.767 ± 0.074, respectively. In comparison, the combined models achieved AUC of 0.801 ± 0.075 to 0.813 ± 0.078 and Harrell’s C-indices of 0.779 ± 0.066 to 0.796 ± 0.069. The results showed that our models were robust to variation in tumor delineation with the coefficient of variation ranging from 4.8% to 6.4% and from 6.7% to 9.3% for AUC and Harrell’s C-index, respectively.

Conclusion

Our results demonstrated that using CT-based radiomic features together with clinical features provided superior NPC prognostic prediction than using either clinical or radiomic features alone.

Hypovascular Cellular Tumor in Primary Central Nervous System Lymphoma is Associated with Treatment Resistance: Tumor Habitat Analysis Using Physiologic MRI

BACKGROUND AND PURPOSE

The microenvironment of lymphomas is known to be highly variable and closely associated with treatment resistance and survival. We tried to develop a physiologic MR imaging-based spatial habitat analysis to identify regions associated with treatment resistance to facilitate the prediction of tumor response after initial chemotherapy in patients with primary central nervous system lymphoma.

MATERIALS AND METHODS

Eighty-one patients with pathologically confirmed primary central nervous system lymphoma were enrolled. Pretreatment physiologic MR imaging was performed, and K-means clustering was used to separate voxels into 3 spatial habitats according to ADC and CBV values. Associations of spatial habitats and clinical and conventional imaging predictors with time to progression were analyzed using Cox proportional hazards modeling. The performance of statistically significant predictors for time to progression was assessed using the concordance probability index.

RESULTS

The 3 spatial habitats of hypervascular cellular tumor, hypovascular cellular tumor, and hypovascular hypocellular tumor were identified. A large hypovascular cellular habitat was most significantly associated with short time to progression (hazard ratio, 2.83; P = . 017). The presence of an atypical finding (hazard ratio, 4.41; P = . 016), high performance score (hazard ratio, 5.82; P = . 04), and high serum lactate dehydrogenase level (hazard ratio, 1.01; P = .013) was significantly associated with time to progression. A predictive model constructed using the habitat score and other imaging parameters showed a concordance probability index for prediction of time to progression of 0.70 (95% CI, 0.54-0.87).

CONCLUSIONS

A hypovascular cellular tumor habitat is associated with treatment resistance in primary central nervous system lymphoma, and its assessment may refine prechemotherapy imaging-based response prediction for patients with primary central nervous system lymphoma.

An investigation of the conformity, feasibility, and expected clinical benefits of multiparametric MRI-guided dose painting radiotherapy in glioblastoma

Background

New technologies developed to improve survival outcomes for glioblastoma (GBM) continue to have limited success. Recently, image-guided dose painting (DP) radiotherapy has emerged as a promising strategy to increase local control rates. In this study, we evaluate the practical application of a multiparametric MRI model of glioma infiltration for DP radiotherapy in GBM by measuring its conformity, feasibility, and expected clinical benefits against standard of care treatment.

Methods

Maps of tumor probability were generated from perfusion/diffusion MRI data from 17 GBM patients via a previously developed model of GBM infiltration. Prescriptions for DP were linearly derived from tumor probability maps and used to develop dose optimized treatment plans. Conformity of DP plans to dose prescriptions was measured via a quality factor. Feasibility of DP plans was evaluated by dose metrics to target volumes and critical brain structures. Expected clinical benefit of DP plans was assessed by tumor control probability. The DP plans were compared to standard radiotherapy plans.

Results

The conformity of the DP plans was >90%. Compared to the standard plans, DP (1) did not affect dose delivered to organs at risk; (2) increased mean and maximum dose and improved minimum dose coverage for the target volumes; (3) reduced minimum dose within the radiotherapy treatment margins; (4) improved local tumor control probability within the target volumes for all patients.

Conclusions

A multiparametric MRI model of GBM infiltration can enable conformal, feasible, and potentially beneficial dose painting radiotherapy plans.

OUP accepted manuscript

Background

Materials and Methods

Results

Conclusion

Effects of tracer uptake time in non-small cell lung cancer 18F-FDG PET radiomics

Positron emission tomography (PET) radiomics applied to oncology allows the measurement of intra-tumoral heterogeneity. This quantification can be affected by image protocols hence there is an increased interest in understanding how radiomic expression on PET images is affected by different imaging conditions. To address that, this study explores how radiomic features are affected by changes in ¹⁸F-FDG uptake time, image reconstruction, lesion delineation, and radiomics binning settings. Methods: Ten non-small cell lung cancer (NSCLC) patients underwent ¹⁸F-FDG PET scans on two consecutive days. On each day, scans were obtained at 60min and 90min post-injection and reconstructed following EARL version 1 (EARL1) and with point-spread-function resolution modelling (PSF-EARL2). Lesions were delineated using thresholds at SUV=4.0, 40% of SUV_max, and with a contrast-based isocontour. PET image intensity was discretized with both fixed bin width (FBW) and fixed bin number (FBN) before the calculation of the radiomic features. Repeatability of features was measured with intraclass correlation (ICC), and the change in feature value over time was calculated as a function of its repeatability. Features were then classified on use-case scenarios based on their repeatability and susceptibility to tracer uptake time. Results: With PSF-EARL2 reconstruction, 40% of SUV_max lesion delineation, and FBW intensity discretization, most features (94%) were repeatable at both uptake times (ICC>0.9), 39% being classified for dual-time-point use-case for being sensitive to changes in uptake time, 39% were classified for cross-sectional studies with unclear dependency on time, 20% classified for cross-sectional use while being robust to tracer uptake time changes, and 6% were discarded for poor repeatability. EARL1 images had one less repeatable feature than PSF-EARL2 (Neighborhood Gray-Level Different Matrix Coarseness), the contrast-based delineation had the poorest repeatability of the delineation methods with 45% features being discarded, and FBN resulted in lower repeatability than FBW (45% and 6% features were discarded, respectively). Conclusion: Repeatability was maximized with PSF-EARL2 reconstruction, lesion delineation at 40% of SUV_max, and FBW intensity discretization. Based on their susceptibility to tracer uptake time, radiomic features were classified into specific NSCLC PET radiomics use-cases.

Delta radiomics analysis of Magnetic Resonance guided radiotherapy imaging data can enable treatment response prediction in pancreatic cancer

Background

Magnetic Resonance Image guided Stereotactic body radiotherapy (MRgRT) is an emerging technology that is increasingly used in treatment of visceral cancers, such as pancreatic adenocarcinoma (PDAC). Given the variable response rates and short progression times of PDAC, there is an unmet clinical need for a method to assess early RT response that may allow better prescription personalization. We hypothesize that quantitative image feature analysis (radiomics) of the longitudinal MR scans acquired before and during MRgRT may be used to extract information related to early treatment response.

Methods

Histogram and texture radiomic features (n = 73) were extracted from the Gross Tumor Volume (GTV) in 0.35T MRgRT scans of 26 locally advanced and borderline resectable PDAC patients treated with 50 Gy RT in 5 fractions. Feature ratios between first (F1) and last (F5) fraction scan were correlated with progression free survival (PFS). Feature stability was assessed through region of interest (ROI) perturbation.

Results

Linear normalization of image intensity to median kidney value showed improved reproducibility of feature quantification. Histogram skewness change during treatment showed significant association with PFS (p = 0.005, HR = 2.75), offering a potential predictive biomarker of RT response. Stability analyses revealed a wide distribution of feature sensitivities to ROI delineation and was able to identify features that were robust to variability in contouring.

Conclusions

This study presents a proof-of-concept for the use of quantitative image analysis in MRgRT for treatment response prediction and providing an analysis pipeline that can be utilized in future MRgRT radiomic studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01957-5.

Texture analysis imaging "what a clinical radiologist needs to know"

Texture analysis has arisen as a tool to explore the amount of data contained in images that cannot be explored by humans visually. Radiomics is a method that extracts a large number of features from radiographic medical images using data-characterisation algorithms. These features, termed radiomic features, have the potential to uncover disease characteristics. The goal of both radiomics and texture analysis is to go beyond size or human-eye based semantic descriptors, to enable the non-invasive extraction of quantitative radiological data to correlate them with clinical outcomes or pathological characteristics. In the latest years there has been a flourishing sub-field of radiology where texture analysis and radiomics have been used in many settings. It is difficult for the clinical radiologist to cope with such amount of data in all the different radiological sub-fields and to identify the most significant papers. The aim of this review is to provide a tool to better understand the basic principles underlining texture analysis and radiological data mining and a summary of the most significant papers of the latest years.

Overall Survival Prognostic Modelling of Non-small Cell Lung Cancer Patients Using Positron Emission Tomography/Computed Tomography Harmonised Radiomics Features: The Quest for the Optimal Machine Learning Algorithm

AIMS

Despite the promising results achieved by radiomics prognostic models for various clinical applications, multiple challenges still need to be addressed. The two main limitations of radiomics prognostic models include information limitation owing to single imaging modalities and the selection of optimum machine learning and feature selection methods for the considered modality and clinical outcome. In this work, we applied several feature selection and machine learning methods to single-modality positron emission tomography (PET) and computed tomography (CT) and multimodality PET/CT fusion to identify the best combinations for different radiomics modalities towards overall survival prediction in non-small cell lung cancer patients.

MATERIALS AND METHODS

A PET/CT dataset from The Cancer Imaging Archive, including subjects from two independent institutions (87 and 95 patients), was used in this study. Each cohort was used once as training and once as a test, followed by averaging of the results. ComBat harmonisation was used to address the centre effect. In our proposed radiomics framework, apart from single-modality PET and CT models, multimodality radiomics models were developed using multilevel (feature and image levels) fusion. Two different methods were considered for the feature-level strategy, including concatenating PET and CT features into a single feature set and alternatively averaging them. For image-level fusion, we used three different fusion methods, namely wavelet fusion, guided filtering-based fusion and latent low-rank representation fusion. In the proposed prognostic modelling framework, combinations of four feature selection and seven machine learning methods were applied to all radiomics modalities (two single and five multimodalities), machine learning hyper-parameters were optimised and finally the models were evaluated in the test cohort with 1000 repetitions via bootstrapping. Feature selection and machine learning methods were selected as popular techniques in the literature, supported by open source software in the public domain and their ability to cope with continuous time-to-event survival data. Multifactor ANOVA was used to carry out variability analysis and the proportion of total variance explained by radiomics modality, feature selection and machine learning methods was calculated by a bias-corrected effect size estimate known as ω².

RESULTS

Optimum feature selection and machine learning methods differed owing to the applied radiomics modality. However, minimum depth (MD) as feature selection and Lasso and Elastic-Net regularized generalized linear model (glmnet) as machine learning method had the highest average results. Results from the ANOVA test indicated that the variability that each factor (radiomics modality, feature selection and machine learning methods) introduces to the performance of models is case specific, i.e. variances differ regarding different radiomics modalities and fusion strategies. Overall, the greatest proportion of variance was explained by machine learning, except for models in feature-level fusion strategy.

CONCLUSION

The identification of optimal feature selection and machine learning methods is a crucial step in developing sound and accurate radiomics risk models. Furthermore, optimum methods are case specific, differing due to the radiomics modality and fusion strategy used.

New method for quantification of intratumoral heterogeneity: a feasibility study on Ktrans maps from preclinical DCE-MRI

OBJECT

To develop new imaging biomarkers of therapeutic efficacy through the quantification of intratumoral microvascular heterogeneity.

MATERIALS AND METHODS

The described method was a combination of non-supervised clustering and extraction of intratumoral complexity features (ICF): number of non-connected objects, volume fraction. It was applied to a set of 3D DCE-MRI K^trans maps acquired previously on tumor bearing mice prior to and on day 4 of anti-angiogenic treatment. Evolutions of ICF were compared to conventional summary statistics (CSS) and to heterogeneity related whole tumor texture features (TF) on treated (n = 9) and control (n = 6) mice.

RESULTS

Computed optimal number of clusters per tumor was 4. Several intratumoral features extracted from the clusters were able to monitor a therapy effect. Whereas no feature significantly changed for the control group, 6 features significantly changed for the treated group (4 ICF, 2 CSS). Among these, 5 also significantly differentiated the two groups (3 ICF, 2 CSS). TF failed in demonstrating differences within and between the two groups.

DISCUSSION

ICF are potential imaging biomarkers for anti-angiogenic therapy assessment. The presented method may be expected to have advantages with respect to texture analysis-based methods regarding interpretability of results and setup of standardized image analysis protocols.

[Artificial intelligence and medical imaging]

The use of artificial intelligence methods for image recognition is one of the most developed branches of the AI field and these technologies are now commonly used in our daily lives. In the field of medical imaging, approaches based on artificial intelligence are particularly promising, with numerous applications and a strong interest in the search for new biomarkers. Here, we will present the general methods used in these approaches as well as the potential areas of application.

Early prediction of tumour-response to radiotherapy in NSCLC patients

Objective. In this study we developed an automatic method to predict tumour volume and shape in weeks 3 and 4 of radiotherapy (RT), using cone-beam computed tomography (CBCT) scans acquired up to week 2, allowing identification of large tumour changes.Approach. 240 non-small cell lung cancer (NSCLC) patients, treated with 55 Gy in 20 fractions, were collected. CBCTs were rigidly registered to the planning CT. Intensity values were extracted in each voxel of the planning target volume across all CBCT images from days 1, 2, 3, 7 and 14. For each patient and in each voxel, four regression models were fitted to voxel intensity; applying linear, Gaussian, quadratic and cubic methods. These models predicted the intensity value for each voxel in weeks 3 and 4, and the tumour volume found by thresholding. Each model was evaluated by computing the root mean square error in pixel value and structural similarity index metric (SSIM) for all patients. Finally, the sensitivity and specificity to predict a 30% change in volume were calculated for each model.Main results. The linear, Gaussian, quadratic and cubic models achieved a comparable similarity score, the average SSIM for all patients was 0.94, 0.94, 0.90, 0.83 in week 3, respectively. At week 3, a sensitivity of 84%, 53%, 90% and 88%, and specificity of 99%, 100%, 91% and 42% were observed for the linear, Gaussian, quadratic and cubic models respectively. Overall, the linear model performed best at predicting those patients that will benefit from RT adaptation. The linear model identified 21% and 23% of patients in our cohort with more than 30% tumour volume reduction to benefit from treatment adaptation in weeks 3 and 4 respectively.Significance. We have shown that it is feasible to predict the shape and volume of NSCLC tumours from routine CBCTs and effectively identify patients who will respond to treatment early.

Quantitative analysis of shear wave elastic heterogeneity for prediction of lymphovascular invasion in breast cancer

OBJECTIVE

To evaluate the correlation between elastic heterogeneity (EH) and lymphovascular invasion (LVI) in breast cancers and assess the clinical value of using EH to predict LVI pre-operatively.

METHODS

This retrospective study consisted of 376 patients with breast cancers that had undergone shear wave elastography (SWE) with virtual touch tissue imaging quantification between June 2017 and June 2018. The EH was determined as the difference between the averaged three highest and three lowest shear wave value. Clinicalpathological parameters including histological type and grades, LVI, axillary lymph node status and molecular markers (estrogen receptor, progesterone receptor, human epidermal growth factor receptor 2 and Ki-67) were reviewed and recorded. Relationship EH and clinicalpathological parameters was investigated respectively. The diagnostic performance of EH in distinguishing LVI or not was analyzed.

RESULTS

At multivariate regression analysis, only EH (p = 0.017) was positively correlated with LVI in all tumors. EH (p = 0.003) and Ki-67 (p = 0.025) were positively correlated with LVI in tumors ≤ 2 cm. None of clinicalpathological parameters were correlated with LVI in tumors > 2 cm (p > 0.05 for all). Using EH to predict LVI in tumors ≤ 2 cm, the sensitivity and negative predictive value were 93 and 89% respectively.

CONCLUSION

EH has the potential to be served as an imaging biomarker to predict LVI in breast cancer especially for tumors ≤ 2 cm.

ADVANCES IN KNOWLEDGE

There was no association between LVI and other most commonly used elastic features such as SWV_mean and SWV_max. Elastic heterogeneity is an independent predictor of LVI, so it can provide additional prognostic information for routine preoperative breast cancer assessment.For tumors ≤ 2cm, using EH value higher than 1.36 m/s to predict LVI involvement, the sensitivity and negative predictive value can reach to 93% and 89%, respectively, suggesting that breast cancer with negative EH value was more likely to be absent of LVI.

MRI-based delta-radiomics are predictive of pathological complete response after neoadjuvant chemoradiotherapy in locally advanced rectal cancer

RATIONALE AND OBJECTIVES

To investigate the capability of delta-radiomics to predict pathological complete response (pCR) after neoadjuvant chemoradiotherapy (nCRT) in patients with locally advanced rectal cancer (LARC).

MATERIALS AND METHODS

This retrospective study enrolled 165 consecutive patients with LARC (training set, n = 116; test set, n = 49) who received nCRT before surgery. All patients underwent pre- and post-nCRT MRI examination from which radiomics features were extracted. A delta-radiomics feature was defined as the percentage change in a radiomics feature from pre- to post-nCRT MRI. A data reduction and feature selection process including the least absolute shrinkage and selection operator algorithm was performed for building T2-weighted imaging (T2WI) and diffusion-weighted imaging (DWI) delta-radiomics signature. Logistic regression was used to build a T2WI and DWI combined radiomics model. Receiver operating characteristic analysis was performed to assess diagnostic performance. Delong method was used to compare the performance of delta-radiomics model with that of magnetic resonance tumor regression grade (mrTRG).

RESULTS

Twenty-seven of 165 patients (16.4%) achieved pCR. T2WI and DWI delta-radiomics signature, and the combined model showed good predictive performance for pCR. The combined model achieved the highest areas under the receiver operating characteristic curves of 0.91 (95% confidence interval: 0.85-0.98) and 0.91 (95% confidence interval: 0.83-0.99) in the training and test sets, respectively (significantly greater than those for mrTRG; training set, p < 0.001; test set, p = 0.04).

CONCLUSION

MRI-based delta-radiomics can help predict pCR after nCRT in patients with LARC with better performance than mrTRG.

Use of machine learning to select texture features in investigating the effects of axial loading on T2-maps from magnetic resonance imaging of the lumbar discs

BACKGROUND

Recent advances in texture analysis and machine learning offer new opportunities to improve the application of imaging to intervertebral disc biomechanics. This study employed texture analysis and machine learning on MRIs to investigate the lumbar disc's response to loading.

METHODS

Thirty-five volunteers (30 (SD 11) yrs.) with and without chronic back pain spent 20 min lying in a relaxed unloaded supine position, followed by 20 min loaded in compression, and then 20 min with traction applied. T₂-weighted MR images were acquired during the last 5 min of each loading condition. Custom image analysis software was used to segment discs from adjacent tissues semi-automatically and segment each disc into the nucleus, anterior and posterior annulus automatically. A grey-level, co-occurrence matrix with one to four pixels offset in four directions (0°, 45°, 90° and 135°) was then constructed (320 feature/tissue). The Random Forest Algorithm was used to select the most promising classifiers. Linear mixed-effect models and Cohen's d compared loading conditions.

FINDINGS

All statistically significant differences (p < 0.001) were observed in the nucleus and posterior annulus in the 135° offset direction at the L4-5 level between lumbar compression and traction. Correlation (P_2-Offset, P_4-Offset) and information measure of correlation 1 (P_3-Offset, P_4-Offset) detected significant changes in the nucleus. Statistically significant changes were also observed for homogeneity (P_2-Offset, P_3-Offset), contrast (P_2-Offset), and difference variance (P_4-Offset) of the posterior annulus.

INTERPRETATION

MRI textural features may have the potential of identifying the disc's response to loading, particularly in the nucleus and posterior annulus, which appear most sensitive to loading.

LEVEL OF EVIDENCE

Diagnostic: individual cross-sectional studies with consistently applied reference standard and blinding.

Decoding intra-tumoral spatial heterogeneity on radiological images using the Hilbert curve

Background

Current intra-tumoral heterogeneous feature extraction in radiology is limited to the use of a single slice or the region of interest within a few context-associated slices, and the decoding of intra-tumoral spatial heterogeneity using whole tumor samples is rare. We aim to propose a mathematical model of space-filling curve-based spatial correspondence mapping to interpret intra-tumoral spatial locality and heterogeneity.

Methods

A Hilbert curve-based approach was employed to decode and visualize intra-tumoral spatial heterogeneity by expanding the tumor volume to a two-dimensional (2D) matrix in voxels while preserving the spatial locality of the neighboring voxels. The proposed method was validated using three-dimensional (3D) volumes constructed from lung nodules from the LIDC-IDRI dataset, regular axial plane images, and 3D blocks.

Results

Dimensionality reduction of the Hilbert volume with a single regular axial plane image showed a sparse and scattered pixel distribution on the corresponding 2D matrix. However, for 3D blocks and lung tumor inside the volume, the dimensionality reduction to the 2D matrix indicated regular and concentrated squares and rectangles. For classification into benign and malignant masses using lung nodules from the LIDC-IDRI dataset, the Inception-V4 indicated that the Hilbert matrix images improved accuracy (85.54% vs. 73.22%, p < 0.001) compared to the original CT images of the test dataset.

Conclusions

Our study indicates that Hilbert curve-based spatial correspondence mapping is promising for decoding intra-tumoral spatial heterogeneity of partial or whole tumor samples on radiological images. This spatial-locality-preserving approach for voxel expansion enables existing radiomics and convolution neural networks to filter structured and spatially correlated high-dimensional intra-tumoral heterogeneity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13244-021-01100-8.

Subregional radiomics analysis for the detection of the EGFR mutation on thoracic spinal metastases from lung cancer

The present study intended to use radiomic analysis of spinal metastasis subregions to detect epidermal growth factor receptor (EGFR) mutation. In total, 94 patients with thoracic spinal metastasis originated from primary lung adenocarcinoma (2017-2020) were studied. All patients underwent T1-weighted (T1W) and T2 fat-suppressed (T2FS) MRI scans. The spinal metastases (tumor region) were subdivided into phenotypically consistent subregions based on patient- and population-level clustering: Three subregions, S1, S2 and S3, and the total tumor region. Radiomics features were extracted from each subregion and from the whole tumor region as well. Least shrinkage and selection operator (LASSO) regression were used for feature selection and radiomics signature definition. Detection performance of S3 was better than all other regions using T1W (AUCs, S1 versus S2 versus S3 versus whole tumor, 0.720 versus 0.764 versus 0.786 versus 0.758) and T2FS (AUCs, S1 versus S2 versus S3 versus whole tumor, 0.791 versus 0.708 versus 0.838 versus 0.797) MRI. The multi-regional radiomics signature derived from the joint of inner subregion S3 from T1W and T2FS MRI achieved the best detection capabilities with AUCs of 0.879 (ACC = 0.774, SEN = 0.838, SPE = 0.840) and 0.777 (ACC = 0.688, SEN = 0.947, SPE = 0.615) in the training and test sets, respectively. Our study revealed that MRI-based radiomic analysis of spinal metastasis subregions has the potential to detect the EGFR mutation in patients with primary lung adenocarcinoma.

Habitat Imaging-Based 18F-FDG PET/CT Radiomics for the Preoperative Discrimination of Non-small Cell Lung Cancer and Benign Inflammatory Diseases

Purpose

To propose and evaluate habitat imaging-based ¹⁸F-fluorodeoxyglucose (¹⁸F-FDG) positron emission tomography/computed tomography (PET/CT) radiomics for preoperatively discriminating non-small cell lung cancer (NSCLC) and benign inflammatory diseases (BIDs).

Methods

Three hundred seventeen ¹⁸F-FDG PET/CT scans were acquired from patients who underwent aspiration biopsy or surgical resection. All volumes of interest (VOIs) were semiautomatically segmented. Each VOI was separated into variant subregions, namely, habitat imaging, based on our adapted clustering-based habitat generation method. Radiomics features were extracted from these subregions. Three feature selection methods and six classifiers were applied to construct the habitat imaging-based radiomics models for fivefold cross-validation. The radiomics models whose features extracted by conventional habitat-based methods and nonhabitat method were also constructed. For comparison, the performances were evaluated in the validation set in terms of the area under the receiver operating characteristic curve (AUC). Pairwise t-test was applied to test the significant improvement between the adapted habitat-based method and the conventional methods.

Results

A total of 1,858 radiomics features were extracted. After feature selection, habitat imaging-based ¹⁸F-FDG PET/CT radiomics models were constructed. The AUC of the adapted clustering-based habitat radiomics was 0.7270 ± 0.0147, which showed significantly improved discrimination performance compared to the conventional methods (p <.001). Furthermore, the combination of features extracted by our adaptive habitat imaging-based method and non-habitat method showed the best performance than the other combinations.

Conclusion

Habitat imaging-based ¹⁸F-FDG PET/CT radiomics shows potential as a biomarker for discriminating NSCLC and BIDs, which indicates that the microenvironmental variations in NSCLC and BID can be captured by PET/CT.

Multi-level multi-modality (PET and CT) fusion radiomics: prognostic modeling for non-small cell lung carcinoma

We developed multi-modality radiomic models by integrating information extracted from¹⁸F-FDG PET and CT images using feature- and image-level fusions, toward improved prognosis for non-small cell lung carcinoma (NSCLC) patients. Two independent cohorts of NSCLC patients from two institutions (87 and 95 patients) were cycled as training and testing datasets. Fusion approaches were applied at two levels, namely feature- and image-levels. For feature-level fusion, radiomic features were extracted individually from CT and PET images and concatenated. Alternatively, radiomic features extracted separately from CT and PET images were averaged. For image-level fusion, wavelet fusion was utilized and tuned with two parameters, namely CT weight and Wavelet Band Pass Filtering Ratio. Clinical and combined clinical + radiomic models were developed. Gray level discretization was performed at 3 different levels (16, 32 and 64) and 225 radiomics features were extracted. Overall survival (OS) was considered as the endpoint. For feature reduction, correlated (redundant) features were excluded using Spearman's correlation, and best combination of top ten features with highest concordance-indices (via univariate Cox model) were selected in each model for further multivariate Cox model. Moreover, prognostic score's median, obtained from the training cohort, was used intact in the testing cohort as a threshold to classify patients into low- versus high-risk groups, and log-rank test was applied to assess differences between the Kaplan-Meier curves. Overall, while models based on feature-level fusion strategy showed limited superiority over single-modalities, image-level fusion strategy significantly outperformed both single-modality and feature-level fusion strategies. As such, the clinical model (C-index = 0.656) outperformed all models from single-modality and feature-level strategies, but was outperformed by certain models from image-level fusion strategy. Our findings indicated that image-level fusion multi-modality radiomics models outperformed single-modality, feature-level fusion, and clinical models for OS prediction of NSCLC patients.

Multiparametric magnetic resonance imaging-derived radiomics for the prediction of disease-free survival in early-stage squamous cervical cancer

OBJECTIVE

To conduct multiparametric magnetic resonance imaging (MRI)-derived radiomics based on multi-scale tumor region for predicting disease-free survival (DFS) in early-stage squamous cervical cancer (ESSCC).

METHODS

A total of 191 ESSCC patients (training cohort, n = 135; validation cohort, n = 56) from March 2016 to September 2019 were retrospectively recruited. Radiomics features were derived from the T2-weighted imaging (T2WI), contrast-enhanced T1-weighted imaging (CET1WI), diffusion-weighted imaging (DWI), and apparent diffusion coefficient (ADC) map for each patient. DFS-related radiomics features were selected in 3 target tumor volumes (VOI^entire, VOI^+5 mm, and VOI^-5 mm) to build 3 rad-scores using the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Logistic regression was applied to build combined model incorporating rad-scores with clinical risk factors and compared with clinical model alone. Kaplan-Meier analysis was used to further validate prognostic value of selected clinical and radiomics characteristics.

RESULTS

Three radiomics scores all showed favorable performances in DFS prediction. Rad-score (VOI^+5 mm) performed best with a C-index of 0.750 in the training set and 0.839 in the validation set. Combined model was constructed by incorporating age categorized by 55, Federation of Gynecology and Obstetrics (Figo) stage, and lymphovascular space invasion with rad-score (VOI^+5 mm). Combined model performed better than clinical model in DFS prediction in both the training set (C-index 0.815 vs 0.709; p = 0.024) and the validation set (C-index 0.866 vs 0.719; p = 0.001).

CONCLUSION

Multiparametric MRI-derived radiomics based on multi-scale tumor region can aid in the prediction of DFS for ESSCC patients, thereby facilitating clinical decision-making.

KEY POINTS

• Three radiomics scores based on multi-scale tumor region all showed favorable performances in DFS prediction. Rad-score (VOI^+5 mm) performed best with favorable C-index values. • Combined model incorporating multiparametric MRI-based radiomics with clinical risk factors performed significantly better in DFS prediction than the clinical model. • Combined model presented as a nomogram can be easily used to predict survival, thereby facilitating clinical decision-making.

Assessment of Intratumor Heterogeneity in Parametric Dynamic Contrast-Enhanced MR Images: A Comparative Study of Novel and Established Methods

Intratumor heterogeneity is associated with aggressive disease and poor survival rates in several types of cancer. A novel method for assessing intratumor heterogeneity in medical images, named the spatial gradient method, has been developed in our laboratory. In this study, we measure intratumor heterogeneity in K^trans maps derived by dynamic contrast-enhanced magnetic resonance imaging using the spatial gradient method, and we compare the performance of the novel method with that of histogram analyses and texture analyses using the Haralick method. K^trans maps of 58 untreated and sunitinib-treated pancreatic ductal adenocaricoma (PDAC) xenografts from two PDAC models were investigated. Intratumor heterogeneity parameters derived by the spatial gradient method were sensitive to tumor line differences as well as sunitinib-induced changes in intratumor heterogeneity. Furthermore, the parameters provided additional information to the median value and were not severely affected by imaging noise. The parameters derived by histogram analyses were insensitive to spatial heterogeneity and were strongly correlated to the median value, and the Haralick features were severely influenced by imaging noise and did not differentiate between untreated and sunitinib-treated tumors. The spatial gradient method was superior to histogram analyses and Haralick features for assessing intratumor heterogeneity in K^trans maps of untreated and sunitinib-treated PDAC xenografts, and can possibly be used to assess intratumor heterogeneity in other medical images and to evaluate effects of other treatments as well.

Predicting Overall Survival Time in Glioblastoma Patients Using Gradient Boosting Machines Algorithm and Recursive Feature Elimination Technique

Despite advances in tumor treatment, the inconsistent response is a major challenge among glioblastoma multiform (GBM) that lead to different survival time. Our aim was to integrate multimodal MRI with non-supervised and supervised machine learning methods to predict GBM patients' survival time. To this end, we identified different compartments of the tumor and extracted their features. Next, we applied Random Forest-Recursive Feature Elimination (RF-RFE) to identify the most relevant features to feed into a GBoost machine. This study included 29 GBM patients with known survival time. RF-RFE GBoost model was evaluated to assess the survival prediction performance using optimal features. Furthermore, overall survival (OS) was analyzed using univariate and multivariate Cox regression analyses, to evaluate the effect of ROIs and their features on survival. The results showed that a RF-RFE Gboost machine was able to predict survival time with 75% accuracy. The results also revealed that the rCBV in the low perfusion area was significantly different between groups and had the greatest effect size in terms of the rate of change of the response variable (survival time). In conclusion, not only integration of multi-modality MRI but also feature selection method can enhance the classifier performance.

Image intensity histograms as imaging biomarkers: application to immune-related colitis

Purpose.To investigate image intensity histograms as a potential source of useful imaging biomarkers in both a clinical example of detecting immune-related colitis (irColitis) in¹⁸F-FDG PET/CT images of immunotherapy patients and an idealized case of classifying digital reference objects (DRO).Methods.Retrospective analysis of bowel¹⁸F-FDG uptake in N = 40 patients receiving immune checkpoint inhibitors was conducted. A CNN trained to segment the bowel was used to generate the histogram of bowel¹⁸F-FDG uptake, and percentiles of the histogram were considered as potential metrics for detecting inflammation associated with irColitis. A model of the colon was also considered using cylindrical DRO. Classification of DRO with different intensity distributions was undertaken under varying geometry and noise settings.Results.The most predictive biomarker of irColitis was the 95th percentile of the bowel SUV histogram (SUV_95%). Patients later diagnosed with irColitis had a significantly higher increase in SUV_95%from baseline to first on-treatment PET than patients who did not experience irColitis (p = 0.02). An increase in SUV_95%> + 40% separated pre-irColitis change from normal variability with a sensitivity of 75% and specificity of 88%. Furthermore, histogram percentiles were ideal metrics for classifying 'hot center' and 'cold center' DRO, and were robust to varying DRO geometry and noise, and to the presence of spoiler volumes unrelated to the detection task.Conclusions.The 95th percentile of the bowel SUV histogram was the optimal metric for detecting irColitis on¹⁸F-FDG PET/CT. Image intensity histograms are a promising source of imaging biomarkers for clinical tasks.

Radiomics as a New Frontier of Imaging for Cancer Prognosis: A Narrative Review

The evaluation of the efficacy of different therapies is of paramount importance for the patients and the clinicians in oncology, and it is usually possible by performing imaging investigations that are interpreted, taking in consideration different response evaluation criteria. In the last decade, texture analysis (TA) has been developed in order to help the radiologist to quantify and identify parameters related to tumor heterogeneity, which cannot be appreciated by the naked eye, that can be correlated with different endpoints, including cancer prognosis. The aim of this work is to analyze the impact of texture in the prediction of response and in prognosis stratification in oncology, taking into consideration different pathologies (lung cancer, breast cancer, gastric cancer, hepatic cancer, rectal cancer). Key references were derived from a PubMed query. Hand searching and clinicaltrials.gov were also used. This paper contains a narrative report and a critical discussion of radiomics approaches related to cancer prognosis in different fields of diseases.

Development and Validation of a 18F-FDG PET-Based Radiomic Model for Evaluating Hypermetabolic Mediastinal-Hilar Lymph Nodes in Non-Small-Cell Lung Cancer

Background

Accurate evaluation of lymph node (LN) status is critical for determining the treatment options in patients with non-small cell lung cancer (NSCLC). This study aimed to develop and validate a ¹⁸F-FDG PET-based radiomic model for the identification of metastatic LNs from the hypermetabolic mediastinal–hilar LNs in NSCLC.

Methods

We retrospectively reviewed 259 patients with hypermetabolic LNs who underwent pretreatment ¹⁸F-FDG PET/CT and were pathologically confirmed as NSCLC from two centers. Two hundred twenty-eight LNs were allocated to a training cohort (LN = 159) and an internal validation cohort (LN = 69) from one center (7:3 ratio), and 60 LNs were enrolled to an external validation cohort from the other. Radiomic features were extracted from LNs of PET images. A PET radiomics signature was constructed by multivariable logistic regression after using the least absolute shrinkage and selection operator (LASSO) method with 10-fold cross-validation. The PET radiomics signature (model 1) and independent predictors from CT image features and clinical data (model 2) were incorporated into a combined model (model 3). A nomogram was plotted for the complex model, and the performance of the nomogram was assessed by its discrimination, calibration, and clinical usefulness.

Results

The area under the curve (AUC) values of model 1 were 0.820, 0.785, and 0.808 in the training, internal, and external validation cohorts, respectively, showing good diagnostic efficacy for lymph node metastasis (LNM). Furthermore, model 2 was able to discriminate metastatic LNs in the training (AUC 0.780), internal (AUC 0.794), and external validation cohorts (AUC 0.802), respectively. Model 3 showed optimal diagnostic performance among the three cohorts, with an AUC of 0.874, 0.845, and 0.841, respectively. The nomogram based on the model 3 showed good discrimination and calibration.

Conclusions

Our study revealed that PET radiomics signature, especially when integrated with CT imaging features, showed the ability to identify true and false positives of mediastinal–hilar LNM detected by PET/CT in patients with NSCLC, which would help clinicians to make individual treatment decisions.

Non-Invasive Ultrasonic Description of Tumor Evolution

Simple Summary

During tumor evolution, heterogeneous structural and functional changes occur in the tumor microenvironment. These complex changes have pro- or anti-tumorigenesis effects and have an impact on therapy efficiency. Therefore, the tumor microenvironment needs to be non-invasively characterized over time. The aim of this preclinical work is to compare the sensitivity of modifications occurring during tumor evolution of volume, immunohistochemistry and non-invasive quantitative ultrasound parameters (Shear Wave Elastography and dynamic Contrast-Enhanced Ultrasound) and to study the link between them. The complementary evaluation over time of multiple morphological and functional parameters during tumor growth underlines the need to integrate histological, morphological, functional, and, ultimately, genomic information into models that can consider the temporal and spatial variability of features to better understand tumor evolution.

Abstract

Purpose: There is a clinical need to better non-invasively characterize the tumor microenvironment in order to reveal evidence of early tumor response to therapy and to better understand therapeutic response. The goals of this work are first to compare the sensitivity to modifications occurring during tumor growth for measurements of tumor volume, immunohistochemistry parameters, and emerging ultrasound parameters (Shear Wave Elastography (SWE) and dynamic Contrast-Enhanced Ultrasound (CEUS)), and secondly, to study the link between the different parameters. Methods: Five different groups of 9 to 10 BALB/c female mice with subcutaneous CT26 tumors were imaged using B-mode morphological imaging, SWE, and CEUS at different dates. Whole-slice immunohistological data stained for the nuclei, T lymphocytes, apoptosis, and vascular endothelium from these tumors were analyzed. Results: Tumor volume and three CEUS parameters (Time to Peak, Wash-In Rate, and Wash-Out Rate) significantly changed over time. The immunohistological parameters, CEUS parameters, and SWE parameters showed intracorrelation. Four immunohistological parameters (the number of T lymphocytes per mm² and its standard deviation, the percentage area of apoptosis, and the colocalization of apoptosis and vascular endothelium) were correlated with the CEUS parameters (Time to Peak, Wash-In Rate, Wash-Out Rate, and Mean Transit Time). The SWE parameters were not correlated with the CEUS parameters nor with the immunohistological parameters. Conclusions: US imaging can provide additional information on tumoral changes. This could help to better explore the effect of therapies on tumor evolution, by studying the evolution of the parameters over time and by studying their correlations.

The exploration of quantitative intra-tumoral metabolic heterogeneity in dual-time 18F-FDG PET/CT of pancreatic cancer

PURPOSE

We aimed to analyze the change of quantitative intra-tumoral metabolic heterogeneity consisting of texture features and conventional metabolic parameters of pancreatic cancer (PC) in dual-time 2-deoxy-2(¹⁸F) fluoro-D-glucose (¹⁸F-FDG) positron emission tomography-computed tomography (PET/CT).

METHODS

A retrospective analysis was conducted considering the texture features and conventional metabolic parameters in dual-time ¹⁸F-FDG PET/CT scans of PC patients. Features were extracted based on spatial distribution of ¹⁸F-FDG uptake in image. Firstly, the texture features and the conventional metabolic parameters of the delayed scan were both compared with that of the early scan. Statistically different data was defined among them. Secondly, the study evaluated the correlations between retention index (RI) of the texture features and the conventional metabolic parameters. Finally, the variation of texture features in dual-time PET/CT of resectable PC patients and unresectable PC patients was calculated separately.

RESULTS

In total, 183 PC patients were analyzed retrospectively in this research. The conventional metabolic parameters were all statistically different between the early and delayed scans except for metabolic tumor volume (MTV). In the radiomics, there were 59 textural features. Nineteen of 59 texture features were statistically different between the early and delayed scans. Features that were more than 10% different during two scans were observed in a substantial percentage of patients. Weak correlations were only found between MTV, TLG (Total lesion glycolysis), SUVpeak and the RI of some texture features in early or delayed scans. There were obviously fewer features with significant difference in resectable PC group than in unresectable PC group. Most features showing the difference in unresectable group while no significant difference in resectable group.

CONCLUSIONS

This study investigated the change and inner correlations of quantitative tumoral metabolic heterogeneity in the dual-time ¹⁸F-FDG-PET/CT scan of PC patients. Some features displayed the difference between dual-time scans. Conventional metabolic parameters were weakly related to the change of texture feature. The change of texture feature in resectable PC group was different from that in unresectable PC group. This result is potential to provide more information for the image evaluation of PC.

Radiomics Study for Predicting the Expression of PD-L1 and Tumor Mutation Burden in Non-Small Cell Lung Cancer Based on CT Images and Clinicopathological Features

Background

The present study compared the predictive performance of pretreatment computed tomography (CT)-based radiomics signatures and clinicopathological and CT morphological factors for ligand programmed death-ligand 1 (PD-L1) expression level and tumor mutation burden (TMB) status and further explored predictive models in patients with advanced-stage non-small cell lung cancer (NSCLC).

Methods

A total of 120 patients with advanced-stage NSCLC were enrolled in this retrospective study and randomly assigned to a training dataset or validation dataset. Here, 462 radiomics features were extracted from region-of-interest (ROI) segmentation based on pretreatment CT images. The least absolute shrinkage and selection operator (LASSO) and logistic regression were applied to select radiomics features and develop combined models with clinical and morphological factors for PD-L1 expression and TMB status prediction. Ten-fold cross-validation was used to evaluate the accuracy, and the predictive performance of these models was assessed using receiver operating characteristic (ROC) and area under the curve (AUC) analyses.

Results

The PD-L1-positive expression level correlated with differentiation degree (p = 0.005), tumor shape (p = 0.006), and vascular convergence (p = 0.007). Stage (p = 0.023), differentiation degree (p = 0.017), and vacuole sign (p = 0.016) were associated with TMB status. Radiomics signatures showed good performance for predicting PD-L1 and TMB with AUCs of 0.730 and 0.759, respectively. Predictive models that combined radiomics signatures with clinical and morphological factors dramatically improved the predictive efficacy for PD-L1 (AUC = 0.839) and TMB (p = 0.818). The results were verified in the validation datasets.

Conclusions

Quantitative CT-based radiomics features have potential value in the classification of PD-L1 expression levels and TMB status. The combined model further improved the predictive performance and provided sufficient information for the guiding of immunotherapy in clinical practice, and it deserves further analysis.

Subregional Radiomics Analysis of PET/CT Imaging with Intratumor Partitioning: Application to Prognosis for Nasopharyngeal Carcinoma

PURPOSE

This work aims to identify intratumoral habitats with distinct heterogeneity based on 2-deoxy-2-[¹⁸F]fluro-D-glucose positron emission tomography (PET)/X-ray computed tomography (CT) imaging, and to develop a subregional radiomics approach to predict progression-free survival (PFS) in patients with nasopharyngeal carcinoma (NPC).

PROCEDURES

In total, 128 NPC patients (85 vs. 43 for primary vs. validation cohorts) who underwent pre-treatment PET/CT scan were enrolled retrospectively. Each tumor was partitioned into several phenotypically consistent subregions based on individual- and population-level clustering. For each subregion, 202 radiomics features were extracted to construct imaging biomarker for prognosis via Cox's proportional hazard model combined with forward stepwise feature selection. Relevance of imaging biomarkers and clinicopathological factors were assessed by multivariate Cox regression analysis and Spearman's correlation analysis. To investigate whether imaging biomarkers could provide complementary prognosis information beyond existing predictors, a scoring system was further developed for risk stratification and compared with AJCC staging system.

RESULTS

Three subregions (denoted as S1, S2, and S3) were discovered with distinct PET/CT imaging characteristics in the two cohorts. The prognostic performance of imaging biomarker S3 outperformed the whole tumor (C-index, 0.69 vs. 0.58; log-rank test, p < 0.001 vs. p = 0.552). Imaging biomarker S3 and AJCC stage were identified as independent predictors (p = 0.011 and 0.042, respectively) after adjusting for clinicopathological factors. The scoring system outperformed the traditional AJCC staging system (log-rank test, p < 0.0001 vs. p = 0.0002 in primary cohort and p = 0.0021 vs. p = 0.0277 in validation cohort, respectively).

CONCLUSIONS

Subregional radiomics analysis of PET/CT imaging has the potential to predict PFS in patients with NPC, which also provides complementary prognostic information for traditional predictors.

Tumor habitat analysis by magnetic resonance imaging distinguishes tumor progression from radiation necrosis in brain metastases after stereotactic radiosurgery

OBJECTIVES

The identification of viable tumor after stereotactic radiosurgery (SRS) is important for future targeted therapy. This study aimed to determine whether tumor habitat on structural and physiologic MRI can distinguish viable tumor from radiation necrosis of brain metastases after SRS.

METHOD

Multiparametric contrast-enhanced T1- and T2-weighted imaging, apparent diffusion coefficient (ADC), and cerebral blood volume (CBV) were obtained from 52 patients with 69 metastases, showing enlarging enhancing masses after SRS. Voxel-wise clustering identified three structural MRI habitats (enhancing, solid low-enhancing, and nonviable) and three physiologic MRI habitats (hypervascular cellular, hypovascular cellular, and nonviable). Habitat-based predictors for viable tumor or radiation necrosis were identified by logistic regression. Performance was validated using the area under the curve (AUC) of the receiver operating characteristics curve in an independent dataset with 24 patients.

RESULTS

None of the physiologic MRI habitats was indicative of viable tumor. Viable tumor was predicted by a high-volume fraction of solid low-enhancing habitat (low T2-weighted and low CE-T1-weighted values; odds ratio [OR] 1.74, p <.001) and a low-volume fraction of nonviable tissue habitat (high T2-weighted and low CE-T1-weighted values; OR 0.55, p <.001). Combined structural MRI habitats yielded good discriminatory ability in both development (AUC 0.85, 95% confidence interval [CI]: 0.77-0.94) and validation sets (AUC 0.86, 95% CI:0.70-0.99), outperforming single ADC (AUC 0.64) and CBV (AUC 0.58) values. The site of progression matched with the solid low-enhancing habitat (72%, 8/11).

CONCLUSION

Solid low-enhancing and nonviable tissue habitats on structural MRI can help to localize viable tumor in patients with brain metastases after SRS.

KEY POINTS

• Structural MRI habitats helped to differentiate viable tumor from radiation necrosis. • Solid low-enhancing habitat was most helpful to find viable tumor. • Providing spatial information, the site of progression matched with solid low-enhancing habitat.

Multiblock Discriminant Analysis of Integrative 18F-FDG-PET/CT Radiomics for Predicting Circulating Tumor Cells in Early-Stage Non-small Cell Lung Cancer Treated With Stereotactic Body Radiation Therapy

PURPOSE

The main objective of the present study was to integrate 18F-FDG-PET/CT radiomics with multiblock discriminant analysis for predicting circulating tumor cells (CTCs) in early-stage non-small cell lung cancer (ES-NSCLC) treated with stereotactic body radiation therapy (SBRT).

METHODS

Fifty-six patients with stage I NSCLC treated with SBRT underwent 18F-FDG-PET/CT imaging pre-SBRT and post-SBRT (median, 5 months; range, 3-10 months). CTCs were assessed via a telomerase-based assay before and within 3 months after SBRT and dichotomized at 5 and 1.3 CTCs/mL. Pre-SBRT, post-SBRT, and delta PET/CT radiomics features (n = 1548 × 3/1562 × 3) were extracted from gross tumor volume. Seven feature blocks were constructed including clinical parameters (n = 12). Multiblock data integration was performed using block sparse partial least squares-discriminant analysis (sPLS-DA) referred to as Data Integration Analysis for Biomarker Discovery Using Latent Components (DIABLO) for identifying key signatures by maximizing common information between different feature blocks while discriminating CTC levels. Optimal input blocks were identified using a pairwise combination method. DIABLO performance for predicting pre-SBRT and post-SBRT CTCs was evaluated using combined AUC (area under the curve, averaged across different blocks) analysis with 20 × 5-fold cross-validation (CV) and compared with that of concatenation-based sPLS-DA that consisted of combining all features into 1 block. CV prediction scores between 1 class versus the other were compared using the Wilcoxon rank sum test.

RESULTS

For predicting pre-SBRT CTCs, DIABLO achieved the best performance with combined pre-SBRT PET radiomics and clinical feature blocks, showing CV AUC of 0.875 (P = .009). For predicting post-SBRT CTCs, DIABLO achieved the best performance with combined post-SBRT CT and delta CT radiomics feature blocks, showing CV AUCs of 0.883 (P = .001). In contrast, all single-block sPLS-DA models could not attain CV AUCs higher than 0.7.

CONCLUSIONS

Multiblock integration with discriminant analysis of 18F-FDG-PET/CT radiomics has the potential for predicting pre-SBRT and post-SBRT CTCs. Radiomics and CTC analysis may complement and together help guide the subsequent management of patients with ES-NSCLC.

The Effect of Heterogenous Subregions in Glioblastomas on Survival Stratification: A Radiomics Analysis Using the Multimodality MRI

Intratumor heterogeneity is partly responsible for the poor prognosis of glioblastoma (GBM) patients. In this study, we aimed to assess the effect of different heterogeneous subregions of GBM on overall survival (OS) stratification. A total of 105 GBM patients were retrospectively enrolled and divided into long-term and short-term OS groups. Four MRI sequences, including contrast-enhanced T1-weighted imaging (T1C), T1, T2, and FLAIR, were collected for each patient. Then, 4 heterogeneous subregions, i.e. the region of entire abnormality (rEA), the regions of contrast-enhanced tumor (rCET), necrosis (rNec) and edema/non-contrast-enhanced tumor (rE/nCET), were manually drawn from the 4 MRI sequences. For each subregion, 50 radiomics features were extracted. The stratification performance of 4 heterogeneous subregions, as well as the performances of 4 MRI sequences, was evaluated both alone and in combination. Our results showed that rEA was superior in stratifying long-and short-term OS. For the 4 MRI sequences used in this study, the FLAIR sequence demonstrated the best performance of survival stratification based on the manual delineation of heterogeneous subregions. Our results suggest that heterogeneous subregions of GBMs contain different prognostic information, which should be considered when investigating survival stratification in patients with GBM.

Characterizing intra-tumor regions on quantitative ultrasound parametric images to predict breast cancer response to chemotherapy at pre-treatment

The efficacy of quantitative ultrasound (QUS) multi-parametric imaging in conjunction with unsupervised classification algorithms was investigated for the first time in characterizing intra-tumor regions to predict breast tumor response to chemotherapy before the start of treatment. QUS multi-parametric images of breast tumors were generated using the ultrasound radiofrequency data acquired from 181 patients diagnosed with locally advanced breast cancer and planned for neo-adjuvant chemotherapy followed by surgery. A hidden Markov random field (HMRF) expectation maximization (EM) algorithm was applied to identify distinct intra-tumor regions on QUS multi-parametric images. Several features were extracted from the segmented intra-tumor regions and tumor margin on different parametric images. A multi-step feature selection procedure was applied to construct a QUS biomarker consisting of four features for response prediction. Evaluation results on an independent test set indicated that the developed biomarker coupled with a decision tree model with adaptive boosting (AdaBoost) as the classifier could predict the treatment response of patient at pre-treatment with an accuracy of 85.4% and an area under the receiver operating characteristic (ROC) curve (AUC) of 0.89. In comparison, the biomarkers consisted of the features derived from the entire tumor core (without consideration of the intra-tumor regions), and the entire tumor core and the tumor margin could predict the treatment response of patients with an accuracy of 74.5% and 76.4%, and an AUC of 0.79 and 0.76, respectively. Standard clinical features could predict the therapy response with an accuracy of 69.1% and an AUC of 0.6. Long-term survival analyses indicated that the patients predicted by the developed model as responders had a significantly better survival compared to the non-responders. Similar findings were observed for the two response cohorts identified at post-treatment based on standard clinical and pathological criteria. The results obtained in this study demonstrated the potential of QUS multi-parametric imaging integrated with unsupervised learning methods in identifying distinct intra-tumor regions in breast cancer to characterize its responsiveness to chemotherapy prior to the start of treatment.

Target Heterogeneity in Oncology: The Best Predictor for Differential Response to Radioligand Therapy in Neuroendocrine Tumors and Prostate Cancer

Simple Summary

In the era of precision medicine, novel targets have emerged on the surface of cancer cells, which have been exploited for the purpose of radioligand therapy. However, there have been variations in the way these receptors are expressed, especially in prostate cancers and neuroendocrine tumors. This variable expression of receptors across the grades of cancers led to the concept of ‘target heterogeneity’, which has not just impacted therapeutic decisions but also their outcomes. Radiopharmaceuticals targeting receptors need to be used when there are specific indicators—either clinical, radiological, or at molecular level—warranting their use. In addition, response to these radioligands can be assessed using different techniques, whereby we can prognosticate further outcomes. We shall also discuss, in this review, the conventional as well as novel approaches of detecting heterogeneity in prostate cancers and neuroendocrine tumors.

Abstract

Tumor or target heterogeneity (TH) implies presence of variable cellular populations having different genomic characteristics within the same tumor, or in different tumor sites of the same patient. The challenge is to identify this heterogeneity, as it has emerged as the most common cause of ‘treatment resistance’, to current therapeutic agents. We have focused our discussion on ‘Prostate Cancer’ and ‘Neuroendocrine Tumors’, and looked at the established methods for demonstrating heterogeneity, each with its advantages and drawbacks. Also, the available theranostic radiotracers targeting PSMA and somatostatin receptors combined with targeted systemic agents, have been described. Lu-177 labeled PSMA and DOTATATE are the ‘standard of care’ radionuclide therapeutic tracers for management of progressive treatment-resistant prostate cancer and NET. These approved therapies have shown reasonable benefit in treatment outcome, with improvement in quality of life parameters. Various biomarkers and predictors of response to radionuclide therapies targeting TH which are currently available and those which can be explored have been elaborated in details. Imaging-based features using artificial intelligence (AI) need to be developed to further predict the presence of TH. Also, novel theranostic tools binding to newer targets on surface of cancer cell should be explored to overcome the treatment resistance to current treatment regimens.

Predictive Radiomic Models for the Chemotherapy Response in Non-Small-Cell Lung Cancer based on Computerized-Tomography Images

The heterogeneity and complexity of non-small cell lung cancer (NSCLC) tumors mean that NSCLC patients at the same stage can have different chemotherapy prognoses. Accurate predictive models could recognize NSCLC patients likely to respond to chemotherapy so that they can be given personalized and effective treatment. We propose to identify predictive imaging biomarkers from pre-treatment CT images and construct a radiomic model that can predict the chemotherapy response in NSCLC. This single-center cohort study included 280 NSCLC patients who received first-line chemotherapy treatment. Non-contrast CT images were taken before and after the chemotherapy, and clinical information were collected. Based on the Response Evaluation Criteria in Solid Tumors and clinical criteria, the responses were classified into two categories: response (n = 145) and progression (n = 135), then all data were divided into two cohorts: training cohort (224 patients) and independent test cohort (56 patients). In total, 1629 features characterizing the tumor phenotype were extracted from a cube containing the tumor lesion cropped from the pre-chemotherapy CT images. After dimensionality reduction, predictive models of the chemotherapy response of NSCLC with different feature selection methods and different machine-learning classifiers (support vector machine, random forest, and logistic regression) were constructed. For the independent test cohort, the predictive model based on a random-forest classifier with 20 radiomic features achieved the best performance, with an accuracy of 85.7% and an area under the receiver operating characteristic curve of 0.941 (95% confidence interval, 0.898–0.982). Of the 20 selected features, four were first-order statistics of image intensity and the others were texture features. For nine features, there were significant differences between the response and progression groups (p < 0.001). In the response group, three features, indicating heterogeneity, were overrepresented and one feature indicating homogeneity was underrepresented. The proposed radiomic model with pre-chemotherapy CT features can predict the chemotherapy response of patients with non-small cell lung cancer. This radiomic model can help to stratify patients with NSCLC, thereby offering the prospect of better treatment.

Artificial intelligence in tumor subregion analysis based on medical imaging: A review

Medical imaging is widely used in the diagnosis and treatment of cancer, and artificial intelligence (AI) has achieved tremendous success in medical image analysis. This paper reviews AI-based tumor subregion analysis in medical imaging. We summarize the latest AI-based methods for tumor subregion analysis and their applications. Specifically, we categorize the AI-based methods by training strategy: supervised and unsupervised. A detailed review of each category is presented, highlighting important contributions and achievements. Specific challenges and potential applications of AI in tumor subregion analysis are discussed.

A subregion-based positron emission tomography/computed tomography (PET/CT) radiomics model for the classification of non-small cell lung cancer histopathological subtypes

Background

This study classifies lung adenocarcinoma (ADC) and squamous cell carcinoma (SCC) using subregion-based radiomics features extracted from positron emission tomography/computed tomography (PET/CT) images.

Methods

In this study, the standard ¹⁸F-fluorodeoxyglucose (FDG) PET/CT images of 150 patients with lung ADC and 100 patients with SCC were retrospectively collected from the PET Center of the First Affiliated Hospital, College of Medicine, Zhejiang University. First, the 3D feature vector of each tumor voxel (whose basis is PET value, CT value, and CT local dominant orientation) was extracted. Using K-means individual clustering and population clustering, each tumor was divided into 4 subregions that reflect intratumoral regional heterogeneity. Next, based on each subregion, 385 radiomics features were extracted. Clinical features including age, gender, and smoking history were included. Thus, there were a total of 1,543 features extracted from PET/CT images and clinical reports. Statistical tests were then used to eliminate irrelevant and redundant features, and the recursive feature elimination (RFE) algorithm was used to select the best feature subset to classify SCC and ADC. Finally, 7 types of classifiers were tested to achieve the optimized model for the classification: support vector machine (SVM) with linear kernel, SVM with radial basis function kernel (SVM-RBF), random forest, logistic regression, Gaussian process classifier, linear discriminant analysis, and the AdaBoost classifier. Furthermore, 5-fold cross-validation was applied to obtain the sensitivity, specificity, accuracy, and area under the curve (AUC) for performance evaluation.

Results

Our model exhibited the best performance with the subregion radiomics features and SVM-RBF classifier, with a 5-fold cross-validation sensitivity, specificity, accuracy, and AUC of 0.8538, 0.8758, 0.8623, and 0.9155, respectively. The interquartile range feature from subregion 2 of CT and the gender feature from the clinical reports are the 2 optimized features that achieved the highest comprehensive score.

Conclusions

Our proposed model showed that SCC and ADC could be classified successfully using PET/CT images, which could be a promising tool to assist radiologists or medical physicists during diagnosis. The subregion-based method illustrated that non-small cell lung cancer (NSCLC) depicts intratumoral regional heterogeneity on both CT and PET images. By defining these heterogeneities through a subregion-based method, the diagnostic performance was improved. The 3D feature vector (whose basis is PET value, CT value, and CT local dominant orientation) showed superiority in reflecting NSCLC intratumoral regional heterogeneity.

Application of radiomics in hepatocellular carcinoma: A review

Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer with low 5-year survival rate. The high molecular heterogeneity in HCC poses huge challenges for clinical practice or trial design and has become a major barrier to improving the management of HCC. However, current clinical practice based on single bioptic or archived tumor tissue has been deficient in identifying useful biomarkers. The concept of radiomics was first proposed in 2012 and is different from the traditional imaging analysis based on the qualitative or semi-quantitative analysis by radiologists. Radiomics refers to high-throughput extraction of large amounts number of high-dimensional quantitative features from medical images through machine learning or deep learning algorithms. Using the radiomics method could quantify tumoral phenotypes and heterogeneity, which may provide benefits in clinical decision-making at a lower cost. Here, we review the workflow and application of radiomics in HCC.

Mask-Guided Convolutional Neural Network for Breast Tumor Prognostic Outcome Prediction on 3D DCE-MR Images

In this proof-of-concept work, we have developed a 3D-CNN architecture that is guided by the tumor mask for classifying several patient-outcomes in breast cancer from the respective 3D dynamic contrast-enhanced MRI (DCE-MRI) images. The tumor masks on DCE-MRI images were generated using pre- and post-contrast images and validated by experienced radiologists. We show that our proposed mask-guided classification has a higher accuracy than that from either the full image without tumor masks (including background) or the masked voxels only. We have used two patient outcomes for this study: (1) recurrence of cancer after 5 years of imaging and (2) HER2 status, for comparing accuracies of different models. By looking at the activation maps, we conclude that an image-based prediction model using 3D-CNN could be improved by even a conservatively generated mask, rather than overly trusting an unguided, blind 3D-CNN. A blind CNN may classify accurately enough, while its attention may really be focused on a remote region within 3D images. On the other hand, only using a conservatively segmented region may not be as good for classification as using full images but forcing the model's attention toward the known regions of interest.

Protocol for Creating Antibodies with Complex Fluorescence Spectra

Fluorescent antibodies are a workhorse of biomedical science, but fluorescence multiplexing has been notoriously difficult due to spectral overlap between fluorophores. We recently established proof-of-principal for fluorescence Multiplexing using Spectral Imaging and Combinatorics (MuSIC), which uses combinations of existing fluorophores to create unique spectral signatures for increased multiplexing. However, a method for labeling antibodies with MuSIC probes has not yet been developed. Here, we present a method for labeling antibodies with MuSIC probes. We conjugate a DBCO-Peg5-NHS ester linker to antibodies and a single-stranded DNA "docking strand" to the linker and, finally, hybridize two MuSIC-compatible, fluorescently labeled oligos to the docking strand. We validate the labeling protocol with spin-column purification and absorbance measurements. We demonstrate the approach using (i) Cy3, (ii) Tex615, and (iii) a Cy3-Tex615 combination as three different MuSIC probes attached to three separate batches of antibodies. We created single-, double-, and triple-positive beads that are analogous to single cells by incubating MuSIC probe-labeled antibodies with protein A beads. Spectral flow cytometry experiments demonstrate that each MuSIC probe can be uniquely distinguished, and the fraction of beads in a mixture with different staining patterns are accurately inferred. The approach is general and might be more broadly applied to cell-type profiling or tissue heterogeneity studies in clinical, biomedical, and drug discovery research.

Hypoxia Supports Differentiation of Terminally Exhausted CD8 T Cells

Hypoxia, angiogenesis, and immunosuppression have been proposed to be interrelated events that fuel tumor progression and impair the clinical effectiveness of anti-tumor therapies. Here we present new mechanistic data highlighting the role of hypoxia in fine-tuning CD8 T cell exhaustion in vitro, in an attempt to reconcile seemingly opposite evidence regarding the impact of hypoxia on functional features of exhausted CD8 T cells. Focusing on the recently characterized terminally-differentiated and progenitor exhausted CD8 T cells, we found that both hypoxia and its regulated mediator, vascular endothelial growth factor (VEGF)-A, promote the differentiation of PD-1+ TIM-3+ CXCR5+ terminally exhausted-like CD8 T cells at the expense of PD-1+ TIM-3- progenitor-like subsets without affecting tumor necrosis factor (TNF)-α and interferon (IFN)-γ production or granzyme B (GZMB) expression by these subpopulations. Interestingly, hypoxia accentuated the proangiogenic secretory profile in exhausted CD8 T cells. VEGF-A was the main factor differentially secreted by exhausted CD8 T cells under hypoxic conditions. In this sense, we found that VEGF-A contributes to generation of terminally exhausted CD8 T cells during in vitro differentiation. Altogether, our findings highlight the reciprocal regulation between hypoxia, angiogenesis, and immunosuppression, providing a rational basis to optimize synergistic combinations of antiangiogenic and immunotherapeutic strategies, with the overarching goal of improving the efficacy of these treatments.

The Role of Radiomics in Lung Cancer: From Screening to Treatment and Follow-Up

Purpose

Lung cancer represents the first cause of cancer-related death in the world. Radiomics studies arise rapidly in this late decade. The aim of this review is to identify important recent publications to be synthesized into a comprehensive review of the current status of radiomics in lung cancer at each step of the patients’ care.

Methods

A literature review was conducted using PubMed/Medline for search of relevant peer-reviewed publications from January 2012 to June 2020

Results

We identified several studies at each point of patient’s care: detection and classification of lung nodules (n=16), determination of histology and genomic (n=10) and finally treatment outcomes predictions (=23). We reported the methodology of those studies and their results and discuss the limitations and the progress to be made for clinical routine applications.

Conclusion

Promising perspectives arise from machine learning applications and radiomics based models in lung cancers, yet further data are necessary for their implementation in daily care. Multicentric collaboration and attention to quality and reproductivity of radiomics studies should be further consider.

Preoperative prediction of regional lymph node metastasis of colorectal cancer based on 18F-FDG PET/CT and machine learning

PURPOSE

To establish and validate a regional lymph node (LN) metastasis prediction model of colorectal cancer (CRC) based on 18F-FDG PET/CT and radiomic features using machine-learning methods.

METHODS

A total of 199 colorectal cancer patients underwent pre-therapy diagnostic 18F-FDG PET/CT scans and CRC radical surgery. The Chang-Gung Image Texture Analysis toolbox (CGITA) was used to extract 70 PET radiomic features reflecting 18F-FDG uptake heterogeneity of tumors. The least absolute shrinkage and selection operator (LASSO) algorithm was used to select radiomic features and develop a radiomic signature score (Rad-score). The training set was used to establish five machine-learning prediction models and the test set was used to test the efficacy of the models. The effectiveness of the models was compared by ROC analysis.

RESULTS

The CRC patients were divided into a training set (n = 144) and a test set (n = 55). Two radiomic features were selected to build the Rad-score. Five machine-learning algorithms including logistic regression, support vector machine (SVM), random forest, neural network and eXtreme gradient boosting (XGBoost) were used to established models. Among the five machine-learning models, logistic regression (AUC 0.866, 95% CI 0.808-0.925) and XGBoost (AUC 0.903, 95% CI 0.855-0.951) models performed the best. In the training set, the AUC of these two models were significantly higher than that of the LN metastasis status reported by 18F-FDG PET/CT for differentiating positive and negative regional LN metastases in CRC (all p < 0.05). Good efficacy of the above two models was also achieved in the test set. We created a nomogram based on the logistic regression model that visualized the results and provided an easy-to-use method for predicting regional LN metastasis in patients with CRC.

CONCLUSION

In this study, five machine-learning models for preoperative prediction of regional LN metastasis of CRC based on 18F-FDG PET/CT and PET-based radiomic features were successfully developed and validated. Among them, the logistic regression and XGBoost models performed the best, with higher efficacy than 18F-FDG PET/CT in both the training and test sets.

Metabolic Classification and Intervention Opportunities for Tumor Energy Dysfunction

A comprehensive view of cell metabolism provides a new vision of cancer, conceptualized as tissue with cellular-altered metabolism and energetic dysfunction, which can shed light on pathophysiological mechanisms. Cancer is now considered a heterogeneous ecosystem, formed by tumor cells and the microenvironment, which is molecularly, phenotypically, and metabolically reprogrammable. A wealth of evidence confirms metabolic reprogramming activity as the minimum common denominator of cancer, grouping together a wide variety of aberrations that can affect any of the different metabolic pathways involved in cell physiology. This forms the basis for a new proposed classification of cancer according to the altered metabolic pathway(s) and degree of energy dysfunction. Enhanced understanding of the metabolic reprogramming pathways of fatty acids, amino acids, carbohydrates, hypoxia, and acidosis can bring about new therapeutic intervention possibilities from a metabolic perspective of cancer.

Identification of Tumor-Specific MRI Biomarkers Using Machine Learning (ML)

The identification of reliable and non-invasive oncology biomarkers remains a main priority in healthcare. There are only a few biomarkers that have been approved as diagnostic for cancer. The most frequently used cancer biomarkers are derived from either biological materials or imaging data. Most cancer biomarkers suffer from a lack of high specificity. However, the latest advancements in machine learning (ML) and artificial intelligence (AI) have enabled the identification of highly predictive, disease-specific biomarkers. Such biomarkers can be used to diagnose cancer patients, to predict cancer prognosis, or even to predict treatment efficacy. Herein, we provide a summary of the current status of developing and applying Magnetic resonance imaging (MRI) biomarkers in cancer care. We focus on all aspects of MRI biomarkers, starting from MRI data collection, preprocessing and machine learning methods, and ending with summarizing the types of existing biomarkers and their clinical applications in different cancer types.