Increased MYB alternative promoter usage is associated with relapse in acute lymphoblastic leukemia.

Therapy-resistant disease is a major cause of death in patients with acute lymphoblastic leukemia (ALL). Activation of the MYB oncogene is associated with ALL and leads to uncontrolled neoplastic cell proliferation and blocked differentiation. Here, we used RNA-seq to study the clinical significance of MYB expression and MYB alternative promoter (TSS2) usage in 133 pediatric ALLs. RNA-seq revealed that all cases analyzed overexpressed MYB and demonstrated MYB TSS2 activity. qPCR analyses confirmed the expression of the alternative MYB promoter also in seven ALL cell lines. Notably, high MYB TSS2 activity was significantly associated with relapse (p = 0.007). Moreover, cases with high MYB TSS2 usage showed evidence of therapy-resistant disease with increased expression of ABC multidrug resistance transporter genes (e.g., ABCA2, ABCB5, and ABCC10) and enzymes catalyzing drug degradation (e.g., CYP1A2, CYP2C9, and CYP3A5). Elevated MYB TSS2 activity was further associated with augmented KRAS signaling (p < 0.05) and decreased methylation of the conventional MYB promoter (p < 0.01). Taken together, our results suggest that MYB alternative promoter usage is a novel potential prognostic biomarker for relapse and therapy resistance in pediatric ALL.

CD99 expression is strongly associated with clinical outcome in children with B-cell precursor acute lymphoblastic leukaemia.

Our study aimed to determine the expression pattern and clinical relevance of CD99 in paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL). Our findings demonstrate that high expression levels of CD99 are mainly found in high-risk BCP-ALL, e.g. BCR-ABL1 and CRLF2Re/Hi, and that high CD99 mRNA levels are strongly associated with a high frequency of relapse, high proportion of positive for minimal residual disease at day 29 and poor overall survival in paediatric cohorts, which indicate that CD99 is a potential biomarker for BCP-ALL.

CopyNumber450kCancer: baseline correction for accurate copy number calling from the 450k methylation array.

UNLABELLED: The Illumina Infinium HumanMethylation450 BeadChip (450k) is widely used for the evaluation of DNA methylation levels in large-scale datasets, particularly in cancer. The 450k design allows copy number variant (CNV) calling using existing bioinformatics tools. However, in cancer samples, numerous large-scale aberrations cause shifting in the probe intensities and thereby may result in erroneous CNV calling. Therefore, a baseline correction process is needed. We suggest the maximum peak of probe segment density to correct the shift in the intensities in cancer samples. AVAILABILITY AND IMPLEMENTATION: CopyNumber450kCancer is implemented as an R package. The package with examples can be downloaded at http://cran.r-project.org CONTACT: nour.marzouka@medsci.uu.se SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The mutational landscape in pediatric acute lymphoblastic leukemia deciphered by whole genome sequencing.

Genomic characterization of pediatric acute lymphoblastic leukemia (ALL) has identified distinct patterns of genes and pathways altered in patients with well-defined genetic aberrations. To extend the spectrum of known somatic variants in ALL, we performed whole genome and transcriptome sequencing of three B-cell precursor patients, of which one carried the t(12;21)ETV6-RUNX1 translocation and two lacked a known primary genetic aberration, and one T-ALL patient. We found that each patient had a unique genome, with a combination of well-known and previously undetected genomic aberrations. By targeted sequencing in 168 patients, we identified KMT2D and KIF1B as novel putative driver genes. We also identified a putative regulatory non-coding variant that coincided with overexpression of the growth factor MDK. Our results contribute to an increased understanding of the biological mechanisms that lead to ALL and suggest that regulatory variants may be more important for cancer development than recognized to date. The heterogeneity of the genetic aberrations in ALL renders whole genome sequencing particularly well suited for analysis of somatic variants in both research and diagnostic applications.

Refining risk prediction in pediatric acute lymphoblastic leukemia through DNA methylation profiling.

