Breast Cancer: Radiogenomic Biomarker Reveals Associations among Dynamic Contrast-enhanced MR Imaging, Long Noncoding RNA, and Metastasis

Magnetic resonance imaging and molecular features associated with tumor-infiltrating lymphocytes in breast cancer

Background

We sought to investigate associations between dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) features and tumor-infiltrating lymphocytes (TILs) in breast cancer, as well as to study if MRI features are complementary to molecular markers of TILs.

Methods

In this retrospective study, we extracted 17 computational DCE-MRI features to characterize tumor and parenchyma in The Cancer Genome Atlas cohort (n = 126). The percentage of stromal TILs was evaluated on H&E-stained histological whole-tumor sections. We first evaluated associations between individual imaging features and TILs. Multiple-hypothesis testing was corrected by the Benjamini-Hochberg method using false discovery rate (FDR). Second, we implemented LASSO (least absolute shrinkage and selection operator) and linear regression nested with tenfold cross-validation to develop an imaging signature for TILs. Next, we built a composite prediction model for TILs by combining imaging signature with molecular features. Finally, we tested the prognostic significance of the TIL model in an independent cohort (I-SPY 1; n = 106).

Results

Four imaging features were significantly associated with TILs (P < 0.05 and FDR < 0.2), including tumor volume, cluster shade of signal enhancement ratio (SER), mean SER of tumor-surrounding background parenchymal enhancement (BPE), and proportion of BPE. Among molecular and clinicopathological factors, only cytolytic score was correlated with TILs (ρ = 0.51; 95% CI, 0.36–0.63; P = 1.6E-9). An imaging signature that linearly combines five features showed correlation with TILs (ρ = 0.40; 95% CI, 0.24–0.54; P = 4.2E-6). A composite model combining the imaging signature and cytolytic score improved correlation with TILs (ρ = 0.62; 95% CI, 0.50–0.72; P = 9.7E-15). The composite model successfully distinguished low vs high, intermediate vs high, and low vs intermediate TIL groups, with AUCs of 0.94, 0.76, and 0.79, respectively. During validation (I-SPY 1), the predicted TILs from the imaging signature separated patients into two groups with distinct recurrence-free survival (RFS), with log-rank P = 0.042 among triple-negative breast cancer (TNBC). The composite model further improved stratification of patients with distinct RFS (log-rank P = 0.0008), where TNBC with no/minimal TILs had a worse prognosis.

Conclusions

Specific MRI features of tumor and parenchyma are associated with TILs in breast cancer, and imaging may play an important role in the evaluation of TILs by providing key complementary information in equivocal cases or situations that are prone to sampling bias.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-1039-2) contains supplementary material, which is available to authorized users.

Relationships Between Human-Extracted MRI Tumor Phenotypes of Breast Cancer and Clinical Prognostic Indicators Including Receptor Status and Molecular Subtype

PURPOSE

The purpose of this study was to investigate if human-extracted MRI tumor phenotypes of breast cancer could predict receptor status and tumor molecular subtype using MRIs from The Cancer Genome Atlas project.

MATERIALS AND METHODS

Our retrospective interpretation study utilized the analysis of HIPAA-compliant breast MRI data from The Cancer Imaging Archive. One hundred and seven preoperative breast MRIs of biopsy proven invasive breast cancers were analyzed by 3 fellowship-trained breast-imaging radiologists. Each study was scored according to the Breast Imaging Reporting and Data System lexicon for mass and nonmass features. The Spearman rank correlation was used for association analysis of continuous variables; the Kruskal-Wallis test was used for associating continuous outcomes with categorical variables. The Fisher-exact test was used to assess correlations between categorical image-derived features and receptor status. Prediction of estrogen receptor (ER), progesterone receptor, human epidermal growth factor receptor, and molecular subtype were performed using random forest classifiers.

RESULTS

ER+ tumors were associated with the absence of rim enhancement (P = 0.019, odds ratio [OR] 5.5), heterogeneous internal enhancement (P = 0.02, OR 6.5), peritumoral edema (P = 0.0001, OR 10.0), and axillary adenopathy (P = 0.04, OR 4.4). ER+ tumors were smaller than ER- tumors (23.7 mm vs 29.2 mm, P = 0.02, OR 8.2). All of these variables except the lack of axillary adenopathy were also associated with progesterone receptor+ status. Luminal A tumors (n = 57) were smaller compared to nonLuminal A (21.8 mm vs 27.5 mm, P = 0.035, OR 7.3) and lacked peritumoral edema (P = 0.001, OR 6.8). Basal like tumors were associated with heterogeneous internal enhancement (P = 0.05, OR 10.1), rim enhancement (P = 0.05, OR6.9), and perituomral edema (P = 0.0001, OR 13.8).

CONCLUSIONS

Human extracted MRI tumor phenotypes may be able to differentiate those tumors with a more favorable clinical prognosis from their more aggressive counterparts.

Lung Cancer Radiogenomics: The Increasing Value of Imaging in Personalized Management of Lung Cancer Patients

Radiogenomics provide a large-scale data analytical framework that aims to understand the broad multiscale relationships between the complex information encoded in medical images (including computational, quantitative, and semantic image features) and their underlying clinical, therapeutic, and biological associations. As such it is a powerful and increasingly important tool for both clinicians and researchers involved in the imaging, evaluation, understanding, and management of lung cancers. Herein we provide an overview of the growing field of lung cancer radiogenomics and its applications.

Precision Medicine and Radiogenomics in Breast Cancer: New Approaches toward Diagnosis and Treatment

Precision medicine is medicine optimized to the genotypic and phenotypic characteristics of an individual and, when present, his or her disease. It has a host of targets, including genes and their transcripts, proteins, and metabolites. Studying precision medicine involves a systems biology approach that integrates mathematical modeling and biology genomics, transcriptomics, proteomics, and metabolomics. Moreover, precision medicine must consider not only the relatively static genetic codes of individuals, but also the dynamic and heterogeneous genetic codes of cancers. Thus, precision medicine relies not only on discovering identifiable targets for treatment and surveillance modification, but also on reliable, noninvasive methods of identifying changes in these targets over time. Imaging via radiomics and radiogenomics is poised for a central role. Radiomics, which extracts large volumes of quantitative data from digital images and amalgamates these together with clinical and patient data into searchable shared databases, potentiates radiogenomics, which is the combination of genetic and radiomic data. Radiogenomics may provide voxel-by-voxel genetic information for a complete, heterogeneous tumor or, in the setting of metastatic disease, set of tumors and thereby guide tailored therapy. Radiogenomics may also quantify lesion characteristics, to better differentiate between benign and malignant entities, and patient characteristics, to better stratify patients according to risk for disease, thereby allowing for more precise imaging and screening. This report provides an overview of precision medicine and discusses radiogenomics specifically in breast cancer. ^© RSNA, 2018.

Receptor-based Surrogate Subtypes and Discrepancies with Breast Cancer Intrinsic Subtypes: Implications for Image Biomarker Development

Purpose To determine the concordance and accuracy of imaging surrogates of immunohistochemical (IHC) markers and the molecular classification of breast cancer. Materials and Methods A total of 3050 patients from 17 public breast cancer data sets containing IHC marker receptor status (estrogen receptor/progesterone receptor/human epidermal growth factor receptor 2 [HER2]) and their molecular classification (basal-like, HER2-enriched, luminal A or B) were analyzed. Diagnostic accuracy and concordance as measured with the κ statistic were calculated between the IHC and molecular classifications. Simulations were performed to assess the relationship between accuracy of imaging-based IHC markers to predict molecular classification. A simulation was performed to examine effects of misclassification of molecular type on patient survival. Results Accuracies of intrinsic subtypes based on IHC subtype were 71.7% (luminal A), 53.7% (luminal B), 64.8% (HER2-enriched), and 81.7% (basal-like). The κ agreement was fair (κ = 0.36) for luminal A and HER2-enriched subtypes, good (κ = 0.65) for the basal-like subtype, and poor (κ = 0.09) for the luminal B subtypes. Introduction of image misclassification by simulation lowered image-true subtype accuracies and κ values. Simulation analysis showed that misclassification caused survival differences between luminal A and basal-like subtypes to decrease. Conclusion There is poor concordance between triple-receptor status and intrinsic molecular subtype in breast cancer, arguing against their use in the design of prognostic genomic-based image biomarkers. © RSNA, 2018 Online supplemental material is available for this article.

