Radiogenomic signatures reveal multiscale intratumour heterogeneity associated with biological functions and survival in breast cancer

The cancer stem cells characteristics analysis of LGR5 + cells that influence lung cancer risk by using single cell RNA-seq analysis

Lung adenocarcinoma (LUAD) is the most popular lung cancer type with highly mortality. We performed a single cell RNA-seq analysis to explore characteristic of cancer stem cells in LUAD. We downloaded the single cell RNA-seq data (GSE149655) from the GEO database, the scRNA-seq analysis was performed by using the "Seurat" and "harmony" R package. The FindMarkers function and "ClusterProlifer" package was used for differentially expressed genes (DEGs) and function enrichment analysis. The protein-protein interaction and transcriptional regulatory network were performed by STRING and ChIPBase database. Immunohistochemistry tests to be used to observe differences in the expression of specific genes in LUAD and paracancerous tissue samples. BEAS-2B and A549 cells was used for vitro assay and the qRT-PCR, western blotting, wound healing, trans-well assays, EdU tests, and flow cytometry were performed. A total of 9 cell clusters were obtained after scRNA-seq analysis, in which the cancer stem cells had higher proportion in LUAD samples. Subsequently, function enrichment analysis revealed that the amino sugar and nucleotide sugar metabolism and DNA replication pathways were activated in cancer stem cells (CSCs), which were further sub-divided into 3 subtypes, the LGR5 + stem cell is a major contributor to cancer progression, its hub genes, such as HLA-DPB1, CD74, CTSH and HLA-DRB5 mediated the unique transcriptional state. In addition, the marker genes of three CSCs were also overexpressed in LUAD cells and the CXCL3 played an important role in mediating cell proliferation, apoptosis, migration and invasion of tumor. We performed a scRNA-seq analysis and identified the LGR5 + stem cell as a major contributor in LUAD progression, our findings are expected to provide new insights into the pathogenesis of LUAD.

CT-based Machine Learning Radiomics Modeling: Survival Prediction and Mechanism Exploration in Ovarian Cancer Patients

RATIONALE AND OBJECTIVES

To create a radiomics model based on computed tomography (CT) to predict overall survival in ovarian cancer patients. To combine Rad-score with genomic data to explore the association between gene expression and Rad-score.

MATERIALS AND METHODS

Imaging and clinical data from 455 patients with ovarian cancer were retrospectively analyzed. Patients were categorized into training cohort, validation cohort and test cohort. Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) methods were utilized to identify characteristics and develop the Rad-score. Radiomics models were developed and evaluated for predictive efficacy and clinical incremental value. Application of genomic data from the cancer genome atlas (TCGA) to reveal differential genes in different Rad-score groups. Screening hub genes and exploring the functions of hub genes through bioinformatics analysis and machine learning.

RESULTS

Prognostic models based on FIGO, tumor residual disease and Rad-score were developed. The receiver operating characteristic (ROC) curves showed that the 1, 3, and 5 year area under curves (AUCs) of the model were in the training group (0.816, 0.865 and 0.862, respectively), validation group (0.845, 0.877, 0.869, respectively) and test group (0.899, 0.906 and 0.869, respectively) had good predictive accuracy. Calibration curves showed good agreement between observations and predictions. Decision curve analysis revealed a high net benefit of the clinical-radiomics model. The clinical impact curve (CIC) showed good clinical applicability of the clinical-radiomics model. Analysis of sequencing data from the TCGA database revealed EMP1 as a hub gene for radiomics modeling. It revealed that its biological function may be associated with extracellular matrix organization and focal adhesion.

CONCLUSION

Prognostic models based on FIGO, Tumor residual disease, and Rad-score can effectively predict the overall survival (OS) of ovarian cancer patients. Rad-score may enable prognostic prediction of ovarian cancer patients by revealing the expression level of EMP1 and its biological function.

MVNMF: Multiview nonnegative matrix factorization for radio-multigenomic analysis in breast cancer prognosis

Radiogenomic research provides a deeper understanding of breast cancer biology by investigating the correlations between imaging phenotypes and genetic data. However, current radiogenomic research primarily focuses on the correlation between imaging phenotypes and single-genomic data (e.g., gene expression data), overlooking the potential of multi-genomics data to unveil more nuances in cancer characterization. To this end, we propose a multiview nonnegative matrix factorization (MVNMF) method for the radio-multigenomic analysis that identifies dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) features associated with multi-genomics data, including DNA copy number alterations, mutations, and mRNAs, each of which is independently predictive of cancer outcomes. MVNMF incorporates subspace learning and multiview regularization into a unified model to simultaneously select features and explore correlations. Subspace learning is utilized to identify representative radiomic features crucial for tumor analysis, while multiview regularization enables the learning of the correlation between the identified radiomic features and multi-genomics data. Experimental results showed that, for overall survival prediction in breast cancer, MVNMF classified patients into two distinct groups characterized by significant differences in survival (p = 0.0012). Furthermore, it achieved better performance with a C-index of 0.698 compared to the method without considering any genomics data (C-index = 0.528). MVNMF is an effective framework for identifying radiomic features linked to multi-genomics data, which improves its predictive power and provides a better understanding of the biological mechanisms underlying observed phenotypes. MVNMF offers a novel framework for prognostic prediction in breast cancer, with the potential to catalyze further radiogenomic/radio-multigenomic studies.

A bibliometric analysis of nasopharyngeal carcinoma radiomics: trends and insights

Background

Nasopharyngeal carcinoma (NPC) is a malignant tumor characterized by distinct geographic and pathological features. Enhancing diagnostic accuracy and timeliness in NPC is crucial for clinical implications. Radiomics has demonstrated significant potential in the clinical management of NPC. Nonetheless, a paucity of bibliometric studies has systematically examined the existing literature in th is domain. The objective of this study was to assess the current landscape and project future trends in NPC research.

Methods

This study conducted a search on English-language literature concerning the application of radiomics within the field of nasopharyngeal carcinoma (NPC) research from January 2015 to July 1, 2024, utilizing the Web of Science Core Collection (WoSCC) database. Bibliometric and visual analyses were performed using VOSviewer and CiteSpace software on publications related to countries/regions, authors, journals, references, and keywords.

Results

A total of 311 documents were retrieved, yielding 229 eligible documents after screening, comprising 209 articles and 20 reviews. Annual publications showed an upward trend, while citations revealed a generally declining trend. Notably, China contributed the most publications (n=175). Tian Jie and Dong Di each published 13 papers, and Zhang B was the most frequently co-cited author. Frontiers in Oncology published the most articles (n=25), and the International Journal of Radiation Oncology Biology Physics had the highest citation count (n=331). Sun Yat-sen University led institutional publications (n=39). The radiomics research in NPC focuses on survival prediction, texture analysis, and distant metastasis, and may guide future research directions.

Conclusion

The application of radiomics in NRC is growing annually, as indicated by bibliometric analysis. Radiomics has enhanced the precision of preoperative diagnosis, prediction, and prognosis in NRC. Bibliometric findings offer insights into radiomics research trends. However, creating extensive NPC datasets and bridging the research-to-clinical gap pose significant challenges. Future research should focus on these areas to advance the development.

An MRI Radiogenomic Signature to Characterize the Transcriptional Heterogeneity Associated with Prognosis and Biological Functions in Glioblastoma

BACKGROUND

The study sought to establish a radiogenomic signature to evaluate the transcriptional heterogeneity that reflects the prognosis and tumour-related biological functions of patients with glioblastoma.

METHODS

Transcriptional subclones were identified via fully unsupervised deconvolution of RNA sequencing. A genomic prognostic risk score was developed from transcriptional subclone proportions in the development dataset (n = 532) and independently verified in the testing dataset (n = 225). Multimodal magnetic resonance imaging (MRI) analysis involved feature extraction from three distinct anatomical regions across four imaging sequences. Key features were selected to construct a radiogenomic signature predictive of the genomic risk score in the radiogenomic dataset (n = 99), with subsequent survival analysis conducted in the image testing dataset (n = 233).

RESULTS

A total of 8 transcriptional subclones were identified, of which the metabolic pathway subclone and spinocerebellar ataxia subclone were independent risk factors for overall survival. The genomic risk score effectively differentiated patient subgroups with divergent survival outcomes in both development (p < 0.001) and testing datasets (p = 0.0003). Nineteen radiomic features were selected to construct a radiogenomic signature, with these features being linked to hallmark cancer pathways and the malignant behaviours of cancer cells. The radiogenomic signature predicted overall survival in the image testing dataset (hazard ratios (HR) = 1.67, p = 0.011).

CONCLUSIONS

A prognostic radiogenomic signature was established and verified to characterize transcriptional subclones with underlying biological functions in glioblastoma.

Radiogenomic Landscape of Metastatic Endocrine-Positive Breast Cancer Resistant to Aromatase Inhibitors

Breast cancer poses a significant global health challenge and includes various subtypes, such as endocrine-positive, HER2-positive, and triple-negative. Endocrine-positive breast cancer, characterized by estrogen and progesterone receptors, is commonly treated with aromatase inhibitors. However, resistance to these inhibitors can hinder patient outcomes due to genetic and epigenetic alterations, mutations in the estrogen receptor 1 gene, and changes in signaling pathways. Radiogenomics combines imaging techniques like MRI and CT scans with genomic profiling methods to identify radiographic biomarkers associated with resistance. This approach enhances our understanding of resistance mechanisms and metastasis patterns, linking them to specific genomic profiles and common metastasis sites like the bone and brain. By integrating radiogenomic data, personalized treatment strategies can be developed, improving predictive and prognostic capabilities. Advancements in imaging and genomic technologies offer promising avenues for enhancing radiogenomic research. A thorough understanding of resistance mechanisms is crucial for developing effective treatment strategies, making radiogenomics a valuable integrative approach in personalized medicine that aims to improve clinical outcomes for patients with metastatic endocrine-positive breast cancer.