Acute lymphoblastic leukemia (ALL) is the most prevalent cancer in children, and despite considerable progress in treatment outcomes, relapses still pose significant risks of mortality and long-term complications. To address this challenge, we employed a supervised machine learning technique, specifically random survival forests, to predict the risk of relapse and mortality using array-based DNA methylation data from a cohort of 763 pediatric ALL patients treated in Nordic countries. The relapse risk predictor (RRP) was constructed based on 16 CpG sites, demonstrating c-indexes of 0.667 and 0.677 in the training and test sets, respectively. The mortality risk predictor (MRP), comprising 53 CpG sites, exhibited c-indexes of 0.751 and 0.754 in the training and test sets, respectively. To validate the prognostic value of the predictors, we further analyzed two independent cohorts of Canadian (n = 42) and Nordic (n = 384) ALL patients. The external validation confirmed our findings, with the RRP achieving a c-index of 0.667 in the Canadian cohort, and the RRP and MRP achieving c-indexes of 0.529 and 0.621, respectively, in an independent Nordic cohort. The precision of the RRP and MRP models improved when incorporating traditional risk group data, underscoring the potential for synergistic integration of clinical prognostic factors. The MRP model also enabled the definition of a risk group with high rates of relapse and mortality. Our results demonstrate the potential of DNA methylation as a prognostic factor and a tool to refine risk stratification in pediatric ALL. This may lead to personalized treatment strategies based on epigenetic profiling.

Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment.

BACKGROUND: Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing. RESULTS: We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792-1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools. CONCLUSION: Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

Presence of FLT3-ITD and high BAALC expression are independent prognostic markers in childhood acute myeloid leukemia.

Mutation status of FLT3, NPM1, CEBPA, and WT1 genes and gene expression levels of ERG, MN1, BAALC, FLT3, and WT1 have been identified as possible prognostic markers in acute myeloid leukemia (AML). We have performed a thorough prognostic evaluation of these genetic markers in patients with pediatric AML enrolled in the Nordic Society of Pediatric Hematology and Oncology (NOPHO) 1993 or NOPHO 2004 protocols. Mutation status and expression levels were analyzed in 185 and 149 patients, respectively. Presence of FLT3-internal tandem duplication (ITD) was associated with significantly inferior event-free survival (EFS), whereas presence of an NPM1 mutation in the absence of FLT3-ITD correlated with significantly improved EFS. Furthermore, high levels of ERG and BAALC transcripts were associated with inferior EFS. No significant correlation with survival was seen for mutations in CEBPA and WT1 or with gene expression levels of MN1, FLT3, and WT1. In multivariate analysis, the presence of FLT3-ITD and high BAALC expression were identified as independent prognostic markers of inferior EFS. We conclude that analysis of the mutational status of FLT3 and NPM1 at diagnosis is important for prognostic stratification of patients with pediatric AML and that determination of the BAALC gene expression level can add valuable information.

Multimodal classification of molecular subtypes in pediatric acute lymphoblastic leukemia.

Genomic analyses have redefined the molecular subgrouping of pediatric acute lymphoblastic leukemia (ALL). Molecular subgroups guide risk-stratification and targeted therapies, but outcomes of recently identified subtypes are often unclear, owing to limited cases with comprehensive profiling and cross-protocol studies. We developed a machine learning tool (ALLIUM) for the molecular subclassification of ALL in retrospective cohorts as well as for up-front diagnostics. ALLIUM uses DNA methylation and gene expression data from 1131 Nordic ALL patients to predict 17 ALL subtypes with high accuracy. ALLIUM was used to revise and verify the molecular subtype of 281 B-cell precursor ALL (BCP-ALL) cases with previously undefined molecular phenotype, resulting in a single revised subtype for 81.5% of these cases. Our study shows the power of combining DNA methylation and gene expression data for resolving ALL subtypes and provides a comprehensive population-based retrospective cohort study of molecular subtype frequencies in the Nordic countries.

DNA methylation for subtype classification and prediction of treatment outcome in patients with childhood acute lymphoblastic leukemia.

