Quantitative imaging biomarkers in the clinical development of targeted therapeutics: current and future perspectives

Improving performance in colorectal cancer histology decomposition using deep and ensemble machine learning

In routine colorectal cancer management, histologic samples stained with hematoxylin and eosin are commonly used. Nonetheless, their potential for defining objective biomarkers for patient stratification and treatment selection is still being explored. The current gold standard relies on expensive and time-consuming genetic tests. However, recent research highlights the potential of convolutional neural networks (CNNs) to facilitate the extraction of clinically relevant biomarkers from these readily available images. These CNN-based biomarkers can predict patient outcomes comparably to golden standards, with the added advantages of speed, automation, and minimal cost. The predictive potential of CNN-based biomarkers fundamentally relies on the ability of CNNs to accurately classify diverse tissue types from whole slide microscope images. Consequently, enhancing the accuracy of tissue class decomposition is critical to amplifying the prognostic potential of imaging-based biomarkers. This study introduces a hybrid deep transfer learning and ensemble machine learning model that improves upon previous approaches, including a transformer and neural architecture search baseline for this task. We employed a pairing of the EfficientNetV2 architecture with a random forest classification head. Our model achieved 96.74% accuracy (95% CI: 96.3%-97.1%) on the external test set and 99.89% on the internal test set. Recognizing the potential of these models in the task, we have made them publicly available.

Monitoring imatinib decreasing pericyte coverage and HIF-1α level in a colorectal cancer model by an ultrahigh-field multiparametric MRI approach

BACKGROUND

Excessive pericyte coverage promotes tumor growth, and a downregulation may solve this dilemma. Due to the double-edged sword role of vascular pericytes in tumor microenvironment (TME), indiscriminately decreasing pericyte coverage by imatinib causes poor treatment outcomes. Here, we optimized the use of imatinib in a colorectal cancer (CRC) model in high pericyte-coverage status, and revealed the value of multiparametric magnetic resonance imaging (mpMRI) at 9.4T in monitoring treatment-related changes in pericyte coverage and the TME.

METHODS

CRC xenograft models were evaluated by histological vascular characterizations and mpMRI. Mice with the highest pericyte coverage were treated with imatinib or saline; then, vascular characterizations, tumor apoptosis and HIF-1α level were analyzed histologically, and alterations in the expression of Bcl-2/bax pathway were assessed through qPCR. The effects of imatinib were monitored by dynamic contrast-enhanced (DCE)-, diffusion-weighted imaging (DWI)- and amide proton transfer chemical exchange saturation transfer (APT CEST)-MRI at 9.4T.

RESULTS

The DCE- parameters provided a good histologic match the tumor vascular characterizations. In the high pericyte coverage status, imatinib exhibited significant tumor growth inhibition, necrosis increase and pericyte coverage downregulation, and these changes were accompanied by increased vessel permeability, decreased microvessel density (MVD), increased tumor apoptosis and altered gene expression of apoptosis-related Bcl-2/bax pathway. Strategically, a 4-day imatinib effectively decreased pericyte coverage and HIF-1α level, and continuous treatment led to a less marked decrease in pericyte coverage and re-elevated HIF-1α level. Correlation analysis confirmed the feasibility of using mpMRI parameters to monitor imatinib treatment, with DCE-derived Ve and Ktrans being most correlated with pericyte coverage, Ve with vessel permeability, AUC with microvessel density (MVD), DWI-derived ADC with tumor apoptosis, and APT CEST-derived MTRasym at 1 µT with HIF-1α.

CONCLUSIONS

These results provided an optimized imatinib regimen to achieve decreasing pericyte coverage and HIF-1α level in the high pericyte-coverage CRC model, and offered an ultrahigh-field multiparametric MRI approach for monitoring pericyte coverage and dynamics response of the TME to treatment.

Fractals in the Neurosciences: A Translational Geographical Approach

The chapter presents three new fractal indices (fractal fragmentation index, fractal tentacularity index, and fractal anisotropy index) and normalized Kolmogorov complexity with proven applicability in geographic research, developed by the authors, and the possibility of their future use in neuroscience. The research demonstrates the relevance of fractal analysis in different fields and the basic concepts and principles of fractal geometry being sufficient for the development of models relevant to the studied reality. Also, the research highlighted the need to continue interdisciplinary research based on known fractal indicators, as well as the development of new analysis methods with the translational potential between fields.

Improved accuracy in colorectal cancer tissue decomposition through refinement of established deep learning solutions

Hematoxylin and eosin-stained biopsy slides are regularly available for colorectal cancer patients. These slides are often not used to define objective biomarkers for patient stratification and treatment selection. Standard biomarkers often pertain to costly and slow genetic tests. However, recent work has shown that relevant biomarkers can be extracted from these images using convolutional neural networks (CNNs). The CNN-based biomarkers predicted colorectal cancer patient outcomes comparably to gold standards. Extracting CNN-biomarkers is fast, automatic, and of minimal cost. CNN-based biomarkers rely on the ability of CNNs to recognize distinct tissue types from microscope whole slide images. The quality of these biomarkers (coined 'Deep Stroma') depends on the accuracy of CNNs in decomposing all relevant tissue classes. Improving tissue decomposition accuracy is essential for improving the prognostic potential of CNN-biomarkers. In this study, we implemented a novel training strategy to refine an established CNN model, which then surpassed all previous solutions . We obtained a 95.6% average accuracy in the external test set and 99.5% in the internal test set. Our approach reduced errors in biomarker-relevant classes, such as Lymphocytes, and was the first to include interpretability methods. These methods were used to better apprehend our model's limitations and capabilities.

Quantifying the spatial clustering characteristics of radiographic emphysema explains variability in pulmonary function

Quantitative assessment of emphysema in CT scans has mostly focused on calculating the percentage of lung tissue that is deemed abnormal based on a density thresholding strategy. However, this overall measure of disease burden discards virtually all the spatial information encoded in the scan that is implicitly utilized in a visual assessment. This simplification is likely grouping heterogenous disease patterns and is potentially obscuring clinical phenotypes and variable disease outcomes. To overcome this, several methods that attempt to quantify heterogeneity in emphysema distribution have been proposed. Here, we compare three of those: one based on estimating a power law for the size distribution of contiguous emphysema clusters, a second that looks at the number of emphysema-to-emphysema voxel adjacencies, and a third that applies a parametric spatial point process model to the emphysema voxel locations. This was done using data from 587 individuals from Phase 1 of COPDGene that had an inspiratory CT scan and plasma protein abundance measurements. The associations between these imaging metrics and visual assessment with clinical measures (FEV[Formula: see text], FEV[Formula: see text]-FVC ratio, etc.) and plasma protein biomarker levels were evaluated using a variety of regression models. Our results showed that a selection of spatial measures had the ability to discern heterogeneous patterns among CTs that had similar emphysema burdens. The most informative quantitative measure, average cluster size from the point process model, showed much stronger associations with nearly every clinical outcome examined than existing CT-derived emphysema metrics and visual assessment. Moreover, approximately 75% more plasma biomarkers were found to be associated with an emphysema heterogeneity phenotype when accounting for spatial clustering measures than when they were excluded.

Prognostic prediction based on histopathologic features of tumor microenvironment in colorectal cancer

Purpose

To automatically quantify colorectal tumor microenvironment (TME) in hematoxylin and eosin stained whole slide images (WSIs), and to develop a TME signature for prognostic prediction in colorectal cancer (CRC).

Methods

A deep learning model based on VGG19 architecture and transfer learning strategy was trained to recognize nine different tissue types in whole slide images of patients with CRC. Seven of the nine tissue types were defined as TME components besides background and debris. Then 13 TME features were calculated based on the areas of TME components. A total of 562 patients with gene expression data, survival information and WSIs were collected from The Cancer Genome Atlas project for further analysis. A TME signature for prognostic prediction was developed and validated using Cox regression method. A prognostic prediction model combined the TME signature and clinical variables was also established. At last, gene-set enrichment analysis was performed to identify the significant TME signature associated pathways by querying Gene Ontology database and Kyoto Encyclopedia of Genes and Genomes database.

