Skip to main content

Predictors of pulmonary metastases on chest computed tomography in children and adolescents with osteosarcoma—tips for qualifying patients for thoracotomy

Abstract

Background

Osteosarcoma is the most common primary malignant bone tumour in children and adolescents. Lungs are the most frequent and often the only site of metastatic disease. The presence of pulmonary metastases is a significant unfavourable prognostic factor. Thoracotomy is strongly recommended in these patients, while computed tomography (CT) remains the gold imaging standard. The purpose of our study was to create tools for the CT-based qualification for thoracotomy in osteosarcoma patients in order to reduce the rate of useless thoracotomies.

Methods

Sixty-four osteosarcoma paediatric patients suspected of lung metastases on CT and their first-time thoracotomies (n = 100) were included in this retrospective analysis. All CT scans were analysed using a compartmental evaluation method based on the number and size of nodules. Calcification and location of lung lesions were also analysed. Inter-observer reliability between two experienced radiologists was assessed. The CT findings were then correlated with the histopathological results of thoracotomies. Various multivariate predictive models (logistic regression, classification tree and random forest) were built and predictors of lung metastases were identified.

Results

All applied models proved that calcified nodules on the preoperative CT scan best predict the presence of pulmonary metastases. The rating of the operated lung on the preoperative CT scan, dependent on the number and size of nodules, and the total number of nodules on this scan were also found to be important predictors. All three models achieved a relatively high sensitivity (72–92%), positive predictive value (81–90%) and accuracy (74–79%). The positive predictive value of each model was higher than of the qualification for thoracotomy performed at the time of treatment. Inter-observer reliability was at least substantial for qualitative variables and excellent for quantitative variables.

Conclusions

The multivariate models built and tested in our study may be useful in the qualification of osteosarcoma patients for metastasectomy through thoracotomy and may contribute to reducing the rate of unnecessary invasive procedures in the future.

Peer Review reports

Background

Osteosarcoma (OS) is the most common primary malignant bone tumour in children and adolescents [1]. OS is characterised by rapid hematogenous spread, with the lung being the most common site [2,3,4]. According to the literature, 10–25% of OS patients present with detectable metastases at the time of initial diagnosis, of which 85–90% have lung metastases [1, 5,6,7,8,9]. Tumour cells in OS metastases produce bone and this potential may be apparent on imaging. Moreover, the recurrence of OS is predominantly located in the lung (~ 80% of cases) [1]. For the assessment of pulmonary metastases, chest computed tomography (CT) has remained for years the gold imaging standard, as also suggested by the Children’s Oncology Group (COG) [6, 10,11,12]. With the continuous advances in multi-row-detector computed tomography (MDCT) scanners, the sensitivity in detecting small lung nodules has improved. Nonetheless, the distinction between malignant and benign pulmonary lesions on CT scans in paediatric patients with sarcomas is still below expectations [6, 11, 13,14,15], and the correct classification of pulmonary nodules on imaging as either malignant or benign is a clinical dilemma even for radiologists experienced in the field.

The presence of metastases has a significant impact on survival in OS patients [2, 6, 16,17,18]. Therefore, early diagnosis and appropriate treatment is an interdisciplinary challenge for the entire team involved, including the oncologist, surgeon, radiologist and pathologist. All OS metastases must be resected completely, regardless of their number and site, if the patient is treated with curative intent [5, 8, 19,20,21,22]. For pulmonary metastases, thoracotomy (TT) with manual exploration of the lung, which may improve survival, is the most strongly recommended technique and should be considered whenever feasible, even when repeated procedures are required [1, 5, 23,24,25,26].

On the other hand, TT is an additional burden for the oncologic patient. Moreover, a substantial percentage of paediatric patients with malignant solid primary tumours undergo invasive thoracic procedures that reveal only benign nodules [13, 14, 27]. We therefore assumed that a proper radiological qualification of OS patients for TT is highly significant. While both in the literature and in multi-institutional clinical trial protocols various criteria of lung metastases have been proposed [5, 13, 15, 27,28,29], to date no CT-based guidelines for qualifying patients with OS for TT have been established.

The purpose of this study was to construct multivariate models revealing tomographic predictors of lung metastases in children and adolescents with OS. An appropriate model will predict which TTs are justified (by confirming the presence of at least one pulmonary metastasis) and which are not. This in turn may contribute to reducing the rate of useless TTs.

Methods

