Integration of high-resolution methylome and transcriptome analyses to dissect epigenomic changes in childhood acute lymphoblastic leukemia

Commentary on "Epigenome-wide analysis across the development span of pediatric acute lymphoblastic leukemia: backtracking to birth"

VTRNA2-1 is a polymorphically imprinted locus. The proportion of individuals with a maternally imprinted VTRNA2-1 locus is consistently approximately 75% in populations of European origin, with the remaining circa 25% having a non-methylated VTRNA2-1 locus. Recently, VTRNA2-1 hypermethylation at birth was suggested to be a precursor of paediatric acute lymphoblastic leukaemia with biomarker potential. The results presented by Ghantous et al. [1] allow for an alternative interpretation to what the authors discussed, and we argue that the observed methylation difference at birth is due to an uneven distribution of imprinted and non-methylated individuals among the cases and controls, with all individuals presenting normative physiological VTRNA2-1 methylation levels. In addition, the notable interindividual variation arising from the polymorphic imprinting in VTRNA2-1 methylation levels calls into question the validity of VTRNA2-1 methylation as a biomarker.

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.

Prediction of disease-free survival of N1/2 non-small cell lung cancer after adjuvant chemotherapy by the biomarker RPMB

No molecular biomarkers have been proven applicable in clinical practice to identify patients who can benefit from adjuvant chemotherapy in non-small cell lung cancer (NSCLC). In this study, we established a biomarker, RPMB, short for promotor methylation burden of DNA repair genes (DRGs), to identify the subgroup of patients who might benefit from adjuvant chemotherapy in NSCLC. Methylation profiles of 828 NSCLC primary tumors and their clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. The RPMB for each patient after radical resection was calculated and its correlation with the prognosis of NSCLC was extensively investigated. DRGs of NSCLC were much more hypomethylated than the other genes (all p＜0.001). RPMB was defined as the ratio of methylated DRGs to the total number of all the DRGs. Patients with higher RPMB values tended to be nonsmokers, had adenocarcinoma, were female and had peripheral tumors. Subgroup analysis of forest plot among different clinical factors showed that high RPMB was significantly correlated to better disease-free survival (DFS) in pathologic N-positive patients after adjuvant chemotherapy (HR = 0.404, n = 62, p = 0.034). Notably, more superior DFS was exhibited in high RPMB NSCLCs with N1 nodal stage compared with those with low RPMB values (HR = 0.348, n = 47, p = 0.043). High RPMB might be used as a potential predictor to identify suitable N-positive NSCLC patients who can benefit from adjuvant chemotherapy after radical surgery.

MicroRNAs and the Diagnosis of Childhood Acute Lymphoblastic Leukemia: Systematic Review, Meta-Analysis and Re-Analysis with Novel Small RNA-Seq Tools

Simple Summary

MicroRNAs (miRNAs) have been under the spotlight for the last three decades. These non-coding RNAs seem to be dynamic regulators of mRNA stability and translation, in addition to interfering with transcription. Circulating miRNAs play a critical role in cell-to-cell interplay; therefore, they can serve as disease biomarkers. Meta-analysis of published data revealed that the CC genotype of rs4938723 in pri-miR-34b/c and the TT genotype of rs543412 in miR-100 confer protection against acute lymphoblastic leukemia (ALL) in children. Reanalysis of small RNA-seq data with novel tools identified significantly overexpressed members of the miR-128, miR-181, miR-130 and miR-17 families and significantly lower expression of miR-30, miR-24-2 and miR143~145 clusters, miR-574 and miR-618 in pediatric T-ALL cases compared with controls. Inconsistencies in methodology and study designs in most published material preclude reproducibility, and further cohort studies need to be conducted in order to empower novel tools, such as ALLSorts and RNAseqCNV.

Abstract

MicroRNAs (miRNAs) have been implicated in childhood acute lymphoblastic leukemia (ALL) pathogenesis. We performed a systematic review and meta-analysis of miRNA single-nucleotide polymorphisms (SNPs) in childhood ALL compared with healthy children, which revealed (i) that the CC genotype of rs4938723 in pri-miR-34b/c and the TT genotype of rs543412 in miR-100 confer protection against ALL occurrence in children; (ii) no significant association between rs2910164 genotypes in miR-146a and childhood ALL; and (iii) SNPs in DROSHA, miR-449b, miR-938, miR-3117 and miR-3689d-2 genes seem to be associated with susceptibility to B-ALL in childhood. A review of published literature on differential expression of miRNAs in children with ALL compared with controls revealed a significant upregulation of the miR-128 family, miR-130b, miR-155, miR-181 family, miR-210, miR-222, miR-363 and miR-708, along with significant downregulation of miR-143 and miR-148a, seem to have a definite role in childhood ALL development. MicroRNA signatures among childhood ALL subtypes, along with differential miRNA expression patterns between B-ALL and T-ALL cases, were scrutinized. With respect to T-ALL pediatric cases, we reanalyzed RNA-seq datasets with a robust and sensitive pipeline and confirmed the significant differential expression of hsa-miR-16-5p, hsa-miR-19b-3p, hsa-miR-92a-2-5p, hsa-miR-128-3p (ranked first), hsa-miR-130b-3p and -5p, hsa-miR-181a-5p, -2-3p and -3p, hsa-miR-181b-5p and -3p, hsa-miR-145-5p and hsa-miR-574-3p, as described in the literature, along with novel identified miRNAs.

