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.

Domestic radon exposure and risk of childhood cancer: a prospective census-based cohort study

Background: In contrast with established evidence linking high doses of ionizing radiation with childhood cancer, research on low-dose ionizing radiation and childhood cancer has produced inconsistent results. Objective: We investigated the association between domestic radon exposure and childhood cancers, particularly leukemia and central nervous system (CNS) tumors. Methods: We conducted a nationwide census-based cohort study including all children < 16 years of age living in Switzerland on 5 December 2000, the date of the 2000 census. Follow-up lasted until the date of diagnosis, death, emigration, a child’s 16th birthday, or 31 December 2008. Domestic radon levels were estimated for each individual home address using a model developed and validated based on approximately 45,000 measurements taken throughout Switzerland. Data were analyzed with Cox proportional hazard models adjusted for child age, child sex, birth order, parents’ socioeconomic status, environmental gamma radiation, and period effects. Results: In total, 997 childhood cancer cases were included in the study. Compared with children exposed to a radon concentration below the median (< 77.7 Bq/m3), adjusted hazard ratios for children with exposure ≥ the 90th percentile (≥ 139.9 Bq/m3) were 0.93 (95% CI: 0.74, 1.16) for all cancers, 0.95 (95% CI: 0.63, 1.43) for all leukemias, 0.90 (95% CI: 0.56, 1.43) for acute lymphoblastic leukemia, and 1.05 (95% CI: 0.68, 1.61) for CNS tumors. Conclusions: We did not find evidence that domestic radon exposure is associated with childhood cancer, despite relatively high radon levels in Switzerland.

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). ^

Panhandle PCR strategy to amplify MLL genomic breakpoints in treatment-related leukemias

Panhandle PCR amplifies genomic DNA with known 5′ and unknown 3′ sequences from a template with an intrastrand loop schematically shaped like a pan with a handle. We used panhandle PCR to clone MLL genomic breakpoints in two pediatric treatment-related leukemias. The karyotype in a case of treatment-related acute lymphoblastic leukemia showed the t(4;11)(q21;q23). Panhandle PCR amplified the translocation breakpoint at position 2158 in intron 6 in the 5′ MLL breakpoint cluster region (bcr). The karyotype in a case of treatment-related acute myeloid leukemia was normal, but Southern blot analysis showed a single MLL gene rearrangement. Panhandle PCR amplified the breakpoint at position 1493 in MLL intron 6. Screening of somatic cell hybrid and radiation hybrid DNAs by PCR and reverse transcriptase-PCR analysis of the leukemic cells indicated that panhandle PCR identified a fusion of MLL intron 6 with a previously uncharacterized sequence in MLL intron 1, consistent with a partial duplication. In both cases, the breakpoints in the MLL bcr were in Alu repeats, and there were Alu repeats in proximity to the breakpoints in the partner DNAs, suggesting that Alu sequences were relevant to these rearrangements. This study shows that panhandle PCR is an effective method for cloning MLL genomic breakpoints in treatment-related leukemias. Analysis of additional pediatric cases will determine whether breakpoint distribution deviates from the predilection for 3′ distribution in the bcr that has been found in adult cases.

The t(10;11)(p13;q14) in the U937 cell line results in the fusion of the AF10 gene and CALM, encoding a new member of the AP-3 clathrin assembly protein family.

The translocation t(10;11)(p13;q14) is a recurring chromosomal abnormality that has been observed in patients with acute lymphoblastic leukemia as well as acute myeloid leukemia. We have recently reported that the monocytic cell line U937 has a t(10;11)(p13;q14) translocation. Using a combination of positional cloning and candidate gene approach, we cloned the breakpoint and were able to show that AF10 is fused to a novel gene that we named CALM (Clathrin Assembly Lymphoid Myeloid leukemia gene) located at 11q14. AF10, a putative transcription factor, had recently been cloned as one of the fusion partners of MLL. CALM has a very high homology in its N-terminal third to the murine ap-3 gene which is one of the clathrin assembly proteins. The N-terminal region of ap-3 has been shown to bind to clathrin and to have a high-affinity binding site for phosphoinositols. The identification of the CALM/AF10 fusion gene in the widely used U937 cell line will contribute to our understanding of the malignant phenotype of this line.

Selective repression of transcriptional activators by Pbx1 does not require the homeodomain.

