ObjectiveThe objective of this study was to non-invasively differentiate the degree of malignancy in two murine breast cancer models based on identification of distinct tissue characteristics in a metastatic and non-metastatic tumor model using a multiparametric Magnetic Resonance Imaging (MRI) approach.MethodsThe highly metastatic 4T1 breast cancer model was compared to the non-metastatic 67NR model. Imaging was conducted on a 9.4 T small animal MRI. The protocol was used to characterize tumors regarding their structural composition, including heterogeneity, intratumoral edema and hemorrhage, as well as endothelial permeability using apparent diffusion coefficient (ADC), T1/T2 mapping and dynamic contrast-enhanced (DCE) imaging. Mice were assessed on either day three, six or nine, with an i.v. injection of the albumin-binding contrast agent gadofosveset. Ex vivo validation of the results was performed with laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS), histology, immunhistochemistry and electron microscopy.ResultsSignificant differences in tumor composition were observed over time and between 4T1 and 67NR tumors. 4T1 tumors showed distorted blood vessels with a thin endothelial layer, resulting in a slower increase in signal intensity after injection of the contrast agent. Higher permeability was further reflected in higher Ktrans values, with consecutive retention of gadolinium in the tumor interstitium visible in MRI. 67NR tumors exhibited blood vessels with a thicker and more intact endothelial layer, resulting in higher peak enhancement, as well as higher maximum slope and area under the curve, but also a visible wash-out of the contrast agent and thus lower Ktrans values. A decreasing accumulation of gadolinium during tumor progression was also visible in both models in LA-ICP-MS. Tissue composition of 4T1 tumors was more heterogeneous, with intratumoral hemorrhage and necrosis and corresponding higher T1 and T2 relaxation times, while 67NR tumors mainly consisted of densely packed tumor cells. Histogram analysis of ADC showed higher values of mean ADC, histogram kurtosis, range and the 90th percentile (p90), as markers for the heterogenous structural composition of 4T1 tumors. Principal component analysis (PCA) discriminated well between the two tumor models.ConclusionsMultiparametric MRI as presented in this study enables for the estimation of malignant potential in the two studied tumor models via the assessment of certain tumor features over time.