Long-term effects of local radiotherapy on growth and vertebral features in children with high-risk neuroblastoma

doi:10.21203/rs.3.rs-3734349/v1

Download PDF

Research Article

Long-term effects of local radiotherapy on growth and vertebral features in children with high-risk neuroblastoma

https://doi.org/10.21203/rs.3.rs-3734349/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background To evaluate the effects of local radiotherapy (RT) on growth, we evaluated the chronological growth profiles and vertebral features of children with high-risk neuroblastoma.

Methods Thirty-eight children who received local photon or proton beam therapy to the abdomen or retroperitoneumbetween January 2014 and September 2019 were included. Simple radiography of the thoracolumbar spine was performed before and every year after RT. The height and vertical length of the irradiated vertebral bodies (VBs) compared with the unirradiated VBs (vertebral body ratio, VBR) were analyzed using the linear mixed model. Shape feature analysis was performed to compare the irradiated and unirradiated vertebrae.

Results The follow-up was a median of 53.5 months (range, 21-81 months) after RT. A decline in height z-scores was mainly found in the early phase after treatment. In the linear mixed model with height, the initial height (fixed, p < 0.001), sex (time interaction, p = 0.008), endocrine dysfunction (time interaction, 0.019), and age at diagnosis (fixed and time interaction, both p = 0.002) were significant. Unlike the trend in height, the change in VBR (ΔVBR) decreased gradually (p < 0.001). The ΔVBR in the group that received more than 30 Gy decreased more than in the group that received smaller doses. In the shape feature analysis, the irradiated VBs changed to a more irregular surface that were neither round nor rectangular.

Conclusion The irradiated VBs in children were gradually restricted compared to the unirradiated VBs in long-term follow-up, and higher RT doses were significantly affected. Radiation-induced irregular features of VBs were observed.

Growth restriction

Radiotherapy

Proton beam therapy

Neuroblastoma

Neuroblastoma is the most common extracranial solid tumor in children [1, 2]. The median age at diagnosis is approximately two years, with most cases diagnosed before five years old [3]. It is an embryonic cancer that arises from neural crest stem cells, and the common primary sites are the adrenal medulla and paraspinal ganglia in the abdomen and pelvis [4]. Its clinical course is heterogeneous, with treatment based on risk. Treatment for high-risk neuroblastoma has recently progressed with high-dose chemotherapy, autologous stem cell transplantation (HDCT/auto-SCT), and external beam radiotherapy (RT) [2, 5]. The clinical outcomes have improved gradually, and the five-year overall survival rate has been reported to be over 70% [6]. However, toxicity after treatment is unavoidable [6]. Some children have experienced growth retardation due to endocrine dysfunction or the effect of RT. In our protocol, total body irradiation (TBI) has been replaced by metaiodobenzylguanidine (MIBG) therapy, with RT used as localized therapy for primary sites and/or metastatic lesions. Growth retardation has improved following the protocol change [7]. However, a decreasing trend in the z-score of height has been reported.

We aimed to evaluate the effect of chronological local RT on growth and spinal features on X-rays in children with neuroblastoma.

Patients

We retrospectively reviewed the medical records of children with high-risk neuroblastoma who were treated with local RT between January 2014 and September 2019. The high-risk group was defined as having at least one of the following risk factors: International Neuroblastoma Staging System (INSS) stage 4 with an age at diagnosis of one year or older and N-myc gene amplification regardless of stage [8]. Children were irradiated in the abdomen or retroperitoneum as a primary site, and simple radiography of the thoracolumbar (T-L) spine was taken before and every year after RT to check the growth and alignment of the spine. The exclusion criteria were: (1) boys over 12 years and girls over ten years and (2) the mediastinal origin of neuroblastoma irradiated out of the T-L spine. Ultimately, 38 children were included.

Treatment of neuroblastoma

The treatment for high-risk neuroblastoma was based on our protocol [8, 9]. Briefly, treatment consisted of nine cycles of induction chemotherapy with surgery after six cycles of chemotherapy, tandem HDCT/auto-SCT, local RT, differentiation therapy with 13-cis-retinoic acid, and immunotherapy with interleukin-2. Most patients in this study were treated using the up-to-date protocol SMC NB-2014 (NCT02771743), which modified the doses of the first HDCT and HD-MIBG according to the tumor status before each HDCT/auto-SCT based on ¹²³I-MIBG or ¹⁸FDG-PET uptake. In the NB-2009 protocol, TBI was not routinely performed because of toxicity and growth retardation [6].