Despite improvements in the prognosis of childhood acute lymphoblastic leukemia (ALL), subgroups of patients would benefit from alternative treatment approaches. Our aim was to identify genes with DNA methylation profiles that could identify such groups. We determined the methylation levels of 1320 CpG sites in regulatory regions of 416 genes in cells from 401 children diagnosed with ALL. Hierarchical clustering of 300 CpG sites distinguished between T-lineage ALL and B-cell precursor (BCP) ALL and between the main cytogenetic subtypes of BCP ALL. It also stratified patients with high hyperdiploidy and t(12;21) ALL into 2 subgroups with different probability of relapse. By using supervised learning, we constructed multivariate classifiers by external cross-validation procedures. We identified 40 genes that consistently contributed to accurate discrimination between the main subtypes of BCP ALL and gene sets that discriminated between subtypes of ALL and between ALL and controls in pairwise classification analyses. We also identified 20 individual genes with DNA methylation levels that predicted relapse of leukemia. Thus, methylation analysis should be explored as a method to improve stratification of ALL patients. The genes highlighted in our study are not enriched to specific pathways, but the gene expression levels are inversely correlated to the methylation levels.

Minimal residual disease assessment in childhood acute lymphoblastic leukaemia: a Swedish multi-centre study comparing real-time polymerase chain reaction and multicolour flow cytometry.

Minimal residual disease (MRD) assessment is a powerful prognostic factor for determining the risk of relapse in childhood acute lymphoblastic leukaemia (ALL). In this Swedish multi-centre study of childhood ALL diagnosed between 2002 and 2006, the MRD levels were analysed in 726 follow-up samples in 228 children using real-time quantitative polymerase chain reaction (RQ-PCR) of rearranged immunoglobulin/T-cell receptor genes and multicolour flow cytometry (FCM). Using an MRD threshold of 0·1%, which was the sensitivity level reached in all analyses, the concordance between RQ-PCR and FCM MRD values at day 29 was 84%. In B-cell precursor ALL, an MRD level of ≥0·1% at day 29 predicted a higher risk of bone marrow relapse (BMR) with both methods, although FCM was a better discriminator. However, considering the higher median MRD values achieved with RQ-PCR, a higher MRD cut-off (≥0·2%) improved the predictive capacity of RQ-PCR. In T-ALL, RQ-PCR was notably superior to FCM in predicting risk of BMR. That notwithstanding, MRD levels of ≥0·1%, detected by either method at day 29, could not predict isolated extramedullary relapse. In conclusion, the concordance between RQ-PCR and FCM was high and hence both methods are valuable clinical tools for identifying childhood ALL cases with increased risk of BMR.

RAG1 co-expression signature identifies ETV6-RUNX1-like B-cell precursor acute lymphoblastic leukemia in children.

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) can be classified into subtypes according to the genetic aberrations they display. For instance, the translocation t(12;21)(p13;q22), representing the ETV6-RUNX1 fusion gene (ER), is present in a quarter of BCP-ALL cases. However, around 10% of the cases lack classifying chromosomal abnormalities (B-other). In pediatric ER BCP-ALL, rearrangement mediated by RAG (recombination-activating genes) has been proposed as the predominant driver of oncogenic rearrangement. Herein we analyzed almost 1600 pediatric BCP-ALL samples to determine which subtypes express RAG. We demonstrate that RAG1 mRNA levels are especially high in the ETV6-RUNX1 (ER) subtype and in a subset of B-other samples. We also define 31 genes that are co-expressed with RAG1 (RAG1-signature) in the ER subtype, a signature that also identifies this subset of B-other samples. Moreover, this subset also shares leukemia and pro-B gene expression signatures as well as high levels of the ETV6 target genes (BIRC7, WBP1L, CLIC5, ANGPTL2) with the ER subtype, indicating that these B-other cases are the recently identified ER-like subtype. We validated our results in a cohort where ER-like has been defined, which confirmed expression of the RAG1-signature in this recently described subtype. Taken together, our results demonstrate that the RAG1-signature identifies the ER-like subtype. As there are no definitive genetic markers to identify this novel subtype, the RAG1-signature represents a means to screen for this leukemia in children.

