Tumor heterogeneity assessed by texture analysis on contrast-enhanced CT in lung adenocarcinoma: association with pathologic grade

Subsolid Nodules: Significance and Current Understanding

Subsolid nodules are heterogeneously appearing and behaving entities, commonly encountered incidentally and in high-risk populations. Accurate characterization of subsolid nodules, and application of evolving surveillance guidelines, facilitates evidence-based and multidisciplinary patient-centered management.

MRI Radiomics Features of Adenohypophysis Determine the Activation of Hypothalamic-Pituitary-Gonadal Axis in Peri-Puberty Children

BACKGROUND

The status of the hypothalamic-pituitary-gonadal (HPG) axis is important for assessing the onset of physiological or pathological puberty. The reference standard gonadotropin-releasing hormone (GnRH) stimulation test requires hospital admission and repeated blood samples. A simple noninvasive method would be beneficial.

OBJECTIVES

To explore a noninvasive method for evaluating HPG axis activation in children using an MRI radiomics model.

STUDY TYPE

Retrospective.

POPULATION

Two hundred thirty-nine children (83 male; 3.6-14.6 years) with hypophysial MRI and GnRH stimulation tests, randomly divided a training set (168 children) and a test set (71 children).

FIELD STRENGTH/SEQUENCE

3.0 T, 3D isotropic fast spin echo (CUBE) T1-weighted imaging (T1WI) sequences.

ASSESSMENT

Radiomics features were extracted from sagittal 3D CUBE T1WI, and imaging signatures were generated using the least absolute shrinkage and selection operator (LASSO) with 10-fold cross-validation. Diagnostic performance for differential diagnosis of HPG status was compared between a radiomics model and MRI features (adenohypophyseal height [aPH] and volume [aPV]).

STATISTICAL TESTS

Receiver operating characteristic (ROC) and decision curve analysis (DCA). A P value <0.05 was considered statistically significant.

RESULTS

Eight hundred fifty-one radiomics features were extracted and reduced to 10 by the LASSO method in the training cohort. The radiomics model based on CUBE T1WI showed good performance in assessment of HPG axis activation with an area under the ROC curve (AUC) of 0.81 (95% CI: 0.71, 0.91) in the test set. The AUC of the radiomics model was significantly higher than that of aPH (0.81 vs. 0.65) but there was no significant difference compared to aPV (0.81 vs. 0.78, P = 0.58). In DCA analysis, the radiomics signature showed higher net benefit over the aPV and aPH models.

DATA CONCLUSIONS

The MRI radiomics model has potential to assess HPG axis activation status noninvasively, potentially providing valuable information in the diagnosis of patients with pathological puberty onset.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 2.

Multidimensional progressive single-cell sequencing reveals cell microenvironment composition and cancer heterogeneity in lung cancer

Despite substantial advances in cancer biology and treatment, the clinical outcomes of patients with lung cancer remain unsatisfactory. The tumor microenvironment (TME) is a potential target. Using single-cell RNA sequencing, we could distinguish eight distinct cell types in the lung cancer microenvironment, demonstrating substantial intratumoral heterogeneity in 19 different lung cancer tumor samples. Through the re-dimensional grouping of cancer-associated fibroblasts (CAFs), myeloid cells, epithelial cells, natural killer (NK) cells, and T cells, the difference in the TME of lung cancer was revealed. We discovered SFTPB, SFN, and KRT8 as possible predictive biomarkers for lung cancer by assessing the gene expression patterns in epithelial cells. Examining cell-to-cell communications showed a robust association between the quantity of matrix CAFs, epithelial cells, and macrophages in the thrombospondin signaling pathway. Additionally, we found that the amyloid precursor protein signaling pathway primarily originated from the matrix, and inflammatory cancer-associated endothelial and fibroblast cells showed a co-expression relationship with myeloid cells and B cells. Through cell-to-cell correlation analysis, we found positive regulation between NK cells, regulatory T cells, GZMB-CD8 T cells, and GZMK-CD8 T cells, which could play a role in developing immune TMEs. These findings support studies on cancer heterogeneity and add to our understanding of lung cancer's cellular microenvironment.

Differentiation Between Heterogeneous GGN and Part-Solid Nodule Using 2 D Grayscale Histogram Analysis of Thin-Section CT Image

INTRODUCTION/BACKGROUND

To evaluate cases of surgically resected pulmonary adenocarcinoma (Ad) with heterogenous ground-glass nodules (HGGNs) or part-solid nodules (PSNs) and to clarify the differences between them, and between invasive adenocarcinoma (IVA) and minimally invasive adenocarcinoma (MIA) + adenocarcinoma in situ (AIS) using grayscale histogram analysis of thin-section computed tomography (TSCT).

MATERIALS AND METHODS

241 patients with pulmonary Ad were retrospectively classified into HGGNs and PSNs on TSCT by three thoracic radiologists. Sixty HGGNs were classified into 17 IVAs, 26 MIAs, and 17 AISs. 181 PSNs were classified into 114 IVAs, 55 MIAs, and 12 AISs.

RESULTS

We found significant differences in area (P = 0.0024), relative size of solid component (P <0.0001), circumference (P <0.0001), mean CT value (P <0.0001), standard deviation of the CT value (P <0.0001), maximum CT value (P <0.0001), skewness (P <0.0001), kurtosis (P <0.0001), and entropy (P <0.0001) between HGGNs and PSNs. In HGGNs, we found significant differences in relative size of solid component (P <0.0001), mean CT value (P = 0.0005), standard deviation of CT value (P = 0.0071), maximum CT value (P = 0.0237), and skewness (P = 0.0027) between IVAs and MIA+AIS lesions. In PSNs, we found significant differences in area (P = 0.0029), relative size of solid component (P = 0.0003), circumference (P = 0.0004), mean CT value (P = 0.0011), skewness (P = 0.0009), and entropy (P = 0.0002) between IVAs and the MIA+AIS lesions.

CONCLUSION