Local RT was administered mainly to the primary tumor site approximately six weeks after the second HDCT/auto-SCT. The total dose in RT per fraction for patients with residual tumors was 15.0 Gy in 1.5 Gy per fraction for patients with residual tumors and 21.6 ~ 36.0 Gy in 1.8 Gy per fraction for patients with residual tumors before RT. The target volume was defined as the preoperative/postchemotherapy primary tumor bed, including residual tumor, if present, with a 1.0 to 1.5 cm margin. If a vertebral body (VB) was partially included in the target volume, the body was encompassed entirely up to at least 15 Gy because of concerns about inhomogeneous growth. From January 2016, our institute started proton beam therapy (PBT) using a posterior beam, which results in less exposure to the anterior part of the body. We previously used three-dimensional conformal radiotherapy (3D-CRT) with anterior-posterior parallel opposite fields.

Analysis of growth and vertebrae

The growth profiles of the children were regularly recorded at the time of diagnosis and at every visit. We analyzed height, weight, and body mass index (BMI) at the time of diagnosis (t_i), immediately prior to RT (pre-RT or t_r), and at annual follow-ups with T-L spine X-ray (t₁-t₅). The parameters are expressed as z-scores or percentiles of the same age based on the 2017 Korean National Growth Charts [10].

To evaluate the effect of RT on the spine, the mid-vertical lengths of irradiated VBs (iVBs) and unirradiated VBs (uVBs) were measured on the lateral view of the T-L spine X-ray. The length of iVB was measured for all the VBs in the RT target volume. The length of uVBs was measured at the same number of vertebrae above the RT field, excluding at least one vertebra directly above the RT field. The vertebral body ratio (VBR) at a certain time was defined as the length per iVBs divided by the length per uVB. The change in VBR (ΔVBR, %) at time t_i was defined as [(VBR at t_i -VBR at t₀) / VBR at t₀] × 100.

Shape feature analysis

To compare the sectional shapes of iVB and uVB, two-dimensional shape feature extraction was performed on anteroposterior (AP) and lateral views of simple X-rays to the T-L spine (Supplementary Fig. S1). For 14 patients who were followed up for five years after RT, one radiation oncologist contoured the VBs of L1 and T7, which were representatives of iVB and uVB, respectively, on the X-ray films taken before RT (pre-RT) and at the five-year follow-up (post-RT). MIM (version 6.4.9, MIM software, OH, US) was used for contouring. The films and contours were transferred to an in-house program as DICOM files. The following feature parameters were extracted: (1) area, perimeter, perimeter to the area, and sphericity from a contour; (2) convex area, convex perimeter, convexity, solidity, and elongation from a convex hull, which defines the smallest convex region, including a contour; and (3) eccentricity and rectangularity from a minimum bounding rectangle [11].

Statistical analysis

Linear mixed model analysis was performed on repeatedly measured data, such as height, weight, BMI, and VBR. The pre-RT and five-year follow-up shape features were compared using the paired t-test. Statistical analyses were performed using SPSS version 27 (SPSS Inc., IBM Company, Chicago, IL, USA). Statistical significance was set at p < 0.05. Bonferroni correction was applied for multiple comparisons.

Patient characteristics

The patient characteristics are summarized in the Supplementary Table S1. The median age at diagnosis was 36.5 months (range: 0-128 months). The study included 23 males (60.5%). Regarding the histologic type, most patients had neuroblastoma (84.2%), except six with ganglioneuroblastoma. Twenty-five (65.8%) patients underwent surgery on the primary tumors, and most patients were treated with MIBG (84.2%).

Local RT was performed at a median age of 50.5 months (13–142 months). The median interval from diagnosis to RT was 14 months (7–17 months). There was a median number of five iVBs (3–9), and the RT doses were 15, 25.2, 30.6, and 36.0 Gy in 17 (44.7%), five (13.2%), five (13.2%), and 11 (28.9%), respectively. PBT was performed in 25 patients (65.8%).

Events and complications related to growth

