Influence of glutathione S-transferase gene polymorphisms on busulfan pharmacokinetics and outcome of hematopoietic stem-cell transplantation in thalassemia pediatric patients

Haplotype Inference Using Long-Read Nanopore Sequencing: Application to GSTA1 Promoter

Recovering true haplotypes can have important clinical consequences. The laboratory process is difficult and is, therefore, most often done through inference. In this paper, we show that when using the Oxford nanopore sequencing technology, we could recover the true haplotypes of the GSTA1 promoter region. Eight LCL cell lines with potentially ambiguous haplotypes were used to characterize the efficacy of Oxford nanopore sequencing to phase the correct GSTA1 promoter haplotypes. The results were compared to Sanger sequencing and inferred haplotypes in the 1000 genomes project. The average read length was 813 bp out of a total PCR length of 1336 bp. The best coverage of sequencing was in the middle of the PCR product and decreased to 50% at the PCR ends. SNPs separated by less than 200 bp showed > 90% of correct haplotypes, while at the distance of 1089 bp, this proportion still exceeded 58%. The number of cycles influences the generation of hybrid haplotypes but not extension or annealing time. The results demonstrate that this long sequencing reads methodology, can accurately determine the haplotypes without the need for inference. The technology proved to be robust but the success of phasing nonetheless depends on the distances and frequencies of SNPs.

Pharmacokinetic Modeling and Simulation with Pharmacogenetic Insights Support the Relevance of Therapeutic Drug Monitoring for Myeloablative Busulfan Dosing in Adult HSCT

Therapeutic drug monitoring (TDM) of busulfan (Bu) is well-established in pediatric hematopoietic stem cell transplantation (HSCT), but its use in adults is limited due to a lack of clear recommendations and scarcity of evidence regarding its utility. GSTA1 promoter variants are reported to affect Bu clearance in both adults and pediatric patients. This study aimed to evaluate the value of preemptive genotyping GSTA1 and body composition (obesity) in individualizing Bu dosing in adults, through pharmacokinetic (PK) modeling and simulations. A population pharmacokinetic (PopPK) model was developed and validated with data from 60 adults who underwent HSCT. Simulations assessed different dosing scenarios based on body size metrics and GSTA1 genotypes. Due to the limited number of obese patients in the cohort, the effect of obesity on Bu pharmacokinetics (PK) was evaluated in silico using a physiologically-based pharmacokinetic (PBPK) model and relevant virtual populations from Simcyp software. Patients with at least 1 GSTA1*B haplotype had 17% lower clearance on average. PopPK simulations indicated that adjusting doses based on genotype increased the probability of achieving the target exposure (3.7 to 5.5 mg.h/L) from 53% to 60 % in GSTA1*A homozygous patients, and from 50% to 61% in *B carriers. Still, Approximately 40% of patients would not achieve this therapeutic window without TDM. A 2-sample optimal design was validated for routine model-based Bu first dose AUC0-∞ estimation, and the model was implemented in the Tucuxi user-friendly TDM software. PBPK simulations confirmed body surface area-based doses of 29 to 31 mg/m2/6h as the most appropriate, regardless of obesity status. This study emphasizes the importance of individualized Bu dosing strategies in adults to achieve therapeutic targets. Preemptive genotyping alone may not have a significant clinical impact, and routine TDM may be necessary for optimal transplantation outcomes.

Treosulfan Exposure Predicts Thalassemia-Free Survival in Patients with Beta Thalassemia Major Undergoing Allogeneic Hematopoietic Cell Transplantation

A toxicity-reduced conditioning regimen with treosulfan, fludarabine, and thiotepa in patients with high-risk β-thalassemia major has significantly improved hematopoietic stem cell transplantation (HCT) outcomes. However, complications resulting from regimen-related toxicities (RRTs), mixed chimerism, and graft rejection remain a challenge. We evaluated the dose-exposure-response relationship of treosulfan and its active metabolite S, S-EBDM, in a uniform cohort of patients with β-thalassemia major to identify whether therapeutic drug monitoring (TDM) and dose adjustment of treosulfan is feasible. Plasma treosulfan/S, S-EBDM levels were measured in 77 patients using a validated liquid chromatography with tandem mass spectrometry method, and the pharmacokinetic parameters were estimated using nlmixr2. The influence of treosulfan and S, S-EBDM exposure, and GSTA1/NQO1 polymorphisms on graft rejection, RRTs, chimerism status, and 1-year overall survival (OS), and thalassemia-free survival (TFS) were assessed. We observed that treosulfan exposure was lower in patients with graft rejection than those without (1,655 vs. 2,037 mg•h/L, P = 0.07). Pharmacodynamic modeling analysis to identify therapeutic cutoff revealed that treosulfan exposure ≥1,660 mg•hour/L was significantly associated with better 1-year TFS (97% vs. 81%, P = 0.02) and a trend to better 1-year OS (90% vs. 69%, P = 0.07). Further, multivariate analysis adjusting for known pre-HCT risk factors also revealed treosulfan exposure <1,660 mg•h/L (hazard ratio (HR) = 3.23; 95% confidence interval (CI) = 1.12-9.34; P = 0.03) and GSTA1*B variant genotype (HR = 3.75; 95% CI = 1.04-13.47; P = 0.04) to be independent predictors for inferior 1-year TFS. We conclude that lower treosulfan exposure increases the risk of graft rejection and early transplant-related mortality affecting TFS. As no RRTs were observed with increasing treosulfan exposure, TDM-based dose adjustment could be feasible and beneficial.