Profiling chromatin accessibility in pediatric acute lymphoblastic leukemia identifies subtype-specific chromatin landscapes and gene regulatory networks

Acute lymphoblastic leukemia (ALL) is a hematopoietic malignancy comprised of molecular subtypes largely characterized by aneuploidy or recurring chromosomal rearrangements. Despite extensive information on the ALL transcriptome and methylome, there is limited understanding of the ALL chromatin landscape. We therefore mapped accessible chromatin in 24 primary ALL cell biospecimens comprising three common molecular subtypes (DUX4/ERG, ETV6-RUNX1 and hyperdiploid) from patients treated at St. Jude Children's Research Hospital. Our findings highlight extensive chromatin reprogramming in ALL, including the identification ALL subtype-specific chromatin landscapes that are additionally modulated by genetic variation. Chromatin accessibility differences between ALL and normal B-cells implicate the activation of B-cell repressed chromatin domains and detail the disruption of normal B-cell development in ALL. Among ALL subtypes, we uncovered roles for basic helix-loop-helix, homeodomain and activator protein 1 transcription factors in promoting subtype-specific chromatin accessibility and distinct gene regulatory networks. In addition to chromatin subtype-specificity, we further identified over 3500 DNA sequence variants that alter the ALL chromatin landscape and contribute to inter-individual variability in chromatin accessibility. Collectively, our data suggest that subtype-specific chromatin landscapes and gene regulatory networks impact ALL biology and contribute to transcriptomic differences among ALL subtypes.

Isoform-specific effects of transcription factor TCFL5 on the pluripotency-related genes SOX2 and KLF4 in colorectal cancer development

Colorectal cancer (CRC) is a very common life‐threatening malignancy. Transcription factor‐like 5 (TCFL5) has been suggested to be involved in CRC. Here, we describe the expression of four alternative transcripts of TCFL5 and their relevance in CRC. Complete deletion of all isoforms drastically decreased pro‐tumoural properties such as spheroids formation and in vivo tumour growth, although increased migration in CRC cell lines. Overexpression of the two main isoforms, TCFL5_E8 and CHA, had opposite effects: TCFL5_E8 reduced proliferation and spheroids formation, while CHA increased them. TCFL5_E8 reduced in vivo tumour formation, while CHA had no effect. In addition, TCFL5_E8 and CHA have different roles in the regulation of the pluripotency‐related genes SOX2 and KLF4. Both isoforms bind directly to their promoters; however, TCFL5_E8 induced SOX2 and reduced KLF4 mRNA levels, whereas CHA did the opposite. Together, our results show that TCFL5 plays an important role in the development of CRC, being however isoform‐specific. This work also points to the need to analyse separately TCFL5 isoforms in cancer, due to their different and opposite functions.

Single cell characterization of B-lymphoid differentiation and leukemic cell states during chemotherapy in ETV6-RUNX1-positive pediatric leukemia identifies drug-targetable transcription factor activities

BACKGROUND

Tight regulatory loops orchestrate commitment to B cell fate within bone marrow. Genetic lesions in this gene regulatory network underlie the emergence of the most common childhood cancer, acute lymphoblastic leukemia (ALL). The initial genetic hits, including the common translocation that fuses ETV6 and RUNX1 genes, lead to arrested cell differentiation. Here, we aimed to characterize transcription factor activities along the B-lineage differentiation trajectory as a reference to characterize the aberrant cell states present in leukemic bone marrow, and to identify those transcription factors that maintain cancer-specific cell states for more precise therapeutic intervention.

METHODS

We compared normal B-lineage differentiation and in vivo leukemic cell states using single cell RNA-sequencing (scRNA-seq) and several complementary genomics profiles. Based on statistical tools for scRNA-seq, we benchmarked a workflow to resolve transcription factor activities and gene expression distribution changes in healthy bone marrow lymphoid cell states. We compared these to ALL bone marrow at diagnosis and in vivo during chemotherapy, focusing on leukemias carrying the ETV6-RUNX1 fusion.

RESULTS

We show that lymphoid cell transcription factor activities uncovered from bone marrow scRNA-seq have high correspondence with independent ATAC- and ChIP-seq data. Using this comprehensive reference for regulatory factors coordinating B-lineage differentiation, our analysis of ETV6-RUNX1-positive ALL cases revealed elevated activity of multiple ETS-transcription factors in leukemic cells states, including the leukemia genome-wide association study hit ELK3. The accompanying gene expression changes associated with natural killer cell inactivation and depletion in the leukemic immune microenvironment. Moreover, our results suggest that the abundance of G1 cell cycle state at diagnosis and lack of differentiation-associated regulatory network changes during induction chemotherapy represent features of chemoresistance. To target the leukemic regulatory program and thereby overcome treatment resistance, we show that inhibition of ETS-transcription factors reduced cell viability and resolved pathways contributing to this using scRNA-seq.

CONCLUSIONS

Our data provide a detailed picture of the transcription factor activities characterizing both normal B-lineage differentiation and those acquired in leukemic bone marrow and provide a rational basis for new treatment strategies targeting the immune microenvironment and the active regulatory network in leukemia.