Neoadjuvant Chemotherapy and Surgery for Breast Cancer: Preoperative MRI Features Associated with Local Recurrence

Purpose To investigate the MRI and clinical-pathologic features associated with local-regional recurrence (LRR) in patients who had undergone breast-conserving surgery (BCS) following neoadjuvant chemotherapy (NAC). Materials and Methods In this retrospective, single-institution study between October 2003 and September 2015, 548 consecutive women, consisting of 468 down-staged and 80 preplanned BCS patients (mean age, 45.7 years; range, 22-75 years), underwent preoperative MRI and BCS following NAC. The rate and site of LRR, preoperative MRI features including Breast Imaging Reporting and Data System lexicon, and clinical-pathologic features (age, stage, tumor subtype, histologic grade, lymphovascular invasion, adjuvant chemotherapy, and endocrine therapy) were analyzed with the Cox proportional hazards model to identify independent factors associated with LRR-free survival (LRFS). Results Of the 548 women, 23 (4.2%) had LRR at a median follow-up of 23.1 months. In Cox regression analysis, younger age (ie, ≤ 40 years) (hazard ratio = 2.932; 95% confidence interval [CI]: 1.233, 6.969; P = .015) or the presence of nonmass enhancement on preoperative MR images (hazard ratio = 3.220; 95% CI: 1.274, 8.140; P = .014) was associated with worse LRFS. LRR was more frequently observed in the same quadrant as the original tumor in the down-staged BCS group than in the preplanned BCS group (80.0% [16 of 20] vs 33.3% [one of three]; P = .021). Conclusion Age of 40 years or younger or the presence of nonmass enhancement on preoperative MR images tends to be associated with worse local-regional recurrence-free survival, and local-regional recurrence frequently occurs in the same quadrant as the original tumor in breast cancer patients who undergo breast-conserving surgery following neoadjuvant chemotherapy. © RSNA, 2018.

Radiomics analysis of apparent diffusion coefficient in cervical cancer: A preliminary study on histological grade evaluation

BACKGROUND

The role of apparent diffusion coefficient (ADC)-based radiomics features in evaluating histopathological grade of cervical cancer is unresolved.

PURPOSE

To determine if there is a difference between radiomics features derived from center-slice 2D versus whole-tumor volumetric 3D for ADC measurements in patients with cervical cancer regarding tumor histopathological grade, and systematically assess the impact of the b value on radiomics analysis in ADC quantifications.

STUDY TYPE

Prospective.

SUBJECTS

In all, 160 patients with histopathologically confirmed squamous cell carcinoma of uterine cervix.

FIELD STRENGTH/SEQUENCE

Conventional and diffusion-weighted MR images (b values = 0, 800, 1000 s/mm² ) were acquired on a 3.0T MR scanner.

ASSESSMENT

Regions of interest (ROIs) were drawn manually along the margin of tumor on each slice, and then the center slice of the tumor was selected with naked eyes in the course of whole-tumor segmentation. A total of 624 radiomics features were derived from T₂ -weighted images and ADC maps. We randomly selected 50 cases and did the reproducibility analysis.

STATISTICAL TESTS

Parameters were compared using Wilcoxon signed rank test, Bland-Altman analysis, t-test, and least absolute shrinkage and selection operator (LASSO) regression with crossvalidation.

RESULTS

In all, 95 radiomics features were insensitive to ROI variation among T₂ images, ADC map of b800, and ADC map of b1000 (P > 0.0002). There was a significant statistical difference between the performances of 2D center-slice and 3D whole-tumor radiomics models in both ADC feature sets of b800 and b1000 (P < 0.0001, P < 0.0001). Compared with ADC features of b800 (0.3758 ± 0.0118), the model of b1000 ADC features appeared to be slightly lower in overall misclassification error (0.3642 ± 0.0162) (P = 0.0076).

DATA CONCLUSION

Several radiomics features extracted from T₂ images and ADC maps were highly reproducible. Whole-tumor volumetric 3D radiomics analysis had a better performance than using the 2D center-slice of tumor in stratifying the histological grade of cervical cancer. A b value of 1000 s/mm² is suggested as the optimal parameter in pelvic DWI scans.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:280-290.

Non-invasive tumor genotyping using radiogenomic biomarkers, a systematic review and oncology-wide pathway analysis

With targeted treatments playing an increasing role in oncology, the need arises for fast non-invasive genotyping in clinical practice. Radiogenomics is a rapidly evolving field of research aimed at identifying imaging biomarkers useful for non-invasive genotyping. Radiogenomic genotyping has the advantage that it can capture tumor heterogeneity, can be performed repeatedly for treatment monitoring, and can be performed in malignancies for which biopsy is not available. In this systematic review of 187 included articles, we compiled a database of radiogenomic associations and unraveled networks of imaging groups and gene pathways oncology-wide. Results indicated that ill-defined tumor margins and tumor heterogeneity can potentially be used as imaging biomarkers for 1p/19q codeletion in glioma, relevant for prognosis and disease profiling. In non-small cell lung cancer, FDG-PET uptake and CT-ground-glass-opacity features were associated with treatment-informing traits including EGFR-mutations and ALK-rearrangements. Oncology-wide gene pathway analysis revealed an association between contrast enhancement (imaging) and the targetable VEGF-signalling pathway. Although the need of independent validation remains a concern, radiogenomic biomarkers showed potential for prognosis prediction and targeted treatment selection. Quantitative imaging enhanced the potential of multiparametric radiogenomic models. A wealth of data has been compiled for guiding future research towards robust non-invasive genomic profiling.

Association of preoperative breast MRI features with locoregional recurrence after breast conservation therapy

Background Locoregional recurrence (LRR) following breast conservation therapy (BCT) is associated with an increased risk of distant metastasis and death in patients with breast cancer. Purpose To investigate whether preoperative breast magnetic resonance imaging (MRI) features are associated with the risk of LRR in patients undergoing BCT. Material and Methods A total of 3781 women with primary invasive breast cancer underwent preoperative MRI and BCT between 2003 and 2013. Forty-eight patients who developed LRR comprised the LRR cohort and one-to-one matching (age, tumor stage, grade, and axillary nodal status) of each patient to a control participant was performed in patients who did not develop recurrence. Three readers independently reviewed MR images of the index cancer and the presence of multifocal disease was assessed. Χ² analysis was used to compare imaging and clinical features between LRR and control cohorts, with multivariate logistic regression analysis used to identify independent features. Results Significant differences were found in the proportion of multifocal disease ( P = 0.001), background parenchymal enhancement level ( P = 0.007), and breast cancer molecular subtype ( P = 0.01) between LRR and control cohorts. Multivariate analysis showed that multifocal disease (odds ratio [OR] = 11.9; 95% confidence interval [CI] = 1.4-102.5; P = 0.02) and human epidermal growth factor receptor 2-positive subtype (OR = 12.7; 95% CI = 1.3-127.6; P = 0.03) were both independently associated with LRR. Conclusion Multifocal disease on preoperative breast MRI may indicate an increased risk of LRR in patients treated with BCT.

The crucial role of multiomic approach in cancer research and clinically relevant outcomes

