Improving tumour heterogeneity MRI assessment with histograms. Br J Cancer. 2014;111:2205-2213. [PMID: 25268373 PMCID: PMC4264439 DOI: 10.1038/bjc.2014.512]

Prospective longitudinal analysis of imaging-based spatiotemporal tumor habitats in glioblastoma, IDH-wild type: implication in patient outcome using multiparametric physiologic MRI

BACKGROUND

Physiologic MRI-based tumor habitat analysis has the potential to predict patient outcomes by identifying the spatiotemporal habitats of glioblastoma. This study aims to prospectively validate the cut-off for tumor progression obtained from tumor habitat analysis based on physiologic MRI in ascertaining time-to-progression (TTP) and the site of progression in glioblastoma patients following concurrent chemoradiotherapy (CCRT).

METHODS

In this prospective study (ClinicalTrials.gov ID: NCT02613988), we will recruit patients with IDH-wild type glioblastoma who underwent CCRT and obtained immediate post-operative and three serial post-CCRT MRI scans within a three-month interval, conducted using diffusion-weighted imaging and dynamic susceptibility contrast imaging. Voxels from cerebral blood volume and apparent diffusion coefficient maps will be grouped using k-means clustering into three spatial habitats (hypervascular cellular, hypovascular cellular, and nonviable tissue). The spatiotemporal habitats of the tumor will be evaluated by comparing changes in each habitat between the serial MRI scans (post-operative and post-CCRT #1, #2, and #3). Associations between spatiotemporal habitats and TTP will be analyzed using cox proportional hazard modeling. The site of progression will be matched with spatiotemporal habitats.

DISCUSSION

The perfusion- and diffusion-derived tumor habitat in glioblastoma is expected to stratify TTP and may serve as an early predictor for tumor progression in patients with IDH wild-type glioblastoma.

TRIAL REGISTRATION

ClinicalTrials.gov ID: NCT02613988.

Histogram analysis of intravoxel incoherent motion imaging: Correlation with molecular prognostic factors and combined subtypes of breast cancer

PURPOSE

To look for links between diffusion and IVIM parameters and different molecular subtypes and prognostic factors through histogram analysis.

MATERIALS AND METHODS

A total of 139 patients with breast cancer who had pre-operative MRI examinations were enrolled in this retrospective study. Histograms of the diffusion and IVIM parameters were analyzed for the whole tumor, and an association was investigated between the parameters and the different molecular prognostic factors and subtypes using the nonparametric test, Spearman's rank correlation, and receiver operating characteristic (ROC) curve.

RESULTS

The histogram metrics of the diffusion and IVIM parameters were significantly different for molecular prognostic factors such as human epidermal receptor factor-2 (HER2), progesterone receptor, estrogen receptor, and ki-67. All histogram metrics displayed a poor correlation with all groups (r = -0.28-0.29). There were significant differences in the histogram metrics for the Luminal B-HER2 (-) vs. HER2-positive (non-luminal) subtypes in the mean and 10th percentile D, with the area under the curves (AUCs) of 0.742 and 0.700, respectively, and for the Luminal A and HER2-positive (non-luminal) subtypes in the 90th percentile and entropy of D*, with AUCs of 0.769 and 0.727, respectively.

CONCLUSION

The histogram metrics of IVIM parameters exhibited links with breast cancer prognosis factors and combined subtypes.

CT Texture Analysis of Perihilar Cholangiocarcinoma-Associations With Tumor Grading, Tumor Markers and Clinical Outcome

BACKGROUND

Texture analysis derived from computed tomography (CT) may provide clinically relevant imaging biomarkers associated with tumor histopathology. Perihilar cholangiocarcinoma is a malignant disease with an overall poor prognosis.

AIMS

The present study sought to elucidate possible associations between texture features derived from CT images with grading, tumor markers, and survival in extrahepatic, perihilar cholangiocarcinomas tumors.

METHODS

This retrospective study included 22 patients (10 females, 45%) with a mean age of 71.8 ± 8.7 years. Texture analysis was performed using the free available Mazda software. All tumors were histopathologically confirmed. Survival and clinical parameters were used as primary study outcomes.

RESULTS

In discrimination analysis, "S(1,1)SumVarnc" was statistically significantly different between patients with long-term survival and nonlong-term survival (mean 275.8 ± 32.6 vs. 239.7 ± 26.0, p = 0.01). The first-order parameter "skewness" was associated with the tumor marker "carcinoembryonic antigen" (CEA) (r = -0.7, p = 0.01). A statistically significant correlation of the texture parameter "S(5,0)SumVarnc" with tumor grading was identified (r = -0.6, p < 0.01). Several other texture features correlated with tumor markers CA-19-9 and AFP, as well as with T and N stage of tumors.

CONCLUSION

Several texture features derived from CT images were associated with tumor characteristics and survival in patients with perihilar cholangiocarcinomas. CT texture features could be used as valuable novel imaging markers in clinical routine.

Super-resolution reconstruction for early cervical cancer magnetic resonance imaging based on deep learning

This study aims to develop a super-resolution (SR) algorithm tailored specifically for enhancing the image quality and resolution of early cervical cancer (CC) magnetic resonance imaging (MRI) images. The proposed method is subjected to both qualitative and quantitative analyses, thoroughly investigating its performance across various upscaling factors and assessing its impact on medical image segmentation tasks. The innovative SR algorithm employed for reconstructing early CC MRI images integrates complex architectures and deep convolutional kernels. Training is conducted on matched pairs of input images through a multi-input model. The research findings highlight the significant advantages of the proposed SR method on two distinct datasets at different upscaling factors. Specifically, at a 2× upscaling factor, the sagittal test set outperforms the state-of-the-art methods in the PSNR index evaluation, second only to the hybrid attention transformer, while the axial test set outperforms the state-of-the-art methods in both PSNR and SSIM index evaluation. At a 4× upscaling factor, both the sagittal test set and the axial test set achieve the best results in the evaluation of PNSR and SSIM indicators. This method not only effectively enhances image quality, but also exhibits superior performance in medical segmentation tasks, thereby providing a more reliable foundation for clinical diagnosis and image analysis.

Impact of radiomics features, pulmonary emphysema score and muscle mass on the rate of pneumothorax and chest tube insertion in CT-guided lung biopsies

Iatrogenic pneumothorax is a relevant complication of computed tomography (CT)-guided percutaneous lung biopsy. The aim of the present study was to analyze the prognostic significance of texture analysis, emphysema score and muscle mass derived from CT-imaging to predict postinterventional pneumothorax after CT-guided lung biopsy. Consecutive patients undergoing CT-guided percutaneous lung biopsy between 2012 and 2021 were analyzed. Multivariate logistic regression analysis included clinical risk factors and CT-imaging features to detect associations with pneumothorax development. Overall, 479 patients (178 females, mean age 65 ± 11.7 years) underwent CT-guided percutaneous lung biopsy of which 180 patients (37.5%) developed pneumothorax including 55 patients (11.5%) requiring chest tube placement. Risk factors associated with pneumothorax were chronic-obstructive pulmonary disease (COPD) (p = 0.03), age (p = 0.02), total lung capacity (p < 0.01) and residual volume (p = 0.01) as well as interventional parameters needle length inside the lung (p < 0.001), target lesion attached to pleura (p = 0.04), and intervention duration (p < 0.001). The combined model demonstrated a prediction accuracy of the occurrence of pneumothorax with an AUC of 0.78 [95%CI: 0.70-0.86] with a resulting sensitivity 0.80 and a specificity of 0.66. In conclusion, radiomics features of the target lesion and the lung lobe CT-emphysema score are predictive for the occurrence of pneumothorax and need for chest insertion after CT-guided lung biopsy.

Whole-liver histogram analysis of hepatocyte-specific contrast-enhanced magnetic resonance imaging for predicting progression in patients with cirrhosis

Background

Liver cirrhosis, as the terminal phase of chronic liver disease fibrosis, is associated with high morbidity and mortality. Traditional methods for assessing liver function, such as clinical scoring systems, offer only a global evaluation and may not accurately reflect regional liver function variations. This study aimed at evaluating the diagnostic potential of whole-liver histogram analysis of gadobenate dimeglumine (Gd-BOPTA)-enhanced magnetic resonance imaging (MRI) for predicting the progression of cirrhosis.

Methods

In this retrospective study, 265 consecutive patients with cirrhosis admitted to the Department of Radiology, Shandong Provincial Hospital Affiliated to Shandong First Medical University from August 2012 to September 2019 were enrolled. After the exclusion criteria were applied, 117 patients (84 males and 33 females) were divided into Child-Pugh A cirrhosis (n=43), Child-Pugh B cirrhosis (n=49), and Child-Pugh C cirrhosis (n=25). After correction for liver signal intensity with the spleen was completed, 19 histogram features of the whole liver were extracted and modeled to evaluate liver function, with the Child-Pugh class being incorporated as a clinical parameter. Receiver operating characteristic (ROC) curves were used to assess the diagnosis capability and determine the optimal cutoffs after a mean follow-up of 42.3±19.1 (range, 8-93) months. The association between significant histogram features and the cumulative incidence of hepatic insufficiency was analyzed with the adjusted Kaplan-Meier curve model.

Results

Among 117 patients (12%), 14 developed hepatic insufficiency through a period of follow-up. Five features, including the median (P<0.01), 90th percentile (P<0.01), root mean squared (P<0.01), mean (P<0.01), and 10th percentile (P<0.05), were significantly different between the groups with and without hepatic insufficiency according to the Kruskal-Wallis test; in the ROC curve analysis, the area under the curve (AUC) of these features was 0.723 [95% confidence interval (CI): 0.653-0.793], 0.722 (95% CI: 0.652-0.792), 0.722 (95% CI: 0.652-0.792), 0.721 (95% CI: 0.651-0.791), and 0.674 (95% CI: 0.600-0.748) after correction, respectively (all P values <0.05). Median, 90th percentile, root mean squared, and mean were found to be significant factors in predicting liver insufficiency. The adjusted Kaplan-Meier curves revealed that patients with a feature level less than the cutoff, as compared to those with a level above the cutoff, showed a statistically shorter progression-free survival and higher incidences of hepatic insufficiency for significant features of median (cutoff =26.001; 21.28% versus 5.71%; P=0.02), 90th percentile (cutoff =86.263; 20.41% versus 5.88%; P<0.01), root mean squared (cutoff =1,028.477; 19.15% versus 7.14%; P=0.049), and mean (cutoff =27.484; 19.15% versus 7.14%; P=0.049). Patients with a 10th percentile less than -39.811 also showed a higher cumulative incidence of hepatic insufficiency than did those with a value higher than the cutoff (0.18% versus 7.46%; P=0.22).

Conclusions

Whole-liver histogram analysis of Gd-BOPTA-enhanced MRI may serve as a noninvasive analytical method to predict hepatic insufficiency in patients with cirrhosis.

Impact of quantitative CT texture analysis on the outcome of CT-guided bone biopsy

Texture analysis can provide new imaging-based biomarkers. Texture analysis derived from computed tomography (CT) might be able to better characterize patients undergoing CT-guided percutaneous bone biopsy. The present study evaluated this and correlated texture features with bioptic outcome in patients undergoing CT-guided bone biopsy. Overall, 123 patients (89 female patients, 72.4 %) were included into the present study. All patients underwent CT-guided percutaneous bone biopsy with an 11 Gauge coaxial needle. Clinical parameters and quantitative imaging features were investigated. Random forest classifier was used to predict a positive biopsy result. Overall, 69 patients had osteolytic metastasis (56.1 %) and 54 had osteoblastic metastasis (43.9 %). The overall positive biopsy rate was 72 %. The developed radiomics model demonstrated a prediction accuracy of a positive biopsy result with an AUC of 0.75 [95 %CI 0.65 - 0.85]. In a subgroup of breast cancer patients, the model achieved an AUC of 0.85 [95 %CI 0.73 - 0.96]. In the subgroup of non-breast cancer patients, the signature achieved an AUC of 0.80 [95 %CI 0.60 - 0.99]. Quantitative CT imaging findings comprised of conventional and texture features can aid to predict the bioptic result of CT-guided bone biopsies. The developed radiomics signature aids in clinical decision-making, and could identify patients at risk for a negative biopsy.

Nomogram of magnetic resonance imaging (MRI) histogram analysis to predict telomerase reverse transcriptase promoter mutation status in glioblastoma

Background

Telomerase reverse transcriptase promoter (pTERT) status is a strong biomarker to diagnose and predict the prognosis of glioblastoma (GBM). In this study, we explored the predictive value of preoperative magnetic resonance imaging (MRI) histogram analysis in the form of nomogram for evaluating pTERT mutation status in GBM.

Methods

The clinical and imaging data of 181 patients with GBM at our hospital between November 2018 and April 2023 were retrospectively assessed. We used the molecular sequencing results to classify the datasets into pTERT mutations (C228T and C250T) and pTERT-wildtype groups. FireVoxel software was used to extract preoperative T1-weighted contrast-enhanced (T1C) histogram parameters of GBM patients. The T1C histogram parameters were compared between groups. Univariate and multivariate logistic regression analyses were used to construct the nomogram, and the predictive efficacy of model was evaluated using calibration and decision curves. Receiver operating characteristic curve was used to assess model performance.

Results

Patient age and percentage of unenhanced tumor area showed statistically significant differences between the pTERT mutation and pTERT-wildtype groups (P<0.001). Among the T1C histogram features, the maximum, standard deviation (SD), variance, coefficient of variation (CV), skewness, 5th, 10th, 25th, 95th and 99th percentiles were statistically significantly different between groups (P=0.000-0.040). Multivariate logistic regression analysis showed that age, percentage of unenhanced tumor area, SD and CV were independent risk factors for predicting pTERT mutation status in GBM patients. The logistic regression model based on these four features showed a better sample predictive performance, and the area under the curve (AUC) [95% confidence interval (CI)], accuracy, sensitivity, specificity were 0.842 (0.767-0.917), 0.796, 0.820, and 0.729, respectively. There were no significant differences in the T1C histogram parameters between the C228T and C250T groups (P=0.055-0.854).

Conclusions

T1C histogram parameters can be used to evaluate pTERT mutations status in GBM. A nomogram based on conventional MRI features and T1C histogram parameters is a reliable tool for the pTERT mutation status, allowing for non-invasive radiological prediction before surgery.

The Predictive Effect of Quantitative Analysis of Signal Intensity Heterogeneity on T2-Weighted MR Images for High-intensity Focused Ultrasound Treatment of Uterine Fibroids

RATIONALE AND OBJECTIVES

To investigate whether the quantitative index of signal intensity (SI) heterogeneity on T2-weighted (T2W) magnetic resonance images can predict the difficulty and efficacy of high-intensity focused ultrasound (HIFU) ablation for uterine fibroids.

MATERIALS AND METHODS

The standard deviation (SD) of T2W image (T2WI) SI was used to quantify SI heterogeneity. The correlation between SD and the non-perfused volume ratio (NPVR) in 575 patients undergoing HIFU treatment was retrospectively analyzed, and the efficacy of SD in predicting NPVR was discussed. Three classifications were made based on the SD, and the ablation difficulty and ablation effect of different grades were compared. A total of 65 cases from another center were used as an external validation set to verify the classification performance of SD.

RESULTS

The SD of SI was negatively correlated with NPVR (r = -0.460, p < 0.001). The predictive efficiency of SD for the ablation effect was higher than that of the scaled signal intensity (0.767 vs. 0.701, p = 0.006). Univariate and multivariate logistic regression analyses showed that SD was an independent predictor of ablation effect. Based on SD, the three classifications were divided into SD I: SD < 101.0, SD II: 101.0 ≤ SD < 138.7, and SD III: SD≥ 138.7. The treatment time, sonication time, treatment intensity, and total energy of SD I were lower than those of SD II and III (p < 0.05).

CONCLUSION

The heterogeneity of T2WI SI of uterine fibroids is negatively correlated with NPVR. The SD of SI can be used to predict the ablation difficulty and ablation effect of HIFU.

Correlation Analysis of Apparent Diffusion Coefficient Histogram Parameters and Clinicopathologic Features for Prognosis Prediction in Uveal Melanoma

Purpose

To investigate the correlation between apparent diffusion coefficient (ADC) histograms and high-risk clinicopathologic features related to uveal melanoma (UM) prognosis.

Methods

This retrospective study included 53 patients with UM who underwent diffusion-weighted imaging (DWI) between August 2015 and March 2024. Axial DWI was performed with a single-shot spin-echo echo-planar imaging sequence. ADC histogram parameters of ADCmean, ADC50%, interquartile range (IQR), skewness, kurtosis, and entropy were obtained from DWI. The relationships between histogram parameters and high-risk clinicopathological characteristics including tumor size, preoperative retinal detachment, histological subtypes, Ki-67 index, and chromosome status, were analyzed by Spearman correlation analysis, Mann-Whitney U test, or Kruskal-Wallis test.

Results

A total of 53 patients (mean ± SD age, 55 ± 15 years; 22 men) were evaluated. The largest basal diameter (LBD) was correlated with kurtosis (r = 0.311, P = 0.024). Tumor prominence (TP) was correlated with entropy (r = 0.581, P < 0.001) and kurtosis (r = 0.273, P = 0.048). Additionally, significant correlations were identified between the Ki-67 index and ADCmean (r = -0.444, P = 0.005), ADC50% (r = -0.487, P = 0.002), and skewness (r = 0.394, P = 0.014). Finally, entropy was correlated with monosomy 3 (r = 0.541, P = 0.017).

Conclusions

The ADC histograms provided valuable insights into high-risk clinicopathologic features of UM and hold promise in the early prediction of UM prognosis.

A Combination of Amide Proton Transfer, Tumor Blood Flow, and Apparent Diffusion Coefficient Histogram Analysis Is Useful for Differentiating Malignant from Benign Intracranial Tumors in Young Patients: A Preliminary Study

PURPOSE

To evaluate the amide proton transfer (APT), tumor blood flow (TBF), and apparent diffusion coefficient (ADC) combined diagnostic value for differentiating intracranial malignant tumors (MTs) from benign tumors (BTs) in young patients, as defined by the 2021 World Health Organization classification of central nervous system tumors.

METHODS

Fifteen patients with intracranial MTs and 10 patients with BTs aged 0-30 years underwent MRI with APT, pseudocontinuous arterial spin labeling (pCASL), and diffusion-weighted imaging. All tumors were evaluated through the use of histogram analysis and the Mann-Whitney U test to compare 10 parameters for each sequence between the groups. The diagnostic performance was evaluated using receiver operating characteristic (ROC) curve analysis.

RESULTS

The APT maximum, mean, 10th, 25th, 50th, 75th, and 90th percentiles were significantly higher in MTs than in BTs; the TBF minimum (min) was significantly lower in MTs than in BTs; TBF kurtosis was significantly higher in MTs than in BTs; the ADC min, 10th, and 25th percentiles were significantly lower in MTs than in BTs (all p < 0.05). The APT 50th percentile (0.900), TBF min (0.813), and ADC min (0.900) had the highest area under the curve (AUC) values of the parameters in each sequence. The AUC for the combination of these three parameters was 0.933.

CONCLUSIONS

The combination of APT, TBF, and ADC evaluated through histogram analysis may be useful for differentiating intracranial MTs from BTs in young patients.

MRI for Differentiation between HPV-Positive and HPV-Negative Oropharyngeal Squamous Cell Carcinoma: A Systematic Review

Human papillomavirus (HPV) is an important risk factor for oropharyngeal squamous cell carcinoma (OPSCC). HPV-positive (HPV+) cases are associated with a different pathophysiology, microstructure, and prognosis compared to HPV-negative (HPV-) cases. This review aimed to investigate the potential of magnetic resonance imaging (MRI) to discriminate between HPV+ and HPV- tumours and predict HPV status in OPSCC patients. A systematic literature search was performed on 15 December 2022 on EMBASE, MEDLINE ALL, Web of Science, and Cochrane according to PRISMA guidelines. Twenty-eight studies (n = 2634 patients) were included. Five, nineteen, and seven studies investigated structural MRI (e.g., T1, T2-weighted), diffusion-weighted MRI, and other sequences, respectively. Three out of four studies found that HPV+ tumours were significantly smaller in size, and their lymph node metastases were more cystic in structure than HPV- ones. Eleven out of thirteen studies found that the mean apparent diffusion coefficient was significantly higher in HPV- than HPV+ primary tumours. Other sequences need further investigation. Fourteen studies used MRI to predict HPV status using clinical, radiological, and radiomics features. The reported areas under the curve (AUC) values ranged between 0.697 and 0.944. MRI can potentially be used to find differences between HPV+ and HPV- OPSCC patients and predict HPV status with reasonable accuracy. Larger studies with external model validation using independent datasets are needed before clinical implementation.

radioGWAS links radiome to genome to discover driver genes with somatic mutations for heterogeneous tumor image phenotype in pancreatic cancer

Addressing the significant level of variability exhibited by pancreatic cancer necessitates the adoption of a systems biology approach that integrates molecular data, biological properties of the tumors, medical images, and clinical features of the patients. In this study, a comprehensive multi-omics methodology was employed to examine a distinctive collection of patient dataset containing rapid autopsy tumor and normal tissue samples as well as longitudinal imaging with a focus on pancreatic cancer. By performing a whole exome sequencing analysis on tumor and normal tissues to identify somatic gene variants and a radiomic feature analysis to tumor CT images, the genome-wide association approach established a connection between pancreatic cancer driver genes and relevant radiomic features, enabling a thorough and quantitative assessment of the heterogeneity of pancreatic tumors. The significant association between sets of genes and radiomic features revealed the involvement of genes in shaping tumor morphological heterogeneity. Some results of the association established a connection between the molecular level mechanism and their outcomes at the level of tumor structural heterogeneity. Because tumor structure and tumor structural heterogeneity are related to the patients' overall survival, patients who had pancreatic cancer driver gene mutations with an association to a certain radiomic feature have been observed to experience worse survival rates than cases without these somatic mutations. Furthermore, the association analysis has revealed potential gene mutations and radiomic feature candidates that warrant further investigation in future research endeavors.

Radiomics analysis for prediction of lymph node metastasis after neoadjuvant chemotherapy based on pretreatment MRI in patients with locally advanced cervical cancer

Background

This study aims to develop and validate a pretreatment MRI-based radiomics model to predict lymph node metastasis (LNM) following neoadjuvant chemotherapy (NACT) in patients with locally advanced cervical cancer (LACC).

Methods

Patients with LACC who underwent NACT from two centers between 2013 and 2022 were enrolled retrospectively. Based on the lymph node (LN) status determined in the pathology reports after radical hysterectomy, patients were categorized as LN positive or negative. The patients from center 1 were assigned as the training set while those from center 2 formed the validation set. Radiomics features were extracted from pretreatment sagittal T2-weighted imaging (Sag-T2WI), axial diffusion-weighted imaging (Ax-DWI), and the delayed phase of dynamic contrast-enhanced sagittal T1-weighted imaging (Sag-T1C) for each patient. The K-best and least absolute shrinkage and selection operator (LASSO) methods were employed to reduce dimensionality, and the radiomics features strongly associated with LNM were selected and used to construct three single-sequence models. Furthermore, clinical variables were incorporated through multivariate regression analysis and fused with the selected radiomics features to construct the clinical-radiomics combined model. The diagnostic performance of the models was assessed using receiver operating characteristic (ROC) curve analysis. The clinical utility of the models was evaluated by the area under the ROC curve (AUC) and decision curve analysis (DCA).

Results

A total of 282 patients were included, comprising 171 patients in the training set, and 111 patients in the validation set. Compared to the Sag-T2WI model (AUC, 95%CI, training set, 0.797, 0.722-0.782; validation set, 0.648, 0.521-0.776) and the Sag-T1C model (AUC, 95%CI, training set, 0.802, 0.723-0.882; validation set, 0.630, 0.505-0.756), the Ax-DWI model exhibited the highest diagnostic performance with AUCs of 0.855 (95%CI, 0.791-0.919) in training set, and 0.753 (95%CI, 0.638-0.867) in validation set, respectively. The combined model, integrating selected features from three sequences and FIGO stage, surpassed predictive ability compared to the single-sequence models, with AUC of 0.889 (95%CI, 0.833-0.945) and 0.859 (95%CI, 0.781-0.936) in the training and validation sets, respectively.

Conclusions

The pretreatment MRI-based radiomics model, integrating radiomics features from three sequences and clinical variables, exhibited superior performance in predicting LNM following NACT in patients with LACC.

Novel 3D morphological characteristics for congenital biliary dilatation diagnosis: a case-control study

BACKGROUND

Congenital biliary dilatation (CBD) necessitates the timely removal of dilated bile ducts. Accurate differentiation between CBD and secondary biliary dilatation (SBD) is crucial for treatment decisions, and identification of CBD with intrahepatic involvement is vital for surgical planning and supportive care. This study aimed to develop quantitative models based on bile duct morphology to distinguish CBD from SBD and further identify CBD with intrahepatic involvement.

MATERIALS AND METHODS

The retrospective study included 131 CBD and 209 SBD patients between December 2014 and December 2021 for model development, internal validation, and testing. A separate cohort of 15 CBD and 34 SBD patients between January 2022 and December 2022 was recruited for temporally-independent validation. Quantitative shape-based (Shape) and diameter-based (Diam) morphological characteristics of bile ducts were extracted to build a CBD diagnosis model to distinguish CBD from SBD and an intrahepatic involvement identification model to classify CBD with/without intrahepatic involvement. The diagnostic performance of the models was compared with that of experienced hepatobiliary surgeons.

RESULTS

The CBD diagnosis model using clinical, Shape, and Diam characteristics showed good performance with an AUROC of 0.942 (95% CI: 0.890-0.994), AUPRC of 0.917 (0.855-0.979), accuracy of 0.891, sensitivity of 0.950, and F1-score of 0.864. The model outperformed two experienced surgeons in accuracy, sensitivity, and F1-score. The intrahepatic involvement identification model using clinical, Shape, and Diam characteristics yielded outstanding performance with an AUROC of 0.944 (0.879-1.000), AUPRC of 0.982 (0.947-1.000), accuracy of 0.932, sensitivity of 0.971, and F1-score of 0.957. The models demonstrated generalizable performance on the temporally-independent validation cohort.

CONCLUSIONS

This study developed two robust quantitative models for distinguishing CBD from SBD and identifying CBD with intrahepatic involvement, respectively, based on morphological characteristics of the bile ducts, showing great potential in risk stratification and surgical planning of CBD.

Precision Phenotypic Profiling and Capture of Circulating Tumor Cells via a Vertical Laminar Flow-Stacked Microfluidic Chip

The heterogeneity of circulating tumor cells has a significant impact on the diagnosis, treatment, and monitoring of cancer. Research on the subtypes of circulating tumor cells can bring better treatment outcomes for cancer patients. Here, we proposed a microfluidic chip for the magnetic capture of subtypes of circulating tumor cells from the whole blood and phenotypic profiling by stacking laminar flow vertically. Circulating tumor cells were sorted and captured by the three-dimensional regulation of both magnetic fields in the vertical direction and flow fields in the lateral direction. Using EpCAM-magnetic beads, we achieved sorting and sectional capture of target cells in whole blood and analyzed the surface expression levels of the captured cells, confirming the functionality of the microfluidic chip in sorting and capturing subtypes of circulating tumor cells. This microfluidic chip can also aid in the subsequent subtype analysis of other rare cells.

Spatial heterogeneity of peri-tumoural lipid composition in postmenopausal patients with oestrogen receptor positive breast cancer

Deregulation of lipid composition in adipose tissue adjacent to breast tumour is observed in ex vivo and animal models. Novel non-invasive magnetic resonance imaging (MRI) allows rapid lipid mapping of the human whole breast. We set out to elucidate the spatial heterogeneity of peri-tumoural lipid composition in postmenopausal patients with oestrogen receptor positive (ER +) breast cancer. Thirteen participants (mean age, 62 ± [SD] 6 years) with ER + breast cancer and 13 age-matched postmenopausal healthy controls were scanned on MRI. The number of double bonds in triglycerides was computed from MRI images to derive lipid composition maps of monounsaturated, polyunsaturated, and saturated fatty acids (MUFA, PUFA, SFA). The spatial heterogeneity measures (mean, median, skewness, entropy and kurtosis) of lipid composition in the peri-tumoural region and the whole breast of participants and in the whole breast of controls were computed. The Ki-67 proliferative activity marker and CD163 antibody on tumour-associated macrophages were assessed histologically. Mann Whitney U or Wilcoxon tests and Spearman's coefficients were used to assess group differences and correlations, respectively. For comparison against the whole breast in participants, peri-tumoural MUFA had a lower mean (median (IQR), 0.40 (0.02), p < .001), lower median (0.42 (0.02), p < .001), a negative skewness with lower magnitude (- 1.65 (0.77), p = .001), higher entropy (4.35 (0.64), p = .007) and lower kurtosis (5.13 (3.99), p = .001). Peri-tumoural PUFA had a lower mean (p < .001), lower median (p < .001), a positive skewness with higher magnitude (p = .005) and lower entropy (p = .002). Peri-tumoural SFA had a higher mean (p < .001), higher median (p < .001), a positive skewness with lower magnitude (p < .001) and lower entropy (p = .012). For comparison against the whole breast in controls, peri-tumoural MUFA had a negative skewness with lower magnitude (p = .01) and lower kurtosis (p = .009), however there was no difference in PUFA or SFA. CD163 moderately correlated with peri-tumoural MUFA skewness (rs = - .64), PUFA entropy (rs = .63) and SFA skewness (rs = .59). There was a lower MUFA and PUFA while a higher SFA, and a higher heterogeneity of MUFA while a lower heterogeneity of PUFA and SFA, in the peri-tumoural region in comparison with the whole breast tissue. The degree of lipid deregulation was associated with inflammation as indicated by CD163 antibody on macrophages, serving as potential marker for early diagnosis and response to therapy.

Radiomics analysis of R2* maps to predict early recurrence of single hepatocellular carcinoma after hepatectomy

Objectives

This study aimed to evaluate the effectiveness of radiomics analysis with R2* maps in predicting early recurrence (ER) in single hepatocellular carcinoma (HCC) following partial hepatectomy.

Methods

We conducted a retrospective analysis involving 202 patients with surgically confirmed single HCC having undergone preoperative magnetic resonance imaging between 2018 and 2021 at two different institutions. 126 patients from Institution 1 were assigned to the training set, and 76 patients from Institution 2 were assigned to the validation set. A least absolute shrinkage and selection operator (LASSO) regularization was conducted to operate a logistic regression, then features were identified to construct a radiomic score (Rad-score). Uni- and multi-variable tests were used to assess the correlations of clinicopathological features and Rad-score with ER. We then established a combined model encompassing the optimal Rad-score and clinical-pathological risk factors. Additionally, we formulated and validated a predictive nomogram for predicting ER in HCC. The nomogram's discrimination, calibration, and clinical utility were thoroughly evaluated.

Results

Multivariable logistic regression revealed the Rad-score, microvascular invasion (MVI), and α fetoprotein (AFP) level > 400 ng/mL as significant independent predictors of ER in HCC. We constructed a nomogram based on these significant factors. The areas under the receiver operator characteristic curve of the nomogram and precision-recall curve were 0.901 and 0.753, respectively, with an F1 score of 0.831 in the training set. These values in the validation set were 0.827, 0.659, and 0.808.

Conclusion

The nomogram that integrates the radiomic score, MVI, and AFP demonstrates high predictive efficacy for estimating the risk of ER in HCC. It facilitates personalized risk classification and therapeutic decision-making for HCC patients.

Diffusion histogram profiles predict molecular features of grade 4 in histologically lower-grade adult diffuse gliomas following WHO classification 2021

OBJECTIVES

In the latest World Health Organization classification 2021, grade 4 adult diffuse gliomas can be diagnosed with several molecular features even without histological evidence of necrosis or microvascular proliferation. We aimed to explore whole tumor histogram-derived apparent diffusion coefficient (ADC) histogram profiles for differentiating between the presence (Mol-4) and absence (Mol-2/3) of grade 4 molecular features in histologically lower-grade gliomas.

METHODS

Between June 2019 and October 2022, 184 adult patients with diffuse gliomas underwent MRI. After excluding 121 patients, 18 (median age, 64.5 [range, 37-84 years]) Mol-4 and 45 (median 40 [range, 18-73] years) Mol-2/3 patients with histologically lower-grade gliomas were enrolled. Whole tumor volume-of-interest-derived ADC histogram profiles were calculated and compared between the two groups. Stepwise logistic regression analysis with Akaike's information criterion using the ADC histogram profiles with p values < 0.01 and age at diagnosis was used to identify independent variables for predicting the Mol-4 group.

RESULTS

The 90th percentile (p < 0.001), median (p < 0.001), mean (p < 0.001), 10th percentile (p = 0.014), and entropy (p < 0.001) of normalized ADC were lower, and kurtosis (p < 0.001) and skewness (p = 0.046) were higher in the Mol-4 group than in the Mol-2/3 group. Multivariate logistic regression analysis revealed that the entropy of normalized ADC and age at diagnosis were independent predictive parameters for the Mol-4 group with an area under the curve of 0.92.

CONCLUSION

ADC histogram profiles may be promising preoperative imaging biomarkers to predict molecular grade 4 among histologically lower-grade adult diffuse gliomas.

CLINICAL RELEVANCE STATEMENT

This study highlighted the diagnostic usefulness of ADC histogram profiles to differentiate histologically lower grade adult diffuse gliomas with the presence of molecular grade 4 features and those without.

KEY POINTS

• ADC histogram profiles to predict molecular CNS WHO grade 4 status among histologically lower-grade adult diffuse gliomas were evaluated. • Entropy of ADC and age were independent predictive parameters for molecular grade 4 status. • ADC histogram analysis is useful for predicting molecular grade 4 among histologically lower-grade gliomas.

Preoperative CT histogram analysis to predict the expression of Ki-67 in solid pseudopapillary tumours of the pancreas

AIM

To explore the value of preoperative computed tomography (CT) histogram features in predicting the expression status of Ki-67 in patients with solid pseudopapillary pancreatic tumours (SPTP).

MATERIALS AND METHODS

This retrospective study analysed venous phase CT images of 39 patients with SPTP confirmed at surgery and histopathology and measured using the Ki-67 proliferation index from November 2015 to February 2022. According to the Ki-67 proliferation index, they were divided into high expression (Ki-67 ≥ 4%) and low expression (Ki-67 < 4%) groups. The histogram features of quantitative parameters were extracted using MaZda software, and the quantitative parameters of CT histograms were compared between groups. The receiver operating characteristic (ROC) curves of the patients were plotted according to the parameters, with statistically significant differences. The area under the curve (AUC), sensitivity, and specificity were calculated, and the effectiveness of the histogram parameters in predicting Ki-67 expression was analysed and evaluated.

RESULTS

In total, 27 SPTP patients were enrolled, including 11 with high expression of Ki-67 and 16 with low expression. Comparative analysis of the Ki-67 high- and low-expression groups revealed a statistically significant in necrosis and variance (p<0.05). ROC curve analysis showed that the AUC of necrosis and variance predicting Ki-67 expression status were 0.753 and 0.841, the sensitivities were 81.8% and 81.3%, and the specificities were 68.7% and 81.8%, respectively.

CONCLUSION

Preoperative CT histogram features help predict Ki-67 expression status in patients with SPTP.

Quantitative assessment of inflammation and evaluation of spatial heterogeneity for non-alcoholic fatty liver disease in mice based on iron-adjustive T1

Background

A sensitive and non-invasive method is necessary to diagnose non-alcoholic fatty liver disease (NAFLD). We explored the iron-adjustive T1 (aT1) ability to quantify the degree of liver inflammation and evaluate the spatial heterogeneity.

Methods

Male C57BL/6J mice were randomly categorized as the NAFLD model (n=40), NAFLD-related liver cirrhosis model (n=20), and normal mice (n=10). T1 and T2* maps were acquired using a 3.0T scanner of magnetic resonance imaging (MRI) and aT1 maps through post-processing corrected iron's effect on T1 using T2*. Pathological changes in the left and right liver lobes were assessed using the Non-alcoholic Steatohepatitis-Clinical Research Network scoring system, though hepatic ballooning lesion were rare in models. Spearman's and partial correlation analyses were used to evaluate correlations, and the receiver operating characteristic curve was used to analyze the diagnostic performance.

Results

aT1 was highly correlated with NAFLD activity score (NAS) (r=0.747, P<0.001) but not with the fibrosis stage when adjusted by NAS (r=-0.135, P=0.147). The area under the curve (AUC) of the aT1 value distinguishing groups with 0< NAS <4 and NAS ≥4 was 0.802. On analyzing the histogram features of aT1, the entropy, interquartile range, range, and variance were significantly different between the groups with 0< NAS <4 and NAS ≥4 (P<0.05). The entropy was the risk factor of NAS ≥4.

Conclusions

aT1 could help evaluate the inflammatory activity in NAFLD mice unaffected by mild fibrosis, and the higher the degree of inflammation, the higher the heterogeneity of the aT1 map.

A computational approach for analysis of intratumoral heterogeneity and standardized uptake value in PET/CT images1

BACKGROUND

By providing both functional and anatomical information from a single scan, digital imaging technologies like PET/CT and PET/MRI hybrids are gaining popularity in medical imaging industry. In clinical practice, the median value (SUVmed) receives less attention owing to disagreements surrounding what defines a lesion, but the SUVmax value, which is a semi-quantitative statistic used to analyse PET and PET/CT images, is commonly used to evaluate lesions.

OBJECTIVE

This study aims to build an image processing technique with the purpose of automatically detecting and isolating lesions in PET/CT images, as well as measuring and assessing the SUVmed.

METHODS

The pictures are separated into their respective lesions using mathematical morphology and the crescent region, which are both part of the image processing method. In this research, a total of 18 different pictures of lesions were evaluated.

RESULTS

The findings of the study reveal that the threshold is satisfied by both the SUVmax and the SUVmed for most of the lesion types. However, in six instances, the SUVmax and SUVmed values are found to be in different courts.

CONCLUSION

The new information revealed by this study needs to be further investigated to determine if it has any practical value in diagnosing and monitoring lesions. However, results of this study suggest that SUVmed should receive more attention in the evaluation of lesions in PET and CT images.

Histogram analysis of mono-exponential, bi-exponential and stretched-exponential diffusion-weighted MR imaging in predicting consistency of meningiomas

BACKGROUND

The consistency of meningiomas is critical to determine surgical planning and has a significant impact on surgical outcomes. Our aim was to compare mono-exponential, bi-exponential and stretched exponential MR diffusion-weighted imaging in predicting the consistency of meningiomas before surgery.

METHODS

Forty-seven consecutive patients with pathologically confirmed meningiomas were prospectively enrolled in this study. Two senior neurosurgeons independently evaluated tumour consistency and classified them into soft and hard groups. A volume of interest was placed on the preoperative MR diffusion images to outline the whole tumour area. Histogram parameters (mean, median, 10th percentile, 90th percentile, kurtosis, skewness) were extracted from 6 different diffusion maps including ADC (DWI), D*, D, f (IVIM), alpha and DDC (SEM). Comparisons between two groups were made using Student's t-Test or Mann-Whitney U test. Parameters with significant differences between the two groups were included for Receiver operating characteristic analysis. The DeLong test was used to compare AUCs.

RESULTS

DDC, D* and ADC 10th percentile were significantly lower in hard tumours than in soft tumours (P ≤ 0.05). The alpha 90th percentile was significantly higher in hard tumours than in soft tumours (P < 0.02). For all histogram parameters, the alpha 90th percentile yielded the highest AUC of 0.88, with an accuracy of 85.10%. The D* 10th percentile had a relatively higher AUC value, followed by the DDC and ADC 10th percentile. The alpha 90th percentile had a significantly greater AUC value than the ADC 10th percentile (P ≤ 0.05). The D* 10th percentile had a significantly greater AUC value than the ADC 10th percentile and DDC 10th percentile (P ≤ 0.03).

CONCLUSION

Histogram parameters of Alpha and D* may serve as better imaging biomarkers to aid in predicting the consistency of meningioma.

Presurgical MRI-Based Radiomics Models for Predicting Cerebellar Mutism Syndrome in Children With Posterior Fossa Tumors

BACKGROUND

Current studies have indicated that tumoral morphologic features are associated with cerebellar mutism syndrome (CMS), but the radiomics application in CMS is scarce.

PURPOSE

To develop a model for CMS discrimination based on multiparametric MRI radiomics in patients with posterior fossa tumors.

STUDY TYPE

Retrospective.

POPULATION

A total of 218 patients (males 132, females 86) with posterior fossa tumors, 169 of which were included in the MRI radiomics analysis. The MRI radiomics study cohort (169) was split into training (119) and testing (50) sets with a ratio of 7:3.

FIELD/SEQUENCE

All the MRI were acquired under 1.5/3.0 T scanners. T2-weighted image (T2W), T1-weighted (T1W), fluid attenuated inversion recovery (FLAIR), diffusion-weighted imaging (DWI).

ASSESSMENT

Apparent diffusion coefficient (ADC) maps were generated from DWI. Each MRI dataset generated 1561 radiomics characteristics. Feature selection was performed with univariable logistic analysis, correlation analysis, and least absolute shrinkage and selection operator (LASSO) penalized logistic regression. Significant clinical features were selected with multivariable logistic analysis and used to constructed the clinical model. Radiomics models (based on T1W, T2W, FLAIR, DWI, ADC) were constructed with selected radiomics features. The mix model was based on the multiparametric MRI radiomics features.

STATISTICAL TEST

Multivariable logistic analysis was utilized during clinical features selection. Models' performance was evaluated using the area under the receiver operating characteristic (AUC) curve. Interobserver variability was assessed using Cohen's kappa. Significant threshold was set as P < 0.05.

RESULTS

Sex (aOR = 3.72), tumor location (aOR = 2.81), hydrocephalus (aOR = 2.14), and tumor texture (aOR = 5.08) were significant features in the multivariable analysis and were used to construct the clinical model (AUC = 0.79); totally, 33 radiomics features were selected to construct radiomics models (AUC = 0.63-0.93). Seven of the 33 radiomics features were selected for the mix model (AUC = 0.93).

DATA CONCLUSION

Multiparametric MRI radiomics may be better at predicting CMS than single-parameter MRI models and clinical model.

EVIDENCE LEVEL

4.

TECHNICAL EFFICACY

2.

MRI radiomics model for predicting TERT promoter mutation status in glioblastoma

BACKGROUND AND PURPOSE

The presence of TERT promoter mutations has been associated with worse prognosis and resistance to therapy for patients with glioblastoma (GBM). This study aimed to determine whether the combination model of different feature selections and classification algorithms based on multiparameter MRI can be used to predict TERT subtype in GBM patients.

METHODS

A total of 143 patients were included in our retrospective study, and 2553 features were obtained. The datasets were randomly divided into training and test sets in a ratio of 7:3. The synthetic minority oversampling technique was used to achieve data balance. The Pearson correlation coefficients were used for dimension reduction. Three feature selections and five classification algorithms were used to model the selected features. Finally, 10-fold cross validation was applied to the training dataset.

RESULTS

A model with eight features generated by recursive feature elimination (RFE) and linear discriminant analysis (LDA) showed the greatest diagnostic performance (area under the curve values for the training, validation, and testing sets: 0.983, 0.964, and 0.926, respectively), followed by relief and random forest (RF), analysis of variance and RF. Furthermore, the relief was the optimal feature selection for separately evaluating those five classification algorithms, and RF was the most preferable algorithm for separately assessing the three feature selectors. ADC entropy was the parameter that made the greatest contribution to the discrimination of TERT mutations.

CONCLUSIONS

Radiomics model generated by RFE and LDA mainly based on ADC entropy showed good performance in predicting TERT promoter mutations in GBM.

Tumor radiogenomics in gliomas with Bayesian layered variable selection

We propose a statistical framework to analyze radiological magnetic resonance imaging (MRI) and genomic data to identify the underlying radiogenomic associations in lower grade gliomas (LGG). We devise a novel imaging phenotype by dividing the tumor region into concentric spherical layers that mimics the tumor evolution process. MRI data within each layer is represented by voxel-intensity-based probability density functions which capture the complete information about tumor heterogeneity. Under a Riemannian-geometric framework these densities are mapped to a vector of principal component scores which act as imaging phenotypes. Subsequently, we build Bayesian variable selection models for each layer with the imaging phenotypes as the response and the genomic markers as predictors. Our novel hierarchical prior formulation incorporates the interior-to-exterior structure of the layers, and the correlation between the genomic markers. We employ a computationally-efficient Expectation-Maximization-based strategy for estimation. Simulation studies demonstrate the superior performance of our approach compared to other approaches. With a focus on the cancer driver genes in LGG, we discuss some biologically relevant findings. Genes implicated with survival and oncogenesis are identified as being associated with the spherical layers, which could potentially serve as early-stage diagnostic markers for disease monitoring, prior to routine invasive approaches. We provide a R package that can be used to deploy our framework to identify radiogenomic associations.

Histogram analysis of quantitative susceptibility mapping and apparent diffusion coefficient for identifying isocitrate dehydrogenase genotypes and tumor subtypes of adult-type diffuse gliomas

Background

Accurate preoperative identification of isocitrate dehydrogenase (IDH) genotypes and tumor subtypes is highly important for proper treatment planning and prognosis evaluation in patients with glioma. This study aimed to differentiate IDH genotypes and tumor subtypes of adult-type diffuse gliomas using histogram features of quantitative susceptibility mapping (QSM) and apparent diffusion coefficient (ADC).

Methods

This prospective study enrolled patients with suspected gliomas between March 2019 and January 2022 in a random series. Histogram features of QSM and ADC were extracted from the tumor parenchyma. The Mann-Whitney U test was used to compare the difference in histogram features between different IDH genotypes and among tumor subtypes. Receiver operating characteristic (ROC) curves were constructed to assess the corresponding diagnostic performance.

Results

This study included 47 patients with histopathologically confirmed adult-type diffuse gliomas. Totals of seven QSM features including 10th percentile (P10), 90th percentile (P90), interquartile range (IQR), maximum, mean absolute deviation (MAD), root mean squared (RMS), and variance, and five ADC features including P10, mean, median, RMS, and skewness exhibited significant differences between different IDH genotypes (P<0.05 for all), with the IQR of QSM demonstrating the highest area under curve (AUC) of 0.774 [95% confidence interval (CI): 0.635-0.913]. For separating tumor subtypes, the IQR of QSM also showed the highest AUC of 0.745 (95% CI: 0.566-0.924) for glioblastoma (GBM) versus astrocytoma and 0.848 (95% CI: 0.706-0.989) for GBM versus oligodendroglioma, but none of the features could discriminate astrocytoma from oligodendroglioma. The combination of the IQR of QSM, P10 of ADC, and age achieved the highest AUC of 0.910 (95% CI: 0.826-0.994) for IDH genotypes, and 0.939 (95% CI: 0.859-1.000) and 0.967 (95% CI: 0.904-1.000) for GBM versus astrocytoma and GBM versus oligodendroglioma, respectively.

Conclusions

QSM and ADC histogram features may serve as potential imaging markers for noninvasively assessing IDH genotypes and tumor subtypes of adult-type diffuse gliomas. Combining significant features may enhance the diagnostic performance substantially.

The statistical challenge of analysing changes in dual energy computed tomography (DECT) urate volumes in people with gout

BACKGROUND

Dual energy computed tomography (DECT) allows direct visualization of monosodium urate crystal deposition in gout. However, DECT urate volume data are often highly skewed (mostly small volumes with the remainder considerably larger), making statistical analyses challenging in longitudinal research. The aim of this study was to explore the ability of various analysis methods to normalise DECT urate volume data and determine change in DECT urate volumes over time.

METHODS

Simulated datasets containing baseline and year 1 DECT urate volumes for 100 people with gout were created from two randomised controlled trials. Five methods were used to transform the DECT urate volume data prior to analysis: log-transformation, Box-Cox transformation, log(X-(min(X)-1)) transformation; inverse hyperbolic sine transformation, and rank order. Linear regression analyses were undertaken to determine the change in DECT urate volume between baseline and year 1. Cohen's d were calculated as a measure of effect size for each data treatment method. These analyses were then tested in a validation clinical trial dataset containing baseline and year 1 DECT urate volumes from 91 people with gout.

RESULTS

No data treatment method successfully normalised the distribution of DECT urate volumes. For both simulated and validation data sets, significant reductions in DECT urate volumes were observed between baseline and Year 1 across all data treatment methods and there were no significant differences in Cohen's d effect sizes.

CONCLUSIONS

Normalising highly skewed DECT urate volume data is challenging. Adopting commonly used transformation techniques may not significantly improve the ability to determine differences in measures of central tendency when comparing the change in DECT urate volumes over time.

Differentiation of Hepatocellular Carcinoma from Intrahepatic Cholangiocarcinoma through MRI Radiomics

The purpose of this study was to investigate the efficacy of magnetic resonance imaging (MRI) radiomics in differentiating hepatocellular carcinoma (HCC) from intrahepatic cholangiocarcinoma (ICC). The clinical and MRI data of 129 pathologically confirmed HCC patients and 48 ICC patients treated at the Affiliated Hospital of North Sichuan Medical College between April 2016 and December 2021 were retrospectively analyzed. The patients were randomly divided at a ratio of 7:3 into a training group of 124 patients (90 with HCC and 34 with ICC) and a validation group of 53 patients (39 with HCC and 14 with ICC). Radiomic features were extracted from axial fat suppression T2-weighted imaging (FS-T2WI) and axial arterial-phase (AP) and portal-venous-phase (PVP) dynamic-contrast-enhanced MRI (DCE-MRI) sequences, and the corresponding datasets were generated. The least absolute shrinkage and selection operator (LASSO) method was used to select the best radiomic features. Logistic regression was used to establish radiomic models for each sequence (FS-T2WI, AP and PVP models), a clinical model for optimal clinical variables (C model) and a joint radiomics model (JR model) integrating the radiomics features of all the sequences as well as a radiomics-clinical model combining optimal radiomic features and clinical risk factors (RC model). The performance of each model was evaluated using the area under the receiver operating characteristic curve (AUC). The AUCs of the FS-T2WI, AP, PVP, JR, C and RC models for distinguishing HCC from ICC were 0.693, 0.863, 0.818, 0.914, 0.936 and 0.977 in the training group and 0.690, 0.784, 0.727, 0.802, 0.860 and 0.877 in the validation group, respectively. The results of this study suggest that MRI-based radiomics may help noninvasively differentiate HCC from ICC. The model integrating the radiomics features and clinical risk factors showed a further improvement in performance.

radioGWAS: link radiome to genome to discover driver genes with somatic mutations for heterogeneous tumor image phenotype in pancreatic cancer

Addressing the significant level of variability exhibited by pancreatic cancer necessitates the adoption of a systems biology approach that integrates molecular data, biological properties of the tumors, and clinical features of the patients. In this study, a comprehensive multi-omics methodology was employed to examine a distinctive collection patient dataset containing rapid autopsy tumor and normal tissue samples as well as longitudinal imaging with a focus on pancreatic cancer. By performing a whole exome sequencing analysis on tumor and normal tissues to identify somatic gene variants and a radiomics feature analysis to tumor CT images, the genome-wide association approach established a connection between pancreatic cancer driver genes and relevant radiomics features, enabling a thorough and quantitative assessment of the heterogeneity of pancreatic tumors. The significant association between sets of genes and radiomics features revealed the involvement of genes in shaping tumor morphological heterogeneity. Some results of the association established a connection between the molecular level mechanism and their outcomes at the level of tumor structural heterogeneity. Because tumor structure and tumor structural heterogeneity are related to the patients' overall survival, patients who had pancreatic cancer driver gene mutations with an association to a certain radiomics feature have been observed to experience worse survival rates than cases without these somatic mutations. Furthermore, the outcome of the association analysis has revealed potential gene mutations and radiomics feature candidates that warrant further investigation in future research endeavors.

Characterization of Microvascular Invasion in Hepatocellular Carcinoma Using Computational Modeling of Interstitial Fluid Pressure and Velocity

BACKGROUND

Most solid tumors show increased interstitial fluid pressure (IFP), and this increased IFP is an obstacle to treatment. A noninvasive model for measuring IFP in hepatocellular carcinoma (HCC) is an unresolved issue.

PURPOSE

To develop a noninvasive model to measure IFP and interstitial fluid velocity (IFV) in HCC and to characterize the microvascular invasion (MVI) status by using this model.

STUDY TYPE

Retrospective.

POPULATION

A total of 97 HCC patients (mean age 57.6 ± 10.9 years, 77.3% males), 53 of them with MVI and 44 of them without MVI.

FIELD STRENGTH/SEQUENCE

A 3-T, three-dimensional spoiled gradient-recalled echo.

ASSESSMENT

MVI was defined as microscopic vascular invasion of small vessels within the peritumoral liver tissue. The volumes of interest (VOIs) were manually delineated and enclosed the tumor lesion and healthy liver parenchyma, respectively. The extended Tofts model (ETM) was used to estimate permeability parameters from all the VOIs. Subsequently, the continuity partial differential equation (PDE) was implemented and IFP and IFV were acquired.

STATISTICAL TESTS

Wilcoxon signed-ranks tests, histogram analysis, Mann-Whitney U test, Fisher's exact test, least absolute shrinkage and selection operator (LASSO) logistic regression, receiver operating characteristic (ROC) curve analysis with the area under the curve (AUC), Youden index, DeLong test, and Benjamini-Hochberg correction. A P value <0.05 was considered statistically significant.

RESULTS

The HCC lesions exhibited elevated IFP and reduced IFV. There were no significant differences in any measured demographic and clinical features between the MVI-positive and MVI-negative groups, except for tumor size. Nine IFP histogram analysis-derived parameters and seven IFV histogram analysis-derived parameters could be used to characterize the MVI status. LASSO regression selected five features: IFP maximum, IFP 10th percentile, IFP 90th percentile, IFV SD, and IFV 10th percentile. The combination of these features showed the highest AUC (0.781) and specificity (77.3%).

DATA CONCLUSION

A noninvasive IFP and IFV measurement model for HCC was developed. Specific IFP- and IFV-derived parameters exhibited significant association with the MVI status.

EVIDENCE LEVEL

3.

TECHNICAL EFFICACY

Stage 2.

Radiogenomic-based multiomic analysis reveals imaging intratumor heterogeneity phenotypes and therapeutic targets

Intratumor heterogeneity (ITH) profoundly affects therapeutic responses and clinical outcomes. However, the widespread methods for assessing ITH based on genomic sequencing or pathological slides, which rely on limited tissue samples, may lead to inaccuracies due to potential sampling biases. Using a newly established multicenter breast cancer radio-multiomic dataset (n = 1474) encompassing radiomic features extracted from dynamic contrast-enhanced magnetic resonance images, we formulated a noninvasive radiomics methodology to effectively investigate ITH. Imaging ITH (IITH) was associated with genomic and pathological ITH, predicting poor prognosis independently in breast cancer. Through multiomic analysis, we identified activated oncogenic pathways and metabolic dysregulation in high-IITH tumors. Integrated metabolomic and transcriptomic analyses highlighted ferroptosis as a vulnerability and potential therapeutic target of high-IITH tumors. Collectively, this work emphasizes the superiority of radiomics in capturing ITH. Furthermore, we provide insights into the biological basis of IITH and propose therapeutic targets for breast cancers with elevated IITH.

Intratumoral Heterogeneity of Fibrosarcoma Xenograft Models: Whole-Tumor Histogram Analysis of DWI and IVIM

RATIONAL AND OBJECTIVE

To explore the correlations of histogram parameters from diffusion-weighted imaging (DWI) and intravoxel incoherent motion (IVIM) with the heterogeneous features in a nude mouse model of fibrosarcoma.

MATERIALS AND METHODS

A total of 44 fibrosarcoma xenograft models were established by inoculating HT-1080 cells on the right thigh of mice and subjected tumors to DWI and IVIM imaging with 3.0 T MRI. Whole-tumor histogram parameters were calculated on apparent diffusion coefficient (ADC), pure diffusion coefficient (D), pseudo-diffusion coefficient (D*), and perfusion fraction (f). Heterogeneous features, including necrosis rate, cell density, Ki-67 labeling index (LI), and microvascular density (MVD) were measured. Intraclass correlation coefficients (ICC), Pearson or Spearman correlation tests, and receiver operating characteristics (ROC) were performed.

RESULTS

The 90th percentile, skewness and kurtosis of ADC and D histograms showed correlations with necrosis rate, and the highest correlation coefficient was found for D90th (r = 0.485). ADC and D histogram parameters showed correlations with cell density and Ki-67 LI; D90th showed the highest correlation coefficient with cell density (r = -0.504); and Dmedian showed the most significant correlation with Ki-67 LI (r = -0.525). D*skewness, D*kurtosis, D*90th, fmean, and fmedian showed correlations with MVD. ADC90th, ADCskewness, ADCkurtosis, D90th, and Dskewness showed significant differences between the low necrosis and high necrosis groups, and the combination model showed the best diagnostic ability (AUC = 0.882), with 97% sensitivity, and 72.7% specificity.

CONCLUSION

Whole-tumor histogram parameters of DWI and IVIM were correlated with heterogeneous features in nude murine models of fibrosarcoma.

Differentiating renal cell carcinoma and oncocytoma with volumetric MRI histogram analysis

OBJECTIVE

In this study, the utility of histogram parameters derived from diffusion-weighted imaging for differentiate renal cell carcinoma (RCC) from oncocytoma was investigated.

METHODS

This research tracked 126 individuals who were diagnosed with RCC and oncocytoma through histopathological analysis, using magnetic resonance imaging (MRI) assessments from 2015 to 2023. We observed various attributes of these patients, including demographic details, surgical records, pre-surgery MRI results, MRI apparent diffusion coefficient (ADC) histogram analysis, and post-surgery histopathological outcomes. Calculations of ADC measurements such as mean, minimum, and maximum in conjunction with the 5th, 10th, 25th, 50th, 75th, 90th, and 95th quantile points were made. In addition, we also noted the skewness, kurtosis, and variance of these data points.

RESULTS

The focus group for this investigation consisted of 75 male and 51 female patients. Out of these, 82 were diagnosed with RCC and 44 with oncocytoma. All ADC parameters including ADCmin, ADCmedian, ADCmean, and ADCmax, including the 5th, 10th, 25th, 50th, 75th, 90th, and 95th quantile divisions among the oncocytoma cohort were observed to be higher than the corresponding ones in the RCC group. A statistically meaningful difference was discovered between the minimum ADC value along with the 5th ranking of ADC measurements (p<0.001), in addition to mean of ADC (p=0.050), and the 10th (p=0.002) and 25th (p=0.015) quantiles of ADC data. When considering the region below the curve (AUC) in ROC analysis, the value of ADCmin was recorded as 0.739, with a sensitivity of 75.0%, and specificity of 68.2%.

CONCLUSION

To distinguish oncocytoma from RCC, it may be useful to conduct a whole-tumor histogram and textural analysis of ADC values.

Whole-tumor histogram analysis of apparent diffusion coefficient maps with machine learning algorithms for predicting histologic grade of sinonasal squamous cell carcinoma: a preliminary study

PURPOSE

Accurate histologic grade assessment is helpful for clinical decision making and prognostic assessment of sinonasal squamous cell carcinoma (SNSCC). This research aimed to explore whether whole-tumor histogram analysis of apparent diffusion coefficient (ADC) maps with machine learning algorithms can predict histologic grade of SNSCC.

METHODS

One hundred and forty-seven patients with pathologically diagnosed SNSCC formed this retrospective study. Sixty-six patients were low-grade (grade I/II) and eighty-one patients were high-grade (grade III). Eighteen histogram features were obtained from quantitative ADC maps. Additionally, the mean ADC value and clinical features were analyzed for comparison with histogram features. Machine learning algorithms were applied to build the best diagnostic model for predicting histological grade. The receiver operating characteristic (ROC) curve was used to evaluate the performance of each model prediction, and the area under the ROC curve (AUC) were analyzed.

RESULTS

The histogram model based on three features (10th Percentile, Mean, and 90th Percentile) with support vector machine (SVM) classifier demonstrated excellent diagnostic performance, with an AUC of 0.947 on the testing dataset. The AUC of the histogram model was similar to that of the mean ADC value model (0.947 vs 0.957; P = 0.7029). The poor diagnostic performance of the clinical model (AUC = 0.692) was improved by the combined model incorporating histogram features or mean ADC value (P < 0.05).

CONCLUSION

ADC histogram analysis improved the projection of SNSCC histologic grade, compared with clinical model. The complex histogram model had comparable but not better performance than mean ADC value model.

The Effectiveness of Volumetric MRI Histogram Analysis in Renal Cell Carcinoma

RATIONALE AND OBJECTIVES

This study investigated the utility of histogram parameters derived from diffusion-weighted imaging (DWI) for evaluating renal cell carcinoma (RCC) grading prior to surgery.

MATERIALS AND METHODS

This retrospective study included 88 patients who were histopathologically diagnosed with RCC and underwent magnetic resonance imaging (MRI) examinations. The patients were divided into two groups as well-differentiated (Group 1) and poorly differentiated (Group 2). Demographic data, preoperative MRI findings, MRI apparent diffusion coefficient (ADC) histogram analyzes, operation types, postoperative histopathological data and cancer stages of the patients were recorded. The histogram parameters of ADC values, comprising the mean, minimum, maximum, 5th, 10th, 25th, 50th, 75th, 90th, and 95th percentiles, as well as skewness, kurtosis, and variance, were calculated.

RESULTS

The study included 59 males and 29 women with an average age of 56.21 ± 1.33 years. There were 52 patients in Group 1 and 36 patients in Group 2. The ADCmin, ADCmean, ADCmax, 5th, 10th, 25th, 50th, 75th, 90th, and 95th percentiles of ADC values of the poorly differentiated group were all lower than those of the well-differentiated group. ADCmin and the 5th percentile of ADC values, as well as ADCmean and the 10th, 25th, 50th, and 75th percentiles of ADC values, showed a statistically significant difference (p < 0.05). The AUC, sensitivity, and specificity of the ADCmin value were 0.703, 56.3%, and 75.7%, respectively.

CONCLUSION

The present study indicated that histogram parameters generated from DWI were capable of differentiating between high-grade and low-grade RCC.

Radiomics model based on contrast-enhanced CT texture features for pretreatment prediction of overall survival in esophageal neuroendocrine carcinoma

Background

Limited studies have observed the prognostic value of CT images for esophageal neuroendocrine carcinoma (NEC) due to rare incidence and low treatment experience in clinical. In this study, the pretreatment enhanced CT texture features and clinical characteristics were investigated to predict the overall survival of esophageal NEC.

Methods

This retrospective study included 89 patients with esophageal NEC. The training and testing cohorts comprised 61 (70%) and 28 (30%) patients, respectively. A total of 402 radiomics features were extracted from the tumor region that segmented pretreatment venous phase CT images. The least absolute shrinkage and selection operator (LASSO) Cox regression was applied to feature dimension reduction, feature selection, and radiomics signature construction. A radiomics nomogram was constructed based on the radiomics signature and clinical risk factors using a multivariable Cox proportional regression. The performance of the nomogram for the pretreatment prediction of overall survival (OS) was evaluated for discrimination and calibration.

Results

Only the enhancement degree was an independent factor in clinical variable influenced OS. The radiomics signatures demonstrated good predictability for prognostic status discrimination. The radiomics nomogram integrating texture signatures was slightly superior to the nomogram derived from the combined model with a C-index of 0.844 (95%CI: 0.783-0.905) and 0.847 (95% CI: 0.782-0.912) in the training set, and 0.805 (95%CI: 0.707-0.903) and 0.745 (95% CI: 0.639-0.851) in the testing set, respectively.

Conclusion

The radiomics nomogram based on pretreatment CT radiomics signature had better prognostic power and predictability of the overall survival in patients with esophageal NEC than the model using combined variables.

Prostate Region-Wise Imaging Biomarker Profiles for Risk Stratification and Biochemical Recurrence Prediction

BACKGROUND

Identifying prostate cancer (PCa) patients with a worse prognosis and a higher risk of biochemical recurrence (BCR) is essential to guide treatment choices. Here, we aimed to identify possible imaging biomarker (perfusion/diffusion + radiomic features) profiles extracted from MRIs that were able to discriminate patients according to their risk or the occurrence of BCR 10 years after diagnosis, as well as to evaluate their predictive value with or without clinical data.

METHODS

Patients with localized PCa receiving neoadjuvant androgen deprivation therapy and radiotherapy were retrospectively evaluated. Imaging features were extracted from MRIs for each prostate region or for the whole gland. Univariate and multivariate analyses were conducted.

RESULTS

128 patients (mean [range] age, 71 [50-83] years) were included. Prostate region-wise imaging biomarker profiles mainly composed of radiomic features allowed discriminating risk groups and patients experiencing BCR. Heterogeneity-related radiomic features were increased in patients with worse prognosis and with BCR. Overall, imaging biomarkers profiles retained good predictive ability (AUC values superior to 0.725 in most cases), which generally improved when clinical data were included (particularly evident for the prediction of the BCR, with AUC values ranging from 0.841 to 0.877 for combined models and sensitivity values above 0.960) and when models were built per prostate region vs. the whole gland.

CONCLUSIONS

Prostate region-aware imaging profiles enable identification of patients with worse prognosis and with a higher risk of BCR, retaining higher predictive values when combined with clinical variables.

MRI-based breast cancer radiogenomics using RNA profiling: association with subtypes in a single-center prospective study

BACKGROUND

There are few prospective studies on the correlations between MRI features and whole RNA-sequencing data in breast cancer according to molecular subtypes. The purpose of our study was to explore the association between genetic profiles and MRI phenotypes of breast cancer and to identify imaging markers that influences the prognosis and treatment according to subtypes.

METHODS

From June 2017 to August 2018, MRIs of 95 women with invasive breast cancer were prospectively analyzed, using the breast imaging-reporting and data system and texture analysis. Whole RNA obtained from surgical specimens was analyzed using next-generation sequencing. The association between MRI features and gene expression profiles was analyzed in the entire tumor and subtypes. Gene networks, enriched functions, and canonical pathways were analyzed using Ingenuity Pathway Analysis. The P value for differential expression was obtained using a parametric F test comparing nested linear models and adjusted for multiple testing by reporting Q value.

RESULTS

In 95 participants (mean age, 53 years ± 11 [standard deviation]), mass lesion type was associated with upregulation of CCL3L1 (sevenfold) and irregular mass shape was associated with downregulation of MIR421 (sixfold). In estrogen receptor-positive cancer with mass lesion type, CCL3L1 (21-fold), SNHG12 (11-fold), and MIR206 (sevenfold) were upregulated, and MIR597 (265-fold), MIR126 (12-fold), and SOX17 (fivefold) were downregulated. In triple-negative breast cancer with increased standard deviation of texture analysis on precontrast T1-weighted imaging, CLEC3A (23-fold), SRGN (13-fold), HSPG2 (sevenfold), KMT2D (fivefold), and VMP1 (fivefold) were upregulated, and IGLC2 (73-fold) and PRDX4 (sevenfold) were downregulated (all, P < 0.05 and Q < 0.1). Gene network and functional analysis showed that mass type estrogen receptor-positive cancers were associated with cell growth, anti-estrogen resistance, and poor survival.

CONCLUSION

MRI characteristics are associated with the different expressions of genes related to metastasis, anti-drug resistance, and prognosis, depending on the molecular subtypes of breast cancer.

The genomic signature of resistance to platinum-containing neoadjuvant therapy based on single-cell data

BACKGROUND

Neoadjuvant chemotherapy (NACT) becomes the first-line option for advanced tumors, while patients who are not sensitive to it may not benefit. Therefore, it is important to screen patients suitable for NACT.

METHODS

Single-cell data of lung adenocarcinoma (LUAD) and esophageal squamous carcinoma (ESCC) before and after cisplatin-containing (CDDP) NACT and cisplatin IC50 data of tumor cell lines were analyzed to establish a CDDP neoadjuvant chemotherapy score (NCS). Differential analysis, GO, KEGG, GSVA and logistic regression models were performed by R. Survival analysis were applied to public databases. siRNA knockdown in A549, PC9, TE1 cell lines, qRT-PCR, western-blot, cck8 and EdU experiments were used for further verification in vitro.

RESULTS

485 genes were expressed differentially in tumor cells before and after neoadjuvant treatment for LUAD and ESCC. After combining the CDDP-associated genes, 12 genes, CAV2, PHLDA1, DUSP23, VDAC3, DSG2, SPINT2, SPATS2L, IGFBP3, CD9, ALCAM, PRSS23, PERP, were obtained and formed the NCS score. The higher the score, the more sensitive the patients were to CDDP-NACT. The NCS divided LUAD and ESCC into two groups. Based on differentially expressed genes, a model was constructed to predict the high and low NCS. CAV2, PHLDA1, ALCAM, CD9, IGBP3 and VDAC3 were significantly associated with prognosis. Finally, we demonstrated that the knockdown of CAV2, PHLDA1 and VDAC3 in A549, PC9 and TE1 significantly increased the sensitivity to cisplatin.

CONCLUSIONS

NCS scores and related predictive models for CDDP-NACT were developed and validated to assist in selecting patients who might benefit from it.

Identification of triple-negative breast cancer and androgen receptor expression based on histogram and texture analysis of dynamic contrast-enhanced MRI

BACKGROUND

Triple-negative breast cancer (TNBC) is highly malignant and has a poor prognosis due to the lack of effective therapeutic targets. Androgen receptor (AR) has been investigated as a possible therapeutic target. This study quantitatively assessed intratumor heterogeneity by histogram analysis of pharmacokinetic parameters and texture analysis on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) to discriminate TNBC from non-triple-negative breast cancer (non-TNBC) and to identify AR expression in TNBC.

