Integration of genomic and gene expression data of childhood ALL without known aberrations identifies subgroups with specific genetic hallmarks

Angiotensin-converting Enzyme Insertion/Deletion Polymorphism (rs4646994) and Susceptibility to Acute Lymphoblastic Leukemia: A Case–control Study

BACKGROUND: Angiotensin-converting enzyme (ACE) stimulates the proliferation of bone marrow hematopoietic progenitors and thought to be involved in pathological neoplastic hematopoiesis and leukemogenesis. AIM: This study aimed to investigate the association between ACE gene I/D polymorphism and the risk of acute lymphoblastic leukemia (ALL). MATERIALS AND METHODS: A total of 96 subjects were recruited for the study, 48 patients with ALL, and 48 apparently healthy volunteers as a control group. Genomic DNA was extracted from peripheral leukocytes and ACE I/D polymorphism was analysed using allele-specific polymerase chain reaction. RESULTS: In both study groups, the ACE D/D polymorphic genotype was the most frequent (52.1% and 54.2%, respectively), followed by the ID genotype (47.9% and 45.8% respectively), while the II genotype was completely absent in both study groups. The distribution of the polymorphic genotypes among the study groups was not significantly different (p = 0. 0.398). The frequency of the D allele was 0.76 in the patients and 0.77 in the control group, while the frequency of I allele was 0.24 in the patients and 0.23 in the control group. No deviation from Hardy–Weinberg equilibrium was observed (χ2 = 4.24, df = 1, p = 0.12). CONCLUSION: ACE I/D polymorphism is not associated with susceptibility to ALL among the Sudanese population.

NT5E gene and CD38 protein as potential prognostic biomarkers for childhood B-acute lymphoblastic leukemia

The risk stratification of B-acute lymphoblastic leukemia (B-ALL) is based on clinical and biological factors. However, B-ALL has significant biological and clinical heterogeneity and 50% of B-ALL patients do not have defined prognostic markers. In this sense, the identification of new prognostic biomarkers is necessary. Considering different cohorts of childhood B-ALL patients, gene (DPP4/CD38/ENTPD1/NT5E) and protein (CD38/CD39/CD73) expressions of ectonucleotidases were analyzed in silico and ex vivo and the association with prognosis was established. In univariate analyses, expression of NT5E was significantly associated with worse progression-free survival (PFS) in bone marrow (BM) samples. In multivariate analyses, Kaplan-Meier analysis, and log-rank test, higher NT5E expression predicted unfavorable PFS in BM samples. Considering minimal residual disease (MRD), higher levels of cellularity were associated with the high NT5E expression at day 8 of induction therapy. In addition, we observed that white blood cells (WBC) of childhood B-ALL patients had more CD38 compared to the same cell population of healthy donors (HD). In fact, MRD > 0.1% patients had higher CD38 protein expression on WBC in comparison to HD. Noteworthy, we observed higher CD38 expression on WBC than blasts in MRD > 0.1% patients. We suggest that NT5E gene and CD38 protein expression, of the ectonucleotidases family, could provide interesting prognostic biomarkers for childhood B-ALL.

Identifying novel genetic alterations in pediatric acute lymphoblastic leukemia based on copy number analysis

Copy number variations (CNVs) analysis may reveal molecular biomarkers and provide information on the pathogenesis of acute lymphoblastic leukemia (ALL). We investigated the gene copy number in childhood ALL by microarray and select three new recurrent CNVs to evaluate by real-time PCR assay: DMBT1, KIAA0125 and PRDM16 were selected due to high frequency of CNVs in ALL samples and based on their potential biological functions in carcinogenesis described in the literature. DBMT1 deletion was associated with patients with chromosomal translocations and is a potential tumor suppressor; KIAA0125 and PRDM16 may act as an oncogene despite having a paradoxical behavior in carcinogenesis. This study reinforces that microarrays/aCGH is it is a powerful tool for detection of genomic aberrations, which may be used in the risk stratification.

Feedforward regulation of Myc coordinates lineage-specific with housekeeping gene expression during B cell progenitor cell differentiation

The differentiation of self-renewing progenitor cells requires not only the regulation of lineage- and developmental stage–specific genes but also the coordinated adaptation of housekeeping functions from a metabolically active, proliferative state toward quiescence. How metabolic and cell-cycle states are coordinated with the regulation of cell type–specific genes is an important question, because dissociation between differentiation, cell cycle, and metabolic states is a hallmark of cancer. Here, we use a model system to systematically identify key transcriptional regulators of Ikaros-dependent B cell–progenitor differentiation. We find that the coordinated regulation of housekeeping functions and tissue-specific gene expression requires a feedforward circuit whereby Ikaros down-regulates the expression of Myc. Our findings show how coordination between differentiation and housekeeping states can be achieved by interconnected regulators. Similar principles likely coordinate differentiation and housekeeping functions during progenitor cell differentiation in other cell lineages.

The human body is made from billions of cells comprizing many specialized cell types. All of these cells ultimately come from a single fertilized oocyte in a process that has two key features: proliferation, which expands cell numbers, and differentiation, which diversifies cell types. Here, we have examined the transition from proliferation to differentiation using B lymphocytes as an example. We find that the transition from proliferation to differentiation involves changes in the expression of genes, which can be categorized into cell-type–specific genes and broadly expressed “housekeeping” genes. The expression of many housekeeping genes is controlled by the gene regulatory factor Myc, whereas the expression of many B lymphocyte–specific genes is controlled by the Ikaros family of gene regulatory proteins. Myc is repressed by Ikaros, which means that changes in housekeeping and tissue-specific gene expression are coordinated during the transition from proliferation to differentiation.

Regulatory Network and Prognostic Effect Investigation of PIP4K2A in Leukemia and Solid Cancers