Results

The deep learning model achieved an accuracy of 94.2% for tissue type recognition. The developed TME signature was found significantly associated to progression-free survival. The clinical combined model achieved a concordance index of 0.714. Gene-set enrichment analysis revealed the TME signature associated genes were enriched in neuroactive ligand-receptor interaction pathway.

Conclusion

The TME signature was proved to be a prognostic factor and the associated biologic pathways would be beneficial to a better understanding of TME in CRC patients.

Directionally Sensitive Fractal Radiomics Compatible With Irregularly Shaped Magnetic Resonance Tumor Regions of Interest: Association With Osteosarcoma Chemoresistance

BACKGROUND

Computational analysis of routinely acquired MRI has potential to improve the tumor chemoresistance prediction and to provide decision support in precision medicine, which may extend patient survival. Most radiomic analytical methods are compatible only with rectangular regions of interest (ROIs) and irregular tumor shape is therefore an important limitation. Furthermore, the currently used analytical methods are not directionally sensitive.

PURPOSE

To implement a tumor analysis that is directionally sensitive and compatible with irregularly shaped ROIs.

STUDY TYPE

Retrospective.

SUBJECTS

A total of 54 patients with histopathologic diagnosis of primary osteosarcoma on tubular long bones and with prechemotherapy MRI.

FIELD STRENGTH/SEQUENCE

A 1.5 T, T2-weighted-short-tau-inversion-recovery-fast-spin-echo.

ASSESSMENT

A model to explore associations with osteosarcoma chemo-responsiveness included MRI data obtained before OsteoSa MAP neoadjuvant cytotoxic chemotherapy. Osteosarcoma morphology was analyzed in the MRI data by calculation of the nondirectional two-dimensional (2D) and directional and nondirectional one-dimensional (1D) Higuchi dimensions (Dh). MAP chemotherapy response was assessed by histopathological necrosis.

STATISTICAL TESTS

The area under the receiver operating characteristic (ROC) curve (AUC) evaluated the association of the calculated features with the actual chemoresponsiveness, using tumor histopathological necrosis (95%) as the endpoint. Least absolute shrinkage and selection operator (LASSO) machine learning and multivariable regression were used for feature selection. Significance was set at <0.05.

RESULTS

The nondirectional 1D Dh reached an AUC of 0.88 in association with the 95% tumor necrosis, while the directional 1D analysis along 180 radial lines significantly improved this association according to the Hanley/McNeil test, reaching an AUC of 0.95. The model defined by variable selection using LASSO reached an AUC of 0.98. The directional analysis showed an optimal predictive range between 90° and 97° and revealed structural osteosarcoma anisotropy manifested by its directionally dependent textural properties.

DATA CONCLUSION

Directionally sensitive radiomics had superior predictive performance in comparison to the standard nondirectional image analysis algorithms with AUCs reaching 0.95 and full compatibility with irregularly shaped ROIs.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 1.

High frequency ultrasound vibrational shear wave elastography for preclinical research

Preclinical evaluation of novel therapies using models of cancer is an important tool in cancer research, where imaging can provide non-invasive tools to characterise the internal structure and function of tumours. The short propagation paths when imaging tumours and organs in small animals allow the use of high frequencies for both ultrasound and shear waves, providing the opportunity for high-resolution shear wave elastography and hence its use for studying the heterogeneity of tissue elasticity, where heterogeneity may be a predictor of tissue response. Here we demonstrate vibrational shear wave elastography (VSWE) using a mechanical actuator to produce high frequency (up to 1000 Hz) shear waves in preclinical tumours, an alternative to the majority of preclinical ultrasound SWE studies where an acoustic radiation force impulse is required to create a relatively low-frequency broad-band shear-wave pulse. We implement VSWE with a high frequency (17.8 MHz) probe running a focused line-by-line ultrasound imaging sequence which as expected was found to offer improved detection of 1000 Hz shear waves over an ultrafast planar wave imaging sequence in a homogenous tissue-mimicking phantom. We test the VSWE in anex vivotumour xenograft, demonstrating the ability to detect shear waves up to 10 mm from the contactor position at 1000 Hz. By reducing the kernel size used for shear wave speed estimation to 1 mm we are able to produce shear wave speed images with spatial resolution of this order. Finally, we present VSWE data from xenograft tumoursin vivo, demonstrating the feasibility of the technique in mice under isoflurane sedation. Mean shear wave speeds in the tumours are in good agreements with those reported by previous authors. Characterising the frequency dependence of shear wave speed demonstrates the potential to quantify the viscoelastic properties of tumoursin vivo.

Predictive biomarkers for personalized medicine in breast cancer

Breast cancer is one of the most frequent malignancies among women worldwide. Based on clinical and molecular features of breast tumors, patients are treated with chemotherapy, hormonal therapy and/or radiotherapy and more recently with immunotherapy or targeted therapy. These different therapeutic options have markedly improved patient outcomes. However, further improvement is needed to fight against resistance to treatment. In the rapidly growing area of research for personalized medicine, predictive biomarkers - which predict patient response to therapy - are essential tools to select the patients who are most likely to benefit from the treatment, with the aim to give the right therapy to the right patient and avoid unnecessary overtreatment. The search for predictive biomarkers is an active field of research that includes genomic, proteomic and/or machine learning approaches. In this review, we describe current strategies and innovative tools to identify, evaluate and validate new biomarkers. We also summarize current predictive biomarkers in breast cancer and discuss companion biomarkers of targeted therapy in the context of precision medicine.

Interpretable Diagnosis for Whole-Slide Melanoma Histology Images Using Convolutional Neural Network

At present, deep learning-based medical image diagnosis had achieved high performance in several diseases. However, the black-box nature of the convolutional neural network (CNN) limits their role in diagnosis. In this study, a novel interpretable diagnosis pipeline using the CNN model was proposed. Furthermore, a sizeable melanoma database that contains 841 digital whole-slide images (WSIs) was built to train and evaluate the model. The model achieved strong melanoma classification ability (0.962 areas under the receiver operating characteristic, 0.887 sensitivity, and 0.925 specificity). Moreover, the proposed model outperformed the existing schemes in terms of accuracy that is 20 pathologists (0.933 vs 0.732 accuracy). Finally, the gradient-weighted class activation mapping (Grad-CAM) method was used to show the inner logic of the proposed model and its feasibility to improve diagnosis process in healthcare. The mechanism of feature heat maps which is visualized through a saliency mapping has demonstrated that features learned or extracted by the proposed model are compatible with the accepted pathological features. Conclusively, the proposed model provides a rapid and accurate diagnosis by locating the distinctive features of melanoma to build doctors' trust in the CNNs' diagnosis results.

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

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

Methods used for the measurement of blood-brain barrier integrity

The blood-brain barrier (BBB) comprises the interface between blood, brain and cerebrospinal fluid. Its primary function, which is mainly carried out by tight junctions, is to stabilize the tightly controlled microenvironment of the brain. To study the development and maintenance of the BBB, as well as various roles their intrinsic mechanisms that play in neurological disorders, suitable measurements are required to demonstrate integrity and functional changes at the interfaces between the blood and brain tissue. Markers and plasma proteins with different molecular weight (MW) are used to measure the permeability of BBB. In addition, the expression changes of tight-junction proteins form the basic structure of BBB, and imaging modalities are available to study the disruption of BBB. In the present review, above mentioned methods are depicted in details, together with the pros and cons as well as the differences between these methods, which maybe benefit research studies focused on the detection of BBB breakdown.

Transformational Role of Medical Imaging in (Radiation) Oncology

Simple Summary

