ORIGINAL RESEARCH article

Front. Oncol., 10 September 2020

Sec. Hematologic Malignancies

Volume 10 - 2020 | https://doi.org/10.3389/fonc.2020.01673

Prognostic Nomogram for Childhood Acute Lymphoblastic Leukemia: A Comprehensive Analysis of 673 Patients

  • 1. The Affiliated Hospital of Southwest Jiaotong University, Chengdu, China

  • 2. Department of Hematology, Peking University Cancer Hospital and Institute, Beijing, China

  • 3. The Center of Gastrointestinal and Minimally Invasive Surgery, The Third People’s Hospital of Chengdu, Chengdu, China

  • 4. Key Laboratory of Carcinogenesis and Translational Research, Ministry of Education, Peking University Cancer Hospital and Institute, Beijing, China

  • 5. State Key Laboratory of Molecular Oncology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Beijing, China

  • 6. Medical Research Center, The Third People’s Hospital of Chengdu, The Affiliated Hospital of Southwest Jiaotong University, The Second Chengdu Hospital Affiliated to Chongqing Medical University, Chengdu, China

Abstract

Objective:

Despite that the survival rate in childhood acute lymphoblastic leukemia (cALL) is excellent, subsets of high-risk patients with cALL still have high relapse rates, and the cure rate is well below that for which we should aim. The present study aims to construct a prognostic nomogram to better inform clinical practitioners and improve risk stratification for clinical trials.

Methods:

The developed nomogram was based on the therapeutically applicable research to generate effective treatment (TARGET) database. With this database, we obtained 673 cALL patients with complete clinical information. We identified and integrated significant prognostic factors to build the nomogram model by univariate and multivariate Cox analysis. The predictive accuracy and discriminative ability of the nomogram were determined by the concordance index (C-index), calibration curve, and area under the receiver operating characteristic (ROC) curve (AUC) of ROC analysis. Internal validations were assessed by the bootstrapping validation.

Results:

In the multivariate analysis of the primary cohort, the independent factors for survival were ETV6 RUNX1 fusion status, karyotype, minimal residual disease (MRD) at day 29, and DNA index, which were all integrated into the nomogram. The calibration curve for the probability of survival showed good agreement between the prediction by the nomogram and the actual observation. The C-index of the nomogram for predicting survival was 0.754 (95% CI, 0.715–0.793), and the AUCs for 3-, 5-, and 7-year survival were 0.775, 0.776, and 0.772, respectively.

Conclusion:

We comprehensively evaluated the risk of clinical factors associated with prognosis and carried out risk stratification. The nomogram proposed in this study objectively and accurately predicted the prognosis of children with ALL.

Introduction

Acute lymphoblastic leukemia (ALL) is the most common cancer in children and represents approximately one quarter of all cancers among persons younger than 15 years (1). The cure rate of ALL is increasing (2), but approximately 15–25% of patients will relapse after recovery, which is the leading cause of death in childhood acute lymphoblastic leukemia (cALL) patients (3, 4). Therefore, it is particularly important to determine the factors that affect the prognosis of cALL.

Many reported factors are related to the prognosis of cALL, such as ETV6 RUNX1 fusion status (5), mixed-lineage leukemia (MLL) gene rearrangement (6), TCF3 PBX1 fusion status (7), BCR-ABL1 fusion status (8), trisomy of leukemic cell chromosomes 4 and 10 (TRISOMY 4 10 status), or all usual combinations of trisomies (chromosomes 4, 10, 17, and 18) (9, 10), minimal residual disease (MRD) status in bone marrow of day 29 (11), the percent of blasts in bone marrow aspirate at day 29 of induction therapy (BMA blasts day 29) (12), cALL with Down syndrome (DS) (5), hypodiploid (DNA index <0.8) (13), and white blood cell (WBC) count at diagnosis (1416). However, no study has used the whole series of factors as a prognostic model to predict the prognosis of cALL.

Therapeutically Applicable Research to Generate Effective Treatment (TARGET) is a dynamically updated database of the National Cancer Institute (NCI)’s Office of Cancer Genomics (OCG), whose mission is to advance the molecular understanding of cancers with the goal of improving patient outcomes. cALL is one of the projects in the TARGET program, which includes phase I (B-ALL), phase II (B-ALL, T-ALL), and phase III (ALL). We used a complete sample of all relevant clinical features in phases 1 and 2, which included 1,842 ALL patients.

Currently, nomograms have been developed for the majority of cancer types. The use of nomograms has compared favorably to the traditional staging systems for many cancers (1719), and thus, they have been proposed as an alternative or even as a new standard (2022). To our knowledge, this study is the first attempt to establish a prognostic nomogram for cALL based on 1,842 cALL samples in the TARGET database.

Materials and Methods

Patients

Clinical parameters associated with childhood ALL patients up to June 10, 2019 were downloaded from the NCI TARGET database1. A total of 1,842 patients were included in the first and second phases of the data set. Of these subjects, patients with cALL who had information of main observation indicators were target subjects of this study, which include survival time, survival status, age at diagnosis, gender, ETV6 RUNX1 fusion status, TRISOMY 4 10 status, MLL status, TCF3 PBX1 fusion status, karyotype, BCR-ABL1 fusion status, central nervous system (CNS) status at diagnosis, BMA blasts day 8, BMA blasts day 29, cell of origin, Down syndrome status, MRD at day 29, and DNA index. The specific screening process is presented in Figure 1. The clinicopathological characteristics of patients in the cohorts are listed in Table 1. All source data are presented in Supplementary Table S1.