Population Pharmacokinetic Modeling for Twice-Daily Intravenous Busulfan in a Large Cohort of Pediatric Patients Undergoing Hematopoietic Stem Cell Transplantation-A 10-Year Single-Center Experience

Reaching target exposure of busulfan-based conditioning prior to hematopoietic stem cell transplantation is vital for favorable therapy outcomes. Yet, a wide inter-patient and inter-occasion variability in busulfan exposure has been reported, especially in children. We aimed to identify factors associated with the variability of busulfan pharmacokinetics in 124 consecutive patients transplanted at the University Children's Hospital Zurich between October 2010 and February 2020. Clinical data and busulfan plasma levels after twice-daily intravenous administration were analyzed retrospectively by population pharmacokinetic modeling. The volume of distribution correlated with total body water. The elimination rate constant followed an age-dependent maturation function, as previously suggested, and correlated with the levels of serum albumin. Acute lymphoblastic leukemia reduced busulfan clearance by 20%. Clearance significantly decreased by 17% on average from the start to the third day of busulfan administration, in agreement with other studies. An average reduction of 31% was found in patients with hemophagocytic lymphohistiocytosis and X-linked lymphoproliferative disease. In conclusion, we demonstrate that in addition to known factors, underlying disease and serum albumin significantly impact busulfan pharmacokinetics in pediatric patients; yet, substantial unexplained variability in some patients remained. Thus, we consider repeated pharmacokinetic assessment essential to achieve the desired target exposure in twice-daily busulfan administration.

Pharmacogenetic Aspects of Drug Metabolizing Enzymes and Transporters in Pediatric Medicine: Study Progress, Clinical Practice and Future Perspectives

As the activity of certain drug metabolizing enzymes or transporter proteins can vary with age, the effect of ontogenetic and genetic variation on the activity of these enzymes is critical for the accurate prediction of treatment outcomes and toxicity in children. This makes pharmacogenetic research in pediatrics particularly important and urgently needed, but also challenging. This review summarizes pharmacogenetic studies on the effects of genetic polymorphisms on pharmacokinetic parameters and clinical outcomes in pediatric populations for certain drugs, which are commonly prescribed by clinicians across multiple therapeutic areas in a general hospital, organized from those with the most to the least pediatric evidence among each drug category. We also further discuss the research status of the gene-guided dosing regimens and clinical implementation of pediatric pharmacogenetics. More and more drug-gene interactions are demonstrated to have clinical validity for children, and pharmacogenomics in pediatrics have shown evidence-based benefits to enhance the efficacy and precision of existing drug dosing regimens in several therapeutic areas. However, the most important limitation to the implementation is the lack of high-quality, rigorous pediatric prospective clinical studies, so adequately powered interventional clinical trials that support incorporation of pharmacogenetics into the care of children are still needed.

Impact of Glutathione S-Transferase Polymorphisms on Busulfan Pharmacokinetics and Outcomes of Hematopoietic Stem Cell Transplantation

BACKGROUND

Busulfan (Bu) is an alkylating drug used in many preparative regimens before hematopoietic stem cell transplantation (HSCT). It is conjugated in the liver mainly by glutathione S-transferase isoenzyme A1-1 ( GSTA1 ). Genetic polymorphisms in these isoenzymes may affect the pharmacokinetics of Bu and the clinical outcomes of HSCT. This study aimed to assess the impact of glutathione S-transferase ( GST ) genetic polymorphisms on the clearance of Bu and the clinical outcomes of patients undergoing HSCT.

METHODS

This single-center retrospective study included patients who received IV Bu before HSCT at Sultan Qaboos University Hospital (SQUH), Oman from January 2003 to October 2016. Genotyping for polymorphisms was performed for GSTM1 , GSTT1 , GSTA1 , and GSTP1 . Each GST polymorphism was analyzed for its impact on Bu clearance and HSCT outcomes.

RESULTS

A total of 135 patients were included. The mean Bu clearance was 3.7 ± 0.98 mL/min/kg. Patients with GSTA1 A-513G heterozygosity (AG) were found to have a higher incidence of graft loss ( P = 0.006). Homozygous double null of GSTM1 and GSTT1 was associated with a higher incidence of acute graft versus host disease ( P = 0.04). Double non-null GSTM1 and GSTT1 and non-null GSTM1 increased the risk of mortality ( P = 0.034 and 0.021, respectively).

CONCLUSIONS

GST genotyping before HSCT may predict HSCT outcomes. The results of this preliminary retrospective study need to be confirmed in a larger prospective study.

External Evaluation of Population Pharmacokinetic Models of Busulfan in Chinese Adult Hematopoietic Stem Cell Transplantation Recipients

Objective: Busulfan (BU) is a bi-functional DNA-alkylating agent used in patients undergoing hematopoietic stem cell transplantation (HSCT). Over the last decades, several population pharmacokinetic (pop PK) models of BU have been established, but external evaluation has not been performed for almost all models. The purpose of the study was to evaluate the predictive performance of published pop PK models of intravenous BU in adults using an independent dataset from Chinese HSCT patients, and to identify the best model to guide personalized dosing.

