Magnetic resonance imaging of tumor response to stroma-modifying pegvorhyaluronidase alpha (PEGPH20) therapy in early-phase clinical trials.

Pre-clinical and clinical studies have shown that PEGPH20 depletes intratumoral hyaluronic acid (HA), which is linked to high interstitial fluid pressures and poor distribution of chemotherapies. 29 patients with metastatic advanced solid tumors received quantitative magnetic resonance imaging (qMRI) in 3 prospective clinical trials of PEGPH20: HALO-109-101 (NCT00834704), HALO-109-102 (NCT01170897), and HALO-109-201 (NCT01453153). Apparent Diffusion Coefficient of water (ADC), T1, ktrans, vp, ve, and iAUC maps were computed from qMRI acquired at baseline and ≥ 1 time point post-PEGPH20. Tumor ADC and T1 decreased, while iAUC, ktrans, vp, and ve increased, on day 1 post-PEGPH20 relative to baseline values. This is consistent with HA depletion leading to a decrease in tumor extracellular water content and an increase in perfusion, permeability, extracellular matrix space, and vascularity. Baseline parameter values predictive of pharmacodynamic responses were: ADC > 1.46 × 10-3 mm2/s (Balanced Accuracy (BA) = 72%, p < 0.01), T1 > 0.54 s (BA = 82%, p < 0.01), iAUC < 9.2 mM-s (BA = 76%, p < 0.05), ktrans < 0.07 min-1 (BA = 72%, p = 0.2), ve < 0.17 (BA = 68%, p < 0.01), and vp < 0.02 (BA = 60%, p < 0.01). A low ve at baseline was moderately predictive of response in any parameter (BA = 65.6%, p < 0.01 averaged across patients). These qMRI biomarkers are potentially useful for guiding patient pre-selection and post-treatment follow-up in future clinical studies of PEGPH20 and other tumor stroma-modifying anti-cancer therapies.

Tumor Response to Stroma-Modifying Therapy: Magnetic Resonance Imaging Findings in Early-Phase Clinical Trials of Pegvorhyaluronidase alpha (PEGPH20).

Pre-clinical and clinical studies have shown that PEGPH20 depletes intratumoral hyaluronic acid (HA), which is linked to high interstitial fluid pressures and poor distribution of chemotherapies. 29 patients with metastatic advanced solid tumors received quantitative magnetic resonance imaging (qMRI) in 3 prospective clinical trials of PEGPH20, HALO-109-101 (NCT00834704), HALO-109-102 (NCT01170897), and HALO-109-201 (NCT01453153). Apparent Diffusion Coefficient of water (ADC), T1, ktrans, vp, ve, and iAUC maps were computed from qMRI acquired at baseline and ≥ 1 time point post-PEGPH20. Tumor ADC and T1 decreased, while iAUC, ktrans, vp, and ve increased, on day 1 post-PEGPH20 relative to baseline values. This is consistent with HA depletion leading to a decrease in tumor water content and an increase in perfusion, permeability, extracellular matrix space, and vascularity. Baseline parameter values that were predictive of pharmacodynamic responses were: ADC > 1.46×10-3 mm2/s (Balanced Accuracy (BA) = 72%, p < 0.01), T1 > 0.54s (BA = 82%, p < 0.01), iAUC < 9.2 mM-s (BA = 76%, p < 0.05), ktrans<0.07min-1 (BA = 72%, p = 0.2), ve<0.17 (BA = 68%, p < 0.01), and vp<0.02 (BA = 60%, p < 0.01). Further, ve<0.39 at baseline was moderately predictive of response in any parameter (BA = 65.6%, p < 0.01 averaged across patients). These qMRI biomarkers are potentially useful for guiding patient pre-selection and post-treatment follow-up in future clinical studies of PEGPH20 and other tumor stroma-modifying anti-cancer therapies.

Pretherapy Ferumoxytol-enhanced MRI to Predict Response to Liposomal Irinotecan in Metastatic Breast Cancer.

Purpose To investigate ferumoxytol (FMX)-enhanced MRI as a pretreatment predictor of response to liposomal irinotecan (nal-IRI) for thoracoabdominal and brain metastases in women with metastatic breast cancer (mBC). Materials and Methods In this phase 1 expansion trial (ClinicalTrials.gov identifier, NCT01770353; 27 participants), 49 thoracoabdominal (19 participants; mean age, 48 years ± 11 [SD]) and 19 brain (seven participants; mean age, 54 years ± 8) metastases were analyzed on MR images acquired before, 1-4 hours after, and 16-24 hours after FMX administration. In thoracoabdominal metastases, tumor transverse relaxation rate (R*2) was normalized to the mean R*2 in the spleen (rR*2), and the tumor histogram metric rR*2,N, representing the average of rR*2 in voxels above the nth percentile, was computed. In brain metastases, a novel compartmentation index was derived by applying the MRI signal equation to phantom-calibrated coregistered FMX-enhanced MRI brain scans acquired before, 1-4 hours after, and 16-24 hours after FMX administration. The fraction of voxels with an FMX compartmentation index greater than 1 was computed over the whole tumor (FCIGT1) and from voxels above the 90th percentile R*2 (FCIGT1 R*2,90). Results rR*2,90 computed from pretherapy MRI performed 16-24 hours after FMX administration, without reference to calibration phantoms, predicted response to nal-IRI in thoracoabdominal metastases (accuracy, 74%). rR*2,90 performance was robust to the inclusion of some peritumoral tissue within the tumor region of interest. FCIGT1 R*2,90 provided 79% accuracy on cross-validation in prediction of response in brain metastases. Conclusion This first in-human study focused on mBC suggests that FMX-enhanced MRI biologic markers can be useful for pretherapy prediction of response to nal-IRI in patients with mBC. Keywords: MRI Contrast Agent, MRI, Breast, Head/Neck, Tumor Response, Experimental Investigations, Brain/Brain Stem Clinical trial registration no. NCT01770353 Supplemental material is available for this article. © RSNA, 2023 See also commentary by Daldrup-Link in this issue.

Assessment of fully automatic segmentation of pulmonary artery and aorta on noncontrast CT with optimal surface graph cuts.

PURPOSE: Accurate segmentation of the pulmonary arteries and aorta is important due to the association of the diameter and the shape of these vessels with several cardiovascular diseases and with the risk of exacerbations and death in patients with chronic obstructive pulmonary disease. We propose a fully automatic method based on an optimal surface graph-cut algorithm to quantify the full 3D shape and the diameters of the pulmonary arteries and aorta in noncontrast computed tomography (CT) scans. METHODS: The proposed algorithm first extracts seed points in the right and left pulmonary arteries, the pulmonary trunk, and the ascending and descending aorta by using multi-atlas registration. Subsequently, the centerlines of the pulmonary arteries and aorta are extracted by a minimum cost path tracking between the extracted seed points, with a cost based on a combination of lumen intensity similarity and multiscale medialness in three planes. The centerlines are refined by applying the path tracking algorithm to curved multiplanar reformatted scans and are then smoothed and dilated nonuniformly according to the extracted local vessel radius from the medialness filter. The resulting coarse estimates of the vessels are used as initialization for a graph-cut segmentation. Once the vessels are segmented, the diameters of the pulmonary artery (PA) and the ascending aorta (AA) and the P A : A A ratio are automatically calculated both in a single axial slice and in a 10 mm volume around the automatically extracted PA bifurcation level. The method is evaluated on noncontrast CT scans from the Danish Lung Cancer Screening Trial (DLCST). Segmentation accuracy is determined by comparing with manual annotations on 25 CT scans. Intraclass correlation (ICC) between manual and automatic diameters, both measured in axial slices at the PA bifurcation level, is computed on an additional 200 CT scans. Repeatability of the automated 3D volumetric diameter and P A : A A ratio calculations (perpendicular to the vessel axis) are evaluated on 118 scan-rescan pairs with an average in-between time of 3 months. RESULTS: We obtained a Dice segmentation overlap of 0.94 ± 0.02 for pulmonary arteries and 0.96 ± 0.01 for the aorta, with a mean surface distance of 0.62 ± 0.33 mm and 0.43 ± 0.07 mm, respectively. ICC between manual and automatic in-slice diameter measures was 0.92 for PA, 0.97 for AA, and 0.90 for the P A : A A ratio, and for automatic diameters in 3D volumes around the PA bifurcation level between scan and rescan was 0.89, 0.95, and 0.86, respectively. CONCLUSION: The proposed automatic segmentation method can reliably extract diameters of the large arteries in non-ECG-gated noncontrast CT scans such as are acquired in lung cancer screening.

Increasing Accuracy of Optimal Surfaces Using Min-Marginal Energies.

Optimal surface methods are a class of graph cut methods posing surface estimation as an n-ary ordered labeling problem. They are used in medical imaging to find interacting and layered surfaces optimally and in low order polynomial time. Representing continuous surfaces with discrete sets of labels, however, leads to discretization errors and, if graph representations are made dense, excessive memory usage. Limiting memory usage and computation time of graph cut methods are important and graphs that locally adapt to the problem has been proposed as a solution. Min-marginal energies computed using dynamic graph cuts offer a way to estimate solution uncertainty and these uncertainties have been used to decide where graphs should be adapted. Adaptive graphs, however, introduce extra parameters, complexity, and heuristics. We propose a way to use min-marginal energies to estimate continuous solution labels that does not introduce extra parameters and show empirically on synthetic and medical imaging datasets that it leads to improved accuracy. The increase in accuracy was consistent and in many cases comparable with accuracy otherwise obtained with graphs up to eight times denser, but with proportionally less memory usage and improvements in computation time.

Cooperative carotid artery centerline extraction in MRI.

Centerline extraction of the carotid artery in MRI is important to analyze the artery geometry and to provide input for further processing such as registration and segmentation. The centerline of the artery bifurcation is often extracted by means of two independent minimum cost paths ranging from the common to the internal and the external carotid artery. Often the cost is not well defined at the artery bifurcation, leading to centerline errors. To solve this problem, we developed a method to cooperatively extract both centerlines, where in the cost to extract each centerline, we integrate a constraint region derived from the estimated position of the neighbor centerline. This method avoids that both centerlines follow the same cheapest path after the bifurcation, which is a common error when the paths are extracted independently. We show that this method results in less error compared to extracting them independently: 10 failed centerlines Vs. 3 failures in a data set of 161 arteries with manual annotations. Additionally, we show that the new method improves the non-cooperative approach in 28 cases (p < 0.0001) in a data set of 3,904 arteries.

Maximization of regional probabilities using Optimal Surface Graphs: Application to carotid artery segmentation in MRI.

PURPOSE: We present a segmentation method that maximizes regional probabilities enclosed by coupled surfaces using an Optimal Surface Graph (OSG) cut approach. This OSG cut determines the globally optimal solution given a graph constructed around an initial surface. While most methods for vessel wall segmentation only use edge information, we show that maximizing regional probabilities using an OSG improves the segmentation results. We applied this to automatically segment the vessel wall of the carotid artery in magnetic resonance images. METHODS: First, voxel-wise regional probability maps were obtained using a Support Vector Machine classifier trained on local image features. Then, the OSG segments the regions which maximizes the regional probabilities considering smoothness and topological constraints. RESULTS: The method was evaluated on 49 carotid arteries from 30 subjects. The proposed method shows good accuracy with a Dice wall overlap of 74.1 ± 4.3%, and significantly outperforms a published method based on an OSG using only surface information, the obtained segmentations using voxel-wise classification alone, and another published artery wall segmentation method based on a deformable surface model. Intraclass correlations (ICC) with manually measured lumen and wall volumes were similar to those obtained between observers. Finally, we show a good reproducibility of the method with ICC = 0.86 between the volumes measured in scans repeated within a short time interval. CONCLUSIONS: In this work, a new segmentation method that uses both an OSG and regional probabilities is presented. The method shows good segmentations of the carotid artery in MRI and outperformed another segmentation method that uses OSG and edge information and the voxel-wise segmentation using the probability maps.

Carotid Artery Wall Segmentation in Multispectral MRI by Coupled Optimal Surface Graph Cuts.

We present a new three-dimensional coupled optimal surface graph-cut algorithm to segment the wall of the carotid artery bifurcation from Magnetic Resonance (MR) images. The method combines the search for both inner and outer borders into a single graph cut and uses cost functions that integrate information from multiple sequences. Our approach requires manual localization of only three seed points indicating the start and end points of the segmentation in the internal, external, and common carotid artery. We performed a quantitative validation using images of 57 carotid arteries. Dice overlap of 0.86 ± 0.06 for the complete vessel and 0.89 ± 0.05 for the lumen compared to manual annotation were obtained. Reproducibility tests were performed in 60 scans acquired with an interval of 15 ± 9 days, showing good agreement between baseline and follow-up segmentations with intraclass correlations of 0.96 and 0.74 for the lumen and complete vessel volumes respectively.

Relation between wall shear stress and carotid artery wall thickening MRI versus CFD.

Wall shear stress (WSS), a parameter associated with endothelial function, is calculated by computational fluid dynamics (CFD) or phase-contrast (PC) MRI measurements. Although CFD is common in WSS (WSSCFD) calculations, PC-MRI-based WSS (WSSMRI) is more favorable in population studies; since it is straightforward and less time consuming. However, it is not clear if WSSMRI and WSSCFD show similar associations with vascular pathology. Our aim was to test the associations between wall thickness (WT) of the carotid arteries and WSSMRI and WSSCFD. The subjects (n=14) with an asymptomatic carotid plaque who underwent MRI scans two times within 4 years of time were selected from the Rotterdam Study. We compared WSSCFD and WSSMRI at baseline and follow-up. Baseline WSSMRI and WSSCFD values were divided into 3 categories representing low, medium and high WSS tertiles. WT of each tertile was compared by a one-way ANOVA test. The WSSMRI and WSSCFD were 0.50±0.13Pa and 0.73±0.25Pa at baseline. Although WSSMRI was underestimated, a significant regression was found between WSSMRI and WSSCFD (r(2)=0.71). No significant difference was found between baseline and follow-up WSS by CFD and MRI-based calculations. The WT at baseline was 1.36±0.16mm and did not change over time. The WT was 1.55±0.21mm in low, 1.33±0.20mm in medium and 1.21±0.21mm in the high WSSMRI tertiles. Similarly, the WT was 1.49±0.21mm in low, 1.33±0.20mm in medium and 1.26±0.21mm in high WSSCFD tertiles. We found that WSSMRI and WSSCFD were inversely related with WT. WSSMRI and WSSCFD patterns were similar although MRI-based calculations underestimated WSS.