Quantitative evaluations using grayscale histogram analysis can clearly distinguish between HGGNs and PSNs, and may be useful for estimating the pathology of such lesions.

Use of machine learning to assess the prognostic utility of radiomic features for in-hospital COVID-19 mortality

As portable chest X-rays are an efficient means of triaging emergent cases, their use has raised the question as to whether imaging carries additional prognostic utility for survival among patients with COVID-19. This study assessed the importance of known risk factors on in-hospital mortality and investigated the predictive utility of radiomic texture features using various machine learning approaches. We detected incremental improvements in survival prognostication utilizing texture features derived from emergent chest X-rays, particularly among older patients or those with a higher comorbidity burden. Important features included age, oxygen saturation, blood pressure, and certain comorbid conditions, as well as image features related to the intensity and variability of pixel distribution. Thus, widely available chest X-rays, in conjunction with clinical information, may be predictive of survival outcomes of patients with COVID-19, especially older, sicker patients, and can aid in disease management by providing additional information.

A comprehensive nomogram combining CT-based radiomics with clinical features for differentiation of benign and malignant lung subcentimeter solid nodules

Objective

To establish a nomogram based on non-enhanced computed tomography(CT) imaging radiomics and clinical features for use in predicting the malignancy of sub-centimeter solid nodules (SCSNs).

Materials and methods

Retrospective analysis was performed of records for 198 patients with SCSNs that were surgically resected and examined pathologically at two medical institutions between January 2020 and June 2021. Patients from Center 1 were included in the training cohort (n = 147), and patients from Center 2 were included in the external validation cohort (n = 52). Radiomic features were extracted from chest CT images. The least absolute shrinkage and selection operator (LASSO) regression model was used for radiomic feature extraction and computation of radiomic scores. Clinical features, subjective CT findings, and radiomic scores were used to build multiple predictive models. Model performance was examined by evaluating the area under the receiver operating characteristic curve (AUC). The best model was selected for efficacy evaluation in a validation cohort, and column line plots were created.

Results

Pulmonary malignant nodules were significantly associated with vascular alterations in both the training (p &lt; 0.001) and external validation (p &lt; 0.001) cohorts. Eleven radiomic features were selected after a dimensionality reduction to calculate the radiomic scores. Based on these findings, three prediction models were constructed: subjective model (Model 1), radiomic score model (Model 2), and comprehensive model (Model 3), with AUCs of 0.672, 0.888, and 0.930, respectively. The optimal model with an AUC of 0.905 was applied to the validation cohort, and decision curve analysis indicated that the comprehensive model column line plot was clinically useful.

Conclusion

Predictive models constructed based on CT-based radiomics with clinical features can help clinicians diagnose pulmonary nodules and guide clinical decision making.

The adding value of contrast-enhanced CT radiomics: Differentiating tuberculosis from non-tuberculous infectious lesions presenting as solid pulmonary nodules or masses

Purpose

To compare the value of contrast-enhanced CT (CECT) and non-contrast-enhanced CT (NCECT) radiomics models in differentiating tuberculosis (TB) from non-tuberculous infectious lesions (NTIL) presenting as solid pulmonary nodules or masses, and develop a combine radiomics model (RM).

Materials and methods

This study was a retrospective analysis of 101 lesions in 95 patients, including 49 lesions (from 45 patients) in the TB group and 52 lesions (from 50 patients) in the NTIL group. Lesions were randomly divided into training and test sets in the ratio of 7:3. Conventional imaging features were used to construct a conventional imaging model (IM). Radiomics features screening and NCECT or CECT RM construction were carried out by correlation analysis and gradient boosting decision tree, and logistic regression. Finally, conventional IM, NCECT RM, and CECT RM were used for combine RM construction. Additionally, we recruited three radiologists for independent diagnosis. The differential diagnostic performance of each model was assessed using the areas under the receiver operating characteristic curve (AUCs).

Results

The CECT RM (training AUC, 0.874; test AUC, 0.796) outperformed the conventional IM (training AUC, 0.792; test AUC, 0.708), the NCECT RM (training AUC, 0.835; test AUC, 0.704), and three radiologists. The diagnostic efficacy of the combine RM (training AUC, 0.922; test AUC, 0.833) was best in the training and test sets.

Conclusions

The diagnostic efficacy of the CECT RM was superior to that of the NCECT RM in identifying TB from NTIL presenting as solid pulmonary nodules or masses. The combine RM had the best performance and may outperform expert radiologists.

Computed tomography radiomics-based distinction of invasive adenocarcinoma from minimally invasive adenocarcinoma manifesting as pure ground-glass nodules with bubble-like signs

BACKGROUND

To explore an effective model based on radiomics features extracted from nonenhanced computed tomography (CT) images to distinguish invasive adenocarcinoma (IAC) from minimally invasive adenocarcinoma (MIA) presenting as pure ground-glass nodules (pGGNs) with bubble-like (B-pGGNs) signs.

PATIENTS AND METHODS

We retrospectively reviewed 511 nodules (MIA, n = 288; IAC, n = 223) between November 2012 and June 2018 from almost all pGGNs pathologically confirmed MIA or IAC. Eventually, a total of 109 B-pGGNs (MIA, n = 55; IAC, n = 54) from 109 patients fulfilling the criteria were randomly assigned to the training and test cluster at a ratio of 7:3. The gradient boosting decision tree (GBDT) method and logistic regression (LR) analysis were applied to feature selection (radiomics, semantic, and conventional CT features). LR was performed to construct three models (the conventional, radiomics and combined model). The performance of the predictive models was evaluated using the area under the curve (AUC).

RESULTS

The radiomics model had good AUCs of 0.947 in the training cluster and of 0.945 in the test cluster. The combined model produced an AUC of 0.953 in the training cluster and of 0.945 in the test cluster. The combined model yielded no performance improvement (vs. the radiomics model). The rad_score was the only independent predictor of invasiveness.

CONCLUSION

The radiomics model showed excellent predictive performance in discriminating IAC from MIA presenting as B-pGGNs and may provide a necessary reference for extending clinical practice.