Methods: The external evaluation methods included prediction-based diagnostics, simulation-based diagnostics, and Bayesian forecasting. In prediction-based diagnostics, the relative prediction error (PE%) was calculated by comparing the population predicted concentration (PRED) with the observations. Simulation-based diagnostics included the prediction- and variability-corrected visual predictive check (pvcVPC) and the normalized prediction distribution error (NPDE). Bayesian forecasting was executed by giving prior one to four observations. The factors influencing the model predictability, including the impact of structural models, were assessed.

Results: A total of 440 concentrations (110 patients) were obtained for analysis. Based on prediction-based diagnostics and Bayesian forecasting, preferable predictive performance was observed in the model developed by Huang et al. The median PE% was -1.44% which was closest to 0, and the maximum F₂₀ of 57.27% and F₃₀ of 72.73% were achieved. Bayesian forecasting demonstrated that prior concentrations remarkably improved the prediction precision and accuracy of all models, even with only one prior concentration.

Conclusion: This is the first study to comprehensively evaluate published pop PK models of BU. The model built by Huang et al. had satisfactory predictive performance, which can be used to guide individualized dosage adjustment of BU in Chinese patients.

Impact of busulfan pharmacokinetics on outcome in adult patients receiving an allogeneic hematopoietic cell transplantation

Busulfan (Bu) is widely used in conditioning regimens before allogeneic hematopoietic cell transplantation, with variable metabolism due to interindividual differences of pharmacokinetics (PK). The purpose of this study was to correlate pharmacokinetics and clinical outcomes. Lower-AUC, in range-AUC and higher-AUC were defined as ±25% of the targeted Bu-AUC. In 2019, we changed Bu dosing from 4×/day (Bu-4) to 1×/day (Bu-1) for ease of application. AUC-target range was reached in 46% of patients; 40% were in low-AUC and 14% in high-AUC. Among all toxicities, viral and fungal infections were significantly more frequent in high-AUC compared with low-AUC (20% vs. 8%; p = 0.01 and 37% vs. 17%; p = 0.03). Bu-1 showed lower PK values (66% vs. 36% of Bu-4 in low-AUC; p < 0.01) and higher incidence of mucositis (p = 0.02). Long-term outcomes at 2 years showed a higher non-relapse mortality (NRM) (p < 0.01) and higher relative risk of death in the high-AUC group compared to the other groups. Cumulative incidence of relapse and acute/chronic GvHD were not significantly different. The optimal cut-off in Bu-AUC associated with low NRM was 969 µmol/l*min (ROC AUC 0.67, sensitivity 0.86 and specificity 0.47) for Bu-4. In conclusion, low-AUC BU-PK seems of benefit regarding NRM and survival.

A potential implication of UDP-glucuronosyltransferase 2B10 in the detoxification of drugs used in pediatric hematopoietic stem cell transplantation setting: an in silico investigation

Background

Sinusoidal occlusion syndrome (SOS) is a potentially severe complication following hematopoietic stem cell transplantation (HSCT) in pediatric patients. Treatment related risk factors such as intensity of conditioning, hepatotoxic co-medication and patient related factors such as genetic variants predispose individuals to develop SOS. The variant allele for SNP rs17146905 in UDP-glucuronosyl transferase 2B10 (UGT2B10) gene was correlated with the occurrence of SOS in an exome-wide association study. UGT2B10 is a phase II drug metabolizing enzyme involved in the N-glucuronidation of tertiary amine containing drugs.

Methods

To shed light on the functionality of UGT2B10 enzyme in the metabolism of drugs used in pediatric HSCT setting, we performed in silico screening against custom based library of putative ligands. First, a list of potential substrates for in silico analysis was prepared using a systematic consensus-based strategy. The list comprised of drugs and their metabolites used in pediatric HSCT setting. The three-dimensional structure of UGT2B10 was not available from the Research Collaboratory Structural Bioinformatics - Protein Data Bank (RCSB - PDB) repository and thus we predicted the first human UGT2B10 3D model by using multiple template homology modeling with MODELLER Version 9.2 and molecular docking calculations with AutoDock Vina Version 1.2 were implemented to quantify the estimated binding affinity between selected putative substrates or ligands and UGT2B10. Finally, we performed molecular dynamics simulations using GROMACS Version 5.1.4 to confirm the potential UGT2B10 ligands prioritized after molecular docking (exhibiting negative free binding energy).

Results

Four potential ligands for UGT2B10 namely acetaminophen, lorazepam, mycophenolic acid and voriconazole n-oxide intermediate were identified. Other metabolites of voriconazole satisfied the criteria of being possible ligands of UGT2B10. Except for bilirubin and 4-Hydroxy Voriconazole, all the ligands (particularly voriconazole and hydroxy voriconazole) are oriented in substrate binding site close to the co-factor UDP (mean ± SD; 0.72 ± 0.33 nm). Further in vitro screening of the putative ligands prioritized by in silico pipeline is warranted to understand the nature of the ligands either as inhibitors or substrates of UGT2B10.

Conclusions

These results may indicate the clinical and pharmacological relevance UGT2B10 in pediatric HSCT setting. With this systematic computational methodology, we provide a rational-, time-, and cost-effective way to identify and prioritize the interesting putative substrates or inhibitors of UGT2B10 for further testing in in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00402-5.

Placental transcriptome sequencing combined with bioinformatics predicts potential genes and circular RNAs associated with hemoglobin Bart's hydrops fetalis syndrome