This retrospective study was approved by the Bioethics Committee of the Institute of Mother and Child in Warsaw, Poland, with a waiver of parental and patient informed consent. The final study cohort was constituted by sixty-four high-grade OS paediatric patients, who were diagnosed and treated in accordance with the trial protocol of the European and American Osteosarcoma Study Group (EURAMOS) 1 [5]. They were selected from a larger group of patients with primary bone sarcomas (n = 112), after excluding patients with other bone sarcomas and patients with osteosarcoma who did not fulfil our inclusion criteria. All included osteosarcoma patients underwent their first-time metastasectomy through TT for presumed pulmonary metastases on chest CT between 2012 and 2018. Out of the total 139 first-time and subsequent TTs performed during the 7-year time period, only 100 first-time procedures (right- and/or left-sided) were included in the study, i.e., one or maximum two TTs per patient, in order to eliminate the influence of postoperative sequelae on CT evaluation. Figure 1 illustrates the study flow chart. In bilateral lesions, staged surgery was performed. Manual palpation of the lung was carried out through lateral TT and wedge resections were performed whenever possible. All lesions removed were studied histologically and the results are available in the archives of the Department of Pathomorphology. For each included TT, we reviewed two spiral chest CT scans queried on our picture archiving and communication system (PACS): the baseline CT scan (CT1), i.e., the first examination which revealed lung nodules on the eventually operated side, and the preoperative CT scan (CT2), i.e., the examination directly preceding the TT. The CT1 scan was performed before the start of treatment (for patients suspected of pulmonary metastases at diagnosis), during treatment (for patients suspected of disease progression), or during follow-up after treatment (for patients suspected of pulmonary recurrence). Nearly all CT scans were performed on MDCT scanners, of which most on Philips Brilliance 64 (64-detector row) with iDose reconstruction algorithms. Intravenous contrast was not administered routinely.

CT scans were retrospectively and independently reviewed by two radiologists with 15 years of experience at our institution (the national reference centre for children with bone and soft tissue tumours) who were blinded to the outcome of each TT. All available tools, including multiplanar reformatted reconstructions (MPR) and the maximum intensity projection (MIP) technique, were used on a dedicated workstation (Philips Extended Brilliance Workspace). Reader 1 reviewed CT scans for all included TTs (n = 100), and his findings were then correlated with the histopathological (HP) status of the TTs. Reader 2 reviewed CT scans for 50 randomly selected TTs, for which the inter-observer reliability between the two readers was tested. On all CT1 scans, both lungs together as well as the operated lung separately were assessed, while on the CT2 scans, the operated lung was assessed once again, each time using a compartmental evaluation method based on the number of nodules and their maximum axial diameter. For the combined assessment of both lungs, scoring criteria were adopted from the EURAMOS-1 protocol: “A”—at least one nodule of ≥ 10 mm or three or more nodules of ≥ 5 mm; “B”—scans not meeting score “A” criteria. We applied the following rating criteria for the operated lung: “a”—single nodule of > 10 mm or more than one nodule of > 5 mm; “b”—solitary nodule of 5–10 mm or multiple nodules of 3–5 mm; “c”—solitary nodule of < 5 mm or several nodules of < 3 mm. On both CT1 and CT2 scans, we also recorded the presence of calcified nodules, while on CT2 scans the total number of nodules, their maximum axial diameter and their location, which was determined using the depth from the nearest pleural surface (< 10 mm = peripheral including subpleural and pleural; ≥10 mm = central), were also noted.

HP status was presented on a three-point scale (0-1-2), where “0” was the presence of solely nonmetastatic nodules in the resected material, “1”—at most nonviable metastases (minimum 1 lesion of complete regression), and “2”—at least 1 viable metastasis.

An analysis of CT1 and CT2 characteristics was performed in order to find variables differentiating the two TT samples: with confirmed metastases (HP group “1 + 2”) and without metastases (HP group “0”). For qualitative variables, the two-sample z-test for proportions was used, while for quantitative variables, the Wilcoxon rank sum test was used. Three various multivariate predictive models were constructed (logistic regression, decision tree and random forest) and predictors of the HP status of TTs were identified. To examine which features were associated with TT outcome in the logistic regression, a multivariate model was created in which all predictors with univariate p values of < 0.20 were entered into the model. The model was reduced until only those predictors with p values of < 0.05 (considered statistically significant) remained. In the decision tree and random forest models predictors were selected by the algorithms used in these models. Each model was tested for sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) and accuracy, using leave-one-out cross-validation (for the logistic regression and decision tree models) and the bootstrap method (for the random forest model). The PPV of the models was compared with the PPV of the qualification for TT at the time of the treatment of the analysed patients. To assess the inter-observer reliability, Cohen’s kappa coefficient (κ) for qualitative variables and Lin’s concordance correlation coefficient (rc) for quantitative variables were used. Pearson’s correlation coefficient (r) was used to test the correlation between the quantitative variables. Statistical analyses were performed using the R software environment (version 4.0.5) under the GNU General Public License.

Fig. 1
figure 1

Study flow chart - osteosarcoma patients and their thoracotomies

Results

Patient and TT characteristics

