Glutathione S-transferase homozygous deletions and relapse in childhood acute lymphoblastic leukemia: a novel study design in a large Italian AIEOP cohort

GSTM1 and GSTT1 double null genotypes determining cell fate and proliferation as potential risk factors of relapse in children with hematological malignancies after hematopoietic stem cell transplantation

PURPOSE

This study aimed to retrospectively evaluate the genetic association of null variants of glutathione S-transferases GSTM1 and GSTT1 with relapse incidence in children with hematological malignancies (HMs) undergoing busulfan (BU)- containing allogeneic hematopoietic stem cell transplantation (HSCT) and to assess the impact of these variants on BU-induced cytotoxicity on the immortalized lymphoblastoid cell lines (LCLs) and tumor THP1 GST gene-edited cell models.

METHODS

GSTM1- and GSTT1-null alleles were genotyped using germline DNA from whole blood prior to a conditioning BU-based regimen. Association of GSTM1- and GSTT1-null variants with relapse incidence was analyzed using multivariable competing risk analysis. BU-induced cell death studies were conducted in GSTs- null and non-null LCLs and CRISPR-Cas9 gene-edited THP1 leukemia cell lines.

RESULTS

Carrying GSTM1/GSTT1 double null genotype was found to be an independent risk factor for post-HSCT relapse in 86 children (adjusted HR: 6.52 [95% Cl, 2.76-15.42; p = 1.9 × 10^-5]). BU-induced cell death preferentially in THP1^{GSTM1(non-null)} and LCLs^{GSTM1(non-null)} as shown by decreased viability, increased necrosis and levels of the oxidized form of glutathione compared to null cells, while GSTT1 non-null cells showed increased baseline proliferation.

CONCLUSION

The clinical association suggests that GSTM1/GSTT1 double null genotype could serve as genetic stratification biomarker for the high risk of post-HSCT relapse. Functional studies have indicated that GSTM1 status modulates BU-induced cell death. On the other hand, GSTT1 is proposed to be involved in baseline cell proliferation.

Sampling strategies to evaluate the prognostic value of a new biomarker on a time-to-event end-point

BACKGROUND

The availability of large epidemiological or clinical data storing biological samples allow to study the prognostic value of novel biomarkers, but efficient designs are needed to select a subsample on which to measure them, for parsimony and economical reasons. Two-phase stratified sampling is a flexible approach to perform such sub-sampling, but literature on stratification variables to be used in the sampling and power evaluation is lacking especially for survival data.

METHODS

We compared the performance of different sampling designs to assess the prognostic value of a new biomarker on a time-to-event endpoint, applying a Cox model weighted by the inverse of the empirical inclusion probability.

RESULTS

Our simulation results suggest that case-control stratified (or post stratified) by a surrogate variable of the marker can yield higher performances than simple random, probability proportional to size, and case-control sampling. In the presence of high censoring rate, results showed an advantage of nested case-control and counter-matching designs in term of design effect, although the use of a fixed ratio between cases and controls might be disadvantageous. On real data on childhood acute lymphoblastic leukemia, we found that optimal sampling using pilot data is greatly efficient.

CONCLUSIONS

Our study suggests that, in our sample, case-control stratified by surrogate and nested case-control yield estimates and power comparable to estimates obtained in the full cohort while strongly decreasing the number of patients required. We recommend to plan the sample size and using sampling designs for exploration of novel biomarker in clinical cohort data.

A hPSC-based platform to discover gene-environment interactions that impact human β-cell and dopamine neuron survival

Common disorders, including diabetes and Parkinson’s disease, are caused by a combination of environmental factors and genetic susceptibility. However, defining the mechanisms underlying gene-environment interactions has been challenging due to the lack of a suitable experimental platform. Using pancreatic β-like cells derived from human pluripotent stem cells (hPSCs), we discovered that a commonly used pesticide, propargite, induces pancreatic β-cell death, a pathological hallmark of diabetes. Screening a panel of diverse hPSC-derived cell types we extended this observation to a similar susceptibility in midbrain dopamine neurons, a cell type affected in Parkinson’s disease. We assessed gene-environment interactions using isogenic hPSC lines for genetic variants associated with diabetes and Parkinson’s disease. We found GSTT1^−/− pancreatic β-like cells and dopamine neurons were both hypersensitive to propargite-induced cell death. Our study identifies an environmental chemical that contributes to human β-cell and dopamine neuron loss and validates a novel hPSC-based platform for determining gene-environment interactions.

