Heterogeneity in DCE-MRI parametric maps: a biomarker for treatment response?

Habitat escalated adaptive therapy (HEAT): a phase 2 trial utilizing radiomic habitat-directed and genomic-adjusted radiation dose (GARD) optimization for high-grade soft tissue sarcoma

BACKGROUND

Soft tissue sarcomas (STS), have significant inter- and intra-tumoral heterogeneity, with poor response to standard neoadjuvant radiotherapy (RT). Achieving a favorable pathologic response (FPR ≥ 95%) from RT is associated with improved patient outcome. Genomic adjusted radiation dose (GARD), a radiation-specific metric that quantifies the expected RT treatment effect as a function of tumor dose and genomics, proposed that STS is significantly underdosed. STS have significant radiomic heterogeneity, where radiomic habitats can delineate regions of intra-tumoral hypoxia and radioresistance. We designed a novel clinical trial, Habitat Escalated Adaptive Therapy (HEAT), utilizing radiomic habitats to identify areas of radioresistance within the tumor and targeting them with GARD-optimized doses, to improve FPR in high-grade STS.

METHODS

Phase 2 non-randomized single-arm clinical trial includes non-metastatic, resectable high-grade STS patients. Pre-treatment multiparametric MRIs (mpMRI) delineate three distinct intra-tumoral habitats based on apparent diffusion coefficient (ADC) and dynamic contrast enhanced (DCE) sequences. GARD estimates that simultaneous integrated boost (SIB) doses of 70 and 60 Gy in 25 fractions to the highest and intermediate radioresistant habitats, while the remaining volume receives standard 50 Gy, would lead to a > 3 fold FPR increase to 24%. Pre-treatment CT guided biopsies of each habitat along with clip placement will be performed for pathologic evaluation, future genomic studies, and response assessment. An mpMRI taken between weeks two and three of treatment will be used for biological plan adaptation to account for tumor response, in addition to an mpMRI after the completion of radiotherapy in addition to pathologic response, toxicity, radiomic response, disease control, and survival will be evaluated as secondary endpoints. Furthermore, liquid biopsy will be performed with mpMRI for future ancillary studies.

DISCUSSION

This is the first clinical trial to test a novel genomic-based RT dose optimization (GARD) and to utilize radiomic habitats to identify and target radioresistance regions, as a strategy to improve the outcome of RT-treated STS patients. Its success could usher in a new phase in radiation oncology, integrating genomic and radiomic insights into clinical practice and trial designs, and may reveal new radiomic and genomic biomarkers, refining personalized treatment strategies for STS.

TRIAL REGISTRATION

NCT05301283.

TRIAL STATUS

The trial started recruitment on March 17, 2022.

Perfusion Imaging of the Musculoskeletal System

Perfusion imaging is the aspect of functional imaging, which is most applicable to the musculoskeletal system. In this review, the anatomy and physiology of bone perfusion is briefly outlined as are the methods of acquiring perfusion data on MR imaging. The current clinical indications of perfusion related to the assessment of soft tissue and bone tumors, synovitis, osteoarthritis, avascular necrosis, Keinbock's disease, diabetic foot, osteochondritis dissecans, and Paget's disease of bone are reviewed. Challenges and opportunities related to perfusion imaging of the musculoskeletal system are also briefly addressed.

Mapping Spatiotemporal Heterogeneity in Tumor Profiles by Integrating High-Throughput Imaging and Omics Analysis

Intratumoral heterogeneity associates with more aggressive disease progression and worse patient outcomes. Understanding the reasons enabling the emergence of such heterogeneity remains incomplete, which restricts our ability to manage it from a therapeutic perspective. Technological advancements such as high-throughput molecular imaging, single-cell omics, and spatial transcriptomics allow recording of patterns of spatiotemporal heterogeneity in a longitudinal manner, thus offering insights into the multiscale dynamics of its evolution. Here, we review the latest technological trends and biological insights from molecular diagnostics as well as spatial transcriptomics, both of which have witnessed burgeoning growth in the recent past in terms of mapping heterogeneity within tumor cell types as well as the stromal constitution. We also discuss ongoing challenges, indicating possible ways to integrate insights across these methods to have a systems-level spatiotemporal map of heterogeneity in each tumor and a more systematic investigation of the implications of heterogeneity for patient outcomes.

Fractal Dimension Analysis of Pixel Dynamic Contrast Enhanced-Magnetic Resonance Imaging Pharmacokinetic Parameters for Discrimination of Benign and Malignant Breast Lesions

PURPOSE

Breast cancer is the most frequent cancer in women worldwide. However, its diagnosis mostly depends on visual examination of radiologic images, leading to an overdiagnosis with substantial costs. Therefore, a quantitative approach such as dynamic contrast enhanced (DCE)-magnetic resonance imaging (MRI) through pharmacokinetic (PK) modeling is required for reliable analysis. As PK parameters lack information on parameter heterogeneity, texture-based analysis is required to quantify PK parameter heterogeneity. Therefore, this study focused on determining the usefulness of fractal dimension (FD) as a potential imaging biomarker of tumor heterogeneity for discriminating benign and malignant breast lesions.

METHODS

Parametric maps for PK parameters, extravasation rate of contrast agent from blood plasma to extravascular extracellular space (K^trans) and volume fraction of extravascular extracellular space (v_e), were generated for the regions of interest (ROIs) under the standard model using 18 lesions. Then, tumor ROI and pixel DCE-MRI time-course data were analyzed to extract pixel values of K^trans and v_e. For each ROI, FD values of K^trans and v_e were computed using the blanket method.

RESULTS

The FD values of K^trans for benign and malignant lesions varied from 2.96 to 3.49 and from 2.37 to 3.16, respectively, whereas FD values of v_e for benign and malignant lesions varied from 3.01 to 5.15 and 2.42 to 3.44, respectively. There were significant differences in FD values derived from K^trans parametric maps (P = .0053) and v_e parametric maps (P = .0271) between benign and malignant lesions according to the statistical analysis.

CONCLUSION

Incorporating texture heterogeneity changes in breast lesions captured by FD with quantitative DCE-MRI parameters generated under the standard model is a potential marker for prediction of malignant lesions.

Assessment of Early Response to Lung Cancer Chemotherapy by Semiquantitative Analysis of Dynamic Contrast-Enhanced MRI

Objective

To evaluate the early chemotherapy response in patients with lung cancer using semiquantitative analysis of dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI).

Methods

Twenty-two patients with lung cancer treated with chemotherapy were subjected to DCE-MRI at two time points: before starting treatment and after one week of therapy. The image data were collected by DCE-MRI, and the semiquantitative parameters including positive enhancement integral (PEI), signal enhancement ratio (SER), maximum slope of increase (MSI), and time to peak (TTP) were calculated. After chemotherapy, the parameters and relevant variations between the responders and nonresponders were compared with Mann–Whitney U tests. Student's t-test for paired samples was used to evaluate the temporal changes between pre- and posttreatment images.

Results

The patients were categorized as 13 responders and 9 nonresponders based on the tumor response evaluation. After chemotherapy, the PEI, SER, and MSI were significantly increased in responders compared with the pretreatment values (P < 0.05), while no obvious decrease in TTP was observed (P > 0.05). However, 9 nonresponders showed no significant changes in PEI, SER, MSI, and TTP values, as compared with those of pretreatment (P > 0.05). Moreover, the increase of PEI was more dramatically in responders than in nonresponders (P < 0.05), but no significantly differences were observed in SER, MSI, and TTP (P > 0.05).

