Comparison of Voxel-Wise Tumor Perfusion Changes Measured With Dynamic Contrast-Enhanced (DCE) MRI and Volumetric DCE CT in Patients With Metastatic Brain Cancer Treated with Radiosurgery

Spatial correlation between in vivo imaging and immunohistochemical biomarkers: A methodological study

In this study, we present a method that enables voxel-by-voxel comparison of in vivo imaging to immunohistochemistry (IHC) biomarkers. As a proof of concept, we investigated the spatial correlation between dynamic contrast enhanced (DCE-)CT parameters and IHC biomarkers Ki-67 (proliferation), HIF-1α (hypoxia), and CD45 (immune cells). 54 whole-mount tumor slices of 15 laryngeal and hypopharyngeal carcinomas were immunohistochemically stained and digitized. Heatmaps of biomarker positivity were created and registered to DCE-CT parameter maps. The adiabatic approximation to the tissue homogeneity model was used to fit the following DCE parameters: Ktrans (transfer constant), Ve (extravascular and extracellular space), and Vi (intravascular space). Both IHC and DCE maps were downsampled to 4 × 4 × 3 mm[3] voxels. The mean values per tumor were used to calculate the between-subject correlations between parameters. For the within-subject (spatial) correlation, values of all voxels within a tumor were compared using the repeated measures correlation (rrm). No between-subject correlations were found between IHC biomarkers and DCE parameters, whereas we found multiple significant within-subject correlations: Ve and Ki-67 (rrm = -0.17, P < .001), Ve and HIF-1α (rrm = -0.12, P < .001), Ktrans and CD45 (rrm = 0.13, P < .001), Vi and CD45 (rrm = 0.16, P < .001), and Vi and Ki-67 (rrm = 0.08, P = .003). The strongest correlation was found between IHC biomarkers Ki-67 and HIF-1α (rrm = 0.35, P < .001). This study shows the technical feasibility of determining the 3 dimensional spatial correlation between histopathological biomarker heatmaps and in vivo imaging. It also shows that between-subject correlations do not reflect within-subject correlations of parameters.

Comparison of computed tomography perfusion and magnetic resonance dynamic susceptibility contrast perfusion-weighted imaging in canine brain

Brain perfusion allows for the evaluation of cerebral hemodynamics, particularly in brain infarcts and tumors. Computed tomography (CT) perfusion (CTP) provides reliable data; however, it has a limited scan field of view and radiation exposure. Magnetic resonance (MR) perfusion provides detailed imaging of small structures and a wide scan field of view. However, no study has compared CTP and MR perfusion and assessed the correlation between the perfusion parameters measured using CTP and MR perfusion. The aim of the present study was to assess the correlation and agreement of the cerebral perfusion derived from dynamic susceptibility contrast (DSC)-MRI and CTP in dogs. In this crossover design study, the cerebral blood volume (CBV), cerebral blood flow (CBF), mean transit time, and time to peak were measured in the temporal cerebral cortex, caudate nucleus, thalamus, piriform lobe, and hippocampus using CTP and DSC-MRI in six healthy beagle dogs and a dog with a pituitary tumor. On the color map of healthy beagles, blood vessels and the perivascular brain parenchyma appeared as red-green, indicating high perfusion, and the areas distant from the vessels appeared as green-blue, indicating low perfusion levels in CTP and DSC-MRI. CTP parameters were highest in the piriform lobe (CBF = 121.11 ± 12.78 mL/100 g/min and CBV = 8.70 ± 2.04 mL/100 g) and lowest in the thalamus (CBF = 63.75 ± 25.24 mL/100 g/min and CBV = 4.02 ± 0.55 mL/100 g). DSC-MRI parameters were also highest in the piriform lobe (CBF = 102.31 ± 14.73 mL/100 g/min and CBV = 3.17 ± 1.23 mL/100 g) and lowest in the thalamus (CBF = 37.73 ± 25.11 mL/100 g/min and CBV = 0.81 ± 0.44 mL/100 g) although there was no statistical correlation in the quantitative perfusion parameters between CTP and DSC-MRI. In a dog with a pituitary tumor, the color map of the tumor appeared as a red scale, indicating high perfusion and higher CBF and CBV on CTP (149 mL/100 g and 20 mL/100 g/min) and on DSC-MRI (116.3 mL/100 g and 15.32 mL/100 g/min) compared to those measured in healthy dogs. These findings indicate that DSC-MRI and CTP maps exhibit comparability and interchangeability in the assessment of canine brain perfusion.