Promoter Methylation of DNA Repair Genes Predicts Disease-free Survival of Gastric Adenocarcinoma after Adjuvant Radiotherapy

The relentless debate on postoperative adjuvant radiotherapy in gastric adenocarcinoma (GA) has been lasting for decades. In this study, a new biomarker, named promoter methylation burden of DNA repair genes (RPMB), was established to identify the subgroup of patients who might benefit from adjuvant radiotherapy. Methylation profiles of 397 GA tumor samples were downloaded from The Cancer Genome Atlas (TCGA). RPMB for a patient was defined as the ratio of methylated DNA repair genes to the number of all DNA repair genes. Subgroup analyses in term of overall survival (OS) and disease-free survival (DFS) indicated that most of the subgroups favored the high-RMPB group. Kaplan-Meier analysis showed that overall the patients with high RPMB after R0 resection had a significantly better clinical outcome regarding DFS (hazard ratio [HR] = 0.013, p = 0.042). Additionally, high-RPMB patients, who underwent adjuvant radiotherapy with both ≥T2 tumor and positive lymph nodes, showed superior DFS in comparison with the low-RPMB group (HR = 5.35 × 10⁻¹⁰, n = 26, p = 0.010). RPMB might be considered as a promising biomarker for decision-making with regard to postoperative adjuvant radiotherapy for GA patients.

Clinicopathological features and prognostic value of SOX11 in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia is marked by aberrant transcriptional features that alter cell differentiation, self-renewal, and proliferative features. We sought to identify the transcription factors exhibiting altered and subtype-specific expression patterns in B-ALL and report here that SOX11, a developmental and neuronal transcription factor, is aberrantly expressed in the ETV6-RUNX1 and TCF3-PBX1 subtypes of acute B-cell leukemias. We show that a high expression of SOX11 leads to alterations of gene expression that are typically associated with cell adhesion, migration, and differentiation. A high expression is associated with DNA hypomethylation at the SOX11 locus and a favorable outcome. The results indicate that SOX11 expression marks a group of patients with good outcomes and thereby prompts further study of its use as a biomarker.

DNA methyltransferase 1-mediated CpG methylation of the miR-150-5p promoter contributes to fibroblast growth factor receptor 1-driven leukemogenesis

MicroRNA-150-5p (miR-150-5p) plays a complex role in normal early hematopoietic development and is also implicated in the development of various different leukemias. We have reported previously that, in myeloid and lymphoid malignancies associated with dysregulated fibroblast growth factor receptor 1 (FGFR1) activities, miR-150-5p is down-regulated compared with healthy cells. Here, using murine cells, we found that this down-regulation is accompanied by CpG methylation of the miR-150-5p promoter region. Of note, analysis of human acute lymphoblastic leukemia (ALL) cohorts also revealed an inverse relationship between miR-150-5p expression and disease progression. We also found that the DNA methyltransferase 1 (DNMT1) enzyme is highly up-regulated in FGFR1-driven leukemias and lymphomas and that FGFR1 inhibition reduces DNMT1 expression. DNMT1 knockdown in stem cell leukemia/lymphoma (SCLL) cells increased miR-150-5p levels and reduced levels of the MYB proto-oncogene transcription factor, a key regulator of leukemogenesis. FGFR1 directly activates the MYC proto-oncogene basic helix-loop-helix transcription factor, which, as we show here, binds and activates the DNMT1 promoter. MYC knockdown decreased DNMT1 expression, which, in turn, increased miR-150-5p expression. One of the known targets of miR-150-5p is MYB, and treatment of leukemic cells with the MYB inhibitor mebendazole dose-dependently increased apoptosis and reduced cell viability. Moreover, mebendazole treatment of murine xenografts models of FGFR1-driven leukemias enhanced survival. These findings provide evidence that MYC activates MYB by up-regulating DNMT1, which silences miR-150-5p and promotes SCLL progression. We propose that inclusion of mebendazole in a combination therapy with FGFR1 inhibitors may be a valuable option to manage SCLL.

The study of METTL3 and METTL14 expressions in childhood ETV6/RUNX1-positive acute lymphoblastic leukemia

BACKGROUND

This study was aimed to explore the METTL3 and METTL14 expressions in children with ETV6/RUNX1(E/R)-positive acute lymphoblastic leukemia (ALL) and investigate the relation between the METTL3 and METTL14 expressions with clinical features.

METHODS

Thirty-seven newly diagnosed E/R-positive ALL children and six controls were included in this study. Real-time quantitative polymerase chain reaction (RT-PCR) was used to detect the mRNA expression level of METTL3 and METTL14.

RESULTS

Among the 37 cases, 51.35% (n = 19) were boys and 48.65% (n = 18) were girls and the median age was 4.72 (1.72-11.99) years. Among the six controls, 50% (n = 3) were boys and 50% (n = 3) were girls and the median age was 5.24 (1.53-13.17) years. The expression level of METTL3 and METTL14 in E/R-positive ALL patients were lower than in controls (p < .05). Although failed to achieve statistical significance, the expression level of METTL3 and METTL14 in relapse patients were lower than nonrelapse patients (p = .171, p = .150, respectively).

CONCLUSION

The reduced levels of METTL3 and METTL14 suggest a possible role in the pathogenesis and course of E/R-positive ALL. METTL3 and METTL14 may become new prognostic factors, and rationalize specific treatment intensification in possible E/R-positive relapse patients.