Onboard, imaging techniques have brought about a huge transformation in the ability to deliver targeted radiation therapies. Each generation of these technologies enables us to better visualize where to deliver lethal doses of radiation and thus allows the shrinking of necessary geometric margins leading to reduced toxicities. Alongside improvements in treatment delivery, advances in medical imaging have also allowed us to better define the volumes we wish to target. The development of imaging techniques that can capture aspects of the tumor’s biology before, during and after therapy is transforming how treatment can be delivered. Technological changes have further made these biological imaging techniques available in real-time providing the opportunity to monitor a patient’s response to treatment closely and often before any volume changes are visible on conventional radiological images. Here we discuss the development of robust quantitative imaging biomarkers and how they can personalize therapy towards meaningful clinical endpoints.

Abstract

Onboard, real-time, imaging techniques, from the original megavoltage planar imaging devices, to the emerging combined MRI-Linear Accelerators, have brought a huge transformation in the ability to deliver targeted radiation therapies. Each generation of these technologies enables lethal doses of radiation to be delivered to target volumes with progressively more accuracy and thus allows shrinking of necessary geometric margins, leading to reduced toxicities. Alongside these improvements in treatment delivery, advances in medical imaging, e.g., PET, and MRI, have also allowed target volumes themselves to be better defined. The development of functional and molecular imaging is now driving a conceptually larger step transformation to both better understand the cancer target and disease to be treated, as well as how tumors respond to treatment. A biological description of the tumor microenvironment is now accepted as an essential component of how to personalize and adapt treatment. This applies not only to radiation oncology but extends widely in cancer management from surgical oncology planning and interventional radiology, to evaluation of targeted drug delivery efficacy in medical oncology/immunotherapy. Here, we will discuss the role and requirements of functional and metabolic imaging techniques in the context of brain tumors and metastases to reliably provide multi-parametric imaging biomarkers of the tumor microenvironment.

Theranostics: The Role of Quantitative Nuclear Medicine Imaging

Theranostics is a precision medicine discipline that integrates diagnostic nuclear medicine imaging with radionuclide therapy in a manner that provides both a tumor phenotype and personalized therapy to patients with cancer using radiopharmaceuticals aimed at the same target-specific biological pathway or receptor. The aim of quantitative nuclear medicine imaging is to plan the alpha or beta-emitting therapy based on an accurate 3-dimensional representation of the in-vivo distribution of radioactivity concentration within the tumor and normal organs/tissues in a noninvasive manner. In general, imaging may be either based on positron emission tomography (PET) or single photon emission computed tomography (SPECT) invariably in combination with X-ray CT (PET/CT; SPECT/CT) or, to a much lesser extent, MRI. PET and SPECT differ in terms of the radionuclides and physical processes that give rise to the emission of high energy photons, as well as the sets of technologies involved in their detection. Using a variety of standardized quantitative parameters, system calibration, patient preparation, imaging acquisition and reconstruction protocols, and image analysis protocols, an accurate quantification of the tracer distribution can be obtained, which helps prescribe the therapeutic dose for each patient.

Development and Validation of a Machine Learning Model to Explore Tyrosine Kinase Inhibitor Response in Patients With Stage IV EGFR Variant-Positive Non-Small Cell Lung Cancer

IMPORTANCE

An end-to-end efficacy evaluation approach for identifying progression risk after epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitor (TKI) therapy in patients with stage IV EGFR variant-positive non-small cell lung cancer (NSCLC) is lacking.

OBJECTIVE

To propose a clinically applicable large-scale bidirectional generative adversarial network for predicting the efficacy of EGFR-TKI therapy in patients with NSCLC.

DESIGN, SETTING, AND PARTICIPANTS

This diagnostic/prognostic study enrolled 465 patients from January 1, 2010, to August 1, 2017, with follow-up from February 1, 2010, to June 1, 2020. A deep learning (DL) semantic signature to predict progression-free survival (PFS) was constructed in the training cohort, validated in 2 external validation and 2 control cohorts, and compared with the radiomics signature.

EXPOSURES

An end-to-end bidirectional generative adversarial network framework was designed to predict the progression risk in patients with NSCLC.

MAIN OUTCOMES AND MEASURES

The primary end point was PFS, considering the time from the initiation of therapy to the date of recurrence, confirmed disease progression, or death.

RESULTS