AIM

Hemoglobin Bart's hydrops fetalis syndrome (BHFS) is the most severe form of α-thalassemia. Histological alternations can be observed in placenta, but placental transcriptome profile and circular RNAs have not been studied in this disease. The aim of this study was to define the placental transcriptional changes and find relevant circular RNAs in BHFS.

METHODS

We performed high-throughput RNA sequencing to detect placental samples from fetuses affected by BHFS (n = 5) and normal fetuses (NF, n = 5), quantitative reverse transcription polymerase chain reaction (RT-qPCR), and Sanger sequencing to validate the differentially expressed circRNAs and their potentially related miRNAs (BHFS, n = 22; NF, n = 11). Bioinformatics methods were performed for further analysis.

RESULTS

Our results showed 152 differentially expressed genes (DEGs), 112 circRNAs, and 45 microRNAs that were differentially expressed. DEGs were found to be involved in Gene Ontology terms related to gas transport, cell adhesion, oxidative stress, organ development, hemopoiesis, and others. RT-qPCR results showed that hsa_circ_0003961 and hsa_circ_0006687 were upregulated (p < 0.05). The competing endogenous RNA and co-expression networks showed that hsa_circ_0003961 and hsa_circ_0006687 were connected with 3 miRNAs and some DEGs, including cell adhesion genes (e.g., CLDN19), hemoglobin related genes (e.g., SOX6 and HBZ) and angiogenesis related genes (e.g., EPHB2). Downregulations of hsa-miR-1299 and hsa-miR-625-5p in ceRNA network were also validated by RT-qPCR. Gene set enrichment analysis results for the two circRNAs showed that some gene sets associated with cell adhesion, hematopoietic system and apoptosis were significantly enriched.

CONCLUSIONS

Our study characterized the placental transcriptome of BHFS. The circRNAs hsa_circ_0003961 and hsa_circ_0006687 in placenta may be relevant to BHFS.

Influence of GSTM1, GSTT1, and GSTP1 genetic polymorphisms on disorders in transplant patients: a systematic review

The glutathione-S-transferase (GST) enzymes are phase II isoenzymes responsible for protection against free radicals and xenobiotics. Since these proteins are described as polymorphic, polymorphisms in genes that encode them may alter enzymatic function and contribute to oxidative stress. In this context, such polymorphisms were already associated with several diseases and multiple therapeutic outcomes. A systematic review was performed to evaluate studies regarding the association between polymorphisms in three genes encoding enzymes of the GST family - GSTM1, GSTT1, and GSTP1 - and disorders in transplant patients. A total of 125 articles on which inclusion and exclusion criteria were applied were identified at PubMed database. Thirty-two studies met the target criteria and were included in the review. The mechanisms by which GST genotypes influence the development of disorders in transplant patients differ by disorder: they may participate in it by decreasing metabolism of drugs administered to patients undergoing transplantation, then exposing them to greater toxicity; by decreasing the repair ability against oxidative stress; or by encoding proteins that may be recognized as foreign, setting of an alloimmune reaction. Although some results are better established - such as GSTM1 null genotype's role in the development of toxicity events in transplant patients - others require further evidences, as GST influence on the development of pulmonary decline and posttransplant diabetes mellitus (PTDM). The importance of investigating these associations lies in a personalized medicine, in which the high-risk genotype patient has its treatment individualized and its care for prophylaxis and surveillance increased, potentially reducing this population's morbimortality.

Is Busulfan Clearance Different in Patients With Sickle Cell Disease? Let's Clear Up That Case With Some Controls

In busulfan-based conditioning regimen for hematopoietic stem cell transplantation in children, accurate a priori determination of the first dose is important because of its narrow therapeutic window. Sickle cell disease (SCD) influences pharmacokinetics of the commonly used drugs by affecting organs responsible for drug metabolism and elimination. This pharmacokinetics study assesses the influence of SCD on the metabolic pathway of busulfan that is mainly metabolized in the liver. In this retrospective cross-sectional case-control study, 16 patients with SCD were matched to 50 patients without SCD on known busulfan clearance's covariates (glutathione-S-transferase alpha1 polymorphisms, age, weight). Clearance of the first dose of busulfan was not significantly different independently of genetic or anthropometric factors in patients with or without SCD.

Neuroprotective Efficacy of Prophylactic Antiepileptic Therapies in Busulfan Conditioning Regimen

Busulfan (BU) is commonly used in the anticancer treatment and bone marrow/hematopoietic stem cell transplantation (BM/HSCT) conditioning regimen. This study aimed to evaluate the neuroprotective efficacy of antiepileptic therapies in the BU conditioning regimen. Patients who underwent BM/HSCT between January 1, 2009, and January 1, 2014, were evaluated retrospectively. Diagnosis, demographic characteristics, neurological side effects, electroencephalography (EEG) abnormalities, and antiepileptic drugs were determined. Levetiracetam was given to 49% and phenytoin to 23% of patients as antiepileptic prophylaxis. After BU treatment, seizures developed in three patients, posterior reversible leukoencephalopathy in two patients, and peripheral neuropathy in two patients. Forty-five patients who underwent EEG before and after BU treatment were identified. The effects of antiepileptic prophylaxis over the electroencephalographic abnormalities after the BU conditioning regimen were found to be not different (p = 0.491). Neurological side effects and electroencephalographic abnormalities following BU conditioning regimen under antiepileptic prophylaxis are comparable.