A feature fusion method based on radiomic features and revised deep features for improving tumor prediction in ultrasound images

BACKGROUND

Radiomic features and deep features are both vitally helpful for the accurate prediction of tumor information in breast ultrasound. However, whether integrating radiomic features and deep features can improve the prediction performance of tumor information is unclear.

METHODS

A feature fusion method based on radiomic features and revised deep features was proposed to predict tumor information. Radiomic features were extracted from the tumor region on ultrasound images, and the optimal radiomic features were subsequently selected based on Gini score. Revised deep features, which were extracted using the revised CNN models integrating prior information, were combined with radiomic features to build a logistic regression classifier for tumor prediction. The performance was evaluated using area under the receiver operating characteristic (ROC) curve (AUC).

RESULTS

The results showed that the proposed feature fusion method (AUC = 0.9845) obtained better prediction performance than that based on radiomic features (AUC = 0.9796) or deep features (AUC = 0.9342).

CONCLUSIONS

Our results demonstrate that the proposed feature fusion framework integrating the radiomic features and revised deep features is an efficient method to improve the prediction performance of tumor information.

From Images to Genes: Radiogenomics Based on Artificial Intelligence to Achieve Non-Invasive Precision Medicine in Cancer Patients

With the increasing demand for precision medicine in cancer patients, radiogenomics emerges as a promising frontier. Radiogenomics is originally defined as a methodology for associating gene expression information from high-throughput technologies with imaging phenotypes. However, with advancements in medical imaging, high-throughput omics technologies, and artificial intelligence, both the concept and application of radiogenomics have significantly broadened. In this review, the history of radiogenomics is enumerated, related omics technologies, the five basic workflows and their applications across tumors, the role of AI in radiogenomics, the opportunities and challenges from tumor heterogeneity, and the applications of radiogenomics in tumor immune microenvironment. The application of radiogenomics in positron emission tomography and the role of radiogenomics in multi-omics studies is also discussed. Finally, the challenges faced by clinical transformation, along with future trends in this field is discussed.

Immune response and mesenchymal transition of papillary thyroid carcinoma reflected in ultrasonography features assessed by radiologists and deep learning

INTRODUCTION

Ultrasonography (US) features of papillary thyroid cancers (PTCs) are used to select nodules for biopsy due to their association with tumor behavior. However, the molecular biological mechanisms that lead to the characteristic US features of PTCs are largely unknown.

OBJECTIVES

This study aimed to investigate the molecular biological mechanisms behind US features assessed by radiologists and three convolutional neural networks (CNN) through transcriptome analysis.

METHODS

Transcriptome data from 273 PTC tissue samples were generated and differentially expressed genes (DEGs) were identified according to US feature. Pathway enrichment analyses were also conducted by gene set enrichment analysis (GSEA) and ClusterProfiler according to assessments made by radiologists and three CNNs - CNN1 (ResNet50), CNN2 (ResNet101) and CNN3 (VGG16). Signature gene scores for PTCs were calculated by single-sample GSEA (ssGSEA).

RESULTS

Individual suspicious US features consistently suggested an upregulation of genes related to immune response and epithelial-mesenchymal transition (EMT). Likewise, PTCs assessed as positive by radiologists and three CNNs showed the coordinate enrichment of similar gene sets with abundant immune and stromal components. However, PTCs assessed as positive by radiologists had the highest number of DEGs, and those assessed as positive by CNN3 had more diverse DEGs and gene sets compared to CNN1 or CNN2. The percentage of PTCs assessed as positive or negative concordantly by radiologists and three CNNs was 85.6% (231/273) and 7.1% (3/273), respectively.

CONCLUSION

US features assessed by radiologists and CNNs revealed molecular biologic features and tumor microenvironment in PTCs.

Radioproteomics stratifies molecular response to antifibrotic treatment in pulmonary fibrosis

Antifibrotic therapy with nintedanib is the clinical mainstay in the treatment of progressive fibrosing interstitial lung disease (ILD). High-dimensional medical image analysis, known as radiomics, provides quantitative insights into organ-scale pathophysiology, generating digital disease fingerprints. Here, we performed an integrative analysis of radiomic and proteomic profiles (radioproteomics) to assess whether changes in radiomic signatures can stratify the degree of antifibrotic response to nintedanib in (experimental) fibrosing ILD. Unsupervised clustering of delta radiomic profiles revealed 2 distinct imaging phenotypes in mice treated with nintedanib, contrary to conventional densitometry readouts, which showed a more uniform response. Integrative analysis of delta radiomics and proteomics demonstrated that these phenotypes reflected different treatment response states, as further evidenced on transcriptional and cellular levels. Importantly, radioproteomics signatures paralleled disease- and drug-related biological pathway activity with high specificity, including extracellular matrix (ECM) remodeling, cell cycle activity, wound healing, and metabolic activity. Evaluation of the preclinical molecular response-defining features, particularly those linked to ECM remodeling, in a cohort of nintedanib-treated fibrosing patients with ILD, accurately stratified patients based on their extent of lung function decline. In conclusion, delta radiomics has great potential to serve as a noninvasive and readily accessible surrogate of molecular response phenotypes in fibrosing ILD. This could pave the way for personalized treatment strategies and improved patient outcomes.

Research Trends and Evolution in Radiogenomics (2005-2023): Bibliometric Analysis

BACKGROUND

Radiogenomics is an emerging technology that integrates genomics and medical image-based radiomics, which is considered a promising approach toward achieving precision medicine.

OBJECTIVE

The aim of this study was to quantitatively analyze the research status, dynamic trends, and evolutionary trajectory in the radiogenomics field using bibliometric methods.

METHODS

The relevant literature published up to 2023 was retrieved from the Web of Science Core Collection. Excel was used to analyze the annual publication trend. VOSviewer was used for constructing the keywords co-occurrence network and the collaboration networks among countries and institutions. CiteSpace was used for citation keywords burst analysis and visualizing the references timeline.

RESULTS

A total of 3237 papers were included and exported in plain-text format. The annual number of publications showed an increasing annual trend. China and the United States have published the most papers in this field, with the highest number of citations in the United States and the highest average number per item in the Netherlands. Keywords burst analysis revealed that several keywords, including "big data," "magnetic resonance spectroscopy," "renal cell carcinoma," "stage," and "temozolomide," experienced a citation burst in recent years. The timeline views demonstrated that the references can be categorized into 8 clusters: lower-grade glioma, lung cancer histology, lung adenocarcinoma, breast cancer, radiation-induced lung injury, epidermal growth factor receptor mutation, late radiotherapy toxicity, and artificial intelligence.

CONCLUSIONS

The field of radiogenomics is attracting increasing attention from researchers worldwide, with the United States and the Netherlands being the most influential countries. Exploration of artificial intelligence methods based on big data to predict the response of tumors to various treatment methods represents a hot spot research topic in this field at present.

CAM3.0: determining cell type composition and expression from bulk tissues with fully unsupervised deconvolution

MOTIVATION

Complex tissues are dynamic ecosystems consisting of molecularly distinct yet interacting cell types. Computational deconvolution aims to dissect bulk tissue data into cell type compositions and cell-specific expressions. With few exceptions, most existing deconvolution tools exploit supervised approaches requiring various types of references that may be unreliable or even unavailable for specific tissue microenvironments.

RESULTS

We previously developed a fully unsupervised deconvolution method-Convex Analysis of Mixtures (CAM), that enables estimation of cell type composition and expression from bulk tissues. We now introduce CAM3.0 tool that improves this framework with three new and highly efficient algorithms, namely, radius-fixed clustering to identify reliable markers, linear programming to detect an initial scatter simplex, and a smart floating search for the optimum latent variable model. The comparative experimental results obtained from both realistic simulations and case studies show that the CAM3.0 tool can help biologists more accurately identify known or novel cell markers, determine cell proportions, and estimate cell-specific expressions, complementing the existing tools particularly when study- or datatype-specific references are unreliable or unavailable.

AVAILABILITY AND IMPLEMENTATION

The open-source R Scripts of CAM3.0 is freely available at https://github.com/ChiungTingWu/CAM3/(https://github.com/Bioconductor/Contributions/issues/3205). A user's guide and a vignette are provided.

Deconvolution-Based Pharmacokinetic Analysis to Improve the Prediction of Pathological Information of Breast Cancer