FIGURE 1

TABLE 1

CharacteristicsTarget-ALL (N = 673)
AliveDead
Survival status, n (%)553 (82.2)120 (17.8)
Age (years)
Median5.5
Range1.0–15.0
Gender, n (%)
Male293 (53.0)66 (55.0)
Female260 (47.0)54 (45.0)
ETV6 RUNX1 Fusion Status, n (%)
Positive106 (19.2)4 (5.0)
Negative447 (80.8)116 (95.0)
CNS Status, n (%)
CNS1452 (81.7)89 (74.2)
CNS277 (13.9)28 (23.3)
CNS324 (4.4)3 (2.5)
TRISOMY 4 10 Status, n (%)
Positive92 (16.6)4 (3.3)
Negative461 (83.4)116 (96.7)
Karyotype, n (%)
No trisomies in 4, 10, 17, 18367 (66.3)108 (89.9)
4, 10, 17, 18 have only one trisomy43 (7.8)5 (4.2)
Double trisomies28 (5.1)2 (1.7)
Triple trisomies115 (20.8)5 (4.2)
MLL Status, n (%)
Positive23 (4.2)9 (7.5)
Negative530 (95.8)111 (92.5)
TCF3 PBX1 Status, n (%)
Positive38 (6.9)15 (12.5)
Negative515 (93.1)105 (87.5)
WBC at Diagnosis (×103/mcL)
Median27
Range0.7–1306
BCR-ABL1 Status, n (%)
Negative543 (98.2)115 (95.8)
Positive10 (1.8)5 (4.2)
BMA Blasts Day 8 (%)
Median9
Range0–97
BMA Blasts Day 29 (%)
Median0
Range0–50
MRD of 29 day (%)
Median0
Range0–43
Cell of Origin, n (%)
B Cell ALL500 (90.4)115 (95.9)
T Cell ALL53 (9.6)5 (4.1)
Down Syndrome, n (%)
Yes5 (0.9)5 (4.2)
No548 (99.1)115 (95.8)
DNA Index
Median1
Range0–1.922

The clinicopathological characteristics of patients in the cohort.

WBC at Diagnosis, WBC count at diagnosis; CNS Status, CNS status at diagnosis.

Diagnosis

Age refers to the age at diagnosis. The WBC count at diagnosis was the absolute peripheral WBC count (in ×103/mcL). CNS status at diagnosis was determined according to the status of CNS leukemia (CNSL) at the time of diagnosis. Diagnosis and typing were according to Pediatric Acute Lymphoblastic Leukemia, Version 2.2020, NCCN Clinical Practice Guidelines in Oncology (23). MRD status was determined by flow cytometry in bone marrow on day 29. BMA blasts day 29 and day 8 represent the percent of blasts in bone marrow aspirate at day 29 or day 8 of induction therapy. The presence or absence of ETV6 RUNX1 fusion status, MLL status, TCF3 PBX1 status, and BCR-ABL1 fusion status were detected by fluorescence in situ hybridization (FISH), PCR, or cytogenetics. The type of cell origin is determined by bone marrow examination. Karyotype analysis was determined by chromosomal banding technique to find abnormal chromosome number of leukemia cells and structural changes such as translocation, inversion, and deletion (23).

Grouping of Clinical Characteristics

Age was divided into two groups based on the cutoff value of 10 years old (24, 25), and WBC count at diagnosis was divided into three groups: “<50,” “50 to 100,” and “≥100” (3, 26, 27) according to the WBC count (in ×103/mcL). Next, CNS status at diagnosis was divided into three groups: CNS1, CNS2, and CNS3. MRD day 29, which means minimum residual disease status at day 29 of induction therapy, was divided into four groups: <0.01, 0.01–0.1%, 0.1–1%, and >1% (28, 29). BMA blasts day 29, which represents the percent of blasts in bone marrow aspirate at day 29 of induction therapy, was divided into two groups with a 5% cutoff. Besides, BMA blasts day 8 was divided into two groups according to whether it is >20%. The karyotypes are divided into four groups according to the trisomy of chromosomes 4, 10, 17, and 18, which include “no trisomies in 4, 10, 17, 18,” “4, 10, 17, 18 have only one trisomy,” “double trisomies (DT),” and “triple trisomies (TT).” Cell of origin contained two subtypes, which included B cell ALL and T cell ALL. Finally, the DNA index, whose number represents the ratio of the DNA content or chromosome number in a tumor sample compared to a normal diploid sample, was divided into two groups: ≤0.8 and >0.8. ETV6 RUNX1 fusion status, MLL status, TCF3 PBX1 status, and BCR-ABL1 fusion status were divided into positive and negative according to whether its corresponding fusion gene were positive or negative.

Statistical Analyses

All data including demographic and disease characteristics were expressed as count (%). Statistical analysis was performed using the R software (Version 3.6.1)2.