Identification of survival‑associated key genes and long non‑coding RNAs in glioblastoma multiforme by weighted gene co‑expression network analysis

Glioblastoma multiforme (GBM) is the most common and aggressive primary malignant brain tumour. However, the causes of GBM are not clear, and the prognosis remains poor. The aim of the present study was to elucidate the key coding genes and long non‑coding RNAs (lncRNAs) associated with the survival time of GBM patients by obtaining the RNA expression profiles from the Chinese Glioma Genome Atlas database and conducting weighted gene co‑expression network analysis. Modules associated with overall survival (OS) were identified, and Gene Ontology and pathway enrichment analyses were performed. The hub genes of these modules were validated via survival analysis, while the biological functions of crucial lncRNAs were also analysed in the publicly available data. The results identified a survival‑associated module with 195 key genes. Among them, 33 key genes were demonstrated to be associated with OS, and the majority of these were involved in extracellular matrix‑associated and tyrosine kinase receptor signalling pathways. Furthermore, LOC541471 was identified as an OS‑associated lncRNA, and was reported to be involved in the oxidative phosphorylation of GBM with pleckstrin‑2. These findings may significantly enhance our understanding on the aetiology and underlying molecular events of GBM, while the identified candidate genes may serve as novel prognostic markers and potential therapeutic targets for GBM.

Very long intergenic non-coding RNA transcripts and expression profiles are associated to specific childhood acute lymphoblastic leukemia subtypes

Very long intergenic non-coding RNAs (vlincRNAs) are a novel class of long transcripts (~50 kb to 1 Mb) with cell type- or cancer-specific expression. We report the discovery and characterization of 256 vlincRNAs from a cohort of 64 primary childhood pre-B and pre-T acute lymphoblastic leukemia (cALL) samples, of which 61% are novel and specifically expressed in cALL. Validation was performed in 35 pre-B and pre-T cALL primary samples. We show that their expression is cALL immunophenotype and molecular subtype-specific and correlated with epigenetic modifications on their promoters, much like protein-coding genes. While the biological functions of these vlincRNAs are still unknown, our results suggest they could play a role in cALL etiology or progression.

Pathogenesis of ETV6/RUNX1-positive childhood acute lymphoblastic leukemia and mechanisms underlying its relapse

ETV6/RUNX1 (E/R) is the most common fusion gene in childhood acute lymphoblastic leukemia (ALL). Multiple lines of evidence imply a “two-hit” model for the molecular pathogenesis of E/R-positive ALL, whereby E/R rearrangement is followed by a series of secondary mutations that trigger overt leukemia. The cellular framework in which E/R arises and the maintenance of a pre-leukemic condition by E/R are fundamental to the mechanism that underlies leukemogenesis. Accordingly, a variety of studies have focused on the relationship between the clones giving rise to the primary and recurrent E/R-positive ALL. We review here the most recent insights into the pathogenic mechanisms underlying E/R-positive ALL, as well as the molecular abnormalities prevailing at relapse.

NGS-based methylation profiling differentiates TCF3-HLF and TCF3-PBX1 positive B-cell acute lymphoblastic leukemia

Aim: To determine whether methylation differences between mostly fatal TCF3-HLF and curable TCF3-PBX1 pediatric acute lymphoblastic leukemia subtypes can be associated with differential gene expression and remission. Materials & methods: Five (extremely rare) TCF3-HLF versus five (very similar) TCF3-PBX1 patients were sampled before and after remission and analyzed using reduced representation bisulfite sequencing and RNA-sequencing. Results: We identified 7000 differentially methylated CpG sites between subtypes, of which 78% had lower methylation levels in TCF3-HLF. Gene expression was negatively correlated with CpG sites in 23 genes. KBTBD11 clearly differed in methylation and expression between subtypes and before and after remission in TCF3-HLF samples. Conclusion: KBTBD11 hypomethylation may be a promising potential target for further experimental validation especially for the TCF3-HLF subtype.

Epigenetics in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.

GATA2 regulates the erythropoietin receptor in t(12;21) ALL

The t(12;21) (p13;q22) chromosomal translocation resulting in the ETV6/RUNX1 fusion gene is the most frequent structural cytogenetic abnormality in children with acute lymphoblastic leukemia (ALL). The erythropoietin receptor (EPOR), usually associated with erythroid progenitor cells, is highly expressed in ETV6/RUNX1 positive cases compared to other B-lineage ALL subtypes. Gene expression analysis of a microarray database and direct quantitative analysis of patient samples revealed strong correlation between EPOR and GATA2 expression in ALL, and higher expression of GATA2 in t(12;21) patients. The mechanism of EPOR regulation was mainly investigated using two B-ALL cell lines: REH, which harbor and express the ETV6/RUNX1 fusion gene; and NALM-6, which do not. Expression of EPOR was increased in REH cells compared to NALM-6 cells. Moreover, of the six GATA family members only GATA2 was differentially expressed with substantially higher levels present in REH cells. GATA2 was shown to bind to the EPOR 5'-UTR in REH, but did not bind in NALM-6 cells. Overexpression of GATA2 led to an increase in EPOR expression in REH cells only, indicating that GATA2 regulates EPOR but is dependent on the cellular context. Both EPOR and GATA2 are hypomethylated and associated with increased mRNA expression in REH compared to NALM-6 cells. Decitabine treatment effectively reduced methylation of CpG sites in the GATA2 promoter leading to increased GATA2 expression in both cell lines. Although Decitabine also reduced an already low level of methylation of the EPOR in NALM-6 cells there was no increase in EPOR expression. Furthermore, EPOR and GATA2 are regulated post-transcriptionally by miR-362 and miR-650, respectively. Overall our data show that EPOR expression in t(12;21) B-ALL cells, is regulated by GATA2 and is mediated through epigenetic, transcriptional and post-transcriptional mechanisms, contingent upon the genetic subtype of the disease.

