MiR-128-3p as blood based liquid biopsy biomarker in childhood acute lymphoblastic leukemia.

BACKGROUND: Minimal residual disease (MRD) is one of the most valuable independent prognostic factors in acute lymphoblastic leukemia (ALL). Bone marrow (BM) aspiration, however, is an invasive process. Previous studies have shown that microRNAs (miR) and extracellular vesicle (EV)-related miRs show different expression profiles at the presence of malignant cells compared to healthy controls. In our previous project, we have reported that two miRs previously described to be overexpressed in blasts were significantly decreased over the first week of the therapy of patients with ALL in the platelet free plasma fraction (PFP) of peripheral blood samples (PB). The aim of the current study was to assess the relation between day 15 flow cytometry (FC) MRD and expression of miR-128-3p and miR-222-3p miRs in exosome-enriched fraction (EEF) of PFP to evaluate whether their expression in EEF correlates with day 15 FC MRD more precisely. METHODS: PB was collected from 13 patients diagnosed with pediatric pre-B ALL at 4 time points. Expression of miR-128-3p and miR-222-3p was measured by qPCR in PFP and EEF. RESULTS: Positive correlation was found between changes of miR-128-3p expression in EEF or PFP by day 8 of chemotherapy and day 15 FC MRD (rEEF = 0.99, pEEF = 1.13E-9 and rPFP = 0.99, pPFP = 4.75E-9, respectively). Furthermore, the decrease of miR-128-3p in EEF by day 15 of treatment also showed a positive correlation with day 15 FC MRD (rEEF = 0.96; pEEF = 4.89E-5). CONCLUSION: Our results show that circulating miRs are potential biomarkers of ALL MRD, asmiR-128-3p level both in PFP and EEF predicts day 15 FC MRD. In addition, the assessment of the EEF gave a more promising result.

Co-Detection of VEGF-A and Its Regulator, microRNA-181a, May Indicate Central Nervous System Involvement in Pediatric Leukemia.