Hematopoietic stem cell transplantation for thalassemia major using HLA fully-matched and mismatched donor grafts

BACKGROUND

Until now, allogeneic hematopoietic stem cell transplantation (allo-HSCT) is the only effective method to cure Thalassemia major. However, it has not been determined whether similar results can be obtained with the same conditioning regimen for both fully-matched and mismatched donors grafts. We hypothesized that using modified NF-08-TM conditioning regimen could achieve similar results for both fully and mismatched donors grafts.

METHODS

This retrospective cohort study included patients with β-thalassemia major who underwent HSCT with modified NF-08-TM conditioning regimen at Guangzhou Women and Children's Medical Centre between January 2013 and January 2019.

RESULTS

Among the 257 patients (172 males) included in this study, 3 had two transplantations. Totally 193 and 67 had fully-matched and mismatched donors were examined, respectively. The median follow-up was 29 months; 6-year overall survival (OS), thalassemia-free survival (TFS), graft rejection (GR) and transplantation-related mortality (TRM) were 92.08%, 90.89%, 1.24% and 8.01%, respectively. Multivariate analysis showed that human leukocyte antigen (HLA) compatibility between patient and donor was not independently associated with OS, TFS, GR or TRM. Mismatched donor graft transplantation for β-thalassemia major is associated with similar survival outcomes and incidences of complications (except for acute GVHD) to fully-matched donor graft transplantation based on modified NF-08-TM conditioning regimen.

CONCLUSIONS

In conclusion, based on the modified NF-08-TM regimen, certain mismatched donor transplantations for thalassemia major patients could achieve similar results as fully-matched donor transplantations.

Genetic Predictors for Sinusoidal Obstruction Syndrome-A Systematic Review

Sinusoidal obstruction syndrome (SOS) is a potentially life-threatening complication after hematopoietic stem cell transplantation (HSCT) or antineoplastic treatment without HSCT. Genetic variants were investigated for their association with SOS, but the evidence is inconclusive. We performed a systematic literature review to identify genes, gene variants, and methods of association analyses of genetic markers with SOS. We identified 23 studies after HSCT and 4 studies after antineoplastic treatment without HSCT. One study (4%) performed whole-exome sequencing (WES) and replicated the analysis in an independent cohort, 26 used a candidate-gene approach. Three studies included >200 participants (11%), and six were of high quality (22%). Variants in 34 genes were tested in candidate gene studies after HSCT. Variants in GSTA1 were associated with SOS in three studies, MTHFR in two, and CPS1, CTH, CYP2B6, GSTM1, GSTP1, HFE, and HPSE in one study each. UGT2B10 and LNPK variants were identified in a WES analysis. After exposure to antineoplastic agents without HSCT, variants in six genes were tested and only GSTM1 was associated with SOS. There was a substantial heterogeneity of populations within and between studies. Future research should be based on sufficiently large homogenous samples, adjust for covariates, and replicate findings in independent cohorts.

The analysis of GSTA1 promoter genetic and functional diversity of human populations

GSTA1 encodes a member of a family of enzymes that function to add glutathione to target electrophilic compounds, including carcinogens, therapeutic drugs, environmental toxins, and products of oxidative stress. GSTA1 has several functional SNPs within its promoter region that are responsible for a change in its expression by altering promoter function. This study aims to investigate distributions of GSTA1 promoter haplotypes across different human populations and to assess their impact on the expression of GSTA1. PHASE 2.1.1 was used to infer haplotypes and diplotypes of six GSTA1 promoter SNPs on 2501 individuals from 26 populations classified by the 1000 Genomes Project into five super-populations that included Africa (N = 660), America (N = 347), East Asia (N = 504), Europe (N = 502), and South Asia (N = 488). We used pairwise FST analysis to compare sub-populations and luciferase reporter assay (LRA) to evaluate the impact of each SNP on activation of transcription and interaction with other SNPs. The distributions of GSTA1 promoter haplotypes and diplotypes were significantly different among the different human populations. Three new promoter haplotypes were found in the African super-population. LRA demonstrated that SNPs at -52 and -69 has the most impact on GSTA1 expression, however other SNPs have a significant impact on transcriptional activity. Based on LRA, a new model of cis-elements interaction is presented. Due to the significant differences in GSTA1 diplotype population frequencies, future pharmacogenomics or disease-related studies would benefit from the inclusion of the complete GSTA1 promoter haplotype based on the newly proposed metabolic grouping derived from the LRA results.

Total Body Irradiation Forever? Optimising Chemotherapeutic Options for Irradiation-Free Conditioning for Paediatric Acute Lymphoblastic Leukaemia