Progression-Free Survival Prediction in Small Cell Lung Cancer Based on Radiomics Analysis of Contrast-Enhanced CT

Purposes and Objectives

The aim of this study was to predict the progression-free survival (PFS) in patients with small cell lung cancer (SCLC) by radiomic signature from the contrast-enhanced computed tomography (CT).

Methods

A total of 186 cases with pathological confirmed small cell lung cancer were retrospectively assembled. First, 1,218 radiomic features were automatically extracted from tumor region of interests (ROIs) on the lung window and mediastinal window, respectively. Then, the prognostic and robust features were selected by machine learning methods, such as (1) univariate analysis based on a Cox proportional hazard (CPH) model, (2) redundancy removing using the variance inflation factor (VIF), and (3) multivariate importance analysis based on random survival forests (RSF). Finally, PFS predictive models were established based on RSF, and their performances were evaluated using the concordance index (C-index) and the cumulative/dynamic area under the curve (C/D AUC).

Results

In total, 11 radiomic features (6 for mediastinal window and 5 for lung window) were finally selected, and the predictive model constructed from them achieved a C-index of 0.7531 and a mean C/D AUC of 0.8487 on the independent test set, better than the predictions by single clinical features (C-index = 0.6026, mean C/D AUC = 0.6312), and single radiomic features computed in lung window (C-index = 0.6951, mean C/D AUC = 0.7836) or mediastinal window (C-index = 0.7192, mean C/D AUC = 0.7964).

Conclusion

The radiomic features computed from tumor ROIs on both lung window and mediastinal window can predict the PFS for patients with SCLC by a high accuracy, which could be used as a useful tool to support the personalized clinical decision for the diagnosis and patient management of patients with SCLC.

Analysis of KRAS Mutation Status Prediction Model for Colorectal Cancer Based on Medical Imaging

This study retrospectively included some patients with colorectal cancer diagnosed by histopathology, to explore the feasibility of CT medical image texture analysis in predicting KRAS gene mutations in patients with colorectal cancer. Before any surgical procedure, all patients received an enhanced CT scan of the abdomen and pelvis, as well as genetic testing. To define patient groups, divide all patients into test and validation sets based on the order of patient enrollment. A radiologist took a look at the plain axial CT image of the tumor, as well as the portal vein CT image, at the corresponding level. The physician points the computer's cursor to the relevant area in the image, and TexRAD software programs together texture parameters based on various spatial scale factors, also known as total mean, total variance, statistical entropy, overall total average, mean total, positive mean, skewness value, kurtosis value, and general skewness. Using the same method again two weeks later, the observer and another physician measured the image of each patient again to see if the method was consistent between observers. With regard to clinical information, the KRAS gene mutation group and the wild group of participants in the test set and validation set each had values for the texture parameter. In a study of patients with colorectal cancer, the results demonstrated that CT texture parameters were correlated with the presence of the KRAS gene mutation. The best CT prediction model includes the values of the medium texture image's slope and the other CT fine texture image's value of entropy, the medium texture image's slope and kurtosis, and the medium texture image's mean and the other CT fine texture image's value of entropy. Regardless of the training set or the validation set, patients with and without KRAS gene mutations did not differ significantly in clinical characteristics. This method can be used to identify mutations in the KRAS gene in patients with colorectal cancer, making it practical to implement CT medical image texture analysis technology for that purpose.

CT Texture Analysis of Pulmonary Neuroendocrine Tumors-Associations with Tumor Grading and Proliferation

Texture analysis derived from computed tomography (CT) might be able to provide clinically relevant imaging biomarkers and might be associated with histopathological features in tumors. The present study sought to elucidate the possible associations between texture features derived from CT images with proliferation index Ki-67 and grading in pulmonary neuroendocrine tumors. Overall, 38 patients (n = 22 females, 58%) with a mean age of 60.8 ± 15.2 years were included into this retrospective study. The texture analysis was performed using the free available Mazda software. All tumors were histopathologically confirmed. In discrimination analysis, "S(1,1)SumEntrp" was significantly different between typical and atypical carcinoids (mean 1.74 ± 0.11 versus 1.79 ± 0.14, p = 0.007). The correlation analysis revealed a moderate positive association between Ki-67 index with the first order parameter kurtosis (r = 0.66, p = 0.001). Several other texture features were associated with the Ki-67 index, the highest correlation coefficient showed "S(4,4)InvDfMom" (r = 0.59, p = 0.004). Several texture features derived from CT were associated with the proliferation index Ki-67 and might therefore be a valuable novel biomarker in pulmonary neuroendocrine tumors. "Sumentrp" might be a promising parameter to aid in the discrimination between typical and atypical carcinoids.

Texture features of primary tumor on 18F-FDG PET images in non-small cell lung cancer: The relationship between imaging and histopathological parameters

PURPOSE

The aims of this study were to evaluate the relationships between textural features of the primary tumor on FDG PET images and clinical-histopathological parameters which are useful in predicting prognosis in newly diagnosed non-small cell lung cancer (NSCLC) patients.

METHODS

PET/CT images of ninety (90) patients with NSCLC prior to surgery were analyzed retrospectively. All patients had resectable tumors. From the images we acquired data related to metabolism (SUVmax, MTV, TLG) and texture features of primary tumors. Histopathological tumor types and subgroups, degree of Ki-67 expression and necrosis rates of the primary tumor, mediastinal lymph node (MLN) status and nodal stages were recorded.

RESULTS

Among the two histologic tumor types (adenocarcinoma and squamous cell carcinoma) significant differences were present regarding metabolic parameters, Ki-67 index with higher values and kurtosis with lower values in the latter group. Textural heterogeneity was found to be higher in poorly differentiated tumors compared to moderately differentiated tumors in patients with adenocarcinoma. While Ki-67 index had significant correlations with metabolic parameters and kurtosis, tumor necrosis rate was only significantly correlated with textural features. By univariate and multivariate analyses of the imaging and histopathological factors examined, only gradient variance was significant predictive factor for the presence of MLN metastasis.