Incorporating cross-voxel exchange for the analysis of dynamic contrast-enhanced imaging data: pre-clinical results

Tumours exhibit abnormal interstitial structures and vasculature function often leading to impaired and heterogeneous drug delivery. The disproportionate spatial accumulation of a drug in the interstitium is determined by several microenvironmental properties (blood vessel distribution and permeability, gradients in the interstitial fluid pressure). Predictions of tumour perfusion are key determinants of drug delivery and responsiveness to therapy. Pharmacokinetic models allow for the quantification of tracer perfusion based on contrast enhancement measured with non-invasive imaging techniques. An advanced cross-voxel exchange model (CVXM) was recently developed to provide a comprehensive description of tracer extravasation as well as advection and diffusion based on cross-voxel tracer kinetics (Sinnoet al2021). Transport parameters were derived from DCE-MRI of twenty TS-415 human cervical carcinoma xenografts by using CVXM. Tracer velocity flows were measured at the tumour periphery (mean 1.78-5.82μm.s^-1) pushing the contrast outward towards normal tissue. These elevated velocity measures and extravasation rates explain the heterogeneous distribution of tracer across the tumour and its accumulation at the periphery. Significant values for diffusivity were deduced across the tumours (mean 152-499μm².s^-1). CVXM resulted in generally smaller values for the extravasation parameterKext(mean 0.01-0.04 min^-1) and extravascular extracellular volume fractionve(mean 0.05-0.17) compared to the standard Tofts parameters, suggesting that Toft model underestimates the effects of inter-voxel exchange. The ratio of Tofts' extravasation parameters over CVXM's was significantly positively correlated to the cross-voxel diffusivity (P< 0.0001) and velocity (P= 0.0005). Tofts' increasedvemeasurements were explained using Sinnoet al(2021)'s theoretical work. Finally, a scan time of 15 min renders informative estimations of the transport parameters. However, a duration as low as 7.5 min is acceptable to recognize the spatial variation of transport parameters. The results demonstrate the potential of utilizing CVXM for determining metrics characterizing the exchange of tracer between the vasculature and the tumour tissue. Like for many earlier models, additional work is strongly recommended, in terms of validation, to develop more confidence in the results, motivating future laboratory work in this regard.

Assessment of brain tumors by magnetic resonance dynamic susceptibility contrast perfusion-weighted imaging and computed tomography perfusion: a comparison study

PURPOSE

To investigate the association and agreement between magnetic resonance dynamic susceptibility contrast perfusion-weighted imaging (DSC-PWI) and computed tomography perfusion (CTP) in determining vascularity and permeability of primary and secondary brain tumors.

MATERIAL AND METHODS

DSC-PWI and CTP studies from 97 patients with high-grade glioma, low-grade glioma and solitary brain metastasis were retrospectively reviewed. Normalized cerebral blood flow (nCBF), cerebral blood volume (nCBV), capillary transfer constant (nK2) and permeability surface area product (nPS) values were obtained. Variables among groups were compared, and correlation and agreement between DSC-PWI and CTP were tested.

RESULTS

All DSC-PWI and CTP parameters were higher in high-grade than in low-grade gliomas (p < 0.01 and p < 0.001). Metastases had greater DSC-PWI nCBV (p < 0.05), nCTP-CBF (p < 0.05), nCTP-CBV (p < 0.01) and nCTP-PS (p < 0.0001) than low-grade gliomas and more elevated nCTP-PS (p < 0.01) than high-grade gliomas. The correlation was strong between DSC-PWI nCBF and CTP nCBF (r = 0.79; p < 0.00001) and between DSC-PWI nCBV and CTP nCBV (r = 0.83; p < 0.00001), weaker between DSC-PWI nK2 and CTP nPS (r = 0.29; p < 0.01). Bland-Altman plots indicated that the agreement was strong between DSC-PWI nCBF and CTP nCBF, good between DSC-PWI nCBV and CTP nCBV and poorer between DSC-PWI nK2 and CTP nPS.