Cancer with heavily economic and social burden is the hot point in the field of medical research. Some remarkable achievements have been made; however, the exact mechanisms of tumor initiation and development remain unclear. Cancer is a complex, whole-body disease that involves multiple abnormalities in the levels of DNA, RNA, protein, metabolite and medical imaging. Biological omics including genomics, transcriptomics, proteomics, metabolomics and radiomics aims to systematically understand carcinogenesis in different biological levels, which is driving the shift of cancer research paradigm from single parameter model to multi-parameter systematical model. The rapid development of various omics technologies is driving one to conveniently get multi-omics data, which accelerates predictive, preventive and personalized medicine (PPPM) practice allowing prediction of response with substantially increased accuracy, stratification of particular patients and eventual personalization of medicine. This review article describes the methodology, advances, and clinically relevant outcomes of different "omics" technologies in cancer research, and especially emphasizes the importance and scientific merit of integrating multi-omics in cancer research and clinically relevant outcomes.

Radiomics Based on Adapted Diffusion Kurtosis Imaging Helps to Clarify Most Mammographic Findings Suspicious for Cancer

Purpose To evaluate a radiomics model of Breast Imaging Reporting and Data System (BI-RADS) 4 and 5 breast lesions extracted from breast-tissue-optimized kurtosis magnetic resonance (MR) imaging for lesion characterization by using a sensitivity threshold similar to that of biopsy. Materials and Methods This institutional study included 222 women at two independent study sites (site 1: training set of 95 patients; mean age ± standard deviation, 58.6 years ± 6.6; 61 malignant and 34 benign lesions; site 2: independent test set of 127 patients; mean age, 58.2 years ± 6.8; 61 malignant and 66 benign lesions). All women presented with a finding suspicious for cancer at x-ray mammography (BI-RADS 4 or 5) and an indication for biopsy. Before biopsy, diffusion-weighted MR imaging (b values, 0-1500 sec/mm²) was performed by using 1.5-T imagers from different MR imaging vendors. Lesions were segmented and voxel-based kurtosis fitting adapted to account for fat signal contamination was performed. A radiomics feature model was developed by using a random forest regressor. The fixed model was tested on an independent test set. Conventional interpretations of MR imaging were also assessed for comparison. Results The radiomics feature model reduced false-positive results from 66 to 20 (specificity 70.0% [46 of 66]) at the predefined sensitivity of greater than 98.0% [60 of 61] in the independent test set, with BI-RADS 4a and 4b lesions benefiting from the analysis (specificity 74.0%, [37 of 50]; 60.0% [nine of 15]) and BI-RADS 5 lesions showing no added benefit. The model significantly improved specificity compared with the median apparent diffusion coefficient (P < .001) and apparent kurtosis coefficient (P = .02) alone. Conventional reading of dynamic contrast material-enhanced MR imaging provided sensitivity of 91.8% (56 of 61) and a specificity of 74.2% (49 of 66). Accounting for fat signal intensity during fitting significantly improved the area under the curve of the model (P = .001). Conclusion A radiomics model based on kurtosis diffusion-weighted imaging performed by using MR imaging machines from different vendors allowed for reliable differentiation between malignant and benign breast lesions in both a training and an independent test data set. ^© RSNA, 2018 Online supplemental material is available for this article.

Dysregulation of long non-coding RNA in breast cancer: an overview of mechanism and clinical implication

Long non-coding RNAs (lncRNAs), which occupy nearly 98% of genome, have crucial roles in cancer development, including breast cancer. Breast cancer is a disease with high incidence. Despite of recent progress in understanding the molecular mechanisms and combined therapy strategies, the functions and mechanisms of lncRNAs in breast cancer remains unclear. This review presents the currently basic knowledge and research approaches of lncRNAs. We also highlight the latest advances of seven classic lncRNAs and three novel lncRNAs in breast cancer, elucidating their mechanisms and possible therapeutic targets. Additionally, association between lncRNA and specific molecular subtype of breast cancer is reported. Lastly, we briefly delineate the potential roles of lncRNAs in clinical applications as biomarkers and treatment targets.

Radiogenomics of hepatocellular carcinoma: multiregion analysis-based identification of prognostic imaging biomarkers by integrating gene data-a preliminary study

Our objective was to identify prognostic imaging biomarkers for hepatocellular carcinoma in contrast-enhanced computed tomography (CECT) with biological interpretations by associating imaging features and gene modules. We retrospectively analyzed 371 patients who had gene expression profiles. For the 38 patients with CECT imaging data, automatic intra-tumor partitioning was performed, resulting in three spatially distinct subregions. We extracted a total of 37 quantitative imaging features describing intensity, geometry, and texture from each subregion. Imaging features were selected after robustness and redundancy analysis. Gene modules acquired from clustering were chosen for their prognostic significance. By constructing an association map between imaging features and gene modules with Spearman rank correlations, the imaging features that significantly correlated with gene modules were obtained. These features were evaluated with Cox's proportional hazard models and Kaplan-Meier estimates to determine their prognostic capabilities for overall survival (OS). Eight imaging features were significantly correlated with prognostic gene modules, and two of them were associated with OS. Among these, the geometry feature volume fraction of the subregion, which was significantly correlated with all prognostic gene modules representing cancer-related interpretation, was predictive of OS (Cox p  =  0.022, hazard ratio  =  0.24). The texture feature cluster prominence in the subregion, which was correlated with the prognostic gene module representing lipid metabolism and complement activation, also had the ability to predict OS (Cox p  =  0.021, hazard ratio  =  0.17). Imaging features depicting the volume fraction and textural heterogeneity in subregions have the potential to be predictors of OS with interpretable biological meaning.

Towards precision medicine: from quantitative imaging to radiomics

Radiology (imaging) and imaging-guided interventions, which provide multi-parametric morphologic and functional information, are playing an increasingly significant role in precision medicine. Radiologists are trained to understand the imaging phenotypes, transcribe those observations (phenotypes) to correlate with underlying diseases and to characterize the images. However, in order to understand and characterize the molecular phenotype (to obtain genomic information) of solid heterogeneous tumours, the advanced sequencing of those tissues using biopsy is required. Thus, radiologists image the tissues from various views and angles in order to have the complete image phenotypes, thereby acquiring a huge amount of data. Deriving meaningful details from all these radiological data becomes challenging and raises the big data issues. Therefore, interest in the application of radiomics has been growing in recent years as it has the potential to provide significant interpretive and predictive information for decision support. Radiomics is a combination of conventional computer-aided diagnosis, deep learning methods, and human skills, and thus can be used for quantitative characterization of tumour phenotypes. This paper discusses the overview of radiomics workflow, the results of various radiomics-based studies conducted using various radiological images such as computed tomography (CT), magnetic resonance imaging (MRI), and positron-emission tomography (PET), the challenges we are facing, and the potential contribution of radiomics towards precision medicine.

High-dimensional regression analysis links magnetic resonance imaging features and protein expression and signaling pathway alterations in breast invasive carcinoma

BACKGROUND

Imaging features derived from MRI scans can be used for not only breast cancer detection and measuring disease extent, but can also determine gene expression and patient outcomes. The relationships between imaging features, gene/protein expression, and response to therapy hold potential to guide personalized medicine. We aim to characterize the relationship between radiologist-annotated tumor phenotypic features (based on MRI) and the underlying biological processes (based on proteomic profiling) in the tumor.

METHODS

Multiple-response regression of the image-derived, radiologist-scored features with reverse-phase protein array expression levels generated association coefficients for each combination of image-feature and protein in the RPPA dataset. Significantly-associated proteins for features were analyzed with Ingenuity Pathway Analysis software. Hierarchical clustering of the results of the pathway analysis determined which features were most strongly correlated with pathway activity and cellular functions.

RESULTS

Each of the twenty-nine imaging features was found to have a set of significantly correlated molecules, associated biological functions, and pathways.

CONCLUSIONS