METHODS

This retrospective study included 99 patients with histopathologically proven breast cancer (TNBC: 36, non-TNBC: 63) who underwent breast DCE-MRI before surgery. The pharmacokinetic parameters of DCE-MRI (K^trans, K_ep and V_e) and their corresponding texture parameters were calculated. The independent t-test, or Mann-Whitney U-test was used to compare quantitative parameters between TNBC and non-TNBC groups, and AR-positive (AR+) and AR-negative (AR-) TNBC groups. The parameters with significant difference between two groups were further involved in logistic regression analysis to build a prediction model for TNBC. The ROC analysis was conducted on each independent parameter and the TNBC predicting model for evaluating the discrimination performance. The area under the ROC curve (AUC), sensitivity and specificity were derived.

RESULTS

The binary logistic regression analysis revealed that K_{ep_Range} (p = 0.032) and V_{e_SumVariance} (p = 0.005) were significantly higher in TNBC than in non-TNBC. The AUC of the combined model for identifying TNBC was 0.735 (p < 0.001) with a cut-off value of 0.268, and its sensitivity and specificity were 88.89% and 52.38%, respectively. The value of K_{ep_Compactness2} (p = 0.049), K_{ep_SphericalDisproportion} (p = 0.049), and V_{e_GlcmEntropy} (p = 0.008) were higher in AR + TNBC group than in AR-TNBC group.