Germline variants of PIP4K2A impact susceptibility of acute lymphoblastic leukemia (ALL) through inducing its overexpression. Although limited reports suggested the oncogenic role of PIP4K2A in cancers, regulatory network and prognostic effect of this gene remains poorly understood in tumorigenesis and leukemogenesis. In this study, we conducted genome-wide gene expression association analyses in pediatric B-ALL cohorts to discover expression associated genes and pathways, which is followed by the bioinformatics analyses to investigate the prognostic role of PIP4K2A and its related genes in multiple cancer types. 214 candidates were identified to be significantly associated with PIP4K2A expression in ALL patients, with known cancer-related genes rankings the top (e.g., RAC2, RBL2, and TFDP1). These candidates do not only tend to be clustered in the same types of leukemia, but can also separate the patients into novel molecular subtypes. PIP4K2A is noticed to be frequently overexpressed in multiple other types of leukemia and solid cancers from cancer cohorts including TCGA, and associated with its candidates in subtype-specific and cancer-specific manners. Interestingly, the association status varied in tumors compared to their matched normal tissues. Moreover, PIP4K2A and its related candidates exhibit stage-independent prognostic effects in multiple cancers, mostly with its lower expression significantly associated with longer overall survival (p < 0.05). Our findings reveal the transcriptional regulatory network of PIP4K2A in leukemia, and suggest its potentially important role on molecular subtypes of multiple cancers and subsequent treatment outcomes.

Regulatory network of GATA3 in pediatric acute lymphoblastic leukemia

GATA3 polymorphisms were reported to be significantly associated with susceptibility of pediatric B-lineage acute lymphoblastic leukemia (ALL), by impacting on GATA3 expression. We noticed that ALL-related GATA3 polymorphism located around in the tissue-specific enhancer, and significantly associated with GATA3 expression. Although the regulatory network of GATA3 has been well reported in T cells, the functional status of GATA3 is poorly understood in B-ALL. We thus conducted genome-wide gene expression association analyses to reveal expression associated genes and pathways in nine independent B-ALL patient cohorts. In B-ALL patients, 173 candidates were identified to be significantly associated with GATA3 expression, including some reported GATA3-related genes (e.g., ITM2A) and well-known tumor-related genes (e.g., STAT4). Some of the candidates exhibit tissue-specific and subtype-specific association with GATA3. Through overexpression and down-regulation of GATA3 in leukemia cell lines, several reported and novel GATA3 regulated genes were validated. Moreover, association of GATA3 expression and its targets can be impacted by SNPs (e.g., rs4894953), which locate in the potential GATA3 binding motif. Our findings suggest that GATA3 may be involved in multiple tumor-related pathways (e.g., STAT/JAK pathway) in B-ALL to impact leukemogenesis through epigenetic regulation.

Mutations in TP53 and JAK2 are independent prognostic biomarkers in B-cell precursor acute lymphoblastic leukaemia

BACKGROUND

In B-cell precursor acute lymphoblastic leukaemia (B-ALL), the identification of additional genetic alterations associated with poor prognosis is still of importance. We determined the frequency and prognostic impact of somatic mutations in children and adult cases with B-ALL treated with Spanish PETHEMA and SEHOP protocols.

METHODS

Mutational status of hotspot regions of TP53, JAK2, PAX5, LEF1, CRLF2 and IL7R genes was determined by next-generation deep sequencing in 340 B-ALL patients (211 children and 129 adults). The associations between mutation status and clinicopathological features at the time of diagnosis, treatment outcome and survival were assessed. Univariate and multivariate survival analyses were performed to identify independent prognostic factors associated with overall survival (OS), event-free survival (EFS) and relapse rate (RR).

RESULTS

A mutation rate of 12.4% was identified. The frequency of adult mutations was higher (20.2% vs 7.6%, P=0.001). TP53 was the most frequently mutated gene (4.1%), followed by JAK2 (3.8%), CRLF2 (2.9%), PAX5 (2.4%), LEF1 (0.6%) and IL7R (0.3%). All mutations were observed in B-ALL without ETV6-RUNX1 (P=0.047) or BCR-ABL1 fusions (P<0.0001). In children, TP53mut was associated with lower OS (5-year OS: 50% vs 86%, P=0.002) and EFS rates (5-year EFS: 50% vs 78.3%, P=0.009) and higher RR (5-year RR: 33.3% vs 18.6% P=0.037), and was independently associated with higher RR (hazard ratio (HR)=4.5; P=0.04). In adults, TP53mut was associated with a lower OS (5-year OS: 0% vs 43.3%, P=0.019) and a higher RR (5-year RR: 100% vs 61.4%, P=0.029), whereas JAK2mut was associated with a lower EFS (5-year EFS: 0% vs 30.6%, P=0.035) and a higher RR (5-year RR: 100% vs 60.4%, P=0.002). TP53mut was an independent risk factor for shorter OS (HR=2.3; P=0.035) and, together with JAK2mut, also were independent markers of poor prognosis for RR (TP53mut: HR=5.9; P=0.027 and JAK2mut: HR=5.6; P=0.036).

CONCLUSIONS

TP53mut and JAK2mut are potential biomarkers associated with poor prognosis in B-ALL patients.

Aberrantly expressed LGR4 empowers Wnt signaling in multiple myeloma by hijacking osteoblast-derived R-spondins

The unrestrained growth of tumor cells is generally attributed to mutations in essential growth control genes, but tumor cells are also affected by, or even addicted to, signals from the microenvironment. As therapeutic targets, these extrinsic signals may be equally significant as mutated oncogenes. In multiple myeloma (MM), a plasma cell malignancy, most tumors display hallmarks of active Wnt signaling but lack activating Wnt-pathway mutations, suggesting activation by autocrine Wnt ligands and/or paracrine Wnts emanating from the bone marrow (BM) niche. Here, we report a pivotal role for the R-spondin/leucine-rich repeat-containing G protein-coupled receptor 4 (LGR4) axis in driving aberrant Wnt/β-catenin signaling in MM. We show that LGR4 is expressed by MM plasma cells, but not by normal plasma cells or B cells. This aberrant LGR4 expression is driven by IL-6/STAT3 signaling and allows MM cells to hijack R-spondins produced by (pre)osteoblasts in the BM niche, resulting in Wnt (co)receptor stabilization and a dramatically increased sensitivity to auto- and paracrine Wnts. Our study identifies aberrant R-spondin/LGR4 signaling with consequent deregulation of Wnt (co)receptor turnover as a driver of oncogenic Wnt/β-catenin signaling in MM cells. These results advocate targeting of the LGR4/R-spondin interaction as a therapeutic strategy in MM.

