Genetic markers, birth characteristics, and childhood leukemia risk

The cause for childhood acute lymphoblastic leukemia (ALL) remains unknown, but male gender is a risk factor, and among ethnicities, Hispanics have the highest risk. In this dissertation, we explored correlations among genetic polymorphisms, birth characteristics, and the risk of childhood ALL in a multi-ethnic sample in 161 cases and 231 controls recruited contemporaneously (2007-2012) in Houston, TX. We first examined three lymphoma risk markers, since lymphoma and ALL both stem from lymphoid cells. Of these, rs2395185 showed a risk association in non-Hispanic White males (OR=2.8, P=0.02; P interaction=0.03 for gender), but not in Hispanics. We verified previously known risk associations to validate the case-control sample. Mutations of HFE (C282Y, H63D) were genotyped to test whether iron-regulatory gene (IRG) variants known to elevate iron levels increase childhood ALL risk. Being positive for either polymorphism yielded only a modestly elevated OR in males, which increased to 2.96 (P=0.01) in the presence of a particular transferrin receptor (TFRC) genotype for rs3817672 (Pinteraction=0.04). SNP rs3817672 itself showed an ethnicity-specific association (P interaction=0.02 for ethnicity). We then examined additional IRG SNPs (rs422982, rs855791, rs733655), which showed risk associations in males (ORs=1.52 to 2.60). A polygenic model based on the number of polymorphic alleles in five IRG SNPs revealed a linear increase in risk (OR=2.00 per incremental change; P=0.002). Having three or more alleles compared with none was associated with increased risk in males (OR=4.12; P=0.004). Significant risk associations with childhood ALL was found with birth length (OR=1.18 per inch, P=0.04), high birth weight (>4,000g) (OR=1.93, P=0.01), and with gestational age (OR=1.10 per week, P=0.04). We observed a negative correlation between HFE SNP rs9366637 and gestational age (P=0.005), again, stronger in males ( P=0.001) and interacting with TFRC (P interaction=0.05). Our results showed that (i) ALL risk markers do not show universal associations across ethnicities or between genders, (ii) IRG SNPs modify ALL risk presumably by their effects on iron levels, (iii) a negative correlation between an HFE SNP and gestational age exists, which implicates an iron-related mechanism. The results suggest that currently unregulated supplemental iron intake may have implications on childhood ALL development.

A Case-Only Genome-wide Association Study of Gender- and Age-specific Risk Markers for Childhood Leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.

A case-only genome-wide association study of gender- and age-specific risk markers for childhood leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. ^ Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5 ) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. ^ The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.^

Interactions of different arsenic species with thols: Chemical and biological implications

Arsenic has been classified as a group I carcinogen. It has been ranked number one in the CERCLA priority list of hazardous substances due to its frequency, toxicity and potential for human exposure. Paradoxically, arsenic has been employed as a successful chemotherapeutic agent for acute promyelocytic leukemia and has found some success in multiple myeloma. Since arsenic toxicity and efficacy is species dependent, a speciation method, based on the complementary use of reverse phase and cation exchange chromatography, was developed. Inductively coupled plasma mass spectrometer (ICP-MS), as an element specific detector, and electrospray ionization mass spectrometer (ESI-MS), as a molecule specific detector, were employed. Low detection limits in the µg. L−1 range on the ICP-MS and mg. L−1 range on the ESI-MS were obtained. The developed methods were validated against each other through the use of a Deming plot. With the developed speciation method, the effects of both pH on the stability of As species and reduced glutathione (GSH) concentration on the formation and stability of arsenic glutathione complexes were studied. To identify arsenicals in multiple myeloma (MM) cell lines post arsenic trioxide (ATO) and darinaparsin (DAR) incubation, an extraction method based on the use of ultrasonic probe was developed. Extraction tools and solvents were evaluated and the effect of GSH concentration on the quantitation of arsenic glutathione (As-GSH) complexes in MM cell extracts was studied. The developed method was employed for the identification of metabolites in DAR incubated cell lines where the effect of extraction pH, DAR incubation concentration and incubation time on the relative distribution of the As metabolites was assessed. A new arsenic species, dimethyarsinothioyl glutathione (DMMTA V-GS), a pentavalent thiolated arsenical, was identified in the cell extracts through the use of liquid chromatography tandem mass spectrometry. The formation of the new metabolite in the extracts was dependent on the decomposition of s-dimethylarsino glutathione (DMA(GS)). These results have major implications in both the medical and toxicological fields of As because they involve the metabolism of a chemotherapeutic agent and the role sulfur compounds play in this mechanism.

Towards the prediction of mutations in genomic sequences

Bio-systems are inherently complex information processing systems. Furthermore, physiological complexities of biological systems limit the formation of a hypothesis in terms of behavior and the ability to test hypothesis. More importantly the identification and classification of mutation in patients are centric topics in today's cancer research. Next generation sequencing (NGS) technologies can provide genome-wide coverage at a single nucleotide resolution and at reasonable speed and cost. The unprecedented molecular characterization provided by NGS offers the potential for an individualized approach to treatment. These advances in cancer genomics have enabled scientists to interrogate cancer-specific genomic variants and compare them with the normal variants in the same patient. Analysis of this data provides a catalog of somatic variants, present in tumor genome but not in the normal tissue DNA. In this dissertation, we present a new computational framework to the problem of predicting the number of mutations on a chromosome for a certain patient, which is a fundamental problem in clinical and research fields. We begin this dissertation with the development of a framework system that is capable of utilizing published data from a longitudinal study of patients with acute myeloid leukemia (AML), who's DNA from both normal as well as malignant tissues was subjected to NGS analysis at various points in time. By processing the sequencing data at the time of cancer diagnosis using the components of our framework, we tested it by predicting the genomic regions to be mutated at the time of relapse and, later, by comparing our results with the actual regions that showed mutations (discovered at relapse time). We demonstrate that this coupling of the algorithm pipeline can drastically improve the predictive abilities of searching a reliable molecular signature. Arguably, the most important result of our research is its superior performance to other methods like Radial Basis Function Network, Sequential Minimal Optimization, and Gaussian Process. In the final part of this dissertation, we present a detailed significance, stability and statistical analysis of our model. A performance comparison of the results are presented. This work clearly lays a good foundation for future research for other types of cancer.^