We interrogated the pathway alterations represented by the protein expression associated with each imaging feature. Our study demonstrates the relationships between biological processes (via proteomic measurements) and MRI features within breast tumors.

A predictive model for distinguishing radiation necrosis from tumour progression after gamma knife radiosurgery based on radiomic features from MR images

OBJECTIVES

To develop a model using radiomic features extracted from MR images to distinguish radiation necrosis from tumour progression in brain metastases after Gamma Knife radiosurgery.

METHODS

We retrospectively identified 87 patients with pathologically confirmed necrosis (24 lesions) or progression (73 lesions) and calculated 285 radiomic features from four MR sequences (T1, T1 post-contrast, T2, and fluid-attenuated inversion recovery) obtained at two follow-up time points per lesion per patient. Reproducibility of each feature between the two time points was calculated within each group to identify a subset of features with distinct reproducible values between two groups. Changes in radiomic features from one time point to the next (delta radiomics) were used to build a model to classify necrosis and progression lesions.

RESULTS

A combination of five radiomic features from both T1 post-contrast and T2 MR images were found to be useful in distinguishing necrosis from progression lesions. Delta radiomic features with a RUSBoost ensemble classifier had an overall predictive accuracy of 73.2% and an area under the curve value of 0.73 in leave-one-out cross-validation.

CONCLUSIONS

Delta radiomic features extracted from MR images have potential for distinguishing radiation necrosis from tumour progression after radiosurgery for brain metastases.

KEY POINTS

• Some radiomic features showed better reproducibility for progressive lesions than necrotic ones • Delta radiomic features can help to distinguish radiation necrosis from tumour progression • Delta radiomic features had better predictive value than did traditional radiomic features.

Biomarkers and Imaging of Breast Cancer

OBJECTIVE

The goals of this review are to provide background information on the definitions and applications of the general term "biomarker" and to highlight the specific roles of breast imaging biomarkers in research and clinical breast cancer care. A search was conducted of the main electronic biomedical databases (PubMed, Cochrane, Embase, MEDLINE [Ovid], Scopus, and Web of Science). The search was focused on review literature in general radiology and biomedical sciences and on reviews and primary research articles on biomarkers in breast imaging over the 15 years ending in June 2017. The keywords included "biomarker," "trial endpoints," "breast imaging," "breast cancer," "radiomics," and "precision medicine" in the titles and abstracts of the papers.

CONCLUSION

Clinical breast care and breast cancer-related research rely on imaging biomarkers for decision support. In the era of precision medicine and big data, the practice of radiology is likely to change. A closer integration of breast imaging with related biomedical fields and the creation of large integrated and shareable databases of clinical, molecular, and imaging biomarkers should allow the field to continue guiding breast cancer care and research.

Background, current role, and potential applications of radiogenomics

With the genomic revolution in the early 1990s, medical research has been driven to study the basis of human disease on a genomic level and to devise precise cancer therapies tailored to the specific genetic makeup of a tumor. To match novel therapeutic concepts conceived in the era of precision medicine, diagnostic tests must be equally sufficient, multilayered, and complex to identify the relevant genetic alterations that render cancers susceptible to treatment. With significant advances in training and medical imaging techniques, image analysis and the development of high-throughput methods to extract and correlate multiple imaging parameters with genomic data, a new direction in medical research has emerged. This novel approach has been termed radiogenomics. Radiogenomics aims to correlate imaging characteristics (ie, the imaging phenotype) with gene expression patterns, gene mutations, and other genome-related characteristics and is designed to facilitate a deeper understanding of tumor biology and capture the intrinsic tumor heterogeneity. Ultimately, the goal of radiogenomics is to develop imaging biomarkers for outcome that incorporate both phenotypic and genotypic metrics. Due to the noninvasive nature of medical imaging and its ubiquitous use in clinical practice, the field of radiogenomics is rapidly evolving and initial results are encouraging. In this article, we briefly discuss the background and then summarize the current role and the potential of radiogenomics in brain, liver, prostate, gynecological, and breast tumors.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2017;47:604-620.

Gene expression profiling of calcifications in breast cancer

We investigated the gene expression profiles of calcifications in breast cancer. Gene expression analysis of surgical specimen was performed using Affymetrix GeneChip® Human Gene 2.0 ST arrays in 168 breast cancer patients. The mammographic calcifications were reviewed by three radiologists and classified into three groups according to malignancy probability: breast cancers without suspicious calcifications; breast cancers with low-to-intermediate suspicious calcifications; and breast cancers with highly suspicious calcifications. To identify differentially expressed genes (DEGs) between these three groups, a one-way analysis of variance was performed with post hoc comparisons with Tukey's honest significant difference test. To explore the biological significance of DEGs, we used DAVID for gene ontology analysis and BioLattice for clustering analysis. A total of 2551 genes showed differential expression among the three groups. ERBB2 genes are up-regulated in breast cancers with highly suspicious calcifications (fold change 2.474, p < 0.001). Gene ontology analysis revealed that the immune, defense and inflammatory responses were decreased in breast cancers with highly suspicious calcifications compared to breast cancers without suspicious calcifications (p from 10^-23 to 10^-8). The clustering analysis also demonstrated that the immune system is associated with mammographic calcifications (p < 0.001). Our study showed calcifications in breast cancers are associated with high levels of mRNA expression of ERBB2 and decreased immune system activity.

Non-Small Cell Lung Cancer Radiogenomics Map Identifies Relationships between Molecular and Imaging Phenotypes with Prognostic Implications

Purpose To create a radiogenomic map linking computed tomographic (CT) image features and gene expression profiles generated by RNA sequencing for patients with non-small cell lung cancer (NSCLC). Materials and Methods A cohort of 113 patients with NSCLC diagnosed between April 2008 and September 2014 who had preoperative CT data and tumor tissue available was studied. For each tumor, a thoracic radiologist recorded 87 semantic image features, selected to reflect radiologic characteristics of nodule shape, margin, texture, tumor environment, and overall lung characteristics. Next, total RNA was extracted from the tissue and analyzed with RNA sequencing technology. Ten highly coexpressed gene clusters, termed metagenes, were identified, validated in publicly available gene-expression cohorts, and correlated with prognosis. Next, a radiogenomics map was built that linked semantic image features to metagenes by using the t statistic and the Spearman correlation metric with multiple testing correction. Results RNA sequencing analysis resulted in 10 metagenes that capture a variety of molecular pathways, including the epidermal growth factor (EGF) pathway. A radiogenomic map was created with 32 statistically significant correlations between semantic image features and metagenes. For example, nodule attenuation and margins are associated with the late cell-cycle genes, and a metagene that represents the EGF pathway was significantly correlated with the presence of ground-glass opacity and irregular nodules or nodules with poorly defined margins. Conclusion Radiogenomic analysis of NSCLC showed multiple associations between semantic image features and metagenes that represented canonical molecular pathways, and it can result in noninvasive identification of molecular properties of NSCLC. Online supplemental material is available for this article.

Heterogeneous Enhancement Patterns of Tumor-adjacent Parenchyma at MR Imaging Are Associated with Dysregulated Signaling Pathways and Poor Survival in Breast Cancer