CONCLUSIONS

Textural features had significant associations with histologic tumor types, degree of pathological differentiation, tumor proliferation and necrosis rates. Texture analysis has potential to differentiate tumor types and subtypes and to predict MLN metastasis in patients with NSCLC.

CT texture analysis of tonsil cancer: Discrimination from normal palatine tonsils

The purposes of the study were to determine whether there are differences in texture analysis parameters between tonsil cancers and normal tonsils, and to correlate texture analysis with 18F-FDG PET/CT to investigate the relationship between texture analysis and metabolic parameters. Sixty-four patients with squamous cell carcinoma of the palatine tonsil were included. A ROI was drawn, including all slices, to involve the entire tumor. The contralateral normal tonsil was used for comparison with the tumors. Texture analysis parameters, mean, standard deviation (SD), entropy, mean positive pixels, skewness, and kurtosis were obtained using commercially available software. Parameters were compared between the tumor and the normal palatine tonsils. Comparisons were also performed among early tonsil cancer, advanced tonsil cancer, and normal tonsils. An ROC curve analysis was performed to assess discrimination of tumor from normal tonsils. Correlation between texture analysis and 18F-FDG PET/CT was performed. Compared to normal tonsils, the tumors showed a significantly lower mean, higher SD, higher entropy, lower skewness, and higher kurtosis on most filters (p<0.001). On comparisons among normal tonsils, early cancers, and advanced tonsil cancers, SD and entropy showed significantly higher values on all filters (p<0.001) between early cancers and normal tonsils. The AUC from the ROC analysis was 0.91, obtained from the entropy. A mild correlation was shown between texture parameters and metabolic parameters. The texture analysis parameters, especially entropy, showed significant differences in contrast-enhanced CT results between tumor and normal tonsils, and between early tonsil cancers and normal tonsils. Texture analysis can be useful as an adjunctive tool for the diagnosis of tonsil cancers.

Reproducibility of CT-based radiomic features against image resampling and perturbations for tumour and healthy kidney in renal cancer patients

Computed Tomography (CT) is widely used in oncology for morphological evaluation and diagnosis, commonly through visual assessments, often exploiting semi-automatic tools as well. Well-established automatic methods for quantitative imaging offer the opportunity to enrich the radiologist interpretation with a large number of radiomic features, which need to be highly reproducible to be used reliably in clinical practice. This study investigates feature reproducibility against noise, varying resolutions and segmentations (achieved by perturbing the regions of interest), in a CT dataset with heterogeneous voxel size of 98 renal cell carcinomas (RCCs) and 93 contralateral normal kidneys (CK). In particular, first order (FO) and second order texture features based on both 2D and 3D grey level co-occurrence matrices (GLCMs) were considered. Moreover, this study carries out a comparative analysis of three of the most commonly used interpolation methods, which need to be selected before any resampling procedure. Results showed that the Lanczos interpolation is the most effective at preserving original information in resampling, where the median slice resolution coupled with the native slice spacing allows the best reproducibility, with 94.6% and 87.7% of features, in RCC and CK, respectively. GLCMs show their maximum reproducibility when used at short distances.

Texture features of primary tumor on 18F-FDG PET images in non-small cell lung cancer: The relationship between imaging and histopathological parameters

OBJECTIVES

The aims of this study were to evaluate the relationships between textural features of the primary tumor on FDG PET images and clinical-histopathological parameters which are useful in predicting prognosis in newly diagnosed non-small cell lung cancer (NSCLC) patients.

MATERIAL AND METHODS

PET/CT images of ninety (90) patients with NSCLC prior to surgery were analyzed retrospectively. All patients had resectable tumors. From the images we acquired data related to metabolism (SUVmax, metabolic tumor volume [MTV] and total lesion glycolysis [TLG]) and texture features of primary tumors. Histopathological tumor types and subgroups, degree of Ki-67 expression and necrosis rates of the primary tumor, mediastinal lymph node (MLN) status and nodal stages were recorded.

RESULTS

Among the 2histologic tumor types (adenocarcinoma and squamous cell carcinoma) significant differences were present regarding metabolic parameters, Ki-67 index with higher values and kurtosis with lower values in the latter group. Textural heterogeneity was found to be higher in poorly differentiated tumors compared to moderately differentiated tumors in patients with adenocarcinoma. While Ki-67 index had significant correlations with metabolic parameters and kurtosis, tumor necrosis rate was only significantly correlated with textural features. By univariate and multivariate analyses of the imaging and histopathological factors examined, only gradient variance was significant predictive factor for the presence of MLN metastasis.

CONCLUSIONS

Textural features had significant associations with histologic tumor types, degree of pathological differentiation, tumor proliferation and necrosis rates. Texture analysis has potential to differentiate tumor types and subtypes and to predict MLN metastasis in patients with NSCLC.

The Biological Meaning of Radiomic Features

An earlier incorrect version appeared online. This article was corrected on February 10, 2021.

Contrast-enhanced computed tomography prior to percutaneous transthoracic needle biopsy reduces the incidence of hemorrhage

Background

Hemorrhage is the second most common complication of percutaneous transthoracic needle biopsy (PTNB), and at present, there is no effective prevention strategy. Contrast-enhanced computed tomography (CECT) has the advantage of clearly visualizing blood supply within the lesion and aiding in the imaging of blood vessels, which can reduce hemorrhage complicating PTNB. As no large-sample studies were evaluating whether CECT could reduce hemorrhage, we conducted the present retrospective study.

Methods

From November 2011 to February 2016, 1,282 biopsies at Jinling Hospital were retrospectively reviewed; 555 underwent CECT, and 727 underwent non-contrast computed tomography (CT). Factors associated with hemorrhage were defined, and hemorrhage rates were compared between the 2 groups.

Results