The median follow-ups were 66.5 (range, 38–96) and 53.5 (range, 21–81) months from diagnosis and RT, respectively. One child (2.6%) died from osteosarcoma, and six patients (15.8%) experienced disease progression (four distant metastases and two in-field recurrences with distant metastasis). Endocrine dysfunction, hypothyroidism, growth hormone deficiency (GHD), and sex hormone deficiency were observed in 16 (42.1%), 14 (36.8%), three (7.9%), and one patient (2.6%), respectively. Two children were diagnosed with hypothyroidism and GHD. All patients were treated with appropriate therapies, including hormone supplementation.

Abnormalities in spine arrangement were observed in three patients (two kyphotic and one scoliotic; Fig. 1). The features in two patients were already shown before RT (Figs. 1A and B and 1E and F). The patient with scoliosis (Fig. 1E and F) seemed to be progressive. None of the patients had any symptoms that required orthopedic treatment.

Growth profiles

Repeated measured z-scores for height, weight, and BMI are shown in Fig. 2, and the linear mixed model analysis results are shown in Table 1. Overall, the trends in the three parameters were significant (p < 0.001, p = 0.001, and p < 0.001, respectively). The height was initially 0.007 and decreased gradually. However, the decrease from t₁ to t₅ was not statistically significant. The z-score of the weight started at -0.320, decreased at t_r (p < 0.001), and then recovered at t_r and t₁.

Table 1

Z-scores of the height, weight, and BMI from diagnosis (timeline t_i), and a comparison between two time points of each year. (A) The analysis used the linear mixed model, and the p-values for subgroup comparison were adjusted using Bonferroni correction. (B) Factors related to the z-score of height by the linear mixed model
(A)	Timeline		t_i	t_r		t₁		t₂		t₃		t₄		t₅		p for time
Height (est. z-score)			0.007	-0.401		-0.647		-0.624		-0.788		-0.793		-0.777		< 0.001
Comparison (p-value)	T_i		-	< 0.001		< 0.001		< 0.001		< 0.001		< 0.001		0.001
	T_r		-	-		0.005		0.341		0.029		0.166		0.724
	T₁		-	-		-		1.000		1.000		1.000		1.000
	T₂		-	-		-		-		0.818		1.000		1.000
	T₃		-	-		-		-		-		1.000		1.000
	T₄		-	-		-		-		-		-		1.000
Weight (est. z-score)			-0.320	-0.694		-0.562		-0.443		-0.393		-0.376		-0.384		0.001
Comparison (p-value)	T_i		-	< 0.001		0.583		1.000		1.000		1.000		1.000
	T_r		-	-		1.000		0.480		0.684		1.000		1.000
	T₁		-	-		-		1.000		1.000		1.000		1.000
	T₂		-	-		-		-		1.000		1.000		1.000
	T₃		-	-		-		-		-		1.000		1.000
	T₄											-		1.000
BMI (est. z-score)			-0.604	-0.806		-0.394		-0.222		-0.108		-0.103		-0.161		< 0.001
Comparison (p-value)	T_i		-	1.000		1.000		0.197		0.075		0.186		0.872
	T_r		-	-		0.002		< 0.001		0.001		0.005		0.068
	T₁		-	-		-		1.000		1.000		1.000		1.000
	T₂		-	-		-		-		1.000		1.000		1.000
	T₃		-	-		-		-		-		1.000		1.000
	T₄		-	-		-		-		-		-		1.000
(B) Variables					Univariable (p-values)						Multivariable (p-values)
(B) Variables					Time		Variable		Interaction		Time		Variable		Interaction
Height at diagnosis (%)		< 30 vs. ≥30			< 0.001		< 0.001		0.155		< 0.001		< 0.001		0.119
Weight at diagnosis (%)		< 5 vs. ≥5			< 0.001		0.019		0.056		< 0.001		0.059		0.607
Age at diagnosis (months)		< 36 vs. ≥36			< 0.001		0.163		< 0.001		< 0.001		0.002		0.002
Sex		F vs. M			< 0.001		0.103		0.005		< 0.001		0.158		0.008
Surgery		Yes vs. No			< 0.001		0.158		0.907
MIBG therapy		Yes vs. No			< 0.001		0.269		0.666
No. of irradiated VB		3–4 vs. ≥5			< 0.001		0.567		0.024		< 0.001		0.943		0.274
RT dose, Gy		15.0-25.2 vs. 30.6–36.0			< 0.001		0.432		0.899
RT type		Photon vs. Proton			< 0.001		0.627		0.853
Endocrine dysfunction		Yes vs. No			< 0.001		0.422		0.004		< 0.001		0.187		0.019
Progression		Yes vs. No			0.001		0.857		0.304
BMI body mass index, MIBG metaiodobenzylguanidine, VB vertebral body, RT radiotherapy