Total-body irradiation (TBI) based conditioning prior to allogeneic hematopoietic stem cell transplantation (HSCT) is generally regarded as the gold-standard for children >4 years of age with acute lymphoblastic leukaemia (ALL). Retrospective studies in the 1990's suggested better survival with irradiation, confirmed in a small randomised, prospective study in the early 2000's. Most recently, this was reconfirmed by the early results of the large, randomised, international, phase III FORUM study published in 2020. But we know survivors will suffer a multitude of long-term sequelae after TBI, including second malignancies, neurocognitive, endocrine and cardiometabolic effects. The drive to avoid TBI directs us to continue optimising irradiation-free, myeloablative conditioning. In chemotherapy-based conditioning, the dominant myeloablative effect is provided by the alkylating agents, most commonly busulfan or treosulfan. Busulfan with cyclophosphamide is a long-established alternative to TBI-based conditioning in ALL patients. Substituting fludarabine for cyclophosphamide reduces toxicity, but may not be as effective, prompting the addition of a third agent, such as thiotepa, melphalan, and now clofarabine. For busulfan, it's wide pharmacokinetic (PK) variability and narrow therapeutic window is well-known, with widespread use of therapeutic drug monitoring (TDM) to individualise dosing and control the cumulative busulfan exposure. The development of first-dose selection algorithms has helped achieve early, accurate busulfan levels within the targeted therapeutic window. In the future, predictive genetic variants, associated with differing busulfan exposures and toxicities, could be employed to further tailor individualised busulfan-based conditioning for ALL patients. Treosulfan-based conditioning leads to comparable outcomes to busulfan-based conditioning in paediatric ALL, without the need for TDM to date. Future PK evaluation and modelling may optimise therapy and improve outcome. More recently, the addition of clofarabine to busulfan/fludarabine has shown encouraging results when compared to TBI-based regimens. The combination shows activity in ALL as well as AML and deserves further evaluation. Like busulfan, optimization of chemotherapy conditioning may be enhanced by understanding not just the PK of clofarabine, fludarabine, treosulfan and other agents, but also the pharmacodynamics and pharmacogenetics, ideally in the context of a single disease such as ALL.

Busulfan Pharmacokinetics in Adenosine Deaminase-Deficient Severe Combined Immunodeficiency Gene Therapy

The pharmacokinetics of low-dose busulfan (BU) were investigated as a nonmyeloablative conditioning regimen for autologous gene therapy (GT) in pediatric subjects with adenosine deaminase-deficient severe combined immunodeficiency disease (ADA SCID). In 3 successive clinical trials, which included either γ-retroviral (γ-RV) or lentiviral (LV) vectors, subjects were conditioned with BU using different dosing nomograms. The first cohort received BU doses based on body surface area (BSA), the second cohort received doses based on actual body weight (ABW), and in the third cohort, therapeutic drug monitoring (TDM) was used to target a specific area under the concentration-time curve (AUC). Neither BSA-based nor ABW-based dosing achieved a consistent cumulative BU AUC; in contrast, TDM-based dosing led to more consistent AUC. BU clearance increased as subject age increased from birth to 18 months. However, weight and age alone were insufficient to accurately predict the dose that would consistently achieve a target AUC. Furthermore, various clinical, laboratory, and genetic factors (eg, genotypes for glutathione-S-transferase isozymes known to participate in BU metabolism) were analyzed, but no single finding predicted subjects with rapid versus slow clearance. Analysis of BU AUC and the postengraftment vector copy number (VCN) in granulocytes, a surrogate marker of the level of engrafted gene-modified hematopoietic stem and progenitor cells (HSPCs), demonstrated gene marking at levels sufficient for therapeutic benefit in the subjects who had achieved the target BU AUC. Although many factors determine the ultimate engraftment following GT, this work demonstrates that the BU AUC correlated with the eventual level of engrafted gene-modified HSPCs within a vector group (γ-RV versus LV), with significantly higher levels of granulocyte VCN in the recipients of LV-modified grafts compared to recipients of γ-RV-transduced grafts. Taken together, these findings provide insight into low-dose BU pharmacokinetics in the unique setting of autologous GT for ADA SCID, and these dosing principles may be applied to future GT trials using low-dose BU to open the bone marrow niche.

Population pharmacokinetics of busulfan in Saudi pediatric patients undergoing hematopoietic stem cell transplantation

Background Busulfan is an antineoplastic drug that is used widely as part of a conditioning regimen in pediatric patients undergoing hematopoietic stem cell transplantation. It has a narrow therapeutic index and highly variable pharmacokinetics; therefore therapeutic drug monitoring is recommended to optimize busulfan dosing. Objective To study the population pharmacokinetics of busulfan in Saudi pediatric patients to optimize its dosing. Settings King Abdullah Specialist Children's Hospital in Riyadh, Saudi Arabia. Methods This pharmacokinetic observational study was conducted between January 2016 and December 2018. All pediatric patients receiving IV busulfan and undergoing routine therapeutic drug monitoring were included. Population pharmacokinetics modeling was conducted using Monolix2019R1. Pharmacokinetic data of busulfan in children. Results The study included 59 patients and 513 samples. The mean ± SD age was 6.10 ± 3.17 years, and the dose administered was 0.994 ± 0.15 mg/kg. The mean ± SD Cmax and area under the curve (AUC) were 900.60 ± 402.8 ng/mL and 1031.14 ± 300.75 µM min, respectively. Based on our simulations, the European Medicines Agency recommended dose were adequate for most patient's groups to achieve the conventional target of an AUC_0-tau of 900-1350 µM min. For patients in the lower weight group < 9 kg, higher doses were need at 1.2 mg/kg. With regards to the newly proposed target of AUC 78-101 mg h/mL, all of the doses we tested had low probability of achieving it. Conclusions Most of our patients had less than a proportional increase in busulfan concentration suggesting autoinduction. The high interindividual variability and autoinduction make dose adjustments challenging and AUC at steady state difficult to predict from the first dose. One approach to improve dose predictions is to use Bayesian dosing software. Based on our simulations, the European Medicines Agency recommended doses were adequate for most patient groups, except those in the lower (< 9 kg) and higher weight groups (> 34 kg).