We found that pre-biopsy CECT was associated with a reduced incidence of biopsy-related hemorrhage compared to non-contrast CT (16.4% vs. 23.1%, P=0.003). Propensity score matching (PSM) analysis also showed that the incidence of hemorrhage in the CECT group was lower than that of the non-contrast CT group at a ratio of 1:1 (P=0.039), 1:2 (P=0.028), or 1:3 (P=0.013). In the multivariate analysis, CECT before PTNB was found to be significantly associated with a reduced risk of hemorrhage [odds ratio (OR): 0.671, 95% confidence interval (CI): 0.499–0.902, P=0.008]. Puncture position, lesion size, depth of needle tract, and the number of punctures were also found to be associated with hemorrhage (all P<0.05).

Conclusions

Compared with non-contrast CT, CECT significantly reduced the risk of post-biopsy pulmonary hemorrhage, which suggests that CECT should be performed before PTNB.

Usefulness of texture analysis for grading pancreatic neuroendocrine tumors on contrast-enhanced computed tomography and apparent diffusion coefficient maps

PURPOSE

To determine whether texture analysis of contrast-enhanced computed tomography (CECT) and apparent diffusion coefficient (ADC) maps could predict tumor grade (G1 vs G2-3) in patients with pancreatic neuroendocrine tumor (PNET).

MATERIALS AND METHODS

Thirty-three PNETs (22 G1 and 11 G2-3) were retrospectively reviewed. Fifty features were individually extracted from the arterial and portal venous phases of CECT and ADC maps by two radiologists. Diagnostic performance was assessed by receiver operating characteristic curves while inter-observer agreement was determined by calculating intraclass correlation coefficients (ICCs).

RESULTS

G2-G3 tumors were significantly larger than G1. Seventeen features significantly differed among the two readers on univariate analysis, with ICCs > 0.6; the largest area under the curve (AUC) for features of each CECT phase and ADC map was log-sigma 1.0 joint-energy = 0.855 for the arterial phase, log-sigma 1.5 kurtosis = 0.860 for the portal venous phase, and log-sigma 1.0 correlation = 0.847 for the ADC map. The log-sigma 1.5 kurtosis of the portal venous phase showed the largest AUC in the CECT and ADC map, and its sensitivity, specificity, and accuracy were 95.5%, 72.7%, and 87.9%, respectively.

CONCLUSION

Texture analysis may aid in differentiating between G1 and G2-3 PNET.

Predicting the invasiveness of lung adenocarcinomas appearing as ground-glass nodule on CT scan using multi-task learning and deep radiomics

Background

Due to different treatment method and prognosis of different subtypes of lung adenocarcinomas appearing as ground-glass nodules (GGNs) on computed tomography (CT) scan, it is important to classify invasive adenocarcinomas from non-invasive adenocarcinomas. The purpose of this paper is to build and evaluate the performance of deep learning networks on the differentiation the invasiveness of lung adenocarcinoma appearing as GGNs.

Methods

This retrospective study included 886 GGNs from 794 pathological confirmed patients with lung adenocarcinoma for training and testing the proposed networks. Three deep learning networks, namely XimaNet (deep learning-based classification model), XimaSharp (classification and nodule segmentation model), and Deep-RadNet (deep learning and radiomics combined classification model, i.e., deep radiomics) were built. Three classification tasks, namely task 1: classification of AAH/AIS and MIA, task 2: classification of MIA and IAC, and task 3: classification of non-invasive adenocarcinomas and invasive adenocarcinomas (AAH/AIS&MIA and IAC) were conducted to evaluate the model performance. The Z-test was used to compare the model performance.

Results

The AUC for classification of AAH/AIS with MIA were 0.891, 0.841 and 0.779 for Deep-RadNet, XimaNet and XimaSharp respectively. The AUC for classification of MIA with IAC were 0.889, 0.785 and 0.778 for three networks and AUC for classification of AAH/AIS&MIA with IAC were 0.941, 0.892 and 0.827 respectively. The performance of deep_RadNet was better than the other two models with the Z-test (P<0.05).

Conclusions

Deep-RadNet with the visual heat map could evaluate the invasiveness of GGNs accurately and intuitively, providing a theoretical basis for individualized and accurate medical treatment of patients with GGNs.

Differentiating Peripherally-Located Small Cell Lung Cancer From Non-small Cell Lung Cancer Using a CT Radiomic Approach

Lung cancer can be classified into two main categories: small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC), which are different in treatment strategy and survival probability. The lung CT images of SCLC and NSCLC are similar such that their subtle differences are hardly visually discernible by the human eye through conventional imaging evaluation. We hypothesize that SCLC/NSCLC differentiation could be achieved via computerized image feature analysis and classification in feature space, as termed a radiomic model. The purpose of this study was to use CT radiomics to differentiate SCLC from NSCLC adenocarcinoma. Patients with primary lung cancer, either SCLC or NSCLC adenocarcinoma, were retrospectively identified. The post-diagnosis pre-treatment lung CT images were used to segment the lung cancers. Radiomic features were extracted from histogram-based statistics, textural analysis of tumor images and their wavelet transforms. A minimal-redundancy-maximal-relevance method was used for feature selection. The predictive model was constructed with a multilayer artificial neural network. The performance of the SCLC/NSCLC adenocarcinoma classifier was evaluated by the area under the receiver operating characteristic curve (AUC). Our study cohort consisted of 69 primary lung cancer patients with SCLC (n = 35; age mean ± SD = 66.91± 9.75 years), and NSCLC adenocarcinoma (n = 34; age mean ± SD = 58.55 ± 11.94 years). The SCLC group had more male patients and smokers than the NSCLC group (P < 0.05). Our SCLC/NSCLC classifier achieved an overall performance of AUC of 0.93 (95% confidence interval = [0.85, 0.97]), sensitivity = 0.85, and specificity = 0.85). Adding clinical data such as smoking history could improve the performance slightly. The top ranking radiomic features were mostly textural features. Our results showed that CT radiomics could quantitatively represent tumor heterogeneity and therefore could be used to differentiate primary lung cancer subtypes with satisfying results. CT image processing with the wavelet transformation technique enhanced the radiomic features for SCLC/NSCLC classification. Our pilot study should motivate further investigation of radiomics as a non-invasive approach for early diagnosis and treatment of lung cancer.