A total of 342 patients with stage IV EGFR variant-positive NSCLC receiving EGFR-TKI therapy met the inclusion criteria. Of these, 145 patients from 2 of the hospitals (n = 117 and 28) formed a training cohort (mean [SD] age, 61 [11] years; 87 [60.0%] female), and the patients from 2 other hospitals comprised 2 external validation cohorts (validation cohort 1: n = 101; mean [SD] age, 57 [12] years; 60 [59.4%] female; and validation cohort 2: n = 96, mean [SD] age, 58 [9] years; 55 [57.3%] female). Fifty-six patients with advanced-stage EGFR variant-positive NSCLC (mean [SD] age, 52 [11] years; 26 [46.4%] female) and 67 patients with advanced-stage EGFR wild-type NSCLC (mean [SD] age, 54 [10] years; 10 [15.0%] female) who received first-line chemotherapy were included. A total of 90 (26%) receiving EGFR-TKI therapy with a high risk of rapid disease progression were identified (median [range] PFS, 7.3 [1.4-32.0] months in the training cohort, 5.0 [0.6-34.6] months in validation cohort 1, and 6.4 [1.8-20.1] months, in validation cohort 2) using the DL semantic signature.The PFS decreased by 36% (hazard ratio, 2.13; 95% CI, 1.30-3.49; P < .001) compared with that in other patients (median [range] PFS, 11.5 [1.5-64.2] months in the training cohort, 10.9 [1.1-50.5] in validation cohort 1, and 8.9 [0.8-40.6] months in validation cohort 2. No significant differences were observed when comparing the PFS of high-risk patients receiving EGFR-TKI therapy with the chemotherapy cohorts (median PFS, 6.9 vs 4.4 months; P = .08). In terms of predicting the tumor progression risk after EGFR-TKI therapy, clinical decisions based on the DL semantic signature led to better survival outcomes than those based on radiomics signature across all risk probabilities by the decision curve analysis.

CONCLUSIONS AND RELEVANCE

This diagnostic/prognostic study provides a clinically applicable approach for identifying patients with stage IV EGFR variant-positive NSCLC who are not likely to benefit from EGFR-TKI therapy. The end-to-end DL-derived semantic features eliminated all manual interventions required while using previous radiomics methods and have a better prognostic performance.

Integrative Models of Histopathological Image Features and Omics Data Predict Survival in Head and Neck Squamous Cell Carcinoma

Background

Both histopathological image features and genomics data were associated with survival outcome of cancer patients. However, integrating features of histopathological images, genomics and other omics for improving prognosis prediction has not been reported in head and neck squamous cell carcinoma (HNSCC).

Methods

A dataset of 216 HNSCC patients was derived from the Cancer Genome Atlas (TCGA) with information of clinical characteristics, genetic mutation, RNA sequencing, protein expression and histopathological images. Patients were randomly assigned into training (n = 108) or validation (n = 108) sets. We extracted 593 quantitative image features, and used random forest algorithm with 10-fold cross-validation to build prognostic models for overall survival (OS) in training set, then compared the area under the time-dependent receiver operating characteristic curve (AUC) in validation set.

Results

In validation set, histopathological image features had significant predictive value for OS (5-year AUC = 0.784). The histopathology + omics models showed better predictive performance than genomics, transcriptomics or proteomics alone. Moreover, the multi-omics model incorporating image features, genomics, transcriptomics and proteomics reached the maximal 1-, 3-, and 5-year AUC of 0.871, 0.908, and 0.929, with most significant survival difference (HR = 10.66, 95%CI: 5.06–26.8, p < 0.001). Decision curve analysis also revealed a better net benefit of multi-omics model.

Conclusion

The histopathological images could provide complementary features to improve prognostic performance for HNSCC patients. The integrative model of histopathological image features and omics data might serve as an effective tool for survival prediction and risk stratification in clinical practice.

18F-FDG PET and DCE kinetic modeling and their correlations in primary NSCLC: first voxel-wise correlative analysis of human simultaneous [18F]FDG PET-MRI data

OBJECTIVES

To decipher the correlations between PET and DCE kinetic parameters in non-small-cell lung cancer (NSCLC), by using voxel-wise analysis of dynamic simultaneous [18F]FDG PET-MRI.

MATERIAL AND METHODS

Fourteen treatment-naïve patients with biopsy-proven NSCLC prospectively underwent a 1-h dynamic [18F]FDG thoracic PET-MRI scan including DCE. The PET and DCE data were normalized to their corresponding T1-weighted MR morphological space, and tumors were masked semi-automatically. Voxel-wise parametric maps of PET and DCE kinetic parameters were computed by fitting the dynamic PET and DCE tumor data to the Sokoloff and Extended Tofts models respectively, by using in-house developed procedures. Curve-fitting errors were assessed by computing the relative root mean square error (rRMSE) of the estimated PET and DCE signals at the voxel level. For each tumor, Spearman correlation coefficients (rs) between all the pairs of PET and DCE kinetic parameters were estimated on a voxel-wise basis, along with their respective bootstrapped 95% confidence intervals (n = 1000 iterations).

RESULTS

Curve-fitting metrics provided fit errors under 20% for almost 90% of the PET voxels (median rRMSE = 10.3, interquartile ranges IQR = 8.1; 14.3), whereas 73.3% of the DCE voxels showed fit errors under 45% (median rRMSE = 31.8%, IQR = 22.4; 46.6). The PET-PET, DCE-DCE, and PET-DCE voxel-wise correlations varied according to individual tumor behaviors. Beyond this wide variability, the PET-PET and DCE-DCE correlations were mainly high (absolute rs values > 0.7), whereas the PET-DCE correlations were mainly low to moderate (absolute rs values < 0.7). Half the tumors showed a hypometabolism with low perfused/vascularized profile, a hallmark of hypoxia, and tumor aggressiveness.

CONCLUSION

A dynamic "one-stop shop" procedure applied to NSCLC is technically feasible in clinical practice. PET and DCE kinetic parameters assessed simultaneously are not highly correlated in NSCLC, and these correlations showed a wide variability among tumors and patients. These results tend to suggest that PET and DCE kinetic parameters might provide complementary information. In the future, this might make PET-MRI a unique tool to characterize the individual tumor biological behavior in NSCLC.

ePAD: An Image Annotation and Analysis Platform for Quantitative Imaging

Medical imaging is critical for assessing the response of patients to new cancer therapies. Quantitative lesion assessment on images is time-consuming, and adopting new promising quantitative imaging biomarkers of response in clinical trials is challenging. The electronic Physician Annotation Device (ePAD) is a freely available web-based zero-footprint software application for viewing, annotation, and quantitative analysis of radiology images designed to meet the challenges of quantitative evaluation of cancer lesions. For imaging researchers, ePAD calculates a variety of quantitative imaging biomarkers that they can analyze and compare in ePAD to identify potential candidates as surrogate endpoints in clinical trials. For clinicians, ePAD provides clinical decision support tools for evaluating cancer response through reports summarizing changes in tumor burden based on different imaging biomarkers. As a workflow management and study oversight tool, ePAD lets clinical trial project administrators create worklists for users and oversee the progress of annotations created by research groups. To support interoperability of image annotations, ePAD writes all image annotations and results of quantitative imaging analyses in standardized file formats, and it supports migration of annotations from various propriety formats. ePAD also provides a plugin architecture supporting MATLAB server-side modules in addition to client-side plugins, permitting the community to extend the ePAD platform in various ways for new cancer use cases. We present an overview of ePAD as a platform for medical image annotation and quantitative analysis. We also discuss use cases and collaborations with different groups in the Quantitative Imaging Network and future directions.

Application of in vivo MR methods in the study of breast cancer metabolism

In the last two decades, various in vivo MR methodologies have been evaluated for their potential in the study of cancer metabolism. During malignant transformation, metabolic alterations occur, leading to morphological and functional changes. Among various MR methods, in vivo MRS has been extensively used in breast cancer to study the metabolism of cells, tissues or whole organs. It provides biochemical information at the metabolite level. Altered choline, phospholipid and energy metabolism has been documented using proton (¹ H), phosphorus (³¹ P) and carbon (¹³ C) isotopes. Increased levels of choline-containing compounds, phosphomonoesters and phosphodiesters in breast cancer, which are indicative of altered choline and phospholipid metabolism, have been reported using in vivo, in vitro and ex vivo NMR studies. These changes are reversed on successful therapy, which depends on the treatment regimen given. Monitoring the various tumor intermediary metabolic pathways using nuclear spin hyperpolarization of ¹³ C-labeled substrates by dynamic nuclear polarization has also been recently reported. Furthermore, the utility of various methods such as diffusion, dynamic contrast and perfusion MRI have also been evaluated to study breast tumor metabolism. Parameters such as tumor volume, apparent diffusion coefficient, volume transfer coefficient and extracellular volume ratio are estimated. These parameters provide information on the changes in tumor microstructure, microenvironment, abnormal vasculature, permeability and grade of the tumor. Such changes seen during cancer progression are due to alterations in the tumor metabolism, leading to changes in cell architecture. Due to architectural changes, the tissue mechanical properties are altered; this can be studied using magnetic resonance elastography, which measures the elastic properties of tissues. Moreover, these structural MRI methods can be used to investigate the effect of therapy-induced changes in tumor characteristics. This review discusses the potential of various in vivo MR methodologies in the study of breast cancer metabolism.

Low perfusion compartments in glioblastoma quantified by advanced magnetic resonance imaging and correlated with patient survival

BACKGROUND AND PURPOSE

Glioblastoma exhibits profound intratumoral heterogeneity in perfusion. Particularly, low perfusion may induce treatment resistance. Thus, imaging approaches that define low perfusion compartments are crucial for clinical management.

MATERIALS AND METHODS

A total of 112 newly diagnosed glioblastoma patients were prospectively recruited for maximal safe resection. The apparent diffusion coefficient (ADC) and relative cerebral blood volume (rCBV) were calculated from diffusion and perfusion imaging, respectively. Based on the overlapping regions of lowest rCBV quartile (rCBVL) with the lowest ADC quartile (ADCL) and highest ADC quartile (ADCH) in each tumor, two low perfusion compartments (ADCH-rCBVL and ADCL-rCBVL) were identified for volumetric analysis. Lactate and macromolecule/lipid levels were determined from multivoxel MR spectroscopic imaging. Progression-free survival (PFS) and overall survival (OS) were analyzed using Kaplan-Meier's and multivariate Cox regression analyses, to evaluate the effects of compartment volume and lactate level on survival.

RESULTS

Two compartments displayed higher lactate and macromolecule/lipid levels compared to contralateral normal-appearing white matter (each P < 0.001). The proportion of the ADCL-rCBVL compartment in the contrast-enhancing tumor was associated with a larger infiltration on FLAIR (P < 0.001, rho = 0.42). The minimally invasive phenotype displayed a lower proportion of the ADCL-rCBVL compartment than the localized (P = 0.031) and diffuse phenotypes (not significant). Multivariate Cox regression showed higher lactate level in the ADCL-rCBVL compartment was associated with worsened survival (PFS: HR 2.995, P = 0.047; OS: HR 4.974, P = 0.005).

CONCLUSIONS

Our results suggest that the ADCL-rCBVL compartment may potentially indicate a clinically measurable resistant compartment.

Imaging tumour hypoxia with oxygen-enhanced MRI and BOLD MRI

Hypoxia is known to be a poor prognostic indicator for nearly all solid tumours and also is predictive of treatment failure for radiotherapy, chemotherapy, surgery and targeted therapies. Imaging has potential to identify, spatially map and quantify tumour hypoxia prior to therapy, as well as track changes in hypoxia on treatment. At present no hypoxia imaging methods are available for routine clinical use. Research has largely focused on positron emission tomography (PET)-based techniques, but there is gathering evidence that MRI techniques may provide a practical and more readily translational alternative. In this review we focus on the potential for imaging hypoxia by measuring changes in longitudinal relaxation [R1; termed oxygen-enhanced MRI or tumour oxygenation level dependent (TOLD) MRI] and effective transverse relaxation [R2*; termed blood oxygenation level dependent (BOLD) MRI], induced by inhalation of either 100% oxygen or the radiosensitising hyperoxic gas carbogen. We explain the scientific principles behind oxygen-enhanced MRI and BOLD and discuss significant studies and their limitations. All imaging biomarkers require rigorous validation in order to translate into clinical use and the steps required to further develop oxygen-enhanced MRI and BOLD MRI into decision-making tools are discussed.

Mass Spectrometry-Based Biomarkers in Drug Development

Advances in mass spectrometry, proteomics, protein bioanalytical approaches, and biochemistry have led to a rapid evolution and expansion in the area of mass spectrometry-based biomarker discovery and development. The last decade has also seen significant progress in establishing accepted definitions, guidelines, and criteria for the analytical validation, acceptance and qualification of biomarkers. These advances have coincided with a decreased return on investment for pharmaceutical research and development and an increasing need for better early decision making tools. Empowering development teams with tools to measure a therapeutic interventions impact on disease state and progression, measure target engagement and to confirm predicted pharmacodynamic effects is critical to efficient data-driven decision making. Appropriate implementation of a biomarker or a combination of biomarkers can enhance understanding of a drugs mechanism, facilitate effective translation from the preclinical to clinical space, enable early proof of concept and dose selection, and increases the efficiency of drug development. Here we will provide descriptions of the different classes of biomarkers that have utility in the drug development process as well as review specific, protein-centric, mass spectrometry-based approaches for the discovery of biomarkers and development of targeted assays to measure these markers in a selective and analytically precise manner.

Predicting survival from colorectal cancer histology slides using deep learning: A retrospective multicenter study

BACKGROUND

For virtually every patient with colorectal cancer (CRC), hematoxylin-eosin (HE)-stained tissue slides are available. These images contain quantitative information, which is not routinely used to objectively extract prognostic biomarkers. In the present study, we investigated whether deep convolutional neural networks (CNNs) can extract prognosticators directly from these widely available images.

METHODS AND FINDINGS

We hand-delineated single-tissue regions in 86 CRC tissue slides, yielding more than 100,000 HE image patches, and used these to train a CNN by transfer learning, reaching a nine-class accuracy of >94% in an independent data set of 7,180 images from 25 CRC patients. With this tool, we performed automated tissue decomposition of representative multitissue HE images from 862 HE slides in 500 stage I-IV CRC patients in the The Cancer Genome Atlas (TCGA) cohort, a large international multicenter collection of CRC tissue. Based on the output neuron activations in the CNN, we calculated a "deep stroma score," which was an independent prognostic factor for overall survival (OS) in a multivariable Cox proportional hazard model (hazard ratio [HR] with 95% confidence interval [CI]: 1.99 [1.27-3.12], p = 0.0028), while in the same cohort, manual quantification of stromal areas and a gene expression signature of cancer-associated fibroblasts (CAFs) were only prognostic in specific tumor stages. We validated these findings in an independent cohort of 409 stage I-IV CRC patients from the "Darmkrebs: Chancen der Verhütung durch Screening" (DACHS) study who were recruited between 2003 and 2007 in multiple institutions in Germany. Again, the score was an independent prognostic factor for OS (HR 1.63 [1.14-2.33], p = 0.008), CRC-specific OS (HR 2.29 [1.5-3.48], p = 0.0004), and relapse-free survival (RFS; HR 1.92 [1.34-2.76], p = 0.0004). A prospective validation is required before this biomarker can be implemented in clinical workflows.

CONCLUSIONS

In our retrospective study, we show that a CNN can assess the human tumor microenvironment and predict prognosis directly from histopathological images.

Plasma Tie2 is a tumor vascular response biomarker for VEGF inhibitors in metastatic colorectal cancer

Oncological use of anti-angiogenic VEGF inhibitors has been limited by the lack of informative biomarkers. Previously we reported circulating Tie2 as a vascular response biomarker for bevacizumab-treated ovarian cancer patients. Using advanced MRI and circulating biomarkers we have extended these findings in metastatic colorectal cancer (n = 70). Bevacizumab (10 mg/kg) was administered to elicit a biomarker response, followed by FOLFOX6-bevacizumab until disease progression. Bevacizumab induced a correlation between Tie2 and the tumor vascular imaging biomarker, Ktrans (R:-0.21 to 0.47) implying that Tie2 originated from the tumor vasculature. Tie2 trajectories were independently associated with pre-treatment tumor vascular characteristics, tumor response, progression free survival (HR for progression = 3.01, p = 0.00014; median PFS 248 vs. 348 days p = 0.0008) and the modeling of progressive disease (p < 0.0001), suggesting that Tie2 should be monitored clinically to optimize VEGF inhibitor use. A vascular response is defined as a 30% reduction in Tie2; vascular progression as a 40% increase in Tie2 above the nadir. Tie2 is the first, validated, tumor vascular response biomarker for VEGFi.

Role of texture analysis in breast MRI as a cancer biomarker: A review

Breast cancer is a known heterogeneous disease. Current clinically utilized histopathologic biomarkers may undersample tumor heterogeneity, resulting in higher rates of misdiagnosis for breast cancer. MRI can provide a whole-tumor sampling of disease burden and is widely utilized in clinical care. Texture analysis can provide a localized description of breast cancer, with particular emphasis on quantifying breast lesion heterogeneity. The object of this review is to provide an overview of texture analysis applications towards breast cancer diagnosis, prognosis, and treatment response evaluation and review the role of image-based texture features as noninvasive prognostic and predictive biomarkers. Level of Evidence: 5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:927-938.

Delivering Functional Imaging on the MRI-Linac: Current Challenges and Potential Solutions

Magnetic resonance imaging (MRI) is a highly versatile imaging modality that can be used to measure features of the tumour microenvironment including cell death, proliferation, metabolism, angiogenesis, and hypoxia. Mapping and quantifying these pathophysiological features has the potential to alter the use of adaptive radiotherapy planning. Although these methods are available for use on diagnostic machines, several challenges exist for implementing these functional MRI methods on the MRI-linear accelerators (linacs). This review considers these challenges and potential solutions.

Novel biomarkers for patient stratification in colorectal cancer: A review of definitions, emerging concepts, and data

Colorectal cancer (CRC) treatment has become more personalised, incorporating a combination of the individual patient risk assessment, gene testing, and chemotherapy with surgery for optimal care. The improvement of staging with high-resolution imaging has allowed more selective treatments, optimising survival outcomes. The next step is to identify biomarkers that can inform clinicians of expected prognosis and offer the most beneficial treatment, while reducing unnecessary morbidity for the patient. The search for biomarkers in CRC has been of significant interest, with questions remaining on their impact and applicability. The study of biomarkers can be broadly divided into metabolic, molecular, microRNA, epithelial-to-mesenchymal-transition (EMT), and imaging classes. Although numerous molecules have claimed to impact prognosis and treatment, their clinical application has been limited. Furthermore, routine testing of prognostic markers with no demonstrable influence on response to treatment is a questionable practice, as it increases cost and can adversely affect expectations of treatment. In this review we focus on recent developments and emerging biomarkers with potential utility for clinical translation in CRC. We examine and critically appraise novel imaging and molecular-based approaches; evaluate the promising array of microRNAs, analyze metabolic profiles, and highlight key findings for biomarker potential in the EMT pathway.

PET/MRI Hybrid Systems

Over the last decade, the combination of PET and MRI in one system has proven to be highly successful in basic preclinical research, as well as in clinical research. Nowadays, PET/MRI systems are well established in preclinical imaging and are progressing into clinical applications to provide further insights into specific diseases, therapeutic assessments, and biological pathways. Certain challenges in terms of hardware had to be resolved concurrently with the development of new techniques to be able to reach the full potential of both combined techniques. This review provides an overview of these challenges and describes the opportunities that simultaneous PET/MRI systems can exploit in comparison with stand-alone or other combined hybrid systems. New approaches were developed for simultaneous PET/MRI systems to correct for attenuation of 511 keV photons because MRI does not provide direct information on gamma photon attenuation properties. Furthermore, new algorithms to correct for motion were developed, because MRI can accurately detect motion with high temporal resolution. The additional information gained by the MRI can be employed to correct for partial volume effects as well. The development of new detector designs in combination with fast-decaying scintillator crystal materials enabled time-of-flight detection and incorporation in the reconstruction algorithms. Furthermore, this review lists the currently commercially available systems both for preclinical and clinical imaging and provides an overview of applications in both fields. In this regard, special emphasis has been placed on data analysis and the potential for both modalities to evolve with advanced image analysis tools, such as cluster analysis and machine learning.

A New Approach to Predict Progression-free Survival in Stage IV EGFR-mutant NSCLC Patients with EGFR-TKI Therapy

Purpose: We established a CT-derived approach to achieve accurate progression-free survival (PFS) prediction to EGFR tyrosine kinase inhibitors (TKI) therapy in multicenter, stage IV EGFR-mutated non-small cell lung cancer (NSCLC) patients.Experimental Design: A total of 1,032 CT-based phenotypic characteristics were extracted according to the intensity, shape, and texture of NSCLC pretherapy images. On the basis of these CT features extracted from 117 stage IV EGFR-mutant NSCLC patients, a CT-based phenotypic signature was proposed using a Cox regression model with LASSO penalty for the survival risk stratification of EGFR-TKI therapy. The signature was validated using two independent cohorts (101 and 96 patients, respectively). The benefit of EGFR-TKIs in stratified patients was then compared with another stage-IV EGFR-mutant NSCLC cohort only treated with standard chemotherapy (56 patients). Furthermore, an individualized prediction model incorporating the phenotypic signature and clinicopathologic risk characteristics was proposed for PFS prediction, and also validated by multicenter cohorts.Results: The signature consisted of 12 CT features demonstrated good accuracy for discriminating patients with rapid and slow progression to EGFR-TKI therapy in three cohorts (HR: 3.61, 3.77, and 3.67, respectively). Rapid progression patients received EGFR TKIs did not show significant difference with patients underwent chemotherapy for progression-free survival benefit (P = 0.682). Decision curve analysis revealed that the proposed model significantly improved the clinical benefit compared with the clinicopathologic-based characteristics model (P < 0.0001).Conclusions: The proposed CT-based predictive strategy can achieve individualized prediction of PFS probability to EGFR-TKI therapy in NSCLCs, which holds promise of improving the pretherapy personalized management of TKIs. Clin Cancer Res; 24(15); 3583-92. ©2018 AACR.

A Novel Unsupervised Segmentation Approach Quantifies Tumor Tissue Populations Using Multiparametric MRI: First Results with Histological Validation

Purpose

We aimed to precisely estimate intra-tumoral heterogeneity using spatially regularized spectral clustering (SRSC) on multiparametric MRI data and compare the efficacy of SRSC with the previously reported segmentation techniques in MRI studies.

Procedures

Six NMRI nu/nu mice bearing subcutaneous human glioblastoma U87 MG tumors were scanned using a dedicated small animal 7T magnetic resonance imaging (MRI) scanner. The data consisted of T2 weighted images, apparent diffusion coefficient maps, and pre- and post-contrast T2 and T2* maps. Following each scan, the tumors were excised into 2–3-mm thin slices parallel to the axial field of view and processed for histological staining. The MRI data were segmented using SRSC, K-means, fuzzy C-means, and Gaussian mixture modeling to estimate the fractional population of necrotic, peri-necrotic, and viable regions and validated with the fractional population obtained from histology.

Results

While the aforementioned methods overestimated peri-necrotic and underestimated viable fractions, SRSC accurately predicted the fractional population of all three tumor tissue types and exhibited strong correlations (r_necrotic = 0.92, r_{peri-necrotic} = 0.82 and r_viable = 0.98) with the histology.

Conclusions

The precise identification of necrotic, peri-necrotic and viable areas using SRSC may greatly assist in cancer treatment planning and add a new dimension to MRI-guided tumor biopsy procedures.

Electronic supplementary material

The online version of this article (doi:10.1007/s11307-016-1009-y) contains supplementary material, which is available to authorized users.

Quantitative radiomic profiling of glioblastoma represents transcriptomic expression

Quantitative imaging biomarkers have increasingly emerged in the field of research utilizing available imaging modalities. We aimed to identify good surrogate radiomic features that can represent genetic changes of tumors, thereby establishing noninvasive means for predicting treatment outcome. From May 2012 to June 2014, we retrospectively identified 65 patients with treatment-naïve glioblastoma with available clinical information from the Samsung Medical Center data registry. Preoperative MR imaging data were obtained for all 65 patients with primary glioblastoma. A total of 82 imaging features including first-order statistics, volume, and size features, were semi-automatically extracted from structural and physiologic images such as apparent diffusion coefficient and perfusion images. Using commercially available software, NordicICE, we performed quantitative imaging analysis and collected the dataset composed of radiophenotypic parameters. Unsupervised clustering methods revealed that the radiophenotypic dataset was composed of three clusters. Each cluster represented a distinct molecular classification of glioblastoma; classical type, proneural and neural types, and mesenchymal type. These clusters also reflected differential clinical outcomes. We found that extracted imaging signatures does not represent copy number variation and somatic mutation. Quantitative radiomic features provide a potential evidence to predict molecular phenotype and treatment outcome. Radiomic profiles represents transcriptomic phenotypes more well.

Whole-body MRI quantitative biomarkers are associated significantly with treatment response in patients with newly diagnosed symptomatic multiple myeloma following bortezomib induction

OBJECTIVES

To evaluate whole-body MRI (WB-MRI) parameters significantly associated with treatment response in multiple myeloma (MM).

METHODS

Twenty-one MM patients underwent WB-MRI at diagnosis and after two cycles of chemotherapy. Scans acquired at 3.0 T included T2, diffusion-weighted-imaging (DWI) and mDixon pre- and post-contrast. Twenty focal lesions (FLs) matched on DWI and post-contrast mDixon were selected for each time point. Estimated tumour volume (eTV), apparent diffusion coefficient (ADC), enhancement ratio (ER) and signal fat fraction (sFF) were derived. Clinical treatment response to chemotherapy was assessed using conventional criteria. Significance of temporal parameter change was assessed by the paired t test and receiver operating characteristics/area under the curve (AUC) analysis was performed. Parameter repeatability was assessed by interclass correlation (ICC) and Bland-Altman analysis of 10 healthy volunteers scanned at two time points.

RESULTS

Fifteen of 21 patients responded to treatment. Of 254 FLs analysed, sFF (p < 0.0001) and ADC (p = 0.001) significantly increased in responders but not non-responders. eTV significantly decreased in 19/21 cases. Focal lesion sFF was the best discriminator of treatment response (AUC 1.0). Bone sFF repeatability was excellent (ICC 0.98) and better than bone ADC (ICC 0.47).

CONCLUSION

WB-MRI derived focal lesion sFF shows promise as an imaging biomarker of treatment response in newly diagnosed MM.

KEY POINTS

• Bone signal fat fraction using mDixon is a robust quantifiable parameter • Fat fraction and ADC significantly increase in myeloma lesions responding to treatment • Bone lesion fat fraction is the best discriminator of myeloma treatment response.

The potential of hyperpolarized 13C magnetic resonance spectroscopy to monitor the effect of combretastatin based vascular disrupting agents

BACKGROUND

Targeting tumor vasculature with vascular disrupting agents (VDAs) results in substantial cell death that precede tumor shrinkage. Here, we investigate the potential of hyperpolarized magnetic resonance spectroscopy (HPMRS) to monitor early metabolic changes associated with VDA treatment.

METHODS

Mice bearing C3H mammary carcinomas were treated with the VDAs combretastatin-A4-phosphate (CA4P) or the analog OXi4503, and HPMRS was performed following [1-¹³C]pyruvate administration. Similarly, treated mice were positron emission tomography (PET) scanned following administration of the glucose analog FDG. Finally, metabolic imaging parameters were compared to tumor regrowth delay and measures of vascular damage, derived from dynamic contrast-agent enhanced magnetic resonance imaging (DCE-MRI) and histology.

RESULTS

VDA-treatment impaired tumor perfusion (histology and DCE-MRI), reduced FDG uptake, increased necrosis, and slowed tumor growth. HPMRS, revealed that the [1-¹³C]pyruvate-to-[1-¹³C]lactate conversion remained unaltered, whereas [1-¹³C]lactate-to-[¹³C]bicarbonate (originating from respiratory CO₂) ratios increased significantly following treatment.

CONCLUSIONS

DCE-MRI and FDG-PET revealed loss of vessel functionality, impaired glucose delivery and reduced metabolic activity prior to cell death. [1-¹³C]lactate-to-[¹³C]bicarbonate ratios increased significantly during treatment, indicating a decline in respiratory activity driven by the onset of hypoxia. HPMRS is promising for early detection of metabolic stress inflicted by VDAs, which cannot easily be inferred based on blood flow measurements.

Cancer Metabolism and Tumor Heterogeneity: Imaging Perspectives Using MR Imaging and Spectroscopy

Cancer cells reprogram their metabolism to maintain viability via genetic mutations and epigenetic alterations, expressing overall dynamic heterogeneity. The complex relaxation mechanisms of nuclear spins provide unique and convertible tissue contrasts, making magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) pertinent imaging tools in both clinics and research. In this review, we summarized MR methods that visualize tumor characteristics and its metabolic phenotypes on an anatomical, microvascular, microstructural, microenvironmental, and metabolomics scale. The review will progress from the utilities of basic spin-relaxation contrasts in cancer imaging to more advanced imaging methods that measure tumor-distinctive parameters such as perfusion, water diffusion, magnetic susceptibility, oxygenation, acidosis, redox state, and cell death. Analytical methods to assess tumor heterogeneity are also reviewed in brief. Although the clinical utility of tumor heterogeneity from imaging is debatable, the quantification of tumor heterogeneity using functional and metabolic MR images with development of robust analytical methods and improved MR methods may offer more critical roles of tumor heterogeneity data in clinics. MRI/MRS can also provide insightful information on pharmacometabolomics, biomarker discovery, disease diagnosis and prognosis, and treatment response. With these future directions in mind, we anticipate the widespread utilization of these MR-based techniques in studying in vivo cancer biology to better address significant clinical needs.