PBX1 is a homeobox-containing gene identified as the chromosome 1 participant of the t(1;19) chromosomal translocation of childhood pre-B-cell acute lymphoblastic leukemia. This translocation produces a fusion gene encoding the chimeric oncoprotein E2A-Pbx1, which can induce both acute myeloid and T-lymphoid leukemia in mice. The binding of Pbx1 to DNA is weak; however, both Pbx1 and E2A-Pbx1 exhibit tight binding to specific DNA motifs in conjunction with certain other homeodomain proteins, and E2A-Pbx1 activates transcription through these motifs, whereas Pbx1 does not. In this report, we investigate potential transcriptional functions of Pbx1, using transient expression assays. While no segments of Pbx1 activated transcription, an internal domain of Pbx1 repressed transcription induced by the activation domain of Sp1, but not by the activation domains of VP16 or p53. This Pbx1 domain, which lies upstream of the homeodomain and is highly conserved among Pbx proteins, is thus predicted to bind a specific transcription factor. Surprisingly, the repression activity of Pbx1 did not require homeodomain-dependent DNA binding. Thus, Pbx1 may be able to alter gene transcription by both DNA-binding-dependent and DNA-binding-independent mechanisms.

Membrane-associated CD19-LYN complex is an endogenous p53-independent and Bc1-2-independent regulator of apoptosis in human B-lineage lymphoma cells.

CD19 receptor is expressed at high levels on human B-lineage lymphoid cells and is physically associated with the Src protooncogene family protein-tyrosine kinase Lyn. Recent studies indicate that the membrane-associated CD19-Lyn receptor-enzyme complex plays a pivotal role for survival and clonogenicity of immature B-cell precursors from acute lymphoblastic leukemia patients, but its significance for mature B-lineage lymphoid cells (e.g., B-lineage lymphoma cells) is unknown. CD19-associated Lyn kinase can be selectively targeted and inhibited with B43-Gen, a CD19 receptor-specific immunoconjugate containing the naturally occurring protein-tyrosine kinase inhibitor genistein (Gen). We now present experimental evidence that targeting the membrane-associated CD19-Lyn complex in vitro with B43-Gen triggers rapid apoptotic cell death in highly radiation-resistant p53-Bax- Ramos-BT B-lineage lymphoma cells expressing high levels of Bcl-2 protein without affecting the Bcl-2 expression level. The therapeutic potential of this membrane-directed apoptosis induction strategy was examined in a scid mouse xenograft model of radiation-resistant high-grade human B-lineage lymphoma. Remarkably, in vivo treatment of scid mice challenged with an invariably fatal number of Ramos-BT cells with B43-Gen at a dose level < 1/10 the maximum tolerated dose resulted in 70% long-term event-free survival. Taken together, these results provide unprecedented evidence that the membrane-associated anti-apoptotic CD19-Lyn complex may be at least as important as Bcl-2/Bax ratio for survival of lymphoma cells.

Tumorigenic activity of the BCR-ABL oncogenes is mediated by BCL2.

BCR-ABL is a chimeric oncogene generated by translocation of sequences from the c-abl protein-tyrosine kinase gene on chromosome 9 into the BCR gene on chromosome 22. Alternative chimeric proteins, p210BCR-ABL and p190BCR-ABL, are produced that are characteristic of chronic myelogenous leukemia and acute lymphoblastic leukemia, respectively. Their role in the etiology of human leukemia remains to be defined. Transformed murine hematopoietic cells can be used as a model of BCR-ABL function since these cells can be made growth factor independent and tumorigenic by the action of the BCR-ABL oncogene. We show that the BCR-ABL oncogenes prevent apoptotic death in these cells by inducing a Bcl-2 expression pathway. Furthermore, BCR-ABL-expressing cells revert to factor dependence and nontumorigenicity after Bcl-2 expression is suppressed. These results help to explain the ability of BCR-ABL oncogenes to synergize with c-myc in cell transformation.

Les oncogènes NUP98-PHF23 et NUP98-HOXD13 confèrent un potentiel aberrant d’auto-renouvellement aux progéniteurs thymiques

L’initiation de la leucémogénèse dans la leucémie aigue lymphoblastique (LAL)-T résulte de l’activation aberrante de facteurs de transcription de la lignée lymphocytaire T. Nous démontrons que les gènes de fusion NUP98-PHF23 (NP23) et NUP98-HOXD13 (NHD13) reprogramment les thymocytes normaux en cellules souches pré-leucémiques (CS-préL) possédant un potentiel aberrant d’auto-renouvellement. Basé sur des essais de clonalité performés sur des thymocytes transplantés en série, nous avons découvert que cette population est hiérarchisée similairement aux cellules souches hématopoïétiques normales. Ces CS-préL dévoilent un enrichissement du compartiment de précurseurs thymiques immatures KIT+ où les deux oncogènes, NP23 et NHD13, activent des gènes impliqués dans l’autorenouvellement, incluant Hoxa9, Hoxa10, Lyl1 et Hhex. De plus, l’activité d’autorenouvellement est abrogée par les ARN interférents contre Lyl1 et Hhex, indiquant leur implication fonctionnelle en aval de NP23 et NHD13. Puisque ces gènes sont aussi activés en aval de trois autres oncogènes dans la LAL-T, SCL/TAL1, LMO1 et LMO2, nous concluons que les niveaux d’activation de Lyl1 et Hhex fixent le seuil de reprogrammation des thymocytes normaux en CS-préL. Malgré l'efficacité des traitements de chimiothérapie actuels à diminuer la masse tumorale, les CS-préL sont épargnées, pouvant mener à des rechutes. Nos résultats répondent à ce besoin et proposent de nouvelles avenues permettant de cibler les CS-préL du compartiment de thymocytes immatures dans la LAL-T.

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.