Purpose To identify the molecular basis of quantitative imaging characteristics of tumor-adjacent parenchyma at dynamic contrast material-enhanced magnetic resonance (MR) imaging and to evaluate their prognostic value in breast cancer. Materials and Methods In this institutional review board-approved, HIPAA-compliant study, 10 quantitative imaging features depicting tumor-adjacent parenchymal enhancement patterns were extracted and screened for prognostic features in a discovery cohort of 60 patients. By using data from The Cancer Genome Atlas (TCGA), a radiogenomic map for the tumor-adjacent parenchymal tissue was created and molecular pathways associated with prognostic parenchymal imaging features were identified. Furthermore, a multigene signature of the parenchymal imaging feature was built in a training cohort (n = 126), and its prognostic relevance was evaluated in two independent cohorts (n = 879 and 159). Results One image feature measuring heterogeneity (ie, information measure of correlation) was significantly associated with prognosis (false-discovery rate < 0.1), and at a cutoff of 0.57 stratified patients into two groups with different recurrence-free survival rates (log-rank P = .024). The tumor necrosis factor signaling pathway was identified as the top enriched pathway (hypergeometric P < .0001) among genes associated with the image feature. A 73-gene signature based on the tumor profiles in TCGA achieved good association with the tumor-adjacent parenchymal image feature (R² = 0.873), which stratified patients into groups regarding recurrence-free survival (log-rank P = .029) and overall survival (log-rank P = .042) in an independent TCGA cohort. The prognostic value was confirmed in another independent cohort (Gene Expression Omnibus GSE 1456), with log-rank P = .00058 for recurrence-free survival and log-rank P = .0026 for overall survival. Conclusion Heterogeneous enhancement patterns of tumor-adjacent parenchyma at MR imaging are associated with the tumor necrosis signaling pathway and poor survival in breast cancer. ^© RSNA, 2017 Online supplemental material is available for this article.

Unsupervised Clustering of Quantitative Image Phenotypes Reveals Breast Cancer Subtypes with Distinct Prognoses and Molecular Pathways

Purpose: To identify novel breast cancer subtypes by extracting quantitative imaging phenotypes of the tumor and surrounding parenchyma and to elucidate the underlying biologic underpinnings and evaluate the prognostic capacity for predicting recurrence-free survival (RFS).Experimental Design: We retrospectively analyzed dynamic contrast-enhanced MRI data of patients from a single-center discovery cohort (n = 60) and an independent multicenter validation cohort (n = 96). Quantitative image features were extracted to characterize tumor morphology, intratumor heterogeneity of contrast agent wash-in/wash-out patterns, and tumor-surrounding parenchyma enhancement. On the basis of these image features, we used unsupervised consensus clustering to identify robust imaging subtypes and evaluated their clinical and biologic relevance. We built a gene expression-based classifier of imaging subtypes and tested their prognostic significance in five additional cohorts with publically available gene expression data but without imaging data (n = 1,160).Results: Three distinct imaging subtypes, that is, homogeneous intratumoral enhancing, minimal parenchymal enhancing, and prominent parenchymal enhancing, were identified and validated. In the discovery cohort, imaging subtypes stratified patients with significantly different 5-year RFS rates of 79.6%, 65.2%, 52.5% (log-rank P = 0.025) and remained as an independent predictor after adjusting for clinicopathologic factors (HR, 2.79; P = 0.016). The prognostic value of imaging subtypes was further validated in five independent gene expression cohorts, with average 5-year RFS rates of 88.1%, 74.0%, 59.5% (log-rank P from <0.0001 to 0.008). Each imaging subtype was associated with specific dysregulated molecular pathways that can be therapeutically targeted.Conclusions: Imaging subtypes provide complimentary value to established histopathologic or molecular subtypes and may help stratify patients with breast cancer. Clin Cancer Res; 23(13); 3334-42. ©2017 AACR.

Associations between gene expression profiles of invasive breast cancer and Breast Imaging Reporting and Data System MRI lexicon

Purpose

To evaluate whether the Breast Imaging Reporting and Data System (BI-RADS) MRI lexicon could reflect the genomic information of breast cancers and to suggest intuitive imaging features as biomarkers.

Methods

Matched breast MRI data from The Cancer Imaging Archive and gene expression profile from The Cancer Genome Atlas of 70 invasive breast cancers were analyzed. Magnetic resonance images were reviewed according to the BI-RADS MRI lexicon of mass morphology. The cancers were divided into 2 groups of gene clustering by gene set enrichment an alysis. Clinicopathologic and imaging characteristics were compared between the 2 groups.

Results

The luminal subtype was predominant in the group 1 gene set and the triple-negative subtype was predominant in the group 2 gene set (55 of 56, 98.2% vs. 9 of 14, 64.3%). Internal enhancement descriptors were different between the 2 groups; heterogeneity was most frequent in group 1 (27 of 56, 48.2%) and rim enhancement was dominant in group 2 (10 of 14, 71.4%). In group 1, the gene sets related to mammary gland development were overexpressed whereas the gene sets related to mitotic cell division were overexpressed in group 2.

Conclusion

We identified intuitive imaging features of breast MRI associated with distinct gene expression profiles using the standard imaging variables of BI-RADS. The internal enhancement pattern on MRI might reflect specific gene expression profiles of breast cancers, which can be recognized by visual distinction.

3D printed pathological sectioning boxes to facilitate radiological–pathological correlation in hepatectomy cases

Radiogenomics promises to identify tumour imaging features indicative of genomic or proteomic aberrations that can be therapeutically targeted allowing precision personalised therapy. An accurate radiological–pathological correlation is critical to the process of radiogenomic characterisation of tumours. An accurate correlation, however, is difficult to achieve with current pathological sectioning techniques which result in sectioning in non-standard planes. The purpose of this work is to present a technique to standardise hepatic sectioning to facilitateradiological–pathological correlation. We describe a process in which three-dimensional (3D)-printed specimen boxes based on preoperative cross-sectional imaging (CT and MRI) can be used to facilitate pathological sectioning in standard planes immediately on hepatic resection enabling improved tumour mapping. We have applied this process in 13 patients undergoing hepatectomy and have observed close correlation between imaging and gross pathology in patients with both unifocal and multifocal tumours.

Computer-aided diagnosis program for classifying the risk of hepatocellular carcinoma on MR images following liver imaging reporting and data system (LI-RADS)

PURPOSE

To develop and evaluate a computer-aided diagnosis (CAD) program for liver lesions on magnetic resonance (MR) images for classification of the risk of hepatocellular carcinoma (HCC) following the liver imaging reporting and data system (LI-RADS).

MATERIALS AND METHODS

Liver MR images from 41 patients with hyperenhancing liver lesions categorized as LR 3, 4, and 5 were evaluated by two radiologists. The major LI-RADS features of each index liver lesion were recorded, including size (maximum transverse diameter), presence of hyperenhancement, washout appearance, and capsule appearance. A CAD program was implemented to register MR images at different contrast-enhancement phases, segment liver lesions, extract lesion features, and classify lesions according to LI-RADS. The LI-RADS features quantified by CAD were compared with those assessed by radiologists using the intraclass correlation coefficient (ICC) and receiver operator curve (ROC) analyses. The LI-RADS categorization between CAD and radiologists was evaluated using the weighted Cohen's kappa coefficient.

RESULTS

The mean and standard deviation of the lesion diameters were 21 ± 11 mm (range, 7-70 mm) by radiologists and 22 ± 11 mm (range, 8-72 mm) by CAD (ICC, 0.96-0.97). The area under the curve (AUC) for the washout assessment by CAD was 0.79-0.93 with sensitivity 0.69-0.82 and specificity 0.79-1. The AUC for the capsule assessment by CAD was 0.79-0.9 with sensitivity 0.75-0.9 and specificity 0.82-0.96. The classifications by the radiologists and CAD coincided in 76-83% lesions (k = 0.57-0.71), while the agreements between radiologists were in 78% lesions (k = 0.59).

CONCLUSION

We developed a CAD program for liver lesions on MR images and showed a substantial agreement in the LI-RADS-based classification of the risk of HCCs between the CAD and radiologists.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:710-722.

Multiregional Radiogenomic Assessment of Prostate Microenvironments with Multiparametric MR Imaging and DNA Whole-Exome Sequencing of Prostate Glands with Adenocarcinoma