CONCLUSION

Histogram and texture analysis of breast lesions on DCE-MRI showed potential to identify TNBC, and the specific features can be possible predictors of AR expression, enhancing the ability to individualize the treatment of patients with TNBC.

Artificial intelligence-driven radiomics study in cancer: the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients' anatomy. However, the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians. Moreover, some potentially useful quantitative information in medical images, especially that which is not visible to the naked eye, is often ignored during clinical practice. In contrast, radiomics performs high-throughput feature extraction from medical images, which enables quantitative analysis of medical images and prediction of various clinical endpoints. Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis, demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine. However, radiomics remains in a developmental phase as numerous technical challenges have yet to be solved, especially in feature engineering and statistical modeling. In this review, we introduce the current utility of radiomics by summarizing research on its application in the diagnosis, prognosis, and prediction of treatment responses in patients with cancer. We focus on machine learning approaches, for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling. Furthermore, we introduce the stability, reproducibility, and interpretability of features, and the generalizability and interpretability of models. Finally, we offer possible solutions to current challenges in radiomics research.

The value of whole-volume apparent diffusion coefficient histogram analysis in preoperatively distinguishing intracranial solitary fibrous tumor and transitional meningioma

Purpose

To investigate the value of whole-volume apparent diffusion coefficient (ADC) histogram analysis in preoperatively distinguishing intracranial solitary fibrous tumors (SFT) from transitional meningiomas (TM), thereby assisting the establishment of the treatment protocol.