Mutational patterns and clonal evolution from diagnosis to relapse in pediatric acute lymphoblastic leukemia.

The mechanisms driving clonal heterogeneity and evolution in relapsed pediatric acute lymphoblastic leukemia (ALL) are not fully understood. We performed whole genome sequencing of samples collected at diagnosis, relapse(s) and remission from 29 Nordic patients. Somatic point mutations and large-scale structural variants were called using individually matched remission samples as controls, and allelic expression of the mutations was assessed in ALL cells using RNA-sequencing. We observed an increased burden of somatic mutations at relapse, compared to diagnosis, and at second relapse compared to first relapse. In addition to 29 known ALL driver genes, of which nine genes carried recurrent protein-coding mutations in our sample set, we identified putative non-protein coding mutations in regulatory regions of seven additional genes that have not previously been described in ALL. Cluster analysis of hundreds of somatic mutations per sample revealed three distinct evolutionary trajectories during ALL progression from diagnosis to relapse. The evolutionary trajectories provide insight into the mutational mechanisms leading relapse in ALL and could offer biomarkers for improved risk prediction in individual patients.

Refined detection and phasing of structural aberrations in pediatric acute lymphoblastic leukemia by linked-read whole-genome sequencing.

Structural chromosomal rearrangements that can lead to in-frame gene-fusions are a leading source of information for diagnosis, risk stratification, and prognosis in pediatric acute lymphoblastic leukemia (ALL). Traditional methods such as karyotyping and FISH struggle to accurately identify and phase such large-scale chromosomal aberrations in ALL genomes. We therefore evaluated linked-read WGS for detecting chromosomal rearrangements in primary samples of from 12 patients diagnosed with ALL. We assessed the effect of input DNA quality on phased haplotype block size and the detectability of copy number aberrations and structural variants in the ALL genomes. We found that biobanked DNA isolated by standard column-based extraction methods was sufficient to detect chromosomal rearrangements even at low 10x sequencing coverage. Linked-read WGS enabled precise, allele-specific, digital karyotyping at a base-pair resolution for a wide range of structural variants including complex rearrangements and aneuploidy assessment. With use of haplotype information from the linked-reads, we also identified previously unknown structural variants, such as a compound heterozygous deletion of ERG in a patient with the DUX4-IGH fusion gene. We conclude that linked-read WGS allows detection of important pathogenic variants in ALL genomes at a resolution beyond that of traditional karyotyping and FISH.

Overexpression of chromatin remodeling and tyrosine kinase genes in iAMP21-positive acute lymphoblastic leukemia.

Intrachromosomal amplification of chromosome 21 (iAMP21) is a cytogenetic subtype associated with relapse and poor prognosis in pediatric B-cell precursor acute lymphoblastic leukemia (BCP ALL). The biology behind the high relapse risk is unknown and the aim of this study was to further characterize the genomic and transcriptional landscape of iAMP21. Using DNA arrays and sequencing, we could identify rearrangements and aberrations characteristic for iAMP21. RNA sequencing revealed that only half of the genes in the minimal region of amplification (20/45) were differentially expressed in iAMP21. Among them were the top overexpressed genes (p < 0.001) in iAMP21 vs. BCP ALL without iAMP21 and three candidate genes could be identified, the tyrosine kinase gene DYRK1A and chromatin remodeling genes CHAF1B and SON. While overexpression of DYRK1A and CHAF1B is associated with poor prognosis in malignant diseases including myeloid leukemia, this is the first study to show significant correlation with iAMP21-positive ALL.

The proliferative history shapes the DNA methylome of B-cell tumors and predicts clinical outcome.

We report a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. Differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

In vitro cellular drug sensitivity at diagnosis is correlated to minimal residual disease at end of induction therapy in childhood acute lymphoblastic leukemia.