Our cohort consisted of 64 patients with OS, 46 (72%) males and 18 (28%) females, with the primary tumour predominantly located in the peripheral skeleton: 34 (53%) in the femur, 15 (23%) in the tibia and 9 (14%) in the humerus; the mean age at the time of the biopsy of the primary tumour was 13.7 years (range: 5–19 years). Among their 100 first-time TTs included in this study, there were 72 staged bilateral metastasectomies (in 36 patients) and 28 unilateral metastasectomies, with 55 right-sided and 45 left-sided procedures. The mean time from the initial diagnosis to the TT was 12.9 months (range: 3–133 months). Table 1 presents the TT characteristics.

Table 1 Thoracotomy characteristics

CT1 scans

The combined rating of both lungs was “B” for 73 (73%) TTs, while the score of the operated lung only was “b” for 43 (43%) TTs. Calcified nodules were recorded on CT1 scans of 41 (41%) TTs. The mean time from the CT1 scan to the TT was 7.5 months (range: 1–33 months). Table 2 presents the characteristics of the CT1 scans.

Table 2 Characteristics of CT1 scans

CT2 scans

The total number of nodules detected and analysed on all CT2 scans was 630, while the mean number of nodules per TT was 6.3 (range: 1–32). The score of the operated lung was “b” for 51 (51%) TTs. Calcified nodules were found in 58 (58%) TTs; the mean number of calcified nodules was 1.3 (range: 0–10). Most lesions were located peripherally (mean 5.1, of which 2.2 sub-/pleurally). The mean time from the CT2 scan to the TT was 13.9 days (range: 0–84 days); the median was 7.5 days. Table 3 presents the characteristics of the CT2 scans.

Table 3 Characteristics of CT2 scans

HP outcome

The total number of nodules resected on all TTs was 422, while the mean number of nodules per TT was 4.2 (range: 1–19), of which 2.3 (54%) were nonmetastatic, 1.2 (29%) were viable metastases and 0.7 (17%) were nonviable metastases. The HP status of the TT was “2” for 54 (54%) TTs with at least 1 viable metastasis, “1” for 18 (18%) TTs with only nonviable metastases and “0” for 28 (28%) TTs confirming nonmetastatic nodules only. Among the latter, fibrosis, reactive changes, inflammation, atelectasis, intrapulmonary lymph nodes, congestion, haemorrhage, granulomas, normal lung parenchyma, focal pleural thickening and adhesions were the most common. Table 4 presents the HP outcome of the TTs.

Table 4 HP outcome of TTs

In the two-sample z-test for proportions and the Wilcoxon rank sum test variables differentiating at a significance level of 0.05, two TT samples (HP group “1 + 2” and HP group “0”) were identified (Table 5). They include: the CT1 rating of both lungs (p = 0.04168 for score “A”), the presence of calcified nodules on the CT1 scan (p = 0.04276), the CT2 rating of the operated lung (p = 0.00334 for score “a”; p = 0.00019 for score “c”), the presence of calcified nodules on CT2 scan (p = 0.00001), the mean number of calcified nodules on the CT2 scan (p = 0.00005), the CT1-CT2 rating of the operated lung (p = 0.03011 for score “a-a”; p = 0.00026 for score “c-c”) and the mean number of nodules of > 10.0 mm on the CT2 scan (p = 0.00381).

Table 5 Analysis of variables differentiating thoracotomies with and without metastases

Logistic regression model

Univariate analysis

The results of the univariate analysis of CT1 and CT2 variables as potential predictors of metastases at TT are provided in Table 6. The features of the CT1 scans, including score “A” for both lungs, score “a” for the operated lung and the presence of calcified nodules, as well as the features of the CT2 scans, such as score “a” or “b” for the operated lung, the number of calcified nodules and total number of nodules, were all associated with a higher risk of malignancy (HP status “1” or “2”). The odds ratio (OR) for score “a” of the operated lung on the CT2 scan was almost 40 (p = 0.000927).

The significance of the number of nodules of < 3.0 mm on the CT2 scan is due to the significance of the total number of nodules on this scan, because these variables show a strong positive correlation (p < 0.00001; r = 0.91).

Table 6 Potential predictors of metastases (univariate analysis)

Multivariate analysis

In the multivariate analysis that included all variables showing a univariate trend, the number of calcified nodules on the CT2 scan remained the only significant predictor of malignancy. With an increase in the number of these lesions by 1, the risk of metastases at TT (HP status “1” or “2”) increased more than threefold (Table 7). Using this model, it has been also found that the predicted probability of metastases in the operated lung was 93% for two calcified nodules on the CT2 scan, 98% for three such foci and almost 100% for four or more (Figs. 2 and 3).

Table 7 Predictors of metastases (multivariate analysis)
Fig. 2
figure 2

Predicted probability of lung metastases in patients with osteosarcoma (OS)

Fig. 3
figure 3