Pharmacokinetic (PK) parameters, revealing changes in the tumor microenvironment, are related to the pathological information of breast cancer. Tracer kinetic models (e.g., Tofts-Kety model) with a nonlinear least square solver are commonly used to estimate PK parameters. However, the method is sensitive to noise in images. To relieve the effects of noise, a deconvolution (DEC) method, which was validated on synthetic concentration-time series, was proposed to accurately calculate PK parameters from breast dynamic contrast-enhanced magnetic resonance imaging. A time-to-peak-based tumor partitioning method was used to divide the whole tumor into three tumor subregions with different kinetic patterns. Radiomic features were calculated from the tumor subregion and whole tumor-based PK parameter maps. The optimal features determined by the fivefold cross-validation method were used to build random forest classifiers to predict molecular subtypes, Ki-67, and tumor grade. The diagnostic performance evaluated by the area under the receiver operating characteristic curve (AUC) was compared between the subregion and whole tumor-based PK parameters. The results showed that the DEC method obtained more accurate PK parameters than the Tofts method. Moreover, the results showed that the subregion-based Ktrans (best AUCs = 0.8319, 0.7032, 0.7132, 0.7490, 0.8074, and 0.6950) achieved a better diagnostic performance than the whole tumor-based Ktrans (AUCs = 0.8222, 0.6970, 0.6511, 0.7109, 0.7620, and 0.5894) for molecular subtypes, Ki-67, and tumor grade. These findings indicate that DEC-based Ktrans in the subregion has the potential to accurately predict molecular subtypes, Ki-67, and tumor grade.

Double Transformer Super-Resolution for Breast Cancer ADC Images

Diffusion-weighted imaging (DWI) has been extensively explored in guiding the clinic management of patients with breast cancer. However, due to the limited resolution, accurately characterizing tumors using DWI and the corresponding apparent diffusion coefficient (ADC) is still a challenging problem. In this paper, we aim to address the issue of super-resolution (SR) of ADC images and evaluate the clinical utility of SR-ADC images through radiomics analysis. To this end, we propose a novel double transformer-based network (DTformer) to enhance the resolution of ADC images. More specifically, we propose a symmetric U-shaped encoder-decoder network with two different types of transformer blocks, named as UTNet, to extract deep features for super-resolution. The basic backbone of UTNet is composed of a locally-enhanced Swin transformer block (LeSwin-T) and a convolutional transformer block (Conv-T), which are responsible for capturing long-range dependencies and local spatial information, respectively. Additionally, we introduce a residual upsampling network (RUpNet) to expand image resolution by leveraging initial residual information from the original low-resolution (LR) images. Extensive experiments show that DTformer achieves superior SR performance. Moreover, radiomics analysis reveals that improving the resolution of ADC images is beneficial for tumor characteristic prediction, such as histological grade and human epidermal growth factor receptor 2 (HER2) status.

AI-powered interpretable imaging phenotypes noninvasively characterize tumor microenvironment associated with diverse molecular signatures and survival in breast cancer

BACKGROUND AND OBJECTIVES

Tumor microenvironment (TME) is a determining factor in decision-making and personalized treatment for breast cancer, which is highly intra-tumor heterogeneous (ITH). However, the noninvasive imaging phenotypes of TME are poorly understood, even invasive genotypes have been largely known in breast cancer.

METHODS

Here, we develop an artificial intelligence (AI)-driven approach for noninvasively characterizing TME by integrating the predictive power of deep learning with the explainability of human-interpretable imaging phenotypes (IMPs) derived from 4D dynamic imaging (DCE-MRI) of 342 breast tumors linked to genomic and clinical data, which connect cancer phenotypes to genotypes. An unsupervised dual-attention deep graph clustering model (DGCLM) is developed to divide bulk tumor into multiple spatially segregated and phenotypically consistent subclusters. The IMPs ranging from spatial heterogeneity to kinetic heterogeneity are leveraged to capture architecture, interaction, and proximity between intratumoral subclusters.

RESULTS

We demonstrate that our IMPs correlate with well-known markers of TME and also can predict distinct molecular signatures, including expression of hormone receptor, epithelial growth factor receptor and immune checkpoint proteins, with the performance of accuracy, reliability and transparency superior to recent state-of-the-art radiomics and 'black-box' deep learning methods. Moreover, prognostic value is confirmed by survival analysis accounting for IMPs.

CONCLUSIONS

Our approach provides an interpretable, quantitative, and comprehensive perspective to characterize TME in a noninvasive and clinically relevant manner.

Efficient estimation of pharmacokinetic parameters from breast dynamic contrast-enhanced MRI based on a convolutional neural network for predicting molecular subtypes

Objective. Tracer kinetic models allow for estimating pharmacokinetic (PK) parameters, which are related to pathological characteristics, from breast dynamic contrast-enhanced magnetic resonance imaging. However, existing tracer kinetic models subject to inaccuracy are time-consuming for PK parameters estimation. This study aimed to accurately and efficiently estimate PK parameters for predicting molecular subtypes based on convolutional neural network (CNN).Approach. A CNN integrating global and local features (GL-CNN) was trained using synthetic data where known PK parameters map was used as the ground truth, and subsequently used to directly estimate PK parameters (volume transfer constantKtransand flux rate constantKep) map. The accuracy assessed by the peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and concordance correlation coefficient (CCC) was compared between the GL-CNN and Tofts-based PK parameters in synthetic data. Radiomic features were calculated from the PK parameters map in 208 breast tumors. A random forest classifier was constructed to predict molecular subtypes using a discovery cohort (n= 144). The diagnostic performance evaluated on a validation cohort (n= 64) using the area under the receiver operating characteristic curve (AUC) was compared between the GL-CNN and Tofts-based PK parameters.Main results. The average PSNR (48.8884), SSIM (0.9995), and CCC (0.9995) between the GL-CNN-basedKtransmap and ground truth were significantly higher than those between the Tofts-basedKtransmap and ground truth. The GL-CNN-basedKtransobtained significantly better diagnostic performance (AUCs = 0.7658 and 0.8528) than the Tofts-basedKtransfor luminal B and HER2 tumors. The GL-CNN method accelerated the computation by speed approximately 79 times compared to the Tofts method for the whole breast of all patients.Significance. Our results indicate that the GL-CNN method can be used to accurately and efficiently estimate PK parameters for predicting molecular subtypes.

Components of the tumor immune microenvironment based on m-IHC correlate with prognosis and subtype of triple-negative breast cancer

BACKGROUND AND AIM

The spatial distribution and interactions of cells in the tumor immune microenvironment (TIME) might be related to the different responses of triple-negative breast cancer (TNBC) to immunomodulators. The potential of multiplex IHC (m-IHC) in evaluating the TIME has been reported, but the efficacy is insufficient. We aimed to research whether m-IHC results could be used to reflect the TIME, and thus to predict prognosis and complement the TNBC subtyping system.

METHODS

The clinical, imaging, and prognosis data for 86 TNBC patients were retrospectively reviewed. CD3, CD4, CD8, Foxp3, PD-L1, and Pan-CK markers were stained by m-IHC. Particular cell spatial distributions and interactions in the TIME were evaluated with the HALO multispectral analysis platform. Then, we calculated the prognostic value of components of the TIME and their correlations with TNBC transcriptomic subtypes and MRI radiomic features reflecting TNBC subtypes.

RESULTS

The components of the TIME score were established by m-IHC and demonstrated positive prognostic value for TNBC (p = 0.0047, 0.039, <0.0001 for DMFS, RFS, and OS). The score was calculated from several indicators, including Treg% in the tumor core (TC) or stromal area (SA), PD-L1+ cell% in the SA, CD3 + cell% in the TC, and PD-L1+ /CD8+ cells in the invasive margin and SA. According to the TNBC subtyping system, a few TIME indicators were significantly different in different subtypes and significantly correlated with MRI radiomic features reflecting TNBC subtypes.

CONCLUSION

We demonstrated that the m-IHC-based quantitative score and indicators related to the spatial distribution and interactions of cells in the TIME can aid in the accurate diagnosis of TNBC in terms of prognosis and classification.

Radiogenomic analysis of cellular tumor-stroma heterogeneity as a prognostic predictor in breast cancer

BACKGROUND

The tumor microenvironment and intercellular communication between solid tumors and the surrounding stroma play crucial roles in cancer initiation, progression, and prognosis. Radiomics provides clinically relevant information from radiological images; however, its biological implications in uncovering tumor pathophysiology driven by cellular heterogeneity between the tumor and stroma are largely unknown. We aimed to identify radiogenomic signatures of cellular tumor-stroma heterogeneity (TSH) to improve breast cancer management and prognosis analysis.

METHODS

This retrospective multicohort study included five datasets. Cell subpopulations were estimated using bulk gene expression data, and the relative difference in cell subpopulations between the tumor and stroma was used as a biomarker to categorize patients into good- and poor-survival groups. A radiogenomic signature-based model utilizing dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) was developed to target TSH, and its clinical significance in relation to survival outcomes was independently validated.

RESULTS

The final cohorts of 1330 women were included for cellular TSH biomarker identification (n = 112, mean age, 57.3 years ± 14.6) and validation (n = 886, mean age, 58.9 years ± 13.1), radiogenomic signature of TSH identification (n = 91, mean age, 55.5 years ± 11.4), and prognostic (n = 241) assessments. The cytotoxic lymphocyte biomarker differentiated patients into good- and poor-survival groups (p < 0.0001) and was independently validated (p = 0.014). The good survival group exhibited denser cell interconnections. The radiogenomic signature of TSH was identified and showed a positive association with overall survival (p = 0.038) and recurrence-free survival (p = 3 × 10-4).

CONCLUSION

Radiogenomic signatures provide insights into prognostic factors that reflect the imbalanced tumor-stroma environment, thereby presenting breast cancer-specific biological implications and prognostic significance.

Triple-negative breast cancer: predictive model of early recurrence based on MRI features

AIM

To develop an integrated model based on preoperative magnetic resonance imaging (MRI) features for predicting early recurrence in patients with triple-negative breast cancer (TNBC).

MATERIALS AND METHODS