The prognostic value of the 16 clinical characteristics was first calculated in the univariate Cox analysis; clinical features with a P < 0.1 in the multivariate Cox regression analysis were used to construct the nomogram via the “rms,” “survival,” and “foreign” packages of R (R version 3.6.1). Hazard ratio (HR) and 95% confidence interval (CI) were calculated. The performance of the nomogram was measured by the C-index and assessed by comparing the nomogram-predicted estimates versus the observed Kaplan–Meier estimates of survival probability (R package “rms”) (30). Based on the regression coefficients of the multivariate Cox regression analysis, a risk score composed of six clinical features in the nomogram was calculated, and the patients were divided into two groups by taking the corresponding median risk score as the cutoff point.

Kaplan–Meier curves and the log-rank test were used to compare the survival outcomes of the two groups with the R packages “survminer” and “survival” (R version 3.6.1). Receiver operating characteristic (ROC) curve analysis was employed to compare prediction concerning the accuracy and precision with the R package “survivalROC” (R version 3.6.1). Bootstrapping validation (1,000 bootstrap resamples) was used to calculate a relatively corrected C-index of the nomogram (31, 32). A P < 0.05 was considered significant.

Results

Clinicopathological Characteristics of the Patients

The cohort included 673 cALL patients with complete clinical information of main observation indicators. The median of age at diagnosis was 5.5 years (range, 1.0–15.0 years). Three-hundred fifty-nine (53.3%) patients in the cohort were male, while 314 patients were female. The median of survival time was 3,016 days (range, 28–5,598 days). Survival status showed 120 deaths and 553 survivals.

Independent Prognostic Factors in the Primary Cohort

In the primary cohort, we performed a univariate Cox regression analysis for each clinical factor (Table 2) and screened factors with P < 0.1, which included age, WBC count at diagnosis, CNS status at diagnosis, ETV6 RUNX1 fusion status, TRISOMY 4 10 status, karyotype, TCF3 PBX1 fusion status, BCR-ABL1 fusion status, MRD day 29, BMA blasts day 8, Down syndrome, and DNA index. Then, these factors were analyzed in a multivariate Cox regression model (Table 3). The results show that the patient and disease variables significantly associated with survival in multivariable modeling included positive ETV6 RUNX1 fusion (HR = 0.293, 95% CI = 0.124–0.695, P = 0.005), triple trisomies (TT) (HR = 0.211, 95% CI = 0.049–0.907, P = 0.036), MRD day 29 with 0.01–0.1% (HR = 1.942, 95% CI = 1.113–3.330, P = 0.016), MRD day 29 with 0.1–1% (HR = 2.140, 95% CI = 1.262–3.628, P = 0.005), MRD day 29 > 1% (HR = 2.508, 95% CI = 1.345–4.675, P = 0.004), and DNA index ≤0.8 (HR = 3.617, 95% CI = 1.233–10.612, P = 0.019).

TABLE 2

Clinical characteristicsHRLower 95% CIUpper 95% CIP-value
Age
≥10 years vs. <10 years1.5571.0712.2640.021*
Gender
Male vs. Female1.0760.7491.5440.693
ETV6 RUNX1 Fusion Status
Positive vs. Negative0.2330.1020.5290.0005***
TRISOMY 4 10 Status
Positive vs. Negative0.1900.0700.5150.001**
MLL Status
Positive vs. Negative1.7010.8593.3700.128
WBC (×103/mcL)
<50
50 to 1001.8061.1472.8450.011*
>1002.4201.4713.981<0.001***
Karyotype
No trisomies in 4, 10, 17, 18
Only trisomies 4 or 10 or 17 or 180.4640.1891.1370.093
Double trisomies (DT)0.2880.0711.1670.081
Triple trisomies (TT)0.1720.0700.4220.0001***
CNS Status
CNS1
CNS21.7391.1372.6620.011*
CNS30.6160.1951.9530.411
TCF3 PBX1 Status
Positive vs. Negative1.6920.9792.9240.060
BCR ABL1 Status
Positive vs. Negative2.2240.9085.4480.080
MRD Day 29
<0.01%
0.01 to 0.1%1.9731.1763.3090.01**
0.1 to 1%2.8901.8094.6158.88e-06***
>1%3.9592.2766.8841.10e-06***
BMA Blasts Day 8
≥20% vs. <20%1.7601.2282.5210.002**
BMA Blasts Day 29
≥5% vs. <5%2.0740.8475.0810.111
Down Syndrome
Yes vs. No3.1441.2837.7040.012*
Cell origin
T Cell ALL vs. B Cell ALL0.4830.1971.1840.112
DNA Index
≤0.8 vs. >0.84.6151.69912.5402.71e-03***

The results of univariate Cox analysis.

HR, hazard ratio; CI, confidence interval; WBC, WBC count at diagnosis; CNS Status, CNS status at diagnosis; * P < 0.05; ** P < 0.01; *** P < 0.001.

TABLE 3