Pulmonary adenocarcinoma appearing as ground-glass opacity nodules identified using non-enhanced and contrast-enhanced CT texture analysis: A retrospective analysis

The present study aimed to investigate the ability of CT-based texture analysis to differentiate invasive adenocarcinoma (IA) from pre-invasive lesions (PIL) or minimally IA (MIA) appearing as ground-glass opacity (GGO) nodules, and to further compare the performance of non-enhanced CT (NECT) images with that of contrast-enhanced CT (CECT) images. A total of 77 patients with GGO nodules and surgically confirmed pulmonary adenocarcinoma were included in the present retrospective study. Each GGO nodule was manually segmented and its texture features were extracted from NECT and CECT images using in-house developed software coded in MATLAB (MathWorks). The independent-samples t-test was used to select the texture features with statistically significant differences between IA and MIA/PIL. Multivariate logistic regression and receiver operating characteristics (ROC) curve analyses were performed to identify predictive features. Of the 77 GGO nodules, 12 were atypical adenomatous hyperplasia or adenocarcinoma in situ (15.6%), 36 were MIA (46.8%) and 29 were IA (37.7%). IA and MIA/PIL exhibited significant differences in most histogram features and gray-level co-occurrence matrix features (P<0.05). Multivariate logistic regression and ROC curve analyses revealed that smaller energy and higher entropy were significant differentiators of IA from MIA and PIL, irrespective of whether NECT images [area under the curve (AUC): 0.839, 0.859] or CECT images (AUC: 0.818, 0.820) are used. Texture analysis of CT images, regardless of whether NECT or CECT is used, has the potential to distinguish IA from PIL or MIA, particularly the parameters of energy and entropy. Furthermore, NECT images were simpler to obtain and no contrast agent was required; thus, analysis with NECT may be a preferred choice.

Radiomics predicts survival of patients with advanced non-small cell lung cancer undergoing PD-1 blockade using Nivolumab

Immune checkpoint blockade is an emerging anticancer strategy, and Nivolumab is a human mAb to PD-1 that is used in the treatment of a number of different malignancies, including non-small cell lung cancer (NSCLC), kidney cancer, urothelial carcinoma and melanoma. Although the use of Nivolumab prolongs survival in a number of patients, this treatment is hampered by high cost. Therefore, the identification of predictive markers of response to treatment in patients is required. In this context, PD-1/PDL1 blockade antitumor effects occur through the reactivation of a pre-existing immune response, and the efficacy of these effects is strictly associated with the presence of necrosis, hypoxia and inflammation at the tumour sites. It has been indicated that these events can be evaluated by specific assessments using a computed tomography (CT) texture analysis (TA) or radiomics. Therefore, a retrospective study was performed, which aimed to evaluate the potential use of this analysis in the identification of patients with NSCLC who may benefit from Nivolumab treatment. A retrospective analysis was performed of 59 patients with metastatic NSCLC who received Nivolumab treatment between January 2015 and July 2017 at Siena University Hospital (35 patients, training dataset), Catanzaro University Hospital and Reggio Calabria Grand Metropolitan Hospital, Italy (24 patients, validation dataset). Pre- and post-contrast CT sequences were used to contour the gross tumour volume (GTV) of the target lesions prior to Nivolumab treatment. The impact of variations on contouring was analysed using two delineations, which were performed on each patient, and the TA parameters were tested for reliability using the Intraclass Coefficient Correlation method (ICC). All analyses for the current study were performed using LifeX Software^©. Imaging, clinical and pathological parameters were correlated with progression free survival and overall survival (OS) using Kaplan Meier analysis. An external validation testing was performed for the TA Score using the validation dataset. A total of 59 patients were included in the analysis of the present study. The reliability ICC analysis of 14 TA parameters indicated a highly reproducibility (ICC >0.70, single measure) in 12 (85%) pre- contrast and 13 (93%) post-contrast exams. A specific cut-off was detected for each of the following parameters: volume (score 1 >36 ml), histogram entropy (score 1 > 1.30), compacity (score 1 <3), gray level co-occurrence matrix (GLCM)-entropy (score 1 >1.80), GLCM-Dissimilarity (score 1 >5) and GLCM-Correlation (score 1<0.54). The global texture score allowed the classification of two subgroups of Low (Score 0-1; 36 patients; 61%) and High Risk patients (Score >1; 23 patients; 39%) that respectively, showed a median OS of 26 (mean +/- SD: 18 +/- 1.98 months; 95% CI 14-21 months) and 5 months (mean +/- SD: 6 +/- 0.99 months; 95% CI: 4-8 months; P=0.002). The current study indicated that TA parameters can identify patients that will benefit from PD-1 blockage by defining the radiological settings that are potentially suggestive of an active immune response. These results require further confirmation in prospective trials.

Preoperative CT texture features predict prognosis after curative resection in pancreatic cancer

Pancreatic cancer is a lethal disease, and resistance to chemotherapy is a critical factor influencing the postoperative prognosis. Tumour heterogeneity is an important indicator of chemoresistance. Therefore, we analysed tumour heterogeneity in preoperative computed tomography scans by performing texture analysis using the grey-level run-length matrix and analysed the correlation of survival with the value obtained in these analyses. We analysed 116 consecutive patients who underwent curative resection and had preoperative contrast-enhanced computed tomography data available for analysis. A region of interest was drawn on all slices with a visible tumour and normal pancreas on the arterial phase computed tomography scans; the correlation of pathological characteristics with grey-level run-length matrix features was analysed. We then performed Kaplan–Meier survival curve analysis among pancreatic cancer patients. The grey-level non-uniformity values in grey-level run-length matrix features for tumours were higher than those for normal pancreas. High grey-level non-uniformity values represent a non-uniform texture, i.e., heterogeneity. Grey-level run-length matrix features showed that recurrence-free survival was shorter in the group with high grey-level non-uniformity 135 values (p = 0.025). Our analyses of the correlation between pathological outcomes and grey-level run-length matrix features in pancreatic cancer patients showed that grey-level non-uniformity values were powerful prognostic indicators.