Diseases such as diabetes and Parkinson's manifest based on interactions between genes and environment. Here, the authors find among a panel of cell types that propargite, a common pesticide, induces pancreatic β-cell and dopamine neuron death and that loss of the gene GSTT1 confers hypersensitivity.

Pharmacogenomic markers of glucocorticoid response in the initial phase of remission induction therapy in childhood acute lymphoblastic leukemia

Background

Response to glucocorticoid (GC) monotherapy in the initial phase of remission induction treatment in childhood acute lymphoblastic leukemia (ALL) represents important biomarker of prognosis and outcome. We aimed to study variants in several pharmacogenes (NR3C1, GSTs and ABCB1) that could contribute to improvement of GC response through personalization of GC therapy.

Methods

Retrospective study enrolling 122 ALL patients was carried out to analyze variants of NR3C1 (rs33389, rs33388 and rs6198), GSTT1 (null genotype), GSTM1 (null genotype), GSTP1 (rs1695 and rs1138272) and ABCB1 (rs1128503, rs2032582 and rs1045642) genes using PCR-based methodology. The marker of GC response was blast count per microliter of peripheral blood on treatment day 8. We carried out analysis in which cut-off value for GC response was 1000 (according to Berlin-Frankfurt-Munster [BFM] protocol), as well as 100 or 0 blasts per microliter.

Results

Carriers of rare NR3C1 rs6198 GG genotype were more likely to have blast count over 1000, than the non-carriers (p = 0.030). NR3C1 CAA (rs33389-rs33388-rs6198) haplotype was associated with blast number below 1000 (p = 0.030). GSTP1 GC haplotype carriers were more likely to have blast number below 1000 (p = 0.036), below 100 (p = 0.028) and to be blast negative (p = 0.054), while GSTP1 GT haplotype and rs1138272 T allele carriers were more likely to be blasts positive (p = 0.034 and p = 0.024, respectively). ABCB1 CGT (rs1128503-rs2032582-rs1045642) haplotype carriers were more likely to be blast positive (p = 0.018).

Conclusions

Our results have shown that NR3C1 rs6198 variant and GSTP1 rs1695-rs1138272 haplotype are the most promising pharmacogenomic markers of GC response in ALL patients.

Association Between Combined Presence of Hepatitis C Virus and Polymorphisms in Different Genes With Toxicities of Methotrexate and 6-Mercaptopurine in Children With Acute Lymphoblastic Leukemia

BACKGROUND

The aim of the present study is to determine the correlation of hepatitis C virus (HCV) infection and polymorphisms in different genes with toxicity of either methotrexate (MTX) or 6-mercaptopurine (6-MP) administered to children with acute lymphoblastic leukemia (ALL).

PROCEDURE

One hundred children with low-risk ALL, who were treated according to the St. Jude Total therapy XV, were recruited. The recruited children were receiving MTX and 6-MP during maintenance phase. Patients were excluded from the study if they had other types of leukemia. Genotyping analyses for the thiopurine methyltransferase (TPMT), methylenetetrahydrofolate reductase (MTHFR), and glutathione S-transferase (GST) genes were performed using a combination of polymerase chain reaction (PCR) and PCR-RFLP (where RFLP is restriction fragment length polymorphism) protocols. Relevant clinical data on adverse drug reactions were collected objectively (blinded to genotypes) from the patient medical records.

RESULTS

There was a significant correlation between the combined presence of HCV and TPMT*3B G460A gene polymorphisms and grades 2-4 hepatotoxicity as aspartate aminotransferase (AST) elevation (P < 0.04). The same observation was seen when comparing either the presence of HCV alone or the presence of the gene polymorphism alone. A significant association between the combined presence of HCV and MTHFR C677T polymorphism and grades 2-4 hepatotoxicity as alanine aminotransferase (ALT), AST, and alkaline phosphatase (ALP) elevation was observed (P values <0.001, 0.02, and 0.001, respectively). The presence of HCV infection had a significant negative effect on hepatic transaminases.

CONCLUSIONS

The present data support a role for combining analysis of genetic variation in drug-metabolizing enzymes and the presence of HCV in the assessment of specific drugs toxicities in multiagent chemotherapeutic treatment regimens.

Crude incidence in two-phase designs in the presence of competing risks

Background

In many studies, some information might not be available for the whole cohort, some covariates, or even the outcome, might be ascertained in selected subsamples. These studies are part of a broad category termed two-phase studies. Common examples include the nested case-control and the case-cohort designs. For two-phase studies, appropriate weighted survival estimates have been derived; however, no estimator of cumulative incidence accounting for competing events has been proposed. This is relevant in the presence of multiple types of events, where estimation of event type specific quantities are needed for evaluating outcome.