Metabolism Regulates the Fate and Function of T Lymphocytes

Proper balancing of the activities of metabolic pathways to meet the challenge of providing necessary products for biosynthetic and energy demands of the cell is a key requirement for maintaining cell viability and allowing for cell proliferation. Cell metabolism has been found to play a crucial role in numerous cell settings, including in the cells of the immune system, where a successful immune response requires rapid proliferation and successful clearance of dangerous pathogens followed by resolution of the immune response. Additionally, it is now well known that cell metabolism is markedly altered from normal cells in the setting of cancer, where tumor cells rapidly and persistently proliferate. In both settings, alterations to the metabolic profile of the cells play important roles in promoting cell proliferation and survival.

It has long been known that many types of tumor cells and actively proliferating immune cells adopt a metabolic phenotype of aerobic glycolysis, whereby the cell, even under normoxic conditions, imports large amounts of glucose and fluxes it through the glycolytic pathway and produces lactate. However, the metabolic programs utilized by various immune cell subsets have only recently begun to be explored in detail, and the metabolic features and pathways influencing cell metabolism in tumor cells in vivo have not been studied in detail. The work presented here examines the role of metabolism in regulating the function of an important subset of the immune system, the regulatory T cell (Treg) and the role and regulation of metabolism in the context of malignant T cell acute lymphoblastic leukemia (T-ALL). We show that Treg cells, in order to properly function to suppress auto-inflammatory disease, adopt a metabolic program that is characterized by oxidative metabolism and active suppression of anabolic signaling and metabolic pathways. We found that the transcription factor FoxP3, which is highly expressed in Treg cells, drives this phenotype. Perturbing the metabolic phenotype of Treg cells by enforcing increased glycolysis or driving proliferation and anabolic signaling through inflammatory signaling pathways results in a reduction in suppressive function of Tregs.

In our studies focused on the metabolism of T-ALL, we observed that while T-ALL cells use and require aerobic glycolysis, the glycolytic metabolism of T-ALL is restrained compared to that of an antigen activated T cell. The metabolism of T-ALL is instead balanced, with mitochondrial metabolism also being increased. We observed that the pro-anabolic growth mTORC1 signaling pathway was limited in primary T-ALL cells as a result of AMPK pathway activity. AMPK pathway signaling was elevated as a result of oncogene induced metabolic stress. AMPK played a key role in the regulation of T-ALL cell metabolism, as genetic deletion of AMPK in an in vivo murine model of T-ALL resulted in increased glycolysis and anabolic metabolism, yet paradoxically increased cell death and increased mouse survival time. AMPK acts to promote mitochondrial oxidative metabolism in T-ALL through the regulation of Complex I activity, and loss of AMPK reduced mitochondrial oxidative metabolism and resulted in increased metabolic stress. Confirming a role for mitochondrial metabolism in T-ALL, we observed that the direct pharmacological inhibition of Complex I also resulted in a rapid loss of T-ALL cell viability in vitro and in vivo. Taken together, this work establishes an important role for AMPK to both balance the metabolic pathways utilized by T-ALL to allow for cell proliferation and to also promote tumor cell viability by controlling metabolic stress.

Overall, this work demonstrates the importance of the proper coupling of metabolic pathway activity with the function needs of particular types of immune cells. We show that Treg cells, which mainly act to keep immune responses well regulated, adopt a metabolic program where glycolytic metabolism is actively repressed, while oxidative metabolism is promoted. In the setting of malignant T-ALL cells, metabolic activity is surprisingly balanced, with both glycolysis and mitochondrial oxidative metabolism being utilized. In both cases, altering the metabolic balance towards glycolytic metabolism results in negative outcomes for the cell, with decreased Treg functionality and increased metabolic stress in T-ALL. In both cases, this work has generated a new understanding of how metabolism couples to immune cell function, and may allow for selective targeting of immune cell subsets by the specific targeting of metabolic pathways.

Immunologically silent cancer clone transmission from mother to offspring

Rare cases of possible materno-fetal transmission of cancer have been recorded over the past 100 years but evidence for a shared cancer clone has been very limited. We provide genetic evidence for mother to offspring transmission, in utero, of a leukemic cell clone. Maternal and infant cancer clones shared the same unique BCR-ABL1 genomic fusion sequence, indicating a shared, single-cell origin. Microsatellite markers in the infant cancer were all of maternal origin. Additionally, the infant, maternally- derived cancer cells had a major deletion on one copy of chromosome 6p that included deletion of HLA alleles that were not inherited by the infant (i.e., foreign to the infant), suggesting a possible mechanism for immune evasion.