Extensions to Gene Set Enrichment

Motivation: Gene Set Enrichment Analysis (GSEA) has been developed recently to capture moderate but coordinated changes in the expression of sets of functionally related genes. We propose number of extensions to GSEA, which uses different statistics to describe the association between genes and phenotype of interest. We make use of dimension reduction procedures, such as principle component analysis to identify gene sets containing coordinated genes. We also address the problem of overlapping among gene sets in this paper. Results: We applied our methods to the data come from a clinical trial in acute lymphoblastic leukemia (ALL) [1]. We identified interesting gene sets using different statistics. We find that gender may have effects on the gene expression in addition to the phenotype effects. Investigating overlap among interesting gene sets indicate that overlapping could alter the interpretation of the significant results.

Associations between vaccination and childhood cancers in Texas counties

Background. In Dr. Mel Greaves "delayed-infection hypothesis," postponed exposure to common infections increases the likelihood of childhood cancer. Hygienic advancements in developed countries have reduced children's exposure to pathogens and children encounter common infectious agents at an older age with an immune system unable to deal with the foreign antigens. Vaccinations may be considered to be simulated infections as they prompt an antigenic response by the immune system. Vaccinations may regulate the risk of childhood cancer by modulating the immune system. The aim of the study was to determine if children born in Texas counties with higher levels of vaccination coverage were at a reduced risk for childhood cancer.^ Methods. We conducted a case-control study to examine the risk of childhood cancers, specifically leukemia, brain tumors, and non-Hodgkin lymphoma, in relation to vaccination rates in Texas counties. We utilized a multilevel mixed-effects regression model of the individual data from the Texas Cancer Registry (TCR) with group-level exposure data (i.e., the county- and public health region-level vaccination rates).^ Results. Utilizing county-level vaccination rates and controlling for child's sex, birth year, ethnicity, birth weight, and mother's age at child's birth the hepatitis B vaccine revealed negative associations with developing all cancer types (OR = 0.81, 95% CI: 0.67–0.98) and acute lymphoblastic leukemia (ALL) (OR = 0.63, 95% CI: 0.46–0.88). The decreased risk for ALL was also evident for the inactivated polio vaccine (IPV) (OR = 0.67, 95% CI: 0.49–0.92) and 4-3-1-3-3 vaccination series (OR = 0.62, 95% CI: 0.44-0.87). Using public health region vaccine coverage levels, an inverse association between the Haemophilus influenzae type b (Hib) vaccine and ALL (OR: 0.58; 95% CI: 0.42–0.82) was present. Conversely, the measles, mumps, and rubella (MMR) vaccine resulted in a positive association with developing non-Hodgkin lymphoma (OR = 2.81, 95% CI: 1.27–6.22). ^

Exploring ant-based algorithms for gene expression data analysis

Objective: Recently, much research has been proposed using nature inspired algorithms to perform complex machine learning tasks. Ant colony optimization (ACO) is one such algorithm based on swarm intelligence and is derived from a model inspired by the collective foraging behavior of ants. Taking advantage of the ACO in traits such as self-organization and robustness, this paper investigates ant-based algorithms for gene expression data clustering and associative classification. Methods and material: An ant-based clustering (Ant-C) and an ant-based association rule mining (Ant-ARM) algorithms are proposed for gene expression data analysis. The proposed algorithms make use of the natural behavior of ants such as cooperation and adaptation to allow for a flexible robust search for a good candidate solution. Results: Ant-C has been tested on the three datasets selected from the Stanford Genomic Resource Database and achieved relatively high accuracy compared to other classical clustering methods. Ant-ARM has been tested on the acute lymphoblastic leukemia (ALL)/acute myeloid leukemia (AML) dataset and generated about 30 classification rules with high accuracy. Conclusions: Ant-C can generate optimal number of clusters without incorporating any other algorithms such as K-means or agglomerative hierarchical clustering. For associative classification, while a few of the well-known algorithms such as Apriori, FP-growth and Magnum Opus are unable to mine any association rules from the ALL/AML dataset within a reasonable period of time, Ant-ARM is able to extract associative classification rules.