Linear mixed analysis was performed to identify the significant factors related to height (Table 1). From the univariate analysis, the initial height and weight were significant when the time factor was fixed (p < 0.001 and p = 0.019, respectively). Regarding the time-variable interactions, age, sex, number of iVBs, and endocrine dysfunction were significant (p < 0.001, p = 0.005, 0.024, and 0.004, respectively). In the multivariable analysis of the factors with p < 0.1 from the univariable analysis, the initial height was associated with the z-score fixed by time (p < 0.001), and sex and endocrine dysfunction were significant in the time-variable interactions (p = 0.008, and 0.019, respectively). Age at diagnosis was related to the z-score trend in the time-fixed level and time-variable interactions (all p = 0.002).

Trends in the vertebral body ratio

Similar to the growth profiles, the VBR was also analyzed using a linear mixed model (Table 2). The VBR showed a significant declining trend (p < 0.001). The VBR started to decrease significantly from t₂ compared to the previous t_r or t₁ (all p-values < .005 from t₂) (Table 2). In the multivariable analysis based on the variables with p < 0.001 in the univariate analysis, only the RT dose was significant in both time-fixed level and time-variable interactions (p = 0.007 and 0.013). The ΔVBR in the group with 30.6 Gy or 36 Gy decreased more than that in the group with 15 Gy or 25.2 Gy (Fig. 3).

Table 2

Changes of the vertebral body ratio (Δ VBR, %) from pre-RT (timeline t_r) and the comparison between two time points of each year. (A) The analysis used the linear mixed model, and the p-values for subgroup comparison were adjusted using Bonferroni correction. (B) Factors related to ΔVBR by the linear mixed model
(A)	Timeline	t_r	t₁		t₂		t₃	t₄	t₅	p for time
ΔVBR (%)		0	-0.467		-2.203		-4.132	-5.944	-7.793	< 0.001
Comparison (p-value)	T_r	-	1.000		< 0.001		< 0.001	< 0.001	< 0.001
	T₁	-	-		1.000		< 0.001	< 0.001	< 0.001
	T₂	-	-		-		< 0.001	< 0.001	< 0.001
	T₃	-	-		-		-	0.008	< 0.001
	T₄	-	-		-		-	-	0.011
(B) Variables			Univariable (p-values)					Multivariable (p-values)
(B) Variables			Time	Variable		Interaction		Time	Variable	Interaction
Height at diagnosis (%)	< 30 vs. ≥30		< 0.001	0.584		0.947
Weight at diagnosis (%)	< 5 vs. ≥5		< 0.001	0.599		0.100		< 0.001	0.971	0.307
Age at diagnosis	< 36 vs. ≥36		< 0.001	0.548		0.427
Sex	F vs. M		< 0.001	0.164		0.247
Surgery	Yes vs. No		< 0.001	0.450		0.185
MIBG therapy	Yes vs. No		< 0.001	0.804		0.018		< 0.001	0.906	0.066
No. of VB	3–4 vs. ≥5		< 0.001	0.296		0.362
RT dose, Gy	15.0-25.2 vs. 30.6–36.0		< 0.001	0.006		0.006		< 0.001	0.007	0.013
RT type	Photon vs. Proton		< 0.001	0.583		0.371
Endocrine dysfunction	Yes vs. No		< 0.001	0.104		0.114
Progression	Yes vs. No		< 0.001	0.677		0.984
VB vertebral body, VBR vertebral body ratio, MIBG metaiodobenzylguanidine, RT radiotherapy

Change of shape feature

Fourteen patients underwent post-RT X-ray films at five years, and the changes in the shape feature parameters were calculated from the pre-RT X-ray (Table 3). The increase in area in L1 was less than that in T7 in both views (p = 0.002 and p = 0.011, respectively). The perimeter in L1 of the lateral view was not different from that in T7 (p = 0.399), which led to a smaller decrease in the perimeter of the area in L1 (p < 0.001). Sphericity decreased in L1 but increased in T7 (p = 0.023 and p < 0.001). The results imply that the body of L1 grew less but had a more irregular surface than T7. Regarding the parameters of the convex hull, the convex area, convex perimeter, and elongation in L1 increased less than those in T7. Convexity and solidity were not different in L1 and T1 in the AP view but were in the lateral view. Convexity and solidity in L1 in the lateral view decreased (p = 0.001 and 0.017, respectively), indicating that the surface of the contour was more irregular and less round. Regarding the minimum bounding rectangle parameters, most comparisons were significant, except for rectangularity in the AP view. After RT, the shape of L1 had a larger difference in the long and short axes than that of T7 and had changed away from the rectangular shape in the lateral view.