Women with TNBC who underwent breast MRI and surgery between 2009 and 2019 were evaluated retrospectively. Two breast radiologists reviewed MRI images independently based on the Breast Imaging Reporting and Data System Lexicon (BI-RADS), and classified the breast oedema scores on T2-weighted imaging (WI) as no oedema, peritumoural oedema, prepectoral oedema, or subcutaneous oedema. The relationship between disease-free survival (DFS) and MRI features was analysed by Cox regression, and a nomogram model was generated based on the results.

RESULTS

150 patients with TNBC were included and divided into a training cohort (n=78) and validation cohort (n=72). MRI features including subcutaneous oedema and rim enhancement showed a tendency to worsen DFS in univariate analysis. Multivariate analysis showed that subcutaneous oedema (p=0.049, HR [95% confidence interval {CI} = 8.24 [1.01-67.52]) and rim enhancement (p=0.016, HR [95% CI] = 4.38 [1.32-14.54]) were independent predictors for DFS. In the nomogram, the areas under the curves (AUCs) of the training cohort was 0.808, and that of the validation cohort was 0.875.

CONCLUSION

The presence of subcutaneous oedema or rim enhancement on preoperative breast MRI was shown to be a good predictor of poor survival outcomes in patients with TNBC.

Recent advancements in artificial intelligence for breast cancer: Image augmentation, segmentation, diagnosis, and prognosis approaches

Breast cancer is a significant global health burden, with increasing morbidity and mortality worldwide. Early screening and accurate diagnosis are crucial for improving prognosis. Radiographic imaging modalities such as digital mammography (DM), digital breast tomosynthesis (DBT), magnetic resonance imaging (MRI), ultrasound (US), and nuclear medicine techniques, are commonly used for breast cancer assessment. And histopathology (HP) serves as the gold standard for confirming malignancy. Artificial intelligence (AI) technologies show great potential for quantitative representation of medical images to effectively assist in segmentation, diagnosis, and prognosis of breast cancer. In this review, we overview the recent advancements of AI technologies for breast cancer, including 1) improving image quality by data augmentation, 2) fast detection and segmentation of breast lesions and diagnosis of malignancy, 3) biological characterization of the cancer such as staging and subtyping by AI-based classification technologies, 4) prediction of clinical outcomes such as metastasis, treatment response, and survival by integrating multi-omics data. Then, we then summarize large-scale databases available to help train robust, generalizable, and reproducible deep learning models. Furthermore, we conclude the challenges faced by AI in real-world applications, including data curating, model interpretability, and practice regulations. Besides, we expect that clinical implementation of AI will provide important guidance for the patient-tailored management.

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.

Radiogenomics: a key component of precision cancer medicine

Radiogenomics, focusing on the relationship between genomics and imaging phenotypes, has been widely applied to address tumour heterogeneity and predict immune responsiveness and progression. It is an inevitable consequence of current trends in precision medicine, as radiogenomics costs less than traditional genetic sequencing and provides access to whole-tumour information rather than limited biopsy specimens. By providing voxel-by-voxel genetic information, radiogenomics can allow tailored therapy targeting a complete, heterogeneous tumour or set of tumours. In addition to quantifying lesion characteristics, radiogenomics can also be used to distinguish benign from malignant entities, as well as patient characteristics, to better stratify patients according to disease risk, thereby enabling more precise imaging and screening. Here, we have characterised the radiogenomic application in precision medicine using a multi-omic approach. we outline the main applications of radiogenomics in diagnosis, treatment planning and evaluations in the field of oncology with the aim of developing quantitative and personalised medicine. Finally, we discuss the challenges in the field of radiogenomics and the scope and clinical applicability of these methods.

Prediction of prostate tumour hypoxia using pre-treatment MRI-derived radiomics: preliminary findings

PURPOSE

To develop a machine learning (ML) model based on radiomic features (RF) extracted from whole prostate gland magnetic resonance imaging (MRI) for prediction of tumour hypoxia pre-radiotherapy.

MATERIAL AND METHODS

Consecutive patients with high-grade prostate cancer and pre-treatment MRI treated with radiotherapy between 01/12/2007 and 1/08/2013 at two cancer centres were included. Cancers were dichotomised as normoxic or hypoxic using a biopsy-based 32-gene hypoxia signature (Ragnum signature). Prostate segmentation was performed on axial T2-weighted (T2w) sequences using RayStation (v9.1). Histogram standardisation was applied prior to RF extraction. PyRadiomics (v3.0.1) was used to extract RFs for analysis. The cohort was split 80:20 into training and test sets. Six different ML classifiers for distinguishing hypoxia were trained and tuned using five different feature selection models and fivefold cross-validation with 20 repeats. The model with the highest mean validation area under the curve (AUC) receiver operating characteristic (ROC) curve was tested on the unseen set, and AUCs were compared via DeLong test with 95% confidence interval (CI).

RESULTS

195 patients were included with 97 (49.7%) having hypoxic tumours. The hypoxia prediction model with best performance was derived using ridge regression and had a test AUC of 0.69 (95% CI: 0.14). The test AUC for the clinical-only model was lower (0.57), but this was not statistically significant (p = 0.35). The five selected RFs included textural and wavelet-transformed features.

CONCLUSION

Whole prostate MRI-radiomics has the potential to non-invasively predict tumour hypoxia prior to radiotherapy which may be helpful for individualised treatment optimisation.

Imaging in metastatic breast cancer, CT, PET/CT, MRI, WB-DWI, CCA: review and new perspectives

BACKGROUND

Breast cancer is the most frequent cancer in women and remains the second leading cause of death in Western countries. It represents a heterogeneous group of diseases with diverse tumoral behaviour, treatment responsiveness and prognosis. While major progress in diagnosis and treatment has resulted in a decline in breast cancer-related mortality, some patients will relapse and prognosis in this cohort of patients remains poor. Treatment is determined according to tumor subtype; primarily hormone receptor status and HER2 expression. Menopausal status and site of disease relapse are also important considerations in treatment protocols.

MAIN BODY

Staging and repeated evaluation of patients with metastatic breast cancer are central to the accurate assessment of disease extent at diagnosis and during treatment; guiding ongoing clinical management. Advances have been made in the diagnostic and therapeutic fields, particularly with new targeted therapies. In parallel, oncological imaging has evolved exponentially with the development of functional and anatomical imaging techniques. Consistent, reproducible and validated methods of assessing response to therapy is critical in effectively managing patients with metastatic breast cancer.

CONCLUSION

Major progress has been made in oncological imaging over the last few decades. Accurate disease assessment at diagnosis and during treatment is important in the management of metastatic breast cancer. CT (and BS if appropriate) is generally widely available, relatively cheap and sufficient in many cases. However, several additional imaging modalities are emerging and can be used as adjuncts, particularly in pregnancy or other diagnostically challenging cases. Nevertheless, no single imaging technique is without limitation. The authors have evaluated the vast array of imaging techniques - individual, combined parametric and multimodal - that are available or that are emerging in the management of metastatic breast cancer. This includes WB DW-MRI, CCA, novel PET breast cancer-epitope specific radiotracers and radiogenomics.

Analyzing breast cancer invasive disease event classification through explainable artificial intelligence

Introduction

Recently, accurate machine learning and deep learning approaches have been dedicated to the investigation of breast cancer invasive disease events (IDEs), such as recurrence, contralateral and second cancers. However, such approaches are poorly interpretable.

Methods

Thus, we designed an Explainable Artificial Intelligence (XAI) framework to investigate IDEs within a cohort of 486 breast cancer patients enrolled at IRCCS Istituto Tumori "Giovanni Paolo II" in Bari, Italy. Using Shapley values, we determined the IDE driving features according to two periods, often adopted in clinical practice, of 5 and 10 years from the first tumor diagnosis.

Results

Age, tumor diameter, surgery type, and multiplicity are predominant within the 5-year frame, while therapy-related features, including hormone, chemotherapy schemes and lymphovascular invasion, dominate the 10-year IDE prediction. Estrogen Receptor (ER), proliferation marker Ki67 and metastatic lymph nodes affect both frames.

Discussion

Thus, our framework aims at shortening the distance between AI and clinical practice.

Image-based deep learning identifies glioblastoma risk groups with genomic and transcriptomic heterogeneity: a multi-center study

OBJECTIVES

To develop and validate a deep learning imaging signature (DLIS) for risk stratification in patients with multiforme (GBM), and to investigate the biological pathways and genetic alterations underlying the DLIS.

METHODS

The DLIS was developed from multi-parametric MRI based on a training set (n = 600) and validated on an internal validation set (n = 164), an external test set 1 (n = 100), an external test set 2 (n = 161), and a public TCIA set (n = 88). A co-profiling framework based on a radiogenomics analysis dataset (n = 127) using multiscale high-dimensional data, including imaging, transcriptome, and genome, was established to uncover the biological pathways and genetic alterations underpinning the DLIS.

RESULTS

The DLIS was associated with survival (log-rank p < 0.001) and was an independent predictor (p < 0.001). The integrated nomogram incorporating the DLIS achieved improved C indices than the clinicomolecular nomogram (net reclassification improvement 0.39, p < 0.001). DLIS significantly correlated with core pathways of GBM (apoptosis and cell cycle-related P53 and RB pathways, and cell proliferation-related RTK pathway), as well as key genetic alterations (del_CDNK2A). The prognostic value of DLIS-correlated genes was externally confirmed on TCGA/CGGA sets (p < 0.01).

CONCLUSIONS