Late pulmonary recurrence of OS (11 years after the initial diagnosis). (a,b) Axial images from non-enhanced preoperative CT2 scan demonstrate two subpleural nodules with calcifications (arrows): (a) lung reconstruction; (b) mediastinal reconstruction with measurement from region of interest (ROI)—max. 469 Hounsfield units (HU). (c) HP slide of resected specimen at TT: haematoxylin and eosin (HE) staining; OS viable metastasis (conventional type) with partial regression after chemotherapy and central ossification (asterisk)

Decision tree model

For TTs with calcified nodules on the CT2 scan (n = 58, 58%), the predicted (90% probability) status was malignant (HP result “1” or “2”). These TTs are represented by the rightmost leaf of the decision tree illustrated in Fig. 4. Moreover, for TTs without calcified nodules but with a CT2 score for the operated lung of “a” or “b” and at least one central nodule on this scan (n = 14, 14%), the predicted (79% probability) status was malignant. For TTs without calcified nodules, with a CT2 score of “a” or “b”, but without central nodules (n = 11, 11%), the predicted status was benign (HP result “0”), although as many as 45% of TTs with these features were confirmed malignant. For TTs without calcified nodules and with a CT2 score of “c” (n = 17, 17%), the predicted status was benign with a probability of 76%, while the remaining 24% (4/17) of TTs with these features were verified malignant, but all of them had at most foci of complete regression, i.e., nonviable metastases (leftmost leaf of the decision tree, Fig. 4). Figure 5 illustrates the most important predictors of malignancy revealed by the decision tree model: the presence of at least one calcified nodule on the CT2 scan and the rating of the operated lung on this scan.

Fig. 4
figure 4

Schematic diagram of the decision tree model revealing predictors of HP status at TT

Fig. 5
figure 5

Pulmonary recurrence of OS 3.5 years after the initial diagnosis. (ac) Axial images from non-enhanced CT scans demonstrate a solitary peripheral nodule (arrow): (a) baseline CT1 scan rated “c”, lung reconstruction; (b,c) preoperative CT2 scan 8 months later shows moderate progression (score “b”) but prominent calcification of the nodule; (b) lung reconstruction; (c) mediastinal reconstruction with measurement from ROI—max. 647 HU. (d) HP slide of resected specimen at TT: HE staining; OS viable metastasis (conventional type) with ossification (asterisk)

Random forest model

In the random forest model, the number of calcified nodules on the CT2 scan and the total number of nodules on this scan were at the top of the ranking of predictors of TTs with pulmonary metastases (Figs. 6 and 7). The CT2 rating of the operated lung was also positioned high on the list.

Fig. 6
figure 6

Ranking of predictors of TTs with lung metastases obtained from the random forest model

Fig. 7
figure 7

Pulmonary OS metastases at initial diagnosis. (ac) Axial images from non-enhanced preoperative CT2 scan post-processed using the maximum intensity projection technique, displayed on split windows (left—lung reconstruction; right—mediastinal reconstruction), demonstrate multiple small nodules (arrows), with some arranged along vessels and some containing calcifications (ellipses). (d) HP slide of resected specimen at TT: HE staining; intravascular OS viable metastasis (conventional type) with ossifications (asterisks)

Evaluation of predictive models

All three multivariate predictive models proved to have a relatively high sensitivity (range: 72–92%), PPV (range: 81–90%) and accuracy (range: 74–79%). The PPV of each model was 9 to 18 percentage points higher than the PPV of the qualification for TT performed at the time of treatment (Table 8), when only 72% of TTs confirmed viable or nonviable metastases (Table 4).

Table 8 Evaluation of multivariate predictive models and comparison with qualification for TT at the time of treatment

Inter-observer reliability

Inter-reader agreement for the two radiologists reading the CT scans was almost perfect (κ > 0.8) or at least substantial (κ > 0.7) for the qualitative variables and excellent for all the quantitative variables (rc > 0.9).

Discussion

Although CT is so far the best imaging modality for the detection of lung nodules as part of disease staging, it is often not sufficient for differentiating benign from malignant lesions [6, 11, 13,14,15]. Many authors have already attempted to find predictors of malignancy at the level of either the nodule or the patient, often concluding that there are no definite imaging criteria of lung metastases [13,14,15, 27]. In patients with OS, TT with manual exploration of the whole lung is the most recommended metastasectomy technique, while minimally invasive surgery focused on the resection of selected nodules is strongly discouraged [5, 22]. Taking into account the above, and also the retrospective nature of our study, which did not guarantee a reliable correlation of nodules detected on CT and resected at TT, we have chosen an analysis at TT level, instead of a nodule-for-nodule comparison of CT characteristics for pathology. Nevertheless, in order to properly rate the lungs, all 630 nodules detected on the preoperative CT2 scans were first assessed in detail. Our inclusion criteria allowed for a homogeneous cohort, composed of 64 paediatric patients with OS suspected of lung metastases, who underwent a total of 100 TTs after two scans per TT performed with the same MDCT scanner in most of the cases.

