Strengths and weaknesses of gene association studies in childhood acute lymphoblastic leukemia

TS Gene Polymorphisms Correlate with Susceptibility to Acute Lymphocytic Leukemia in Children

Background

Acute lymphocytic leukemia (ALL) in children is a clonal disease of bone marrow hematopoietic stem cells. This study aimed to explore the associations between MTHFR or TS genetic polymorphisms and susceptibility to acute lymphocytic leukemia (ALL) in children.

Material/Methods

This case-control study included 79 ALL patients (case group) and 102 non-ALL patients (control group). Post-PCR genomic DNA sequencing revealed MTHFR C677T and MTHFR A1298C genotypes and TS polymorphisms. The χ² test was used to compare differences in MTHFR and TS polymorphisms (including genotypic and allelic distributions) between groups. Logistic regression analysis was used to determine genetic polymorphisms and ALL risk associations.

Results

The results indicated that TS 3R allele frequency was significantly higher in the case group than in the control group (χ²=7.45, P<0.05). The MTHFR C677T and MTHFR A1298C polymorphisms were not associated with ALL risk. Compared to the TS 2R/2R genotype, subjects carrying TS 2R/3R were twice as likely to develop ALL, and the TS 3R/3R+3R/4R genotype carried a 4-fold higher risk of developing ALL (OR=1.96, CI: 1.14–3.36).

Conclusions

The TS genetic polymorphisms increase the ALL risk. The TS 3R allele was a risk factor for ALL. There were no associations between MTHFR C677T or MTHFR A1298C polymorphisms and ALL susceptibility.

Roles of genetic polymorphisms in the folate pathway in childhood acute lymphoblastic leukemia evaluated by Bayesian relevance and effect size analysis

In this study we investigated whether polymorphisms in the folate pathway influenced the risk of childhood acute lymphoblastic leukemia (ALL) or the survival rate of the patients. For this we selected and genotyped 67 SNPs in 15 genes in the folate pathway in 543 children with ALL and 529 controls. The results were evaluated by gender adjusted logistic regression and by the Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) methods. Bayesian structure based odds ratios for the relevant variables and interactions were also calculated. Altogether 9 SNPs in 8 genes were associated with altered susceptibility to ALL. After correction for multiple testing, two associations remained significant. The genotype distribution of the MTHFD1 rs1076991 differed significantly between the ALL and control population. Analyzing the subtypes of the disease the GG genotype increased only the risk of B-cell ALL (p = 3.52×10⁻⁴; OR = 2.00). The GG genotype of the rs3776455 SNP in the MTRR gene was associated with a significantly reduced risk to ALL (p = 1.21×10⁻³; OR = 0.55), which resulted mainly from the reduced risk to B-cell and hyperdiploid-ALL. The TC genotype of the rs9909104 SNP in the SHMT1 gene was associated with a lower survival rate comparing it to the TT genotype (80.2% vs. 88.8%; p = 0.01). The BN-BMLA confirmed the main findings of the frequentist-based analysis and showed structural interactional maps and the probabilities of the different structural association types of the relevant SNPs especially in the hyperdiploid-ALL, involving additional SNPs in genes like TYMS, DHFR and GGH. We also investigated the statistical interactions and redundancies using structural model properties. These results gave further evidence that polymorphisms in the folate pathway could influence the ALL risk and the effectiveness of the therapy. It was also shown that in gene association studies the BN-BMLA could be a useful supplementary to the traditional frequentist-based statistical method.

Influence of Methylenetetrahydrofolate Reductase C677T, A1298C, and G80A Polymorphisms on the Survival of Pediatric Patients with Acute Lymphoblastic Leukemia

The influence of genic polymorphisms involved in metabolism of chemotherapeutic agents as the methotrexate (MTX) has been studied mainly in acute lymphoblastic leukemia (ALL) of childhood. Advances in treatment may be attributed to identification of prognostic factors added to chemotherapy protocol. The aim of this study was to analyze the association of the C677T, A1298C, and G80A polymorphisms on MTHFR gene and on the overall survival of pediatric patients (n = 126) with lymphoblastic leukemia treated with MTX according to the Brazilian protocol in 187 months. The C677T and G80A polymorphisms were genotyped by PCR-RFLP and A1298C polymorphism by allele-specific PCR. We observed that ALL patients presented rate (dead/alive) of 0.36 for the 677CC genotype, corresponding also to lower overall survival (P = 0.0013); on the other hand, the 677TT genotype showed a better survival (98%). Thus, we believe that patients with 80AA genotype presented a small reduction in MTX plasma level, suggesting that ALL children, carrying the 80AA genotype, showed a high toxicity to MTX (P < 0.0001).

Candidate gene association study in pediatric acute lymphoblastic leukemia evaluated by Bayesian network based Bayesian multilevel analysis of relevance

Background

We carried out a candidate gene association study in pediatric acute lymphoblastic leukemia (ALL) to identify possible genetic risk factors in a Hungarian population.

Methods

The results were evaluated with traditional statistical methods and with our newly developed Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) method. We collected genomic DNA and clinical data from 543 children, who underwent chemotherapy due to ALL, and 529 healthy controls. Altogether 66 single nucleotide polymorphisms (SNPs) in 19 candidate genes were genotyped.

Results

With logistic regression, we identified 6 SNPs in the ARID5B and IKZF1 genes associated with increased risk to B-cell ALL, and two SNPs in the STAT3 gene, which decreased the risk to hyperdiploid ALL. Because the associated SNPs were in linkage in each gene, these associations corresponded to one signal per gene. The odds ratio (OR) associated with the tag SNPs were: OR = 1.69, P = 2.22x10^-7 for rs4132601 (IKZF1), OR = 1.53, P = 1.95x10^-5 for rs10821936 (ARID5B) and OR = 0.64, P = 2.32x10^-4 for rs12949918 (STAT3). With the BN-BMLA we confirmed the findings of the frequentist-based method and received additional information about the nature of the relations between the SNPs and the disease. E.g. the rs10821936 in ARID5B and rs17405722 in STAT3 showed a weak interaction, and in case of T-cell lineage sample group, the gender showed a weak interaction with three SNPs in three genes. In the hyperdiploid patient group the BN-BMLA detected a strong interaction among SNPs in the NOTCH1, STAT1, STAT3 and BCL2 genes. Evaluating the survival rate of the patients with ALL, the BN-BMLA showed that besides risk groups and subtypes, genetic variations in the BAX and CEBPA genes might also influence the probability of survival of the patients.

Conclusions

In the present study we confirmed the roles of genetic variations in ARID5B and IKZF1 in the susceptibility to B-cell ALL. With the newly developed BN-BMLA method several gene-gene, gene-phenotype and phenotype-phenotype connections were revealed. We showed several advantageous features of the new method, and suggested that in gene association studies the BN-BMLA might be a useful supplementary to the traditional frequentist-based statistical method.