Conclusion

Semiquantitative analysis of DCE-MRI could provide a reliable noninvasive method for assessing early chemotherapy response in lung cancer patients.

Bridging the macro to micro resolution gap with angiographic optical coherence tomography and dynamic contrast enhanced MRI

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is emerging as a valuable tool for non-invasive volumetric monitoring of the tumor vascular status and its therapeutic response. However, clinical utility of DCE-MRI is challenged by uncertainty in its ability to quantify the tumor microvasculature (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \mathrm{m}$$\end{document}μm scale) given its relatively poor spatial resolution (mm scale at best). To address this challenge, we directly compared DCE-MRI parameter maps with co-registered micron-scale-resolution speckle variance optical coherence tomography (svOCT) microvascular images in a window chamber tumor mouse model. Both semi and fully quantitative (Toft’s model) DCE-MRI metrics were tested for correlation with microvascular svOCT biomarkers. svOCT’s derived vascular volume fraction (VVF) and the mean distance to nearest vessel (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm{DNV} }$$\end{document}DNV¯) metrics were correlated with DCE-MRI vascular biomarkers such as time to peak contrast enhancement (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=-0.81$$\end{document}r=-0.81 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.83$$\end{document}0.83 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both), the area under the gadolinium-time concentration curve (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=0.50$$\end{document}r=0.50 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.48$$\end{document}-0.48 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{trans}$$\end{document}ktrans (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=0.64$$\end{document}r=0.64 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.61$$\end{document}-0.61 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both). Several other correlated micro–macro vascular metric pairs were also noted. The microvascular insights afforded by svOCT may help improve the clinical utility of DCE-MRI for tissue functional status assessment and therapeutic response monitoring applications.

A review on advances in 18F-FDG PET/CT radiomics standardisation and application in lung disease management

Radiomics analysis quantifies the interpolation of multiple and invisible molecular features present in diagnostic and therapeutic images. Implementation of 18-fluorine-fluorodeoxyglucose positron emission tomography/computed tomography (¹⁸F-FDG PET/CT) radiomics captures various disorders in non-invasive and high-throughput manner. ¹⁸F-FDG PET/CT accurately identifies the metabolic and anatomical changes during cancer progression. Therefore, the application of ¹⁸F-FDG PET/CT in the field of oncology is well established. Clinical application of ¹⁸F-FDG PET/CT radiomics in lung infection and inflammation is also an emerging field. Combination of bioinformatics approaches or textual analysis allows radiomics to extract additional information to predict cell biology at the micro-level. However, radiomics texture analysis is affected by several factors associated with image acquisition and processing. At present, researchers are working on mitigating these interrupters and developing standardised workflow for texture biomarker establishment. This review article focuses on the application of ¹⁸F-FDG PET/CT in detecting lung diseases specifically on cancer, infection and inflammation. An overview of different approaches and challenges encountered on standardisation of ¹⁸F-FDG PET/CT technique has also been highlighted. The review article provides insights about radiomics standardisation and application of ¹⁸F-FDG PET/CT in lung disease management.

Quantification of tumor heterogeneity: from data acquisition to metric generation

Tumors are unique and complex ecosystems, in which heterogeneous cell subpopulations with variable molecular profiles, aggressiveness, and proliferation potential coexist and interact. Understanding how heterogeneity influences tumor progression has important clinical implications for improving diagnosis, prognosis, and treatment response prediction. Several recent innovations in data acquisition methods and computational metrics have enabled the quantification of spatiotemporal heterogeneity across different scales of tumor organization. Here, we summarize the most promising efforts from a common experimental and computational perspective, discussing their advantages, shortcomings, and challenges. With personalized medicine entering a new era of unprecedented opportunities, our vision is that of future workflows integrating across modalities, scales, and dimensions to capture intricate aspects of the tumor ecosystem and to open new avenues for improved patient care.

New method for quantification of intratumoral heterogeneity: a feasibility study on Ktrans maps from preclinical DCE-MRI

OBJECT

To develop new imaging biomarkers of therapeutic efficacy through the quantification of intratumoral microvascular heterogeneity.

MATERIALS AND METHODS

The described method was a combination of non-supervised clustering and extraction of intratumoral complexity features (ICF): number of non-connected objects, volume fraction. It was applied to a set of 3D DCE-MRI K^trans maps acquired previously on tumor bearing mice prior to and on day 4 of anti-angiogenic treatment. Evolutions of ICF were compared to conventional summary statistics (CSS) and to heterogeneity related whole tumor texture features (TF) on treated (n = 9) and control (n = 6) mice.

RESULTS

Computed optimal number of clusters per tumor was 4. Several intratumoral features extracted from the clusters were able to monitor a therapy effect. Whereas no feature significantly changed for the control group, 6 features significantly changed for the treated group (4 ICF, 2 CSS). Among these, 5 also significantly differentiated the two groups (3 ICF, 2 CSS). TF failed in demonstrating differences within and between the two groups.

DISCUSSION

ICF are potential imaging biomarkers for anti-angiogenic therapy assessment. The presented method may be expected to have advantages with respect to texture analysis-based methods regarding interpretability of results and setup of standardized image analysis protocols.

Advanced Imaging and Computational Techniques for the Diagnostic and Prognostic Assessment of Malignant Gliomas

ABSTRACT

Advanced imaging techniques provide a powerful tool to assess the intratumoral and intertumoral heterogeneity of gliomas. Advances in the molecular understanding of glioma subgroups may allow improved diagnostic assessment combining imaging and molecular tumor features, with enhanced prognostic utility and implications for patient treatment. In this article, a comprehensive overview of the physiologic basis for conventional and advanced imaging techniques is presented, and clinical applications before and after treatment are discussed. An introduction to the principles of radiomics and the advanced integration of imaging, clinical outcomes, and genomic data highlights the future potential for this field of research to better stratify and select patients for standard as well as investigational therapies.

Primary Gastro-Intestinal Lymphoma and Gastro-Intestinal Adenocarcinoma: An Initial Study of CT Texture Analysis as Quantitative Biomarkers for Differentiation

BACKGROUND

To explore the potential role of computed tomography (CT) texture analysis and an imaging biomarker in differentiating primary gastro-intestinal lymphoma (PGIL) from gastro-intestinal adenocarcinoma (GIAC).

METHODS

A total of 131 patients with surgical pathologically PGIL and GIAC were enrolled in this study. Histogram parameters of arterial and venous phases extracted from contrast enhanced modified discrete cosine transform (MDCT) images were compared between PGIL and GIAC by Mann-Whitney U tests. The optimal parameters for differentiating these two groups were obtained through receiver operating characteristic (ROC) curves and the area under the curve (AUC) was calculated.

RESULTS

Compared with GIAC, in arterial phase, PGIL had statistically higher 5th, 10th percentiles (p = 0.003 and 0.011) and statistically lower entropy (p = 0.001). In the venous phase, PGIL had statistically lower mean, median, 75th, 90th, 95th percentiles, and entropy (p = 0.036, 0.029, 0.007, 0.001 and 0.001, respectively). For differentiating PGIL from GIAC, V-median + A-5th percentile was an optimal parameter for combined diagnosis (AUC = 0.746, p < 0.0001), and the corresponding sensitivity and specificity were 81.7 and 64.8%, respectively.

CONCLUSION

CT texture analysis could be useful for differential diagnosis of PGIL and GIAC.

The Biological Meaning of Radiomic Features

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

Personalising sarcoma care using quantitative multimodality imaging for response assessment

Over the last decades, technological developments in the field of radiology have resulted in a widespread use of imaging for personalising medicine in oncology, including patients with a sarcoma. New scanner hardware, imaging protocols, image reconstruction algorithms, radiotracers, and contrast media, enabled the assessment of the physical and biological properties of tumours associated with response to treatment. In this context, medical imaging has the potential to select sarcoma patients who do not benefit from (neo-)adjuvant treatment and facilitate treatment adaptation. Due to the biological heterogeneity in sarcomas, the challenge at hand is to acquire a practicable set of imaging features for specific sarcoma subtypes, allowing response assessment. This review provides a comprehensive overview of available clinical data on imaging-based response monitoring in sarcoma patients and future research directions. Eventually, it is expected that imaging-based response monitoring will help to achieve successful modification of (neo)adjuvant treatments and improve clinical care for these patients.

Quantitative tumor heterogeneity MRI profiling improves machine learning-based prognostication in patients with metastatic colon cancer

OBJECTIVES

Intra-tumor heterogeneity has been previously shown to be an independent predictor of patient survival. The goal of this study is to assess the role of quantitative MRI-based measures of intra-tumor heterogeneity as predictors of survival in patients with metastatic colorectal cancer.

METHODS

In this IRB-approved retrospective study, we identified 55 patients with stage 4 colon cancer with known hepatic metastasis on MRI. Ninety-four metastatic hepatic lesions were identified on post-contrast images and manually volumetrically segmented. A heterogeneity phenotype vector was extracted from each lesion. Univariate regression analysis was used to assess the contribution of 110 extracted features to survival prediction. A random forest-based machine learning technique was applied to the feature vector and to the standard prognostic clinical and pathologic variables. The dataset was divided into a training and test set at a ratio of 4:1. ROC analysis and confusion matrix analysis were used to assess classification performance.

RESULTS

Mean survival time was 39 ± 3.9 months for the study population. A total of 22 texture features were associated with patient survival (p < 0.05). The trained random forest machine learning model that included standard clinical and pathological prognostic variables resulted in an area under the ROC curve of 0.83. A model that adds imaging-based heterogeneity features to the clinical and pathological variables resulted in improved model performance for survival prediction with an AUC of 0.94.

CONCLUSIONS

MRI-based texture features are associated with patient outcomes and improve the performance of standard clinical and pathological variables for predicting patient survival in metastatic colorectal cancer.

KEY POINTS

• MRI-based tumor heterogeneity texture features are associated with patient survival outcomes. • MRI-based tumor texture features complement standard clinical and pathological variables for prognosis prediction in metastatic colorectal cancer. • Agglomerative hierarchical clustering shows that patient survival outcomes are associated with different MRI tumor profiles.

Evaluating Heterogeneity of Primary Lung Tumor Using Clinical Routine Magnetic Resonance Imaging and a Tumor Heterogeneity Index

Objective

To improve the assessment of primary tumor heterogeneity in magnetic resonance imaging (MRI) of non-small cell lung cancer (NSCLC), we proposed a method using basic measurements from T1- and T2-weighted MRI.

Methods

One hundred and four NSCLC patients with different T stages were studied. Fifty-two patients were analyzed as training group and another 52 as testing group. The ratios of standard deviation (SD)/mean signal value of primary tumor from T1-weighted (T1WI), T1-enhanced (T1C), T2-weighted (T2WI), and T2 fat suppression (T2fs) images were calculated. In the training group, correlation analyses were performed between the ratios and T stages. Then an ordinal regression model was built to generate the tumor heterogeneous index (THI) for evaluating the heterogeneity of tumor. The model was validated in the testing group.

Results

There were 11, 32, 40, and 21 patients with T1, T2, T3, and T4 disease, respectively. In the training group, the median SD/mean on T1WI, T1C, T2WI, and T2fs sequences was 0.11, 0.19, 0.16, and 0.15 respectively. The SD/mean on T1C (p=0.003), T2WI (p=0.000), and T2fs sequences (p=0.002) correlated significantly with T stages. Patients with more advanced T stage showed higher SD/mean on T2-weighted, T2fs, and T1C sequences. The median THI in the training group was 2.15. THI correlated with T stage significantly (p=0.000). In the testing group, THI was also significantly related to T stages (p=0.001). Higher THI had relevance to more advanced T stage.

Conclusions

The proposed ratio measurements and THI based on MRI can serve as functional radiomic markers that correlated with T stages for evaluating heterogeneity of lung tumors.

A Role for FDG PET Radiomics in Personalized Medicine?

Radiomics describes the extraction of multiple features from medical images, including molecular imaging modalities, that with bioinformatic approaches, provide additional clinically relevant information that may be invisible to the human eye. This information may complement standard radiological interpretation with data that may better characterize a disease or that may provide predictive or prognostic information. Progressing from predefined image features, often describing heterogeneity of voxel intensities within a volume of interest, there is increasing use of machine learning to classify disease characteristics and deep learning methods based on artificial neural networks that can learn features without a priori definition and without the need for preprocessing of images. There have been advances in standardization and harmonization of methods to a level that should support multicenter studies. However, in this relatively early phase of research in the field, there are limited aspects that have been adopted into routine practice. Most of the reports in the molecular imaging field describe radiomic approaches in cancer using 18F-fluorodeoxyglucose positron emission tomography (18F-FDG-PET). In this review, we will describe radiomics in molecular imaging and summarize the pertinent literature in lung cancer where reports are most prevalent and mature.

MR Imaging Texture Analysis in the Abdomen and Pelvis

Add "which is a" before "distribution"? Texture analysis (TA) is a form of radiomics that refers to quantitative measurements of the histogram, distribution and/or relationship of pixel intensities or gray scales within a region of interest on an image. TA can be applied to MR images of the abdomen and pelvis, with the main strength quantitative analysis of pixel intensities and heterogeneity rather than subjective/qualitative analysis. There are multiple limitations of MRTA. Despite these limitations, there is a growing body of literature supporting MRTA. This review discusses application of MRTA to the abdomen and pelvis.

Selection of Fitting Model and Arterial Input Function for Repeatability in Dynamic Contrast-Enhanced Prostate MRI

RATIONALE AND OBJECTIVES

Analysis of dynamic contrast-enhanced (DCE) magnetic resonance imaging is notable for the variability of calculated parameters. The purpose of this study was to evaluate the level of measurement variability and error/variability due to modeling in DCE magnetic resonance imaging parameters.

MATERIALS AND METHODS

Two prostate DCE scans were performed on 11 treatment-naïve patients with suspected or confirmed prostate peripheral zone cancer within an interval of less than two weeks. Tumor-suspicious and normal-appearing regions of interest (ROI) in the prostate peripheral zone were segmented. Different Tofts-Kety based models and different arterial input functions, with and without bolus arrival time (BAT) correction, were used to extract pharmacokinetic parameters. The percent repeatability coefficient (%RC) of fitted model parameters K_trans, v_e, and k_ep was calculated. Paired t-tests comparing parameters in tumor-suspicious ROIs and in normal-appearing tissue evaluated each parameter's sensitivity to pathology.

RESULTS

Although goodness-of-fit criteria favored the four-parameter extended Tofts-Kety model with the BAT correction included, the simplest two-parameter Tofts-Kety model overall yielded the best repeatability scores. The best %RC in the tumor-suspicious ROI was 63% for k_ep, 28% for v_e, and 83% for K_trans . The best p values for discrimination between tissues were p <10^-5 for k_ep and K_trans, and p = 0.11 for v_e. Addition of the BAT correction to the models did not improve repeatability.

CONCLUSION

The parameter k_ep, using an arterial input functions directly measured from blood signals, was more repeatable than K_trans. Both K_trans and k_ep values were highly discriminatory between healthy and diseased tissues in all cases. The parameter v_e had high repeatability but could not distinguish the two tissue types.

Pharmacokinetic modeling of dynamic contrast-enhanced MRI using a reference region and input function tail

PURPOSE

Quantitative analysis of dynamic contrast-enhanced MRI (DCE-MRI) requires an arterial input function (AIF) which is difficult to measure. We propose the reference region and input function tail (RRIFT) approach which uses a reference tissue and the washout portion of the AIF.

METHODS

RRIFT was evaluated in simulations with 100 parameter combinations at various temporal resolutions (5-30 s) and noise levels (σ = 0.01-0.05 mM). RRIFT was compared against the extended Tofts model (ETM) in 8 studies from patients with glioblastoma multiforme. Two versions of RRIFT were evaluated: one using measured patient-specific AIF tails, and another assuming a literature-based AIF tail.

RESULTS

RRIFT estimated the transfer constant K trans and interstitial volume v e with median errors within 20% across all simulations. RRIFT was more accurate and precise than the ETM at temporal resolutions slower than 10 s. The percentage error of K trans had a median and interquartile range of -9 ± 45% with the ETM and -2 ± 17% with RRIFT at a temporal resolution of 30 s under noiseless conditions. RRIFT was in excellent agreement with the ETM in vivo, with concordance correlation coefficients (CCC) of 0.95 for K trans , 0.96 for v e , and 0.73 for the plasma volume v p using a measured AIF tail. With the literature-based AIF tail, the CCC was 0.89 for K trans , 0.93 for v e and 0.78 for v p .

CONCLUSIONS

Quantitative DCE-MRI analysis using the input function tail and a reference tissue yields absolute kinetic parameters with the RRIFT method. This approach was viable in simulation and in vivo for temporal resolutions as low as 30 s.

Molecular Imaging of Deoxycytidine Kinase Activity Using Deoxycytidine-Enhanced CEST MRI

Deoxycytidine kinase (DCK) is a key enzyme for the activation of a broad spectrum of nucleoside-based chemotherapy drugs (e.g., gemcitabine); low DCK activity is one of the most important causes of cancer drug-resistance. Noninvasive imaging methods that can quantify DCK activity are invaluable for assessing tumor resistance and predicting treatment efficacy. Here we developed a "natural" MRI approach to detect DCK activity using its natural substrate deoxycytidine (dC) as the imaging probe, which can be detected directly by chemical exchange saturation transfer (CEST) MRI without any synthetic labeling. CEST MRI contrast of dC and its phosphorylated form, dCTP, successfully discriminated DCK activity in two mouse leukemia cell lines with different DCK expression. This dC-enhanced CEST MRI in xenograft leukemic cancer mouse models demonstrated that DCK(+) tumors have a distinctive dynamic CEST contrast enhancement and a significantly higher CEST contrast than DCK(-) tumors (AUC^{0-60 min} = 0.47 ± 0.25 and 0.20 ± 0.13, respectively; P = 0.026, paired Student t test, n = 4) at 1 hour after the injection of dC. dC-enhanced CEST contrast also correlated well with tumor responses to gemcitabine treatment. This study demonstrates a novel MR molecular imaging approach for predicting cancer resistance using natural, nonradioactive, nonmetallic, and clinically available agents. This method has great potential for pursuing personalized chemotherapy by stratifying patients with different DCK activity. SIGNIFICANCE: A new molecular MRI method that detects deoxycytidine kinase activity using its natural substrate deoxycytidine has great translational potential for clinical assessment of tumor resistance and prediction of treatment efficacy.

Biofabrication strategies for creating microvascular complexity

Design and fabrication of effective biomimetic vasculatures constitutes a relevant and yet unsolved challenge, lying at the heart of tissue repair and regeneration strategies. Even if cell growth is achieved in 3D tissue scaffolds or advanced implants, tissue viability inevitably requires vascularization, as diffusion can only transport nutrients and eliminate debris within a few hundred microns. This engineered vasculature may need to mimic the intricate branching geometry of native microvasculature, referred to herein as vascular complexity, to efficiently deliver blood and recreate critical interactions between the vascular and perivascular cells as well as parenchymal tissues. This review first describes the importance of vascular complexity in labs- and organs-on-chips, the biomechanical and biochemical signals needed to create and maintain a complex vasculature, and the limitations of current 2D, 2.5D, and 3D culture systems in recreating vascular complexity. We then critically review available strategies for design and biofabrication of complex vasculatures in cell culture platforms, labs- and organs-on-chips, and tissue engineering scaffolds, highlighting their advantages and disadvantages. Finally, challenges and future directions are outlined with the hope of inspiring researchers to create the reliable, efficient and sustainable tools needed for design and biofabrication of complex vasculatures.

The use of texture-based radiomics CT analysis to predict outcomes in early-stage non-small cell lung cancer treated with stereotactic ablative radiotherapy

OBJECTIVE:

Stereotactic ablative radiotherapy (SABR) is being increasingly used as a non-invasive treatment for early-stage non-small cell lung cancer (NSCLC). A non-invasive method to estimate treatment outcomes in these patients would be valuable, especially since access to tissue specimens is often difficult in these cases.

METHODS:

We developed a method to predict survival following SABR in NSCLC patients using analysis of quantitative image features on pre-treatment CT images. We developed a Cox Lasso model based on two-dimensional Riesz wavelet quantitative texture features on CT scans with the goal of separating patients based on survival.

RESULTS:

The median log-rank p-value for 1000 cross-validations was 0.030. Our model was able to separate patients based upon predicted survival. When we added tumor size into the model, the p-value lost its significance, demonstrating that tumor size is not a key feature in the model but rather decreases significance likely due to the relatively small number of events in the dataset. Furthermore, running the model using Riesz features extracted either from the solid component of the tumor or from the ground glass opacity (GGO) component of the tumor maintained statistical significance. However, the p-value improved when combining features from the solid and the GGO components, demonstrating that there are important data that can be extracted from the entire tumor.

CONCLUSIONS:

The model predicting patient survival following SABR in NSCLC may be useful in future studies by enabling prediction of survival-based outcomes using radiomics features in CT images.

ADVANCES IN KNOWLEDGE:

Quantitative image features from NSCLC nodules on CT images have been found to significantly separate patient populations based on overall survival (p = 0.04). In the long term, a non-invasive method to estimate treatment outcomes in patients undergoing SABR would be valuable, especially since access to tissue specimens is often difficult in these cases.

Challenges and Promises of PET Radiomics

PURPOSE

Radiomics describes the extraction of multiple, otherwise invisible, features from medical images that, with bioinformatic approaches, can be used to provide additional information that can predict underlying tumor biology and behavior.

METHODS AND MATERIALS

Radiomic signatures can be used alone or with other patient-specific data to improve tumor phenotyping, treatment response prediction, and prognosis, noninvasively. The data describing 18F-fluorodeoxyglucose positron emission tomography radiomics, often using texture or heterogeneity parameters, are increasing rapidly.

RESULTS

In relation to radiation therapy practice, early data have reported the use of radiomic approaches to better define tumor volumes and predict radiation toxicity and treatment response.

CONCLUSIONS

Although at an early stage of development, with many technical challenges remaining and a need for standardization, promise nevertheless exists that PET radiomics will contribute to personalized medicine, especially with the availability of increased computing power and the development of machine-learning approaches for imaging.

Investigating the role of functional imaging in the management of soft-tissue sarcomas of the extremities

Background and purpose

In this work, we validate a texture-based model computed from positron emission tomography (PET) and magnetic resonance imaging (MRI) for the prediction of lung metastases in soft-tissue sarcomas (STS). We explore functional imaging at different treatment time points and evaluate the feasibility of radiotherapy dose painting as a potential treatment strategy for patients with higher metastatic risk.

Materials and methods

We acquired fluorodeoxyglucose (FDG)-PET, fluoromisonidazole (FMISO)-PET, diffusion weighting (DW)-MRI and dynamic contrast enhanced (DCE)-MRI data for 18 patients with extremity STS before, during, and after pre-operative radiotherapy. We tested the lung metastases prediction model using pre-treatment images. We evaluated the feasibility of dose painting using volumetric arc therapy (VMAT) via treatment re-planning with a prescription of 50 Gy to the planning target volume (PTV_50Gy) and boost doses of 60 Gy to the FDG hypermetabolic gross tumour volume (GTV_60Gy) and 65 Gy to the low-perfusion DCE-MRI hypoxic GTV contained within the GTV_60Gy (GTV_65Gy).

Results

The texture-based model for lung metastases prediction reached an area under the curve (AUC), sensitivity, specificity and accuracy of 0.71, 0.75, 0.85 and 0.82, respectively. Dose painting resulted in adequate coverage and homogeneity in the re-planned treatments: D_95% to the PTV_50Gy, GTV_60Gy and GTV_65Gy were 50.0 Gy, 60.3 Gy and 65.4 Gy, respectively.

Conclusions

Textural biomarkers extracted from FDG-PET and MRI could be useful to identify STS patients that might benefit from dose escalation. The feasibility of treatment planning with double boost levels to intratumoural GTV functional sub-volumes was established.

Exploring Applications of Radiomics in Magnetic Resonance Imaging of Head and Neck Cancer: A Systematic Review

Background

Radiomics has been widely investigated for non-invasive acquisition of quantitative textural information from anatomic structures. While the vast majority of radiomic analysis is performed on images obtained from computed tomography, magnetic resonance imaging (MRI)-based radiomics has generated increased attention. In head and neck cancer (HNC), however, attempts to perform consistent investigations are sparse, and it is unclear whether the resulting textural features can be reproduced. To address this unmet need, we systematically reviewed the quality of existing MRI radiomics research in HNC.

Methods

Literature search was conducted in accordance with guidelines established by Preferred Reporting Items for Systematic Reviews and Meta-Analyses. Electronic databases were examined from January 1990 through November 2017 for common radiomic keywords. Eligible completed studies were then scored using a standardized checklist that we developed from Enhancing the Quality and Transparency of Health Research guidelines for reporting machine-learning predictive model specifications and results in biomedical research, defined by Luo et al. (1). Descriptive statistics of checklist scores were populated, and a subgroup analysis of methodology items alone was conducted in comparison to overall scores.

Results

Sixteen completed studies and four ongoing trials were selected for inclusion. Of the completed studies, the nasopharynx was the most common site of study (37.5%). MRI modalities varied with only four of the completed studies (25%) extracting radiomic features from a single sequence. Study sample sizes ranged between 13 and 118 patients (median of 40), and final radiomic signatures ranged from 2 to 279 features. Analyzed endpoints included either segmentation or histopathological classification parameters (44%) or prognostic and predictive biomarkers (56%). Liu et al. (2) addressed the highest number of our checklist items (total score: 48), and a subgroup analysis of methodology checklist items alone did not demonstrate any difference in scoring trends between studies [Spearman’s ρ = 0.94 (p < 0.0001)].

Conclusion

Although MRI radiomic applications demonstrate predictive potential in analyzing diverse HNC outcomes, methodological variances preclude accurate and collective interpretation of data.

Quantitative image analysis for evaluation of tumor response in clinical oncology

The objective, accurate, and standardized evaluation of tumor response to treatment is an indispensable procedure in clinical oncology. Compared to manual measurement, computer-assisted linear measurement can significantly improve the accuracy and reproducibility of tumor burden quantification. For irregular-shaped and infiltrating or diffuse tumors, which are difficult to quantify by linear measurement, computer-assisted volumetric measurement may provide a more objective and sensitive quantification to evaluate tumor response to treatment than linear measurement does. In the evaluation of tumor response to novel oncologic treatments such as targeted therapy, changes in overall tumor size do not necessarily reflect tumor response to therapy due to the presence of internal necrosis or hemorrhages. This leads to a new generation of imaging biomarkers to evaluate tumor response by using texture analysis methods, also called radiomics. Computer-assisted texture analysis technology offers a more comprehensive and in-depth imaging biomarker to evaluate tumor response. The application of computer-assisted quantitative imaging analysis techniques not only reduces the inaccuracy and improves the reliability in tumor burden quantification, but facilitates the development of more comprehensive and intelligent approaches to evaluate treatment response, and hence promotes precision imaging in the evaluation of tumor response in clinical oncology. This article summarizes the state-of-the-art technical developments and clinical applications of quantitative imaging analysis in evaluation of tumor response in clinical oncology.

Characterization of tumor vascular permeability using natural dextrans and CEST MRI

PURPOSE

To investigate the use of natural dextrans as nano-sized chemical exchange saturation transfer (CEST) MRI probes for characterizing size-dependent tumor vascular permeability.

METHODS

Dextrans of different molecular weight (10, 70, 150, and 2000 kD) were characterized for their CEST contrast. Mice (N = 5) bearing CT26 subcutaneous colon tumors were injected intravenously with 10 kD (D10, 6 nm) and 70 kD (D70, 12 nm) dextran at a dose of 375 mg/kg. The CEST-MRI signal in the tumors was assessed before and approximately 40 min after each injection using a dynamic CEST imaging scheme.

RESULTS

All dextrans of different molecular weights have a strong CEST signal with an apparent maximum of approximately 0.9 ppm. The detectability and effects of pH and saturation conditions (B₁ and T_sat ) were investigated. When applied to CT26 tumors, the injection of D10 could produce a significant "dexCEST" enhancement in the majority of the tumor area, whereas the injection of D70 only resulted in an increase in the tumor periphery. Quantitative analysis revealed the differential permeability of CT26 tumors to different size particles, which was validated by fluorescence imaging and immunohistochemistry.

CONCLUSIONS

As a first application, we used 10- and 70-kD dextrans to visualize the spatially variable, size-dependent permeability in the tumor, indicating that nano-sized dextrans can be used for characterizing tumor vascular permeability with dexCEST MRI and, potentially, for developing dextran-based theranostic drug delivery systems. Magn Reson Med 79:1001-1009, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Primary Rectal Cancer: Repeatability of Global and Local-Regional MR Imaging Texture Features

Purpose To assess the day-to-day repeatability of global and local-regional magnetic resonance (MR) imaging texture features derived from primary rectal cancer. Materials and Methods After ethical approval and patient informed consent were obtained, two pretreatment T2-weighted axial MR imaging studies performed prospectively with the same imaging unit on 2 consecutive days in 14 patients with rectal cancer (11 men [mean age, 61.7 years], three women [mean age, 70.0 years]) were analyzed to extract (a) global first-order statistical histogram and model-based fractal features reflecting the whole-tumor voxel intensity histogram distribution and repeating patterns, respectively, without spatial information and (b) local-regional second-order and high-order statistical texture features reflecting the intensity and spatial interrelationships between adjacent in-plane or multiplanar voxels or regions, respectively. Repeatability was assessed for 46 texture features, and mean difference, 95% limits of agreement, within-subject coefficient of variation (wCV), and repeatability coefficient (r) were recorded. Results Repeatability was better for global parameters than for most local-regional parameters. In particular, histogram mean, median, and entropy, fractal dimension mean and standard deviation, and second-order entropy, homogeneity, difference entropy, and inverse difference moment demonstrated good repeatability, with narrow limits of agreement and wCVs of 10% or lower. Repeatability was poorest for the following high-order gray-level run-length (GLRL) gray-level zone size matrix (GLZSM) and neighborhood gray-tone difference matrix (NGTDM) parameters: GLRL intensity variability, GLZSM short-zone emphasis, GLZSM intensity nonuniformity, GLZSM intensity variability, GLZSM size zone variability, and NGTDM complexity, demonstrating wider agreement limits and wCVs of 50% or greater. Conclusion MR imaging repeatability is better for global texture parameters than for local-regional texture parameters, indicating that global texture parameters should be sufficiently robust for clinical practice. Online supplemental material is available for this article.

Application of Texture Analysis to Study Small Vessel Disease and Blood-Brain Barrier Integrity

Objectives

We evaluate the alternative use of texture analysis for evaluating the role of blood–brain barrier (BBB) in small vessel disease (SVD).

Methods

We used brain magnetic resonance imaging from 204 stroke patients, acquired before and 20 min after intravenous gadolinium administration. We segmented tissues, white matter hyperintensities (WMH) and applied validated visual scores. We measured textural features in all tissues pre- and post-contrast and used ANCOVA to evaluate the effect of SVD indicators on the pre-/post-contrast change, Kruskal–Wallis for significance between patient groups and linear mixed models for pre-/post-contrast variations in cerebrospinal fluid (CSF) with Fazekas scores.

Results

Textural “homogeneity” increase in normal tissues with higher presence of SVD indicators was consistently more overt than in abnormal tissues. Textural “homogeneity” increased with age, basal ganglia perivascular spaces scores (p < 0.01) and SVD scores (p < 0.05) and was significantly higher in hypertensive patients (p < 0.002) and lacunar stroke (p = 0.04). Hypertension (74% patients), WMH load (median = 1.5 ± 1.6% of intracranial volume), and age (mean = 65.6 years, SD = 11.3) predicted the pre/post-contrast change in normal white matter, WMH, and index stroke lesion. CSF signal increased with increasing SVD post-contrast.

Conclusion

A consistent general pattern of increasing textural “homogeneity” with increasing SVD and post-contrast change in CSF with increasing WMH suggest that texture analysis may be useful for the study of BBB integrity.

Haralick texture features from apparent diffusion coefficient (ADC) MRI images depend on imaging and pre-processing parameters

In recent years, texture analysis of medical images has become increasingly popular in studies investigating diagnosis, classification and treatment response assessment of cancerous disease. Despite numerous applications in oncology and medical imaging in general, there is no consensus regarding texture analysis workflow, or reporting of parameter settings crucial for replication of results. The aim of this study was to assess how sensitive Haralick texture features of apparent diffusion coefficient (ADC) MR images are to changes in five parameters related to image acquisition and pre-processing: noise, resolution, how the ADC map is constructed, the choice of quantization method, and the number of gray levels in the quantized image. We found that noise, resolution, choice of quantization method and the number of gray levels in the quantized images had a significant influence on most texture features, and that the effect size varied between different features. Different methods for constructing the ADC maps did not have an impact on any texture feature. Based on our results, we recommend using images with similar resolutions and noise levels, using one quantization method, and the same number of gray levels in all quantized images, to make meaningful comparisons of texture feature results between different subjects.

Diagnostic value of combined diffusion-weighted imaging with dynamic contrast enhancement MRI in differentiating malignant from benign bone lesions

AIM

To determine the diagnostic value of combined diffusion-weighted imaging (DWI) with dynamic contrast enhancement magnetic resonance imaging (DCE-MRI) in differentiating malignant from benign bone lesions.

MATERIALS AND METHODS

DWI and DCE-MRI were performed in 36 patients (14 were benign and 22 were malignant). The mean apparent diffusion coefficient (ADC) values and signal enhanced extent (SEE), slope value, and time-signal intensity curve (TIC) type were recorded by two observers. Between-group comparison was made using the independent sample t-test and receiver-operating characteristic (ROC) analysis.

RESULTS

There was a significant difference between the mean ADC value of the benign ([1.75±0.50]×10^-3 mm²/s) and malignant ([1.11±0.47]×10^-3 mm²/s) groups (p=0.001). The threshold ADC value of ≤1.10×10^-3 mm²/s resulted in a sensitivity of 77.3%, a specificity of 92.9%, and an accuracy of 85.1%. A type III curve was found in 23 cases (21 malignant and two benign), a type II curve was seen in six cases (one malignant and five benign), and a type I curve in seven cases (all were benign). The SEE and slope values in the benign and malignant groups were 227.96±172.08, 325.60±125.86 (p=0.058); 0.97±0.67%/s, 3.19±3.20%/s (p=0.016), respectively. ROC analysis showed a sensitivity of 95.5%, a specificity of 85.7%, and an accuracy of 90.6% for malignancy, based on a slope cut-off value of >1.46%/s. Combining ADC and slope values resulted in a sensitivity of 100%, a specificity of 85.7%, and an accuracy of 92.9%.

CONCLUSIONS

Both DWI and DCE-MRI showed promising results for differentiating malignant from benign bone lesions. A combination of DWI and DCE-MRI was the most valuable of the three.

Dynamic contrast enhanced MRI detects changes in vascular transport rate constants following treatment with thermally-sensitive liposomal doxorubicin

Temperature-sensitive liposomal formulations of chemotherapeutics, such as doxorubicin, can achieve locally high drug concentrations within a tumor and tumor vasculature while maintaining low systemic toxicity. Further, doxorubicin delivery by temperature-sensitive liposomes can reliably cure local cancer in mouse models. Histological sections of treated tumors have detected red blood cell extravasation within tumors treated with temperature-sensitive doxorubicin and ultrasound hyperthermia. We hypothesize that the local release of drug into the tumor vasculature and resulting high drug concentration can alter vascular transport rate constants along with having direct tumoricidal effects. Dynamic contrast enhanced MRI (DCE-MRI) coupled with a pharmacokinetic model can detect and quantify changes in such vascular transport rate constants. Here, we set out to determine whether changes in rate constants resulting from intravascular drug release were detectable by MRI. We found that the accumulation of gadoteridol was enhanced in tumors treated with temperature-sensitive liposomal doxorubicin and ultrasound hyperthermia. While the initial uptake rate of the small molecule tracer was slower (k₁=0.0478±0.011s^-1 versus 0.116±0.047s^-1) in treated compared to untreated tumors, the tracer was retained after treatment due to a larger reduction in the rate of clearance (k₂=0.291±0.030s^-1 versus 0.747±0.24s^-1). While DCE-MRI assesses a combination of blood flow and permeability, ultrasound imaging of microvascular flow rate is sensitive only to changes in vascular flow rate; based on this technique, blood flow was not significantly altered 30min after treatment. In summary, DCE-MRI provides a means to detect changes that are associated with treatment by thermally-activated particles and such changes can be exploited to enhance local delivery.

Texture analysis of medical images for radiotherapy applications

The high-throughput extraction of quantitative information from medical images, known as radiomics, has grown in interest due to the current necessity to quantitatively characterize tumour heterogeneity. In this context, texture analysis, consisting of a variety of mathematical techniques that can describe the grey-level patterns of an image, plays an important role in assessing the spatial organization of different tissues and organs. For these reasons, the potentiality of texture analysis in the context of radiotherapy has been widely investigated in several studies, especially for the prediction of the treatment response of tumour and normal tissues. Nonetheless, many different factors can affect the robustness, reproducibility and reliability of textural features, thus limiting the impact of this technique. In this review, an overview of the most recent works that have applied texture analysis in the context of radiotherapy is presented, with particular focus on the assessment of tumour and tissue response to radiations. Preliminary, the main factors that have an influence on features estimation are discussed, highlighting the need of more standardized image acquisition and reconstruction protocols and more accurate methods for region of interest identification. Despite all these limitations, texture analysis is increasingly demonstrating its ability to improve the characterization of intratumour heterogeneity and the prediction of clinical outcome, although prospective studies and clinical trials are required to draw a more complete picture of the full potential of this technique.

Dynamic fractal signature dissimilarity analysis for therapeutic response assessment using dynamic contrast-enhanced MRI

PURPOSE

To develop a dynamic fractal signature dissimilarity (FSD) method as a novel image texture analysis technique for the quantification of tumor heterogeneity information for better therapeutic response assessment with dynamic contrast-enhanced (DCE)-MRI.

METHODS

A small animal antiangiogenesis drug treatment experiment was used to demonstrate the proposed method. Sixteen LS-174T implanted mice were randomly assigned into treatment and control groups (n = 8/group). All mice received bevacizumab (treatment) or saline (control) three times in two weeks, and one pretreatment and two post-treatment DCE-MRI scans were performed. In the proposed dynamic FSD method, a dynamic FSD curve was generated to characterize the heterogeneity evolution during the contrast agent uptake, and the area under FSD curve (AUCFSD) and the maximum enhancement (MEFSD) were selected as representative parameters. As for comparison, the pharmacokinetic parameter K(trans) map and area under MR intensity enhancement curve AUCMR map were calculated. Besides the tumor's mean value and coefficient of variation, the kurtosis, skewness, and classic Rényi dimensions d1 and d2 of K(trans) and AUCMR maps were evaluated for heterogeneity assessment for comparison. For post-treatment scans, the Mann-Whitney U-test was used to assess the differences of the investigated parameters between treatment/control groups. The support vector machine (SVM) was applied to classify treatment/control groups using the investigated parameters at each post-treatment scan day.

RESULTS

The tumor mean K(trans) and its heterogeneity measurements d1 and d2 values showed significant differences between treatment/control groups in the second post-treatment scan. In contrast, the relative values (in reference to the pretreatment value) of AUCFSD and MEFSD in both post-treatment scans showed significant differences between treatment/control groups. When using AUCFSD and MEFSD as SVM input for treatment/control classification, the achieved accuracies were 93.8% and 93.8% at first and second post-treatment scan days, respectively. In comparison, the classification accuracies using d1 and d2 of K(trans) map were 87.5% and 100% at first and second post-treatment scan days, respectively.

CONCLUSIONS

As quantitative metrics of tumor contrast agent uptake heterogeneity, the selected parameters from the dynamic FSD method accurately captured the therapeutic response in the experiment. The potential application of the proposed method is promising, and its addition to the existing DCE-MRI techniques could improve DCE-MRI performance in early assessment of treatment response.

Evaluation of Head and Neck Tumors with Functional MR Imaging

Head and neck cancer is one of the most common cancers worldwide. MR imaging-based diffusion and perfusion techniques enable the noninvasive assessment of tumor biology and physiology, which supplement information obtained from standard structural scans. Diffusion and perfusion MR imaging techniques provide novel biomarkers that can aid monitoring in pretreatment, during treatment, and posttreatment stages to improve patient selection for therapeutic strategies; provide evidence for change of therapy regime; and evaluate treatment response. This review discusses pertinent aspects of the role of diffusion and perfusion MR imaging and computational analysis methods in studying head and neck cancer.

Simulation of heterogeneous molecular delivery in tumours using μCT reconstructions and MRI validation

PRIMARY OBJECTIVE

Utilizing the detailed vascular network obtained from micro-computed tomography (μCT) to establish a mathematical model of the temporal molecular distribution within a murine C3H mammary carcinoma.

PROCEDURES

Female CDF1 mice with a C3H mammary carcinoma on the right rear foot were used in this study. Dynamic information for each tumour was achieved by Dynamic Contrast Enhanced-Magnetic Resonance Imaging (DCE-MRI) on a 16.4 T system. Detailed morphologic information on the tumour vasculature was obtained by ex vivo μCT and compared to CD34 immunohistochemical staining of tissue sections. The reconstructed vascular network served as origin for the diffusion (described by the apparent diffusion coefficient) within the tumour (the restricted volume described by the interstitial volume fraction derived from DCE-MRI). The resulting partial differential equation was solved using Finite-Element and a combined mathematical graph describing molecular distribution within the tumour was obtained.

RESULTS

The established molecular distribution model predicted a heterogeneous distribution throughout the tumour related to the layout of the vascular network. Central tumour section concentration-time curves estimated from the established molecular distribution model were compared with physical measurements obtained by DCE-MRI of the same tumours and showed excellent correlation.

CONCLUSIONS

A mathematical model describing temporal molecular distribution based on detailed vascular network structures was established and compared to DCE-MRI. The improved morphological insight will enhance future studies of heterogeneous tumours.

Differentiation of tumor vasculature heterogeneity levels in small animals based on total hemoglobin concentration using magnetic resonance-guided diffuse optical tomography in vivo

Insight into the vasculature of the tumor in small animals has the potential to impact many areas of cancer research. The heterogeneity of the vasculature of a tumor is directly related to tumor stage and disease progression. In this small scale animal study, we investigated the feasibility of differentiating tumors with different levels of vasculature heterogeneity in vivo using a previously developed hybrid magnetic resonance imaging (MRI) and diffuse optical tomography (DOT) system for small animal imaging. Cross-sectional total hemoglobin concentration maps of 10 Fisher rats bearing R3230 breast tumors are reconstructed using multi-wavelength DOT measurements both with and without magnetic resonance (MR) structural a priori information. Simultaneously acquired MR structural images are used to guide and constrain the DOT reconstruction, while dynamic contrast-enhanced MR functional images are used as the gold standard to classify the vasculature of the tumor into two types: high versus low heterogeneity. These preliminary results show that the stand-alone DOT is unable to differentiate tumors with low and high vascular heterogeneity without structural a priori information provided by a high resolution imaging modality. The mean total hemoglobin concentrations comparing the vasculature of the tumors with low and high heterogeneity are significant (p-value 0.02) only when MR structural a priori information is utilized.

Evaluation of the effect of transcytolemmal water exchange analysis for therapeutic response assessment using DCE-MRI: a comparison study

This study compares the shutter-speed (SS) and the Tofts models as used in assessing therapeutic response in a longitudinal DCE-MRI experiment. Sixteen nu/nu mice with implanted colorectal adenocarcinoma cell line (LS-174T) were randomly assigned into treatment/control groups (n  =  8/group) and received bevacizumab/saline twice weekly (Day1/Day4/Day8). All mice were scanned at one pre- (Day0) and two post-treatment (Day2/Day9) time points using a high spatiotemporal resolution DCE-MRI pulse sequence. The CA extravasation rate constant [Formula: see text] from the Tofts/SS model and the mean intracellular water residence time [Formula: see text] from the SS model were analyzed. A biological subvolume (BV) within the tumor was identified based on the [Formula: see text] intensity distribution, and the SS model parameters within the BV ([Formula: see text] and [Formula: see text]) were analyzed. It is found that [Formula: see text] and [Formula: see text] have a similar spatial distribution in the tumor volume. The Bayesian information criterion results show that the SS model was a better fit for all scans. At Day9, the treatment group had significantly higher tumor mean [Formula: see text] (p  =  0.021), [Formula: see text] (p  =  0.021) and [Formula: see text] (p  = 0.045). When BV from transcytolemmal water exchange analysis was adopted, the treatment group had higher mean [Formula: see text] at both Day2 (p  =  0.038) and Day9 (p  =  0.007). Additionally, at Day9, the treatment group had higher mean [Formula: see text] (p  =  0.045) and higher [Formula: see text] spatial heterogeneity indices (Rényi dimensions) d 1 (p  = 0.010) and d 2 (p  = 0.021). When mean [Formula: see text] and its coefficient of variation (CV) were used to separate treatment/control group samples using supporting vector machine, the accuracy of treatment/control classification was 68.8% at Day2 and 87.5% at Day9; in contrast, the Day2/Day9 accuracy were 62.5%/87.5% using tumor mean [Formula: see text] and its CV and were 50.0%/87.5% using tumor mean [Formula: see text] and its CV, respectively. These results suggest that the SS model parameters outperformed the Tofts model parameters in terms of capturing bevacizumab therapeutic effect in this longitudinal experiment.

Accelerated Brain DCE-MRI Using Iterative Reconstruction With Total Generalized Variation Penalty for Quantitative Pharmacokinetic Analysis: A Feasibility Study

PURPOSE

To investigate the feasibility of using undersampled k-space data and an iterative image reconstruction method with total generalized variation penalty in the quantitative pharmacokinetic analysis for clinical brain dynamic contrast-enhanced magnetic resonance imaging.

METHODS

Eight brain dynamic contrast-enhanced magnetic resonance imaging scans were retrospectively studied. Two k-space sparse sampling strategies were designed to achieve a simulated image acquisition acceleration factor of 4. They are (1) a golden ratio-optimized 32-ray radial sampling profile and (2) a Cartesian-based random sampling profile with spatiotemporal-regularized sampling density constraints. The undersampled data were reconstructed to yield images using the investigated reconstruction technique. In quantitative pharmacokinetic analysis on a voxel-by-voxel basis, the rate constant K^trans in the extended Tofts model and blood flow F_B and blood volume V_B from the 2-compartment exchange model were analyzed. Finally, the quantitative pharmacokinetic parameters calculated from the undersampled data were compared with the corresponding calculated values from the fully sampled data. To quantify each parameter's accuracy calculated using the undersampled data, error in volume mean, total relative error, and cross-correlation were calculated.

RESULTS

The pharmacokinetic parameter maps generated from the undersampled data appeared comparable to the ones generated from the original full sampling data. Within the region of interest, most derived error in volume mean values in the region of interest was about 5% or lower, and the average error in volume mean of all parameter maps generated through either sampling strategy was about 3.54%. The average total relative error value of all parameter maps in region of interest was about 0.115, and the average cross-correlation of all parameter maps in region of interest was about 0.962. All investigated pharmacokinetic parameters had no significant differences between the result from original data and the reduced sampling data.

CONCLUSION

With sparsely sampled k-space data in simulation of accelerated acquisition by a factor of 4, the investigated dynamic contrast-enhanced magnetic resonance imaging pharmacokinetic parameters can accurately estimate the total generalized variation-based iterative image reconstruction method for reliable clinical application.

Texture analysis on parametric maps derived from dynamic contrast-enhanced magnetic resonance imaging in head and neck cancer

AIM: To investigate the merits of texture analysis on parametric maps derived from pharmacokinetic modeling with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) as imaging biomarkers for the prediction of treatment response in patients with head and neck squamous cell carcinoma (HNSCC).