Complete metastasectomy is the best predictor of survival in patients with metastatic OS [25, 28, 30,31,32]. Even small microscopic deposits of regressing metastases may recur if the patient does not have a TT with the removal of all residual disease [18]. We therefore considered viable and nonviable metastases (the latter formed by a calcified osteoid matrix or by a matrix and necrotic material without sarcoma cells [33]), which are distinguishable on pathology, but not on CT, together in a single HP group “1 + 2” for most of our analyses.

The presence of pulmonary nodules in sarcoma patients does not always represent metastatic disease, and there are still far too many TTs confirming benign lesions only. This problem has been already highlighted by many authors, among them: McCarville et al. [14], Picci et al. [27], Brader et al. [13], Kayton et al. [33] and Ciccarese et al. [28]. In our study, there were 28 (28%) useless TTs.

To address the aim of our study, we constructed three various multivariate predictive models—the logistic regression model used by most authors [13, 14, 29, 34, 35] and two more—the decision tree and random forest models. This is a strength of our report, because to the best of our knowledge the latter two models have not been previously used in a similar cohort. We only found one recent (2021) report in which three machine learning models were used in an attempt to predict lung metastases in patients with OS, but the radiomic features were extracted from CT images of the local tumour [36].

All of our multivariate predictive models proved that calcified nodules on the preoperative CT2 scan—either the number of these foci (logistic regression and random forest models) or the mere fact of their occurrence (decision tree model)—have the strongest ability to predict pulmonary metastases. Furthermore, in the logistic regression model, it was found that the predicted probability of metastases in the operated lung is very high for two or more calcified nodules on the CT2 scan, reaching nearly 100% for four such foci. The correlation between HP status (malignant versus benign) and the calcification in OS patients is unique and the opposite of what is observed in the general population [13, 37]. The explanation of this fact is that the osteoid matrix produced by OS cells may mineralise and become both apparent upon imaging [2, 38] and readily palpated at TT as a small “grain of sand” [18]. Our findings are in agreement with those described by Brader et al. [13] and Ciccarese et al. [28]. Interestingly, the EURAMOS-1 [5] and COG [12] definition of metastases does not mention the calcification criterion. In our cohort, calcified nodules were more frequent on the CT2 scans (58% of TTs) than on the CT1 scans (41% of TTs). The broad histomorphological spectrum of OS, dependent on the predominant matrix [2], presumably influences imaging appearance on the baseline scan, while the higher incidence of calcified nodules on the preoperative scan is assumed to be induced by chemotherapy often administered in between the scans. Although both the presence of calcified nodules on the CT1 scan and the number of calcified nodules on the CT2 scan showed a univariate trend, only the latter finally proved to be a significant predictor of malignancy in the multivariate logistic regression model.

Moreover, two of the three multivariate predictive models revealed another significant predictor—the CT2 rating of the operated lung, dependent on the number and size of nodules. This variable is positioned on the second level of the decision tree and relatively high in the ranking obtained from the random forest model. Score “a” (single nodule of > 10 mm or more than one nodule of > 5 mm) and score “b” (solitary nodule of 5–10 mm or multiple nodules of 3–5 mm) are associated with an increased probability of malignancy, as revealed by the decision tree model. Analysing previous reports, we found that regarding the individual nodule size, most authors proposed a threshold for metastases of 5–6 mm [13, 15, 27,28,29, 39]. However, Ciccarese et al. [28], Kusma et al. [29] and Lautz et al. [39] admitted that nodules < 5 mm in diameter still have a reasonable likelihood of malignancy and therefore the proposed values cannot be considered as a strict rule. In turn, the definition provided by the EURAMOS-1 trial [5] and COG [12] concerns the patient (not the individual nodule). Their authors defined “certain” pulmonary metastases at presentation as three or more lesions which are ≥ 5 mm in maximum diameter or as a single lesion ≥ 1 cm. Although it is not straightforward to compare three different scoring systems (focused on the nodule or the patient, or the operated lung), our findings seem to be similar to the results of the cited papers, while the scores “a” and “b” criteria are closest to the EURAMOS-1 and COG definition.

In the random forest model, the total number of nodules on the CT2 scan was positioned very high (second place) in the ranking of predictors. A similar finding was described by Picci et al. [27], who claimed that in their cohort all patients with more than seven nodules were metastatic, as opposed to 31% of patients with only one nodule.

The impact of central nodules on the CT2 scan (positioned at the lowest level of the decision tree model and in fifth place in the ranking of the random forest model) seems to be less significant than that of the predictors discussed above. There are divergent reports in the literature. According to most authors, lung metastases are usually peripherally located [30, 35, 40], but Ciccarese et al. [28] found that they were equally located in the central and in the subpleural region, while benign lesions had a predilection for the subpleural region.

Interestingly, the most significant predictors of pulmonary metastases (calcified nodules, rating of the operated lung and total number of nodules) are all CT2 variables, suggesting that the preoperative CT scan is more important for qualifying patients with OS for metastasectomy than the baseline scan.