Our study offers a biologically interpretable deep learning predictor of survival outcomes in patients with GBM, which is crucial for better understanding GBM patient's prognosis and guiding individualized treatment.

KEY POINTS

• MRI-based deep learning imaging signature (DLIS) stratifies GBM into risk groups with distinct molecular characteristics. • DLIS is associated with P53, RB, and RTK pathways and del_CDNK2A mutation. • The prognostic value of DLIS-correlated pathway genes is externally demonstrated.

Radiogenomic analysis of prediction HER2 status in breast cancer by linking ultrasound radiomic feature module with biological functions

BACKGROUND

Human epidermal growth factor receptor 2 (HER2) overexpressed associated with poor prognosis in breast cancer and HER2 has been defined as a therapeutic target for breast cancer treatment. We aimed to explore the molecular biological information in ultrasound radiomic features (URFs) of HER2-positive breast cancer using radiogenomic analysis. Moreover, a radiomics model was developed to predict the status of HER2 in breast cancer.

METHODS

This retrospective study included 489 patients who were diagnosed with breast cancer. URFs were extracted from a radiomics analysis set using PyRadiomics. The correlations between differential URFs and HER2-related genes were calculated using Pearson correlation analysis. Functional enrichment of the identified URFs-correlated HER2 positive-specific genes was performed. Lastly, the radiomics model was developed based on the URF-module mined from auxiliary differential URFs to assess the HER2 status of breast cancer.

RESULTS

Eight differential URFs (p < 0.05) were identified among the 86 URFs extracted by Pyradiomics. 25 genes that were found to be the most closely associated with URFs. Then, the relevant biological functions of each differential URF were obtained through functional enrichment analysis. Among them, Zone Entropy is related to immune cell activity, which regulate the generation of calcification in breast cancer. The radiomics model based on the Logistic classifier and URF-module showed good discriminative ability (AUC = 0.80, 95% CI).

CONCLUSION

We searched for the URFs of HER2-positive breast cancer, and explored the underlying genes and biological functions of these URFs. Furthermore, the radiomics model based on the Logistic classifier and URF-module relatively accurately predicted the HER2 status in breast cancer.

Radiogenomic association of deep MR imaging features with genomic profiles and clinical characteristics in breast cancer

BACKGROUND

It has been believed that traditional handcrafted radiomic features extracted from magnetic resonance imaging (MRI) of tumors are normally shallow and low-ordered. Recent advancement in deep learning technology shows that the high-order deep radiomic features extracted automatically from tumor images can capture tumor heterogeneity in a more efficient way. We hypothesize that MRI-based deep radiomic phenotypes have significant associations with molecular profiles of breast cancer tumors. We aim to identify deep radiomic features (DRFs) from MRI, evaluate their significance in predicting breast cancer (BC) clinical characteristics and explore their associations with multi-level genomic factors.

METHODS

A denoising autoencoder was built to retrospectively extract 4,096 DRFs from 110 BC patients' MRI. Visualization and clustering were applied to these DRFs. Linear Mixed Effect models were used to test their associations with multi-level genomic features (GFs) (risk genes, gene signatures, and biological pathway activities) extracted from the same patients' mRNA expression profile. A Least Absolute Shrinkage and Selection Operator model was used to identify the most predictive DRFs for each clinical characteristic (tumor size (T), lymph node metastasis (N), estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) status).

RESULTS

Thirty-six conventional radiomic features (CRFs) for 87 of the 110 BC patients provided by a previous study were used for comparison. More than 1,000 DRFs were associated with the risk genes, gene signatures, and biological pathways activities (adjusted P-value < 0.05). DRFs produced better performance in predicting T, N, ER, PR, and HER2 status (AUC > 0.9) using DRFs. These DRFs showed significant powers of stratifying patients, linking to relevant biological and clinical characteristics. As a contrast, only eight risk genes were associated with CRFs. The RFs performed worse in predicting clinical characteristics than DRFs.

CONCLUSIONS

The deep learning-based auto MRI features perform better in predicting BC clinical characteristics, which are more significantly associated with GFs than traditional semi-auto MRI features. Our radiogenomic approach for identifying MRI-based imaging signatures may pave potential pathways for the discovery of genetic mechanisms regulating specific tumor phenotypes and may enable a more rapid innovation of novel imaging modalities, hence accelerating their translation to personalized medicine.

Lung Cancer Subtype Diagnosis by Fusing Image-Genomics Data and Hybrid Deep Networks

Accurate diagnosis of cancer subtypes is crucial for precise treatment, because different cancer subtypes are involved with different pathology and require different therapies. Although deep learning techniques have made great success in computer vision and other fields, they do not work well on Lung cancer subtype diagnosis, due to the distinction of slide images between different cancer subtypes is ambiguous. Furthermore, they often over-fit to high-dimensional genomics data with limited samples, and do not fuse the image and genomics data in a sensible way. In this paper, we propose a hybrid deep network based approach LungDIG for Lung cancer subtype Diagnosis by fusing Image-Genomics data. LungDIG first tiles the tissue slide image into small patches and extracts the patch-level features by fine-tuning an Inception-V3 model. Since the patches may contain some false positives in non-diagnostic regions, it further designs a patch-level feature combination strategy to integrate the extracted patch features and maintain the diversity between different cancer subtypes. At the same time, it extracts the genomics features from Copy Number Variation data by an attention based nonlinear extractor. Next, it fuses the image and genomics features by an attention based multilayer perceptron (MLP) to diagnose cancer subtype. Experiments on TCGA lung cancer data show that LungDIG can not only achieve higher accuracy for cancer subtype diagnosis than state-of-the-art methods, but also have a high authenticity and good interpretability.

A review of radiomics and genomics applications in cancers: the way towards precision medicine

The application of radiogenomics in oncology has great prospects in precision medicine. Radiogenomics combines large volumes of radiomic features from medical digital images, genetic data from high-throughput sequencing, and clinical-epidemiological data into mathematical modelling. The amalgamation of radiomics and genomics provides an approach to better study the molecular mechanism of tumour pathogenesis, as well as new evidence-supporting strategies to identify the characteristics of cancer patients, make clinical decisions by predicting prognosis, and improve the development of individualized treatment guidance. In this review, we summarized recent research on radiogenomics applications in solid cancers and presented the challenges impeding the adoption of radiomics in clinical practice. More standard guidelines are required to normalize radiomics into reproducible and convincible analyses and develop it as a mature field.

Identifying Associations between DCE-MRI Radiomic Features and Expression Heterogeneity of Hallmark Pathways in Breast Cancer: A Multi-Center Radiogenomic Study

BACKGROUND

To investigate the relationship between dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) radiomic features and the expression activity of hallmark pathways and to develop prediction models of pathway-level heterogeneity for breast cancer (BC) patients.

METHODS

Two radiogenomic cohorts were analyzed (n = 246). Tumor regions were segmented semiautomatically, and 174 imaging features were extracted. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were performed to identify significant imaging-pathway associations. Random forest regression was used to predict pathway enrichment scores. Five-fold cross-validation and grid search were used to determine the optimal preprocessing operation and hyperparameters.

RESULTS

We identified 43 pathways, and 101 radiomic features were significantly related in the discovery cohort (p-value < 0.05). The imaging features of the tumor shape and mid-to-late post-contrast stages showed more transcriptional connections. Ten pathways relevant to functions such as cell cycle showed a high correlation with imaging in both cohorts. The prediction model for the mTORC1 signaling pathway achieved the best performance with the mean absolute errors (MAEs) of 27.29 and 28.61% in internal and external test sets, respectively.

CONCLUSIONS

The DCE-MRI features were associated with hallmark activities and may improve individualized medicine for BC by noninvasively predicting pathway-level heterogeneity.

CT-based radiogenomic analysis dissects intratumor heterogeneity and predicts prognosis of colorectal cancer: a multi-institutional retrospective study

BACKGROUND

This study aimed to develop a radiogenomic prognostic prediction model for colorectal cancer (CRC) by investigating the biological and clinical relevance of intratumoural heterogeneity.

METHODS

This retrospective multi-cohort study was conducted in three steps. First, we identified genomic subclones using unsupervised deconvolution analysis. Second, we established radiogenomic signatures to link radiomic features with prognostic subclone compositions in an independent radiogenomic dataset containing matched imaging and gene expression data. Finally, the prognostic value of the identified radiogenomic signatures was validated using two testing datasets containing imaging and survival information collected from separate medical centres.

RESULTS

This multi-institutional retrospective study included 1601 patients (714 females and 887 males; mean age, 65 years ± 14 [standard deviation]) with CRC from 5 datasets. Molecular heterogeneity was identified using unsupervised deconvolution analysis of gene expression data. The relative prevalence of the two subclones associated with cell cycle and extracellular matrix pathways identified patients with significantly different survival outcomes. A radiogenomic signature-based predictive model significantly stratified patients into high- and low-risk groups with disparate disease-free survival (HR = 1.74, P = 0.003). Radiogenomic signatures were revealed as an independent predictive factor for CRC by multivariable analysis (HR = 1.59, 95% CI:1.03-2.45, P = 0.034). Functional analysis demonstrated that the 11 radiogenomic signatures were predominantly associated with extracellular matrix and immune-related pathways.

CONCLUSIONS

The identified radiogenomic signatures might be a surrogate for genomic signatures and could complement the current prognostic strategies.

The dynamic dysregulated network identifies stage-specific markers during lung adenocarcinoma malignant progression and metastasis