Purpose To assess the underlying genomic variation of prostate gland microenvironments of patients with prostate adenocarcinoma in the context of colocalized multiparametric magnetic resonance (MR) imaging and histopathologic assessment of normal and abnormal regions by using whole-exome sequencing. Materials and Methods Six patients with prostate adenocarcinoma who underwent robotic prostatectomy with whole-mount preservation of the prostate were identified, which enabled spatial mapping between preoperative multiparametric MR imaging and the gland. Four regions of interest were identified within each gland, including regions found to be normal and abnormal via histopathologic analysis. Whole-exome DNA sequencing (>50 times coverage) was performed on each of these spatially targeted regions. Radiogenomic analysis of imaging and mutation data were performed with hierarchical clustering, phylogenetic analysis, and principal component analysis. Results Radiogenomic multiparametric MR imaging and whole-exome spatial characterization in six patients with prostate adenocarcinoma (three patients, Gleason score of 3 + 4; and three patients, Gleason score of 4 + 5) was performed across 23 spatially distinct regions. Hierarchical clustering separated histopathologic analysis-proven high-grade lesions from the normal regions, and this reflected concordance between multiparametric MR imaging and resultant histopathologic analysis in all patients. Seventy-seven mutations involving 29 cancer-associated genes across the 23 spatially distinct prostate samples were identified. There was no significant difference in mutation load in cancer-associated genes between regions that were proven to be normal via histopathologic analysis (34 mutations per sample ± 19), mildly suspicious via multiparametric MR imaging (37 mutations per sample ± 21), intermediately suspicious via multiparametric MR imaging (31 mutations per sample ± 15), and high-grade cancer (33 mutations per sample ± 18) (P = .30). Principal component analysis resolved samples from different patients and further classified samples (regardless of histopathologic status) from prostate glands with Gleason score 3 + 4 versus 4 + 5 samples. Conclusion Multiregion spatial multiparametric MR imaging and whole-exome radiogenomic analysis of prostate glands with adenocarcinoma shows a continuum of mutations across regions that were found via histologic analysis to be high grade and normal. ^© RSNA, 2017 Online supplemental material is available for this article.

Radiogenomic Analysis of Oncological Data: A Technical Survey

In the last few years, biomedical research has been boosted by the technological development of analytical instrumentation generating a large volume of data. Such information has increased in complexity from basic (i.e., blood samples) to extensive sets encompassing many aspects of a subject phenotype, and now rapidly extending into genetic and, more recently, radiomic information. Radiogenomics integrates both aspects, investigating the relationship between imaging features and gene expression. From a methodological point of view, radiogenomics takes advantage of non-conventional data analysis techniques that reveal meaningful information for decision-support in cancer diagnosis and treatment. This survey is aimed to review the state-of-the-art techniques employed in radiomics and genomics with special focus on analysis methods based on molecular and multimodal probes. The impact of single and combined techniques will be discussed in light of their suitability in correlation and predictive studies of specific oncologic diseases.

lncRNA MEG3 had anti-cancer effects to suppress pancreatic cancer activity

AIM

The aim of this study was to explain the mechanism of lncRNA MEG 3 in pancreatic cancer.

METHODS

We were collecting 30 pancreatic cancer patients, taking the sample from these patients. We measured the PI3K protein expressions from 30 patients by IHC and WB methods and MEG 3 expression by RT-PCR, and analyzed the relationship between PI3K protein expression and pancreatic cancer patients' clinical pathology and the correlation between lncRNA MEG 3 and PI3K. In the cell experiment, PANC-1 cells were divided into three groups: NC, BL and lncRNA groups, after treatment,we measured cell proliferation rate of 3 groups by MTT methods, evaluated cell apoptosis and cell cycle using flow cytometry, tested the invasion cells and migrate rate of 3 groups by transwell and wound healing assays.

RESULTS

Compared with carcinoma adjacent tissue, The PI3K protein expression of pancreatic cancer tissue were significantly up-regulation (P>0.05). MEG 3 gene expression was negatively correlated with PI3K expression. The MEG 3 was negatively correlated with tumor size, Metastasis and Vascular invasion in pancreatic cancer (P<0.05, respectively). In the cell experiment, The cell proliferation and apoptosis rates of lncRNA group were significantly difference compared with NC group (P<0.05, respectively), and the G1 phase rate of lncRNA group was higher than NC group (P<0.05). The invasion cells and wound healing rate were significantly reduced in lncRNA group than those in NC group (P<0.05, respectively).

CONCLUSION

MEG 3 over-expressing had anti-cancer effects to suppress pancreatic cancer activity by regulation PI3K/AKT/Bcl-2/Bax/Cyclin D1/P53 and PI3K/AKT/MMP-2/MMP-9 signaling pathways.

Identifying relations between imaging phenotypes and molecular subtypes of breast cancer: Model discovery and external validation

PURPOSE

To determine whether dynamic contrast enhancement magnetic resonance imaging (DCE-MRI) characteristics of the breast tumor and background parenchyma can distinguish molecular subtypes (ie, luminal A/B or basal) of breast cancer.

MATERIALS AND METHODS

In all, 84 patients from one institution and 126 patients from The Cancer Genome Atlas (TCGA) were used for discovery and external validation, respectively. Thirty-five quantitative image features were extracted from DCE-MRI (1.5 or 3T) including morphology, texture, and volumetric features, which capture both tumor and background parenchymal enhancement (BPE) characteristics. Multiple testing was corrected using the Benjamini-Hochberg method to control the false-discovery rate (FDR). Sparse logistic regression models were built using the discovery cohort to distinguish each of the three studied molecular subtypes versus the rest, and the models were evaluated in the validation cohort.

RESULTS

On univariate analysis in discovery and validation cohorts, two features characterizing tumor and two characterizing BPE were statistically significant in separating luminal A versus nonluminal A cancers; two features characterizing tumor were statistically significant for separating luminal B; one feature characterizing tumor and one characterizing BPE reached statistical significance for distinguishing basal (Wilcoxon P < 0.05, FDR < 0.25). In discovery and validation cohorts, multivariate logistic regression models achieved an area under the receiver operator characteristic curve (AUC) of 0.71 and 0.73 for luminal A cancer, 0.67 and 0.69 for luminal B cancer, and 0.66 and 0.79 for basal cancer, respectively.

CONCLUSION

DCE-MRI characteristics of breast cancer and BPE may potentially be used to distinguish among molecular subtypes of breast cancer.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2017;46:1017-1027.

Prediction of malignancy by a radiomic signature from contrast agent-free diffusion MRI in suspicious breast lesions found on screening mammography

PURPOSE

To assess radiomics as a tool to determine how well lesions found suspicious on breast cancer screening X-ray mammography can be categorized into malignant and benign with unenhanced magnetic resonance (MR) mammography with diffusion-weighted imaging and T₂ -weighted sequences.

MATERIALS AND METHODS

From an asymptomatic screening cohort, 50 women with mammographically suspicious findings were examined with contrast-enhanced breast MRI (ceMRI) at 1.5T. Out of this protocol an unenhanced, abbreviated diffusion-weighted imaging protocol (ueMRI) including T₂ -weighted, (T₂ w), diffusion-weighted imaging (DWI), and DWI with background suppression (DWIBS) sequences and corresponding apparent diffusion coefficient (ADC) maps were extracted. From ueMRI-derived radiomic features, three Lasso-supervised machine-learning classifiers were constructed and compared with the clinical performance of a highly experienced radiologist: 1) univariate mean ADC model, 2) unconstrained radiomic model, 3) constrained radiomic model with mandatory inclusion of mean ADC.

RESULTS

The unconstrained and constrained radiomic classifiers consisted of 11 parameters each and achieved differentiation of malignant from benign lesions with a .632 + bootstrap receiver operating characteristics (ROC) area under the curve (AUC) of 84.2%/85.1%, compared to 77.4% for mean ADC and 95.9%/95.9% for the experienced radiologist using ceMRI/ueMRI.

CONCLUSION