Methods

We develop a non parametric estimator of the cumulative incidence function of events accounting for possible competing events. It handles a general sampling design by weights derived from the sampling probabilities. The variance is derived from the influence function of the subdistribution hazard.

Results

The proposed method shows good performance in simulations. It is applied to estimate the crude incidence of relapse in childhood acute lymphoblastic leukemia in groups defined by a genotype not available for everyone in a cohort of nearly 2000 patients, where death due to toxicity acted as a competing event. In a second example the aim was to estimate engagement in care of a cohort of HIV patients in resource limited setting, where for some patients the outcome itself was missing due to lost to follow-up. A sampling based approach was used to identify outcome in a subsample of lost patients and to obtain a valid estimate of connection to care.

Conclusions

A valid estimator for cumulative incidence of events accounting for competing risks under a general sampling design from an infinite target population is derived.

Pharmacogenetics and induction/consolidation therapy toxicities in acute lymphoblastic leukemia patients treated with AIEOP-BFM ALL 2000 protocol

Drug-related toxicities represent an important clinical concern in chemotherapy, genetic variants could help tailoring treatment to patient. A pharmacogenetic multicentric study was performed on 508 pediatric acute lymphoblastic leukemia patients treated with AIEOP-BFM 2000 protocol: 28 variants were genotyped by VeraCode and Taqman technologies, deletions of GST-M1 and GST-T1 by multiplex PCR. Toxicities were derived from a central database: 251 patients (49.4%) experienced at least one gastrointestinal (GI) or hepatic (HEP) or neurological (NEU) grade III/IV episode during the remission induction phase: GI occurred in 63 patients (12.4%); HEP in 204 (40.2%) and NEU in 44 (8.7%). Logistic regression model adjusted for sex, risk and treatment phase revealed that ITPA rs1127354 homozygous mutated patients showed an increased risk of severe GI and NEU. ABCC1 rs246240 and ADORA2A rs2236624 homozygous mutated genotypes were associated to NEU and HEP, respectively. These three variants could be putative predictive markers for chemotherapy-related toxicities in AIEOP-BFM protocols.

Role of Pharmacogenetics in Hematopoietic Stem Cell Transplantation Outcome in Children

Hematopoietic stem cell transplantation (HSCT) is an established therapeutic procedure for several congenital and acquired disorders, both malignant and nonmalignant. Despite the great improvements in HSCT clinical practices over the last few decades, complications, such as graft vs. host disease (GVHD) and sinusoidal obstructive syndrome (SOS), are still largely unpredictable and remain the major causes of morbidity and mortality. Both donor and patient genetic background might influence the success of bone marrow transplantation and could at least partially explain the inter-individual variability in HSCT outcome. This review summarizes some of the recent studies on candidate gene polymorphisms in HSCT, with particular reference to pediatric cohorts. The interest is especially focused on pharmacogenetic variants affecting myeloablative and immunosuppressive drugs, although genetic traits involved in SOS susceptibility and transplant-related mortality are also reviewed.

Non-monotonic relationships between arsenic and selenium excretion and its implication on arsenic methylation pattern in a Bangladeshi population

The toxicity of arsenic differs markedly between individuals and populations, which might be related to the metabolism (methylation) of inorganic arsenic (As), as well as the selenium (Se) nutritional status. Urinary excretion of As (u-As) and Se (u-Se) was examined in an adult population (n=128) living in an As-contaminated area in Bangladesh. Although there was a significant negative correlation between u-Se and u-As (median 137; range 49-927 μg/g creatinine), closer examination revealed a non-monotonous relationship between them. A quadratic curve with an axis of As at 155 μg/g Cre gave a better fit, and u-As and u-Se were positively or negatively correlated depending on whether the As concentration was lower or higher than 155 μg As/g Cre, respectively. Likewise, the relationships between the As methylation pattern and glutathione-S-transferase (GST) polymorphism, body mass index (BMI), and u-Se differed depending on the u-As range; i.e., higher or lower than 155 μg/g Cre. Although we did not determine the causal mechanism for these observations, the non-monotonic relationship between As exposure and the variables examined suggested the existence of a threshold at which the handling of As by human body is qualitatively changed. The possible importance of Se nutrition for As toxicity is also discussed.