Brain metastasis occurs in approximately 30% of patients with lung adenocarcinoma (LUAD) and is closely associated with poor prognosis, recurrence, and death. However, dynamic gene regulation and molecular mechanism driving LUAD progression remain poorly understood. In this study, we performed a comprehensive single-cell transcriptome analysis using data from normal, early stage, advanced stage, and brain metastasis LUAD. Our single-cell-level analysis reveals the cellular composition heterogeneity at different stages during LUAD progression. We identified stage-specific risk genes that could contribute to LUAD progression and metastasis by reprogramming immune-related and metabolic-related functions. We constructed an early advanced metastatic dysregulated network and revealed the dynamic changes in gene regulations during LUAD progression. We identified 6 early advanced (HLA-DRB1, HLA-DQB1, SFTPB, SFTPC, PLA2G1B, and FOLR1), 8 advanced metastasis (RPS15, RPS11, RPL13A, RPS24, HLA-DRB5, LYPLA1, KCNJ15, and PSMA3), and 2 common risk genes in different stages (SFTPD and HLA-DRA) as prognostic markers in LUAD. Particularly, decreased expression of HLA-DRA, HLA-DRB1, HLA-DQB1, and HLA-DRB5 refer poor prognosis in LUAD by controlling antigen processing and presentation and T cell activation. Increased expression of PSMA3 and LYPLA1 refer poor prognosis by reprogramming fatty acid metabolism and RNA catabolic process. Our findings will help further understanding the pathobiology of brain metastases in LUAD.

Unsupervised Analysis Based on DCE-MRI Radiomics Features Revealed Three Novel Breast Cancer Subtypes with Distinct Clinical Outcomes and Biological Characteristics

Background: This study aimed to reveal the heterogeneity of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) of breast cancer (BC) and identify its prognosis values and molecular characteristics. Methods: Two radiogenomics cohorts (n = 246) were collected and tumor regions were segmented semi-automatically. A total of 174 radiomics features were extracted, and the imaging subtypes were identified and validated by unsupervised analysis. A gene-profile-based classifier was developed to predict the imaging subtypes. The prognostic differences and the biological and microenvironment characteristics of subtypes were uncovered by bioinformatics analysis. Results: Three imaging subtypes were identified and showed high reproducibility. The subtypes differed remarkably in tumor sizes and enhancement patterns, exhibiting significantly different disease-free survival (DFS) or overall survival (OS) in the discovery cohort (p = 0.024) and prognosis datasets (p ranged from <0.0001 to 0.0071). Large sizes and rapidly enhanced tumors usually had the worst outcomes. Associations were found between imaging subtypes and the established subtypes or clinical stages (p ranged from <0.001 to 0.011). Imaging subtypes were distinct in cell cycle and extracellular matrix (ECM)-receptor interaction pathways (false discovery rate, FDR < 0.25) and different in cellular fractions, such as cancer-associated fibroblasts (p < 0.05). Conclusions: The imaging subtypes had different clinical outcomes and biological characteristics, which may serve as potential biomarkers.

Predicting hormone receptors and PAM50 subtypes of breast cancer from multi-scale lesion images of DCE-MRI with transfer learning technique

BACKGROUND

The recent development of artificial intelligence (AI) technologies coupled with medical imaging data has gained considerable attention, and offers a non-invasive approach for cancer diagnosis and prognosis. In this context, improved breast cancer (BC) molecular characteristics assessment models are foreseen to enable personalized strategies with better clinical outcomes compared to existing screening strategies. And it is a promising approach to developing models for hormone receptors (HR) and subtypes of BC patients from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) data.

METHODS

In this institutional review board-approved study, 174 BC patients with both DCE-MRI and RNA-seq data in the local database were analyzed. Slice images from tumor lesions and multi-scale peri-tumor regions were used as model inputs, and five representative pre-trained transfer learning (TF) networks, such as Inception-v3 and Xception, were employed to establish prediction models. A comprehensive analysis was performed using five-fold cross-validation to avoid overfitting, and accuracy (ACC) and area under the receiver operating characteristic curve (AUROC) to evaluate model performance.

RESULTS

Xception achieved the superior results when using solely tumor regions, with highest AUROCs of 0.844 (95% CI: [0.841, 0.847]) and 0.784 (95% CI: [0.781, 0.788]) for estrogen receptor (ER) and progesterone receptor (PR), respectively, and best ACC of 0.467 (95% CI: [0.462, 0.470]) for PAM50 subtypes. A significant improvement in the model performance was observed when images of the peri-tumor region were included, with optimal results achieved using images of the tumor and the 10 mm peri-tumor regions. Xception-based TF models performed most effectively in predicting ER and PR statuses, with the AUROCs were 0.942 (95% CI: [0.940, 0.944]) and 0.920 (95% CI: [0.917, 0.922]), respectively, whereas for PAM50 subtypes, the Inception-v3-based network yielded the highest ACC as 0.742 (95% CI: [0.738, 0.746]).

CONCLUSIONS

Transfer learning analysis based on DCE-MRI data of tumor and peri-tumor regions was helpful to the non-invasive assessment of molecular characteristics of BC.

Integrative Histology-Genomic Analysis Predicts Hepatocellular Carcinoma Prognosis Using Deep Learning

Cancer prognosis analysis is of essential interest in clinical practice. In order to explore the prognostic power of computational histopathology and genomics, this paper constructs a multi-modality prognostic model for survival prediction. We collected 346 patients diagnosed with hepatocellular carcinoma (HCC) from The Cancer Genome Atlas (TCGA), each patient has 1-3 whole slide images (WSIs) and an mRNA expression file. WSIs were processed by a multi-instance deep learning model to obtain the patient-level survival risk scores; mRNA expression data were processed by weighted gene co-expression network analysis (WGCNA), and the top hub genes of each module were extracted as risk factors. Information from two modalities was integrated by Cox proportional hazard model to predict patient outcomes. The overall survival predictions of the multi-modality model (Concordance index (C-index): 0.746, 95% confidence interval (CI): ±0.077) outperformed these based on histopathology risk score or hub genes, respectively. Furthermore, in the prediction of 1-year and 3-year survival, the area under curve of the model achieved 0.816 and 0.810. In conclusion, this paper provides an effective workflow for multi-modality prognosis of HCC, the integration of histopathology and genomic information has the potential to assist clinical prognosis management.

A hybrid hemodynamic knowledge-powered and feature reconstruction-guided scheme for breast cancer segmentation based on DCE-MRI

Automatically and accurately annotating tumor in dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), which provides a noninvasive in vivo method to evaluate tumor vasculature architectures based on contrast accumulation and washout, is a crucial step in computer-aided breast cancer diagnosis and treatment. However, it remains challenging due to the varying sizes, shapes, appearances and densities of tumors caused by the high heterogeneity of breast cancer, and the high dimensionality and ill-posed artifacts of DCE-MRI. In this paper, we propose a hybrid hemodynamic knowledge-powered and feature reconstruction-guided scheme that integrates pharmacokinetics prior and feature refinement to generate sufficiently adequate features in DCE-MRI for breast cancer segmentation. The pharmacokinetics prior expressed by time intensity curve (TIC) is incorporated into the scheme through objective function called dynamic contrast-enhanced prior (DCP) loss. It contains contrast agent kinetic heterogeneity prior knowledge, which is important to optimize our model parameters. Besides, we design a spatial fusion module (SFM) embedded in the scheme to exploit intra-slices spatial structural correlations, and deploy a spatial-kinetic fusion module (SKFM) to effectively leverage the complementary information extracted from spatial-kinetic space. Furthermore, considering that low spatial resolution often leads to poor image quality in DCE-MRI, we integrate a reconstruction autoencoder into the scheme to refine feature maps in an unsupervised manner. We conduct extensive experiments to validate the proposed method and show that our approach can outperform recent state-of-the-art segmentation methods on breast cancer DCE-MRI dataset. Moreover, to explore the generalization for other segmentation tasks on dynamic imaging, we also extend the proposed method to brain segmentation in DSC-MRI sequence. Our source code will be released on https://github.com/AI-medical-diagnosis-team-of-JNU/DCEDuDoFNet.

Radiogenomics analysis reveals the associations of dynamic contrast-enhanced-MRI features with gene expression characteristics, PAM50 subtypes, and prognosis of breast cancer

BACKGROUND

To investigate reliable associations between dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) features and gene expression characteristics in breast cancer (BC) and to develop and validate classifiers for predicting PAM50 subtypes and prognosis from DCE-MRI non-invasively.

METHODS

Two radiogenomics cohorts with paired DCE-MRI and RNA-sequencing (RNA-seq) data were collected from local and public databases and divided into discovery (n = 174) and validation cohorts (n = 72). Six external datasets (n = 1,443) were used for prognostic validation. Spatial-temporal features of DCE-MRI were extracted, normalized properly, and associated with gene expression to identify the imaging features that can indicate subtypes and prognosis.

RESULTS

Expression of genes including RBP4, MYBL2, and LINC00993 correlated significantly with DCE-MRI features (q-value < 0.05). Importantly, genes in the cell cycle pathway exhibited a significant association with imaging features (p-value < 0.001). With eight imaging-associated genes (CHEK1, TTK, CDC45, BUB1B, PLK1, E2F1, CDC20, and CDC25A), we developed a radiogenomics prognostic signature that can distinguish BC outcomes in multiple datasets well. High expression of the signature indicated a poor prognosis (p-values < 0.01). Based on DCE-MRI features, we established classifiers to predict BC clinical receptors, PAM50 subtypes, and prognostic gene sets. The imaging-based machine learning classifiers performed well in the independent dataset (areas under the receiver operating characteristic curve (AUCs) of 0.8361, 0.809, 0.7742, and 0.7277 for estrogen receptor (ER), human epidermal growth factor receptor 2 (HER2)-enriched, basal-like, and obtained radiogenomics signature). Furthermore, we developed a prognostic model directly using DCE-MRI features (p-value < 0.0001).

