Acute Temporal Changes of MRI-Tracked Tumor Vascular Parameters after Combined Anti-angiogenic and Radiation Treatments in a Rat Glioma Model: Identifying Signatures of Synergism

Radiomics characterization of tissues in an animal brain tumor model imaged using dynamic contrast enhanced (DCE) MRI

Here, we investigate radiomics-based characterization of tumor vascular and microenvironmental properties in an orthotopic rat brain tumor model measured using dynamic-contrast-enhanced (DCE) MRI. Thirty-two immune compromised-RNU rats implanted with human U-251N cancer cells were imaged using DCE-MRI (7Tesla, Dual-Gradient-Echo). The aim was to perform pharmacokinetic analysis using a nested model (NM) selection technique to classify brain regions according to vasculature properties considered as the source of truth. A two-dimensional convolutional-based radiomics analysis was performed on the raw-DCE-MRI of the rat brains to generate dynamic radiomics maps. The raw-DCE-MRI and respective radiomics maps were used to build 28 unsupervised Kohonen self-organizing-maps (K-SOMs). A Silhouette-Coefficient (SC), k-fold Nested-Cross-Validation (k-fold-NCV), and feature engineering analyses were performed on the K-SOMs' feature spaces to quantify the distinction power of radiomics features compared to raw-DCE-MRI for classification of different Nested Models. Results showed that eight radiomics features outperformed respective raw-DCE-MRI in prediction of the three nested models. The average percent difference in SCs between radiomics features and raw-DCE-MRI was: 29.875% ± 12.922%, p < 0.001. This work establishes an important first step toward spatiotemporal characterization of brain regions using radiomics signatures, which is fundamental toward staging of tumors and evaluation of tumor response to different treatments.

Dynamic contrast enhanced (DCE) MRI estimation of vascular parameters using knowledge-based adaptive models

We introduce and validate four adaptive models (AMs) to perform a physiologically based Nested-Model-Selection (NMS) estimation of such microvascular parameters as forward volumetric transfer constant, Ktrans, plasma volume fraction, vp, and extravascular, extracellular space, ve, directly from Dynamic Contrast-Enhanced (DCE) MRI raw information without the need for an Arterial-Input Function (AIF). In sixty-six immune-compromised-RNU rats implanted with human U-251 cancer cells, DCE-MRI studies estimated pharmacokinetic (PK) parameters using a group-averaged radiological AIF and an extended Patlak-based NMS paradigm. One-hundred-ninety features extracted from raw DCE-MRI information were used to construct and validate (nested-cross-validation, NCV) four AMs for estimation of model-based regions and their three PK parameters. An NMS-based a priori knowledge was used to fine-tune the AMs to improve their performance. Compared to the conventional analysis, AMs produced stable maps of vascular parameters and nested-model regions less impacted by AIF-dispersion. The performance (Correlation coefficient and Adjusted R-squared for NCV test cohorts) of the AMs were: 0.914/0.834, 0.825/0.720, 0.938/0.880, and 0.890/0.792 for predictions of nested model regions, vp, Ktrans, and ve, respectively. This study demonstrates an application of AMs that quickens and improves DCE-MRI based quantification of microvasculature properties of tumors and normal tissues relative to conventional approaches.

Characterization of the Response of 9L and U-251N Orthotopic Brain Tumors to 3D Conformal Radiation Therapy

In a study employing MRI-guided stereotactic radiotherapy (SRS) in two orthotopic rodent brain tumor models, the radiation dose yielding 50% survival (the TCD50) was sought. Syngeneic 9L cells, or human U-251N cells, were implanted stereotactically in 136 Fischer 344 rats or 98 RNU athymic rats, respectively. At approximately 7 days after implantation for 9L, and 18 days for U-251N, rats were imaged with contrast-enhanced MRI (CE-MRI) and then irradiated using a Small Animal Radiation Research Platform (SARRP) operating at 220 kV and 13 mA with an effective energy of ∼70 keV and dose rate of ∼2.5 Gy per min. Radiation doses were delivered as single fractions. Cone-beam CT images were acquired before irradiation, and tumor volumes were defined using co-registered CE-MRI images. Treatment planning using MuriPlan software defined four non-coplanar arcs with an identical isocenter, subsequently accomplished by the SARRP. Thus, the treatment workflow emulated that of current clinical practice. The study endpoint was animal survival to 200 days. The TCD50 inferred from Kaplan-Meier survival estimation was approximately 25 Gy for 9L tumors and below 20 Gy, but within the 95% confidence interval in U-251N tumors. Cox proportional-hazards modeling did not suggest an effect of sex, with the caveat of wide confidence intervals. Having identified the radiation dose at which approximately half of a group of animals was cured, the biological parameters that accompany radiation response can be examined.

Dual integrin αvβ3 and αvβ5 blockade attenuates cardiac dysfunction by reducing fibrosis in a rat model of doxorubicin-induced cardiomyopathy

OBJECTIVE

The present study aimed to evaluate the protective role of cilengitide (CGT), an integrin αvβ3 and αvβ5 inhibitor, on doxorubicin (DOX)-induced myocardial fibrosis and cardiac dysfunction in a rat model. Methods. Forty male rats were randomly divided into four groups: DOX (n = 12), intraperitoneal (i.p.) injection of DOX 0.8 ∼ 1.0 mg/kg three times a week for up to 6 weeks, then saline i.p. three times a week for another 3 weeks; CGT (n = 8), CGT 10 mg/kg, i.p. three times a week for 9 weeks; DOX + CGT (n = 12), DOX and CGT co-administration as above for 6 weeks, then CGT alone for another 3 weeks; Control (n = 8), saline i.p. three times a week for 9 weeks. Echocardiography, serum procollagen I C-terminal propeptide (PICP) procollagen III N-terminal propeptide (PIIINP) and C telopeptide type I (CTX-I) were evaluated at baseline and 3, 6 and 9 weeks after initial DOX administration for all surviving rats. The heart tissues were then harvested for myocardial hydroxyproline (HYP) evaluation, qRT-PCR, and western blotting. Results. CGT attenuated DOX-induced eccentric remodeling by improving relative wall thickness at the 9th week. CGT also improved systolic function at the 9th week and diastolic function at the 6th and 9th week. CGT reduced myocardial HYP and serum PICP, PIIINP, CTX-I, and the PICP/PIIINP ratio. RT-PCR and western blot showed that CGT blocked the TGF-β1/SMAD3 pathway and mitigating extracellular matrix turnover. Conclusions. CGT exerted a cardioprotective effect against doxorubicin-induced fibrosis and improved cardiac function.

A computational model of glioma reveals opposing, stiffness-sensitive effects of leaky vasculature and tumor growth on tissue mechanical stress and porosity

A biphasic computational model of a growing, vascularized glioma within brain tissue was developed to account for unique features of gliomas, including soft surrounding brain tissue, their low stiffness relative to brain tissue, and a lack of draining lymphatics. This model is the first to couple nonlinear tissue deformation with porosity and tissue hydraulic conductivity to study the mechanical interaction of leaky vasculature and solid growth in an embedded glioma. The present model showed that leaky vasculature and elevated interstitial fluid pressure produce tensile stress within the tumor in opposition to the compressive stress produced by tumor growth. This tensile effect was more pronounced in softer tissue and resulted in a compressive stress concentration at the tumor rim that increased when tumor was softer than host. Aside from generating solid stress, fluid pressure-driven tissue deformation decreased the effective stiffness of the tumor while growth increased it, potentially leading to elevated stiffness in the tumor rim. A novel prediction of reduced porosity at the tumor rim was corroborated by direct comparison with estimates from our in vivo imaging studies. Antiangiogenic and radiation therapy were simulated by varying vascular leakiness and tissue hydraulic conductivity. These led to greater solid compression and interstitial pressure in the tumor, respectively, the former of which may promote tumor infiltration of the host. Our findings suggest that vascular leakiness has an important influence on in vivo solid stress, stiffness, and porosity fields in gliomas given their unique mechanical microenvironment.

Imaging acute effects of bevacizumab on tumor vascular kinetics in a preclinical orthotopic model of U251 glioma

The effect of a human vascular endothelial growth factor antibody on the vasculature of human tumor grown in rat brain was studied. Using dynamic contrast-enhanced magnetic resonance imaging, the effects of intravenous bevacizumab (Avastin; 10 mg/kg) were examined before and at postadministration times of 1, 2, 4, 8, 12 and 24 h (N = 26; 4-5 per time point) in a rat model of orthotopic, U251 glioblastoma (GBM). The commonly estimated vascular parameters for an MR contrast agent were: (i) plasma distribution volume (v_p ), (ii) forward volumetric transfer constant (K^trans ) and (iii) reverse transfer constant (k_ep ). In addition, extracellular distribution volume (V_D ) was estimated in the tumor (V_D-tumor ), tumor edge (V_D-edge ) and the mostly normal tumor periphery (V_D-peri ), along with tumor blood flow (TBF), peri-tumoral hydraulic conductivity (K) and interstitial flow (Flux) and tumor interstitial fluid pressure (TIFP). Studied as % changes from baseline, the 2-h post-treatment time point began showing significant decreases in v_p , V_D-tumor, V_D-edge and V_D-peri , as well as K, with these changes persisting at 4 and 8 h in v_p , K, V_{D-tumor, -edge} and _-peri (t-tests; p < 0.05-0.01). Decreases in K^trans were observed at the 2- and 4-h time points (p < 0.05), while interstitial volume fraction (v_e ; = K^trans /k_ep ) showed a significant decrease only at the 2-h time point (p < 0.05). Sustained decreases in Flux were observed from 2 to 24 h (p < 0.01) while TBF and TIFP showed delayed responses, increases in the former at 12 and 24 h and a decrease in the latter only at 12 h. These imaging biomarkers of tumor vascular kinetics describe the short-term temporal changes in physical spaces and fluid flows in a model of GBM after Avastin administration.

The impact of initial tumor microenvironment on imaging phenotype

Models of human cancer, to be useful, must replicate human disease with high fidelity. Our focus in this study is rat xenograft brain tumors as a model of human embedded cerebral tumors. A distinguishing signature of such tumors in humans, that of contrast-enhancement on imaging, is often not present when the human cells grow in rodents, despite the xenografts having nearly identical DNA signatures to the original tumor specimen. Although contrast enhancement was uniformly evident in all the human tumors from which the xenografts’ cells were derived, we show that long-term contrast enhancement in the model tumors may be determined conditionally by the tumor microenvironment at the time of cell implantation. We demonstrate this phenomenon in one of two patient-derived orthotopic xenograft (PDOX) models using cancer stem-like cell (CSC)-enriched neurospheres from human tumor resection specimens, transplanted to groups of immune-compromised rats in the presence or absence of a collagen/fibrin scaffolding matrix, Matrigel. The rats were imaged by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and their brains were examined by histopathology. Targeted proteomics of the PDOX tumor specimens grown from CSC implanted with and without Matrigel showed that while the levels of the majority of proteins and post-translational modifications were comparable between contrast-enhancing and non-enhancing tumors, phosphorylation of Fox038 showed a differential expression. The results suggest key proteins determine contrast enhancement and suggest a path toward the development of better animal models of human glioma. Future work is needed to elucidate fully the molecular determinants of contrast-enhancement.

Glioma-Targeted Delivery of a Theranostic Liposome Integrated with Quantum Dots, Superparamagnetic Iron Oxide, and Cilengitide for Dual-Imaging Guiding Cancer Surgery

Herein, a theranostic liposome (QSC-Lip) integrated with superparamagnetic iron oxide nanoparticles (SPIONs) and quantum dots (QDs) and cilengitide (CGT) into one platform is constructed to target glioma under magnetic targeting (MT) for guiding surgical resection of glioma. Transmission electron microscopy and X-ray photoelectron spectroscopy confirm the complete coencapsulation of SPIONs and QDs in liposome. Besides, CGT is also effectively encapsulated into the liposome with an encapsulation efficiency of ∼88.9%. QSC-Lip exhibits a diameter of 100 ± 1.24 nm, zeta potential of -17.10 ± 0.11 mV, and good stability in several mediums. Moreover, each cargo shows a biphasic release pattern from QSC-Lip, a rapid initial release within initial 10 h followed by a sustained release. Cellular uptake of QSC-Lip is significantly enhanced by C₆ cells under MT. In vivo dual-imaging studies show that QSC-Lip not only produces an obvious negative-contrast enhancement effect on glioma by magnetic resonance imaging but also makes tumor emitting fluorescence under MT. The dual-imaging of QSC-Lip guides the accurate resection of glioma by surgery. Besides, CGT is also specifically distributed to glioma after administration of QSC-Lip under MT, resulting in an effective inhibition of tumors. The integrated liposome may be a potential carrier for theranostics of tumor.

Toward a noninvasive estimate of interstitial fluid pressure by dynamic contrast-enhanced MRI in a rat model of cerebral tumor

PURPOSE

This study demonstrates a DCE-MRI estimate of tumor interstitial fluid pressure (TIFP) and hydraulic conductivity in a rat model of glioblastoma, with validation against an invasive wick-in-needle (WIN) technique. An elevated TIFP is considered a mark of aggressiveness, and a decreased TIFP a predictor of response to therapy.

METHODS

The DCE-MRI studies were conducted in 36 athymic rats (controls and posttreatment animals) with implanted U251 cerebral tumors, and with TIFP measured using a WIN method. Using a model selection paradigm and a novel application of Patlak and Logan plots to DCE-MRI data, the MRI parameters required for estimating TIFP noninvasively were estimated. Two models, a fluid-mechanical model and a multivariate empirical model, were used for estimating TIFP, as verified against WIN-TIFP.

RESULTS

Using DCE-MRI, the mean estimated hydraulic conductivity (MRI-K) in U251 tumors was (2.3 ± 3.1) × 10^-5 (mm² /mmHg-s) in control studies. Significant positive correlations were found between WIN-TIFP and MRI-TIFP in both mechanical and empirical models. For instance, in the control group of the fluid-mechanical model, MRI-TIFP was a strong predictor of WIN-TIFP (R²  = 0.76, p < .0001). A similar result was found in the bevacizumab-treated group of the empirical model (R²  = 0.93, p = .014).

CONCLUSION

This research suggests that MRI dynamic studies contain enough information to noninvasively estimate TIFP in this, and possibly other, tumor models, and thus might be used to assess tumor aggressiveness and response to therapy.

Reproducibility and relative stability in magnetic resonance imaging indices of tumor vascular physiology over a period of 24h in a rat 9L gliosarcoma model

PURPOSE

The objective was to study temporal changes in tumor vascular physiological indices in a period of 24h in a 9L gliosarcoma rat model.

METHODS

Fischer-344 rats (N=14) were orthotopically implanted with 9L cells. At 2weeks post-implantation, they were imaged twice in a 24h interval using dynamic contrast enhanced magnetic resonance imaging (DCE-MRI). Data-driven model-selection-based analysis was used to segment tumor regions with varying vascular permeability characteristics. The region with the maximum number of estimable parameters of vascular kinetics was chosen for comparison across the two time points. It provided estimates of three parameters for an MR contrast agent (MRCA): i) plasma volume (v_p), ii) forward volumetric transfer constant (K^trans) and interstitial volume fraction (v_e, ratio of K^trans to reverse transfer constant, k_ep). In addition, MRCA extracellular distribution volume (V_D) was estimated in the tumor and its borders, along with tumor blood flow (TBF) and peritumoral MRCA flux. Descriptors of parametric distributions were compared between the two times. Tumor extent was examined by hematoxylin and eosin (H&E) staining. Picrosirus red staining of secreted collagen was performed as an additional index for 9L cells.

RESULTS

Test-retest differences between population summaries for any parameter were not significant (paired t and Wilcoxon signed rank tests). Bland-Altman plots showed no apparent trends between the differences and averages of the test-retest measures for all indices. The intraclass correlation coefficients showed moderate to almost perfect reproducibility for all of the parameters, except v_p. H&E staining showed tumor infiltration in parenchyma, perivascular space and white matter tracts. Collagen staining was observed along the outer edges of main tumor mass.

CONCLUSION

The data suggest the relative stability of these MR indices of tumor microenvironment over a 24h duration in this gliosarcoma model.