Methods

Preoperative diffusion-weighted imaging datasets of 24 patients with SFT and 28 patients with TM were used to extract whole-volume ADC histogram parameters, including variance, skewness, kurtosis, and mean, as well as 1st (AP1), 10th (AP10), 50th (AP50), 90th (AP90), and 99th (AP99) percentiles of ADC using MaZda software. The independent t-test or Mann-Whitney U test was used to compare the differences between ADC histogram parameters of SFT and TM. Receiver operating characteristic (ROC) curves were generated to evaluate the performance of significant ADC histogram parameters. Spearman's correlation coefficients were calculated to evaluate correlations between these parameters and the Ki-67 expression levels.

Results

SFT exhibited significantly higher variance, and lower AP1 and AP10 (all P < 0.05) than TM. The best diagnostic performance was obtained by variance, with an area under the ROC curve of 0.848 (0.722-0.933). However, there was no significant difference in skewness, kurtosis, mean, or other percentiles of ADC between the two groups (all P > 0.05). Significant correlations were also observed between the Ki-67 proliferation index and variance (r = 0.519), AP1 (r = -0.425), and AP10 (r = -0.372) (all P < 0.05).

Conclusion

Whole-volume ADC histogram analysis is a feasible tool for non-invasive preoperative discrimination between intracranial SFT and TM, with variance being the most promising prospective parameter.