Central nervous system (CNS) involvement is a leading cause of therapy-refractory pediatric acute lymphoblastic leukemia (pALL), which is aggravated by underdiagnosing CNS disease with the currently used cell-based approach of cerebrospinal fluid (CSF) diagnostics. Our study focused on developing novel subcellular CNS leukemia indicators in the CSF and the bone marrow (BM) of patients with pALL. Serial liquid biopsy samples (n = 65) were analyzed by Elisas to measure the level of essential proteins associated with blast cell CNS trafficking, vascular endothelial growth factor A (VEGF-A) and integrin alpha 6 (ITGA6). In CSF samples from early induction chemotherapy, VEGF-A concentration were uniformly elevated in the CNS-positive group compared to those patients without unambiguous meningeal infiltration (9 vs Nine patients, Δc = 17.2 pg/ml, p = 0.016). Expression of miR-181a, a VEGFA-regulating microRNA which showed increased level in CNS leukemia in our previous experiments, was then paralleled with VEGF-A concentration. A slight correlation between the levels of miR-181a and VEGF-A indicators in CSF and BM samples was revealed (n = 46, Pearson's r = 0.36, p = 0.015). After validating in international cohorts, the joint quantification of miR-181a and VEGF-A might provide a novel tool to precisely diagnose CNS involvement and adjust CNS-directed therapy in pALL.

Melanoma susceptibility: an update on genetic and epigenetic findings.

Malignant melanoma is one of the most highly ranked cancers in terms of years of life lost. Hereditary melanoma with its increased familial susceptibility is thought to affect up to 12% of all melanoma patients. In the past, only a few high-penetrance genes associated with familial melanoma, such as CDKN2A and CDK4, have been clinically tested. However, findings now indicate that melanoma is a cancer most likely to develop not only due to high-penetrance variants but also due to polygenic inheritance patterns, leaving no clear division between the hereditary and sporadic development of malignant melanoma. Various pathogenic low-penetrance variants were recently discovered through genome-wide association studies, and are now translated into polygenic risk scores. These can show superior sensitivity rates for the prediction of melanoma susceptibility and related mixed cancer syndromes than risk scores based on phenotypic traits of the patients, with odds ratios of up to 5.7 for patients in risk groups. In addition to describing genetic findings, we also review the first results of epigenetic research showing constitutional methylation changes that alter the susceptibility to cutaneous melanoma and its risk factors.

Pharmacogenetics of the Central Nervous System-Toxicity and Relapse Affecting the CNS in Pediatric Acute Lymphoblastic Leukemia.

Despite improving cure rates in childhood acute lymphoblastic leukemia (ALL), therapeutic side effects and relapse are ongoing challenges. These can also affect the central nervous system (CNS). Our aim was to identify germline gene polymorphisms that influence the risk of CNS events. Sixty single nucleotide polymorphisms (SNPs) in 20 genes were genotyped in a Hungarian non-matched ALL cohort of 36 cases with chemotherapy related acute toxic encephalopathy (ATE) and 544 controls. Five significant SNPs were further analyzed in an extended Austrian-Czech-NOPHO cohort (n = 107 cases, n = 211 controls) but none of the associations could be validated. Overall populations including all nations' matched cohorts for ATE (n = 426) with seizure subgroup (n = 133) and posterior reversible encephalopathy syndrome (PRES, n = 251) were analyzed, as well. We found that patients with ABCB1 rs1045642, rs1128503 or rs2032582 TT genotypes were more prone to have seizures but those with rs1045642 TT developed PRES less frequently. The same SNPs were also examined in relation to ALL relapse on a case-control matched cohort of 320 patients from all groups. Those with rs1128503 CC or rs2032582 GG genotypes showed higher incidence of CNS relapse. Our results suggest that blood-brain-barrier drug transporter gene-polymorphisms might have an inverse association with seizures and CNS relapse.

Investigation of circulating lncRNAs as potential biomarkers in chronic respiratory diseases.

BACKGROUND: In the present study the blood expression level of inflammatory response and autoimmunity associated long non-coding RNAs (lncRNAs) were compared in patients with different chronic respiratory diseases and investigated whether they could be used as biomarkers in these diseases. METHODS: In the discovery cohort, the gene expression level of 84 lncRNAs were measured in the blood of 24 adult patients including healthy controls and patients with asthma and COPD. In the replication cohort the expression of 6 selected lncRNAs were measured in 163 subjects including healthy controls and adults with allergic rhinitis, asthma, COPD and children with asthma. It was evaluated whether these lncRNAs can be used as diagnostic biomarkers for any studied disease. With systems biology analysis the biological functions of the selected lncRNAs were predicted. RESULTS: In the discovery cohort, the mean expression of 27 lncRNAs showed nominally significant differences in at least one comparison. OIP5-AS1, HNRNPU, RP11-325K4.3, JPX, RP11-282O18.3, MZF1-AS1 were selected for measurement in the replication cohort. Three lncRNAs (HNRNPU, RP11-325K4.3, JPX) expressed significantly higher in healthy children than in adult controls. All the mean expression level of the 6 lncRNAs differed significantly between adult allergic rhinitis patients and controls. RP11-325K4.3, HNRNPU and OIP5-AS1 expressed higher in allergic asthma than in non-allergic asthma. COPD and asthma differed in the expression of RP11-325K4.3 from each other. In examining of the lncRNAs as biomarkers the weighted accuracy (WA) values were especially high in the comparison of healthy controls and patients with allergic rhinitis. OIP5-AS1 and JPX achieved 0.98 and 0.9 WA values, respectively, and the combination of the selected lncRNAs also resulted in a high performance (WA = 0.98). Altogether, OIP5-AS1 had the highest discriminative power in case of three out of six comparisons. CONCLUSION: Differences were detected in the expression of circulating lncRNAs in chronic respiratory diseases. Some of these differences might be utilized as biomarkers and also suggest a possible role of these lncRNAs in the pathomechanism of these diseases. The lncRNAs and the associated pathways are potential therapeutic targets in these diseases, but naturally additional studies are needed for the confirmation of these results.

MicroRNA-181a as novel liquid biopsy marker of central nervous system involvement in pediatric acute lymphoblastic leukemia.

BACKGROUND: Refractory central nervous system (CNS) involvement is among the major causes of therapy failure in childhood acute leukemia. Applying contemporary diagnostic methods, CNS disease is often underdiagnosed. To explore more sensitive and less invasive CNS status indicators, we examined microRNA (miR) expressions and extracellular vesicle (EV) characteristics. METHODS: In an acute lymphoblastic leukemia (ALL) discovery cohort, 47 miRs were screened using Custom TaqMan Advanced Low-Density Array gene expression cards. As a validation step, a candidate miR family was further scrutinized with TaqMan Advanced miRNA Assays on serial cerebrospinal fluid (CSF), bone marrow (BM) and peripheral blood samples with different acute leukemia subtypes. Furthermore, small EV-rich fractions were isolated from CSF and the samples were processed for immunoelectron microscopy with anti-CD63 and anti-CD81 antibodies, simultaneously. RESULTS: Regarding the discovery study, principal component analysis identified the role of miR-181-family (miR-181a-5p, miR-181b-5p, miR-181c-5p) in clustering CNS-positive (CNS+) and CNS-negative (CNS‒) CSF samples. We were able to validate miR-181a expression differences: it was about 52 times higher in CSF samples of CNS+ ALL patients compared to CNS‒ cases (n = 8 vs. n = 10, ΔFC = 52.30, p = 1.5E-4), and CNS+ precursor B cell subgroup also had ninefold higher miR-181a levels in their BM (p = 0.04). The sensitivity of CSF miR-181a measurement in ALL highly exceeded those of conventional cytospin in the initial diagnosis of CNS leukemia (90% vs. 54.5%). Pellet resulting from ultracentrifugation of CNS+ CSF samples of ALL patients showed atypical CD63-/CD81- small EVs in high density by immunoelectron microscopy. CONCLUSIONS: After validating in extensive cohorts, quantification of miR-181a or a specific EV subtype might provide novel tools to monitor CNS disease course and further adjust CNS-directed therapy in pediatric ALL.

Circulating microRNAs as minimal residual disease biomarkers in childhood acute lymphoblastic leukemia.

BACKGROUND: Treatment stratification based on bone marrow minimal residual disease (MRD) at set time points has resulted in considerably improved survival in pediatric acute lymphoblastic leukemia (ALL). Treatment response is assessed using bone marrow samples. MicroRNAs (miRs) easily traffic among fluid spaces and are more stable than most other RNA classes. We examined the role of circulating miRs as putative less invasive MRD biomarkers. METHODS: In an exploratory experiment, expression of 46 preselected miRs was studied in platelet-free blood plasma samples of 15 de novo, 5 relapsed ALL patients and 10 controls by Custom TaqMan Array Advanced MicroRNA Card. Based on their high expression in ALL compared to controls, and on the reduction observed along the induction therapy, four miRs were selected for further analyses: miR-128-3p, -181a-5p, -181b-5p and 222-3p. Their expression was measured by qPCR at 4 time points in 27 de novo ALL patients treated in the ALL IC-BFM 2009 study. RESULTS: The expression of all 4 miRs significantly decreased over the first week of therapy (miR-128-3p: log2 fold change - 2.86; adjusted p 3.6 × 10-7; miR-181b-5p: log2 fold change - 1.75; adjusted p 1.48 × 10-2; miR-181a-5p: log2 fold change -1.33; adjusted p 3.12 × 10-2; miR-222-3p: log2 fold change - 1.25; adjusted p 1.66 × 10-2). However, no significant further reduction in miR expression was found after the 8th day of therapy. Measured drop in expression of 2 miRs at day 8 strongly correlated with day 15 bone marrow flow cytometry MRD results (miR-128-3p: Pearson's r = 0.88, adjusted p = 2.71 × 10-4; miR-222-3p: r = 0.81, adjusted p = 2.99 × 10-3). CONCLUSION: In conclusion, these circulating miRs might act as biomarkers of residual leukemia. MiR-128-3p and miR-222-3p in blood predict day 15 flow cytometry MRD results 7 days earlier. Although, their sensitivity falls behind that of bone marrow flow cytometry MRD at day 15.

Possible roles of genetic variations in chemotherapy related cardiotoxicity in pediatric acute lymphoblastic leukemia and osteosarcoma.

BACKGROUND: The treatment of acute lymphoblastic leukemia (ALL) and osteosarcoma (OSC) is very effective: the vast majority of patients recover and survive for decades. However, they still need to face serious adverse effects of chemotherapy. One of these is cardiotoxicity which may lead to progressive heart failure in the long term. Cardiotoxicity is contributed mainly to the use of anthracyclines and might have genetic risk factors. Our goal was to test the association between left ventricular function and genetic variations of candidate genes. METHODS: Echocardiography data from medical records of 622 pediatric ALL and 39 OSC patients were collected from the period 1989-2015. Fractional shortening (FS) and ejection fraction (EF) were determined, 70 single nucleotide polymorphisms (SNPs) in 26 genes were genotyped. Multivariate logistic regression and multi-adjusted general linear model were performed to investigate the influence of genetic polymorphisms on the left ventricular parameters. Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) method was applied to test for the potential interaction of the studied cofactors and SNPs. RESULTS: Our results indicate that variations in ABCC2, CYP3A5, NQO1, SLC22A6 and SLC28A3 genes might influence the left ventricular parameters. CYP3A5 rs4646450 TT was 17% among ALL cases with FS lower than 28, and 3% in ALL patients without pathological FS (p = 5.60E-03; OR = 6.94 (1.76-27.39)). SLC28A3 rs7853758 AA was 12% in ALL cases population, while only 1% among controls (p = 6.50E-03; OR = 11.56 (1.98-67.45)). Patients with ABCC2 rs3740066 GG genotype had lower FS during the acute phase of therapy and 5-10 years after treatment (p = 7.38E-03, p = 7.11E-04, respectively). NQO1 rs1043470 rare T allele was associated with lower left ventricular function in the acute phase and 5-10 years after the diagnosis (p = 4.28E-03 and 5.82E-03, respectively), and SLC22A6 gene rs6591722 AA genotype was associated with lower mean FS (p = 1.71E-03), 5-10 years after the diagnosis. CONCLUSIONS: Genetic variants in transporters and metabolic enzymes might modulate the individual risk to cardiac toxicity after chemotherapy.

HLA-DRB1*07:01-HLA-DQA1*02:01-HLA-DQB1*02:02 haplotype is associated with a high risk of asparaginase hypersensitivity in acute lymphoblastic leukemia.

Hypersensitivity reactions are the most frequent dose-limiting adverse reactions to Escherichia coli-derived asparaginase in pediatric acute lymphoblastic leukemia (ALL) patients. The aim of the present study was to identify associations between sequence-based Human Leukocyte Antigen Class II region alleles and asparaginase hypersensitivity in a Hungarian ALL population. Four-digit typing of HLA-DRB1 and HLA-DQB1 loci was performed in 359 pediatric ALL patients by using next-generation sequencing method. Based on genotypic data of the two loci, haplotype reconstruction was carried out. In order to investigate the possible role of the HLA-DQ complex, the HLA-DQA1 alleles were also inferred. Multivariate logistic regression analysis and a Bayesian network-based approach were applied to identify relevant genetic risk factors of asparaginase hypersensitivity. Patients with HLA-DRB1*07:01 and HLA-DQB1*02:02 alleles had significantly higher risk of developing asparaginase hypersensitivity compared to non-carriers [P=4.56×10-5; OR=2.86 (1.73-4.75) and P=1.85×10-4; OR=2.99 (1.68-5.31); n=359, respectively]. After haplotype reconstruction, the HLA-DRB1*07:01-HLA-DQB1*02:02 haplotype was associated with an increased risk. After inferring the HLA-DQA1 alleles the HLA-DRB1*07:01-HLA-DQA1*02:01-HLA-DQB1*02:02 haplotype was associated with the highest risk of asparaginase hypersensitivity [P=1.22×10-5; OR=5.00 (2.43-10.29); n=257]. Significantly fewer T-cell ALL patients carried the HLA-DQB1*02:02 allele and the associated haplotype than did pre-B-cell ALL patients (6.5%; vs. 19.2%, respectively; P=0.047). In conclusion, we identified a haplotype in the Human Leukocyte Antigen Class II region associated with a higher risk of asparaginase hypersensitivity. Our results confirm that variations in HLA-D region might influence the development of asparaginase hypersensitivity.

Investigation of the Possible Role of the Hippo/YAP1 Pathway in Asthma and Allergy.

PURPOSE: Several lines of evidence indicate that the Hippo/Yes-associated protein 1 (YAP1) pathways might play a role in the pathogenesis of asthma. To investigate the possible role of the Hippo/YAP1 pathway in the pathogenesis of asthma or its phenotypes. METHODS: The levels of gene expressions of the members of the Hippo/YAP1 were compared. The presence of the proteins of the YAP1 and FRMD6 were analyzed with Western blot in induced sputum of 18 asthmatic subjects and 10 control subjects. Fourteen single nucleotide polymorphisms (SNPs) in the YAP1 gene were genotyped in 522 asthmatic subjects and 711 healthy controls. The results were evaluated with traditional frequentist methods and with Bayesian network-based Bayesian multilevel analysis of relevance (BN-BMLA). RESULTS: The mRNA of all the members of the Hippo/YAP1 pathway could be detected in the induced sputum of both controls and cases. A correlation was found between YAP1 mRNA levels and sputum bronchial epithelial cells (r=0.575, P=0.003). The signal for the FRMD6 protein could be detected in all sputum samples while the YAP1 protein could not be detected in the sputum samples, of the healthy controls and severe asthmatics, but it was detectable in mild asthmatics. The rs2846836 SNP of the YAP1 gene was significantly associated with exercise-induced asthma (odds ratio [OR]=2.1 [1.3-3.4]; P=0.004). The distribution of genotypes of rs11225138 and certain haplotypes of the YAP1 gene showed significant differences between different asthma severity statuses. With BN-BMLA, 2 SNPs, genetic variations in the FRMD6 gene proved to be the most relevant to exercise-induced asthma and allergic rhinitis. These 2 SNPs through allergic rhinitis and exercise-induced asthma were in epistatic interaction with each other. CONCLUSIONS: Our results provided additional evidence that the FRMD6/Hippo/YAP1 pathway plays a role in the pathogenesis of asthma. If additional studies can confirm these findings, this pathway can be a potential novel therapeutic target in asthma and other inflammatory airway diseases.

Pharmacogenetic analysis of high-dose methotrexate treatment in children with osteosarcoma.

Inter-individual differences in toxic symptoms and pharmacokinetics of high-dose methotrexate (MTX) treatment may be caused by genetic variants in the MTX pathway. Correlations between polymorphisms and pharmacokinetic parameters and the occurrence of hepato- and myelotoxicity were studied. Single nucleotide polymorphisms (SNPs) of the ABCB1, ABCC1, ABCC2, ABCC3, ABCC10, ABCG2, GGH, SLC19A1 and NR1I2 genes were analyzed in 59 patients with osteosarcoma. Univariate association analysis and Bayesian network-based Bayesian univariate and multilevel analysis of relevance (BN-BMLA) were applied. Rare alleles of 10 SNPs of ABCB1, ABCC2, ABCC3, ABCG2 and NR1I2 genes showed a correlation with the pharmacokinetic values and univariate association analysis. The risk of toxicity was associated with five SNPs in the ABCC2 and NR1I2 genes. Pharmacokinetic parameters were associated with four SNPs of the ABCB1, ABCC3, NR1I2, and GGH genes, and toxicity was shown to be associated with ABCC1 rs246219 and ABCC2 rs717620 using the univariate and BN-BMLA method. BN-BMLA analysis detected relevant effects on the AUC0-48 in the following SNPs: ABCB1 rs928256, ABCC3 rs4793665, and GGH rs3758149. In both univariate and multivariate analyses the SNPs ABCB1 rs928256, ABCC3 rs4793665, GGH rs3758149, and NR1I2 rs3814058 SNPs were relevant. These SNPs should be considered in future dose individualization during treatment.

Pharmacogenetics of anthracyclines.

Anthracyclines constitute a fundamental part of the chemotherapy regimens utilized to treat a number of different malignancies both in pediatric and adult patients. These drugs are one of the most efficacious anticancer agents ever invented. On the other hand, anthracyclines are cardiotoxic. Childhood cancer survivors treated with anthracyclines often undergo cardiac complications which are influenced by genetic variations of the patients. The scientific literature comprises numerous investigations in the subject of the pharmacogenetics of anthracyclines. In this review, we provide a comprehensive overview of this research topic. Genetic variants are proposed targets in the personalized treatment in order to individualize dosing and therefore reduce side effects.

Subgroups of Paediatric Acute Lymphoblastic Leukaemia Might Differ Significantly in Genetic Predisposition to Asparaginase Hypersensitivity.

L-asparaginase (ASP) is a key element in the treatment of paediatric acute lymphoblastic leukaemia (ALL). However, hypersensitivity reactions (HSRs) to ASP are major challenges in paediatric patients. Our aim was to investigate genetic variants that may influence the risk to Escherichia coli-derived ASP hypersensitivity. Sample and clinical data collection was carried out from 576 paediatric ALL patients who were treated according to protocols from the Berlin-Frankfurt-Münster Study Group. A total of 20 single nucleotide polymorphisms (SNPs) in GRIA1 and GALNT10 genes were genotyped. Patients with GRIA1 rs4958351 AA/AG genotype showed significantly reduced risk to ASP hypersensitivity compared to patients with GG genotype in the T-cell ALL subgroup (OR = 0.05 (0.01-0.26); p = 4.70E-04), while no such association was found in pre-B-cell ALL. In the medium risk group two SNPs of GRIA1 (rs2055083 and rs707176) were associated significantly with the occurrence of ASP hypersensitivity (OR = 0.21 (0.09-0.53); p = 8.48E-04 and OR = 3.02 (1.36-6.73); p = 6.76E-03, respectively). Evaluating the genders separately, however, the association of rs707176 with ASP HSRs was confined only to females. Our results suggest that genetic variants of GRIA1 might influence the risk to ASP hypersensitivity, but subgroups of patients can differ significantly in this respect.

Impact of single nucleotide polymorphisms of cytarabine metabolic genes on drug toxicity in childhood acute lymphoblastic leukemia.

BACKGROUND: Cytarabine (cytosine arabinoside, ara-C) is a chemotherapeutical agent used in the treatment of pediatric acute lymphoblastic leukemia (ALL). Adverse drug reactions, such as interpatient variability in sensitivity to ara-C, are considerable and may cause difficulties during chemotherapy. Single nucleotide polymorphisms (SNPs) can play a significant role in modifying nucleoside-drug pharmacokinetics and pharmacodynamics and thus the development of adverse effects. Our aim was to determine whether polymorphisms in genes encoding transporters and enzymes responsible for the metabolism of ara-C are associated with toxicity and clinical outcome in a patient population with childhood ALL. PROCEDURE: We studied 8 SNPs in the CDA, DCK, DCTD, SLC28A3, and SLC29A1 genes in 144 patients with childhood acute lymphoblastic leukemia treated according to ALLIC BFM 1990, 1995 and 2002 protocols. RESULTS: DCK rs12648166 and DCK rs4694362 SNPs were associated with hematologic toxicity (OR = 2.63, CI 95% = 1.37-5.04, P = 0.0036 and OR = 2.53, CI 95% = 1.34-4.80, P = 0.0044, respectively). CONCLUSIONS: Our results indicate that DCK polymorphisms might be important genetic risk factors for hematologic toxicity during ALL treatment with ara-C. Individualized chemotherapy based on genetic profiling may help to optimize ara-C dosing, leading to improvements in clinical outcome and reduced toxicity.

Associations of novel genetic variations in the folate-related and ARID5B genes with the pharmacokinetics and toxicity of high-dose methotrexate in paediatric acute lymphoblastic leukaemia.

High-dose methotrexate (HD-MTX) plays an important role in the consolidation therapy of acute lymphoblastic leukaemia (ALL) in many treatment regimens worldwide. However, there is a large interpatient variability in the pharmacokinetics and toxicity of the drug. We investigated the influence of single nucleotide polymorphisms (SNPs) in genes of the folate metabolic pathway, transporter molecules and transcription proteins on the pharmacokinetics and toxicity of MTX and 7-hydroxy-methotrexate (7-OH-MTX). 63 SNPs of 14 genes were genotyped and a total of 463 HD-MTX courses (administered according to the ALL-BFM 95 and ALL IC-BFM 2002 protocols) were analysed. Haematological, hepatic and renal toxicities, estimated by routine laboratory parameters were evaluated. Random forest and regression trees were used for variable selection and model building. Linear mixed models were established to prove the significance of the selected variables. SNPs (rs4948502, rs4948496, rs4948487) of the ARID5B gene were associated with the serum levels of MTX (P < 0·02), serum levels and area under the curve of 7-OH-MTX (P < 0·02) and with hypoproteinaemia (P = 0·004). SLCO1B1 rs4149056 also showed a significant association with serum MTX levels (P < 0·001). Our findings confirm the association of novel genetic variations in folate-related and ARID5B genes with the serum MTX levels and acute toxicity.

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(-4); 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(-3); 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.

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.

Non-synonymous single nucleotide polymorphisms in genes for immunoregulatory galectins: association of galectin-8 (F19Y) occurrence with autoimmune diseases in a Caucasian population.

BACKGROUND: Galectins are potent immune regulators, with galectin-8 acting as a pro-apoptotic effector on synovial fluid cells and thymocytes and stimulator on T-cells. To set a proof-of-principle example for risk assessment in autoimmunity, and for a mutation affecting physiological galectin sensor functions, a polymorphism in the coding region of the galectin-8 gene (rs2737713; F19Y) was studied for its association with two autoimmune disorders, i.e. rheumatoid arthritis and myasthenia gravis. METHODS: A case-control analysis and a related quantitative trait-association study were performed to investigate the association of this polymorphism in patients (myasthenia gravis 149, rheumatoid arthritis 214 and 134 as primary and repetitive cohorts, respectively) and 365 ethnically matched (Caucasian) healthy controls. Distribution was also investigated in patients grouped according to their antibody status and age at disease onset. Comparative testing for lectin activity was carried out in ELISA/ELLA-based binding tests with both wild-type and F19Y mutant galectin-8 from peripheral blood mononuclear cell lysates of healthy individuals with different genotypes as well as with recombinant wild-type and F19Y mutant galectin-8 proteins. RESULTS: A strong association was found for rheumatoid arthritis, and a mild one with myasthenia gravis. Furthermore, the presence of the sequence deviation also correlated with age at disease onset in the case of rheumatoid arthritis. The F19Y substitution did not appear to affect carbohydrate binding in solid-phase assays markedly. GENERAL SIGNIFICANCE: This is the first report of an association between a galectin-based polymorphism leading to a mutant protein and autoimmune diseases, with evidence for antagonistic pleiotropy.

Asthma endophenotypes and polymorphisms in the histamine receptor HRH4 gene.

BACKGROUND: Histamine as an inflammatory mediator plays an important role in chronic allergic and asthmatic conditions. However, the role of genetic polymorphisms of the histamine receptor HRH4 (histamine receptor H4) gene in asthma susceptibility and endophenotypes has not been studied yet. Our aim was to investigate the possible association between single nucleotide polymorphisms (SNPs) in the HRH4 gene and asthma or some endophenotypes of asthma. METHODS: Twenty-one SNPs of the HRH4 gene were genotyped in 313 asthmatic patients and 360 controls using Sequenom® iPLEX® Gold Genotyping Technology. RESULTS: Genotype distribution of three HRH4 SNPs, namely rs17187619 [p = 0.002; odds ratio, OR (95% confidence interval, CI) = 2.4 (4.1-1.4)], rs527790 [p = 0.0002; OR (95% CI) = 3.3 (6.1-1.8)] and rs487202 [p = 0.00007; OR (95% CI) = 3.5 (6.6-1.9)] differed significantly between patients with or without infection-induced asthma. Haplotypes, which included the rs4800573-rs527790 CC allele combination, were found to be associated with infection-induced asthma [p = 0.0009, OR (95% CI) = 0.5 (0.4-0.8)]. The rs487202-rs574913 CA haplotype was more frequent among patients with infection-induced asthma [p = 0.0006, OR (95% CI) = 1.9 (1.3-2.6)]. None of the SNPs contributed directly to the risk of asthma. CONCLUSIONS: Our results suggest that genetic variation in the HRH4 gene might influence the pathogenesis of infection-induced asthma.

Evaluation of a partial genome screening of two asthma susceptibility regions using bayesian network based bayesian multilevel analysis of relevance.

Genetic studies indicate high number of potential factors related to asthma. Based on earlier linkage analyses we selected the 11q13 and 14q22 asthma susceptibility regions, for which we designed a partial genome screening study using 145 SNPs in 1201 individuals (436 asthmatic children and 765 controls). The results were evaluated with traditional frequentist methods and we applied a new statistical method, called bayesian network based bayesian multilevel analysis of relevance (BN-BMLA). This method uses bayesian network representation to provide detailed characterization of the relevance of factors, such as joint significance, the type of dependency, and multi-target aspects. We estimated posteriors for these relations within the bayesian statistical framework, in order to estimate the posteriors whether a variable is directly relevant or its association is only mediated.With frequentist methods one SNP (rs3751464 in the FRMD6 gene) provided evidence for an association with asthma (OR = 1.43(1.2-1.8); p = 3×10(-4)). The possible role of the FRMD6 gene in asthma was also confirmed in an animal model and human asthmatics.In the BN-BMLA analysis altogether 5 SNPs in 4 genes were found relevant in connection with asthma phenotype: PRPF19 on chromosome 11, and FRMD6, PTGER2 and PTGDR on chromosome 14. In a subsequent step a partial dataset containing rhinitis and further clinical parameters was used, which allowed the analysis of relevance of SNPs for asthma and multiple targets. These analyses suggested that SNPs in the AHNAK and MS4A2 genes were indirectly associated with asthma. This paper indicates that BN-BMLA explores the relevant factors more comprehensively than traditional statistical methods and extends the scope of strong relevance based methods to include partial relevance, global characterization of relevance and multi-target relevance.