Clinical characteristicsHRLower 95% CIUpper 95% CIP-value
Age
≥10 years vs. <10 years1.1640.7841.7260.452
ETV6 RUNX1 Fusion Status
Positive vs. Negative0.2930.1240.6950.005**
TRISOMY 4 10 Status
Positive vs. Negative0.8660.1754.2800.860
WBC (×103/mcL)
<50
50 to 1001.2610.7872.0180.335
>1001.5390.9122.6270.114
CNS Status
CNS1
CNS21.4650.9372.2930.094
CNS30.5280.1641.6940.283
TCF3 PBX1 Status
Positive vs. Negative1.4580.8122.6160.206
BCR ABL1 Status
Positive vs. Negative1.3910.5533.5040.483
MRD Day 29
<0.01%
0.01 to 0.1%1.9421.1333.3300.016*
0.1 to 1%2.1401.2623.6280.005**
>1%2.5081.3454.6750.004**
BMA Blasts Day 8
≥20% vs. <20%1.3540.8902.0580.157
Down Syndrome
Yes vs. No2.3510.9405.8850.068
DNA Index
≤0.8 vs. >0.83.6171.23310.6120.019*
Karyotype
No trisomies in 4, 10, 17, 18
Only trisomies 4 or 10 or 17 or 180.5150.2071.2780.152
Double trisomies (DT)0.3080.0741.2890.107
Triple trisomies (TT)0.2110.0490.907 0.036*

The results of multivariate Cox analysis.

HR, hazard ratio; CI, confidence interval; WBC, WBC count at diagnosis; CNS Status, CNS status at diagnosis; * P < 0.05; ** P < 0.01.

Prognostic Nomogram for Overall Survival

The prognostic nomogram that integrated all significant independent factors for overall survival (OS) in the primary cohort is shown in Figure 2. The nomogram demonstrated that the karyotype had the largest contribution to prognosis, followed by ETV6 RUNX1 fusion status, DNA index, and CNS status. BMA blast day 8 and MRD day 29 showed a moderate effect on survival rate. Each category within these variables was assigned a point on the top scale based on the coefficients from multivariate Cox regression. By summing all of the assigned points for the eight variables and drawing a vertical line between total points and survival probability axis, we were easily able to obtain the estimated probability of 3-, 5-, and 7-year survival. The C-index for OS prediction were 0.754 (95% CI, 0.715–0.793) and 0.731 for the primary cohorts and bootstrapping validation, respectively. The calibration plot for the probability of survival at 3, 5, and 7 years showed an optimal agreement between the predictions by the nomogram and the actual observations (Figure 3).

FIGURE 2

FIGURE 3

ROC and Kaplan–Meier Curve Analyses

We calculated the risk score, which was composed of six clinical features in the nomogram, based on the regression coefficients from the multivariate Cox regression analysis. We divided all patients in the training set into two groups, the high-risk score group (n = 336) and the low-risk score group (n = 337), by taking the corresponding median risk score as the cutoff point. Next, we compared survival predictions with regard to specificity and sensitivity according to the risk scores and clinical characteristics in the nomogram by ROC curve analysis. The results show that the areas under the curve (AUCs) of the risk score for 3-, 5-, and 7-year survival were 0.775, 0.776, and 0.772, respectively, which were higher than those of all of the clinical factors in the nomogram (Figures 4A–C). Patients with a high-risk score had a markedly worse OS than patients with a low-risk score (log-rank test P < 0.0001, Figure 4D).

FIGURE 4

Survival Analysis for Subgroups According to Variable in Nomogram

In general, subtype patients with positive ETV6 RUNX1 fusion gene showed better OS than patients with negative ETV6 RUNX1 fusion gene, and patients in the high-risk group had worse OS than those in the low-risk group in ETV6 RUNX1 fusion-negative segment, while in the ETV6 RUNX1 fusion gene-positive segment, all patients were in the low-risk group (Figure 5A). Besides, the Kaplan–Meier curves of OS suggested that hypodiploid patients had significantly worse OS than patients not. What is more, patients in the high-risk group had worse OS than those in the low-risk group in non-hypodiploid segment, and all patients of the hypodiploid segment were in the high-risk group (Figure 5B). These above results indicated that the risk scores we have obtained are quite accurate and predictive.

FIGURE 5

Validated Nomogram in the Independent and Total Cohorts

To validate the robustness of our nomogram, we reviewed the data from the TARGET database again retrospectively, and the filtering process is presented in Supplementary Figure S1. We ended up with 299 valid cases of data as independent validation queues. Through a similar analysis process, we get supportive results. As Supplementary Figure S2 shows, in the validation cohort, the AUCs of the risk score for 3-, 5-, and 7-year survival were 0.683, 0.723, and 0.737, respectively, which were higher than those of all of the clinical factors in the nomogram (Supplementary Figure S2D). Patients with a high-risk score had a markedly worse OS than patients with a low-risk score (log-rank test P < 0.0001, Supplementary Figure S2A). The C-index of the nomogram for predicting survival was 0.703 (95% CI, 0.640–0.766). Besides, to further confirm the robustness of the model, we also validate the nomogram in the total cohort (N = 299 + 673 = 972). Patients with a high-risk score had a markedly worse OS than patients with a low-risk score (log-rank test P < 0.0001, Supplementary Figure S2C), and the AUCs of the risk score for 3-, 5-, and 7-year survival were 0.753, 0.783, and 0.773 (Supplementary Figure S2B), respectively. The C-index of the nomogram for predicting survival was 0.723 (95% CI, 0.684–0.762).

Discussion

In the present study, we downloaded the data from TARGET database and screened patients with complete information of main observation indicators. Then, we identified independent prognostic factors by univariate and multivariate Cox regression analyses. Significant factors such as ETV6 RUNX1 fusion status, karyotype, MRD day 29, Down syndrome, and DNA index were applied to construct a prognostic nomogram. C-index, Kaplan–Meier analyses, ROC curves and AUC values show that the nomogram objectively and accurately predicted the prognosis of patients with cALL.