Radiomics models based on multisequence MRI for predicting PD-1/PD-L1 expression in hepatocellular carcinoma

The purpose of this study was to explore the effectiveness of radiomics based on multisequence MRI in predicting the expression of PD-1/PD-L1 in hepatocellular carcinoma (HCC). One hundred and eight patients with HCC who underwent contrast-enhanced MRI 2 weeks before surgical resection were enrolled in this retrospective study. Corresponding paraffin sections were collected for immunohistochemistry to detect the expression of PD-1 and PD-L1. All patients were randomly divided into a training cohort and a validation cohort at a ratio of 7:3. Univariate and multivariate analyses were used to select potential clinical characteristics related to PD-1 and PD-L1 expression. Radiomics features were extracted from the axial fat-suppression T2-weighted imaging (FS-T2WI) images and the arterial phase and portal venous phase images from the axial dynamic contrast-enhanced MRI, and the corresponding feature sets were generated. The least absolute shrinkage and selection operator (LASSO) was used to select the optimal radiomics features for analysis. Logistic regression analysis was performed to construct single-sequence and multisequence radiomics and radiomic-clinical models. The predictive performance was judged by the area under the receiver operating characteristic curve (AUC) in the training and validation cohorts. In the whole cohort, PD-1 expression was positive in 43 patients, and PD-L1 expression was positive in 34 patients. The presence of satellite nodules served as an independent predictor of PD-L1 expression. The AUC values of the FS-T2WI, arterial phase, portal venous phase and multisequence models in predicting the expression of PD-1 were 0.696, 0.843, 0.863, and 0.946 in the training group and 0.669, 0.792, 0.800 and 0.815 in the validation group, respectively. The AUC values of the FS-T2WI, arterial phase, portal venous phase, multisequence and radiomic-clinical models in predicting PD-L1 expression were 0.731, 0.800, 0.800, 0.831 and 0.898 in the training group and 0.621, 0.743, 0.771, 0.810 and 0.779 in the validation group, respectively. The combined models showed better predictive performance. The results of this study suggest that a radiomics model based on multisequence MRI has the potential to predict the preoperative expression of PD-1 and PD-L1 in HCC, which could become an imaging biomarker for immune checkpoint inhibitor (ICI)-based treatment.