Genome-Wide DNA Copy Number Analysis of Acute Lymphoblastic Leukemia Identifies New Genetic Markers Associated with Clinical Outcome

Identifying additional genetic alterations associated with poor prognosis in acute lymphoblastic leukemia (ALL) is still a challenge. Aims: To characterize the presence of additional DNA copy number alterations (CNAs) in children and adults with ALL by whole-genome oligonucleotide array (aCGH) analysis, and to identify their associations with clinical features and outcome. Array-CGH was carried out in 265 newly diagnosed ALLs (142 children and 123 adults). The NimbleGen CGH 12x135K array (Roche) was used to analyze genetic gains and losses. CNAs were analyzed with GISTIC and aCGHweb software. Clinical and biological variables were analyzed. Three of the patients showed chromothripsis (cth6, cth14q and cth15q). CNAs were associated with age, phenotype, genetic subtype and overall survival (OS). In the whole cohort of children, the losses on 14q32.33 (p = 0.019) and 15q13.2 (p = 0.04) were related to shorter OS. In the group of children without good- or poor-risk cytogenetics, the gain on 1p36.11 was a prognostic marker independently associated with shorter OS. In adults, the gains on 19q13.2 (p = 0.001) and Xp21.1 (p = 0.029), and the loss of 17p (p = 0.014) were independent markers of poor prognosis with respect to OS. In summary, CNAs are frequent in ALL and are associated with clinical parameters and survival. Genome-wide DNA copy number analysis allows the identification of genetic markers that predict clinical outcome, suggesting that detection of these genetic lesions will be useful in the management of patients newly diagnosed with ALL.

SNP Array in Hematopoietic Neoplasms: A Review

Cytogenetic analysis is essential for the diagnosis and prognosis of hematopoietic neoplasms in current clinical practice. Many hematopoietic malignancies are characterized by structural chromosomal abnormalities such as specific translocations, inversions, deletions and/or numerical abnormalities that can be identified by karyotype analysis or fluorescence in situ hybridization (FISH) studies. Single nucleotide polymorphism (SNP) arrays offer high-resolution identification of copy number variants (CNVs) and acquired copy-neutral loss of heterozygosity (LOH)/uniparental disomy (UPD) that are usually not identifiable by conventional cytogenetic analysis and FISH studies. As a result, SNP arrays have been increasingly applied to hematopoietic neoplasms to search for clinically-significant genetic abnormalities. A large numbers of CNVs and UPDs have been identified in a variety of hematopoietic neoplasms. CNVs detected by SNP array in some hematopoietic neoplasms are of prognostic significance. A few specific genes in the affected regions have been implicated in the pathogenesis and may be the targets for specific therapeutic agents in the future. In this review, we summarize the current findings of application of SNP arrays in a variety of hematopoietic malignancies with an emphasis on the clinically significant genetic variants.

Identification of renin progenitors in the mouse bone marrow that give rise to B-cell leukaemia

The cell of origin and triggering events for leukaemia are mostly unknown. Here we show that the bone marrow contains a progenitor that expresses renin throughout development and possesses a B-lymphocyte pedigree. This cell requires RBP-J to differentiate. Deletion of RBP-J in these renin-expressing progenitors enriches the precursor B-cell gene programme and constrains lymphocyte differentiation, facilitated by H3K4me3 activating marks in genes that control the pre-B stage. Mutant cells undergo neoplastic transformation, and mice develop a highly penetrant B-cell leukaemia with multi-organ infiltration and early death. These renin-expressing cells appear uniquely vulnerable as other conditional models of RBP-J deletion do not result in leukaemia. The discovery of these unique renin progenitors in the bone marrow and the model of leukaemia described herein may enhance our understanding of normal and neoplastic haematopoiesis.

Prognostification of ALL by Cytogenetics

Cytogenetic abnormalities in chromosomal number and structure are common in pediatric ALL and some have prognostic significance. One interesting association between cytogenetic status and treatment response involves the metabolism of methotrexate. Hyperdiploid lymphoblasts accumulate increased amounts of MTX and MTX polyglutamates, and they have higher basal apoptotic rates compared with leukemic cells with lower ploidy and normal cells. These characteristics may contribute to the better outcomes observed for patients with hyperdiploid lymphoblasts. A number of recurrent chromosomal abnormalities have been shown to have prognostic significance, especially in B-precursor ALL. Some chromosomal abnormalities are associated with more favorable outcomes, such as high hyperdiploidy (51-65 chromosomes) and the ETV6-RUNX1 fusion. Others are associated with a poorer prognosis, including the Philadelphia chromosome [t(9;22)], rearrangements of the MLL gene (chromosome 11q23), and intrachromosomal amplification of the AML1 gene (iAMP21).

An advanced fragment analysis-based individualized subtype classification of pediatric acute lymphoblastic leukemia

Pediatric acute lymphoblastic leukemia (ALL) is the most common neoplasm and one of the primary causes of death in children. Its treatment is highly dependent on the correct classification of subtype. Previously, we developed a microarray-based subtype classifier based on the relative expression levels of 62 marker genes, which can predict 7 different ALL subtypes with an accuracy as high as 97% in completely independent samples. Because the classifier is based on gene expression rank values rather than actual values, the classifier enables an individualized diagnosis, without the need to reference the background distribution of the marker genes in a large number of other samples, and also enables cross platform application. Here, we demonstrate that the classifier can be extended from a microarray-based technology to a multiplex qPCR-based technology using the same set of marker genes as the advanced fragment analysis (AFA). Compared to microarray assays, the new assay system makes the convenient, low cost and individualized subtype diagnosis of pediatric ALL a reality and is clinically applicable, particularly in developing countries.

Consequences of Aneuploidy in Cancer: Transcriptome and Beyond

Cancer cells differ from normal healthy cells in multiple aspects ranging from altered cellular signaling through metabolic changes to aberrant chromosome content, so called aneuploidy. The large-scale changes in copy numbers of chromosomes or large chromosomal regions due to aneuploidy alter significantly the gene expression, as several hundreds of genes are gained or lost. Comparison of quantitative genome, transcriptome and proteome data enables dissection of the molecular causes that underlie the gene expression changes observed in cancer cells and provides a new perspective on the molecular consequences of aneuploidy. Here, we will map to what degree aneuploidy affects the expression of genes located on the affected chromosomes. We will also address the effects of aneuploidy on global gene expression in cancer cells as well as whether and how it may contribute to the physiology of cancer cells.