The Role of Histone Protein Modifications and Mutations in Histone Modifiers in Pediatric B-Cell Progenitor Acute Lymphoblastic Leukemia

While cancer has been long recognized as a disease of the genome, the importance of epigenetic mechanisms in neoplasia was acknowledged more recently. The most active epigenetic marks are DNA methylation and histone protein modifications and they are involved in basic biological phenomena in every cell. Their role in tumorigenesis is stressed by recent unbiased large-scale studies providing evidence that several epigenetic modifiers are recurrently mutated or frequently dysregulated in multiple cancers. The interest in epigenetic marks is especially due to the fact that they are potentially reversible and thus druggable. In B-cell progenitor acute lymphoblastic leukemia (BCP-ALL) there is a relative paucity of reports on the role of histone protein modifications (acetylation, methylation, phosphorylation) as compared to acute myeloid leukemia, T-cell ALL, or other hematologic cancers, and in this setting chromatin modifications are relatively less well studied and reviewed than DNA methylation. In this paper, we discuss the biomarker associations and evidence for a driver role of dysregulated global and loci-specific histone marks, as well as mutations in epigenetic modifiers in BCP-ALL. Examples of chromatin modifiers recurrently mutated/disrupted in BCP-ALL and associated with disease outcomes include MLL1, CREBBP, NSD2, and SETD2. Altered histone marks and histone modifiers and readers may play a particular role in disease chemoresistance and relapse. We also suggest that epigenetic regulation of B-cell differentiation may have parallel roles in leukemogenesis.

Integrated methylome and transcriptome analysis reveals novel regulatory elements in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common cancer diagnosed in children under the age of 15. In addition to genetic aberrations, epigenetic modifications such as DNA methylation are altered in cancer and impact gene expression. To identify epigenetic alterations in ALL, genome-wide methylation profiles were generated using the methylated CpG island recovery assay followed by next-generation sequencing. More than 25,000 differentially methylated regions (DMR) were observed in ALL patients with ∼90% present within intronic or intergenic regions. To determine the regulatory potential of the DMR, whole-transcriptome analysis was performed and integrated with methylation data. Aberrant promoter methylation was associated with the altered expression of genes involved in transcriptional regulation, apoptosis, and proliferation. Novel enhancer-like sequences were identified within intronic and intergenic DMR. Aberrant methylation in these regions was associated with the altered expression of neighboring genes involved in cell cycle processes, lymphocyte activation and apoptosis. These genes include potential epi-driver genes, such as SYNE1, PTPRS, PAWR, HDAC9, RGCC, MCOLN2, LYN, TRAF3, FLT1, and MELK, which may provide a selective advantage to leukemic cells. In addition, the differential expression of epigenetic modifier genes, pseudogenes, and non-coding RNAs was also observed accentuating the role of erroneous epigenetic gene regulation in ALL.

An evaluation of methods correcting for cell-type heterogeneity in DNA methylation studies

BACKGROUND

Many different methods exist to adjust for variability in cell-type mixture proportions when analyzing DNA methylation studies. Here we present the result of an extensive simulation study, built on cell-separated DNA methylation profiles from Illumina Infinium 450K methylation data, to compare the performance of eight methods including the most commonly used approaches.

RESULTS

We designed a rich multi-layered simulation containing a set of probes with true associations with either binary or continuous phenotypes, confounding by cell type, variability in means and standard deviations for population parameters, additional variability at the level of an individual cell-type-specific sample, and variability in the mixture proportions across samples. Performance varied quite substantially across methods and simulations. In particular, the number of false positives was sometimes unrealistically high, indicating limited ability to discriminate the true signals from those appearing significant through confounding. Methods that filtered probes had consequently poor power. QQ plots of p values across all tested probes showed that adjustments did not always improve the distribution. The same methods were used to examine associations between smoking and methylation data from a case-control study of colorectal cancer, and we also explored the effect of cell-type adjustments on associations between rheumatoid arthritis cases and controls.

CONCLUSIONS

We recommend surrogate variable analysis for cell-type mixture adjustment since performance was stable under all our simulated scenarios.

Population whole-genome bisulfite sequencing across two tissues highlights the environment as the principal source of human methylome variation

BACKGROUND

CpG methylation variation is involved in human trait formation and disease susceptibility. Analyses within populations have been biased towards CpG-dense regions through the application of targeted arrays. We generate whole-genome bisulfite sequencing data for approximately 30 adipose and blood samples from monozygotic and dizygotic twins for the characterization of non-genetic and genetic effects at single-site resolution.

RESULTS