METHODS: In this retrospective study, 19 HNSCC patients underwent pre- and intra-treatment DCE-MRI scans at a 1.5T MRI scanner. All patients had chemo-radiation treatment. Pharmacokinetic modeling was performed on the acquired DCE-MRI images, generating maps of volume transfer rate (K^trans) and volume fraction of the extravascular extracellular space (v_e). Image texture analysis was then employed on maps of K^trans and v_e, generating two texture measures: Energy (E) and homogeneity.

RESULTS: No significant changes were found for the mean and standard deviation for K^trans and v_e between pre- and intra-treatment (P > 0.09). Texture analysis revealed that the imaging biomarker E of v_e was significantly higher in intra-treatment scans, relative to pretreatment scans (P < 0.04).

CONCLUSION: Chemo-radiation treatment in HNSCC significantly reduces the heterogeneity of tumors.

Volumetric segmentation of ADC maps and utility of standard deviation as measure of tumor heterogeneity in soft tissue tumors

PURPOSE

Determine interobserver concordance of semiautomated three-dimensional volumetric and two-dimensional manual measurements of apparent diffusion coefficient (ADC) values in soft tissue masses (STMs) and explore standard deviation (SD) as a measure of tumor ADC heterogeneity.

RESULTS

Concordance correlation coefficients for mean ADC increased with more extensive sampling. Agreement on the SD of tumor ADC values was better for large regions of interest and multislice methods. Correlation between mean and SD ADC was low, suggesting that these parameters are relatively independent.