Prediction of prognostic factors in breast cancer: A noninvasive method utilizing histogram parameters derived from Adc maps

Objective

The aim of this study is to investigate the relationship between histogram parameters and prognostic factors of breast cancer and to reveal the diagnostic performance of histogram parameters in predicting prognostic factors status.

Materials and methods

Ninety-two patients with a confirmed histopathological diagnosis of breast cancer were included in the study. Magnetic resonance imaging (MRI) was performed using a 1.5T scanner and two different b values were used for diffusion-weighted imaging (DWI) (b values: 0 s/mm², b: 800 s/mm²). For 3D histogram analysis, regions of interest (ROI) were drawn each slice of the lesion on apparent diffusion coefficient (ADC) maps. The following data were derived from the histogram analysis data: percentiles, skewness, kurtosis, and entropy. The relationship between prognostic factors and histogram analysis data was investigated using the Kolmogorov-Smirnov test, Shapiro-Wilk test, skewness-kurtosis test, independent t-test, Mann-Whitney U test, and Kruskal-Wallis test. Receiver operator characteristic (ROC) curve analysis was performed to evaluate the diagnostic performance of the histogram parameters.

Results

ADC_max, kurtosis, and entropy parameters were statistically significantly correlated with tumor diameter (p = 0.002, p = 0.008, and p = 0.001, respectively). There was a significant difference in ADC_90% and ADC_max values, depending on estrogen receptor (ER) and progesterone receptor (PR) status. These values were lower in ER- and PR-positive than ER- and PR-negative patients (p = 0.02 and p = 0.001 vs. p = 0.018, p = 0.008). All ADC percentage values were lower in patients with a positive Ki-67 proliferation index as compared with those with a negative Ki-67 proliferation index (all p = 0.001). The entropy value was high in high-grade lesions and lesions with axillary involvement (p = 0.039 and p = 0.048, respectively). The highest area under the curve (AUC) for ER and PR status was calculated for the ADC_90% value with ROC curve analysis. The highest AUC for Ki-67 proliferation index was found for the ADC_50%.