Impact of genetic polymorphisms determining leukocyte/neutrophil count on chemotherapy toxicity

Nous avons investigué la relation entre les polymorphismes de nucléotides simples (SNPs) chez trois gènes/loci candidats : DARC, CXCL2 et le loci ORMDL3-GSDMA-CSF3 situés sur le chromosome 17q21 et les complications neutropéniques et infectieuses qui en résultent durant la chimiothérapie chez les patients atteints de la leucémie lymphoblastique aigue. Ces loci codent pour certaines composantes du système immunitaire altérant la concentration de chémokines et leur distribution (DARC), stimulant le relâchement et la migration des neutophiles de la moelle épinière (CXCL2) et régulant la prolifération et la survie des granulocytes (G-CSF). Il est possible que des polymorphismes dans ces loci lorsqu’associés à de la chimiothérapie puissent mettre des individus suceptibles à un risque plus élevé de complication reliées à la chimiothérapie. Une sélection des marqueurs SNPs dans ces gènes ont été génotypés chez des enfants traités au CHU Ste-Justine pour une ALL entre 1989 et 2005. Après correction pour tests multiples, un polymorphisme DARC rs3027012 situé dans le 5’UTR a été associé à un compte phagocytaire peu élevé (APC<500 et <1000 cellules/µL, p=0.001 and p=0.0005, respectivement) ainsi qu’une hospitalisation due à une neutropénie (p=0.007) ou due à une infection et/ou neutropénie (p=0.007). Un effet protecteur a été identifié pour la mutation non sense Gly42Asp variant rs12075 (p=0.006). Des polymorphismes sur le chromosome 17q2 étaient associés à une hospitalisation due à une infection (rs3859192, p= 0.004) et à une neutropénie (rs17609240, p=0.006) L’infection était aussi modulée par CXCL2 (rs16850408, p=0.008) Cette étude identifie pour la première fois que les loci modulant le décompte des leucocytes et des neutrophiles pourraient jouer un rôle dans de déclenchement de complications dues à la chimiothérapie et pourraient ainsi servir de marqueurs pour un ajustement et un suivi du traitement.

Impact of genetic polymorphisms determining leukocyte/neutrophil count on chemotherapy toxicity

Nous avons investigué la relation entre les polymorphismes de nucléotides simples (SNPs) chez trois gènes/loci candidats : DARC, CXCL2 et le loci ORMDL3-GSDMA-CSF3 situés sur le chromosome 17q21 et les complications neutropéniques et infectieuses qui en résultent durant la chimiothérapie chez les patients atteints de la leucémie lymphoblastique aigue. Ces loci codent pour certaines composantes du système immunitaire altérant la concentration de chémokines et leur distribution (DARC), stimulant le relâchement et la migration des neutophiles de la moelle épinière (CXCL2) et régulant la prolifération et la survie des granulocytes (G-CSF). Il est possible que des polymorphismes dans ces loci lorsqu’associés à de la chimiothérapie puissent mettre des individus suceptibles à un risque plus élevé de complication reliées à la chimiothérapie. Une sélection des marqueurs SNPs dans ces gènes ont été génotypés chez des enfants traités au CHU Ste-Justine pour une ALL entre 1989 et 2005. Après correction pour tests multiples, un polymorphisme DARC rs3027012 situé dans le 5’UTR a été associé à un compte phagocytaire peu élevé (APC<500 et <1000 cellules/µL, p=0.001 and p=0.0005, respectivement) ainsi qu’une hospitalisation due à une neutropénie (p=0.007) ou due à une infection et/ou neutropénie (p=0.007). Un effet protecteur a été identifié pour la mutation non sense Gly42Asp variant rs12075 (p=0.006). Des polymorphismes sur le chromosome 17q2 étaient associés à une hospitalisation due à une infection (rs3859192, p= 0.004) et à une neutropénie (rs17609240, p=0.006) L’infection était aussi modulée par CXCL2 (rs16850408, p=0.008) Cette étude identifie pour la première fois que les loci modulant le décompte des leucocytes et des neutrophiles pourraient jouer un rôle dans de déclenchement de complications dues à la chimiothérapie et pourraient ainsi servir de marqueurs pour un ajustement et un suivi du traitement.

