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.

Classification and clinical behavior of blastic plasmacytoid dendritic cell neoplasms according to their maturation-associated immunophenotypic profile

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare subtype of leukemia/lymphoma, whose diagnosis can be difficult to achieve due to its clinical and biological heterogeneity, as well as its overlapping features with other hematologic malignancies. In this study we investigated whether the association between the maturational stage of tumor cells and the clinico-biological and prognostic features of the disease, based on the analysis of 46 BPDCN cases classified into three maturation-associated subgroups on immunophenotypic grounds. Our results show that blasts from cases with an immature plasmacytoid dendritic cell (pDC) phenotype exhibit an uncommon CD56- phenotype, coexisting with CD34+ non-pDC tumor cells, typically in the absence of extramedullary (e.g. skin) disease at presentation. Conversely, patients with a more mature blast cell phenotype more frequently displayed skin/extramedullary involvement and spread into secondary lymphoid tissues. Despite the dismal outcome, acute lymphoblastic leukemia-type therapy (with central nervous system prophylaxis) and/or allogeneic stem cell transplantation appeared to be the only effective therapies. Overall, our findings indicate that the maturational profile of pDC blasts in BPDCN is highly heterogeneous and translates into a wide clinical spectrum -from acute leukemia to mature lymphoma-like behavior-, which may also lead to variable diagnosis and treatment.

Caracterização Hematológica e Imunofenotípica em Pacientes com Leucemia Linfoblástica Aguda

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior

Can Peto's paradox be used as the null hypothesis to identify the role of evolution in natural resistance to cancer? A critical review

BACKGROUND: Carcinogenesis affects not only humans but almost all metazoan species. Understanding the rules driving the occurrence of cancers in the wild is currently expected to provide crucial insights into identifying how some species may have evolved efficient cancer resistance mechanisms. Recently the absence of correlation across species between cancer prevalence and body size (coined as Peto's paradox) has attracted a lot of attention. Indeed, the disparity between this null hypothesis, where every cell is assumed to have an identical probability to undergo malignant transformation, and empirical observations is particularly important to understand, due to the fact that it could facilitate the identification of animal species that are more resistant to carcinogenesis than expected. Moreover it would open up ways to identify the selective pressures that may be involved in cancer resistance. However, Peto's paradox relies on several questionable assumptions, complicating the interpretation of the divergence between expected and observed cancer incidences. DISCUSSIONS: Here we review and challenge the different hypotheses on which this paradox relies on with the aim of identifying how this null hypothesis could be better estimated in order to provide a standard protocol to study the deviation between theoretical/theoretically predicted and observed cancer incidence. We show that due to the disproportion and restricted nature of available data on animal cancers, applying Peto's hypotheses at species level could result in erroneous conclusions, and actually assume the existence of a paradox. Instead of using species level comparisons, we propose an organ level approach to be a more accurate test of Peto's assumptions. SUMMARY: The accuracy of Peto's paradox assumptions are rarely valid and/or quantifiable, suggesting the need to reconsider the use of Peto's paradox as a null hypothesis in identifying the influence of natural selection on cancer resistance mechanisms.

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.

Characterization of myeloid or lymphoid acute leukemia by a chemiluminescence assay: Comparison with immunocytochemistry using an antimyeloperoxidase antibody

A simple and sensitive chemiluminescence assay for the demonstration of the activity of intracellular myeloperoxidase (MPO) is described, which is useful for the distinction between myeloid and lymphoid commitment in blasts from acute leukemia patients. When the cut-off point was settled at 13 mV of chemiluminescence all cases of acute myeloid leukemia (AML) were distinguished from those of acute lymphoid leukemia. In addition, this technique was able to demonstrate MPO activity in AML poorly differentiated (FAB-M0) which usually does not stain for MPO in classical cytochemistry preparations and could be negative also by immunocytochemistry with anti-MPO monoclonal antibody. Therefore the method here described presented a higher sensitivity than the immunocytochemistry procedure with anti-MPO.

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.

Farmacogenomica della Clofarabina nel trattamento delle Leucemie Acute pediatriche: identificazione di nuovi bersagli molecolari e del profilo genomico associato all’efficacia terapeutica del farmaco antitumorale

In quest’ultimi decenni si è assistito ad un notevole miglioramento nella terapia delle Leucemie Acute (LA) pediatriche, nonostante tutto si assiste oggi ad una fase di plateau della curva di sopravvivenza e le leucemie continuano a costituire la principale causa di morte pediatrica per malattia. Ulteriori progressi nel trattamento delle LA potrebbero essere ottenuti mediante studi di farmacogenomica che, identificando le componenti genetiche associate alla risposta individuale ai trattamenti farmacologici, consentono il disegno di terapie personalizzate e tumore-specifiche, ad alta efficacia e bassa tossicità per ciascun paziente. Il lavoro svolto è stato, dunque, finalizzato allo studio della farmacogenomica del farmaco antitumorale Clofarabina (CLO) nel trattamento delle LA pediatriche al fine di identificare marcatori genetici predittivi di risposta delle cellule leucemiche al farmaco, delucidare i meccanismi di resistenza cellulare ed individuare nuovi bersagli verso cui indirizzare terapie più mirate ed efficaci. L’analisi in vitro della sensibilità alla CLO di blasti provenienti da pazienti pediatrici affetti da Leucemia Acuta Linfoblastica (LAL) e Mieloide (LAM) ha consentito l’identificazione di due sottopopolazioni di cellule LAL ad immunofenotipo T a diversa sensibilità alla CLO. Mediante DNA-microarrays, si è identificata la “signature” genetica specificamente associata alla diversa risposta delle cellule LAL-T al farmaco. Successivamente, la caratterizzazione funzionale dei geni differenziali e l’analisi dei pathways hanno consentito l’identificazione specifica di potenziali biomarcatori di risposta terapeutica aprendo nuove prospettive per la comprensione dei meccanismi di resistenza cellulare alla CLO e suggerendo un nuovo bersaglio terapeutico per le forme LAL-T a bassa sensibilità al farmaco. In conclusione, nel lavoro svolto si sono identificati set di geni e pathways di rilievo biologico per la risposta delle cellule LAL-T alla CLO suggerendo marcatori genetici in grado di identificare i soggetti eleggibili per il trattamento o verso cui disegnare terapie innovative. Il lavoro è paradigma per l’applicazione della farmacogenomica in altre neoplasie.

Arabinosylguanine: Metabolism, action, and lineage-specific cytotoxicity in human leukemia cells

9-β-D-arabinosylguanine (ara-G), an analogue of deoxyguanosine, has demonstrated T-lymphoblast selective anti-leukemia activity both in vitro and in vivo in cell lines and primary cells and in phase I investigations. The present work was initiated to identify factors that result in this selectivity. ^ The cytotoxicity of ara-G is manifest only after its phosphorylation. Experiments using cell lines transfected to overexpress specific nucleoside kinases demonstrated that the phosphorylation of ara-G to its monophosphate is by both cytoplasmic deoxycytidine kinase and mitochondria) deoxyguanosine kinase. Ara-G monophosphate is converted to its 5′-triphosphate (ara-GTP) in cells by these kinases and then incorporated into DNA. Mechanistic studies demonstrated that incorporation of ara-GTP into DNA was a necessary event for the induction of cell death. ^ Pharmacokinetic and pharmacodynamic studies utilizing three human acute leukemia cell lines, CEM (T-lymphoblastic), Raji (B-lymphoblastic), and ML-1 (myeloid) were performed. CEM cells were most sensitive to ara-G-induced inhibition of colony formation, accumulated ara-GTP at a faster rate and to a greater degree than either Raji or ML-1, but incorporated the lowest number of ara-G molecules into DNA. The position of incorporation was internal and similar in all cell lines. The terminal elimination phase of ara-GTP was >24 h and similar in these cells. Comparisons between inhibition of colony formation and ara-GTP incorporation into DNA demonstrated that while within a cell line there was correlation among these parameters, between cell lines there was no relationship between number of incorporated ara-G molecules and ara-G(TP)-mediated toxicity suggesting that there were additional factors. ^ The expression of membrane bound Fas and Fast was unchanged in all cell lines. In contrast, there was a 2-fold increase in soluble Fast, which was found exclusively in CEM cells. Ara-G-mediated apoptosis in CEM occurred from all phases of the cell cycle and was abrogated partially by Fas antagonist antibodies. These data suggest that Fas-mediated cell death due to the liberation of sFasL may be responsible for the hypersensitivity to ara-G manifested by immature T-cells such as CEM. The role of Fas in ara-G induced death of acute T-lymphoblastic leukemia cells during therapy needs to be tested. ^