Integrated DNA Copy Number and Gene Expression Regulatory Network Analysis of Non-small Cell Lung Cancer Metastasis

Integrative analysis of multi-level molecular profiles can distinguish interactions that cannot be revealed based on one kind of data in the analysis of cancer susceptibility and metastasis. DNA copy number variations (CNVs) are common in cancer cells, and their role in cell behaviors and relationship to gene expression (GE) is poorly understood. An integrative analysis of CNV and genome-wide mRNA expression can discover copy number alterations and their possible regulatory effects on GE. This study presents a novel framework to identify important genes and construct potential regulatory networks based on these genes. Using this approach, DNA copy number aberrations and their effects on GE in lung cancer progression were revealed. Specifically, this approach contains the following steps: (1) select a pool of candidate driver genes, which have significant CNV in lung cancer patient tumors or have a significant association with the clinical outcome at the transcriptional level; (2) rank important driver genes in lung cancer patients with good prognosis and poor prognosis, respectively, and use top-ranked driver genes to construct regulatory networks with the COpy Number and EXpression In Cancer (CONEXIC) method; (3) identify experimentally confirmed molecular interactions in the constructed regulatory networks using Ingenuity Pathway Analysis (IPA); and (4) visualize the refined regulatory networks with the software package Genatomy. The constructed CNV/mRNA regulatory networks provide important insights into potential CNV-regulated transcriptional mechanisms in lung cancer metastasis.

CGBayesNets: conditional Gaussian Bayesian network learning and inference with mixed discrete and continuous data

Bayesian Networks (BN) have been a popular predictive modeling formalism in bioinformatics, but their application in modern genomics has been slowed by an inability to cleanly handle domains with mixed discrete and continuous variables. Existing free BN software packages either discretize continuous variables, which can lead to information loss, or do not include inference routines, which makes prediction with the BN impossible. We present CGBayesNets, a BN package focused around prediction of a clinical phenotype from mixed discrete and continuous variables, which fills these gaps. CGBayesNets implements Bayesian likelihood and inference algorithms for the conditional Gaussian Bayesian network (CGBNs) formalism, one appropriate for predicting an outcome of interest from, e.g., multimodal genomic data. We provide four different network learning algorithms, each making a different tradeoff between computational cost and network likelihood. CGBayesNets provides a full suite of functions for model exploration and verification, including cross validation, bootstrapping, and AUC manipulation. We highlight several results obtained previously with CGBayesNets, including predictive models of wood properties from tree genomics, leukemia subtype classification from mixed genomic data, and robust prediction of intensive care unit mortality outcomes from metabolomic profiles. We also provide detailed example analysis on public metabolomic and gene expression datasets. CGBayesNets is implemented in MATLAB and available as MATLAB source code, under an Open Source license and anonymous download at http://www.cgbayesnets.com.

Frequent and sex-biased deletion of SLX4IP by illegitimate V(D)J-mediated recombination in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) accounts for ∼25% of pediatric malignancies. Of interest, the incidence of ALL is observed ∼20% higher in males relative to females. The mechanism behind the phenomenon of sex-specific differences is presently not understood. Employing genome-wide genetic aberration screening in 19 ALL samples, one of the most recurrent lesions identified was monoallelic deletion of the 5' region of SLX4IP. We characterized this deletion by conventional molecular genetic techniques and analyzed its interrelationships with biological and clinical characteristics using specimens and data from 993 pediatric patients enrolled into trial AIEOP-BFM ALL 2000. Deletion of SLX4IP was detected in ∼30% of patients. Breakpoints within SLX4IP were defined to recurrent positions and revealed junctions with typical characteristics of illegitimate V(D)J-mediated recombination. In initial and validation analyses, SLX4IP deletions were significantly associated with male gender and ETV6/RUNX1-rearranged ALL (both overall P < 0.0001). For mechanistic validation, a second recurrent deletion affecting TAL1 and caused by the same molecular mechanism was analyzed in 1149 T-cell ALL patients. Validating a differential role by sex of illegitimate V(D)J-mediated recombination at the TAL1 locus, 128 out of 1149 T-cell ALL samples bore a deletion and males were significantly more often affected (P = 0.002). The repeatedly detected association of SLX4IP deletion with male sex and the extension of the sex bias to deletion of the TAL1 locus suggest that differential illegitimate V(D)J-mediated recombination events at specific loci may contribute to the consistent observation of higher incidence rates of childhood ALL in boys compared with girls.

IKZF1 status as a prognostic feature in BCR-ABL1-positive childhood ALL

Childhood BCR-ABL1-positive B-cell precursor acute lymphoblastic leukemia (BCP-ALL) has an unfavorable outcome and shows high frequency of IKZF1 deletions. The prognostic value of IKZF1 deletions was evaluated in 2 cohorts of BCR-ABL1-positive BCP-ALL patients, before tyrosine kinase inhibitors (pre-TKI) and after introduction of imatinib (in the European Study for Philadelphia-Acute Lymphoblastic Leukemia [EsPhALL]). In 126/191 (66%) cases an IKZF1 deletion was detected. In the pre-TKI cohort, IKZF1-deleted patients had an unfavorable outcome compared with wild-type patients (4-year disease-free survival [DFS] of 30.0 ± 6.8% vs 57.5 ± 9.4%; P = .01). In the EsPhALL cohort, the IKZF1 deletions were associated with an unfavorable prognosis in patients stratified in the good-risk arm based on early clinical response (4-year DFS of 51.9 ± 8.8% for IKZF1-deleted vs 78.6 ± 13.9% for IKZF1 wild-type; P = .03), even when treated with imatinib (4-year DFS of 55.5 ± 9.5% for IKZF1-deleted vs 75.0 ± 21.7% for IKZF1 wild-type; P = .05). In conclusion, the highly unfavorable outcome for childhood BCR-ABL1-positive BCP-ALL with IKZF1 deletions, irrespective of imatinib exposure, underscores the need for alternative therapies. In contrast, good-risk patients with IKZF1 wild-type responded remarkably well to imatinib-containing regimens, providing a rationale to potentially avoid hematopoietic stem-cell transplantation in this subset of patients.