CONCLUSIONS

Our results identified the DCE-MRI features that are robust and associated with the gene expression in BC and displayed the possibility of using the features to predict clinical receptors and PAM50 subtypes and to indicate BC prognosis.

A framework for deep multitask learning with multiparametric magnetic resonance imaging for the joint prediction of histological characteristics in breast cancer

The clinical management and decision-making process related to breast cancer are based on multiple histological indicators. This study aims to jointly predict the Ki-67 expression level, luminal A subtype and histological grade molecular biomarkers using a new deep multitask learning method with multiparametric magnetic resonance imaging. A multitask learning network structure was proposed by introducing a common-task layer and task-specific layers to learn the high-level features that are common to all tasks and related to a specific task, respectively. A network pretrained with knowledge from the ImageNet dataset was used and fine-tuned with MRI data. Information from multiparametric MR images was fused using the strategy at the feature and decision levels. The area under the receiver operating characteristic curve (AUC) was used to measure model performance. For single-task learning using a single image series, the deep learning model generated AUCs of 0.752, 0.722, and 0.596 for the Ki-67, luminal A and histological grade prediction tasks, respectively. The performance was improved by freezing the first 5 convolutional layers, using 20% shared layers and fusing multiparametric series at the feature level, which achieved AUCs of 0.819, 0.799 and 0.747 for Ki-67, luminal A and histological grade prediction tasks, respectively. Our study showed advantages in jointly predicting correlated clinical biomarkers using a deep multitask learning framework with an appropriate number of fine-tuned convolutional layers by taking full advantage of common and complementary imaging features. Multiparametric image series-based multitask learning could be a promising approach for the multiple clinical indicator-based management of breast cancer.

A novel integrative computational framework for breast cancer radiogenomic biomarker discovery

•

Bayesian tensor factorization is used to integrate multiomics data for radiogenomics analysis.

•

A regression framework is proposed to handle the unmatched data issue in radiogenomics analysis.

•

Deep learning is used to identify prognostic meaningful radiogenomic biomarkers for cancer.

In precise medicine, it is with great value to develop computational frameworks for identifying prognostic biomarkers which can capture both multi-genomic and phenotypic heterogeneity of breast cancer (BC). Radiogenomics is a field where medical images and genomic measurements are integrated and mined to solve challenging clinical problems. Previous radiogenomic studies suffered from data incompleteness, feature subjectivity and low interpretability. For example, the majority of the radiogenomic studies miss one or two of medical imaging data, genomic data, and clinical outcome data, which results in the data incomplete issue. Feature subjectivity issue comes from the extraction of imaging features with significant human involvement. Thus, there is an urgent need to address above-mentioned limitations so that fully automatic and transparent radiogenomic prognostic biomarkers could be identified for BC.

We proposed a novel framework for BC prognostic radiogenomic biomarker identification. This framework involves an explainable DL model for image feature extraction, a Bayesian tensor factorization (BTF) processing for multi-genomic feature extraction, a leverage strategy to utilize unpaired imaging, genomic, and survival outcome data, and a mediation analysis to provide further interpretation for identified biomarkers. This work provided a new perspective for conducting a comprehensive radiogenomic study when only limited resources are given. Compared with baseline traditional radiogenomic biomarkers, the 23 biomarkers identified by the proposed framework performed better in indicating patients’ survival outcome. And their interpretability is guaranteed by different levels of build-in and follow-up analyses.

A Multi-Classification Model for Predicting the Invasiveness of Lung Adenocarcinoma Presenting as Pure Ground-Glass Nodules

Objectives

To establish a multi-classification model for precisely predicting the invasiveness (pre-invasive adenocarcinoma, PIA; minimally invasive adenocarcinoma, MIA; invasive adenocarcinoma, IAC) of lung adenocarcinoma manifesting as pure ground-glass nodules (pGGNs).

Methods

By the inclusion and exclusion criteria, this retrospective study enrolled 346 patients (female, 297, and male, 49; age, 55.79 ± 10.53 (24-83)) presenting as pGGNs from 1292 consecutive patients with pathologically confirmed lung adenocarcinoma. A total of 27 clinical were collected and 1409 radiomics features were extracted by PyRadiomics package on python. After feature selection with L2,1-norm minimization, logistic regression (LR), extra w(ET) and gradient boosting decision tree (GBDT) were used to construct the three-classification model. Then, an ensemble model of the three algorithms based on model ensemble strategy was established to further improve the classification performance.

Results

After feature selection, a hybrid of 166 features consisting of 1 clinical (short-axis diameter, ranked 27th) and 165 radiomics (4 shape, 71 intensity and 90 texture) features were selected. The three most important features are wavelet-HLL_firstorder_Minimum, wavelet-HLL_ngtdm_Busyness and square_firstorder_Kurtosis. The hybrid-ensemble model based on hybrid clinical-radiomics features and the ensemble strategy showed more accurate predictive performance than other models (hybrid-LR, hybrid-ET, hybrid-GBDT, clinical-ensemble and radiomics-ensemble). On the training set and test set, the model can obtain the accuracy values of 0.918 ± 0.022 and 0.841, and its F1-scores respectively were 0.917 ± 0.024 and 0.824.

Conclusion

The multi-classification of invasive pGGNs can be precisely predicted by our proposed hybrid-ensemble model to assist patients in the early diagnosis of lung adenocarcinoma and prognosis.

Deep learning signatures reveal multiscale intratumor heterogeneity associated with biological functions and survival in recurrent nasopharyngeal carcinoma

PURPOSE

How to discriminate different risks of recurrent nasopharyngeal carcinoma (rNPC) patients and guide individual treatment has become of great importance. This study aimed to explore the associations between deep learning signatures and biological functions as well as survival in (rNPC) patients.

METHODS

A total of 420 rNPC patients with PET/CT imaging and follow-up of overall survival (OS) were retrospectively enrolled. All patients were randomly divided into a training set (n = 269) and test set (n = 151) with a 6:4 ratio. We constructed multi-modality deep learning signatures from PET and CT images with a light-weighted deep convolutional neural network EfficienetNet-lite0 and survival loss DeepSurvLoss. An integrated nomogram was constructed incorporating clinical factors and deep learning signatures from PET/CT. Clinical nomogram and single-modality deep learning nomograms were also built for comparison. Furthermore, the association between biological functions and survival risks generated from an integrated nomogram was analyzed by RNA sequencing (RNA-seq).

RESULTS

The C-index of the integrated nomogram incorporating age, rT-stage, and deep learning PET/CT signature was 0.741 (95% CI: 0.688-0.794) in the training set and 0.732 (95% CI: 0.679-0.785) in the test set. The nomogram stratified patients into two groups with high risk and low risk in both the training set and test set with hazard ratios (HR) of 4.56 (95% CI: 2.80-7.42, p < 0.001) and 4.05 (95% CI: 2.21-7.43, p < 0.001), respectively. The C-index of the integrated nomogram was significantly higher than the clinical nomogram and single-modality nomograms. When stratified by sex, N-stage, or EBV DNA, risk prediction of our integrated nomogram was valid in all patient subgroups. Further subgroup analysis showed that patients with a low-risk could benefit from surgery and re-irradiation, while there was no difference in survival rates between patients treated by chemotherapy in the high-risk and low-risk groups. RNA sequencing (RNA-seq) of data further explored the mechanism of high- and low-risk patients from the genetic and molecular level.

CONCLUSION

Our study demonstrated that PET/CT-based deep learning signatures showed satisfactory prognostic predictive performance in rNPC patients. The nomogram incorporating deep learning signatures successfully divided patients into different risks and had great potential to guide individual treatment: patients with a low-risk were supposed to be treated with surgery and re-irradiation, while for high-risk patients, the application of palliative chemotherapy may be sufficient.

Radiomics predicts the prognosis of patients with locally advanced breast cancer by reflecting the heterogeneity of tumor cells and the tumor microenvironment

Background

This study investigated the efficacy of radiomics to predict survival outcome for locally advanced breast cancer (LABC) patients and the association of radiomics with tumor heterogeneity and microenvironment.

Methods

Patients with LABC from 2010 to 2015 were retrospectively reviewed. Radiomics features were extracted from enhanced MRI. We constructed the radiomics score using lasso and assessed its prognostic value. An external validation cohort from The Cancer Imaging Archive was used to assess phenotype reproducibility. Sequencing data from TCGA and our center were applied to reveal genomic landscape of different radiomics score groups. Tumor infiltrating lymphocytes map and bioinformatics methods were applied to evaluate the heterogeneity of tumor microenvironment. Computational histopathology was also applied.

Results

A total of 278 patients were divided into training cohort and validation cohort. Radiomics score was constructed and significantly associated with disease-free survival (DFS) of the patients in training cohort, validation cohort and external validation cohort (p < 0.001, p = 0.014 and p = 0.041, respectively). The radiomics-based nomogram showed better predictive performance of DFS compared with TNM model. Distinct gene expression patterns were identified. Immunophenotype and immune cell composition was different in each radiomics score group. The link between radiomics and computational histopathology was revealed.