All three multivariate predictive models achieved a relatively high sensitivity, PPV and accuracy. If they were applied to our cohort at the time of treatment, the PPV would have increased from 72% to 81–90%. Using these models in the future may therefore contribute to reducing the percentage of unnecessary invasive thoracic procedures in children and adolescents with OS. The selected predictors of pulmonary metastases may serve as helpful tips during patient qualification for TT, rather than absolute indications, because the final decision about metastasectomy should always be made by an interdisciplinary team, taking into account the sine qua non conditions of TT, i.e., a locally controlled primary tumour, the absence of extrapulmonary metastases, completely resectable lung metastases and sufficient postoperative lung function. Further research with the use of applications based on predictive models, computer-aided detection on CT workstations or artificial intelligence will support the development of guidelines and optimise the process of selecting candidates for pulmonary metastasectomy through thoracotomy.

There are some limitations to our study. The first is the retrospective design of our single-centre trial. Our patient cohort was not very large, but it did not differ much from those found in the literature concerning this rare disease. Furthermore, the trial inclusion criteria entering patients at various stages of disease and treatment, as well as the compartmental evaluation method used to rate the lung, limited the assessment of change in the number and size of lesions over time.

Conclusions

The three multivariate predictive models that were constructed and validated in this study revealed the most important predictors of lung metastases on the preoperative CT scan in children and adolescents with OS (Fig. 8), which can contribute to a decreased rate of useless thoracotomies in the future.

Fig. 8
figure 8

Venn diagram illustrating the most important predictors of pulmonary metastases on the preoperative CT scan. “a”—single nodule of > 10 mm or more than one nodule of > 5 mm; “b”—solitary nodule of 5–10 mm or multiple nodules of 3–5 mm

Data availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CI:

Confidence interval

COG:

Children’s Oncology Group

CT:

Computed tomography

CT1:

Baseline CT scan

CT2:

Preoperative CT scan

EURAMOS:

European and American Osteosarcoma Study Group

HE:

Haematoxylin and eosin

HU:

Hounsfield units

HP:

Histopathological

MDCT:

Multi-row-detector computed tomography

MIP:

Maximum intensity projection

MPR:

Multiplanar reformatted reconstructions

NPV:

Negative predictive value

OS:

Osteosarcoma

PACS:

Picture archiving and communication system

PPV:

Positive predictive value

OR:

Odds ratio

ROI:

Region of interest

SD:

Standard deviation

TT:

Thoracotomy