Table 3

Comparison of the five-year change after radiotherapy between irradiated and unirradiated vertebral bodies (L1 and T7) from the shape feature analysis
Shape	Parameter, Δ (%)	AP			LAT
Shape	Parameter, Δ (%)	L1	T7	p-value	L1	T7	p-value
(1) Contour	Area	88.4 ± 31.4	113.2 ± 37.7	0.002	102.7 ± 40.9	125.5 ± 48.3	0.011
	Perimeter	37.4 ± 10.7	43.5 ± 11.7	0.021	46.5 ± 15.7	48.9 ± 14.8	0.399
	Perimeter to area	-26.1 ± 6.4	-31.8 ± 5.8	0.001	-26.5 ± 7.2	-32.8 ± 6.3	< 0.001
	Sphericity¹⁾	-0.4 ± 2.3	1.5 ± 2.5	0.023	-3.2 ± 2.5	0.4 ± 2.6	< 0.001
(2) Convex hull	Convex area	90.0 ± 32.3	114.8 ± 36.8	0.023	104.8 ± 42.0	125.1 ± 47.2	0.047
	Convex perimeter	36.5 ± 11.0	41.8 ± 11.0	0.001	44.2 ± 15.5	49.5 ± 14.6	0.022
	Convexity²⁾	-0.6 ± 1.2	-1.1 ± 1.9	0.457	-1.6 ± 1.5	0.4 ± 1.4	0.001
	Solidity³⁾	-0.8 ± 1.4	-0.8 ± 1.4	0.944	-1.0 ± 1.5	0.1 ± 1.0	0.017
	Elongation⁴⁾	3.1 ± 6.6	9.9 ± 9.3	0.022	-4.0 ± 5.6	4.2 ± 10.1	0.013
(3) Minimum bounding rectangle	Eccentricity⁵⁾	2.5 ± 3.9	11.2 ± 5.2	< 0.001	-3.7 ± 3.9	2.9 ± 6.6	0.009
(3) Minimum bounding rectangle	Rectangularity⁶⁾	-0.9 ± 4.3	-1.0 ± 3.6	0.930	-1.4 ± 3.2	2.4 ± 2.7	0.001
¹⁾Sphericity = perimeter of a round having the same area as the contour/perimeter of the contour; ²⁾Convexity = perimeter of the convex hull/perimeter of contour; ³⁾Solidity = area of the contour/area of the convex hull; ⁴⁾Elongation = the length of the major axis of the convex hull/the length of the minor axis of the convex hull; ⁵⁾Eccentricity = the length of the minor axis of the minimum bounding rectangle/the length of the major axis of the minimum bounding rectangle; and ⁶⁾Rectangularity = the area of the contour/the area of the minimum bounding rectangle
AP anteroposterior view, LAT lateral view

We chronologically analyzed the growth of prepubertal children after treatment with neuroblastoma, including local RT. Similar to previous studies [7, 12], the z-scores of height decreased sharply from the diagnosis, and the slope became smooth from the second year. The z-score of weight decreased less than that of height and recovered slowly after treatment. This means that while the disease and treatment affect the growth restriction temporarily, the loss of height cannot be recovered despite appropriate nutritional status. The growth of pediatric patients with cancer has been reported to be affected by hormonal and non-hormonal factors [13]. Endocrine problems frequently occur in 40–60% of pediatric patients with cancer after treatment, and the main endocrine disturbances are disorders of the hypothalamic-pituitary axis, including GHD and thyroid dysfunction [14]. Chemotherapy and cranial irradiation were the main treatment options. These hormonal problems were found to be significant for growth retardation despite the absence of TBI or cranial RT and adequate supplements. Non-hormonal factors include age at the time of treatment, parental height, physical activity, nutrition, and damage to the growth plate from surgery or RT [13]. Similar to the results of this study, initial growth status and age at diagnosis were also significant factors in the z-score height.