Population pharmacokinetic analysis of intravenous busulfan: GSTA1 genotype is not a predictive factor of initial dose in Chinese adult patients undergoing hematopoietic stem cell transplantation

PURPOSE

This study aimed to develop a population pharmacokinetic (PPK) model to investigate the impact of GSTA1, GSTP1, and GSTM1 genotypes on busulfan pharmacokinetic (PK) variability in Chinese adult patients.

METHODS

Forty-three and 19 adult patients who underwent hematopoietic stem cell transplantation (HSCT) were enrolled for modeling group and validation group, respectively. All patients received twice-daily intravenous busulfan as part of conditioning regimen before HSCT. The PPK model was developed by nonlinear mixed-effect modeling. Covariates investigated were age, sex, actual body weight, body surface area, diagnoses, hepatic function markers, GST genotypes and conditioning regimen.

RESULTS

A total of 488 busulfan concentrations from 43 patients were obtained for the PPK model. The PK of intravenous busulfan was described by one-compartment model with first-order elimination with estimated clearance (CL) of 14.2 L/h and volume of distribution of 64.1 L. Inclusion of GSTA1 genotype as a covariate accounted for 1.1% of the inter-individual variability of busulfan CL (from 17.8% in the basic model to 16.7% in the final model). The accuracy and applicability of the final model were externally validated in the independent group. The difference of busulfan PK between Chinese patients and Caucasian patients existed because of the rarity of haplotype *B in Chinese population.

CONCLUSIONS

Although the GSTA1 genotype-based PPK model of intravenous busulfan was successfully developed and externally validated, the GSTA1 genotype was not considered to be clinically relevant to busulfan CL. We did not suggest the guidance of GSTA1 genotype on initial busulfan dose in Chinese adult patients.

Comparison of Different Conditioning Regimens in Allogeneic Hematopoietic Stem-Cell Transplantation Shows Superiority of Total Body Irradiation-Based Regimen for Younger Patients With Acute Leukemia: A Nationwide Study

BACKGROUND

The optimal the conditioning regimens for allogeneic hematopoietic stem-cell transplantation, especially for East Asian patients, remains unknown.

PATIENTS AND METHODS

We collected and analyzed clinical and survival data of 4255 patients from the Korean National Health Insurance Claims Database.

RESULTS

Between 1562 myeloablative conditioning and 2693 nonmyeloablative conditioning groups, the overall survival was not statistically different. However, in the myeloablative conditioning group, the overall survival of the total body irradiation-based regimen was better than that of chemotherapy-alone regimen (P = .005). In subgroup analysis, the superiority of the total body irradiation-based regimen was especially prominent in acute leukemia (P = .012 for acute myeloid leukemia; P = .005 for acute lymphoblastic leukemia) and for younger patients (< 50 years old vs. ≥ 50 years old, P = .015).

CONCLUSION

Total body irradiation combination might be the best conditioning regimen for young patients undergoing hematopoietic stem-cell transplantation for acute leukemias in Korea.

Personalized pharmacokinetic targeting with busulfan in allogeneic hematopoietic stem cell transplantation in infants with acute lymphoblastic leukemia

Individual busulfan (BU) dosing based on pharmacokinetic (PK) data is preferable for hematopoietic stem cell transplantation (HSCT) conditioning, but information on BU PK in infants is scarce. We report BU PK data on HSCT conditioning for infants with KMT2A-gene-rearrangement-positive acute lymphoblastic leukemia (MLL-r ALL). Infants showed wide variation in BU PK indices, such as clearance (CL) and volume of distribution (V_d) value, which are distributed more widely among those who received oral, rather than intravenous (IV), BU. Because the steady state concentration (C_ss) fluctuates readily in infants, dose re-adjustment based on PK at the initial administration was important even if the initial dose was determined by a PK test. HSCT can be performed safely within the C_ss range of 600-900 ng/mL per dose, although it was difficult to fit within the therapeutic index of BU. The clinical outcome of engraftment, graft-versus-host disease, adverse events, including sinusoidal obstruction syndrome, and survival did not correlate with the BU PK data, which paradoxically suggests that remaining within this C_ss range helped minimize transplant-related toxicities, while securing engraftment in infants with MLL-r ALL.

Effect of glutathione S-transferase genetic polymorphisms on busulfan pharmacokinetics and veno-occlusive disease in hematopoietic stem cell transplantation: A meta-analysis

This meta-analysis was conducted to derive an integrated conclusion about the influence of glutathione S-transferase (GST) genetic polymorphisms on busulfan pharmacokinetic (PK) parameters and veno-occlusive disease (VOD). Studies which analysed the effect of GST genetic polymorphisms on area under the curve (AUC), clearance (CL) or VOD were searched for and selected. A pooled analysis was conducted using Comprehensive Meta-Analysis programme. Nineteen studies were included in this meta-analysis. GSTA1*B and GSTM1 null genotypes significantly decreased CL_IV of busulfan (standardized difference in means (SDM) = -1.103; P = 0.019 and SDM = -0.418; P = 0.002, respectively). GSTA1*B significantly increased AUC_IV of busulfan (SDM = 0.832; P = 0.046), whereas GSTM1 did not (SDM = 0.155; P = 0.478). The PK parameters of oral busulfan did not differ according to GST genotype. GSTA1, GSTM1 and GSTP1 were not significantly associated with VOD occurrence. GSTA1 and GSTM1 genotypes affected CL_IV of busulfan, but only GSTA1 affected AUC_IV . There was no significant difference in the PK parameters of oral busulfan (CL_PO and AUC_PO ) and VOD when only GST genotypes were considered.