Purely invariable CpGs display a bimodal distribution with enrichment of unmethylated CpGs and depletion of fully methylated CpGs in promoter and enhancer regions. Population-variable CpGs account for approximately 15-20 % of total CpGs per tissue, are enriched in enhancer-associated regions and depleted in promoters, and single nucleotide polymorphisms at CpGs are a frequent confounder of extreme methylation variation. Differential methylation is primarily non-genetic in origin, with non-shared environment accounting for most of the variance. These non-genetic effects are mainly tissue-specific. Tobacco smoking is associated with differential methylation in blood with no evidence of this exposure impacting cell counts. Opposite to non-genetic effects, genetic effects of CpG methylation are shared across tissues and thus limit inter-tissue epigenetic drift. CpH methylation is rare, and shows similar characteristics of variation patterns as CpGs.

CONCLUSIONS

Our study highlights the utility of low pass whole-genome bisulfite sequencing in identifying methylome variation beyond promoter regions, and suggests that targeting the population dynamic methylome of tissues requires assessment of understudied intergenic CpGs distal to gene promoters to reveal the full extent of inter-individual variation.

Developmental genes significantly afflicted by aberrant promoter methylation and somatic mutation predict overall survival of late-stage colorectal cancer

Carcinogenesis is an exceedingly complicated process, which involves multi-level dysregulations, including genomics (majorly caused by somatic mutation and copy number variation), DNA methylomics, and transcriptomics. Therefore, only looking into one molecular level of cancer is not sufficient to uncover the intricate underlying mechanisms. With the abundant resources of public available data in the Cancer Genome Atlas (TCGA) database, an integrative strategy was conducted to systematically analyze the aberrant patterns of colorectal cancer on the basis of DNA copy number, promoter methylation, somatic mutation and gene expression. In this study, paired samples in each genomic level were retrieved to identify differentially expressed genes with corresponding genetic or epigenetic dysregulations. Notably, the result of gene ontology enrichment analysis indicated that the differentially expressed genes with corresponding aberrant promoter methylation or somatic mutation were both functionally concentrated upon developmental process, suggesting the intimate association between development and carcinogenesis. Thus, by means of random walk with restart, 37 significant development-related genes were retrieved from a priori-knowledge based biological network. In five independent microarray datasets, Kaplan-Meier survival and Cox regression analyses both confirmed that the expression of these genes was significantly associated with overall survival of Stage III/IV colorectal cancer patients.

Preserved global histone H4 acetylation linked to ETV6-RUNX1 fusion and PAX5 deletions is associated with favorable outcome in pediatric B-cell progenitor acute lymphoblastic leukemia

Epigenetic dysregulation is a hallmark of cancer executed by a number of complex processes the most important of which converge on DNA methylation and histone protein modifications. Epigenetic marks are potentially reversible and thus promising drug targets. In the setting of acute lymphoblastic leukemia (ALL) they have been associated with clinicopathological features including risk of relapse or molecular subgroups of the disease. Here, using immunocytochemistry of bone marrow smears from diagnosis, we studied global histone H4 acetylation, whose loss was previously linked to treatment failure in adults with ALL, in pediatric patients. We demonstrate that preserved global histone H4 acetylation is significantly associated with favorable outcome (RFS, EFS, OS) in children with B cell progenitor (BCP) ALL, recapitulating the findings from adult populations. Further, for the first time we demonstrate differential histone H4 acetylation in molecular subclasses of BCP-ALL including cases with ETV6-RUNX1 fusion gene or PAX5 deletion or deletions in genes linked to B cell development. We conclude global histone H4 acetylation is a prognostic marker and a potential therapeutic target in ALL.

Parsimonious covariate selection for a multicategory ordered response

We propose a flexible continuation ratio (CR) model for an ordinal categorical response with potentially ultrahigh dimensional data that characterizes the unique covariate effects at each response level. The CR model is the logit of the conditional discrete hazard function for each response level given covariates. We propose two modeling strategies, one that keeps the same covariate set for each hazard function but allows regression coefficients to arbitrarily change with response level, and one that allows both the set of covariates and their regression coefficients to arbitrarily change with response. Evaluating a covariate set is accomplished by using the nonparametric bootstrap to estimate prediction error and their robust standard errors that do not rely on proper model specification. To help with interpretation of the selected covariate set, we flexibly estimate the conditional cumulative distribution function given the covariates using the separate hazard function models. The goodness-of-fit of our flexible CR model is assessed with graphical and numerical methods based on the cumulative sum of residuals. Simulation results indicate the methods perform well in finite samples. An application to B-cell acute lymphocytic leukemia data is provided.

Epigenome data release: a participant-centered approach to privacy protection

Large-scale epigenome mapping by the NIH Roadmap Epigenomics Project, the ENCODE Consortium and the International Human Epigenome Consortium (IHEC) produces genome-wide DNA methylation data at one base-pair resolution. We examine how such data can be made open-access while balancing appropriate interpretation and genomic privacy. We propose guidelines for data release that both reduce ambiguity in the interpretation of open-access data and limit immediate access to genetic variation data that are made available through controlled access.

Characterization of CpG sites that escape methylation on the inactive human X-chromosome