Conclusion

Histogram analysis parameters derived from of ADC maps of whole lesions can reflect histopathological features of the tumors. Based on our study, it was concluded that histogram analysis parameters were related to the prognostic factors of the tumor.

Artificial Intelligence in CT and MR Imaging for Oncological Applications

Cancer care increasingly relies on imaging for patient management. The two most common cross-sectional imaging modalities in oncology are computed tomography (CT) and magnetic resonance imaging (MRI), which provide high-resolution anatomic and physiological imaging. Herewith is a summary of recent applications of rapidly advancing artificial intelligence (AI) in CT and MRI oncological imaging that addresses the benefits and challenges of the resultant opportunities with examples. Major challenges remain, such as how best to integrate AI developments into clinical radiology practice, the vigorous assessment of quantitative CT and MR imaging data accuracy, and reliability for clinical utility and research integrity in oncology. Such challenges necessitate an evaluation of the robustness of imaging biomarkers to be included in AI developments, a culture of data sharing, and the cooperation of knowledgeable academics with vendor scientists and companies operating in radiology and oncology fields. Herein, we will illustrate a few challenges and solutions of these efforts using novel methods for synthesizing different contrast modality images, auto-segmentation, and image reconstruction with examples from lung CT as well as abdome, pelvis, and head and neck MRI. The imaging community must embrace the need for quantitative CT and MRI metrics beyond lesion size measurement. AI methods for the extraction and longitudinal tracking of imaging metrics from registered lesions and understanding the tumor environment will be invaluable for interpreting disease status and treatment efficacy. This is an exciting time to work together to move the imaging field forward with narrow AI-specific tasks. New AI developments using CT and MRI datasets will be used to improve the personalized management of cancer patients.