References

  1. Ingley KM, Maleddu A, Grange FL, Gerrand C, Bleyer A, Yasmin E, et al. Current approaches to management of bone sarcoma in adolescent and young adult patients. Pediatr Blood Cancer. 2022;69(2):e29442.

    Article  PubMed  Google Scholar 

  2. WHO Classification of Tumours. Soft Tissue and Bone Tumours. 5th ed. Lyon, France: IARC Press; 2020 2020.

  3. Kempf-Bielack B, Bielack SS, Jürgens H, Branscheid D, Berdel WE, Exner GU, et al. Osteosarcoma relapse after combined modality therapy: an analysis of unselected patients in the Cooperative Osteosarcoma Study Group (COSS). J Clin Oncol. 2005;23(3):559–68.

    Article  PubMed  Google Scholar 

  4. Chi SN, Conklin LS, Qin J, Meyers PA, Huvos AG, Healey JH, et al. The patterns of relapse in osteosarcoma: the Memorial Sloan-Kettering experience. Pediatr Blood Cancer. 2004;42(1):46–51.

    Article  PubMed  Google Scholar 

  5. EURAMOS-1. Clinical trial protocol. Version 4.0 2015 [updated 12.02.2015. https://www.ctu.mrc.ac.uk/euramos/euramos_i_trial.asp.

  6. Chiesa AM, Spinnato P, Miceli M, Facchini G. Radiologic Assessment of Osteosarcoma Lung Metastases: state of the art and recent advances. Cells. 2021;10(3).

  7. Bhattasali O, Vo AT, Roth M, Geller D, Randall RL, Gorlick R, et al. Variability in the reported management of pulmonary metastases in osteosarcoma. Cancer Med. 2015;4(4):523–31.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Costelloe CM, Chuang HH, Daw NC. PET/CT of Osteosarcoma and Ewing Sarcoma. Sem Roentgenol. 2017;52(4):255–68.

    Article  Google Scholar 

  9. Isakoff MS, Bielack SS, Meltzer P, Gorlick R. Osteosarcoma: current treatment and a collaborative pathway to Success. J Clin Oncol. 2015;33(27):3029–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Ritter J, Bielack SS, Osteosarcoma. Ann Oncol. 2010;21(Suppl 7):vii320–5.

    Article  PubMed  Google Scholar 

  11. Carrle D, Bielack S. Osteosarcoma lung metastases detection and principles of multimodal therapy. Cancer Treat Res. 2009;152:165–84.

    Article  PubMed  Google Scholar 

  12. Meyer JS, Nadel HR, Marina N, Womer RB, Brown KL, Eary JF, et al. Imaging guidelines for children with Ewing sarcoma and osteosarcoma: a report from the Children’s Oncology Group Bone Tumor Committee. Pediatr Blood Cancer. 2008;51(2):163–70.

    Article  PubMed  Google Scholar 

  13. Brader P, Abramson SJ, Price AP, Ishill NM, Emily ZC, Moskowitz CS, et al. Do characteristics of pulmonary nodules on computed tomography in children with known osteosarcoma help distinguish whether the nodules are malignant or benign? J Pediatr Surg. 2011;46(4):729–35.

    Article  PubMed  PubMed Central  Google Scholar 

  14. McCarville MB, Lederman HM, Santana VM, Daw NC, Shochat SJ, Li CS, et al. Distinguishing benign from malignant pulmonary nodules with helical chest CT in children with malignant solid tumors. Radiology. 2006;239(2):514–20.

    Article  PubMed  Google Scholar 

  15. Ghosh KM, Lee LH, Beckingsale TB, Gerrand CH, Rankin KS. Indeterminate nodules in osteosarcoma: what’s the follow-up? Br J Cancer. 2018;118(5):634–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ługowska I, Woźniak W, Klepacka T, Michalak E, Szamotulska K. A prognostic evaluation of vascular endothelial growth factor in children and young adults with osteosarcoma. Pediatr Blood Cancer. 2011;57(1):63–8.

    Article  PubMed  Google Scholar 

  17. Yahiro K, Matsumoto Y. Immunotherapy for osteosarcoma. Hum Vaccin Immunother. 2021;17(5):1294–5.

    Article  CAS  PubMed  Google Scholar 

  18. Picci P. Osteosarcoma (osteogenic sarcoma). Orphanet J Rare Dis. 2007;2:6.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Kager L, Zoubek A, Kastner U, Kempf-Bielack B, Potratz J, Kotz R, et al. Skip metastases in osteosarcoma: experience of the Cooperative Osteosarcoma Study Group. J Clin Oncol. 2006;24(10):1535–41.

    Article  PubMed  Google Scholar 

  20. Kager L, Zoubek A, Pötschger U, Kastner U, Flege S, Kempf-Bielack B, et al. Primary metastatic osteosarcoma: presentation and outcome of patients treated on neoadjuvant Cooperative Osteosarcoma Study Group protocols. J Clin Oncol. 2003;21(10):2011–8.

    Article  PubMed  Google Scholar 

  21. Mialou V, Philip T, Kalifa C, Perol D, Gentet J-C, Marec-Berard P, et al. Metastatic osteosarcoma at diagnosis. Cancer. 2005;104(5):1100–9.

    Article  PubMed  Google Scholar 

  22. Gao E, Li Y, Zhao W, Zhao T, Guo X, He W, et al. Necessity of thoracotomy in pulmonary metastasis of osteosarcoma. J Thorac Dis. 2019;11(8):3578–83.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ferrari S, Briccoli A, Mercuri M, Bertoni F, Picci P, Tienghi A, et al. Postrelapse survival in osteosarcoma of the extremities: prognostic factors for long-term survival. J Clin Oncol. 2003;21(4):710–5.

    Article  PubMed  Google Scholar 

  24. Bielack SS, Kempf-Bielack B, Branscheid D, Carrle D, Friedel G, Helmke K, et al. Second and subsequent recurrences of osteosarcoma: presentation, treatment, and outcomes of 249 consecutive cooperative osteosarcoma study group patients. J Clin Oncol. 2009;27(4):557–65.

    Article  PubMed  Google Scholar 

  25. Pastorino U, Palmerini E, Porcu L, Luksch R, Scanagatta P, Meazza C et al. Lung metastasectomy for osteosarcoma in children, adolescents, and young adults: proof of permanent cure. Tumori. 2021:3008916211053048.

  26. de Bree E, Drositis I, Michelakis D, Mavroudis D. Resection of Pulmonary metastases in Osteosarcoma. Is it justified? Hellenic J Surg. 2018;90(6):293–8.

    Article  Google Scholar 

  27. Picci P, Vanel D, Briccoli A, Talle K, Haakenaasen U, Malaguti C, et al. Computed tomography of pulmonary metastases from osteosarcoma: the less poor technique. A study of 51 patients with histological correlation. Ann Oncol. 2001;12(11):1601–4.

    Article  CAS  PubMed  Google Scholar 

  28. Ciccarese F, Bazzocchi A, Ciminari R, Righi A, Rocca M, Rimondi E, et al. The many faces of pulmonary metastases of osteosarcoma: retrospective study on 283 lesions submitted to surgery. Eur J Radiol. 2015;84(12):2679–85.

    Article  PubMed  Google Scholar 

  29. Kusma J, Young C, Yin H, Stanek JR, Yeager N, Aldrink JH. Pulmonary nodule size < 5 mm still warrants Investigation in patients with Osteosarcoma and Ewing Sarcoma. J Pediatr Hematol Oncol. 2017;39(3):184–7.

    Article  CAS  PubMed  Google Scholar 

  30. Ahmed G, Zamzam M, Kamel A, Ahmed S, Salama A, Zaki I, et al. Effect of timing of pulmonary metastasis occurrence on the outcome of metastasectomy in osteosarcoma patients. J Pediatr Surg. 2019;54(4):775–9.

    Article  PubMed  Google Scholar 

  31. Heaton TE, Hammond WJ, Farber BA, Pallos V, Meyers PA, Chou AJ, et al. A 20-year retrospective analysis of CT-based pre-operative identification of pulmonary metastases in patients with osteosarcoma: a single-center review. J Pediatr Surg. 2017;52(1):115–9.

    Article  PubMed  Google Scholar 

  32. Cheung FP, Alam NZ, Wright GM. The past, Present and Future of Pulmonary Metastasectomy: a review article. Ann Thorac Cardiovasc Surg. 2019;25(3):129–41.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kayton ML, Huvos AG, Casher J, Abramson SJ, Rosen NS, Wexler LH, et al. Computed tomographic scan of the chest underestimates the number of metastatic lesions in osteosarcoma. J Pediatr Surg. 2006;41(1):200–6. discussion – 6.

    Article  PubMed  Google Scholar 

  34. Murrell Z, Dasgupta R. What predicts the risk of recurrent lung metastases? J Pediatr Surg. 2013;48(5):1020–4.

    Article  PubMed  Google Scholar 

  35. Murrell Z, Dickie B, Dasgupta R. Lung nodules in pediatric oncology patients: a prediction rule for when to biopsy. J Pediatr Surg. 2011;46(5):833–7.

    Article  PubMed  Google Scholar 

  36. Pereira HM, Leite Duarte ME, Ribeiro Damasceno I, de Oliveira Moura Santos LA, Nogueira-Barbosa MH. Machine learning-based CT radiomics features for the prediction of pulmonary metastasis in osteosarcoma. Br J Radiol. 2021;94(1124):20201391.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Khan A. ACR appropriateness criteria on solitary pulmonary nodule. J Am Coll Radiol. 2007;4(3):152–5.

    Article  PubMed  Google Scholar 

  38. Oliveira GAd, Werlang HZ, Bergoli PM, Frechiani M, Oliveira F. Tomografia computadorizada na análise dos padrões de calcificações nos tumores ósseos da bacia em pediatria: nova abordagem. Radiologia Brasileira. 2006;39:413–8.

    Article  Google Scholar 

  39. Lautz TB, Farooqui Z, Jenkins T, Heaton TE, Doski JJ, Cooke-Barber J, et al. Thoracoscopy vs thoracotomy for the management of metastatic osteosarcoma: a Pediatric Surgical Oncology Research Collaborative Study. Int J Cancer. 2021;148(5):1164–71.

    Article  CAS  PubMed  Google Scholar 

  40. Briccoli A, Rocca M, Salone M, Guzzardella GA, Balladelli A, Bacci G. High grade osteosarcoma of the extremities metastatic to the lung: long-term results in 323 patients treated combining surgery and chemotherapy, 1985–2005. Surg Oncol. 2010;19(4):193–9.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors wish to thank Irena Koch from the Department of Pathomorphology at the Institute of Mother and Child for her assistance in collecting histopathological results.

Funding

This research received no external funding and was supported by the internal grant of the Institute of Mother and Child in Warsaw, Poland (No. 510-14-01).

Author information

Authors and Affiliations

Authors

Contributions

M.D. contributed to the conception and design of the work, acquisition, analysis, and interpretation of data, has drafted the work, and has substantively revised the work. A.D. contributed to the design of the work, acquisition and analysis of data, has drafted the work, and has substantively revised the work. A.O. contributed to the design of the work, analysis, and interpretation of data, and has substantively revised the work. E.M. and K.B. contributed to the acquisition of data and have substantively revised the work. T.K., M.R.-P. and A.R. contributed to the acquisition of data. M.B.-F. contributed to the conception of the work and has substantively revised the work. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Marek Duczkowski.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Bioethics Committee of the Institute of Mother and Child in Warsaw, Poland (protocol code: 43/2019, date of approval: 30 May 2019), which waived the parental and patient consent due to the retrospective nature of the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Duczkowski, M., Duczkowska, A., Olwert, A. et al. Predictors of pulmonary metastases on chest computed tomography in children and adolescents with osteosarcoma—tips for qualifying patients for thoracotomy. BMC Pediatr 24, 382 (2024). https://doi.org/10.1186/s12887-024-04858-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12887-024-04858-0

Keywords