In many whole genome studies of gene expression or modified cytosines, data from probes localized to the X-chromosome are removed from analyses due to gender bias. Previously, we observed population differences in cytosine modifications between Caucasian and African lymphoblastoid cell lines (LCLs) on the autosomes using whole genome arrays to measure modified cytosines. DNA methylation plays a critical role in establishment and maintenance of X-chromosome inactivation in females. Therefore, we reasoned that by investigating cytosine modification patterns specifically on the X-chromosome, we could obtain valuable information about a chromosome that is often disregarded in genome-wide analyses. We investigated population differences in cytosine modification patterns along the X-chromosome between Caucasian and African LCLs and identified novel sites that escape methylation on the inactive X-chromosome (Xi) in females. We characterized the chromatin state of these loci by incorporating the extensive histone modification ChIP-seq data generated by ENCODE. To explore the relationship between DNA and histone modifications further, we hypothesized that BRD4, a protein that binds acetylated histones, could be preventing some sites from becoming de novo methylated. To test this, we treated 4 female LCLs with JQ1, a small molecule inhibitor of BRD4, but found that JQ1 treatment induced minor changes in cytosine modification levels, and the majority of sites escaping methylation on the Xi remained unmethylated. This suggests that other epigenetic mechanisms or transcription factors are likely playing a larger role in protecting these sites from de novo methylation on the Xi.

Whole-genome fingerprint of the DNA methylome during human B cell differentiation

We analyzed the DNA methylome of ten subpopulations spanning the entire B cell differentiation program by whole-genome bisulfite sequencing and high-density microarrays. We observed that non-CpG methylation disappeared upon B cell commitment, whereas CpG methylation changed extensively during B cell maturation, showing an accumulative pattern and affecting around 30% of all measured CpG sites. Early differentiation stages mainly displayed enhancer demethylation, which was associated with upregulation of key B cell transcription factors and affected multiple genes involved in B cell biology. Late differentiation stages, in contrast, showed extensive demethylation of heterochromatin and methylation gain at Polycomb-repressed areas, and genes with apparent functional impact in B cells were not affected. This signature, which has previously been linked to aging and cancer, was particularly widespread in mature cells with an extended lifespan. Comparing B cell neoplasms with their normal counterparts, we determined that they frequently acquire methylation changes in regions already undergoing dynamic methylation during normal B cell differentiation.

Emerging technologies for studying DNA methylation for the molecular diagnosis of cancer

DNA methylation is an epigenetic mechanism that plays a key role in regulating gene expression and other functions. Although this modification is seen in different sequence contexts, the most frequently detected DNA methylation in mammals involves cytosine-guanine dinucleotides. Pathological alterations in DNA methylation patterns are described in a variety of human diseases, including cancer. Unlike genetic changes, DNA methylation is heavily influenced by subtle modifications in the cellular microenvironment. In all cancers, aberrant DNA methylation is involved in the alteration of a large number of oncological pathways with relevant theranostic utility. Several technologies for DNA methylation mapping have been developed recently and successfully applied in cancer studies. The scope of these technologies varies from assessing a single cytosine-guanine locus to genome-wide distribution of DNA methylation. Here, we review the strengths and weaknesses of these approaches in the context of clinical utility for the molecular diagnosis of human cancers.

DNA methylation-based subtype prediction for pediatric acute lymphoblastic leukemia

BACKGROUND

We present a method that utilizes DNA methylation profiling for prediction of the cytogenetic subtypes of acute lymphoblastic leukemia (ALL) cells from pediatric ALL patients. The primary aim of our study was to improve risk stratification of ALL patients into treatment groups using DNA methylation as a complement to current diagnostic methods. A secondary aim was to gain insight into the functional role of DNA methylation in ALL.

RESULTS

We used the methylation status of ~450,000 CpG sites in 546 well-characterized patients with T-ALL or seven recurrent B-cell precursor ALL subtypes to design and validate sensitive and accurate DNA methylation classifiers. After repeated cross-validation, a final classifier was derived that consisted of only 246 CpG sites. The mean sensitivity and specificity of the classifier across the known subtypes was 0.90 and 0.99, respectively. We then used DNA methylation classification to screen for subtype membership of 210 patients with undefined karyotype (normal or no result) or non-recurrent cytogenetic aberrations ('other' subtype). Nearly half (n = 106) of the patients lacking cytogenetic subgrouping displayed highly similar methylation profiles as the patients in the known recurrent groups. We verified the subtype of 20% of the newly classified patients by examination of diagnostic karyotypes, array-based copy number analysis, and detection of fusion genes by quantitative polymerase chain reaction (PCR) and RNA-sequencing (RNA-seq). Using RNA-seq data from ALL patients where cytogenetic subtype and DNA methylation classification did not agree, we discovered several novel fusion genes involving ETV6, RUNX1, and PAX5.

CONCLUSIONS

Our findings indicate that DNA methylation profiling contributes to the clarification of the heterogeneity in cytogenetically undefined ALL patient groups and could be implemented as a complementary method for diagnosis of ALL. The results of our study provide clues to the origin and development of leukemic transformation. The methylation status of the CpG sites constituting the classifiers also highlight relevant biological characteristics in otherwise unclassified ALL patients.

Whole-epigenome analysis in multiple myeloma reveals DNA hypermethylation of B cell-specific enhancers

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

Epigenetic landscape correlates with genetic subtype but does not predict outcome in childhood acute lymphoblastic leukemia