Conclusions

The radiomics score could effectively predict prognosis of LABC after neoadjuvant chemotherapy and radiotherapy. Radiomics revealed heterogeneity of tumor cell and tumor microenvironment and holds great potential to facilitate individualized DFS estimation and guide personalized care.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-022-01516-0.

swCAM: estimation of subtype-specific expressions in individual samples with unsupervised sample-wise deconvolution

MOTIVATION

Complex biological tissues are often a heterogeneous mixture of several molecularly distinct cell subtypes. Both subtype compositions and subtype-specific (STS) expressions can vary across biological conditions. Computational deconvolution aims to dissect patterns of bulk tissue data into subtype compositions and STS expressions. Existing deconvolution methods can only estimate averaged STS expressions in a population, while many downstream analyses such as inferring co-expression networks in particular subtypes require subtype expression estimates in individual samples. However, individual-level deconvolution is a mathematically underdetermined problem because there are more variables than observations.

RESULTS

We report a sample-wise Convex Analysis of Mixtures (swCAM) method that can estimate subtype proportions and STS expressions in individual samples from bulk tissue transcriptomes. We extend our previous CAM framework to include a new term accounting for between-sample variations and formulate swCAM as a nuclear-norm and ℓ2,1-norm regularized matrix factorization problem. We determine hyperparameter values using cross-validation with random entry exclusion and obtain a swCAM solution using an efficient alternating direction method of multipliers. Experimental results on realistic simulation data show that swCAM can accurately estimate STS expressions in individual samples and successfully extract co-expression networks in particular subtypes that are otherwise unobtainable using bulk data. In two real-world applications, swCAM analysis of bulk RNASeq data from brain tissue of cases and controls with bipolar disorder or Alzheimer's disease identified significant changes in cell proportion, expression pattern and co-expression module in patient neurons. Comparative evaluation of swCAM versus peer methods is also provided.

AVAILABILITY AND IMPLEMENTATION

The R Scripts of swCAM are freely available at https://github.com/Lululuella/swCAM. A user's guide and a vignette are provided.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A deep matrix factorization framework for identifying underlying tissue-specific patterns of DCE-MRI: applications for molecular subtype classification in breast cancer

Objective.Breast cancer is heterogeneous in that different angiogenesis and blood flow characteristics could be present within a tumor. The pixel kinetics of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) can assume several distinct signal patterns related to specific tissue characteristics. Identification of the latent, tissue-specific dynamic patterns of intratumor heterogeneity can shed light on the biological mechanisms underlying the heterogeneity of tumors.Approach.To mine this information, we propose a deep matrix factorization-based dynamic decomposition (DMFDE) model specifically designed according to DCE-MRI characteristics. The time-series imaging data were decomposed into tissue-specific dynamic patterns and their corresponding proportion maps. The image pixel matrix and the reference matrix of population-level kinetics obtained by clustering the dynamic signals were used as the inputs. Two multilayer neural network branches were designed to collaboratively project the input matrix into a latent dynamic pattern and a dynamic proportion matrix, which was justified using simulated data. Clinical implications of DMFDE were assessed by radiomics analysis of proportion maps obtained from the tumor/parenchyma region for classifying the luminal A subtype.Main results.The decomposition performance of DMFDE was evaluated by the root mean square error and was shown to be better than that of the conventional convex analysis of mixtures (CAM) method. The predictive model withK = 3, 4, and 5 dynamic proportion maps generated AUC values of 0.780, 0.786 and 0.790, respectively, in distinguishing between luminal A and nonluminal A tumors, which are better than the CAM method (AUC = 0.726). The combination of statistical features from images with different proportion maps has the highest prediction value (AUC = 0.813), which is significantly higher than that based on CAM.Conclusion.This proposed method identified the latent dynamic patterns associated with different molecular subtypes, and radiomics analysis based on the pixel compositions of the uncovered dynamic patterns was able to determine molecular subtypes of breast cancer.

A radiogenomics application for prognostic profiling of endometrial cancer

Prognostication is critical for accurate diagnosis and tailored treatment in endometrial cancer (EC). We employed radiogenomics to integrate preoperative magnetic resonance imaging (MRI, n = 487 patients) with histologic-, transcriptomic- and molecular biomarkers (n = 550 patients) aiming to identify aggressive tumor features in a study including 866 EC patients. Whole-volume tumor radiomic profiling from manually (radiologists) segmented tumors (n = 138 patients) yielded clusters identifying patients with high-risk histological features and poor survival. Radiomic profiling by a fully automated machine learning (ML)-based tumor segmentation algorithm (n = 336 patients) reproduced the same radiomic prognostic groups. From these radiomic risk-groups, an 11-gene high-risk signature was defined, and its prognostic role was reproduced in orthologous validation cohorts (n = 554 patients) and aligned with The Cancer Genome Atlas (TCGA) molecular class with poor survival (copy-number-high/p53-altered). We conclude that MRI-based integrated radiogenomics profiling provides refined tumor characterization that may aid in prognostication and guide future treatment strategies in EC.

Radiogenomic Signatures of Oncotype DX Recurrence Score Enable Prediction of Survival in Estrogen Receptor-Positive Breast Cancer: A Multicohort Study

Background Radiogenomics explores the association between imaging features and genomic assays to uncover relevant prognostic features; however, the prognostic implications of the derived signatures remain unclear. Purpose To identify preoperative radiogenomic signatures of estrogen receptor-positive breast cancer associated with the Oncotype DX recurrence score (RS) and to evaluate whether they are biomarkers for survival and responses to neoadjuvant chemotherapy (NACT). Materials and Methods In this retrospective multicohort study, three data sets were analyzed. The radiogenomic development data set, with preoperative dynamic contrast-enhanced MRI and RS data obtained between January 2016 and October 2019 was used to identify radiogenomic signatures. Prognostic implications of the imaging signatures were assessed by measuring overall survival and recurrence-free survival in the prognostic assessment data set using a multivariable Cox proportional hazards model. The therapeutic implication of the radiogenomic signatures was evaluated by determining their ability to predict the response to NACT using the treatment assessment data set obtained between August 2015 and March 2019. Prediction performance was estimated by using the area under the receiver operating characteristic curve (AUC). Results The final cohorts included a radiogenomic development data set with 130 women (mean age, 52 years ± 10 [standard deviation]), a prognostic assessment data set with 116 women (mean age, 48 years ± 9), and a treatment assessment data set with 135 women (mean age, 50 years ± 11). Radiogenomic signatures (n = 11) of texture and morphologic and statistical features were identified to generate the predicted RS (R² = 0.33, P < .001). A predicted RS greater than 29.9 was associated with poor overall and recurrence-free survival (P = .001 and P = .007, respectively); predicted RS was greater in women with a good NACT response (30.51 ± 6.92 vs 27.35 ± 4.04 [responders vs nonresponders], P = .001). By combining the predicted RS and complementary features, the model achieved improved performance in prediction of the NACT response (AUC, 0.85; P < .001). Conclusion Radiogenomic signatures associated with genomic assays provide markers of prognosis and treatment in estrogen receptor-positive breast cancer. © RSNA, 2021 Online supplemental material is available for this article.

Radiomic Analysis of Pharmacokinetic Heterogeneity Within Tumor Based on the Unsupervised Decomposition of Dynamic Contrast-Enhanced MRI for Predicting Histological Characteristics of Breast Cancer

BACKGROUND

Breast tumor heterogeneity is associated with histological characteristics. However, pharmacokinetic (PK) heterogeneity within tumor might merit further exploration.

PURPOSE

To enhance the predictive power of molecular subtypes, Ki-67, and tumor grade by analyzing PK heterogeneity within tumor based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

STUDY TYPE

Retrospective.

POPULATION

Two hundred and eight biopsy-proven breast cancer patients, randomly divided into a training cohort (N = 144) and a testing cohort (N = 64).

FIELD STRENGTH/SEQUENCE

T₁ -weighted DCE-MRI at 3.0 T.

ASSESSMENT

A convex analysis of mixtures-compartmental modeling decomposition method was used to estimate the PK parameter (i.e., the volume transfer constant K^trans ) in tumor subregions with distinct physiological kinetic patterns, including fast-flow kinetics, slow-flow kinetics, and plasma input. Radiomic features based on the PK parameter were calculated from each tumor subregion.

STATISTICAL TESTS

The training cohort was used to build random forest classifiers based on the optimal features determined by the 5-fold cross-validation method. The performance was assessed on the testing cohort using the area under the receiver operating characteristic curve (AUC). The AUCs derived from the tumor subregion-based PK parameter were compared with those of the original images of the entire tumor using the DeLong test. A P-value of <0.05 was considered statistically significant.

RESULTS

The tumor subregion-based PK parameter, which yielded the highest AUCs of 0.8782, 0.7568, 0.7019, 0.7963, 0.8080, and 0.7375 for luminal A, luminal B, basal-like, human epidermal growth factor receptor 2, Ki-67, and tumor grade, respectively, obtained better diagnostic performance than the original images in the entire tumor (highest AUCs = 0.8612, 0.6191, 0.5593, 0.7704, 0.7494, and 0.6261, respectively). In particular, statistically significant improvement in the diagnostic performance was obtained for luminal B.

DATA CONCLUSION

Radiomic analysis of PK heterogeneity within tumor can enhance the predictive performance of radiomic models compared with that of the entire tumor.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY STAGE: 3.