An international study of intrachromosomal amplification of chromosome 21 (iAMP21): cytogenetic characterization and outcome

Intrachromosomal amplification of chromosome 21 (iAMP21) defines a distinct cytogenetic subgroup of childhood B-cell precursor acute lymphoblastic leukaemia (BCP-ALL). To date, fluorescence in situ hybridisation (FISH), with probes specific for the RUNX1 gene, provides the only reliable detection method (five or more RUNX1 signals per cell). Patients with iAMP21 are older (median age 9 years) with a low white cell count. Previously, we demonstrated a high relapse risk when these patients were treated as standard risk. Recent studies have shown improved outcome on intensive therapy. In view of these treatment implications, accurate identification is essential. Here we have studied the cytogenetics and outcome of 530 iAMP21 patients that highlighted the association of specific secondary chromosomal and genetic changes with iAMP21 to assist in diagnosis, including the gain of chromosome X, loss or deletion of chromosome 7, ETV6 and RB1 deletions. These iAMP21 patients when treated as high risk showed the same improved outcome as those in trial-based studies regardless of the backbone chemotherapy regimen given. This study reinforces the importance of intensified treatment to reduce the risk of relapse in iAMP21 patients. This now well-defined patient subgroup should be recognised by World Health Organisation (WHO) as a distinct entity of BCP-ALL.

A phenolic ester from Aglaia loheri leaves reveals cytotoxicity towards sensitive and multidrug-resistant cancer cells

BACKGROUND

Bioactivity-guided fractionation of extracts of Aglaia loheri Blanco (Meliaceae) yielded a cytotoxic isolate, termed Maldi 531.2[M + H]+. This phenolic ester was further investigated for its in vitro cytotoxicity toward human CCRF-CEM leukemia cells and their multi-drug resistant (MDR) subline, CEM/ADR5000. The intrinsic mitochondrial membrane potential (ΔΨm) and induction of apoptosis by this isolate were evaluated.

METHODS