Although children with acute lymphoblastic leukemia (ALL) generally have a good outcome, some patients do relapse and survival following relapse is poor. Altered DNA methylation is highly prevalent in ALL and raises the possibility that DNA methylation-based biomarkers could predict patient outcome. In this study, genome-wide methylation analysis, using the Illumina Infinium HumanMethylation450 BeadChip platform, was carried out on 52 diagnostic patient samples from 4 genetic subtypes [ETV6-RUNX1, high hyperdiploidy (HeH), TCF3-PBX1 and dic(9;20)(p11-13;q11)] in a 1:1 case-control design with patients who went on to relapse (as cases) and patients achieving long-term remission (as controls). Pyrosequencing assays for selected loci were used to confirm the array-generated data. Non-negative matrix factorization consensus clustering readily clustered samples according to genetic subgroups and gene enrichment pathway analysis suggested that this is in part driven by epigenetic disruption of subtype specific signaling pathways. Multiple bioinformatics approaches (including bump hunting and individual locus analysis) were used to identify CpG sites or regions associated with outcome. However, no associations with relapse were identified. Our data revealed that ETV6-RUNX1 and dic(9;20) subtypes were mostly associated with hypermethylation; conversely, TCF3-PBX1 and HeH were associated with hypomethylation. We observed significant enrichment of the neuroactive ligand-receptor interaction pathway in TCF3-PBX1 as well as an enrichment of genes involved in immunity and infection pathways in ETV6-RUNX1 subtype. Taken together, our results suggest that altered DNA methylation may have differential impacts in distinct ALL genetic subtypes.

Differences in DNA methylation signatures reveal multiple pathways of progression from adenoma to colorectal cancer

BACKGROUND & AIMS

Genetic and epigenetic alterations contribute to the pathogenesis of colorectal cancer (CRC). There is considerable molecular heterogeneity among colorectal tumors, which appears to arise as polyps progress to cancer. This heterogeneity results in different pathways to tumorigenesis. Although epigenetic and genetic alterations have been detected in conventional tubular adenomas, little is known about how these affect progression to CRC. We compared methylomes of normal colon mucosa, tubular adenomas, and colorectal cancers to determine how epigenetic alterations might contribute to cancer formation.

METHODS

We conducted genome-wide array-based studies and comprehensive data analyses of aberrantly methylated loci in 41 normal colon tissue, 42 colon adenomas, and 64 cancers using HumanMethylation450 arrays.

RESULTS

We found genome-wide alterations in DNA methylation in the nontumor colon mucosa and cancers. Three classes of cancers and 2 classes of adenomas were identified based on their DNA methylation patterns. The adenomas separated into classes of high-frequency methylation and low-frequency methylation. Within the high-frequency methylation adenoma class a subset of adenomas had mutant KRAS. Additionally, the high-frequency methylation adenoma class had DNA methylation signatures similar to those of cancers with low or intermediate levels of methylation, and the low-frequency methylation adenoma class had methylation signatures similar to that of nontumor colon tissue. The CpG sites that were differentially methylated in these signatures are located in intragenic and intergenic regions.

CONCLUSIONS

Genome-wide alterations in DNA methylation occur during early stages of progression of tubular adenomas to cancer. These findings reveal heterogeneity in the pathogenesis of colorectal cancer, even at the adenoma step of the process.

Epigenetics and chemoresistance in childhood acute lymphoblastic leukemia

SUMMARY For children with acute lymphoblastic leukemia (ALL) who relapse, prognosis is poor and novel therapeutic strategies are needed. In the last decade, it has become apparent that ALL exhibits unique epigenetic patterns in addition to the well known cytogenetic findings. Furthermore, whole genome sequencing efforts are revealing recurrent mutations in epigenetic modifiers in ALL. Aberrant epigenetic modulation may be involved in leukemic transformation and resistance to chemotherapy. Consequently, compounds that specifically modulate the maintenance of such epigenetic programs may offer new approaches to therapy, including the modulation or prevention of chemoresistance in ALL. In this article, we review some of the most recent findings with regard to epigenetic aberrations in ALL, and discuss therapeutic strategies that are currently in development.

Epigenetic deregulation in pediatric acute lymphoblastic leukemia

Similar to most cancers, genome-wide DNA methylation profiles are commonly altered in pediatric acute lymphoblastic leukemia (ALL); however, recent observations highlight that a large portion of malignancy-associated DNA methylation alterations are not accompanied by related gene expression changes. By analyzing and integrating the methylome and transcriptome profiles of pediatric B-cell ALL cases and primary tissue controls, we report 325 genes hypermethylated and downregulated and 45 genes hypomethylated and upregulated in pediatric B-cell ALL, irrespective of subtype. Repressed cation channel subunits and cAMP signaling activators and transducers are overrepresented, potentially indicating a reduced cellular potential to receive and propagate apoptotic signals. Furthermore, we report specific DNA methylation alterations with concurrent gene expression changes within individual ALL subtypes. The ETV6-RUNX1 translocation was associated with downregulation of ASNS and upregulation of the EPO-receptor, while Hyperdiploid patients (> 50 chr) displayed upregulation of B-cell lymphoma (BCL) members and repression of PTPRG and FHIT. In combination, these data indicate genetically distinct B-cell ALL subtypes contain cooperative epimutations and genome-wide epigenetic deregulation is common across all B-cell ALL subtypes.

Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia

BACKGROUND

Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood.

RESULTS

We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status.

CONCLUSIONS

Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.