Central nervous system involvement at first relapse in patients with acute promyelocytic leukemia treated with all-trans retinoic acid and anthracycline monochemotherapy without intrathecal prophylaxis.

BACKGROUND The prevalence of and risk factors for central nervous system recurrence in patients with acute promyelocytic leukemia are not well established and remain a controversial matter. DESIGN AND METHODS Between 1996 and 2005, 739 patients with newly diagnosed acute promyelocytic leukemia enrolled in two consecutive trials (PETHEMA LPA96 and LPA99) received induction therapy with all-trans retinoic acid and idarubicin. Consolidation therapy comprised three courses of anthracycline monochemotherapy (LPA96), with all-trans retinoic acid and reinforced doses of idarubicin in patients with an intermediate or high risk of relapse (LPA99). Central nervous system prophylaxis was not given. RESULTS Central nervous system relapse was documented in 11 patients. The 5-year cumulative incidence of central nervous system relapse was 1.7% (LPA96 3.2% and LPA99 1.2%; p=0.09). The cumulative incidence was 0%, 0.8%, and 5.5% in low-, intermediate-, and high-risk patients, respectively. Relapse risk score (p=0.0001) and the occurrence of central nervous system hemorrhage during induction (5-year cumulative incidence 18.7%, p=0.006) were independent risk factors for central nervous system relapse. CONCLUSIONS This study shows a low incidence of central nervous system relapse in patients with acute promyelocytic leukemia following therapy with all-trans retinoic acid and anthracycline without specific central nervous system prophylaxis. Central nervous system relapse was significantly associated with high white blood cell counts and prior central nervous system hemorrhage, which emerged as independent prognostic factors.

A non-random chromosome abnormality found in precursor-B lineage acute lymphoblastic leukaemia: dic(9;20)(p1?3;q11).

A comparison of cytogenetical data on acute lymphoblastic leukaemia studied at four large European centres has revealed a non-random dicentric chromosome abnormality: dic(9;20) (p1?3;q11) in 10 patients, nine of whom were children. All had early precursor-B lineage ALL, and eight children had a non-standard risk clinical presentation. The origin of the dicentric chromosome was demonstrated using a range of chromosome banding techniques. This was confirmed by FISH using paints and centromeric probes for chromosomes 9 and 20, together with a number of cosmid probes. The follow-up time of these patients is presently too short and the number of patients too few to determine the prognostic significant of this chromosome abnormality.

Effect of all-trans retinoic acid on newly diagnosed acute promyelocytic leukemia patients: results of a Brazilian center

Thirty-seven patients with acute promyelocytic leukemia (APL) were treated with all-trans retinoic acid (ATRA). Patients received 45 mg m-2 day-1 po of ATRA until complete remission (CR) was achieved, defined as: a) presence of less than 5% blasts in the bone marrow, with b) white blood cells >103/mm3, c) platelets >105/mm3 and d) hemoglobin concentration >8 g/dl, with no blood or platelet transfusions. Thirty-one (83.7%) patients achieved CR by day 50, and 75% of these before day 30. Correction of the coagulopathy, achieved between days 2 and 10 (mean, 3 days), was the first evidence of response to treatment. Only one patient had been previously treated with chemotherapy and three had the microgranular variant M3 form. Dryness of skin and mucosae was the most common side effect observed in 82% of the patients. Thrombosis, hepatotoxicity and retinoid acid syndrome (RAS) were observed in 7 (19%), 6 (16%) and 4 (11%) patients, respectively. Thirteen (35%) patients had to be submitted to chemotherapy due to hyperleukocytosis (above 40 x 103/mm3) and six of these presented with new signs of coagulopathy after chemotherapy. Four (11%) patients died secondarily to intracerebral hemorrhage (IH) and two (5.4%) dropped out of the protocol due to severe ATRA side effects (one RAS and one hepatotoxicity). RAS and IH were related strictly to hyperleukocytosis. The reduced use of platelets and fresh frozen plasma probably lowered the total cost of treatment. We conclude that ATRA is an effective agent for inducing complete remission in APL patients.