Chromatography techniques, mass spectrometry and proton NMR were employed to isolate Maldi 531.2[M + H]+. XTT cell proliferation and viability assay was used for cytotoxic test, and JC-1[5',5',6,6',-tetrachloro-1,1',3,3'-tetraethylbenzimidazoyl carbocyanine iodide was used to assess ΔΨm and initiation of apoptosis; Annexin V/FITC-PI staining was employed to analyse apoptosis.

RESULTS

Maldi 531.2[M + H]+ was cytotoxic towards both CCRF-CEM and CEM/ADR5000 cells with IC50 values of 0.02 and 0.03 μM, respectively. The mitochondrial membrane potential (ΔΨm) of MDR cells was significantly reduced in a dose-dependent manner leading to apoptosis as detected by flow cytometric Annexin V-FITC/ PI staining.

CONCLUSION

Maldi 531.2[M + H]+ may be a potential anti-cancer drug candidate whose mode of action include reduction of the mitochondrial membrane potential and induction of apoptosis.

Micro-Analyzer: automatic preprocessing of Affymetrix microarray data

A current trend in genomics is the investigation of the cell mechanism using different technologies, in order to explain the relationship among genes, molecular processes and diseases. For instance, the combined use of gene-expression arrays and genomic arrays has been demonstrated as an effective instrument in clinical practice. Consequently, in a single experiment different kind of microarrays may be used, resulting in the production of different types of binary data (images and textual raw data). The analysis of microarray data requires an initial preprocessing phase, that makes raw data suitable for use on existing analysis platforms, such as the TIGR M4 (TM4) Suite. An additional challenge to be faced by emerging data analysis platforms is the ability to treat in a combined way those different microarray formats coupled with clinical data. In fact, resulting integrated data may include both numerical and symbolic data (e.g. gene expression and SNPs regarding molecular data), as well as temporal data (e.g. the response to a drug, time to progression and survival rate), regarding clinical data. Raw data preprocessing is a crucial step in analysis but is often performed in a manual and error prone way using different software tools. Thus novel, platform independent, and possibly open source tools enabling the semi-automatic preprocessing and annotation of different microarray data are needed. The paper presents Micro-Analyzer (Microarray Analyzer), a cross-platform tool for the automatic normalization, summarization and annotation of Affymetrix gene expression and SNP binary data. It represents the evolution of the μ-CS tool, extending the preprocessing to SNP arrays that were not allowed in μ-CS. The Micro-Analyzer is provided as a Java standalone tool and enables users to read, preprocess and analyse binary microarray data (gene expression and SNPs) by invoking TM4 platform. It avoids: (i) the manual invocation of external tools (e.g. the Affymetrix Power Tools), (ii) the manual loading of preprocessing libraries, and (iii) the management of intermediate files, such as results and metadata. Micro-Analyzer users can directly manage Affymetrix binary data without worrying about locating and invoking the proper preprocessing tools and chip-specific libraries. Moreover, users of the Micro-Analyzer tool can load the preprocessed data directly into the well-known TM4 platform, extending in such a way also the TM4 capabilities. Consequently, Micro Analyzer offers the following advantages: (i) it reduces possible errors in the preprocessing and further analysis phases, e.g. due to the incorrect choice of parameters or due to the use of old libraries, (ii) it enables the combined and centralized pre-processing of different arrays, (iii) it may enhance the quality of further analysis by storing the workflow, i.e. information about the preprocessing steps, and (iv) finally Micro-Analzyer is freely available as a standalone application at the project web site http://sourceforge.net/projects/microanalyzer/.

Linking genomic lesions with minimal residual disease improves prognostic stratification in children with T-cell acute lymphoblastic leukaemia

Multiple lesions in genes that are involved in cell cycle control, proliferation, survival and differentiation underlie T-cell acute lymphoblastic leukaemia (T-ALL). We translated these biological insights into clinical practice to improve diagnostic work-ups and patient management. Combined interphase fluorescence in situ hybridization (CI-FISH), single nucleotide polymorphism (SNP), and gene expression profiles (GEP) were applied in 51 children with T-ALL who were stratified according to minimal residual disease (MRD) risk categories (AIEOP-BFM ALL2000). CI-FISH identified type A abnormalities in 90% of patients. Distribution of each was in line with the estimated incidence in childhood T-ALL: 37.5% TAL/LMO, 22.5% HOXA, 20% TLX3, 7.5% TLX1, and 2.5% NKX2-1. GEP predictions concurred. SNP detected type B abnormalities in all cases, thus linking type A and B lesions. This approach provided an accurate, comprehensive genomic diagnosis and a complementary GEP-based classification of T-ALL in children. Dissecting primary and secondary lesions within MRD categories could improve prognostic criteria for the majority of patients and be a step towards personalized diagnosis.

Genes commonly deleted in childhood B-cell precursor acute lymphoblastic leukemia: association with cytogenetics and clinical features

In childhood B-cell precursor acute lymphoblastic leukemia, cytogenetics is important in diagnosis and as an indicator of response to therapy, thus playing a key role in risk stratification of patients for treatment. Little is known of the relationship between different cytogenetic subtypes in B-cell precursor acute lymphoblastic leukemia and the recently reported copy number abnormalities affecting significant leukemia associated genes. In a consecutive series of 1427 childhood B-cell precursor acute lymphoblastic leukemia patients, we have determined the incidence and type of copy number abnormalities using multiplex ligation-dependent probe amplification. We have shown strong links between certain deletions and cytogenetic subtypes, including the novel association between RB1 deletions and intrachromosomal amplification of chromosome 21. In this study, we characterized the different copy number abnormalities and show heterogeneity of PAX5 and IKZF1 deletions and the recurrent nature of RB1 deletions. Whole gene losses are often indicative of larger deletions, visible by conventional cytogenetics. An increased number of copy number abnormalities is associated with NCI high risk, specifically deletions of IKZF1 and CDKN2A/B, which occur more frequently among these patients. IKZF1 deletions and rearrangements of CRLF2 among patients with undefined karyotypes may point to the poor risk BCR-ABL1-like group. In conclusion, this study has demonstrated in a large representative cohort of children with B-cell precursor acute lymphoblastic leukemia that the pattern of copy number abnormalities is highly variable according to the primary genetic abnormality.

Phenotype prediction by integrative network analysis of SNP and gene expression microarrays

A long-term goal of biomedical research is to decipher how genetic processes influence disease formation. Ubiquitous and advancing microarray technology can measure millions of DNA structural variants (single-nucleotide polymorphisms, or SNPs) and thousands of gene transcripts (RNA expression microarrays) in cells. Both of these information modalities can be brought to bear on disease etiology. This paper develops a Bayesian network-based approach to integrate SNP and expression microarray data. The network models SNP-gene interactions using a phenotype-centric network. Inferring the network consists of two steps: variable selection and network learning. The learned network illustrates how functionally dependent SNPs and genes influence each other, and also serves as a predictor of the phenotype. The application of the proposed method to a pediatric acute lymphoblastic leukemia dataset demonstrates the feasibility of our approach and its impact on biological investigation and clinical practice.

IKAROS deletions dictate a unique gene expression signature in patients with adult B-cell acute lymphoblastic leukemia

BACKGROUND

Deletions of IKAROS (IKZF1) frequently occur in B-cell precursor acute lymphoblastic leukemia (B-ALL) but the mechanisms by which they influence pathogenesis are unclear. To address this issue, a cohort of 144 adult B-ALL patients (106 BCR-ABL1-positive and 38 B-ALL negative for known molecular rearrangements) was screened for IKZF1 deletions by single nucleotide polymorphism (SNP) arrays; a sub-cohort of these patients (44%) was then analyzed for gene expression profiling.

PRINCIPAL FINDINGS

Total or partial deletions of IKZF1 were more frequent in BCR-ABL1-positive than in BCR-ABL1-negative B-ALL cases (75% vs 58%, respectively, p = 0.04). Comparison of the gene expression signatures of patients carrying IKZF1 deletion vs those without showed a unique signature featured by down-regulation of B-cell lineage and DNA repair genes and up-regulation of genes involved in cell cycle, JAK-STAT signalling and stem cell self-renewal. Through chromatin immunoprecipitation and luciferase reporter assays we corroborated these findings both in vivo and in vitro, showing that Ikaros deleted isoforms lacked the ability to directly regulate a large group of the genes in the signature, such as IGLL1, BLK, EBF1, MSH2, BUB3, ETV6, YES1, CDKN1A (p21), CDKN2C (p18) and MCL1.

CONCLUSIONS

Here we identified and validated for the first time molecular pathways specifically controlled by IKZF1, shedding light into IKZF1 role in B-ALL pathogenesis.

Integrated analysis of gene copy number, copy neutral LOH, and microRNA profiles in adult acute lymphoblastic leukemia

We adopted an integrated analysis of gene copy number alterations (CNAs), copy number neutral loss of heterozygosity (CNN LOH), and microRNA (miRNA) profiling in 21 adult acute lymphoblastic leukemia (ALL) patients. This study revealed the most frequent CNAs to be at chromosomes 9p, 7, and 17 and recurrent CNN LOH at 5p, 9p, and Xq. As for the most differentially expressed miRNAs, they included 8 upregulated and 14 downregulated miRNAs, of which miR-148a at 7p15.2, miR-22 at 17p13.3, miR-223 at Xq12, as well as miR-101-2 at 9p24.1 exhibited recurrent CNAs or CNN LOH. miR-101-2 was recurrently downregulated, and although the related CNN LOH was detected only in BCR-ABL1 negative cases (2/14), deletions of miR-101-2 were observed solely in BCR-ABL1 positive cases (4/7). Finally, BCR-ABL1 positive cases, in contrast to negative ones, were characterized by slightly, but still significantly, higher expression levels of miR-29b.

Digital gene expression profiling of primary acute lymphoblastic leukemia cells

We determined the genome-wide digital gene expression (DGE) profiles of primary acute lymphoblastic leukemia (ALL) cells from 21 patients taking advantage of ‘second-generation' sequencing technology. Patients included in this study represent four cytogenetically distinct subtypes of B-cell precursor (BCP) ALL and T-cell lineage ALL (T-ALL). The robustness of DGE combined with supervised classification by nearest shrunken centroids (NSC) was validated experimentally and by comparison with published expression data for large sets of ALL samples. Genes that were differentially expressed between BCP ALL subtypes were enriched to distinct signaling pathways with dic(9;20) enriched to TP53 signaling, t(9;22) to interferon signaling, as well as high hyperdiploidy and t(12;21) to apoptosis signaling. We also observed antisense tags expressed from the non-coding strand of ∼50% of annotated genes, many of which were expressed in a subtype-specific pattern. Antisense tags from 17 gene regions unambiguously discriminated between the BCP ALL and T-ALL subtypes, and antisense tags from 76 gene regions discriminated between the 4 BCP subtypes. We observed a significant overlap of gene regions with alternative polyadenylation and antisense transcription (P<1 × 10⁻¹⁵). Our study using DGE profiling provided new insights into the RNA expression patterns in ALL cells.

The projection score--an evaluation criterion for variable subset selection in PCA visualization

BACKGROUND

In many scientific domains, it is becoming increasingly common to collect high-dimensional data sets, often with an exploratory aim, to generate new and relevant hypotheses. The exploratory perspective often makes statistically guided visualization methods, such as Principal Component Analysis (PCA), the methods of choice. However, the clarity of the obtained visualizations, and thereby the potential to use them to formulate relevant hypotheses, may be confounded by the presence of the many non-informative variables. For microarray data, more easily interpretable visualizations are often obtained by filtering the variable set, for example by removing the variables with the smallest variances or by only including the variables most highly related to a specific response. The resulting visualization may depend heavily on the inclusion criterion, that is, effectively the number of retained variables. To our knowledge, there exists no objective method for determining the optimal inclusion criterion in the context of visualization.

RESULTS

We present the projection score, which is a straightforward, intuitively appealing measure of the informativeness of a variable subset with respect to PCA visualization. This measure can be universally applied to find suitable inclusion criteria for any type of variable filtering. We apply the presented measure to find optimal variable subsets for different filtering methods in both microarray data sets and synthetic data sets. We note also that the projection score can be applied in general contexts, to compare the informativeness of any variable subsets with respect to visualization by PCA.

CONCLUSIONS

We conclude that the projection score provides an easily interpretable and universally applicable measure of the informativeness of a variable subset with respect to visualization by PCA, that can be used to systematically find the most interpretable PCA visualization in practical exploratory analysis.

Genomic characterization implicates iAMP21 as a likely primary genetic event in childhood B-cell precursor acute lymphoblastic leukemia

Intrachromosomal amplification of chromosome 21 (iAMP21) defines a distinct subgroup of childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) that has a dismal outcome when treated with standard therapy. For improved diagnosis and risk stratification, the initiating genetic events need to be elucidated. To investigate the genetic basis of BCP-ALL, genomes of 94 iAMP21 patients were interrogated by arrays, FISH, and multiplex ligation-dependent probe amplification. Most copy number alterations targeted chromosome 21, reinforcing the complexity of this chromosome. The common region of amplification on chromosome 21 was refined to a 5.1-mb region that included RUNX1, miR-802, and genes mapping to the Down syndrome critical region. Recurrent abnormalities affecting genes in key pathways were identified: IKZF1 (22%), CDKN2A/B (17%), PAX5 (8%), ETV6 (19%), and RB1 (37%). Investigation of clonal architecture provided evidence that these abnormalities, and P2RY8-CRLF2, were secondary to chromosome 21 rearrangements. Patient outcome was uniformly poor with standard therapy irrespective of the presence or absence of these changes. This study has provided evidence that chromosome 21 instability is the only anomaly among those so far investigated that is common to all iAMP21 patients, and therefore the initiating event is likely to be found among the complex structural rearrangements of this abnormal chromosome.