Leukemic cells from 85 children with newly diagnosed precursor B-lineage ALL were tested for in vitro drug sensitivity to a panel of anti-cancer drugs. Minimal residual disease (MRD) was measured by RQ-PCR. There was a significant correlation between MRD day 29 and in vitro sensitivity to prednisolone (p<0.001) and doxorubicin (p=0.017), drugs administered during induction therapy. In patients with t(12;21) (n=20), in vitro sensitivity to doxorubicin was an independent factor for MRD <0.1% (p=0.031; R(2)=0.66). Thus, data show that in vitro drug sensitivity at diagnosis is correlated to cell kill during induction therapy as measured by MRD day 29.

Thioguanine pharmacokinetics in induction therapy of children with acute myeloid leukemia.

We studied the pharmacokinetics of 6-thioguanine (6TG) in 50 children treated for newly diagnosed acute myeloid leukemia, four of them with Down syndrome (DS). They received oral 6TG 100 mg/m2 body surface area twice daily for 4 days. Etoposide, 100 mg/m2/24 h, and cytarabine, 200 mg/m2/24 h, were administered concomitantly by intravenous infusion. On day 5, doxorubicin 75 mg/m2 was given as an 8-h infusion. The concentration of thioguanine nucleotides (TGN) in erythrocytes, the active metabolites of 6TG, was determined by high-performance liquid chromatography. The mean TGN concentration from 72, 95, and 106-h samples was used as a measure of drug exposure for each individual. The median TGN concentration in non-DS children above 2 years of age was 2.30 micromol/mmol Hb (range 0.57-25.3). The TGN concentrations varied widely (30-fold) also after dose normalization. We found no correlation with demographic, clinical, or biochemical parameters, and differences in bioavailability might be the most important explanation to interpatient variability. Children with high TGN concentration tended to have longer treatment interval to the next course, but we found no correlation with our predefined parameters for clinical response, that is, remission and relapse rate. Therefore, 6TG does not seem to be a candidate for therapeutic drug monitoring by TGN measurement, at least not in the setting of short multidrug treatment courses. Children with DS had significantly higher TGN concentrations, indicating that dose reduction might be considered to reach the same drug exposure as in non-DS children.

Vincristine pharmacokinetics in children with Down syndrome.

Children with Down syndrome (DS), who represent about 2% of childhood acute lymphoblastic leukemia, have inferior prognosis compared to non-DS children. For vincristine (and many other anticancer agents) pharmacokinetic data are scant or missing, and there is considerable uncertainty about the optimal dosing of drugs to patients with DS. We studied vincristine pharmacokinetics on treatment day one in six children with DS and compared to 92 non-DS children. No differences were found. Thus, we found no rationale for dose reduction of vincristine in DS children from a strictly pharmacokinetic point of view.

Pharmacokinetics of high-dose methotrexate in infants treated for acute lymphoblastic leukemia.

BACKGROUND: Interfant-99 was an international collaborative treatment protocol for infants with acute lymphoblastic leukemia (ALL). PROCEDURE: We collected data on 103 infants at the time of their first treatment with high-dose methotrexate (HD MTX), 5 g/m(2). Children <6 months of age received two-third of the calculated dose based on body surface area (BSA), children 6-12 months three-fourth of the calculated dose, and children >12 months full dose. RESULTS: The median steady-state MTX concentration at the end of the 24-hr infusion was 57.8 microM (range 9.5-313). The median systemic clearance was 6.22 L/hr/m(2) BSA, and tended to increase with age (P = 0.099). Boys had higher clearance than girls, 6.77 and 5.28 L/hr/m(2) (P = 0.030), and tended to have lower median MTX concentration at 24 hr. Eight infants had MTX levels below 20 microM, a level judged to be sufficient in B-lineage ALL in children >1 year of age. All infants tolerated the dose well enough to receive a second dose of HD MTX without dose reduction. We found no significant effect on disease-free survival for MTX steady-state concentration, MTX clearance, or time to MTX below 0.2 microM. CONCLUSIONS: Our data provide no support for a change in the dosing rules for MTX used in Interfant-99. However, in view of the poor treatment results for infants, one might consider increase in the dose for patients who reach plasma levels below median after the first MTX dose.