A large study showed that cases with trisomy of chromosomes 4, 10, 17, and 18 appear to have the most favorable outcome (10, 3335). Besides, Harris et al. (9) found that, among patients with a DNA index >1.16, patients with trisomies of both chromosomes 4 and 10 had a 4-year EFS of 96.6% (n = 161, SE = 3.8%), whereas patients with neither or only one of these trisomies had a 4-year EFS of 70.4% (n = 73, SE = 11.5%). Convincingly, the Kaplan–Meier curves of OS suggested that patients with TT had significantly better OS than patients not in this study (HR = 0.211, P = 0.036).

The ETV6-RUNX1 fusion gene, which grew out of t(12; 21) (p13; q22) translocation, is the most common chromosome translocation abnormality among cALL. Rubnitz et al. (36, 37) found that the positive frequency of ETV6-RUNX1 in newly diagnosed and recurrent children was 25%, and the 5-year event-free survival (EFS) survival rate of positive children was more than 90%. A study with an average follow-up time of 8 years showed that the 5-year EFS of 244 ETV6-RUNX1-positive ALL children accounted (86 ± 2)%, while that of the ETV6-RUNX1-negative B-ALL children was (72 ± 2)% (38). Obviously, the fusion gene, ETV6-RUNX1, is associated with a favorable prognosis. Similarly, our results showed that the OS rate of ETV6-RUNX1-positive group was significantly higher than that of the group with negative ETV6-RUNX1 (HR = 0.293, P = 0.005).

Hypodiploid acute lymphoblastic leukemia (<44 chromosomes) comprises two subtypes with distinct transcriptional profiles and genetic alterations (33). Numerous studies have shown that hypodiploid (chromosome number ≤44) or DNA index <0.8 is a high-risk type for patients with cALL (13, 39, 40). Low-hypodiploid acute lymphoblastic leukemia has a very poor outcome (39). The frequency increases with age, from being extremely rare in children (<1%), to 5% in adolescents and young adults, and over 10% in adults (41). In our study, hypodiploid accounted for 1.52%, and OS was significantly inferior to non-hypodiploid (HR = 3.617, P = 0.019).

Current research suggests that MRD may be the main cause of relapse (42). MRD refers to the state of trace tumor cells that cannot be detected morphologically in vivo after the induction of remission or bone marrow transplantation in children with leukemia. It is considered to be a more objective and sensitive assessment of the specificity of the clinical treatment response and disease control. In the CCLG-ALL-2008 program of the Chinese Children’s Leukemia Collaborative Group, MRD has been used as an important indicator for risk stratification, and multiple studies have also shown that MRD can be used as an independent prognostic factor (4244). In our study, KM analysis (log-rank test P < 0.0001, Supplementary Figure S4) and univariable COX regression analysis suggest that MRD has a great influence on survival of cALL patients. The presence of day 29 marrow MRD was associated with shorter OS in all risk groups; even patients with 0.01–0.1% day 29 MRD had poor outcome compared with patients negative for MRD patients (80.1 vs. 88.9% 5-year OS). Besides, multivariate COX regression analysis suggests that MRD is an independent risk factor in cALL patients, which was consistent with previous research results (45).

According to the proposed nomogram, we are able to estimate the 3, 5, and 7 years survival rate in patients with cALL, for example, a patient (TARGET-10-PARCVT) with negative TEL-AML1 (corresponds to 79 points), no trisomies in 4, 10, 17, and 18 (corresponds to 100 points), 0.8% of MRD day 29 (corresponds to 49 points), no Down Syndrome (corresponds to 0 points), and hypodiploid (corresponds to 77 points). The calculation according to the proposed nomogram is thus 79 + 100 + 49 + 77 = 305 points, predicting a 5-year survival rate of 17.5% postoperatively. In fact, she died with an OS of 237 days. Another example is that of a patient (TARGET-10-PAMEEK) with negative TEL-AML1 (corresponds to 79 points), no trisomies in 4, 10, 17, and 18 (corresponds to 100 points), 12.6% of MRD day 29 (corresponds to 57 points), no Down syndrome (corresponds to 0 points), and DNA index = 1 (corresponds to 0 points). The calculation according to the proposed nomogram is thus 79 + 100 + 57 = 236 points, predicting a 5-year survival rate of 60.1% postoperatively. In fact, she is still alive with an OS of 60 days, but he could be in danger due to the high value of MRD day 29. If the patient used the scoring system and control his MRD value <0.01%, his total score would be 179, with a 5-year survival rate of about 83.5%. There is a great significance for guiding clinical stratified treatment. Low-risk patients can appropriately reduce the intensity of treatment and do not need to do allogeneic bone marrow transplantation. Besides, high-risk patients need to consider more actively to do transplant and strengthen the consolidation of treatment after induction and remission.

Nomograms are a commonly used tool in oncology that can be used to calculate an individual probability by integrating diverse prognostic and determinant variables according to corresponding clinical characteristics (46, 47). At present, the study on the prognosis of cALL patients focuses on individual factors and lacks a prognostic model covering a comprehensive range of factors. In this study, a prognostic nomogram combining clinical factors was established. The clinical factors in the nomogram are not affected by researchers and can be easily obtained. Moreover, our nomogram had a better predictive accuracy than that of each factor alone. However, the limitation of this study is the lack of external validation. Due to the lack of the number of patients and the corresponding information, it is difficult to obtain relevant resources in public databases or disease centers. Multicenter prospective cohort study may predict patient’s prognosis more accurately.