Mapping transcription mechanisms from multimodal genomic data

BACKGROUND

Identification of expression quantitative trait loci (eQTLs) is an emerging area in genomic study. The task requires an integrated analysis of genome-wide single nucleotide polymorphism (SNP) data and gene expression data, raising a new computational challenge due to the tremendous size of data.

RESULTS

We develop a method to identify eQTLs. The method represents eQTLs as information flux between genetic variants and transcripts. We use information theory to simultaneously interrogate SNP and gene expression data, resulting in a Transcriptional Information Map (TIM) which captures the network of transcriptional information that links genetic variations, gene expression and regulatory mechanisms. These maps are able to identify both cis- and trans- regulating eQTLs. The application on a dataset of leukemia patients identifies eQTLs in the regions of the GART, PCP4, DSCAM, and RIPK4 genes that regulate ADAMTS1, a known leukemia correlate.

CONCLUSIONS

The information theory approach presented in this paper is able to infer the dependence networks between SNPs and transcripts, which in turn can identify cis- and trans-eQTLs. The application of our method to the leukemia study explains how genetic variants and gene expression are linked to leukemia.

Integrative analysis of gene expression and copy number alterations using canonical correlation analysis

Background

With the rapid development of new genetic measurement methods, several types of genetic alterations can be quantified in a high-throughput manner. While the initial focus has been on investigating each data set separately, there is an increasing interest in studying the correlation structure between two or more data sets. Multivariate methods based on Canonical Correlation Analysis (CCA) have been proposed for integrating paired genetic data sets. The high dimensionality of microarray data imposes computational difficulties, which have been addressed for instance by studying the covariance structure of the data, or by reducing the number of variables prior to applying the CCA. In this work, we propose a new method for analyzing high-dimensional paired genetic data sets, which mainly emphasizes the correlation structure and still permits efficient application to very large data sets. The method is implemented by translating a regularized CCA to its dual form, where the computational complexity depends mainly on the number of samples instead of the number of variables. The optimal regularization parameters are chosen by cross-validation. We apply the regularized dual CCA, as well as a classical CCA preceded by a dimension-reducing Principal Components Analysis (PCA), to a paired data set of gene expression changes and copy number alterations in leukemia.

Results