Mechanism of resistance to tyrosine kinase inhibitors in philadelphia-positive acute lymphblastic leukaemia (all): from genetic alterations to impaired RNA editing

The Ph chromosome is the most frequent cytogenetic aberration associated with adult ALL and it represents the single most significant adverse prognostic marker. Despite imatinib has led to significant improvements in the treatment of patients with Ph+ ALL, in the majority of cases resistance developed quickly and disease progressed. Some mechanisms of resistance have been widely described but the full knowledge of contributing factors, driving both the disease and resistance, remains to be defined. The observation of rapid development of lymphoblastic leukemia in mice expressing altered Ikaros (Ik) isoforms represented the background of this study. Ikaros is a zinc finger transcription factor required for normal hemopoietic differentiation and proliferation, particularly in the lymphoid lineages. By means of alternative splicing, Ikaros encodes several proteins that differ in their abilities to bind to a consensus DNA-binding site. Shorter, DNA nonbinding isoforms exert a dominant negative effect, inhibiting the ability of longer heterodimer partners to bind DNA. The differential expression pattern of Ik isoforms in Ph+ ALL patients was analyzed in order to determine if molecular abnormalities involving the Ik gene could associate with resistance to imatinib and dasatinib. Bone marrow and peripheral blood samples from 46 adult patients (median age 55 yrs, 18-76) with Ph+ ALL at diagnosis and during treatment with imatinib (16 pts) or dasatinib (30 pts) were collected. We set up a fast, high-throughput method based on capillary electrophoresis technology to detect and quantify splice variants. 41% Ph+ ALL patients expressed high levels of the non DNA-binding dominant negative Ik6 isoform lacking critical N-terminal zinc-fingers which display abnormal subcellular compartmentalization pattern. Nuclear extracts from patients expressed Ik6 failed to bind DNA in mobility shift assay using a DNA probe containing an Ikaros-specific DNA binding sequence. In 59% Ph+ ALL patients there was the coexistence in the same PCR sample and at the same time of many splice variants corresponded to Ik1, Ik2, Ik4, Ik4A, Ik5A, Ik6, Ik6 and Ik8 isoforms. In these patients aberrant full-length Ikaros isoforms in Ph+ ALL characterized by a 60-bp insertion immediately downstream of exon 3 and a recurring 30-bp in-frame deletion at the end of exon 7 involving most frequently the Ik2, Ik4 isoforms were also identified. Both the insertion and deletion were due to the selection of alternative splice donor and acceptor sites. The molecular monitoring of minimal residual disease showed for the first time in vivo that the Ik6 expression strongly correlated with the BCR-ABL transcript levels suggesting that this alteration could depend on the Bcr-Abl activity. Patient-derived leukaemia cells expressed dominant-negative Ik6 at diagnosis and at the time of relapse, but never during remission. In order to mechanistically demonstrated whether in vitro the overexpression of Ik6 impairs the response to tyrosine kinase inhibitors (TKIs) and contributes to resistance, an imatinib-sensitive Ik6-negative Ph+ ALL cell line (SUP-B15) was transfected with the complete Ik6 DNA coding sequence. The expression of Ik6 strongly increased proliferation and inhibited apoptosis in TKI sensitive cells establishing a previously unknown link between specific molecular defects that involve the Ikaros gene and the resistance to TKIs in Ph+ ALL patients. Amplification and genomic sequence analysis of the exon splice junction regions showed the presence of 2 single nucleotide polymorphisms (SNPs): rs10251980 [A/G] in the exon2/3 splice junction and of rs10262731 [A/G] in the exon 7/8 splice junction in 50% and 36% of patients, respectively. A variant of the rs11329346 [-/C], in 16% of patients was also found. Other two different single nucleotide substitutions not recognized as SNP were observed. Some mutations were predicted by computational analyses (RESCUE approach) to alter cis-splicing elements. In conclusion, these findings demonstrated that the post-transcriptional regulation of alternative splicing of Ikaros gene is defective in the majority of Ph+ ALL patients treated with TKIs. The overexpression of Ik6 blocking B-cell differentiation could contribute to resistance opening a time frame, during which leukaemia cells acquire secondary transforming events that confer definitive resistance to imatinib and dasatinib.

