- Research article
- Open Access
- Open Peer Review
A new relative tumor sizing method in epi-metaphyseal osteosarcoma
BMC Cancervolume 15, Article number: 284 (2015)
The goal of this study was to develop a new method for determining tumor size to predict prognosis with high performance in osteosarcoma.
This study was approved by the institutional review board. We retrospectively reviewed 41 magnetic resonance (MR) images at diagnosis and 57 MR images after neoadjuvant chemotherapy from 59 patients with non-metastatic, high-grade extremity osteosarcoma, who had undergone surgery between October 1994 and October 2009.
A new parameter of tumor axial ratio (TAR) was designed to normalize tumor size by dividing the absolute tumor axial size by the reference bone axial size (RBS) of the affected bone. RBS was defined using anatomical landmarks for each type of bone. Absolute tumor length (ATL), absolute tumor volume (ATV), and relative tumor volume (RTV) were comparatively analyzed. TAR was only significantly decreased after chemotherapy in the survival (P = 0.009) and metastasis-free (P = 0.018) group in the paired t-test. With the Kaplan–Meier method, significant differences in overall survival (log rank P = 0.004) and disease-free survival (Log Rank P = 0.009) were noted between decreased TAR after chemotherapy and increased TAR. After Cox regression analysis, TAR showed an odds ratios of 5.931 for survival (95% Confidence Interval [CI], 1.153–30.513) and 14.144 for metastasis (95% CI, 2.826–70.784), whereas ATL, ATV, and RTV showed no associations with these clinical variables. The AUC value of TAR was 0.713 (95% CI, 0.548 to 0.878) for survival and 0.759 (95% CI, 0.608 to 0.909) for metastasis.
TAR is a novel sizing method with potential as a prognostic tool in osteosarcoma.
The size of a primary tumor is the most fundamental parameter for staging solid tumors. In bone sarcomas, the American Joint Committee on Cancer (AJCC) 7th edition (2010) has established 8 cm in greatest dimension as the cutoff for subclassification of stage A or B . Besides AJCC staging, several studies have suggested various tumor sizing methods and cutoff points with prognostic powers in osteosarcoma [2-8].
Controversy over methods used to measure sizes of bone sarcomas has arisen from two major problems. First, there is concern over the performance of imaging devices. Prior to magnetic resonance imaging (MRI), the only two-dimensional (2D) sizing methods on roentgenological images did not show significant results . However, significant relationships between tumor burden (tumor volume) and prognosis have been reported for tumor volume based on MRI . Second, size heterogeneity of the affected organ is another major concern that is unique to bone sarcomas compared with other cancers. The sizes of different types of bones vary (e.g. humerus, femur, tibia, and fibula). Furthermore, the same type of bone varies in size depending on sex and age. Size heterogeneity is a serious concern because bone sarcomas mainly occur in children and teenagers. For this reason, the need for a better method to determine relative tumor size is important.
The mostly widely used method for adjusting tumor size is to normalize tumor burden according to individual patient size (body surface area [BSA]), as described by Bieling et al. . This method, which uses an ellipsoid formula to calculate tumor volume, was initially applied on plain X-ray images; however modified applications on MRI have also been reported [5,8,11]. In those studies, there were some inherent limitations that are important to consider. First, adjustment by BSA dose not discriminate bias from the different types of bones involved, although it can rule out bias from heterogeneity in individual patient size. Second, tumor volume based on ellipsoid formulas assumes that the shape of tumors is ellipsoidal, but this may not always be true. Advances in MRI have led to more correct measurement of tumor volume. Indeed, reports based on three-dimensional (3D) volumetry of tumors have shown prognostication of tumor burden [4,6]. However, calculating 3D volumetry on MRI is not clinically practical.
The purpose of this study was to develop a new method for determining tumor size with high performance for predicting prognosis. We developed a novel method for determining relative tumor size that focused on the relative axial length of a tumor.
We retrospectively reviewed 41 MR images at diagnosis and 57 after neoadjuvant chemotherapy from 59 patients with nonmetastatic, high-grade extremity osteosarcoma, who had undergone surgery between October 1994 and October 2009 and analyzed them together with other clinical data. Mean follow-up period was 114.7 months (range, 4.8–240 months). This study was done under a protocol approved by Severance Hospital Institutional Review Boardl. Both MRI at diagnosis and after neoadjuvant chemotherapy were available for 39 patients. MRI after neoadjuvant chemotherapy but not MRI at diagnosis was available for 18 patients. MRI at diagnosis but not MRI after neoadjuvant chemotherapy was available for two patients. The clinical characteristics of 59 patients are listed in Table 1. Sixty-five patients (92.9%) had received neoadjuvant chemotherapy. Patients were treated as follows: 22 received combination of intraarterial cisplatin and doxorubicin, while 37 received combination of intraarterial cisplatin, doxorubicin, and ifosfamide. Outcomes of neoadjuvant chemotherapy were not significantly different between doublet and triplet regimens in our cohorts .
Measurement of tumor size and novel parameter for tumor sizing
All parameters for tumor sizing were measured on MR images. Parameters for tumor sizing were defined as follows: absolute tumor length (ATL), the greatest longitudinal tumor length on either coronal or sagittal images; absolute tumor axial size (ATA), the greatest horizontal tumor length in any direction on axial images. Reference bone axial size (RBS) was defined using anatomical landmarks for each type of bone around the epi-metaphyseal area (Table 2). Measurement of all parameters was done with the Centricity Radiology RA1000 program (General Electrics Healthcare, United Kingdom). All parameters mentioned above were measured independently by three orthopaedic surgeons. Absolute tumor volume (ATV) was measured by one radiologist using 3D region of interest (ROI) magnetic resonance volumetry [13,14].
A new parameter for tumor sizing, tumor axial ratio (TAR), was designed to normalize tumor size by dividing ATA by RBS for each affected bone; thus, TAR represented the ratio of tumor axial size to affected bone axial size. The applications of our method are illustrated in Figure 1. For example, as shown in Figure 1a, distal femur with a RBS of 89.5 mm, and an ATA of 84.2 mm, the TAR would be 0.94 (84.2/89.5). Other applications for proximal tibia (Figure 1b), proximal humerus (Figure 1c), and proximal femur (Figure 1d) are also provided.
The reliability on measuring those parameters among different observers was evaluated using the interclass correlation coefficient. Changes in those parameters after neoadjuvant chemotherapy were analyzed by the paired t-test. Overall survival and disease-free survival between the decreased the TAR group after neoadjuvant chemotherapy and increased TAR were compared using the Kaplan–Meier method and log rank test. The correlations between each parameter and prognosis (5-year survival and metastasis) were expressed as the hazard ratio using Cox regression. Prognostication of TAR and ATV was compared by evaluating the area under the receiver operating characteristic curve (AUC).
All statistical analyses were performed using SPSS (version 20.0, SPSS, Inc., Chicago, IL, USA). All P values were two-tailed, and P < 0.05 was considered statistically significant.
New parameters and their reliability
To be accepted as a standard, the reliability of a method should be evaluated. The lengths of ATL, ATA, ASA, and RBS for all MR images enrolled in this study were measured independently by three orthopaedic surgeons and reliability analyses were performed (Table 3). Interclass coefficients (ICCs) for RBS were 0.947 on MR images at diagnosis and 0.965 on MR images after chemotherapy. These were convincing values because RBSs were defined from well-known and widely used anatomical landmarks. ICCs of ATA were 0.725 on MR images at diagnosis and 0.862 on MR images after neoadjuvant chemotherapy. Factors confounding measurement of ATA were peritumoral inflammatory changes and bone edema, which eventually resulted in the lowest ICCs on MR images at diagnosis. The ICC of ATA on MR images after neoadjuvant chemotherapy was increased as bone edema decreased by neoadjuvant chemotherapy. ICCs of TAR were 0.807 on MR images at diagnosis and 0.862 on MR images after neoadjuvant chemotherapy. Although ICCs of TAR were lower than those of ATL, they showed sufficient reliability.
Prognostications of parameters
Changes in parameters for tumor size after neoadjuvant chemotherapy were analyzed by the paired t-test (Table 4). Although there were no significant changes in ATL, ATV, and relative tumor volume (RTV) in all groups, ATA and TAR were significantly reduced in the survival group (P = 0.011 and P = 0.009, respectively) and metastasis-free group (P = 0.016 and P = 0.018, respectively). ATA and TAR in the mortality group and metastasis-positive group were not affected.
Overall survival and disease-free survival in the decreased TAR group after neoadjuvant chemotherapy and increased TAR were compared using the Kaplan–Meier method (Figure 2). Significant differences in overall survival (log rank P = 0.004) and disease-free survival (log rank P =0.009) were noted.
Prognostication of all parameters was analyzed with respect to 5-year survival and metastasis using Cox regression analysis (Table 5). None of the parameters on MR images at diagnosis showed associations with survival and metastasis. However, parameters for tumor axial extent on MR images after neoadjuvant chemotherapy showed associations with these prognostic variables, whereas the longitudinal extent and volume of tumors did not. ATA showed significant associations with prognosis (P < 0.001 for both survival and metastasis); but the odds ratios for ATA (1.070 [95% confidence interval (CI), 1.030–1.112] for survival and 1.073 [95% CI 1.035–1.112] for metastasis) revealed that they appeared to be useless parameters in practice. However, normalization with RBS made ATA strongly useful. TAR also showed significant associations with survival (P = 0.033) and metastasis (P = 0.001). TAR showed an odds ratio of 5.931 for survival (95% CI, 1.153–30.513) and 14.144 for metastasis (95% CI, 2.826–70.784). The AUC value for TAR was 0.713 (95% CI, 0.548–0.878) for survival and 0.759 (95% CI, 0.608–0.909) for metastasis, which was superior to that of ATV (0.588 for survival [95% CI, 0.397–0.778] and 0.609 for metastasis [95% CI, 0.431–0.788]) (Figure 3). Apart from associations with 5-years survival and metastasis, TAR also showed significant associations with histological response to chemotherapy, with an odds ratio of 10.746 (95% CI, 1.650–69.989) in Cox regression analysis, while the other parameters did not. In conclusion, TAR on MR images after neoadjuvant chemotherapy was the only parameter that predicted prognosis among parameters for tumor size.
In consideration of clinical application, we tried to assign cutoff values to TAR for predicting dichotomous outcomes. The cutoff value was 0.85 for survival and metastasis, which was derived from the Youden index. According to dichotomous predictions validated by two-way contingency table analysis, the accuracy, sensitivity, and specificity of TAR for survival were 75.81% (95% CI, 0.625–0.857), 73.91% (95% CI, 0.559–0.872), and 76.92% (95% CI, 0.663–0.848) respectively, while those for metastasis were 76.36% (95% CI, 0.624–0.865), 72.22% (95% CI, 0.509–0.877) and 78.38% (95% C,. 0.680–0.859).
The size of a primary tumor is regarded as a fundamental parameter for staging solid tumors, but there has been no proven standard method or cutoff value with high performance for predicting prognosis in bone sarcomas. Although many methods and cutoff values have been suggested, most studies have reported on the prognostication of their sizing method using a cutoff value optimized in their study populations. In those studies, there are still questions regarding generalization of the cutoff values. This study is believed to be the first to report a linear correlation (as an odds ratio) between parameter of tumor size and prognosis using successive data analysis in bone sarcomas.
Histological subtypes of osteosarcoma show distinct effects on survival and responses to chemotherapy. In this study, histological responses to chemotherapy were dependent on histological subtype (osteoblastic, chondroblastic, fibroblastic and mixed), validated by χ2 test (P < 0.000). However, there was no significant difference among the four subtypes for survival (P = 0.173) and metastasis (P = 0.261) in Cox regression analysis. We also evaluated relationships between parameters for tumor size and histological subtype. According to one-way analysis of variance, among the four subtypes, only longitudinal extent, ATL, was significantly different (P = 0.020 at diagnosis, P = 0.015 after neoadjuvant chemotherapy). Axial extent, TAR (P = 0.238 at diagnosis, P = 0.527 after neoadjuvant chemotherapy) and 3D ROI volume (P = 0.494 at diagnosis, P = 0.112 after neoadjuvant chemotherapy) were not.
The significance of the axial extent of a tumor has been considered. The Enneking staging system is also a concern for the axial extent of a tumor . Spanier et al. classified the amount of local axial extension of osteosarcoma as six grades. Multivariate analysis has shown that only the axial extent of a tumor has a significant effect on disease-free survival . Kim et al. reported that longitudinally growing tumors were associated with better survival than concentrically growing tumors in AJCC IIB osteosarcoma . On the basis of these findings, we focused on the axial extent of tumors. However, the application of a method for determining the axial extent of a tumor has limitations with respect to integration into clinical practice. The anatomical relationships between osteosarcoma and surrounding structures differ among affected bones, so the absolute value for invasion axial depth is not applied equally among other affected bones. Although the need for a method to determine relative tumor size has been raised, there is no established universal standard for normalizing tumor size in different bones. BSA is the most widely used method for normalizing tumor size. This approach may be suitable for normalizing tumor burden but not the axial invasion extent of tumors because it does not reflect axial size and spatial relationships between tumors and affected bones. Therefore, we defined new standards for normalizing tumor size from each affected bone by using well-known anatomical landmarks. This adjustment led to significant correlations with prognosis in our study. Absolute parameters for the axial extent of tumors did not show linear correlations with prognosis in this study, but relative parameters did show linear correlations with prognosis.
Many imaging modalities have been suggested to predict response to neoadjuvant chemotherapy and oncological outcomes. The maximum standardized uptake value (SUVmax) on 18 F-fluorodeoxyglucose positron emission tomography (18 F-FDG PET) has been reported as an indicator [18-20]. 3D volumetry on MRI with greater precisions compared with previous methods based on the ellipsoid formula has reported reduction of tumor volume as an indicator [4,6]. The apparent diffusion coefficient on diffusion-weighted MRI is also an indicator [21,22]. The major concerns with respect to these indicators are that the values measured after neoadjuvant chemotherapy are more convincing than the values measured at diagnosis. In the present study, TAR showed similar results. TAR measured on MR images at diagnosis was not prognostic, but it was prognostic when measured after neoadjuvant chemotherapy.
Our new method has several advantages over previous methods. First, tumor size is normalized according to the size of the affected bone rather than individual patient size. Normalization with BSA does not discriminate differences in the types of bones involved; however, our method is able to distinguish differences due to bone heterogeneity. Second, our method is easy to use. Calculation of tumor volume using ROI-based analysis on MR images is not practical clinically, but our method requires only measuring the lengths of tumors on the 2D plane of MR images. Third, our method focuses on the extent of tumor axial growth, which results in invasions into the surrounding tissue.
There were several limitations to our study. First, our study was a small retrospective study and only 39 of 59 patients had both MR images at diagnosis and after neoadjuvant chemotherapy, which may have led to wide confidence intervals for the associations of TAR with 5-year survival and metastasis, despite the statistically significant P values. This study needs to be validated in larger studies to assess applicability of TAR. Second, the application of our new method was restricted to osteosarcomas located around the epiphysis and metaphysis because RBS was defined by anatomical landmarks around these areas. Thus, our method cannot be applied to diaphyseal and periosteal osteosarcomas. It is also restricted to osteosarcoma located in large joints (shoulder, hip, and knee), although they comprise most cases of osteosarcoma. Third, further optimization of our method by more accurately defining RBS should be considered. Indeed, RBS of the proximal femur was defined as the epiphyseal plate of the femoral head in order to decrease RBS and increase TAR in light of poor prognosis, even though the trochanteric area is more frequently affected than the femoral head.
TAR is a novel relative sizing method with potential as a prognostic tool in osteosarcoma, which can discriminate differences caused by affected bone heterogeneity as well as individual patient size.
American Joint Committee on Cancer
Absolute tumor axial size
Absolute tumor length
Absolute tumor volume
Area under receiver operating characteristic curve
Body surface area
Reference bone axial size
Relative tumor volume
Tumor axial ratio
AJCC. Cancer Staging Manual. 7th ed. New York: Springer; 2010.
Wellings RM, Davies AM, Pynsent PB, Carter SR, Grimer RJ. The value of computed tomographic measurements in osteosarcoma as a predictor of response to adjuvant chemotherapy. Clin Radiol. 1994;49(1):19–23.
Bieling P, Rehan N, Winkler P, Helmke K, Maas R, Fuchs N, et al. Tumor size and prognosis in aggressively treated osteosarcoma. J Clin Oncol. 1996;14(3):848–58.
Shin KH, Moon SH, Suh JS, Yang WI. Tumor volume change as a predictor of chemotherapeutic response in osteosarcoma. Clin Orthop Relat Res. 2000;376:200–8.
Kaste SC, Liu T, Billups CA, Daw NC, Pratt CB, Meyer WH. Tumor size as a predictor of outcome in pediatric non-metastatic osteosarcoma of the extremity. Pediatr Blood Cancer. 2004;43(7):723–8.
Moon SH, Shin KH, Suh JS, Yang WI, Noh JK, Hahn SB. Tumor volume change after chemotheraphy as a predictive factor of disease free survival for osteosarcoma. Yonsei Med J. 2005;46(1):119–24.
Kim MS, Lee SY, Cho WH, Song WS, Koh JS, Lee JA, et al. An examination of the efficacy of the 8 cm maximal tumor diameter cutoff for the subdivision of AJCC stage II osteosarcoma patients. J Surg Oncol. 2008;98(6):427–31.
Lee JA, Kim MS, Kim DH, Lim JS, Yoo JY, Koh JS, et al. Relative tumor burden predicts metastasis-free survival in pediatric osteosarcoma. Pediatr Blood Cancer. 2008;50(2):195–200.
Davis AM, Bell RS, Goodwin PJ. Prognostic factors in osteosarcoma: a critical review. J Clin Oncol. 1994;12(2):423–31.
Bramer JA, van Linge JH, Grimer RJ, Scholten RJ. Prognostic factors in localized extremity osteosarcoma: a systematic review. Eur J Surg Oncol. 2009;35(10):1030–6.
Kim MS, Lee SY, Cho WH, Song WS, Koh JS, Lee JA, et al. Initial tumor size predicts histologic response and survival in localized osteosarcoma patients. J Surg Oncol. 2008;97(5):456–61.
Hong S, Shin SJ, Jung M, Jeong J, Lee YJ, Shin KH, et al. Comparison of long-term outcome between doublet and triplet neoadjuvant chemotherapy in non-metastatic osteosarcoma of the extremity. Oncology. 2011;80(1–2):107–17.
Yeo SG, Kim DY, Kim TH, Jung KH, Hong YS, Chang HJ, et al. Tumor volume reduction rate measured by magnetic resonance volumetry correlated with pathologic tumor response of preoperative chemoradiotherapy for rectal cancer. Int J Radiat Oncol Biol Phys. 2010;78(1):164–71.
Mayr NA, Taoka T, Yuh WT, Denning LM, Zhen WK, Paulino AC, et al. Method and timing of tumor volume measurement for outcome prediction in cervical cancer using magnetic resonance imaging. Int J Radiat Oncol Biol Phys. 2002;52(1):14–22.
Enneking WF, Spanier SS, Goodman MA. A system for the surgical staging of musculoskeletal sarcoma. Clin Orthop Relat Res. 1980;153:106–20.
Spanier SS, Shuster JJ, Vander Griend RA. The effect of local extent of the tumor on prognosis in osteosarcoma. J Bone Joint Surg Am. 1990;72(5):643–53.
Kim MS, Lee SY, Lee TR, Cho WH, Song WS, Cho SH, et al. Prognostic effect of pathologic fracture in localized osteosarcoma: a cohort/case controlled study at a single institute. J Surg Oncol. 2009;100(3):233–9.
Cheon GJ, Kim MS, Lee JA, Lee SY, Cho WH, Song WS, et al. Prediction model of chemotherapy response in osteosarcoma by 18 F-FDG PET and MRI. J Nucl Med. 2009;50(9):1435–40.
Costelloe CM, Macapinlac HA, Madewell JE, Fitzgerald NE, Mawlawi OR, Rohren EM, et al. 18 F-FDG PET/CT as an indicator of progression-free and overall survival in osteosarcoma. J Nucl Med. 2009;50(3):340–7.
Hongtao L, Hui Z, Bingshun W, Xiaojin W, Zhiyu W, Shuier Z, et al. 18 F-FDG positron emission tomography for the assessment of histological response to neoadjuvant chemotherapy in osteosarcomas: a meta-analysis. Surg Oncol. 2012;21(4):e165–170.
Bajpai J, Gamnagatti S, Kumar R, Sreenivas V, Sharma MC, Khan SA, et al. Role of MRI in osteosarcoma for evaluation and prediction of chemotherapy response: correlation with histological necrosis. Pediatr Radiol. 2011;41(4):441–50.
Neubauer H, Evangelista L, Hassold N, Winkler B, Schlegel PG, Kostler H, et al. Diffusion-weighted MRI for detection and differentiation of musculoskeletal tumorous and tumor-like lesions in pediatric patients. World J Pediatr. 2012;8(4):342–9.
The authors would like to thank all the patients enrolled in this study.
The authors declare that they have no competing interests.
SHK carried out the overall study design, data collection, data organization, data analysis/interpretation, writing of all drafts of the manuscript, and has approved final version of the submitted manuscript. KHS participated in study design, data collection, data organization, data analysis/interpretation, writing of all drafts of the manuscript, and has approved final version of the submitted manuscript. EHP participated in data analysis and has approved final version of the submitted manuscript. YJC participated in discussion about study design, data analysis/interpretation, and has approved final version of the submitted manuscript. BKP participated in data collection, data analysis/interpretation, and has approved final version of the submitted manuscript. JSS participated in data collection, data analysis/interpretation, and has approved final version of the submitted manuscript. WIY participated in data collection, data analysis/interpretation, and has approved final version of the submitted manuscript.