In association with RT, the irradiated volume, total dose, or fraction size can affect growth [15]. The number of VBs and TBI were associated with a lower height percentile in our previous study [12]. However, the number of VBs was significant only in the univariate analysis. The estimated VBR change over five years was − 7.793%. For the five-year follow-up after RT in 14 children, the vertical length per iVB and uVB grew from 21.4 mm to 25.3 mm and 16.5 mm to 21.4 mm, respectively. The restricted length of iVB was estimated to be approximately 2.5 mm per VB. It is thought that a loss of growth in a median of five vertebrae would have a limited effect on overall height.

The higher the RT dose, the greater its effect on spinal growth [15–18]. The effect of the RT dose varied according to age at RT and can be critical for growth in children aged zero to two years, even at less than 10 Gy. In children aged two to six years, it showed a substantial effect at > 15 Gy. For prepubertal patients aged six years or older, the effect of RT might cause damage at 35 Gy or higher, and the negative effect may be less but possible at between 15 and 35 Gy. Although the approximate limit of the RT dose is known, the cut-off value may be unclear. Therefore, the European Society for Pediatric Oncology radiotherapy working group recommended a vertebral dose gradient within 5 Gy for prepubertal children aged two years or older, avoiding a dose of more than 20 Gy to seven or more thoracic vertebrae if possible [15]. At our institute, RT targets for paravertebral sites in prepubertal patients include all VBs with at least 15 Gy. We observed that the growth restriction of iVB at more than 30 Gy was more evident than at relatively low doses. Therefore, although the intensity of vertebral damage might differ between patients, bony growth restriction in the RT field seems impossible to avoid after local RT with ≥ 15 Gy. Therefore, we should be careful in deciding the RT dose and irradiated volume.

PBT as an RT modality was not significantly associated with height restriction. PBT has been preferred for pediatric patients with cancer because it has shown potential benefits, such as lower secondary malignancy [19] and better hematologic outcomes after RT [20]. PBT can spare the anterior part of the body when locally irradiated using posterior beams. However, little is known about the effect of PBT compared to X-ray therapy in relation to height restriction. A study that followed patients for a median of 13.9 months after PBT reported a decrease in the growth rate of approximately 23% at 15 Gy and 60% at 30 Gy [21]. However, the effect of PBT was difficult to ascertain in this study because only growth at one point over one to two years was observed, and no patients underwent X-ray therapy. In our patients who underwent PBT or photon therapy and were followed up for a longer time, PBT was found to have a similar effect on height or vertebral growth as X-ray therapy. Nevertheless, a longer follow-up period is required.

A previous study by our institute reported on the shape change of iVBs in the axial view on magnetic resonance imaging [12]. The signal change of the irradiated vertebrae in the T1 weighted image and the roundness in the T2 weighted image decreased at a median of 86 months after RT. In this study, using simple radiography, we analyzed changes in the vertical length and contours of the AP and lateral views, which are thought to be more related to actual growth than the axial view. In addition to the decrease in the vertical length, the surface irregularity of iVB increased compared to uVB. iVB showed atypical features (Fig. 1G). In the lateral view of the radiograph, severe irregularities at the anterior borders and convexity were often observed. Growth plates, as secondary ossification centers, are located at the upper and lower ends of the VB, which is near the cartilage of discs [22]. After the growth plate of the VB was damaged from treatment, the surface appeared to have permanent scars as ossification, although the volume of iVB increased gradually. Finally, iVB exhibited an irregular and bulging shape at the top and bottom. It can be suggested that the primary ossification center in the middle of the VB is less sensitive to RT than the superior and inferior growth plates as secondary ossification centers [15]. Researchers have suggested a direct effect on chondroblasts [23] and microvascular damage [24]. Blood supply might be better in the middle of the bone than at the edges, which might be an advantage for the primary ossification center. Because of the primary center, the radiation dose effect may be complicated but sometimes obvious at higher doses [25].

This study had several limitations. Because this was a retrospective, single-center study with a small sample, detailed statistical verification of the factors was challenging. A simple X-ray film was not obtained every year for some patients, and several missing values were adjusted using a linear mixed model. Although vertical length included loss of length from the VB and abnormal arrangement of the spine [15], spinal deformities were not analyzed due to their rarity. If more comprehensive data are obtained, factors related to RT could be verified for growth restriction and RT-related orthopedic diseases.