TNF-α SNP rs1800629 and risk of relapse in childhood acute lymphoblastic leukemia: relation to immunophenotype

AIM

In the AIEOP-BFM ALL (Associazione Italiana Ematologia Oncologia Pediatrica-Berlin Frankfurt Münster acute lymphoblastic leukemia) 2000 protocol, 70% of relapsed patients had favorable prognostic features and fell within less intensive polychemotherapeutic regimens, suggesting the need for better assessing lower risk stratification.

MATERIALS & METHODS

A novel two-phase study design selected 614 children to be genotyped for TNF-α SNP rs1800629 (-308G>A). A weighted Cox model was applied to evaluate the SNP effect on hazard of relapse, adjusting for immunophenotype, risk group, age and gender and including interaction terms.

RESULTS

Significant interaction was found with immunophenotypes (p = 0.0007, with minor allele genotypes being adverse genetic markers in B-cell acute lymphoblastic leukemia and protective ones in T-cell acute lymphoblastic leukemia), and also with risk protocols (p = 0.0041, with minor allele genotypes as prognostic factor of relapse for standard risk patients [only one T-cell acute lymphoblastic leukemia in the subgroup analyzed]).

CONCLUSION

The presence of at least one A allele in TNF-α SNP rs1800629 should suggest a closer monitoring in B-cell acute lymphoblastic leukemia standard risk patients.

[Structure and functions of glutathione transferases]

Data about classification, nomenclature, structure, substrate specificity and role of many glutathione transferase's isoenzymes in cell functions have been summarised. The enzyme has been discovered more than 50 years ago. This family of proteins is updated continuously. It has very different composition and will have demand for system analysis for many years.

Pharmacogenetics of childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer in developed countries. Although treatment outcome has improved owing to advances in chemotherapy, there is still a group of patients for which therapy fails while some patients experience severe toxicity. In the last few years, several pharmacogenetic studies have been performed to search for markers of outcome and toxicity in pediatric ALL. However, to date, TPMT is the only pharmacogenetic marker in ALL with clinical guidelines for drug dosing. In this article, we will provide an overview of the most important findings carried out in pharmacogenetics for pediatric ALL, such as the interest drawn by methotrexate transporters in the context of methotrexate treatment. Even if most of the studies are centered on coding genes, we will also point to new approaches focusing on noncoding regions and epigenetic variation that could be interesting for consideration in the near future.

Survival estimation in two-phase cohort studies with application to biomarkers evaluation

Two-phase studies are attractive for their economy and efficiency in research settings where large cohorts are available for investigating the prognostic and predictive role of novel genetic and biological factors. In this type of study, information on novel factors is collected only in a convenient subcohort (phase II) drawn from the cohort (phase I) according to a given (optimal) sampling strategy. Estimation of survival in the subcohort needs to account for the design. The Kaplan-Meier method, based on counts of events and of subjects at risk in time, must be applied accounting, with suitable weights, for the sampling probabilities of the subjects in phase II, in order to recover the representativeness of the subcohort for the entire cohort. The authors derived a proper variance estimator of survival by linearization. The proposed method is applied in the context of a two-phase study on childhood acute lymphoblastic leukemia, which was planned in order to evaluate the role of genetic polymorphisms on treatment failure due to relapse. The method has shown satisfactory performance through simulations under different scenarios, including the case-control setting, and proved to be useful for describing results in the clinical example.

Predicting relapse risk in childhood acute lymphoblastic leukaemia

Intensive multi-agent chemotherapy regimens and the introduction of risk-stratified therapy have substantially improved cure rates for children with acute lymphoblastic leukaemia (ALL). Current risk allocation schemas are imperfect, as some children are classified as lower-risk and treated with less intensive therapy relapse, while others deemed higher-risk are probably over-treated. Most cooperative groups previously used morphological clearance of blasts in blood and marrow during the initial phases of chemotherapy as a primary factor for risk group allocation; however, this has largely been replaced by the detection of minimal residual disease (MRD). Other than age and white blood cell count (WBC) at presentation, many clinical variables previously used for risk group allocation are no longer prognostic, as MRD and the presence of sentinel genetic lesions are more reliable at predicting outcome. Currently, a number of sentinel genetic lesions are used by most cooperative groups for risk stratification; however, in the near future patients will probably be risk-stratified using genomic signatures and clustering algorithms, rather than individual genetic alterations. This review will describe the clinical, biological, and response-based features known to predict relapse risk in childhood ALL, including those currently used and those likely to be used in the near future to risk-stratify therapy.