Radiomics: the bridge between medical imaging and personalized medicine

Radiomics, the high-throughput mining of quantitative image features from standard-of-care medical imaging that enables data to be extracted and applied within clinical-decision support systems to improve diagnostic, prognostic, and predictive accuracy, is gaining importance in cancer research. Radiomic analysis exploits sophisticated image analysis tools and the rapid development and validation of medical imaging data that uses image-based signatures for precision diagnosis and treatment, providing a powerful tool in modern medicine. Herein, we describe the process of radiomics, its pitfalls, challenges, opportunities, and its capacity to improve clinical decision making, emphasizing the utility for patients with cancer. Currently, the field of radiomics lacks standardized evaluation of both the scientific integrity and the clinical relevance of the numerous published radiomics investigations resulting from the rapid growth of this area. Rigorous evaluation criteria and reporting guidelines need to be established in order for radiomics to mature as a discipline. Herein, we provide guidance for investigations to meet this urgent need in the field of radiomics.

Cancer heterogeneity and imaging

There is interest in identifying and quantifying tumor heterogeneity at the genomic, tissue pathology and clinical imaging scales, as this may help better understand tumor biology and may yield useful biomarkers for guiding therapy-based decision making. This review focuses on the role and value of using x-ray, CT, MRI and PET based imaging methods that identify, measure and map tumor heterogeneity. In particular we highlight the potential value of these techniques and the key challenges required to validate and qualify these biomarkers for clinical use.