In children with neuroblastoma, local RT minimally affected height growth. However, iVB gradually grew significantly less than uVB, particularly at higher RT doses. We also observed highly irregular features of the vertebrae on the follow-up X-rays after RT.

AP anteroposterior

BMI Body mass index

GH Growth hormone

GHD Growth hormone deficiency

IMPT intensity-modulated proton therapy

IMRT intensity-modulated radiotherapy

INSS International Neuroblastoma Staging System

LAT lateral view; PBT: Proton beam therapy

PT Proton therapy; RT: Radiotherapy

SCT Stem cell transplantation

TBI Total body irradiation

VB Vertebral body

VBR Vertebral body ratio

Acknowledgements

Not applicable.

Authors’ contributions

DHL participated in the design of the study. KY participated in the data analysis, and the drafting of the manuscript. JHA and SYW participated in the statistical analysis. SHJ participated in the data collection and analysis. KWS and JWL participated in data collection. All authors read and approved the final manuscript.

Funding

Not applicable.

Data availability

The datasets generated and/or analyzed during the current study are not publicly available due to concerns related to patient confidentiality and privacy. However, they are available from the corresponding author upon reasonable request.

Ethics approval and consent to participate

This study was conducted in accordance with the Declaration of Helsinki for experiments involving human subjects. Also, this study was approved by the Institutional Review Board of the Samsung Medical Center (No. 2023-01-087), and the requirement for written informed consent was waived due to the study’s retrospective nature.

Consent for publication

Not Applicable.

Competing interests

The authors declare no competing interests