Extramedullary relapse in acute promyelocytic leukemia treated with all-trans retinoic acid and chemotherapy

We analyzed the incidence, presenting features, risk factors of extramedullary (EM) relapse occurring in acute promyelocytic leukemia (APL) treated with all-trans retinoic acid (ATRA) and chemotherapy by using a competing-risk method. In total, 740/ 806 (92%) patients included in three multicenter trials (APL91, APL93 trials and PETHEMA 96) achieved CR, of whom 169 (23%) relapsed, including 10 EM relapses. Nine relapses involved the central nervous system (CNS) and one the skin, of which two were isolated EM relapse. In patients with EM disease, median WBC count was 26950/mm3 (7700-162000). The 3-year cumulative incidence of EM disease at first relapse was 5.0%. Univariate analysis identified age <45 years (P=0.05), bcr3 PML-RARalpha isoform (P= 0.0003) and high WBC counts (> or = 10,000/ mm3) (P<0.0001) as risk factors for EM relapse. In multivariate analysis, only high WBC count remained significant (P= 0.001). Patients with EM relapse had a poorer outcome since median survival from EM relapse was 6.7 months as compared to 26.3 months for isolated BM relapse (P=0.04). In conclusion, EM relapse in APL occurs more frequently in patients with increased WBC counts (> or = 10,000/mm3) and carries a poor prognosis. Whether CNS prophylaxis should be systematically performed in patients with WBC > or = 10,000/mm3 at diagnosis remains to be established.

Venous thromboembolism in patients with acute leukemia: Prevalence, incidence, recurrence of thrombosis, risk factors, and impact on survival

Venous thromboembolism (VTE), including deep vein thrombosis (DVT) and pulmonary embolism (PE), is the third most preventable cardiovascular disease and a growing public health problem in the United States. The incidence of VTE remains high with an annual estimate of more than 600,000 symptomatic events. DVT affects an estimated 2 million American each year with a death toll of 300,000 persons per year from DVT-related PE. Leukemia patients are at high risk for both hemorrhage and thrombosis; however, little is known about thrombosis among acute leukemia patients. The ultimate goal of this dissertation was to obtain deep understanding of thrombotic issue among acute leukemia patients. The dissertation was presented in a format of three papers. First paper mainly looked at distribution and risk factors associated with development of VTE among patients with acute leukemia prior to leukemia treatment. Second paper looked at incidence, risk factors, and impact of VTE on survival of patients with acute lymphoblastic leukemia during treatment. Third paper looked at recurrence and risk factors for VTE recurrence among acute leukemia patients with an initial episode of VTE. Descriptive statistics, Chi-squared or Fisher's exact test, median test, Mann-Whitney test, logistic regression analysis, Nonparametric Estimation Kaplan-Meier with a log-rank test or Cox model were used when appropriate. Results from analyses indicated that acute leukemia patients had a high prevalence, incidence, and recurrent rate of VTE. Prior history of VTE, obesity, older age, low platelet account, presence of Philadelphia positive ALL, use of oral contraceptives or hormone replacement therapy, presence of malignancies, and co-morbidities may place leukemia patients at an increased risk for VTE development or recurrence. Interestingly, development of VTE was not associated with a higher risk of death among hospitalized acute leukemia patients.^

