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.

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.

A complete understanding of the mechanisms involved in the development of pre-B ALL is lacking. In this study we integrated DNA methylation data and gene expression data to elucidate the impact of aberrant intergenic DNA methylation on gene expression in pre-B ALL. We found a subset of differentially methylated intergenic loci that were associated with altered gene expression in pre-B ALL patients. Notably, 84% of these regions were also bound by transcription factors known to play roles in differentiation and B-cell development in a lymphoblastoid cell line. Further, an overall down-regulation of eRNA transcripts was observed in pre-B ALL patients and these transcripts were associated with the down-regulation of putative target genes involved in B-cell migration, proliferation and apoptosis. The identification of novel putative regulatory regions highlights the significance of intergenic DNA sequences and may contribute to the identification of new therapeutic targets for the treatment of pre-B ALL.

The potential effect of trans regulation by AS-lncRNA through DNA/RNA binding is significant as sequence alignment analysis of the 25 most differentially methylated AS-lncRNAs revealed 368 genes containing highly similar sequences with a median nucleotide identity of 90.8% and binding span of 122 base pairs. Regulation of biological processes and anatomical structure development were over represented. ALL classification schemes based on AS-lncRNA methylation can provide new insights into its pathogenesis and treatment.