Kim S, Chung DH. Pediatric solid malignancies: neuroblastoma and Wilms' tumor. Surg Clin North Am. 2006;86:469–87.
Matthay KK, Reynolds CP, Seeger RC, Shimada H, Adkins ES, Haas-Kogan D, et al. Long-term results for children with high-risk neuroblastoma treated on a randomized trial of myeloablative therapy followed by 13-cis-retinoic acid: a children's oncology group study. J Clin Oncol. 2009;27:1007–13.
Davis S, Rogers MA, Pendergrass TW. The incidence and epidemiologic characteristics of neuroblastoma in the United States. Am J Epidemiol. 1987;126:1063–74.
Newman EA, Abdessalam S, Aldrink JH, Austin M, Heaton TE, Bruny J, et al. Update on neuroblastoma. J Pediatr Surg. 2019;54:383–9.
Castleberry RP, Kun LE, Shuster JJ, Altshuler G, Smith IE, Nitschke R, et al. Radiotherapy improves the outlook for patients older than 1 year with Pediatric Oncology Group stage C neuroblastoma. J Clin Oncol. 1991;9:789–95.
Sung KW, Son MH, Lee SH, Yoo KH, Koo HH, Kim JY, et al. Tandem high-dose chemotherapy and autologous stem cell transplantation in patients with high-risk neuroblastoma: results of SMC NB-2004 study. Bone Marrow Transplant. 2013;48:68–73.
Lee JW, Lee S, Cho HW, Ma Y, Yoo KH, Sung KW, et al. Incorporation of high-dose (131)I-metaiodobenzylguanidine treatment into tandem high-dose chemotherapy and autologous stem cell transplantation for high-risk neuroblastoma: results of the SMC NB-2009 study. J Hematol Oncol. 2017;10:108.
Yi ES, Son MH, Hyun JK, Cho HW, Ju HY, Lee JW, et al. Predictors of survival in patients with high-risk neuroblastoma who failed tandem high-dose chemotherapy and autologous stem cell transplantation. Pediatr Blood Cancer. 2020;67:e28066.
Lee JW, Son MH, Cho HW, Ma YE, Yoo KH, Sung KW, et al. Clinical significance of MYCN amplification in patients with high-risk neuroblastoma. Pediatr Blood Cancer. 2018;65:e27257.
Kim JH, Yun S, Hwang SS, Shim JO, Chae HW, Lee YJ, et al. The 2017 Korean National Growth Charts for children and adolescents: development, improvement, and prospects. Korean J Pediatr. 2018;61:135–49.
Mingqiang Y, Kidiyo K, Joseph R. A survey of shape feature extraction techniques. Pattern Recognit. 2008;15:43–90.
Yu JI, Lim DH, Jung SH, Sung KW, Yoo SY, Nam H. The effects of radiation therapy on height and spine MRI characteristics in children with neuroblastoma. Radiother Oncol. 2015;114:384–8.
Diller L, Chow EJ, Gurney JG, Hudson MM, Kadin-Lottick NS, Kawashima TI, et al. Chronic disease in the Childhood Cancer Survivor Study cohort: a review of published findings. J Clin Oncol. 2009;27:2339–55.
Thomas-Teinturier C, Salenave S. [Endocrine sequelae after treatment of pediatric cancer: From childhood to adulthood]. Bull Cancer. 2015;102:612–21.
Hoeben BA, Carrie C, Timmermann B, Mandeville HC, Gandola L, Dieckmann K, et al. Management of vertebral radiotherapy dose in paediatric patients with cancer: consensus recommendations from the SIOPE radiotherapy working group. Lancet Oncol. 2019;20:e155–e66.
Probert JC, Parker BR. The effects of radiation therapy on bone growth. Radiology. 1975;114:155–62.
Vaeth JM, Levitt SH, Jones MD, Holtfreter C. Effects of radiation therapy in survivors of Wilms's tumor. Radiology. 1962;79:560–8.
Probert JC, Parker BR, Kaplan HS. Growth retardation in children after megavoltage irradiation of the spine. Cancer. 1973;32:634–9.
Chung CS, Yock TI, Nelson K, Xu Y, Keating NL, Tarbell NJ. Incidence of second malignancies among patients treated with proton versus photon radiation. Int J Radiat Oncol Biol Phys. 2013;87:46–52.
Yoo GS, Yu JI, Cho S, Han Y, Oh Y, Lim DH, et al. Chronological Analysis of Acute Hematological Outcomes after Proton and Photon Beam Craniospinal Irradiation in Pediatric Brain Tumors. Cancer Res Treat. 2022;54:907–16.
Baba K, Mizumoto M, Oshiro Y, Shimizu S, Nakamura M, Hiroshima Y, et al. An Analysis of Vertebral Body Growth after Proton Beam Therapy for Pediatric Cancer. Cancers (Basel). 2021;13:349.
Zhang Y, Lenart BA, Lee JK, Chen D, Shi P, Ren J, et al. Histological features of endplates of the mammalian spine: from mice to men. Spine (Phila Pa 1976). 2014;39:E312–7.
Kember NF. Cell survival and radiation damage in growth cartilage. Br J Radiol. 1967;40:496–505.
Eifel PJ, Donaldson SS, Thomas PR. Response of growing bone to irradiation: a proposed late effects scoring system. Int J Radiat Oncol Biol Phys. 1995;31:1301–7.
Rao AD, Ladra M, Dunn E, Kumar R, Rao SS, Sehgal S, et al. A Road Map for Important Centers of Growth in the Pediatric Skeleton to Consider During Radiation Therapy and Associated Clinical Correlates of Radiation-Induced Growth Toxicity. Int J Radiat Oncol Biol Phys. 2019;103:669–79.

No competing interests reported.

Supplementarymaterials1.docx
Supplementary Information Table S1: Patient characteristics (N = 38). Fig. S1: A sample of the feature extraction. Contours of T7 (orange color) and L1 (blue color) vertebral bodies on the anteroposterior (A) and lateral (B) views of thoracolumbar X-ray films.

Download PDF

Editorial decision: Revision requested
05 Apr, 2024
Reviews received at journal
04 Apr, 2024
Reviews received at journal
02 Apr, 2024
Reviewers agreed at journal
10 Mar, 2024
Reviewers agreed at journal
09 Mar, 2024
Reviewers agreed at journal
23 Jan, 2024
Reviewers agreed at journal
05 Jan, 2024
Reviewers invited by journal
26 Dec, 2023
Editor assigned by journal
26 Dec, 2023
Editor invited by journal
26 Dec, 2023
Submission checks completed at journal
26 Dec, 2023
First submitted to journal
10 Dec, 2023

You are reading this latest preprint version

Long-term effects of local radiotherapy on growth and vertebral features in children with high-risk neuroblastoma

Status:

Version 1

Abstract

Figures

Background

Methods

Patients

Treatment of neuroblastoma

Analysis of growth and vertebrae

Shape feature analysis

Statistical analysis

Results

Patient characteristics

Events and complications related to growth

Growth profiles

Trends in the vertebral body ratio

Change of shape feature

Discussion

Conclusion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1