- Open Access
Chemotherapy-related hyperbilirubinemia in pediatric acute lymphoblastic leukemia: a genome-wide association study from the AIEOP-BFM ALL study group
Journal of Experimental & Clinical Cancer Research volume 42, Article number: 21 (2023)
Characterization of clinical phenotypes in context with tumor and host genomic information can aid in the development of more effective and less toxic risk-adapted and targeted treatment strategies. To analyze the impact of therapy-related hyperbilirubinemia on treatment outcome and to identify contributing genetic risk factors of this well-recognized adverse effect we evaluated serum bilirubin levels in 1547 pediatric patients with acute lymphoblastic leukemia (ALL) and conducted a genome-wide association study (GWAS).
Patients and methods
Patients were treated in multicenter trial AIEOP-BFM ALL 2000 for pediatric ALL. Bilirubin toxicity was graded 0 to 4 according to the Common Toxicity Criteria (CTC) of the National Cancer Institute. In the GWAS discovery cohort, including 650 of the 1547 individuals, genotype frequencies of 745,895 single nucleotide variants were compared between 435 patients with hyperbilirubinemia (CTC grades 1-4) during induction/consolidation treatment and 215 patients without it (grade 0). Replication analyses included 224 patients from the same trial.
Compared to patients with no (grade 0) or moderate hyperbilirubinemia (grades 1-2) during induction/consolidation, patients with grades 3-4 had a poorer 5-year event free survival (76.6 ± 3% versus 87.7 ± 1% for grades 1-2, P = 0.003; 85.2 ± 2% for grade 0, P < 0.001) and a higher cumulative incidence of relapse (15.6 ± 3% versus 9.0 ± 1% for grades 1-2, P = 0.08; 11.1 ± 1% for grade 0, P = 0.007). GWAS identified a strong association of the rs6744284 variant T allele in the UGT1A gene cluster with risk of hyperbilirubinemia (allelic odds ratio (OR) = 2.1, P = 7 × 10− 8). TT-homozygotes had a 6.5-fold increased risk of hyperbilirubinemia (grades 1-4; 95% confidence interval (CI) = 2.9-14.6, P = 7 × 10− 6) and a 16.4-fold higher risk of grade 3-4 hyperbilirubinemia (95% CI 6.1-43.8, P = 2 × 10− 8). Replication analyses confirmed these associations with joint analysis yielding genome-wide significance (allelic OR = 2.1, P = 6 × 10− 11; 95% CI 1.7-2.7). Moreover, rs6744284 genotypes were strongly linked to the Gilbert’s syndrome-associated UGT1A1*28/*37 allele (r2 = 0.70), providing functional support for study findings. Of clinical importance, the rs6744284 TT genotype counterbalanced the adverse prognostic impact of high hyperbilirubinemia on therapy outcome.
Chemotherapy-related hyperbilirubinemia is a prognostic factor for treatment outcome in pediatric ALL and genetic variation in UGT1A aids in predicting the clinical impact of hyperbilirubinemia.
The systematic evaluation of tumor and host genomic information can help to identify new predictive markers and cancer vulnerabilities and lead to novel effective and less toxic risk-adapted and targeted treatment strategies . During the last five decades, treatment of pediatric patients with acute lymphoblastic leukemia (ALL) has significantly improved and is one of the success stories in clinical oncology . However, although most children diagnosed with ALL can be cured by contemporary clinical protocols [2, 3], benefits in survival come with acute and long-term adverse effects. This may complicate administration of therapy or impact on health and quality of life during and after completion of treatment [2, 3]. Thus, there is continuing medical need to improve risk assessment and tailored treatment for children and adolescents with ALL.
Hepatotoxicity, mainly captured in clinical trials by assessment of hyperbilirubinemia and transaminasemia, is a well-known complication frequently occurring in the early phases of ALL treatment [4, 5]. Surprisingly, detailed reports describing incidence and impact of hepatotoxicity on the overall outcome of pediatric ALL are lacking. A recent study including 1872 pediatric ALL patients treated according to the ALL IC-BFM 2002 protocol reported 934 grade 3 or above events of hepatotoxicity according to the National Cancer Institute (NCI) Common Toxicity Criteria (CTC; for details see Suppl. Methods section)  in 527 individuals (28%) during the entire therapy . A previous single institution study reported that event-free survival (EFS) was rarely influenced by hyperbilirubinemia, whereas treatment modifications including delays were common . Overall, 17% of patients in the latter study had bilirubin levels equal or above CTC grade 3 at least once in the course of therapy . Another study determined such high levels in 10, 2 and 15% of patients during induction, consolidation and maintenance phases, respectively .
High total serum bilirubin levels – with or without transaminasemia – may indicate liver dysfunction upon exposure to various antileukemic drugs including L-asparaginase or antimetabolites [8,9,10]. However, in this context the biological basis of hyperbilirubinemia is poorly understood. In the last decade, genome-wide association studies (GWAS) enabled a refined risk assessment of many diseases – including pediatric ALL – by identifying risk-associated genetic variants [11,12,13].
In the present study, we evaluated the clinical impact of hyperbilirubinemia in a large cohort of pediatric ALL patients and applied a GWAS approach to identify genetic variants influencing chemotherapy-related hyperbilirubinemia.
Patients included were 1 to 18 years of age at diagnosis of ALL and enrolled in the German part of the European AIEOP-BFM ALL 2000 multicenter clinical trial for frontline treatment of pediatric ALL from August 1999 to November 2005 [14,15,16] (for treatment details see Suppl. Table 1). Primary patient selection criteria were the availability of bilirubin toxicity gradings for induction or consolidation – protocols IA/IB – and, for the assembly of our GWAS discovery cohort, the availability of genome-wide germline genotyping information (Suppl. Fig. 1) obtained during a previous project . When available, we also assessed bilirubin gradings related to subsequent treatment phases (extracompartment therapy, high-risk block treatment, and re-intensification protocols II and III). No information was available for maintenance treatment. For further details including statistical analyses, see Suppl. Methods section.
As part of the routine safety management, toxicity was assessed for all treatment elements except for interim maintenance and maintenance phases. Considering 17.1 μmol/L as the upper normal limit (UNL), total bilirubin serum levels were graded according to the CTC of the NCI, version 2  (for details see Suppl. Methods section).
In the discovery cohort, DNA obtained from bone marrow samples in morphological remission was genotyped on Human Omni1-Quad v1 arrays (Illumina, San Diego, CA, USA) as previously described . Genotype information on rs6744284 for replication purposes was derived from a preceding analysis using Affymetrix Genome-wide Human SNP arrays 5.0 (Affymetrix, South San Francisco, CA, USA) .
Out of the discovery cohort, we genotyped 544 patients with additional DNA available for UGT1A1*28/*37 variations using a pre-developed KASP assay (LGC Biosearch Technologies, Hoddesdon, UK).
Genome-wide association study
GWAS conduction, including data pruning, association testing with plink 1.9 (www.cog-genomics.org/plink/1.9/), evaluation and plotting, was realized with an in-house developed R (3.6.0) script using RStudio (1.0.143). A strict quality control was performed prior to association testing. For details on data pruning, association testing, genotype imputation and further statistical analyses see Suppl. Methods section.
Clinical characterization of hyperbilirubinemia
Within the 1547 patients of the AIEOP ALL-BFM 2000 study population with available toxicity information, 540 (34.9%) had normal and 1007 (65.1%) had increased bilirubin levels during induction/consolidation (protocols IA/IB) (Table 1). This included 825 (53.3%) patients with moderate hyperbilirubinemia (grades 1-2; 707 during induction, 575 during consolidation, 355 in both phases) and 182 (11.8%) with high hyperbilirubinemia (grades 3-4; 158 reports during induction, 50 during consolidation, 26 in both phases). Comparing patients without to patients with hyperbilirubinemia (grade 0 vs. grades 1-4), we noticed a larger proportion of older patients (P < 0.001) and more T cell ALL patients (P = 0.022) among those affected. The group of B cell ALL patients exhibiting hyperbilirubinemia contained fewer hyperdiploid patients (P = 0.002). However, no differences with regard to other genetic subgroups were observed (ETV6-RUNX1, BCR-ABL1 and KMT2A-AFF1) (Table 1). Considering the entire course of therapy, 245 patients (16.0%) had high hyperbilirubinemia in at least one treatment element. We determined a median time to protocol day 78 (after completion of induction and consolidation) of 89 ± 11 days (range 64-185 days), analyzing 1453 of 1547 individuals of our study population with available information. Compared to patients with moderate or no hyperbilirubinemia, patients with hyperbilirubinemia grades 3-4 experienced more therapy delays, requiring 91 days (range 64-174 days) to complete induction/consolidation vs. 88 days (range 70-154 days) for grades 1-2 and 89 days (range 70-185 days) for grade 0 (P = 0.002). No alterations of therapy in response to hyperbilirubinemia were noted.
In outcome analyses, patients with high hyperbilirubinemia (grades 3-4) during induction/consolidation fared significantly worse compared to patients with moderate or no hyperbilirubinemia: 5-year EFS 76.7 ± 3% vs. 87.7 ± 1% (P < 0.0001), and vs. 85.2 ± 2% (P = 0.0031), respectively (Fig. 1A). The corresponding 5-year cumulative incidences of relapse (CIR) were 15.6 ± 3% for high, 9.0 ± 1% for moderate hyperbilirubinemia, and 11.1 ± 1% for patients without hyperbilirubinemia (Fig. 1B).
In our study cohort, 1443 (93%) patients were observed with elevated hepatic transaminase activity levels (grades 1-4) during induction/consolidation (Table 1), 765 of which demonstrated grades 3-4. Transaminase levels were positively correlated with hyperbilirubinemia and, in particular, patients with high hyperbilirubinemia were at risk for concurrent high transaminase levels (grades 3-4) compared to the remaining patients (70% vs. 40%, odds ratio (OR) = 21.3, 95% confidence interval (CI) = 5.1-88.2, P = 2.53 × 10− 5) (Suppl. Table 2).
Interestingly, 5-year EFS and CIR did not differ between patients with moderate (grades 1-2), high (grades 3-4) or absent transaminasemia during induction/consolidation (Suppl. Fig. 2A and B). The 5-year EFS of patients with concurrent high hyperbilirubinemia and high transaminasemia was 76.2 ± 4% and 78.1 ± 6% in patients with high hyperbilirubinemia accompanied by moderate or no transaminasemia (grades 0-2, P = 0.63; Suppl. Fig. 2C). The corresponding CIR were 17.5 ± 3% and 11.0 ± 4% (P = 0.13; Suppl. Fig. 2D).
Multivariate analyses including established prognostic factors in AIEOP-BFM trials identified high hyperbilirubinemia as an independent predictor of outcome, while severe transaminasemia (CTC grade 4) did not demonstrate an impact here (Table 2).
Characteristics of the GWAS discovery cohort
The finally pruned GWAS discovery cohort included 650 of 1547 patients with available genome-wide genotyping information (Suppl. Table 3). Both the distribution of hyperbilirubinemia and the clinical characteristics were comparable to those of the entire study population (Table 1 and Suppl. Tables 3 and 4).
Of the 650 patients in this discovery cohort, 215 (33%) patients had normal and 435 (67%) had increased bilirubin levels during induction/consolidation: 367 (56.5%) patients demonstrated moderate hyperbilirubinemia (grades 1-2; 313 during induction, 248 during consolidation, 159 in both phases) and 68 (10.5%) had high hyperbilirubinemia (grades 3-4; 59 reports during induction, 21 during consolidation, 12 in both phases). Pre-treatment hyperbilirubinemia at diagnosis was rare, but more frequent among patients developing chemotherapy-related hyperbilirubinemia during induction/consolidation compared to those not (7.5% (19/252) vs. 1.7% (2/116), P = 0.025). Considering the entire course of therapy, 91 patients (14.0%) had high hyperbilirubinemia in at least one treatment element.
Genome-wide association study
When comparing the 435 patients with hyperbilirubinemia (grades 1-4) during induction/consolidation to the 215 patients with normal bilirubin levels, the five loci most associated with this phenotype were the UGT1A gene cluster, MARK2P5, SULF2, MIR924HG and USH2A (Suppl. Table 5). The strongest associations were observed for variants residing in the UGT1A (UDP glucuronosyltransferase family 1 member A) locus at 2q37. The only variant reaching near genome-wide significance (OR = 2.1, 95% CI = 1.6-2.7, P = 7.3 × 10− 8), rs6744284 was also the index SNV of a 189 kb region of high linkage disequilibrium (LD). Information on the complex UGT1A cluster with its overlapping genes and results of imputation are presented in supplementary material (Suppl. Fig. 3, Suppl. Tables 6 and 7).
To examine whether inclusion of age and immunophenotype would influence allelic association, we compared results from crude and adjusted logistic regression analyses. We did not detect any differences in conferred risk for the variant rs6744284 T allele with reference to the wild-type C allele (unadjusted allelic OR = 2.1, 95% CI = 1.6-2.7, P = 1.8 × 10− 7; adjusted OR = 2.1, 95% CI = 1.6-2.8, P = 1.2 × 10− 7) (Suppl. Table 5).
The genotypic association of rs6744284 with frequency and risk of hyperbilirubinemia during induction/consolidation increased stepwise (Table 3). Compared to wild-type patients (CC; 58% with hyperbilirubinemia grades 1-4), heterozygotes (TC; 71% with hyperbilirubinemia grades 1-4) demonstrated a 1.7-fold higher risk of hyperbilirubinemia, while homozygosity for the T allele (TT; 90% with hyperbilirubinemia grades 1-4) conferred an OR of 6.5 (95% CI = 2.9-14.6, P = 7.0 × 10− 6) (Fig. 2A and Table 3). Inclusion of age and immunophenotype as covariates or as stratifying variables did not significantly alter these results (Suppl. Tables 7, 8 and 9). Notably, TT-homozygotes were at particular risk of developing high hyperbilirubinemia (grades 3-4, OR with reference to CC genotype 16.4, 95% CI = 6.1-43.8, P = 2 × 10− 8).
Independent replication analysis
We performed independent replication analyses in a cohort of 224 ETV6-RUNX1-rearranged pediatric ALL patients (Suppl. Table 10), results of which supported our initial GWAS findings. The allelic OR for hyperbilirubinemia (grades 1-4) during induction/consolidation conferred by the variant rs6744284 T allele versus the wild-type C allele was 2.3 (95% CI = 1.5-3.7, P = 2.4 × 10− 4). Genotypic OR in comparison to wild-type patients (CC) were 2.4 (95% CI = 1.3-4.3, P = 3.8 × 10− 3) for heterozygotes (TC) and 6.1 (95% CI = 1.7-21.6, P = 5.6 × 10− 3) for homozygous variant patients (TT) (Table 3). Similar to our findings in the GWAS discovery cohort, patients possessing the rs6744284 TT genotype were at particular risk of high hyperbilirubinemia (OR with reference to CC genotype 13.6, 95% CI = 2.6-71.8, P = 0.002).
Association testing in the combined discovery and replication cohorts resulted in genome-wide significance. Compared to the rs6744284 wild-type allele, presence of the T allele was associated with an OR of 2.1 (CI = 1.7-2.7) for hyperbilirubinemia (grades 1-4) during induction/consolidation at a significance level of P = 5.7 × 10− 11.
UGT1A rs6744284 genotype in subsequent treatment elements
Similar to initial findings, we observed that the rs6744284 TT genotype was also strongly associated with hyperbilirubinemia during extracompartment therapy (Protocol M, OR = 4.1, 95% CI = 2.2-7.9; P < 0.001), re-intensification (Protocols II and III, OR = 9.1, 95% CI = 4.5-18.6, P < 0.001) and high-risk (HR) block elements (OR = 15.3, 95% CI = 1.8-126.6; P = 0.012) (Table 3, Suppl. Fig. 4). Thus, the effect of rs6744284 on risk of hyperbilirubinemia was not limited to early chemotherapy, but was generalizable to all intensive treatment phases for pediatric ALL.
UGT1A rs6744284 and Gilbert’s syndrome-associated variants
The UGT1A enzyme family is crucial for bilirubin glucuronidation and related impairing genetic alterations form the mechanistic basis of the Gilbert’s syndrome (GS) [18,19,20]. Therefore, we genotyped the GS-related functional genetic variations UGT1A1*28 and *37 [21, 22] in 544 (84%) patients of our discovery cohort with available remission DNA (Suppl. Table 11). Comparable to rs6744284, we observed a strong association with hyperbilirubinemia: the allelic OR for *28/*37 vs. wild-type (*1) was 1.9 (95% CI = 1.4-2.5, P = 5.0 × 10− 6). Genotype-based analyses demonstrated a stepwise increase of frequency and risk of hyperbilirubinemia for the variant alleles. Out of 544 patients 62 (11%) were homozygous for either UGT1A1*28/*28 or *37/*37 – this cannot be differentiated by our assay – and had the highest rate and risk of hyperbilirubinemia (89% compared to 58% for *1/*1; OR in comparison to *1/*1 5.8; 95% CI = 2.5-13.3; P = 3.3 × 10− 5) (Suppl. Table 12). Homozygous variant patients were at particular risk of developing high hyperbilirubinemia (grades 3-4, OR = 12.4, 95% CI = 4.4-34.8, P = 1.9 × 10− 6). The strong interrelationship of rs6744284 with UGT1A1*28/*37 is depicted in Fig. 2B. Extended haplotype analyses including eight additional GS-related variants further documented a strong association with rs6744284 (see Suppl. Table 13 and related additional information). Of note, none of the GS-related variants showed a stronger association with hyperbilirubinemia than rs6744284.
Hyperbilirubinemia, transaminase levels, rs6744284 genotype, treatment delay and outcome in the GWAS discovery cohort (n = 650)
Similar to the patients of the entire study population (n = 1547), patients in our discovery cohort with high hyperbilirubinemia during induction/consolidation tended to take 2 days longer to complete consolidation (P = 0.072) (for details see Suppl. Information, page 25). Of interest, we did not observe significant differences between rs6744284 genotypes: 88 days for TT vs. 89 days for CT and 90 days for CC (P = 0.122).
Consistent with the results obtained for the entire study cohort, outcome analyses of the discovery cohort showed that high hyperbilirubinemia during induction/consolidation was associated with a poor 5-year EFS of 71.8 ± 5%, compared to 87.4 ± 2% and 81.7 ± 3% in patients with moderate and without hyperbilirubinemia, respectively (Fig. 3A). The corresponding 5-year CIR were 19.3 ± 5% for high hyperbilirubinemia, 10.7 ± 2% for moderate, and 13.2 ± 2% for no hyperbilirubinemia (Fig. 4A). Although rs6744284 was strongly associated with high hyperbilirubinemia and the proportion of patients with TT genotype among those with high hyperbilirubinemia was 28% (19/68), there were no differences between rs6744284 genotypes related to EFS or CIR (Figs. 3B and 4B) in the discovery cohort. However, within high hyperbilirubinemic patients those carrying the TT genotype had a better EFS (84.2 ± 8% vs 66.9 ± 7%, P = 0.110) and a lower CIR (5.3 ± 5% vs 24.9 ± 6%, P = 0.039) at 5 years compared to the remaining genotypes (TC, CC) (Figs. 3D and 4D).
In the GWAS discovery cohort 604 (93%) patients were observed with elevated hepatic transaminase levels (grades 1-4) during induction/consolidation (Suppl. Table 3), 315 of which demonstrated high grades 3-4. Similar to the results obtained for the 1443 patients in the entire study cohort with available information, transaminase levels were positively associated with hyperbilirubinemia. Especially patients with high hyperbilirubinemia were at increased risk for concurrent high transaminase levels compared to the remaining patients (75% vs. 41%, OR = 16.8, 95% CI = 2.2-127.1, P = 0.006) (Suppl. Table 14). Of importance, transaminasemia was not associated with rs6744284 genotype (P = 0.74). Five-year EFS and CIR did not differ between patients with no (grade 0), moderate (grades 1-2) or high (grades 3-4) transaminasemia during induction/consolidation (Suppl. Fig. 5A and B).
Patients with high hyperbilirubinemia and concurrent high transaminasemia tended to have a higher 5-year EFS of 74.3 ± 6% compared to 64.7 ± 12% in patients with high hyperbilirubinemia accompanied by moderate or no transaminasemia (grades 0-2, P = 0.480; Suppl. Fig. 5C). Corresponding CIR were 21.8 ± 6% and 11.8 ± 8% (P = 0.260, Suppl. Fig. 5D).
Multivariate analyses including established prognostic factors in AIEOP-BFM trials identified high hyperbilirubinemia as an independent predictor of outcome, while rs6744284 TT genotype demonstrated only a tentative protective effect in these analyses (Table 4). However, in multivariate analysis restricted to patients with prognostically unfavorable high hyperbilirubinemia in induction/consolidation, the rs6744284 TT genotype was associated with a statistically significant 14-fold lower relapse-risk compared to rs6744284 wild-type or heterozygous variant patients (CC or TC) (Suppl. Table 15).
Routine hepatotoxicity monitoring in clinical trials for pediatric ALL is typically performed by grading of elevated bilirubin and liver transaminase levels as absent, mild, moderate, severe or life-threatening/fatal according to the NCI CTC criteria. Although the evaluation of laboratory values for hepatotoxicity is common practice, it can be debated how precise such measurements reflect liver dysfunction. To enhance the phenotypic characterization in context with abnormal laboratory values, alternative classifications integrate additional clinical information (e.g., coagulopathy, impairment of liver function-dependent organs) . Similarly, genetic biomarkers hold the potential to aid strategies directed at improved evaluation of hepatotoxicity. Nonetheless, only a few are currently used in clinical routine to guide a genotype-adapted dosing of specific chemotherapeutic agents and thereby reduce adverse reactions (e.g., TPMT with thiopurines , UGT1A1 with irinotecan [25, 26]). In the present study, we applied an unbiased genome-wide approach and identified genetic variation in the UGT1A gene cluster as a major contributor to hyperbilirubinemia associated with chemotherapy for pediatric ALL.
Genetic variation in UGT1A is well-established to affect enzymatic glucuronidation activity and to modulate the metabolism of endogenous metabolites as well as multiple xenobiotics [27,28,29,30,31,32,33,34,35]. Out of nine functional isoforms, only UGT1A1 is relevant for bilirubin glucuronidation [31, 36]. Various functional UGT1A1 variants result in partial or complete reduction of enzymatic activity and determine the phenotype of heritable diseases of bilirubin metabolism [31, 36, 37]. The most common genetic cause for reduced bilirubin conjugation is an insertion polymorphism in the TATA box of UGT1A1 – the (TA)7 variant allele UGT1A1*28; homozygosity for this allele confers a reduced transcriptional activity of 18 to 33% [22, 38], corresponding to the residual glucuronidation activity of ~ 30% determined in patients with Gilbert’s syndrome (GS) [39, 40]. While other associated risk alleles have been described, e.g. UGT1A1*6 (rs4148323) in Asian populations , UGT1A1*28 is by far the most common cause for the Gilbert’s syndrome in Caucasians and African Americans [22, 37]. Our lead SNV, rs6744284, was closely correlated to all GS-related variants assessable in our investigational setting – including UGT1A1*28 – and was the best predictor of hyperbilirubinemia in our patients. These findings are concordant with reports on multi-SNV haplotypes of UGT1A involved in impaired glucuronidation associated with GS [20, 42] and support a potential diagnostic role for the simple assessment of rs6744284.
With relevance to cancer treatment, UGT1A1*28 or UGT1A1*6 homozygotes are low metabolizers of irinotecan and have an increased risk of severe neutropenia, requiring preventive dose adjustments [25, 43]. Likewise, UGT1A*28 homozygotes are recommended to receive reduced doses of belinostat [44, 45]. The UGT1A1 genotype also influences the pharmacokinetics of other glucuronidation-dependent drugs – including several ones used in ALL treatment (e.g., methotrexate, etoposide, and cyclophosphamide) [28, 32, 33, 35, 46, 47]. For example, low UGT1A1 activity was associated with higher plasma methotrexate and bilirubin levels, suggestive of competitive interactions between the three [7, 48]. Similar to our study, these former investigations in the field of pediatric ALL showed that patients with GS were prone to hyperbilirubinemia throughout all treatment phases [7, 48].
One of the previous studies demonstrated that despite higher bilirubin levels, ALL patients with GS did not experience significant treatment modifications, including delays, or worse therapy outcomes , which is in line with our findings. These observations may explain why universal screening for GS in all patients diagnosed with ALL is not routinely performed so far. Although it was recommended previously to screen at least ALL patients with hyperbilirubinemia for GS , common standard recommendations agreed on between international ALL trial consortia do not exist, and may be promoted through the Ponte di Legno initiative .
In our study, high hyperbilirubinemia was an independent prognostic factor negatively affecting EFS and CIR of patients treated on a modern risk-adapted BFM protocol. To our knowledge, this is the first report from a large pediatric ALL cohort receiving relatively homogenous therapy demonstrating an effect of hyperbilirubinemia on long-term treatment outcome. Importantly, the negative prognostic impact was not determined in patients who demonstrated high hyperbilirubinemia and were homozygous for the variant T allele of our lead SNV rs6744284. This observation could have direct clinical implications by helping to differentiate hyperbilirubinemia conferring a negative prognostic impact from hyperbilirubinemia of a less severe clinical phenotype. Moreover, our findings imply that patients with phenotypically relevant genetic variation in UGT1A/GS can be spared from experiencing treatment modifications, as it was suggested previously [48, 49]. Besides replication in other independent clinically and genetically comprehensively characterized cohorts, functional studies are particularly required to augment our knowledge of hyperbilirubinemia as a treatment-related toxicity in pediatric ALL and to elucidate involved pathomechanisms.
Despite interesting perspectives, there are several limitations associated with our study: 1) We were only able to study total serum bilirubin levels. However, a separate analysis of unconjugated and conjugated bilirubin will be important for an improved understanding of therapy-related hyperbilirubinemia and its prognostic value. 2) No data on additional intake (e.g., vitamin B complex, ursodiol) and/or pharmacokinetics of glucuronidation dependent drugs (e.g., methotrexate), potentially confounding our observations, were available to us. 3) Selection bias for inclusion in our GWAS and/or replication cohorts is immanent to a clinical investigation depending on availability of reported hepatotoxicity and biological material. 4) The clinical importance of our findings could be enhanced by systematic collection of information on the potential long-term burden of hepatotoxicity. 5) Finally, the generalizability of our findings to other therapy protocols for pediatric ALL is limited due to differences in medication and timing between them. All of these issues need to be addressed in future investigations and will help to resolve the limitations of our current observations.
High hyperbilirubinemia acted as an independent prognostic factor of therapy outcome in pediatric ALL patients treated on the AIEOP-BFM ALL 2000 protocol. Further, the rs6744284 genotype reliably predicted hyperbilirubinemia throughout all intensive treatment phases. Thus, both assessment of hyperbilirubinemia and UGT1A genotyping will be useful for complementing toxicity risk profiling and optimizing risk-adapted therapeutic strategies for pediatric ALL.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
- AIEOP-BFM ALL study group:
Associazione Italiana Ematologia Oncologia Pediatrica - Berlin Frankfurt-Muenster Acute Lymphoblastic Leukemia study group
Acute lymphoblasic leukemia
Cumulative incidence of relapse
Central nervous system
Common toxicity criteria of the NCI
Event free survival
- ERG :
ETS Transcription Factor ERG
- ETV6 :
ETS Variant Transcription Factor 6
Genomwide Association Study
- IU/m2 :
International Unit per Square Meter
- MARK2P5 :
Microtubule Affinity Regulating Kinase 2 Pseudogene 5
Microliter (1 μL = 1 × 10− 6 l)
Miligram (1 mg =1 × 10− 3 g)
- MIR924HG :
MIR924 Host Gene
Minimal residual disease
National Cancer Institute
- PO :
- RUNX1 :
RUNX family transcription factor 1
Single nucleotide variant
- SULF2 :
- Suppl. :
- TPMT :
- UGT1A :
UDP Glucuronosyltransferase Family 1 Member A
Upper normal limit /upper limit of normal
- USH2A :
white blood cell
Di Martino MT, Meschini S, Scotlandi K, Riganti C, De Smaele E, Zazzeroni F, et al. From single gene analysis to single cell profiling: a new era for precision medicine. J Exp Clin Cancer Res. 2020;39(1):48.
Pui CH, Yang JJ, Hunger SP, Pieters R, Schrappe M, Biondi A, et al. Childhood acute lymphoblastic leukemia: Progress through collaboration. J Clin Oncol. 2015;33(27):2938–48.
Möricke A, Zimmermann M, Reiter A, Henze G, Schrauder A, Gadner H, et al. Long-term results of five consecutive trials in childhood acute lymphoblastic leukemia performed by the ALL-BFM study group from 1981 to 2000. Leukemia. 2010;24(2):265–84.
Denton CC, Rawlins YA, Oberley MJ, Bhojwani D, Orgel E. Predictors of hepatotoxicity and pancreatitis in children and adolescents with acute lymphoblastic leukemia treated according to contemporary regimens. Pediatr Blood Cancer. 2018;65(3):e26891.
Zawitkowska J, Lejman M, Zaucha-Prazmo A, Drabko K, Plonowski M, Bulsa J, et al. Grade 3 and 4 toxicity profiles during therapy of childhood acute lymphoblastic leukemia. In Vivo. 2019;33(4):1333–9.
National Cancer Institute Common Toxicity Criteria (NCI-CTC), version 2.0, DCTD, NCI, NIH, DHHS; 1998. Available from: https://ctep.cancer.gov/protocolDevelopment/electronic_applications/ctc.htm.
Kishi S, Cheng C, French D, Pei D, Das S, Cook EH, et al. Ancestry and pharmacogenetics of antileukemic drug toxicity. Blood. 2007;109(10):4151–7.
Raetz EA, Salzer WL. Tolerability and efficacy of L-asparaginase therapy in pediatric patients with acute lymphoblastic leukemia. J Pediatr Hematol Oncol. 2010;32(7):554–63.
Schmiegelow K, Nielsen SN, Frandsen TL, Nersting J. Mercaptopurine/methotrexate maintenance therapy of childhood acute lymphoblastic leukemia: clinical facts and fiction. J Pediatr Hematol Oncol. 2014;36(7):503–17.
Schmiegelow K, Muller K, Mogensen SS, Mogensen PR, Wolthers BO, Stoltze UK, et al. Non-infectious chemotherapy-associated acute toxicities during childhood acute lymphoblastic leukemia therapy. F1000Res. 2017;6:444.
Bartram T, Burkhardt B, Wossmann W, Seidemann K, Zimmermann M, Cario G, et al. Childhood acute lymphoblastic leukemia-associated risk-loci IKZF1, ARID5B and CEBPE and risk of pediatric non-Hodgkin lymphoma: a report from the Berlin-Frankfurt-Munster study group. Leuk Lymphoma. 2015;56(3):814–6.
Vijayakrishnan J, Qian M, Studd JB, Yang W, Kinnersley B, Law PJ, et al. Identification of four novel associations for B-cell acute lymphoblastic leukaemia risk. Nat Commun. 2019;10(1):5348.
Migliorini G, Fiege B, Hosking FJ, Ma Y, Kumar R, Sherborne AL, et al. Variation at 10p12.2 and 10p14 influences risk of childhood B-cell acute lymphoblastic leukemia and phenotype. Blood. 2013;122(19):3298–307.
Conter V, Bartram CR, Valsecchi MG, Schrauder A, Panzer-Grumayer R, Möricke A, et al. Molecular response to treatment redefines all prognostic factors in children and adolescents with B-cell precursor acute lymphoblastic leukemia: results in 3184 patients of the AIEOP-BFM ALL 2000 study. Blood. 2010;115(16):3206–14.
Schrappe M, Valsecchi MG, Bartram CR, Schrauder A, Panzer-Grumayer R, Moricke A, et al. Late MRD response determines relapse risk overall and in subsets of childhood T-cell ALL: results of the AIEOP-BFM-ALL 2000 study. Blood. 2011;118(8):2077–84.
Moricke A, Zimmermann M, Valsecchi MG, Stanulla M, Biondi A, Mann G, et al. Dexamethasone vs prednisone in induction treatment of pediatric ALL: results of the randomized trial AIEOP-BFM ALL 2000. Blood. 2016;127(17):2101–12.
Ellinghaus E, Stanulla M, Richter G, Ellinghaus D, te Kronnie G, Cario G, et al. Identification of germline susceptibility loci in ETV6-RUNX1-rearranged childhood acute lymphoblastic leukemia. Leukemia. 2012;26(5):902–9.
Strassburg CP. Hyperbilirubinemia syndromes (Gilbert-Meulengracht, Crigler-Najjar, Dubin-Johnson, and rotor syndrome). Best Pract Res Clin Gastroenterol. 2010;24(5):555–71.
Strassburg CP. Gilbert-Meulengracht's syndrome and pharmacogenetics: is jaundice just the tip of the iceberg? Drug Metab Rev. 2010;42(1):168–81.
Ehmer U, Kalthoff S, Fakundiny B, Pabst B, Freiberg N, Naumann R, et al. Gilbert syndrome redefined: a complex genetic haplotype influences the regulation of glucuronidation. Hepatology. 2012;55(6):1912–21.
Strassburg CP. Pharmacogenetics of Gilbert's syndrome. Pharmacogenomics. 2008;9(6):703–15.
Beutler E, Gelbart T, Demina A. Racial variability in the UDP-glucuronosyltransferase 1 (UGT1A1) promoter: a balanced polymorphism for regulation of bilirubin metabolism? Proc Natl Acad Sci U S A. 1998;95(14):8170–4.
Schmiegelow K, Attarbaschi A, Barzilai S, Escherich G, Frandsen TL, Halsey C, et al. Consensus definitions of 14 severe acute toxic effects for childhood lymphoblastic leukaemia treatment: a Delphi consensus. Lancet Oncol. 2016;17(6):e231–e9.
Relling MV, Schwab M, Whirl-Carrillo M, Suarez-Kurtz G, Pui CH, Stein CM, et al. Clinical Pharmacogenetics implementation consortium guideline for Thiopurine dosing based on TPMT and NUDT15 genotypes: 2018 update. Clin Pharmacol Ther. 2019;105(5):1095–105.
Nelson RS, Seligson ND, Bottiglieri S, Carballido E, Cueto AD, Imanirad I, et al. UGT1A1 guided Ccancer therapy: review of the evidence and considerations for clinical implementation. Cancers. 2021;13(7):1566.
Recommendations from the EGAPP Working Group. Can UGT1A1 genotyping reduce morbidity and mortality in patients with metastatic colorectal cancer treated with irinotecan? Genet Med. 2009;11(1):15–20.
Mackenzie PI, Owens IS, Burchell B, Bock KW, Bairoch A, Belanger A, et al. The UDP glycosyltransferase gene superfamily: recommended nomenclature update based on evolutionary divergence. Pharmacogenet Genomics. 1997;7(4):255–69.
Zahreddine HA, Culjkovic-Kraljacic B, Gasiorek J, Duchaine J, Borden KLB. GLI1-inducible Glucuronidation targets a broad Spectrum of drugs. ACS Chem Biol. 2019;14(3):348–55.
Gessner T, Vaughan LA, Beehler BC, Bartels CJ, Baker RM. Elevated pentose cycle and glucuronyltransferase in daunorubicin-resistant P388 cells. Cancer Res. 1990;50(13):3921–7.
Weenen H, van Maanen JM, de Planque MM, McVie JG, Pinedo HM. Metabolism of 4′-modified analogs of doxorubicin. Unique glucuronidation pathway for 4′-epidoxorubicin. Eur J Cancer Clin Oncol. 1984;20(7):919–26.
Tukey RH, Strassburg CP. Human UDP-glucuronosyltransferases: metabolism, expression, and disease. Annu Rev Pharmacol Toxicol. 2000;40:581–616.
Watanabe Y, Nakajima M, Ohashi N, Kume T, Yokoi T. Glucuronidation of etoposide in human liver microsomes is specifically catalyzed by UDP-glucuronosyltransferase 1A1. Drug Metab Dispos. 2003;31(5):589–95.
Guillemette C, Lévesque É, Rouleau M. Pharmacogenomics of human uridine diphospho-glucuronosyltransferases and clinical implications. Clin Pharmacol Ther. 2014;96(3):324–39.
Ha VH, Jupp J, Tsang RY. Oncology drug dosing in Gilbert syndrome associated with UGT1A1: a summary of the literature. Pharmacotherapy. 2017;37(8):956–72.
Allain EP, Rouleau M, Levesque E, Guillemette C. Emerging roles for UDP-glucuronosyltransferases in drug resistance and cancer progression. Br J Cancer. 2020;122(9):1277–87.
Ritter JK, Crawford JM, Owens IS. Cloning of two human liver bilirubin UDP-glucuronosyltransferase cDNAs with expression in COS-1 cells. J Biol Chem. 1991;266(2):1043–7.
Bosma PJ. Inherited disorders of bilirubin metabolism. J Hepatol. 2003;38(1):107–17.
Bosma PJ, Chowdhury JR, Bakker C, Gantla S, de Boer A, Oostra BA, et al. The genetic basis of the reduced expression of bilirubin UDP-glucuronosyltransferase 1 in Gilbert's syndrome. N Engl J Med. 1995;333(18):1171–5.
Black M, Billing BH. Hepatic bilirubin udp-glucuronyl transferase activity in liver disease and gilbert's syndrome. N Engl J Med. 1969;280(23):1266–71.
Jansen PL, Bosma PJ, Chowdhury JR. Molecular biology of bilirubin metabolism. Prog Liver Dis. 1995;13:125–50.
Kim JY, Cheong HS, Park BL, Kim LH, Namgoong S, Kim JO, et al. Comprehensive variant screening of the UGT gene family. Yonsei Med J. 2014;55(1):232–9.
Lankisch TO, Behrens G, Ehmer U, Mobius U, Rockstroh J, Wehmeier M, et al. Gilbert's syndrome and hyperbilirubinemia in protease inhibitor therapy--an extended haplotype of genetic variants increases risk in indinavir treatment. J Hepatol. 2009;50(5):1010–8.
Hoskins JM, Goldberg RM, Qu P, Ibrahim JG, McLeod HL. UGT1A1*28 genotype and irinotecan-induced neutropenia: dose matters. J Natl Cancer Inst. 2007;99(17):1290–5.
Goey AK, Sissung TM, Peer CJ, Trepel JB, Lee MJ, Tomita Y, et al. Effects of UGT1A1 genotype on the pharmacokinetics, pharmacodynamics, and toxicities of belinostat administered by 48-hour continuous infusion in patients with cancer. J Clin Pharmacol. 2016;56(4):461–73.
Peer CJ, Goey AKL, Sissung TM, Erlich S, Lee M-J, Tomita Y, et al. UGT1A1 genotype-dependent dose adjustment of belinostat in patients with advanced cancers using population pharmacokinetic modeling and simulation. J Clin Pharmacol. 2016;56(4):450–60.
Zahreddine HA, Culjkovic-Kraljacic B, Assouline S, Gendron P, Romeo AA, Morris SJ, et al. The sonic hedgehog factor GLI1 imparts drug resistance through inducible glucuronidation. Nature. 2014;511(7507):90–3.
Kishi S, Yang W, Boureau B, Morand S, Das S, Chen P, et al. Effects of prednisone and genetic polymorphisms on etoposide disposition in children with acute lymphoblastic leukemia. Blood. 2004;103(1):67–72.
Berrueco R, Alonso-Saladrigues A, Martorell-Sampol L, Catala-Temprano A, Ruiz-Llobet A, Toll T, et al. Outcome and toxicities associated to chemotherapy in children with acute lymphoblastic leukemia and Gilbert syndrome. Usefulness of UGT1A1 mutational screening. Pediatr Blood Cancer. 2015;62(7):1195–201.
Nomura A, Maruo Y, Taga T, Takeuchi Y. Contribution of UGT1A1 variations to chemotherapy-induced unconjugated hyperbilirubinemia in pediatric leukemia patients. Pediatr Res. 2016;80(2):252–7.
We are indebted to all patients, parents, nurses, technicians, data managers and doctors who participated in or contributed to AIEOP-BFM ALL 2000.
Open Access funding enabled and organized by Projekt DEAL. TRANSCALL2, ERA-NET TRANSCAN/European Commission under the 7th Framework Programme (FP7), Madeleine-Schickedanz-Kinderkrebsstiftung, Deutsche Krebshilfe, Verein für krebskranke Kinder Hannover e.V., Deutsche José Carreras Leukämie-Stiftung, Wilhelm Sander-Stiftung, Titus Clinician Scientist Program of the Else Kröner-Fresenius-Stiftung. M. Schwab and E. Schaeffeler are funded by the Robert Bosch Stiftung Stuttgart, the German Cancer Consortium (DKTK), the German Cancer Research Center (DKFZ) Partner Site Tübingen, and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2180 – 390900677.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board of Hannover Medical School, Hannover, Germany, and informed consent was obtained from patients and/or their guardians.
None of the authors has a conflict of interest regarding the present study. The data presented here have not been previously presented and have not previously been submitted for publication in any other journal.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Suppl. Table 1.
Treatment details of protocol AIEOP-BFM ALL 2000. Suppl. Table 2. Clinical characteristics of the patients in the study cohort by severity of bilirubin toxicity during induction/consolidation (protocols IA/IB, n = 1547). Suppl. Table 3. Characteristics of the patients in the GWAS discovery cohort by serum bilirubin levels during induction/consolidation (n = 650). Suppl. Table 4. Characteristics of the patients in the GWAS discovery cohort compared to all patients of the study cohort with toxicity information. Suppl. Table 5. Summary of genome-wide association analysis for therapy-related hyperbilirubinemia during induction/consolidation (protocols IA/IB). Suppl. Table 6. Allelic association of hyperbilirubinemia phenotype with the 20 most strongly associated variants around rs6744284 resulting from genotype imputation. Suppl. Table 7. Genotypic association between hyperbilirubinemia phenotype and the 20 most strongly associated SNV around rs6744284 after genotype imputation. Suppl. Table 8. Adjusted genotypic association of rs6744284 with hyperbilirubinemia during protocols IA/IB and later therapeutic elements, including age and immunophenotype as covariates. Suppl. Table 9. Genotypic association of rs6744284 with hyperbilirubinemia phenotype stratified for potential effect modifiers. Suppl. Table 10. Characteristics of acute lymphoblastic leukemia (ALL) patients included in the replication cohort (n = 224). Suppl. Table 11. Characteristics of acute lymphoblastic leukemia (ALL) patients included in subsequent UGT1A1*28/*37 genotyping (n = 544). Suppl. Table 12. Association of identified risk loci and known Gilbert’s syndrome related variants with hyperbilirubinemia. Suppl. Table 13. Correlation of rs6744284 with imputed top SNV and UGT1A variations related to hyperbilirubinemia and the Gilbert’s syndrome. Suppl. information on the correlation analysis of known Gilbert’s syndrome (GS) related variations with rs6744284. Suppl. information on the impact of hyperbilirubinemia on therapy delays in the discovery cohort. Suppl. Table 14. Clinical characteristics of the acute lymphoblastic leukemia (ALL) patients of the discovery cohort according to the severity of bilirubin toxicity during induction/consolidation (protocols IA/IB, n = 650). Suppl. Table 15. Estimated hazard ratios from the multivariable Cox proportional model on the hazard of relapse in patients with high hyperbilirubinemia (≥CTC grade 3) during induction and/or consolidation (n = 68). Suppl. Fig. 1. Consolidated Standards of Reporting Trials (CONSORT) diagram of inclusion criteria for the study population (N = 1547). Suppl. Fig. 2. Estimated 5-year event-free survival (EFS) and cumulative incidence of relapse (CIR) at 5 years in the study cohort by the maximum transaminase levels during protocol IA/IB [%]. Suppl. Fig. 3. Regional plot of association results and recombination rates for the identified risk locus in the UGT1A region (2q37). Suppl. Fig. 4. Total serum bilirubin levels by treatment element and rs6744284 genotype. Suppl. Fig. 5. Estimated 5-year event-free survival (EFS) and cumulative incidence of relapse (CIR) at 5 years in the discovery cohort by the maximum transaminase levels during protocol IA/IB [%].
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.
About this article
Cite this article
Junk, S.V., Schaeffeler, E., Zimmermann, M. et al. Chemotherapy-related hyperbilirubinemia in pediatric acute lymphoblastic leukemia: a genome-wide association study from the AIEOP-BFM ALL study group. J Exp Clin Cancer Res 42, 21 (2023). https://doi.org/10.1186/s13046-022-02585-x
- Childhood acute lymphoblastic leukemia
- Treatment-related toxicity