Comparison of CT- and MRI-Based Quantification of Tumor Heterogeneity and Vascularity for Correlations with Prognostic Biomarkers and Survival Outcomes: A Single-Center Prospective Cohort Study

BACKGROUND

Tumor heterogeneity and vascularity can be noninvasively quantified using histogram and perfusion analyses on computed tomography (CT) and magnetic resonance imaging (MRI). We compared the association of histogram and perfusion features with histological prognostic factors and progression-free survival (PFS) in breast cancer patients on low-dose CT and MRI.

METHODS

This prospective study enrolled 147 women diagnosed with invasive breast cancer who simultaneously underwent contrast-enhanced MRI and CT before treatment. We extracted histogram and perfusion parameters from each tumor on MRI and CT, assessed associations between imaging features and histological biomarkers, and estimated PFS using the Kaplan-Meier analysis.

RESULTS

Out of 54 histogram and perfusion parameters, entropy on T2- and postcontrast T1-weighted MRI and postcontrast CT, and perfusion (blood flow) on CT were significantly associated with the status of subtypes, hormone receptors, and human epidermal growth factor receptor 2 (p < 0.05). Patients with high entropy on postcontrast CT showed worse PFS than patients with low entropy (p = 0.053) and high entropy on postcontrast CT negatively affected PFS in the Ki67-positive group (p = 0.046).

CONCLUSIONS

Low-dose CT histogram and perfusion analysis were comparable to MRI, and the entropy of postcontrast CT could be a feasible parameter to predict PFS in breast cancer patients.

Radiogenomics Reveals Correlation between Quantitative Texture Radiomic Features of Biparametric MRI and Hypoxia-Related Gene Expression in Men with Localised Prostate Cancer

OBJECTIVES

To perform multiscale correlation analysis between quantitative texture feature phenotypes of pre-biopsy biparametric MRI (bpMRI) and targeted sequence-based RNA expression for hypoxia-related genes.

MATERIALS AND METHODS

Images from pre-biopsy 3T bpMRI scans in clinically localised PCa patients of various risk categories (n = 15) were used to extract textural features. The genomic landscape of hypoxia-related gene expression was obtained using post-radical prostatectomy tissue for targeted RNA expression profiling using the TempO-sequence method. The nonparametric Games Howell test was used to correlate the differential expression of the important hypoxia-related genes with 28 radiomic texture features. Then, cBioportal was accessed, and a gene-specific query was executed to extract the Oncoprint genomic output graph of the selected hypoxia-related genes from The Cancer Genome Atlas (TCGA). Based on each selected gene profile, correlation analysis using Pearson's coefficients and survival analysis using Kaplan-Meier estimators were performed.

RESULTS

The quantitative bpMR imaging textural features, including the histogram and grey level co-occurrence matrix (GLCM), correlated with three hypoxia-related genes (ANGPTL4, VEGFA, and P4HA1) based on RNA sequencing using the TempO-Seq method. Further radiogenomic analysis, including data accessed from the cBioportal genomic database, confirmed that overexpressed hypoxia-related genes significantly correlated with a poor survival outcomes, with a median survival ratio of 81.11:133.00 months in those with and without alterations in genes, respectively.

CONCLUSION

This study found that there is a correlation between the radiomic texture features extracted from bpMRI in localised prostate cancer and the hypoxia-related genes that are differentially expressed. The analysis of expression data based on cBioportal revealed that these hypoxia-related genes, which were the focus of the study, are linked to an unfavourable survival outcomes in prostate cancer patients.

CT texture analysis in predicting treatment response and survival in patients with hepatocellular carcinoma treated with transarterial chemoembolization using random forest models

BACKGROUND

Using texture features derived from contrast-enhanced computed tomography (CT) combined with general imaging features as well as clinical information to predict treatment response and survival in patients with hepatocellular carcinoma (HCC) who received transarterial chemoembolization (TACE) treatment.

METHODS

From January 2014 to November 2022, 289 patients with HCC who underwent TACE were retrospectively reviewed. Their clinical information was documented. Their treatment-naïve contrast-enhanced CTs were retrieved and reviewed by two independent radiologists. Four general imaging features were evaluated. Texture features were extracted based on the regions of interest (ROIs) drawn on the slice with the largest axial diameter of all lesions using Pyradiomics v3.0.1. After excluding features with low reproducibility and low predictive value, the remaining features were selected for further analyses. The data were randomly divided in a ratio of 8:2 for model training and testing. Random forest classifiers were built to predict patient response to TACE treatment. Random survival forest models were constructed to predict overall survival (OS) and progress-free survival (PFS).

RESULTS

We retrospectively evaluated 289 patients (55.4 ± 12.4 years old) with HCC treated with TACE. Twenty features, including 2 clinical features (ALT and AFP levels), 1 general imaging feature (presence or absence of portal vein thrombus) and 17 texture features, were included in model construction. The random forest classifier achieved an area under the curve (AUC) of 0.947 with an accuracy of 89.5% for predicting treatment response. The random survival forest showed good predictive performance with out-of-bag error rate of 0.347 (0.374) and a continuous ranked probability score (CRPS) of 0.170 (0.067) for the prediction of OS (PFS).

CONCLUSIONS

Random forest algorithm based on texture features combined with general imaging features and clinical information is a robust method for predicting prognosis in patients with HCC treated with TACE, which may help avoid additional examinations and assist in treatment planning.