Radiomics for Classification of Lung Cancer Histological Subtypes Based on Nonenhanced Computed Tomography

OBJECTIVES

To evaluate the performance of using radiomics method to classify lung cancer histological subtypes based on nonenhanced computed tomography images.

MATERIALS AND METHODS

278 patients with pathologically confirmed lung cancer were collected, including 181 nonsmall cell lung cancer (NSCLC) and 97 small cell lung cancers (SCLC) patients. Among the NSCLC patients, 88 patients were adenocarcinomas (AD) and 93 patients were squamous cell carcinomas (SCC). In total, 1695 quantitative radiomic features (QRF) were calculated from the primary lung cancer tumor in each patient. To build radiomic classification model based on the extracted QRFs, several machine-learning algorithms were applied sequentially. First, unsupervised hierarchical clustering was used to exclude highly correlated QRFs; second, the minimum Redundancy Maximum Relevance feature selection algorithm was employed to select informative and nonredundant QRFs; finally, the Incremental Forward Search and Support Vector Machine classification algorithms were used to combine the selected QRFs and build the model. In our work, to study the phenotypic differences among lung cancer histological subtypes, four classification models were built. They were models of SCLC vs NSCLC, SCLC vs AD, SCLC vs SCC, and AD vs SCC. The performance of the classification models was evaluated by the area under the receiver operating characteristic curve (AUC) estimated by three-fold cross-validation.

RESULTS

The AUC (95% confidence interval) for the model of SCLC vs NSCLC was 0.741(0.678, 0.795). For the models of SCLC vs AD and SCLC vs SCC, the AUCs were 0.822(0.755, 0.875) and 0.665(0.583, 0.738), respectively. The AUC for the model of AD vs SCC was 0.655(0.570, 0.731). Several QRFs ("Law_15," "LoG_Uniformity," "GLCM_Contrast," and "Compactness Factor") that characterize tumor heterogeneity and shape were selected as the significant features to build the models.

CONCLUSION

Our results show that phenotypic differences exist among different lung cancer histological subtypes on nonenhanced computed tomography image.

CT Texture Analysis-Correlations With Histopathology Parameters in Head and Neck Squamous Cell Carcinomas

Introduction: Texture analysis is an emergent imaging technique to quantify heterogeneity in radiological images. It is still unclear whether this technique is capable to reflect tumor microstructure. The present study sought to correlate histopathology parameters with texture features derived from contrast-enhanced CT images in head and neck squamous cell carcinomas (HNSCC).

Materials and Methods: Twenty-eight patients with histopathological proven HNSCC were retrospectively analyzed. In every case EGFR, VEGF, Hif1-alpha, Ki67, p53 expression derived from immunhistochemical specimen were semiautomatically calculated. Furthermore, mean cell count was estimated. Texture analysis was performed on contrast-enhanced CT images as a whole lesion measurement. Spearman's correlation analysis was performed, adjusted with Benjamini-Hochberg correction for multiple tests.

Results: Several texture features correlated with histopathological parameters. After correction only CT texture joint entropy and CT entropy correlation with Hif1-alpha expression remained statistically significant (ρ = −0.60 and ρ = −0.50, respectively).

Conclusions: CT texture joint entropy and CT entropy were associated with Hif1-alpha expression in HNSCC and might be able to reflect hypoxic areas in this entity.

Pancreatic neuroendocrine tumor: prediction of the tumor grade using magnetic resonance imaging findings and texture analysis with 3-T magnetic resonance

Purpose

The purpose of this study was to evaluate the performance of magnetic resonance imaging (MRI) findings and texture parameters for prediction of the histopathologic grade of pancreatic neuroendocrine tumors (PNETs) with 3-T magnetic resonance.

Patients and methods

PNETs are classified into Grade 1 (G1), Grade 2 (G2), and Grade 3 (G3) tumors based on the Ki-67 proliferation index and the mitotic activity. A total of 77 patients with pathologically confirmed PNETs met the inclusion criteria. Texture analysis (TA) was applied to T2-weighted imaging (T2WI) and diffusion-weighted imaging (DWI) maps. Patient demographics, MRI findings, and texture parameters were compared among three different histopathologic subtypes by using Fisher’s exact tests or Kruskal–Wallis test. Then, logistic regression analysis was adopted to predict tumor grades. ROC curves and AUCs were calculated to assess the diagnostic performance of MRI findings and texture parameters in prediction of tumor grades.

Results

There were 31 G1, 29 G2, and 17 G3 patients. Compared with G1, G2/G3 tumors showed higher frequencies of an ill-defined margin, a predominantly solid tumor type, local invasion or metastases, hypo-enhancement at the arterial phase, and restriction diffusion. Four T2-based (inverse difference moment, energy, correlation, and differenceEntropy) and five DWI-based (correlation, contrast, inverse difference moment, maxintensity, and entropy) TA parameters exhibited statistical significance among PNETs (P<0.001). The AUCs of six predicting models on T2WI and DWI ranged from 0.703–0.989.

Conclusion

Our data indicate that MRI findings, including tumor margin, texture, local invasion or metastases, tumor enhancement, and diffusion restriction, as well as texture parameters can aid the prediction of PNETs grading.

CT Enhancement and 3D Texture Analysis of Pancreatic Neuroendocrine Neoplasms