GSTA1 diplotypes affect busulfan clearance and toxicity in children undergoing allogeneic hematopoietic stem cell transplantation: a multicenter study

Busulfan (BU) dose adjustment following therapeutic drug monitoring contributes to better outcome of hematopoietic stem cell transplantation (HSCT). Further improvement could be achieved through genotype-guided BU dose adjustments. To investigate this aspect, polymorphism within glutathione S transferase genes were assessed. Particularly, promoter haplotypes of the glutathione S transferase A1 (GSTA1) were evaluated in vitro, with reporter gene assays and clinically, in a pediatric multi-center study (N =138) through association with BU pharmacokinetics (PK) and clinical outcomes. Promoter activity significantly differed between the GSTA1 haplotypes (p<0.001) supporting their importance in capturing PK variability. Four GSTA1 diplotype groups that significantly correlated with clearance (p=0.009) were distinguished. Diplotypes underlying fast and slow metabolizing capacity showed higher and lower BU clearance (ml/min/kg), respectively. GSTA1 diplotypes with slow metabolizing capacity were associated with higher incidence of sinusoidal obstruction syndrome, acute graft versus host disease and combined treatment-related toxicity (p<0.0005). Among other GST genes investigated, GSTP1 313GG correlated with acute graft versus host disease grade 1-4 (p=0.01) and GSTM1 non-null genotype was associated with hemorrhagic cystitis (p=0.003). This study further strengthens the hypothesis that GST diplotypes/genotypes could be incorporated into already existing population pharmacokinetic models for improving first BU dose prediction and HSCT outcomes. (N^o Clinicaltrials.gov identifier: NCT01257854. Registered 8 December 2010, retrospectively registered).

GSTA1 Genetic Variants and Conditioning Regimen: Missing Key Factors in Dosing Guidelines of Busulfan in Pediatric Hematopoietic Stem Cell Transplantation

Busulfan (Bu) is a key component of conditioning regimens used before hematopoietic stem cell transplantation (SCT) in children. Different predictive methods have been used to calculate the first dose of Bu. To evaluate the necessity of further improvements, we retrospectively analyzed the currently available weight- and age-based guidelines to calculate the first doses in 101 children who underwent allogenic SCT in CHU Sainte-Justine, Montreal, after an intravenous Bu-containing conditioning regimen according to genetic and clinical factors. The measured areas under the curve (AUCs) were within target (900 to 1500 µM/min) in 38.7% of patients after the administration of the first dose calculated based on age and weight, as locally recommended. GSTA1 diplotypes linked to poor Bu metabolism (G3) and fludarabine-containing regimens were the only factors associated with AUC within target (OR, 4.7 [95% CI, 1.1 to 19.8, P = .04]; and OR, 9.9 [95% CI, 1.6 to 61.7, P = .01], respectively). From the 11 methods selected for dose calculation, the percentage of AUCs within the target varied between 16% and 74%. In some models G3 was associated with AUCs within the therapeutic and the toxic range, whereas rapid metabolizers (G1) were correlated with subtherapeutic AUCs when different methods were used. These associations were confirmed by clearance-prediction analysis, in which GSTA1 diplotypes consistently influenced the prediction errors of the methods. These findings suggest that these factors should be considered in Bu dose prediction in addition to the anthropometric data from patients. Furthermore, our data indicated that GSTA1 diplotypes was a factor that should be included in future population pharmacokinetic models, including similar conditioning regiments, to improve the prediction of Bu exposure after its initial dose.

Pharmacometabolomics for predicting variable busulfan exposure in paediatric haematopoietic stem cell transplantation patients

Owing to its narrow therapeutic range and high pharmacokinetic variability, optimal dosing for busulfan is important to minimise overexposure-related systemic toxicity and underexposure-related graft failure. Using global metabolomics, we investigated biomarkers for predicting busulfan exposure. We analysed urine samples obtained before busulfan administration from 59 paediatric patients divided into 3 groups classified by area under the busulfan concentration-time curve (AUC), i.e., low-, medium-, and high-AUC groups. In the high-AUC group, deferoxamine metabolites were detected. Phenylacetylglutamine and two acylcarnitines were significantly lower in the high-AUC group than in the low-AUC group. Deferoxamine, an iron-chelating agent that lowers serum ferritin levels, was detected in the high-AUC group, indicating that those patients had high ferritin levels. Therefore, in a retrospective study of 130 paediatric patients, we confirmed our hypothesis that busulfan clearance (dose/AUC) and serum ferritin level has a negative correlation (r = −0.205, P = 0.019). Ferritin, acylcarnitine, and phenylacetylglutamine are associated with liver damage, including free radical formation, deregulation of hepatic mitochondrial β-oxidation, and hyperammonaemia. Our findings reveal potential biomarkers predictive of busulfan exposure and suggest that liver function may affect busulfan exposure.