Tumour size measurement in a mouse model using high resolution MRI
© Montelius et al. licensee BioMed Central Ltd. 2012
Received: 8 June 2011
Accepted: 15 May 2012
Published: 30 May 2012
Animal models are frequently used to assess new treatment methods in cancer research. MRI offers a non-invasive in vivo monitoring of tumour tissue and thus allows longitudinal measurements of treatment effects, without the need for large cohorts of animals. Tumour size is an important biomarker of the disease development, but to our knowledge, MRI based size measurements have not yet been verified for small tumours (10−2–10−1 g). The aim of this study was to assess the accuracy of MRI based tumour size measurements of small tumours on mice.
2D and 3D T2-weighted RARE images of tumour bearing mice were acquired in vivo using a 7 T dedicated animal MR system. For the 3D images the acquired image resolution was varied. The images were exported to a PC workstation where the tumour mass was determined assuming a density of 1 g/cm3, using an in-house developed tool for segmentation and delineation. The resulting data were compared to the weight of the resected tumours after sacrifice of the animal using regression analysis.
Strong correlations were demonstrated between MRI- and necropsy determined masses. In general, 3D acquisition was not a prerequisite for high accuracy. However, it was slightly more accurate than 2D when small (<0.2 g) tumours were assessed for inter- and intraobserver variation. In 3D images, the voxel sizes could be increased from 1603 μm3 to 2403 μm3 without affecting the results significantly, thus reducing acquisition time substantially.
2D MRI was sufficient for accurate tumour size measurement, except for small tumours (<0.2 g) where 3D acquisition was necessary to reduce interobserver variation. Acquisition times between 15 and 50 minutes, depending on tumour size, were sufficient for accurate tumour volume measurement. Hence, it is possible to include further MR investigations of the tumour, such as tissue perfusion, diffusion or metabolic composition in the same MR session.
KeywordsAnimal model Magnetic resonance Volume determination Cancer
The efficacy of new tumour treatment methods are usually tested in animal models prior to clinical trials(e.g. [1–3]). The end-point studied is often the change in tumour size using mice with subcutaneous xenografts of human tumour tissue (e.g. ). Two or three perpendicular diameters of the tumour are measured and volume is calculated under assumptions of e.g. ellipsoidal tumour shape (e.g. ). These methods are fast, inexpensive, and do not require anaesthesia. Studies on orthotopic and metastatic tumour models have increased [6–9]. Generally, such studies exclude external volume measurements and require non-invasive imaging techniques, such as magnetic resonance imaging (MRI), computed tomography (CT), positron emission tomography or single photon emission computed tomography[10, 11].
MicroCT imaging allows accurate measurement of very small tumour volumes (~10−2 g) in short acquisition time , but the radiation exposure might affect tumour growth and thus confound the interpretation of the treatment results. MRI does not expose the animal to ionizing radiation and gives no known side effects which opens the possibility to follow the same animal by consecutive non-invasive measurements in longitudinal studies. In addition, MRI has intrinsic contrast abilities which can be used to assess tumour tissue characteristics, e.g. diffusion, perfusion and metabolic properties, parameters reflecting cell proliferation, apoptosis, vascularisation, etc. [11, 13, 14, 15, 16, 17, 18, 19, 20, 21].
Multiparametric MR imaging requires optimised MR methods for minimum time consumption. Small tumours require careful attention regarding imaging parameters due to the partial volume effect (PVE). To our knowledge, there is a lack of information on the accuracy and precision of MRI based size measurements on small tumours.
The aim of this study was 1) to develop and assess the accuracy of an MRI based method for in vivo measurement of tumour size in a small animal tumour model, and 2) to investigate the possibility to reduce acquisition time by reducing image resolution, without affecting the accuracy.
Tumour models and experimental setup
The study was performed on 17 tumour bearing female, athymic BALB/c and SCID mice s.c. inoculated in the neck/back region with the human midgut carcinoid cell line, GOT1 (n = 15), or the medullary thyroid carcinoma, GOT2 (n = 2) under anaesthesia using tribromoethanol (Avertin®, Winthrop Laboratories, Surbiton, UK). Animals had free access to water and standard diet.
Within two days after MR imaging, animals were sacrificed by surgical incision of the heart after i.p. injection of sodium pentobarbital (Apoteket Farmaci, Stockholm, Sweden). Gauge block measurements of the length, l, width, w and height, h of the tumours were performed. Assuming ellipsoidal tumours, the volume V = (π/6)·l·w·h, was determined and converted to mass, m GB , assuming a tumour density of 1.0 g/cm3. The tumour was then surgically removed and weighed for reference tumour mass, m T , determination. All procedures were approved by the local Animal ethics committee in Gothenburg.
MRI was performed using a 7 T MR system (Bruker BioSpin MRI GmbH, Ettlingen, Germany; software: ParaVision 5.0) equipped with water cooled gradient coils (maximum gradient strength 400 mT/m). A 72 mm volume coil was used for transmission and an actively decoupled 4 channel array rat brain coil was used for signal receiving (RAPID Biomedical GmbH, Rimpar, Germany). Proper animal positioning was assured using a fast gradient-echo localizer scan. All other imaging experiments were based on a respiration triggered and fat suppressed T2-weighted rapid acquisition with relaxation enhancement (RARE) sequence . A 2D set of axial images was acquired (TR = 3438–4186 ms, TEeff = 30 ms, NSA = 2–5, turbo factor Tf = 4, pixel size: ~150·200 μm2, slice thickness = 700 μm) prior to a 3D volume acquisition (TR = 3000 ms, TEeff = 35–82 ms, NSA = 1–2, Tf = 4–14). In the 3D acquisition, the matrix size was adjusted to achieve 1603 μm3 voxels with minimal field-of-view (FOV) but full tumour coverage. The 3D sequence was repeated twice using voxel sizes of 2003 μm3 and 2403 μm3 but with maintained FOV. The methods are henceforth designated 3D-160, 3D-200 and 3D-240, respectively. Saturation slices were always used in the phase encoding direction.
During the MRI experiments anaesthesia was maintained using a mixture of air and isoflurane (1.5–2.0 %) (Isoba vet., Schering-Plough Animal-health, Farum, Denmark). Body temperature was maintained with a heating pad on the animal bed and a pressure sensitive pad was used for respiratory triggering.
Image processing and tumour mass estimation
The operator determining the masses from the repeated scans with increasing voxel sizes (m 3D-160 , m 3D-200 , and m 3D-240 ) was blinded to images with higher resolution when such existed.
Intra- and interobserver variability
Repeated (n = 10) segmentation and mass calculation was performed to investigate the intraobserver variability. Five image series including a small tumour (m T = 0.01, 2D and 3D-160), a medium sized tumour (m T = 0.10 g, 3D-160) and a large tumour (m T = 0.87 g, 2D and 3D-160) were included in the assessment. The operator was blinded to the 3D images when evaluating the 2D images. Learning effects were minimized by presenting the image series in a random order and allowing at least several hours to pass between sessions. The coefficient of variation (CV) was calculated for each image series.
Two additional observers repeated the process on the five image series to investigate interobserver variations. Their results were compared to the corresponding average values obtained by observer 1 in the intraobserver assessment.
Estimation of influence of partial volume effect
A simulation was set up using MATLAB to quantify the maximum influence that PVE could have on volume estimations based on 3D MRI images of a spherically shaped tumour. Tumour diameter and voxel dimensions were used as input parameters. A 3D grid was used to simulate the image matrix, with grid elements acting as image voxels. A sphere, simulating the tumour, was superimposed on the grid with coinciding origins. Grid elements located inside and outside the sphere were defined as tumour and background elements, respectively. Elements intersected by the rim of the sphere were defined as PVE elements. Each one-element thick section of the grid (representing the image slices) was considered for two extreme cases of segmentation, where PVE elements were classified either as tumour elements, or as background elements, thus overestimating or underestimating the true tumour volume, respectively.
The total influence of PVE on the tumour volume estimation was defined as %, where represents the total volume of tumour elements in the simulation and V is the true tumour volume. Voxel sizes of 1003–2503 μm3 and 3–9 mm (diameter) tumours were simulated, and the resulting data were presented as relative volume errors for different number of voxels-per-tumour-diameter ratios.
Regression analyses were performed using the Microsoft Excel (2003) statistics toolkit. Pearson’s r squared correlation coefficient was used.
Determination of tumour mass
All tumours found at necropsy were well visualized in the MR images. The mass of the resected tumours, m T , was between 0.01 g and 2.28 g. Tumours with m T <0.2 g are henceforth classified as small.
Comparing m 2D to m T (n = 15, one image series was lost, and one did not cover the entire tumour) a strong correlation was obtained (R2 = 0.99) (Figure 3b). Also for small tumours a strong correlation was found (R2 = 0.96), i.e. the 2D method was comparable to the 3D-160 method regarding accuracy.
To exclude the inherent uncertainty of the digital balance, m 2D was compared directly to m 3D-160 . The linear relation in the regression and the correlation approached unity m 2D = 1.01m 3D-160 –0.00, R2 = 1.00, n = 15. For small tumours the corresponding relation was m 2D =1.04 m 3D-160 + 0.00, R2 = 0.99, n = 9.
The analysis of m 3D-160 , m 3D-200 , and m 3D-240 vs. m T did not reveal any obvious voxel size dependence in the range of tumour- and voxel sizes investigated. The acquisition time required for the smallest tumour studied (m T = 0.01 g) was clearly reduced with larger voxels: 55, 34 and 22 minutes for the 1603, 2003 and 2403 μm3 images, respectively.
Intra- and interobserver variability
The coefficient of variation (CV) calculated for the five intraobserver variability assessments. The CV is based on 10 volume calculations performed on each of the five image series
m T (g)
Interobserver variation. The interobserver variation of tumour mass measurements for three observers given as relative deviation (in per cent) from the average value of the mass (n = 10) obtained by observer 1 for three different tumours and two imaging methods studied
Interobserver variation (%)
Influence of partial volume effect on volume determination. The maximum possible influence of the partial volume effect (PVE) on tumour volume estimations, assuming isotropic voxels and a spherical tumour, determined by simulations
Maximum relative volume difference (%)
The present study evaluated the possibilities and requirements for accurate measurement of small tumour sizes in mice using MRI. A linear, highly correlated relationship between weighed and MRI measured tumour masses down to 10−2 g was found. 3D acquisition should be considered when tumour masses of 10−1 gram or less are expected, due to the relative increase of PVE. In very small tumours (10−2 g) image acquisition at high resolution (in our setup 1603 μm3 voxels) should also be considered. The increased acquisition time using high resolution is compensated for by the smaller FOV needed to cover the tumour. Short acquisition times allow either additional MR investigations on the same animal, such as determining tumour diffusion, perfusion or metabolic parameters within one MR session, or a higher animal throughput.
Initially, the 2D method was included only for anatomical reference since gradient performance limited the minimum slice thickness to 700 μm, resulting in substantial PVE. However, the accuracy of volume estimations based on the 2D images was similar to that based on the 3D images (Figures 2 and 3b). This might be due to the possibility to study adjacent image slices, which probably improves the view of the tumour shape and delineation of the tumour border. However, the 2D method has a higher interobserver variability (Table 2), which supports the use of a 3D method for very small tumours, in order to limit the subjectivity in the evaluation.
In the 3D method, the turbo factor (Tf) was adjusted to reduce acquisition time. An increase in Tf results in an increased point spread function (PSF). Computer simulations assuming similar acquisition parameters as those used, and T2 values common at 7 T  showed that the PSF was broadened only by a factor of 1.6 compared to the value for Tf = 1 (data not shown). The minimum TEeff is also affected by the Tf, i.e. the image contrast will vary slightly with Tf. However, tumours were always easily visualized and, altogether, the range of Tfs used in the study might only affect the results to a minor extent.
The most time consuming process in the volume determination was probably when adjusting the threshold value in images where the global segmentation had failed. To reduce the time of analysis one could e.g. calculate an average volume based on two extreme segmentations; one including most of the border, and one excluding it. Such a procedure would, however, overestimate the volume, especially for small voxel-to-diameter ratios, where asymmetry between over- and underestimated volume errors is more pronounced (Table 3). In situations when small tumours require polygon delineation of the tumour border, the decision to assign voxels intersected by the polygon line to the tumour or the background compartment, will require asymmetry consideration since it might have a significant effect on the volume estimation (Table 3).
The tumour density assumption (1.0 g/cm3) might be an underestimation that would account for the fact that the relations between the determined masses (m 3D-160 , m 2D and m GB ) and the weight (m T ) were less than unity. Another contributing factor could be the inherent uncertainty in the digital balance, since errors in the predictor used in regression analysis are typically manifested as a decrease of the slope coefficient towards zero . These two considerations are justified by the fact that the slope was 1.01 when m 3D-160 and m 2D were compared directly to each other, thus excluding the density effect and predictor uncertainties.
Few studies were found in the literature where the accuracy of MRI based tumour size measurement was verified by e.g. comparison with weight after resection. He et al. found a correlation of R2 = 0.96 (n = 7) when comparing volumes of pancreatic tumours in mice from T2-weighted 2D MR images (similar to our 2D method), acquired at 4.7 T, for 0.2–2.0 g tumours . We obtained a similar correlation (R2 = 0.96, n = 9) but for smaller tumours (0.01–0.2 g) since we used six times smaller voxel volume, i.e. the influence of PVE might be comparable. Other groups have reported MRI tumour size measurements in mice but without verification with tumour weight, e.g. [7–9]. One group reported MRI measurements of tumours in mouse pancreas down to 0.14 g at 7.0 T using sequence parameters similar to our 2D method (T2-weighted RARE sequence, 0.015 mm3 voxels), but verified the volume determination by one phantom measurement only .
Using high resolution microCT (7 min acquisition time, voxel size of 813 μm3, voxel volume of 0.0005 mm3) a close correlation was found for 0.02–0.25 g s.c. tumours in mice, verified by weight after resection (R2 = 0.97, n = 20) . Thus, microCT is a fast and accurate method, but the absorbed dose delivered to the animal, and especially to the tumour tissue, is a confounding factor in therapy response assessments.
The generally faster T1-weighted sequences have also been used for tumour imaging (e.g. ). Often however, the tumour and surrounding tissue have similar T1-values, thus requiring use of contrast agents. T1-weighted sequences without contrast agents may be useful for imaging tumours in e.g. the bladder wall, where fluid constitutes the surrounding tissue, due to widely different T1-values in fluid and solid tissue .
This study shows the feasibility of accurate and precise MRI based size measurement of small tumours in mice. Furthermore, short acquisition times that allow additional MR investigations within the same session can be realized by careful selection of MR sequence parameters, e.g. using 2D instead of 3D methods and optimal voxel size.
The authors thank Lilian Karlsson and Ann Wikström for excellent animal handling. This study was supported by grants from the European Commission FP7 Collaborative Project TARCC HEALTH-F2-2007-201962, the Swedish Research Council, the Swedish Cancer Society, BioCARE - a National Strategic Research Program at University of Gothenburg, the King Gustav V Jubilee Clinic Cancer Research Foundation. The work was performed within the EC COST Action BM0607.
- Schmitt A, Bernhardt P, Nilsson O, Ahlman H, Kolby L, Maecke HR, Forssell-Aronsson E: Radiation Therapy of Small Cell Lung Cancer with 177Lu-DOTA-Tyr3-Octreotate in an Animal Model. J Nucl Med. 2004, 45 (9): 1542-1548.PubMedGoogle Scholar
- Masciotti J, Provenzano F, Papa J, Klose A, Hur J, Gu X, Yamashiro D, Kandel J, Hielscher A: Monitoring Tumor Growth and Treatment in Small Animals with Magnetic Resonance and Optical Tomographic Imaging. Proceedings of SPIE: 21 January 2006; San Jose. Edited by: Fred S. 2006, Azar, Dimitris N, Metaxas, 608105-1-Google Scholar
- Kölby L, Bernhardt P, Johanson V, Schmitt A, Ahlman H, Forssell-Aronsson E, Macke H, Nilsson O: Successful receptor-mediated radiation therapy of xenografted human midgut carcinoid tumour. Br J Cancer. 2005, 93 (10): 1144-1151. 10.1038/sj.bjc.6602845.View ArticlePubMedPubMed CentralGoogle Scholar
- Tomayko M, Reynolds C: Determination of subcutaneous tumor size in athymic (nude) mice. Cancer chemotherapy and pharmacology. 1989, 24 (3): 148-154. 10.1007/BF00300234.View ArticlePubMedGoogle Scholar
- Euhus DM, Hudd C, Laregina MC, Johnson FE: Tumor measurement in the nude mouse. J Surg Oncol. 1986, 31 (4): 229-234. 10.1002/jso.2930310402.View ArticlePubMedGoogle Scholar
- He Z, Evelhoch JL, Mohammad RM, Adsay NV, Pettit GR, Vaitkevicius VK, Sarkar FH: Magnetic Resonance Imaging to Measure Therapeutic Response Using an Orthotopic Model of Human Pancreatic Cancer. Pancreas. 2000, 21 (1): 69-76. 10.1097/00006676-200007000-00054.View ArticlePubMedGoogle Scholar
- Mazurchuk R, Glaves D, Raghavan D: Magnetic resonance imaging of response to chemotherapy in orthotopic xenografts of human bladder cancer. Clin Cancer Res. 1997, 3 (9): 1635-1641.PubMedGoogle Scholar
- Grimm J, Potthast A, Wunder A, Moore A: Magnetic resonance imaging of the pancreas and pancreatic tumors in a mouse orthotopic model of human cancer. Int J Cancer. 2003, 106 (5): 806-811. 10.1002/ijc.11281.View ArticlePubMedGoogle Scholar
- Partecke I, Kaeding A, Sendler M, Albers N, Kuhn J-P, Speerforck S, Roese S, Seubert F, Diedrich S, Kuehn S, et al: In vivo imaging of pancreatic tumours and liver metastases using 7 Tesla MRI in a murine orthotopic pancreatic cancer model and a liver metastases model. BMC Cancer. 2011, 11 (1): 40-10.1186/1471-2407-11-40.View ArticlePubMedPubMed CentralGoogle Scholar
- Forssell-Aronsson E, Kjellén E, Mattsson S, Hellström M, Swedish Cancer Society Investigation Group T: Medical Imaging for Improved Tumour Characterization, Delineation and Treatment Verification. Acta Oncologica. 2002, 41 (7): 604-614. 10.1080/028418602321028201.View ArticlePubMedGoogle Scholar
- Weissleder R: Scaling down imaging: molecular mapping of cancer in mice. Nat Rev Cancer. 2002, 2 (1): 11-18. 10.1038/nrc701.View ArticlePubMedGoogle Scholar
- Jensen M, Jorgensen J, Binderup T, Kjaer A: Tumor volume in subcutaneous mouse xenografts measured by microCT is more accurate and reproducible than determined by 18 F-FDG-microPET or external caliper. BMC Medical Imaging. 2008, 8 (1): 16-10.1186/1471-2342-8-16.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorensen AG: Magnetic Resonance As a Cancer Imaging Biomarker. J Clin Oncol. 2006, 24 (20): 3274-3281. 10.1200/JCO.2006.06.6597.View ArticlePubMedGoogle Scholar
- De Schepper AM, De Beuckeleer L, Vandevenne J, Somville J: Magnetic resonance imaging of soft tissue tumors. Eur Radiol. 2000, 10 (2): 213-223. 10.1007/s003300050037.View ArticlePubMedGoogle Scholar
- Low RN: Abdominal MRI advances in the detection of liver tumours and characterisation. The Lancet Oncology. 2007, 8 (6): 525-535. 10.1016/S1470-2045(07)70170-5.View ArticlePubMedGoogle Scholar
- Hakumäki JM, Brindle KM: Techniques: Visualizing apoptosis using nuclear magnetic resonance. Trends Pharmacol Sci. 2003, 24 (3): 146-149. 10.1016/S0165-6147(03)00032-4.View ArticlePubMedGoogle Scholar
- Raman V, Pathak AP, Glunde K, Artemov D, Bhujwalla ZM: Magnetic resonance imaging and spectroscopy of transgenic models of cancer. NMR Biomed. 2007, 20 (3): 186-199. 10.1002/nbm.1136.View ArticlePubMedGoogle Scholar
- Kettunen MI, Brindle KM: Apoptosis detection using magnetic resonance imaging and spectroscopy. Prog Nucl Magn Reson Spectrosc. 2005, 47 (3–4): 175-185.View ArticleGoogle Scholar
- Barrett T, Brechbiel M, Bernardo M, Choyke PL: MRI of tumor angiogenesis. J Magn Reson Imaging. 2007, 26 (2): 235-249. 10.1002/jmri.20991.View ArticlePubMedGoogle Scholar
- Schmuecking M, Boltze C, Geyer H, Salz H, Schilling B, Wendt TG, Kloetzer K-H, Marx C: Dynamic MRI and CAD vs. Choline MRS: Where is the detection level for a lesion characterisation in prostate cancer?. Int J Radiat Biol. 2009, 85 (9): 814-824. 10.1080/09553000903090027.View ArticlePubMedGoogle Scholar
- Pacheco-Torres J, López-Larrubia P, Ballesteros P, Cerdán S: Imaging tumor hypoxia by magnetic resonance methods. NMR Biomed. 2011, 24 (1): 1-16. 10.1002/nbm.1558.View ArticlePubMedGoogle Scholar
- Hennig J, Nauerth A, Friedburg H: RARE imaging: A fast imaging method for clinical MR. Magn Reson Med. 1986, 3 (6): 823-833. 10.1002/mrm.1910030602.View ArticlePubMedGoogle Scholar
- Crémillieux Y, Ding S, Dunn JF: High-resolution in vivo measurements of transverse relaxation times in rats at 7 Tesla. Magn Reson Med. 1998, 39 (2): 285-290. 10.1002/mrm.1910390216.View ArticlePubMedGoogle Scholar
- Frost C, Thompson SG: Correcting for regression dilution bias: comparison of methods for a single predictor variable. Journal of the Royal Statistical Society: Series A (Statistics in Society). 2000, 163 (2): 173-189. 10.1111/1467-985X.00164.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/12/12/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.