To evaluate pancreatic neuroendocrine neoplasms (panNENs) grade prediction by means of qualitative and quantitative CT evaluation, and 3D CT-texture analysis. Patients with histopathologically-proven panNEN, availability of Ki67% values and pre-treatment CT were included. CT images were retrospectively reviewed, and qualitative and quantitative images analysis were done; for quantitative analysis four enhancement-ratios and three permeability-ratios were created. 3D CT-texture imaging analysis was done (Mean Value; Variance; Skewness; Kurtosis; Entropy). Subsequently, these features were compared among the three grading (G) groups. 304 patients affected by panNENs were considered, and 100 patients were included. At qualitative evaluation, frequency of irregular margins was significantly different between tumor G groups. At quantitative evaluation, for all ratios, comparisons resulted statistical significant different between G1 and G3 groups and between G2 and G3 groups. At 3D CT-texture analysis, Kurtosis resulted statistical significant different among three G groups and Entropy resulted statistical significant different between G1 and G3 and between G2 and G3 groups. Quantitative CT evaluation of panNENs can predict tumor grade, discerning G1 from G3 and G2 from G3 tumors. CT-texture analysis can predict panNENs tumor grade, distinguishing G1 from G3 and G2 from G3, and G1 from G2 tumors.

Predicting locally advanced rectal cancer response to neoadjuvant therapy with 18F-FDG PET and MRI radiomics features

PURPOSE

Pathological complete response (pCR) following neoadjuvant chemoradiotherapy or radiotherapy in locally advanced rectal cancer (LARC) is reached in approximately 15-30% of cases, therefore it would be useful to assess if pretreatment of ¹⁸F-FDG PET/CT and/or MRI texture features can reliably predict response to neoadjuvant therapy in LARC.

METHODS

Fifty-two patients were dichotomized as responder (pR+) or non-responder (pR-) according to their pathological tumor regression grade (TRG) as follows: 22 as pR+ (nine with TRG = 1, 13 with TRG = 2) and 30 as pR- (16 with TRG = 3, 13 with TRG = 4 and 1 with TRG = 5). First-order parameters and 21 second-order texture parameters derived from the Gray-Level Co-Occurrence matrix were extracted from semi-automatically segmented tumors on T2w MRI, ADC maps, and PET/CT acquisitions. The role of each texture feature in predicting pR+ was assessed with monoparametric and multiparametric models.

RESULTS

In the mono-parametric approach, PET homogeneity reached the maximum AUC (0.77; sensitivity = 72.7% and specificity = 76.7%), while PET glycolytic volume and ADC dissimilarity reached the highest sensitivity (both 90.9%). In the multiparametric analysis, a logistic regression model containing six second-order texture features (five from PET and one from T2w MRI) yields the highest predictivity in distinguish between pR+ and pR- patients (AUC = 0.86; sensitivity = 86%, and specificity = 83% at the Youden index).

CONCLUSIONS

If preliminary results of this study are confirmed, pretreatment PET and MRI could be useful to personalize patient treatment, e.g., avoiding toxicity of neoadjuvant therapy in patients predicted pR-.

HRCT texture analysis for pure or part-solid ground-glass nodules: distinguishability of adenocarcinoma in situ or minimally invasive adenocarcinoma from invasive adenocarcinoma

PURPOSE

To distinguish between adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC) showing pure or part-solid ground-glass nodules (GGNs) by high-resolution computed tomography (HRCT) texture analysis.

MATERIALS AND METHODS

This retrospective study included 101 consecutive patients with 115 pure or part-solid GGNs ≤ 3 cm diameter, which were surgically resected and pathologically diagnosed with AIS, MIA, or IAC (48 AIS-MIA and 67 IAC) between April 2011 and March 2015. Each tumor was manually segmented on axial CT images, and the following texture features were calculated: volume, mass, mean CT value, variance, skewness, kurtosis, entropy, uniformity, and percentile CT numbers (10th, 25th, 50th, 75th, 90th, 95th percentiles). The differences between AIS-MIA and IAC were statistically evaluated using univariate, multivariate, and receiver operating characteristic analysis.

RESULTS

Compared with IAC, AIS-MIA had significantly greater skewness, kurtosis, and uniformity, whereas in the other parameters, AIS-MIA demonstrated significantly lower values than those of IAC. Multivariate analysis revealed that independent differentiators were the 90th percentile CT numbers (P < 0.001) and entropy (P = 0.005) with an excellent accuracy (area under the curve, 0.90).

CONCLUSIONS

The 90th percentile CT numbers and entropy can accurately distinguish AIS-MIA from IAC.

Predicting survival time of lung cancer patients using radiomic analysis

OBJECTIVES

This study investigates the prediction of Non-small cell lung cancer (NSCLC) patient survival outcomes based on radiomic texture and shape features automatically extracted from tumor image data.

MATERIALS AND METHODS

Retrospective analysis involves CT scans of 315 NSCLC patients from The Cancer Imaging Archive (TCIA). A total of 24 image features are computed from labeled tumor volumes of patients within groups defined using NSCLC subtype and TNM staging information. Spearman's rank correlation, Kaplan-Meier estimation and log-rank tests were used to identify features related to long/short NSCLC patient survival groups. Automatic random forest classification was used to predict patient survival group from multivariate feature data. Significance is assessed at P < 0.05 following Holm-Bonferroni correction for multiple comparisons.

RESULTS

Significant correlations between radiomic features and survival were observed for four clinical groups: (group, [absolute correlation range]): (large cell carcinoma (LCC) [0.35, 0.43]), (tumor size T2, [0.31, 0.39]), (non lymph node metastasis N0, [0.3, 0.33]), (TNM stage I, [0.39, 0.48]). Significant log-rank relationships between features and survival time were observed for three clinical groups: (group, hazard ratio): (LCC, 3.0), (LCC, 3.9), (T2, 2.5) and (stage I, 2.9). Automatic survival prediction performance (i.e. below/above median) is superior for combined radiomic features with age-TNM in comparison to standard TNM clinical staging information (clinical group, mean area-under-the-ROC-curve (AUC)): (LCC, 75.73%), (N0, 70.33%), (T2, 70.28%) and (TNM-I, 76.17%).

CONCLUSION

Quantitative lung CT imaging features can be used as indicators of survival, in particular for patients with large-cell-carcinoma (LCC), primary-tumor-sizes (T2) and no lymph-node-metastasis (N0).