CONCLUSION

Mean ADC of STMs can be determined by volumetric quantification with high interobserver agreement. STM heterogeneity merits further investigation as a potential imaging biomarker that complements other functional magnetic resonance imaging parameters.

Models and methods for analyzing DCE-MRI: a review

PURPOSE

To present a review of most commonly used techniques to analyze dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), discusses their strengths and weaknesses, and outlines recent clinical applications of findings from these approaches.

METHODS

DCE-MRI allows for noninvasive quantitative analysis of contrast agent (CA) transient in soft tissues. Thus, it is an important and well-established tool to reveal microvasculature and perfusion in various clinical applications. In the last three decades, a host of nonparametric and parametric models and methods have been developed in order to quantify the CA's perfusion into tissue and estimate perfusion-related parameters (indexes) from signal- or concentration-time curves. These indexes are widely used in various clinical applications for the detection, characterization, and therapy monitoring of different diseases.

RESULTS

Promising theoretical findings and experimental results for the reviewed models and techniques in a variety of clinical applications suggest that DCE-MRI is a clinically relevant imaging modality, which can be used for early diagnosis of different diseases, such as breast and prostate cancer, renal rejection, and liver tumors.

CONCLUSIONS

Both nonparametric and parametric approaches for DCE-MRI analysis possess the ability to quantify tissue perfusion.

Childhood extracranial neoplasms: the role of imaging in drug development and clinical trials

Cancer is the leading cause of death in children older than 1 year of age and new drugs are necessary to improve outcomes. Imaging is crucial to the drug development process and assessment of therapeutic response. In adults, tumours are often assessed with CT using size criteria. Unfortunately, techniques established in adults are not necessarily applicable in children due to differing pathophysiology, ability to cooperate and increased susceptibility to ionising radiation. MRI, in particular quantitative MRI, has to date not been fully utilised in children with extracranial neoplasms. The specific challenges of imaging in children, the potential for functional imaging techniques to inform upon and their inclusion in clinical trials are discussed.