Conclusion

In conclusion, we comprehensively evaluated the risk of clinical factors associated with prognosis and carried out risk stratification. The nomogram proposed in this study objectively and accurately predicted the prognosis of children with ALL. This nomogram may be a useful tool that can help clinicians develop personalized treatment plans, thereby effectively improving the survival rate of cALL patients.

Statements

Data availability statement

All datasets presented in this study are included in the article/Supplementary Material.

Ethics statement

As the data (TARGET datasets) are publicly available, no ethical approval was required.

Author contributions

RM, YL, and TZ conceived the project and designed the experiments. RM and SH wrote the manuscript. FD, YCZ, and YZ carried out the statistical analysis. YL, RM, and TZ contributed to manuscript revision. All authors provided suggestions during manuscript preparation and read the final version.

Funding

This work was supported by grants from the National Natural Science Foundation of China (81502075) and the Foundation of Science and Technology of Sichuan Province (2019YJ0635). The funders had no role in the study design and implementation.

Acknowledgments

The authors would like to thank the efforts of the National Cancer Institute (NCI)’s Office of Cancer Genomics (OCG).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.01673/full#supplementary-material

Abbreviations

  • ALL

    acute lymphoblastic leukemia

  • AUC

    area under the curve

  • cALL

    childhood acute lymphoblastic leukemia

  • CI

    confidence interval

  • CNS

    central nervous system

  • DCA

    decision curve analysis

  • HR

    hazard ratio

  • MRD

    minimal residual disease

  • NCI

    National Cancer Institute

  • OCG

    Office of Cancer Genomics

  • ROC

    receiver operating characteristic

  • TARGET

    Therapeutically Applicable Research to Generate Effective Treatment

  • WBC

    white blood cell.