CONCLUSION

DSC-PWI and CTP CBF and CBV maps were comparable and interchangeable in the assessment of tumor vascularity, unlike DSC-PWI K2 and CTP PS maps that were more discordant in the analysis of tumor permeability. CTP could be an alternative method to quantify tumor neoangiogenesis when MRI is not available or when the patient does not tolerate it.

Incorporating cross-voxel exchange into the analysis of dynamic contrast-enhanced imaging data: theory, simulations and experimental results

Predictions of tumour perfusion are key determinants of drug delivery and responsiveness to therapy. Pharmacokinetic models allow for the estimation of perfusion properties of tumour tissues but many assume no dispersion associated with tracer transport away from the capillaries and through the tissue. At the level of a voxel, this translates to assuming no cross-voxel tracer exchange, often leading to the misinterpretation of derived perfusion parameters. Tofts model (TM), a compartmental model widely used in oncology, also makes this assumption. A more realistic description is required to quantify kinetic properties of tracers, such as convection and diffusion. We propose a Cross-Voxel Exchange Model (CVXM) for analysing cross-voxel tracer kinetics.In silicodatasets quantifying the roles of convection and diffusion in tracer transport (which TM ignores) were employed to investigate the interpretation of Tofts' perfusion parameters compared to CVXM. TM returned inaccurate values ofKtransandvewhere diffusive and convective mechanisms are pronounced (up to 20% and 300% error respectively). A mathematical equation, developed in this work, predicts and gives the correct physiological interpretation of Tofts've.Finally, transport parameters were derived from dynamic contrast enhanced-magnetic resonance imaging of a TS-415 human cervical carcinoma xenograft by using TM and CVXM. The latter deduced lower values ofKtransandvecompared to TM (lower by up to 63% and 76% respectively). It also allowed the detection of a diffusive flux (mean diffusivity 155μm²s^-1) in the tumour tissue, as well as an increased convective flow at the periphery (mean velocity 2.3μm s^-1detected). The results serve as a proof of concept establishing the feasibility of using CVXM for accurately determining transport metrics that characterize the exchange of tracer between voxels. CVXM needs to be investigated further as its parameters can be linked to the tumour microenvironment properties (permeability, pressure…), potentially leading to enhanced personalized treatment planning.

Multisite Technical and Clinical Performance Evaluation of Quantitative Imaging Biomarkers from 3D FDG PET Segmentations of Head and Neck Cancer Images

Quantitative imaging biomarkers (QIBs) provide medical image-derived intensity, texture, shape, and size features that may help characterize cancerous tumors and predict clinical outcomes. Successful clinical translation of QIBs depends on the robustness of their measurements. Biomarkers derived from positron emission tomography images are prone to measurement errors owing to differences in image processing factors such as the tumor segmentation method used to define volumes of interest over which to calculate QIBs. We illustrate a new Bayesian statistical approach to characterize the robustness of QIBs to different processing factors. Study data consist of 22 QIBs measured on 47 head and neck tumors in 10 positron emission tomography/computed tomography scans segmented manually and with semiautomated methods used by 7 institutional members of the NCI Quantitative Imaging Network. QIB performance is estimated and compared across institutions with respect to measurement errors and power to recover statistical associations with clinical outcomes. Analysis findings summarize the performance impact of different segmentation methods used by Quantitative Imaging Network members. Robustness of some advanced biomarkers was found to be similar to conventional markers, such as maximum standardized uptake value. Such similarities support current pursuits to better characterize disease and predict outcomes by developing QIBs that use more imaging information and are robust to different processing factors. Nevertheless, to ensure reproducibility of QIB measurements and measures of association with clinical outcomes, errors owing to segmentation methods need to be reduced.

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.

Routine or enhanced imaging to differentiate between type 1 and type 2 papillary renal cell carcinoma

AIM

To investigate imaging features and differentiating qualities of type 1 and type 2 papillary renal cell carcinoma (pRCC) by different imaging techniques.

MATERIALS AND METHODS