Polymorphisms in the methylenetetrahydrofolate reductase gene are associated with susceptibility to acute leukemia in adults

Reduction of 5,10-methylenetetrahydrofolate (methyleneTHF), a donor for methylating dUMP to dTMP in DNA synthesis, to 5-methyltetrahydrofolate (methylTHF), the primary methyl donor for methionine synthesis, is catalyzed by 5,10-methylenetetrahydrofolate reductase (MTHFR). A common 677 C → T polymorphism in the MTHFR gene results in thermolability and reduced MTHFR activity that decreases the pool of methylTHF and increases the pool of methyleneTHF. Recently, another polymorphism in MTHFR (1298 A → C) has been identified that also results in diminished enzyme activity. We tested whether carriers of these variant alleles are protected from adult acute leukemia. We analyzed DNA from a case–control study in the United Kingdom of 308 adult acute leukemia patients and 491 age- and sex-matched controls. MTHFR variant alleles were determined by a PCR-restriction fragment length polymorphism assay. The MTHFR 677TT genotype was lower among 71 acute lymphocytic leukemia (ALL) cases compared with 114 controls, conferring a 4.3-fold decrease in risk of ALL [odds ratio (OR = 0.23; 95% CI = 0.06–0.81]. We observed a 3-fold reduction in risk of ALL in individuals with the MTHFR 1298AC polymorphism (OR = 0.33; 95% CI = 0.15–0.73) and a 14-fold decreased risk of ALL in those with the MTHFR 1298CC variant allele (OR = 0.07; 95% CI = 0.00–1.77). In acute myeloid leukemia, no significant difference in MTHFR 677 and 1298 genotype frequencies was observed between 237 cases and 377 controls. Individuals with the MTHFR 677TT, 1298AC, and 1298CC genotypes have a decreased risk of adult ALL, but not acute myeloid leukemia, which suggests that folate inadequacy may play a key role in the development of ALL.