Using the correlation-maximizing methods, regularized dual CCA and PCA+CCA, we show that without pre-selection of known disease-relevant genes, and without using information about clinical class membership, an exploratory analysis singles out two patient groups, corresponding to well-known leukemia subtypes. Furthermore, the variables showing the highest relevance to the extracted features agree with previous biological knowledge concerning copy number alterations and gene expression changes in these subtypes. Finally, the correlation-maximizing methods are shown to yield results which are more biologically interpretable than those resulting from a covariance-maximizing method, and provide different insight compared to when each variable set is studied separately using PCA.

Conclusions

We conclude that regularized dual CCA as well as PCA+CCA are useful methods for exploratory analysis of paired genetic data sets, and can be efficiently implemented also when the number of variables is very large.

Down syndrome acute lymphoblastic leukemia, a highly heterogeneous disease in which aberrant expression of CRLF2 is associated with mutated JAK2: a report from the International BFM Study Group

We report gene expression and other analyses to elucidate the molecular characteristics of acute lymphoblastic leukemia (ALL) in children with Down syndrome (DS). We find that by gene expression DS-ALL is a highly heterogeneous disease not definable as a unique entity. Nevertheless, 62% (33/53) of the DS-ALL samples analyzed were characterized by high expression of the type I cytokine receptor CRLF2 caused by either immunoglobulin heavy locus (IgH@) translocations or by interstitial deletions creating chimeric transcripts P2RY8-CRLF2. In 3 of these 33 patients, a novel activating somatic mutation, F232C in CRLF2, was identified. Consistent with our previous research, mutations in R683 of JAK2 were identified in 10 specimens (19% of the patients) and, interestingly, all 10 had high CRLF2 expression. Cytokine receptor-like factor 2 (CRLF2) and mutated Janus kinase 2 (Jak2) cooperated in conferring cytokine-independent growth to BaF3 pro-B cells. Intriguingly, the gene expression signature of DS-ALL is enriched with DNA damage and BCL6 responsive genes, suggesting the possibility of B-cell lymphocytic genomic instability. Thus, DS confers increased risk for genetically highly diverse ALLs with frequent overexpression of CRLF2, associated with activating mutations in the receptor itself or in JAK2. Our data also suggest that the majority of DS children with ALL may benefit from therapy blocking the CRLF2/JAK2 pathways.

DNA copy-number abnormalities do not occur in infant ALL with t(4;11)/MLL-AF4

The pathogenesis of infant acute lymphoblastic leukemia (ALL) is still not well defined. Short latency to leukemia and very high concordance rate for ALL in Mixed-Lineage Leukemia (MLL)-positive infant twins suggest that the MLL rearrangement itself could be sufficient for overt leukemia. Attempts to generate a suitable mouse model for MLL-AF4-positive ALL did not thoroughly resolve the issue of whether cooperating mutations are required to reduce latency and to generate overt leukemia in vivo. In this study, we applied single-nucleotide polymorphism array technology to perform genomic profiling of 28 infant ALL cases carrying t(4;11) to detect MLL-cooperating aberrations hidden to conventional techniques and to gain new insights into infant ALL pathogenesis. In contrast to pediatric, adolescent and adult ALL cases, the MLL rearrangement in infant ALL is associated with an exceptionally low frequency of copy-number abnormalities, thus confirming the unique nature of this disease. By contrast, additional genetic aberrations are acquired at disease relapse. Small-segmental uniparental disomy traits were frequently detected, mostly constitutional, and widely distributed throughout the genome. It can be argued that the MLL rearrangement as a first hit, rather than inducing the acquisition of additional genetic lesions, has a major role to drive and hasten the onset of leukemia.

Molecular inversion probes reveal patterns of 9p21 deletion and copy number aberrations in childhood leukemia

Childhood leukemia, which accounts for >30% of newly diagnosed childhood malignancies, is one of the leading causes of death for children with cancer. Genome-wide studies using microarray chips to identify copy number changes in human cancer are becoming more common. In this pilot study, 45 pediatric leukemia samples were analyzed for gene copy aberrations using novel molecular inversion probe (MIP) technology. Acute leukemia subtypes included precursor B-cell acute lymphoblastic leukemia (ALL) (n=23), precursor T-cell ALL (n=6), and acute myeloid leukemia (n=14). The MIP analysis identified 69 regions of recurring copy number changes, of which 41 have not been identified with other DNA microarray platforms. Copy number gains and losses were validated in 98% of clinical karyotypes and 100% of fluorescence in situ hybridization studies available. We report unique patterns of copy number loss in samples with 9p21.3 (CDKN2A) deletion in the precursor B-cell ALL patients, compared with the precursor T-cell ALL patients. MIPs represent an attractive technology for identifying novel copy number aberrations, validating previously reported copy number changes, and translating molecular findings into clinically relevant targets for further investigation.

Different molecular mechanisms causing 9p21 deletions in acute lymphoblastic leukemia of childhood

Deletion of chromosome 9p21 is a crucial event for the development of several cancers including acute lymphoblastic leukemia (ALL). Double strand breaks (DSBs) triggering 9p21 deletions in ALL have been reported to occur at a few defined sites by illegitimate action of the V(D)J recombination activating protein complex. We have cloned 23 breakpoint junctions for a total of 46 breakpoints in 17 childhood ALL (9 B- and 8 T-lineages) showing different size deletions at one or both homologous chromosomes 9 to investigate which particular sequences make the region susceptible to interstitial deletion. We found that half of 9p21 deletion breakpoints were mediated by ectopic V(D)J recombination mechanisms whereas the remaining half were associated to repeated sequences, including some with potential for non-B DNA structure formation. Other mechanisms, such as microhomology-mediated repair, that are common in other cancers, play only a very minor role in ALL. Nucleotide insertions at breakpoint junctions and microinversions flanking the breakpoints have been detected at 20/23 and 2/23 breakpoint junctions, respectively, both in the presence of recombination signal sequence (RSS)-like sequences and of other unspecific sequences. The majority of breakpoints were unique except for two cases, both T-ALL, showing identical deletions. Four of the 46 breakpoints coincide with those reported in other cases, thus confirming the presence of recurrent deletion hotspots. Among the six cases with heterozygous 9p deletions, we found that the remaining CDKN2A and CDKN2B alleles were hypermethylated at CpG islands.