From 2007 to 2019, 107 patients with type 1 pRCC (T1-pRCC) and 147 with type 2 pRCC (T2-pRCC) were included in this retrospective study. All patients underwent conventional ultrasound (US); some also underwent contrast-enhanced ultrasound (CEUS), contrast-enhanced computed tomography (CECT), or contrast-enhanced magnetic resonance imaging (CE-MRI). Tumour Fuhrman grade or World Health Organization (WHO)/International Society of Urological Pathology (ISUP) grade (after June 2016) and invasive ranges were recorded. The two types of pRCC were analysed and compared for imaging features including tumour position, size, margin, echo type, and colour Doppler flow imaging (CDFI) using US as well as enhanced features from CEUS, CECT, or CE-MRI.

RESULTS

T2-pRCC showed a higher Fuhrman grade (p<0.001) and greater propensity to invade extrarenal tissue (p<0.001) than T1-pRCC. On US imaging, T2-pRCC was more likely to be a cystic-solid lesion (p<0.001), and colour flow with a higher resistance index (RI; p=0.014) was more easily detected (p=0.001) in T2-pRCC than in T1-pRCC. Within contrast-enhanced examinations, more T2-pRCC lesions had blurred tumour borders (p=0.003), hypervascular characteristics (p=0.003), and heterogeneous enhancement (p<0.001) than those of T1-pRCC.

CONCLUSIONS

T2-pRCC manifests more aggressively than T1-pRCC. T2-pRCC has a higher proportion of hypervascular and heterogeneous enhancement than T1-RCC.

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.

Accuracy and Performance of Functional Parameter Estimation Using a Novel Numerical Optimization Approach for GPU-Based Kinetic Compartmental Modeling

Quantitative kinetic parameters derived from dynamic contrast-enhanced (DCE) data are dependent on signal measurement quality and choice of pharmacokinetic model. However, the fundamental optimization analysis method is equally important and its impact on pharmacokinetic parameters has been mostly overlooked. We examine the effects of those choices on accuracy and performance of parameter estimation using both computer processing unit and graphical processing unit (GPU) numerical optimization implementations and evaluate the improvements offered by a novel optimization approach. A test framework was developed where experimentally derived population-average arterial input function and randomly sampled parameter sets {K_trans, K_ep, V_b, τ} were used to generate known tissue curves. Five numerical optimization algorithms were evaluated: sequential quadratic programming, downhill simplex (Nelder–Mead), pattern search, simulated annealing, and differential evolution. This was combined with various objective function implementation details: delay approximation, discretization and varying sampling rates. Then, impact of noise and CPU/GPU implementation was tested for speed and accuracy. Finally, the optimal method was compared to conventional implementation as applied to clinical DCE computed tomography. Nelder–Mead, differential evolution and sequential quadratic programming produced good results on clean and noisy input data outperforming simulated annealing and pattern search in terms of speed and accuracy in the respective order of 10⁻⁸%, 10⁻⁷%, and ×10⁻⁶%). A novel approach for DCE numerical optimization (infinite impulse response with fractional delay approximation) was implemented on GPU for speed increase of at least 2 orders of magnitude. Applied to clinical data, the magnitude of overall parameter error was <10%.

Unified platform for multimodal voxel-based analysis to evaluate tumour perfusion and diffusion characteristics before and after radiation treatment evaluated in metastatic brain cancer

OBJECTIVE:

Early changes in tumour behaviour following stereotactic radiosurgery) are potential biomarkers of response. To-date quantitative model-based measures of dynamic contrast-enhanced (DCE) and diffusion-weighted (DW) MRI parameters have shown widely variable findings, which may be attributable to variability in image acquisition, post-processing and analysis. Big data analytic approaches are needed for the automation of computationally intensive modelling calculations for every voxel, independent of observer interpretation.

METHODS:

This unified platform is a voxel-based, multimodality architecture that brings complimentary solute transport processes such as perfusion and diffusion into a common framework. The methodology was tested on synthetic data and digital reference objects and consequently evaluated in patients who underwent volumetric DCE-CT, DCE-MRI and DWI-MRI scans before and after treatment. Three-dimensional pharmacokinetic parameter maps from both modalities were compared as well as the correlation between apparent diffusion coefficient (ADC) values and the extravascular, extracellular volume (V_e). Comparison of histogram parameters was done via Bland-Altman analysis, as well as Student's t-test and Pearson's correlation using two-sided analysis.

RESULTS:

System testing on synthetic Tofts model data and digital reference objects recovered the ground truth parameters with mean relative percent error of 1.07 × 10^-7 and 5.60 × 10^-4 respectively. Direct voxel-to-voxel Pearson's analysis showed statistically significant correlations between CT and MR which peaked at Day 7 for K_trans (R = 0.74, p <= 0.0001). Statistically significant correlations were also present between ADC and Ve derived from both DCE-MRI and DCE-CT with highest median correlations found at Day 3 between median ADC and Ve,_MRI values (R = 0.6, p < 0.01) The strongest correlation to DCE-CT measurements was found with DCE-MRI analysis using voxelwise T₁₀ maps (R = 0.575, p < 0.001) instead of assigning a fixed T₁₀ value.

CONCLUSION:

The unified implementation of multiparametric transport modelling allowed for more robust and timely observer-independent data analytics. Utility of a common analysis platform has shown higher correlations between pharmacokinetic parameters obtained from different modalities than has previously been reported.

ADVANCES IN KNOWLEDGE:

Utility of a common analysis platform has shown statistically higher correlations between pharmacokinetic parameters obtained from different modalities than has previously been reported.

Detectability of radiation-induced changes in magnetic resonance biomarkers following stereotactic radiosurgery: A pilot study

Our objective was to investigate direct voxel-wise relationship between dose and early MR biomarker changes both within and in the high-dose region surrounding brain metastases following stereotactic radiosurgery (SRS). Specifically, we examined the apparent diffusion coefficient (ADC) from diffusion-weighted imaging and the contrast transfer coefficient (K^trans) and volume of extracellular extravascular space (v_e) derived from dynamic contrast-enhanced (DCE) MRI data. We investigated 29 brain metastases in 18 patients using 3 T MRI to collect imaging data at day 0, day 3 and day 20 following SRS. The ADC maps were generated by the scanner and K^trans and v_e maps were generated using in-house software for dynamic tracer-kinetic analysis. To enable spatially-correlated voxel-wise analysis, we developed a registration pipeline to register all ADC, K^trans and v_e maps to the planning MRI scan. To interrogate longitudinal changes, we computed absolute ΔADC, ΔK^trans and Δv_e for day 3 and 20 post-SRS relative to day 0. We performed a Kruskall-Wallice test on each biomarker between time points and investigated dose correlations within the gross tumour volume (GTV) and surrounding high dose region > 12 Gy via Spearman’s rho. Only v_e exhibited significant differences between day 0 and 20 (p < 0.005) and day 3 and 20 (p < 0.05) within the GTV following SRS. Strongest dose correlations were observed for ADC within the GTV (rho = 0.17 to 0.20) and weak correlations were observed for ADC and K^trans in the surrounding > 12 Gy region. Both ΔK^trans and Δv_e showed a trend with dose at day 20 within the GTV and > 12 Gy region (rho = -0.04 to -0.16). Weak dose-related decreases in K^trans and v_e within the GTV and high dose region at day 20 most likely reflect underlying vascular responses to radiation. Our study also provides a voxel-wise analysis schema for future MR biomarker studies with the goal of elucidating surrogates for radionecrosis.

The Use of Quantitative Imaging in Radiation Oncology: A Quantitative Imaging Network (QIN) Perspective

Modern radiation therapy is delivered with great precision, in part by relying on high-resolution multidimensional anatomic imaging to define targets in space and time. The development of quantitative imaging (QI) modalities capable of monitoring biologic parameters could provide deeper insight into tumor biology and facilitate more personalized clinical decision-making. The Quantitative Imaging Network (QIN) was established by the National Cancer Institute to advance and validate these QI modalities in the context of oncology clinical trials. In particular, the QIN has significant interest in the application of QI to widen the therapeutic window of radiation therapy. QI modalities have great promise in radiation oncology and will help address significant clinical needs, including finer prognostication, more specific target delineation, reduction of normal tissue toxicity, identification of radioresistant disease, and clearer interpretation of treatment response. Patient-specific QI is being incorporated into radiation treatment design in ways such as dose escalation and adaptive replanning, with the intent of improving outcomes while lessening treatment morbidities. This review discusses the current vision of the QIN, current areas of investigation, and how the QIN hopes to enhance the integration of QI into the practice of radiation oncology.