Genetic heterogeneity in familial acute myelogenous leukemia: Evidence for a second locus at chromosome 16q21-23.2

The identification of genes responsible for the rare cases of familial leukemia may afford insight into the mechanism underlying the more common sporadic occurrences. Here we test a single family with 11 relevant meioses transmitting autosomal dominant acute myelogenous leukemia (AML) and myelodysplasia for linkage to three potential candidate loci. In a different family with inherited AML, linkage to chromosome 21q22.1-22.2 was recently reported; we exclude linkage to 21q22.1-22.2, demonstrating that familial AML is a heterogeneous disease. After reviewing familial leukemia and observing anticipation in the form of a declining age of onset with each generation, we had proposed 9p21-22 and 16q22 as additional candidate loci. Whereas linkage to 9p21-22 can be excluded, the finding of a maximum two-point LOD score of 2.82 with the microsatellite marker D16S522 at a recombination fraction theta = 0 provides evidence supporting linkage to 16q22. Haplotype analysis reveals a 23.5-cM (17.9-Mb) commonly inherited region among all affected family members extending from D16S451 to D1GS289, In order to extract maximum linkage information with missing individuals, incomplete informativeness with individual markers in this interval, and possible deviance from strict autosomal dominant inheritance, we performed nonparametric linkage analysis (NPL) and found a maximum NPL statistic corresponding to a P-value of .00098, close to the maximum conditional probability of linkage expected for a pedigree with this structure. Mutational analysis in this region specifically excludes expansion of the AT-rich minisatellite repeat FRA16B fragile site and the CAG trinucleotide repeat in the E2F-4 transcription factor. The ''repeat expansion detection'' method, capable of detecting dynamic mutation associated with anticipation, more generally excludes large CAG repeat expansion as a cause of leukemia in this family.

Constitutional hypomorphic telomerase mutations in patients with acute myeloid leukemia

Loss-of-function mutations in telomerase complex genes can cause bone marrow failure, dyskeratosis congenita, and acquired aplastic anemia, both diseases that predispose to acute myeloid leukemia. Loss of telomerase function produces short telomeres, potentially resulting in chromosome recombination, end-to-end fusion, and recognition as damaged DNA. We investigated whether mutations in telomerase genes also occur in acute myeloid leukemia. We screened bone marrow samples from 133 consecutive patients with acute myeloid leukemia and 198 controls for variations in TERT and TERC genes. An additional 89 patients from a second cohort, selected based on cytogenetic status, and 528 controls were further examined for mutations. A third cohort of 372 patients and 384 controls were specifically tested for one TERT gene variant. In the first cohort, 11 patients carried missense TERT gene variants that were not present in controls (P<0.0001); in the second cohort, TERT mutations were associated with trisomy 8 and inversion 16. Mutation germ-line origin was demonstrated in 5 patients from whom other tissues were available. Analysis of all 3 cohorts (n = 594) for the most common gene variant (A1062T) indicated a prevalence 3 times higher in patients than in controls (n = 1,110; P = 0.0009). Introduction of TERT mutants into telomerase-deficient cells resulted in loss of enzymatic activity by haploinsufficiency. Inherited mutations in TERT that reduce telomerase activity are risk factors for acute myeloid leukemia. We propose that short and dysfunctional telomeres limit normal stem cell proliferation and predispose for leukemia by selection of stem cells with defective DNA damage responses that are prone to genome instability.