Imaging biomarker roadmap for cancer studies

Imaging biomarkers (IBs) are integral to the routine management of patients with cancer. IBs used daily in oncology include clinical TNM stage, objective response and left ventricular ejection fraction. Other CT, MRI, PET and ultrasonography biomarkers are used extensively in cancer research and drug development. New IBs need to be established either as useful tools for testing research hypotheses in clinical trials and research studies, or as clinical decision-making tools for use in healthcare, by crossing 'translational gaps' through validation and qualification. Important differences exist between IBs and biospecimen-derived biomarkers and, therefore, the development of IBs requires a tailored 'roadmap'. Recognizing this need, Cancer Research UK (CRUK) and the European Organisation for Research and Treatment of Cancer (EORTC) assembled experts to review, debate and summarize the challenges of IB validation and qualification. This consensus group has produced 14 key recommendations for accelerating the clinical translation of IBs, which highlight the role of parallel (rather than sequential) tracks of technical (assay) validation, biological/clinical validation and assessment of cost-effectiveness; the need for IB standardization and accreditation systems; the need to continually revisit IB precision; an alternative framework for biological/clinical validation of IBs; and the essential requirements for multicentre studies to qualify IBs for clinical use.

Decision support systems for personalized and participative radiation oncology

A paradigm shift from current population based medicine to personalized and participative medicine is underway. This transition is being supported by the development of clinical decision support systems based on prediction models of treatment outcome. In radiation oncology, these models 'learn' using advanced and innovative information technologies (ideally in a distributed fashion - please watch the animation: http://youtu.be/ZDJFOxpwqEA) from all available/appropriate medical data (clinical, treatment, imaging, biological/genetic, etc.) to achieve the highest possible accuracy with respect to prediction of tumor response and normal tissue toxicity. In this position paper, we deliver an overview of the factors that are associated with outcome in radiation oncology and discuss the methodology behind the development of accurate prediction models, which is a multi-faceted process. Subsequent to initial development/validation and clinical introduction, decision support systems should be constantly re-evaluated (through quality assurance procedures) in different patient datasets in order to refine and re-optimize the models, ensuring the continuous utility of the models. In the reasonably near future, decision support systems will be fully integrated within the clinic, with data and knowledge being shared in a standardized, dynamic, and potentially global manner enabling truly personalized and participative medicine.

[Quantitative perfusion imaging in magnetic resonance imaging]

CLINICAL/METHODICAL ISSUE

Magnetic resonance imaging (MRI) is recognized for its superior tissue contrast while being non-invasive and free of ionizing radiation. Due to the development of new scanner hardware and fast imaging techniques during the last decades, access to tissue and organ functions became possible. One of these functional imaging techniques is perfusion imaging with which tissue perfusion and capillary permeability can be determined from dynamic imaging data.

STANDARD RADIOLOGICAL METHODS

Perfusion imaging by MRI can be performed by two approaches, arterial spin labeling (ASL) and dynamic contrast-enhanced (DCE) MRI. While the first method uses magnetically labelled water protons in arterial blood as an endogenous tracer, the latter involves the injection of a contrast agent, usually gadolinium (Gd), as a tracer for calculating hemodynamic parameters.

PERFORMANCE

Studies have demonstrated the potential of perfusion MRI for diagnostics and also for therapy monitoring.

ACHIEVEMENTS

The utilization and application of perfusion MRI are still restricted to specialized centers, such as university hospitals. A broad application of the technique has not yet been implemented.

PRACTICAL RECOMMENDATIONS

The MRI perfusion technique is a valuable tool that might come broadly available after implementation of standards on European and international levels. Such efforts are being promoted by the respective professional bodies.

Acute tumour response to a bispecific Ang-2-VEGF-A antibody: insights from multiparametric MRI and gene expression profiling

BACKGROUND

To assess antivascular effects, and evaluate clinically translatable magnetic resonance imaging (MRI) biomarkers of tumour response in vivo, following treatment with vanucizumab, a bispecific human antibody against angiopoietin-2 (Ang-2) and vascular endothelial growth factor-A (VEGF-A).

METHODS

Colo205 colon cancer xenografts were imaged before and 5 days after treatment with a single 10 mg kg(-1) dose of either vanucizumab, bevacizumab (anti-human VEGF-A), LC06 (anti-murine/human Ang-2) or omalizumab (anti-human IgE control). Volumetric response was assessed using T2-weighted MRI, and diffusion-weighted, dynamic contrast-enhanced (DCE) and susceptibility contrast MRI used to quantify tumour water diffusivity (apparent diffusion coefficient (ADC), × 10(6) mm(2) s(-1)), vascular perfusion/permeability (K(trans), min(-1)) and fractional blood volume (fBV, %) respectively. Pathological correlates were sought, and preliminary gene expression profiling performed.

RESULTS

Treatment with vanucizumab, bevacizumab or LC06 induced a significant (P<0.01) cytolentic response compared with control. There was no significant change in tumour ADC in any treatment group. Uptake of Gd-DTPA was restricted to the tumour periphery in all post-treatment groups. A significant reduction in tumour K(trans) (P<0.05) and fBV (P<0.01) was determined 5 days after treatment with vanucizumab only. This was associated with a significant (P<0.05) reduction in Hoechst 33342 uptake compared with control. Gene expression profiling identified 20 human genes exclusively regulated by vanucizumab, 6 of which are known to be involved in vasculogenesis and angiogenesis.

CONCLUSIONS

Vanucizumab is a promising antitumour and antiangiogenic treatment, whose antivascular activity can be monitored using DCE and susceptibility contrast MRI. Differential gene expression in vanucizumab-treated tumours is regulated by the combined effect of Ang-2 and VEGF-A inhibition.

Geometrical Measures Obtained from Pretreatment Postcontrast T1 Weighted MRIs Predict Survival Benefits from Bevacizumab in Glioblastoma Patients

Background

Antiangiogenic therapies for glioblastoma (GBM) such as bevacizumab (BVZ), have been unable to extend survival in large patient cohorts. However, a subset of patients having angiogenesis-dependent tumors might benefit from these therapies. Currently, there are no biomarkers allowing to discriminate responders from non-responders before the start of the therapy.

Methods

40 patients from the randomized GENOM009 study complied the inclusion criteria (quality of images, clinical data available). Of those, 23 patients received first line temozolomide (TMZ) for eight weeks and then concomitant radiotherapy and TMZ. 17 patients received BVZ+TMZ for seven weeks and then added radiotherapy to the treatment. Clinical variables were collected, tumors segmented and several geometrical measures computed including: Contrast enhancing (CE), necrotic, and total volumes; equivalent spherical CE width; several geometric measures of the CE ‘rim’ geometry and a set of image texture measures. The significance of the results was studied using Kaplan-Meier and Cox proportional hazards analysis. Correlations were assessed using Spearman correlation coefficients.

Results

Kaplan-Meier and Cox proportional hazards analysis showed that total, CE and inner volume (p = 0.019, HR = 4.258) and geometric heterogeneity of the CE areas (p = 0.011, HR = 3.931) were significant parameters identifying response to BVZ. The group of patients with either regular CE areas (small geometric heterogeneity, median difference survival 15.88 months, p = 0.011) or those with small necrotic volume (median survival difference 14.50 months, p = 0.047) benefited substantially from BVZ.

Conclusion

Imaging biomarkers related to the irregularity of contrast enhancing areas and the necrotic volume were able to discriminate GBM patients with a substantial survival benefit from BVZ. A prospective study is needed to validate our results.

Differences in iron and manganese concentration may confound the measurement of myelin from R1 and R2 relaxation rates in studies of dysmyelination

A model of dysmyelination, the Long Evans Shaker (les) rat, was used to study the contribution of myelin to MR tissue properties in white matter. A large region of white matter was identified in the deep cerebellum and was used for measurements of the MR relaxation rate constants, R1  = 1/T1 and R2  = 1/T2 , at 7 T. In this study, R1 of the les deep cerebellar white matter was found to be 0.55 ± 0.08 s (-1) and R2 was found to be 15 ± 1 s(-1) , revealing significantly lower R1 and R2 in les white matter relative to wild-type (wt: R1  = 0.69 ± 0.05 s(-1) and R2  = 18 ± 1 s(-1) ). These deviated from the expected ΔR1 and ΔR2 values, given a complete lack of myelin in the les white matter, derived from the literature using values of myelin relaxivity, and we suspect that metals could play a significant role. The absolute concentrations of the paramagnetic transition metals iron (Fe) and manganese (Mn) were measured by a micro-synchrotron radiation X-ray fluorescence (μSRXRF) technique, with significantly greater Fe and Mn in les white matter than in wt (in units of μg [metal]/g [wet weight tissue]: les: Fe concentration,19 ± 1; Mn concentration, 0.71 ± 0.04; wt: Fe concentration,10 ± 1; Mn concentration, 0.47 ± 0.04). These changes in Fe and Mn could explain the deviations in R1 and R2 from the expected values in white matter. Although it was found that the influence of myelin still dominates R1 and R2 in wt rats, there were non-negligible changes in the contribution of the metals to relaxation. Although there are already problems with the estimation of myelin from R1 and R2 changes in disease models with pathology that also affects the relaxation rate constants, this study points to a specific pitfall in the estimation of changes in myelin in diseases or models with disrupted concentrations of paramagnetic transition metals. Copyright © 2016 John Wiley & Sons, Ltd.

Detection of Early Tumor Response to Axitinib in Advanced Hepatocellular Carcinoma by Dynamic Contrast Enhanced Ultrasound

This study aimed to evaluate the utility of dynamic contrast-enhanced ultrasound (DCE-US) in measuring early tumor response of advanced hepatocellular carcinoma to axitinib. Twenty patients were enrolled (aged 18-78 y; median 65). DCE-US was performed with bolus injection and infusion/disruption replenishment. Median overall survival was 7.1 mo (1.8-27.3) and progression free survival was 3.6 mo (1.8-17.4). Fifteen patients completed infusion scans and 12 completed bolus scans at 2 wk. Among the perfusion parameters, fractional blood volume at infusion (INFBV) decreased at 2 wk in 10/15 (16%-81% of baseline, mean 47%) and increased in 5/15 (116%-535%, mean 220%). This was not significantly associated with progression free survival (p = 0.310) or progression at 16 wk (p = 0.849), but was borderline statistically significant (p = 0.050) with overall survival, limited by a small sample size. DCE-US is potentially useful in measuring early tumor response of advanced hepatocellular carcinoma to axitinib, but a larger trial is needed.

Quantitative Computed Tomography Imaging Biomarkers in the Diagnosis and Management of Lung Cancer

Tumor diameter has traditionally been used as a standard metric in terms of diagnosis and prognosis prediction of lung cancer. However, recent advances in imaging techniques and data analyses have enabled novel quantitative imaging biomarkers that can characterize disease status more comprehensively and/or predict tumor behavior more precisely. The most widely used imaging modality for lung tumor assessment is computed tomography. Therefore, we focused on computed tomography imaging biomarkers such as tumor volume and mass, ground-glass opacities, perfusion parameters, as well as texture features in this review. Herein, we first appraised the conventional 1- or 2-dimensional measurement with brief discussion on their limits and then introduced the potential imaging biomarkers with emphasis on the current understanding of their clinical usefulness with respect to the malignancy differentiation, treatment response monitoring, and patient outcome prediction.