In this pilot study we identified two ueMRI radiomics classifiers that performed well in the differentiation of malignant from benign lesions and achieved higher performance than the mean ADC parameter alone. Classification was lower than the almost perfect performance of a highly experienced breast radiologist. The potential of radiomics to provide a training-independent diagnostic decision tool is indicated. A performance reaching the human expert would be highly desirable and based on our results is considered possible when the concept is extended in larger cohorts with further development and validation of the technique.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;46:604-616.

Predictive radiogenomics modeling of EGFR mutation status in lung cancer

Molecular analysis of the mutation status for EGFR and KRAS are now routine in the management of non-small cell lung cancer. Radiogenomics, the linking of medical images with the genomic properties of human tumors, provides exciting opportunities for non-invasive diagnostics and prognostics. We investigated whether EGFR and KRAS mutation status can be predicted using imaging data. To accomplish this, we studied 186 cases of NSCLC with preoperative thin-slice CT scans. A thoracic radiologist annotated 89 semantic image features of each patient’s tumor. Next, we built a decision tree to predict the presence of EGFR and KRAS mutations. We found a statistically significant model for predicting EGFR but not for KRAS mutations. The test set area under the ROC curve for predicting EGFR mutation status was 0.89. The final decision tree used four variables: emphysema, airway abnormality, the percentage of ground glass component and the type of tumor margin. The presence of either of the first two features predicts a wild type status for EGFR while the presence of any ground glass component indicates EGFR mutations. These results show the potential of quantitative imaging to predict molecular properties in a non-invasive manner, as CT imaging is more readily available than biopsies.

Radiomic Profiling of Glioblastoma: Identifying an Imaging Predictor of Patient Survival with Improved Performance over Established Clinical and Radiologic Risk Models

Purpose To evaluate whether radiomic feature-based magnetic resonance (MR) imaging signatures allow prediction of survival and stratification of patients with newly diagnosed glioblastoma with improved accuracy compared with that of established clinical and radiologic risk models. Materials and Methods Retrospective evaluation of data was approved by the local ethics committee and informed consent was waived. A total of 119 patients (allocated in a 2:1 ratio to a discovery [n = 79] or validation [n = 40] set) with newly diagnosed glioblastoma were subjected to radiomic feature extraction (12 190 features extracted, including first-order, volume, shape, and texture features) from the multiparametric (contrast material-enhanced T1-weighted and fluid-attenuated inversion-recovery imaging sequences) and multiregional (contrast-enhanced and unenhanced) tumor volumes. Radiomic features of patients in the discovery set were subjected to a supervised principal component (SPC) analysis to predict progression-free survival (PFS) and overall survival (OS) and were validated in the validation set. The performance of a Cox proportional hazards model with the SPC analysis predictor was assessed with C index and integrated Brier scores (IBS, lower scores indicating higher accuracy) and compared with Cox models based on clinical (age and Karnofsky performance score) and radiologic (Gaussian normalized relative cerebral blood volume and apparent diffusion coefficient) parameters. Results SPC analysis allowed stratification based on 11 features of patients in the discovery set into a low- or high-risk group for PFS (hazard ratio [HR], 2.43; P = .002) and OS (HR, 4.33; P < .001), and the results were validated successfully in the validation set for PFS (HR, 2.28; P = .032) and OS (HR, 3.45; P = .004). The performance of the SPC analysis (OS: IBS, 0.149; C index, 0.654; PFS: IBS, 0.138; C index, 0.611) was higher compared with that of the radiologic (OS: IBS, 0.175; C index, 0.603; PFS: IBS, 0.149; C index, 0.554) and clinical risk models (OS: IBS, 0.161, C index, 0.640; PFS: IBS, 0.139; C index, 0.599). The performance of the SPC analysis model was further improved when combined with clinical data (OS: IBS, 0.142; C index, 0.696; PFS: IBS, 0.132; C index, 0.637). Conclusion An 11-feature radiomic signature that allows prediction of survival and stratification of patients with newly diagnosed glioblastoma was identified, and improved performance compared with that of established clinical and radiologic risk models was demonstrated. (©) RSNA, 2016 Online supplemental material is available for this article.

Noninvasive strategies for breast cancer early detection

Breast cancer screening and presurgical diagnosis are currently based on mammography, ultrasound and more sensitive imaging technologies; however, noninvasive biomarkers represent both a challenge and an opportunity for early detection of cancer. An extensive number of potential breast cancer biomarkers have been discovered by microarray hybridization or sequencing of circulating DNA, noncoding RNA and blood cell RNA; multiplex analysis of immune-related molecules and mass spectrometry-based approaches for high-throughput detection of protein, endogenous peptides, circulating and volatile metabolites. However, their medical relevance and their translation to clinics remain to be exploited. Once they will be fully validated, cancer biomarkers, used in combination with the current and emerging imaging technologies, represent an avenue to a personalized breast cancer diagnosis.

Quantitative MRI radiomics in the prediction of molecular classifications of breast cancer subtypes in the TCGA/TCIA data set

Using quantitative radiomics, we demonstrate that computer-extracted magnetic resonance (MR) image-based tumor phenotypes can be predictive of the molecular classification of invasive breast cancers. Radiomics analysis was performed on 91 MRIs of biopsy-proven invasive breast cancers from National Cancer Institute’s multi-institutional TCGA/TCIA. Immunohistochemistry molecular classification was performed including estrogen receptor, progesterone receptor, human epidermal growth factor receptor 2, and for 84 cases, the molecular subtype (normal-like, luminal A, luminal B, HER2-enriched, and basal-like). Computerized quantitative image analysis included: three-dimensional lesion segmentation, phenotype extraction, and leave-one-case-out cross validation involving stepwise feature selection and linear discriminant analysis. The performance of the classifier model for molecular subtyping was evaluated using receiver operating characteristic analysis. The computer-extracted tumor phenotypes were able to distinguish between molecular prognostic indicators; area under the ROC curve values of 0.89, 0.69, 0.65, and 0.67 in the tasks of distinguishing between ER+ versus ER−, PR+ versus PR−, HER2+ versus HER2−, and triple-negative versus others, respectively. Statistically significant associations between tumor phenotypes and receptor status were observed. More aggressive cancers are likely to be larger in size with more heterogeneity in their contrast enhancement. Even after controlling for tumor size, a statistically significant trend was observed within each size group (P=0.04 for lesions ⩽2 cm; P=0.02 for lesions >2 to ⩽5 cm) as with the entire data set (P-value=0.006) for the relationship between enhancement texture (entropy) and molecular subtypes (normal-like, luminal A, luminal B, HER2-enriched, basal-like). In conclusion, computer-extracted image phenotypes show promise for high-throughput discrimination of breast cancer subtypes and may yield a quantitative predictive signature for advancing precision medicine.

Radiogenomic Analysis Demonstrates Associations between (18)F-Fluoro-2-Deoxyglucose PET, Prognosis, and Epithelial-Mesenchymal Transition in Non-Small Cell Lung Cancer

Purpose To investigate whether non-small cell lung cancer (NSCLC) tumors that express high normalized maximum standardized uptake value (SUVmax) are associated with a more epithelial-mesenchymal transition (EMT)-like phenotype. Materials and Methods In this institutional review board-approved study, a public NSCLC data set that contained fluorine 18 ((18)F) fluoro-2-deoxyglucose positron emission tomography (PET) and messenger RNA expression profile data (n = 26) was obtained, and patients were categorized on the basis of measured normalized SUVmax values. Significance analysis of microarrays was then used to create a radiogenomic signature. The prognostic ability of this signature was assessed in a second independent data set that consisted of clinical and messenger RNA expression data (n = 166). Signature concordance with EMT was evaluated by means of validation in a publicly available cell line data set. Finally, by establishing an in vitro EMT lung cancer cell line model, an attempt was made to substantiate the radiogenomic signature with quantitative polymerase chain reaction, and functional assays were performed, including Western blot, cell migration, glucose transporter, and hexokinase assays (paired t test), as well as pharmacologic assays against chemotherapeutic agents (half-maximal effective concentration). Results Differential expression analysis yielded a 14-gene radiogenomic signature (P < .05, false discovery rate [FDR] < 0.20), which was confirmed to have differences in disease-specific survival (log-rank test, P = .01). This signature also significantly overlapped with published EMT cell line gene expression data (P < .05, FDR < 0.20). Finally, an EMT cell line model was established, and cells that had undergone EMT differentially expressed this signature and had significantly different EMT protein expression (P < .05, FDR < 0.20), cell migration, glucose uptake, and hexokinase activity (paired t test, P < .05). Cells that had undergone EMT also had enhanced chemotherapeutic resistance, with a higher half-maximal effective concentration than that of cells that had not undergone EMT (P < .05). Conclusion Integrative radiogenomic analysis demonstrates an association between increased normalized (18)F fluoro-2-deoxyglucose PET SUVmax, outcome, and EMT in NSCLC. (©) RSNA, 2016 Online supplemental material is available for this article.

Imaging genomics in cancer research: limitations and promises

Recently, radiogenomics or imaging genomics has emerged as a novel high-throughput method of associating imaging features with genomic data. Radiogenomics has the potential to provide comprehensive intratumour, intertumour and peritumour information non-invasively. This review article summarizes the current state of radiogenomic research in tumour characterization, discusses some of its limitations and promises and projects its future directions. Semi-radiogenomic studies that relate specific gene expressions to imaging features will also be briefly reviewed.

Breast MRI radiogenomics: Current status and research implications

Breast magnetic resonance imaging (MRI) radiogenomics is an emerging area of research that has the potential to directly influence clinical practice. Clinical MRI scanners today are capable of providing excellent temporal and spatial resolution, which allows extraction of numerous imaging features via human extraction approaches or complex computer vision algorithms. Meanwhile, advances in breast cancer genetics research has resulted in the identification of promising genes associated with cancer outcomes. In addition, validated genomic signatures have been developed that allow categorization of breast cancers into distinct molecular subtypes as well as predict the risk of cancer recurrence and response to therapy. Current radiogenomics research has been directed towards exploratory analysis of individual genes, understanding tumor biology, and developing imaging surrogates to genetic analysis with the long-term goal of developing a meaningful tool for clinical care. The background of breast MRI radiogenomics research, image feature extraction techniques, approaches to radiogenomics research, and promising areas of investigation are reviewed. J. Magn. Reson. Imaging 2016;43:1269-1278.

Using computer-extracted image phenotypes from tumors on breast magnetic resonance imaging to predict breast cancer pathologic stage

BACKGROUND

The objective of this study was to demonstrate that computer-extracted image phenotypes (CEIPs) of biopsy-proven breast cancer on magnetic resonance imaging (MRI) can accurately predict pathologic stage.

METHODS

The authors used a data set of deidentified breast MRIs organized by the National Cancer Institute in The Cancer Imaging Archive. In total, 91 biopsy-proven breast cancers were analyzed from patients who had information available on pathologic stage (stage I, n = 22; stage II, n = 58; stage III, n = 11) and surgically verified lymph node status (negative lymph nodes, n = 46; ≥ 1 positive lymph node, n = 44; no lymph nodes examined, n = 1). Tumors were characterized according to 1) radiologist-measured size and 2) CEIP. Then, models were built that combined 2 CEIPs to predict tumor pathologic stage and lymph node involvement, and the models were evaluated in a leave-1-out, cross-validation analysis with the area under the receiver operating characteristic curve (AUC) as the value of interest.

RESULTS

Tumor size was the most powerful predictor of pathologic stage, but CEIPs that captured biologic behavior also emerged as predictive (eg, stage I and II vs stage III demonstrated an AUC of 0.83). No size measure was successful in the prediction of positive lymph nodes, but adding a CEIP that described tumor "homogeneity" significantly improved discrimination (AUC = 0.62; P = .003) compared with chance.

CONCLUSIONS

The current results indicate that MRI phenotypes have promise for predicting breast cancer pathologic stage and lymph node status. Cancer 2016;122:748-757. © 2015 American Cancer Society.

The radiogenomic risk score stratifies outcomes in a renal cell cancer phase 2 clinical trial

OBJECTIVES

To characterize a radiogenomic risk score (RRS), a previously defined biomarker, and to evaluate its potential for stratifying radiological progression-free survival (rPFS) in patients with metastatic renal cell carcinoma (mRCC) undergoing pre-surgical treatment with bevacizumab.

METHODOLOGY

In this IRB-approved study, prospective imaging analysis of the RRS was performed on phase II clinical trial data of mRCC patients (n = 41) evaluating whether patient stratification according to the RRS resulted in groups more or less likely to have a rPFS to pre-surgical bevacizumab prior to cytoreductive nephrectomy. Survival times of RRS subgroups were analyzed using Kaplan-Meier survival analysis.

RESULTS

The RRS is enriched in diverse molecular processes including drug response, stress response, protein kinase regulation, and signal transduction pathways (P < 0.05). The RRS successfully stratified rPFS to bevacizumab based on pre-treatment computed tomography imaging with a median progression-free survival of 6 versus >25 months (P = 0.005) and overall survival of 25 versus >37 months in the high and low RRS groups (P = 0.03), respectively. Conventional prognostic predictors including the Motzer and Heng criteria were not predictive in this cohort (P > 0.05).

CONCLUSIONS

The RRS stratifies rPFS to bevacizumab in patients from a phase II clinical trial with mRCC undergoing cytoreductive nephrectomy and pre-surgical bevacizumab.

KEY POINTS

• The RRS SOMA stratifies patient outcomes in a phase II clinical trial. • RRS stratifies subjects into prognostic groups in a discrete or continuous fashion. • RRS is biologically enriched in diverse processes including drug response programs.

The Radiogenomic Risk Score: Construction of a Prognostic Quantitative, Noninvasive Image-based Molecular Assay for Renal Cell Carcinoma

PURPOSE

To evaluate the feasibility of constructing radiogenomic-based surrogates of molecular assays (SOMAs) in patients with clear-cell renal cell carcinoma (CCRCC) by using data extracted from a single computed tomographic (CT) image.

MATERIALS AND METHODS

In this institutional review board approved study, gene expression profile data and contrast material-enhanced CT images from 70 patients with CCRCC in a training set were independently assessed by two radiologists for a set of predefined imaging features. A SOMA for a previously validated CCRCC-specific supervised principal component (SPC) risk score prognostic gene signature was constructed and termed the radiogenomic risk score (RRS). It uses the microarray data and a 28-trait image array to evaluate each CT image with multiple regression of gene expression analysis. The predictive power of the RRS SOMA was then prospectively validated in an independent dataset to confirm its relationship to the SPC gene signature (n = 70) and determination of patient outcome (n = 77). Data were analyzed by using multivariate linear regression-based methods and Cox regression modeling, and significance was assessed with receiver operator characteristic curves and Kaplan-Meier survival analysis.

RESULTS

Our SOMA faithfully represents the tissue-based molecular assay it models. The RRS scaled with the SPC gene signature (R = 0.57, P < .001, classification accuracy 70.1%, P < .001) and predicted disease-specific survival (log rank P < .001). Independent validation confirmed the relationship between the RRS and the SPC gene signature (R = 0.45, P < .001, classification accuracy 68.6%, P < .001) and disease-specific survival (log-rank P < .001) and that it was independent of stage, grade, and performance status (multivariate Cox model P < .05, log-rank P < .001).

CONCLUSION

A SOMA for the CCRCC-specific SPC prognostic gene signature that is predictive of disease-specific survival and independent of stage was constructed and validated, confirming that SOMA construction is feasible.