References

  • 1.

    TorpyJMLynmCGlassRM.JAMA patient page. Acute lymphoblastic leukemia.JAMA. (2009) 301:452.

  • 2.

    SoulierJCortesJ.Introduction to the review series on acute lymphoblastic leukemia.Blood. (2015) 125:39656. 10.1182/blood-2015-05-635300

  • 3.

    HungerSPMullighanCG.Acute lymphoblastic leukemia in children.N Engl J Med. (2015) 373:154152. 10.1056/NEJMra1400972

  • 4.

    PuiCHEvansWE.Treatment of acute lymphoblastic leukemia.N Engl J Med. (2006) 354:16678. 10.1056/NEJMra052603

  • 5.

    BuitenkampTDIzraeliSZimmermannMForestierEHeeremaNAvan den Heuvel-EibrinkMMet alAcute lymphoblastic leukemia in children with Down syndrome: a retrospective analysis from the Ponte di Legno study group.Blood. (2014) 123:707. 10.1182/blood-2013-06-509463

  • 6.

    Sanjuan-PlaABuenoCPrietoCAchaPStamRWMarschalekRet alRevisiting the biology of infant t(4;11)/MLL-AF4+ B-cell acute lymphoblastic leukemia.Blood. (2015) 126:267685. 10.1182/blood-2015-09-667378

  • 7.

    EldforsSKuusanmäkiHKontroMMajumderMMParsonsAEdgrenHet alIdelalisib sensitivity and mechanisms of disease progression in relapsed TCF3-PBX1 acute lymphoblastic leukemia.Leukemia. (2017) 31:517. 10.1038/leu.2016.202

  • 8.

    SchultzKRPrestidgeTCamittaB.Philadelphia chromosome-positive acute lymphoblastic leukemia in children: new and emerging treatment options.Expert Rev Hematol. (2010) 3:73142. 10.1586/ehm.10.60

  • 9.

    HarrisMBShusterJJCarrollALookATBorowitzMJCristWMet alTrisomy of leukemic cell chromosomes 4 and 10 identifies children with B-progenitor cell acute lymphoblastic leukemia with a very low risk of treatment failure: a Pediatric Oncology Group study.Blood. (1992) 79:331624.

  • 10.

    DastugueNSuciuSPlatGSpelemanFCavéHGirardSet alHyperdiploidy with 58-66 chromosomes in childhood B-acute lymphoblastic leukemia is highly curable: 58951 CLG-EORTC results.Blood. (2013) 121:241523. 10.1182/blood-2012-06-437681

  • 11.

    PetitATrinquandAChevretSBalleriniPCayuelaJMGrardelNet alOncogenetic mutations combined with MRD improve outcome prediction in pediatric T-cell acute lymphoblastic leukemia.Blood. (2018) 131:289300. 10.1182/blood-2017-04-778829

  • 12.

    ParkJHRivièreIGonenMWangXSénéchalBCurranKJet alLong-term follow-up of CD19 CAR therapy in acute lymphoblastic leukemia.N Engl J Med. (2018) 378:44959. 10.1056/NEJMoa1709919

  • 13.

    SafaviSPaulssonK.Near-haploid and low-hypodiploid acute lymphoblastic leukemia: two distinct subtypes with consistently poor prognosis.Blood. (2017) 129:4203. 10.1182/blood-2016-10-743765

  • 14.

    O’ConnorDEnshaeiABartramJHancockJHarrisonCJHoughRet alGenotype-specific minimal residual disease interpretation improves stratification in pediatric acute lymphoblastic leukemia.J Clin Oncol. (2018) 36:3443. 10.1200/JCO.2017.74.0449

  • 15.

    ZhengRXiaoxiaPHongmeiZSiweiZTianhuiCHuanminWet alIncidence, mortality and survival of childhood cancer in China during 2000-2010 period: a population-based study.Cancer Lett. (2015) 363:17680.

  • 16.

    WardEDeSantisCRobbinsAKohlerBJemalA.Childhood and adolescent cancer statistics.CA Cancer J Clin. (2014) 64:83103.

  • 17.

    CheungWYRenfroLAKerrDde GramontASaltzLBGrotheyAet alDeterminants of early mortality among 37,568 patients with colon cancer who participated in 25 clinical trials from the adjuvant colon cancer endpoints database.J Clin Oncol. (2016) 34:11829. 10.1200/JCO.2015.65.1158

  • 18.

    WilliamsSBKamatAMSamanthaCRichardJSLaurenceCet alEORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive stage Ta-T1 urothelial bladder cancer patients treated with 1-3 years of maintenance bacillus calmette-guérin.Eur Urol. (2016) 69:e1212. 10.1016/j.eururo.2016.01.053

  • 19.

    TsoutsouPGJeanneret SozziWMatzingerOOzsahinM.Nomograms predicting locoregional recurrence in the subtype era of breast cancer.J Clin Oncol. (2013) 31:6478. 10.1200/JCO.2012.46.8587

  • 20.

    SternbergCN.Are nomograms better than currently available stage groupings for bladder cancer.J Clin Oncol. (2006) 24:381920. 10.1200/JCO.2006.07.1290

  • 21.

    MarianiLMiceliRKattanMWBrennanMFColecchiaMFioreMet alValidation and adaptation of a nomogram for predicting the survival of patients with extremity soft tissue sarcoma using a three-grade system.Cancer. (2005) 103:4028. 10.1002/cncr.20778

  • 22.

    WangLHricakHKattanMWChenHNKuroiwaKEisenbergHFet alPrediction of seminal vesicle invasion in prostate cancer: incremental value of adding endorectal MR imaging to the Kattan nomogram.Radiology. (2007) 242:1828. 10.1148/radiol.2421051254

  • 23.

    BrownPInabaHAnnesleyCBeckJColaceSDallasMet alPediatric acute lymphoblastic leukemia, version 2.2020, NCCN clinical practice guidelines in oncology.J Natl Compr Canc Netw. (2020) 18:81112. 10.6004/jnccn.2020.0001

  • 24.

    ToftNBirgensHAbrahamssonJGriškevièiusLHallböökHHeymanMet alResults of NOPHO ALL2008 treatment for patients aged 1-45 years with acute lymphoblastic leukemia.Leukemia. (2018) 32:60615. 10.1038/leu.2017.265

  • 25.

    OskarssonTSöderhällSArvidsonJForestierEMontgomerySBottaiMet alRelapsed childhood acute lymphoblastic leukemia in the Nordic countries: prognostic factors, treatment and outcome.Haematologica. (2016) 101:6876. 10.3324/haematol.2015.131680

  • 26.

    EguigurenJMSchellMJCristWMKunkelKRiveraGK.Complications and outcome in childhood acute lymphoblastic leukemia with hyperleukocytosis.Blood. (1992) 79:8715.

  • 27.

    LoweEJPuiCHHancockMLGeigerTLKhanRBSandlundJT.Early complications in children with acute lymphoblastic leukemia presenting with hyperleukocytosis.Pediatr Blood Cancer. (2005) 45:105. 10.1002/pbc.20178

  • 28.

    BorowitzMJDevidasMHungerSPBowmanWPCarrollAJCarrollWLet alClinical significance of minimal residual disease in childhood acute lymphoblastic leukemia and its relationship to other prognostic factors: a Children’s Oncology Group study.Blood. (2008) 111:547785. 10.1182/blood-2008-01-132837

  • 29.

    NyvoldCMadsenHORyderLPSeyfarthJSvejgaardAClausenNet alPrecise quantification of minimal residual disease at day 29 allows identification of children with acute lymphoblastic leukemia and an excellent outcome.Blood. (2002) 99:12538. 10.1182/blood.v99.4.1253

  • 30.

    KramerAAZimmermanJE.Assessing the calibration of mortality benchmarks in critical care: the Hosmer-Lemeshow test revisited.Crit Care Med. (2007) 35:20526. 10.1097/01.CCM.0000275267.64078.B0

  • 31.

    ElhaiMMeuneCBoubayaMAvouacJHachullaEBalbir-GurmanAet alMapping and predicting mortality from systemic sclerosis.Ann Rheum Dis. (2017) 76:1897905. 10.1136/annrheumdis-2017-211448

  • 32.

    RenfroLAGoldbergRMGrotheyASobreroAAdamsRSeymourMTet alClinical calculator for early mortality in metastatic colorectal cancer: an analysis of patients from 28 clinical trials in the aide et recherche en cancérologie digestive database.J Clin Oncol. (2017) 35:192937. 10.1200/JCO.2016.71.5771

  • 33.

    MalardFMohtyM.Acute lymphoblastic leukaemia.Lancet. (2020) 395:114662. 10.1016/S0140-6736(19)33018-1

  • 34.

    MoormanAVRichardsSMMartineauMCheungKLRobinsonHMJalaliGRet alOutcome heterogeneity in childhood high-hyperdiploid acute lymphoblastic leukemia.Blood. (2003) 102:275662. 10.1182/blood-2003-04-1128

  • 35.

    HeeremaNASatherHNSenselMGZhangTHutchinsonRJNachmanJBet alPrognostic impact of trisomies of chromosomes 10, 17, and 5 among children with acute lymphoblastic leukemia and high hyperdiploidy (> 50 chromosomes).J Clin Oncol. (2000) 18:187687. 10.1200/JCO.2000.18.9.1876

  • 36.

    RubnitzJEDowningJRPuiCHShurtleffSARaimondiSCEvansWEet alTEL gene rearrangement in acute lymphoblastic leukemia: a new genetic marker with prognostic significance.J Clin Oncol. (1997) 15:11507. 10.1200/JCO.1997.15.3.1150

  • 37.

    RubnitzJEBehmFGWichlanDRyanCSandlundJTRibeiroRCet alLow frequency of TEL-AML1 in relapsed acute lymphoblastic leukemia supports a favorable prognosis for this genetic subgroup.Leukemia. (1999) 13:1921. 10.1038/sj.leu.2401257

  • 38.

    RubnitzJEWichlanDDevidasMShusterJLindaSBKurtzbergJet alProspective analysis of TEL gene rearrangements in childhood acute lymphoblastic leukemia: a Children’s Oncology Group study.J Clin Oncol. (2008) 26:218691. 10.1200/JCO.2007.14.3552

  • 39.

    PuiCHReboraPSchrappeMAttarbaschiABaruchelABassoGet alOutcome of children with hypodiploid acute lymphoblastic leukemia: a retrospective multinational study.J Clin Oncol. (2019) 37:7709. 10.1200/JCO.18.00822

  • 40.

    PuiCHWilliamsDLRaimondiSCRiveraGKLookATDodgeRKet alHypodiploidy is associated with a poor prognosis in childhood acute lymphoblastic leukemia.Blood. (1987) 70:24753.

  • 41.

    HolmfeldtLWeiLDiaz-FloresEWalshMZhangJDingLet alThe genomic landscape of hypodiploid acute lymphoblastic leukemia.Nat Genet. (2013) 45:24252. 10.1038/ng.2532

  • 42.

    BeldjordKChevretSAsnafiVHuguetFBoullandMLLeguayTet alOncogenetics and minimal residual disease are independent outcome predictors in adult patients with acute lymphoblastic leukemia.Blood. (2014) 123:373949. 10.1182/blood-2014-01-547695

  • 43.

    O’ConnorDMoormanAVWadeRHancockJTanRMBartramJet alUse of minimal residual disease assessment to redefine induction failure in pediatric acute lymphoblastic leukemia.J Clin Oncol. (2017) 35:6607. 10.1200/JCO.2016.69.6278

  • 44.

    PietersRde Groot-KrusemanHVan der VeldenVFioccoMvan den BergHde BontEet alSuccessful therapy reduction and intensification for childhood acute lymphoblastic leukemia based on minimal residual disease monitoring: study all10 from the dutch childhood oncology group.J Clin Oncol. (2016) 34:2591601. 10.1200/JCO.2015.64.6364

  • 45.

    VoraAGouldenNWadeRMitchellCHancockJHoughRet alTreatment reduction for children and young adults with low-risk acute lymphoblastic leukaemia defined by minimal residual disease (UKALL 2003): a randomised controlled trial.Lancet Oncol. (2013) 14:199209. 10.1016/S1470-2045(12)70600-9

  • 46.

    HanDSSuhYSKongSHLeeHJChoiYAikouSet alNomogram predicting long-term survival after d2 gastrectomy for gastric cancer.J Clin Oncol. (2012) 30:383440. 10.1200/JCO.2012.41.8343

  • 47.

    WangTSSosaJA.Thyroid gland: can a nomogram predict death in patients with thyroid cancer.Nat Rev Endocrinol. (2013) 9:1923. 10.1038/nrendo.2013.45

Summary

Keywords

childhood acute lymphoblastic leukemia, prognosis, nomogram, therapeutically applicable research to generate effective treatment database, genetic predisposition

Citation

Mao R, Hu S, Zhang Y, Du F, Zhang Y, Liu Y and Zhang T (2020) Prognostic Nomogram for Childhood Acute Lymphoblastic Leukemia: A Comprehensive Analysis of 673 Patients. Front. Oncol. 10:1673. doi: 10.3389/fonc.2020.01673

Received

04 May 2020

Accepted

29 July 2020

Published

10 September 2020

Volume

10 - 2020

Edited by

Annalisa Lonetti, University of Bologna, Italy

Reviewed by

Denis Martin Schewe, University of Kiel, Germany; Valentina Indio, Università degli Studi di Bologna, Italy

Updates

Copyright

*Correspondence: Yanjun Liu, Tongtong Zhang,

These authors have contributed equally to this work

This article was submitted to Hematologic Malignancies, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics