While mutations of splicing and epigenetic factors have been reported in adult AML and related myeloid disorders, relatively few such changes have been identified in pediatric AML. We previously identified a chromatin remodeling helicase, PASG (SMARCA6, HELLS, LSH), by down-regulation in AML cell lines following cytokine withdrawal and identified an alternatively spliced variant lacking a highly conserved (STRAGGLG) domain. To assess the prevalence of this splicing variant (PASGΔ75) in pediatric AML, we tested 167 diagnostic specimens from the TARGET-AML cohort for fractional PASGΔ75 expression (PASGΔ75/PASG) using a discriminatory RT-qPCR assay. These studies demonstrated a broad, continuous distribution of PASGΔ75 with right skew (mean PASGΔ75: 26%, interquartile range: 9% – 41%) that was not significantly associated with cytogenetic class (inv16, t(8;21), MLL, normal) or FAB subtype. For further comparison, specimens were quantized by PASGΔ75 quartile.Given reported associations between loss of PASG function and abnormalities of genomic methylation, we tested 48 AML specimens at the extremes of the PASGΔ75 distribution for total 5-methylcytosine (5-mC) content by liquid chromatography/tandem mass spectrometry. The mean total methylation was significantly lower in the high compared to the low PASGΔ75 groups (mean 5-mC 3.95% vs. 4.22% of cytosine, p=0.015 Mann-Whitney). To identify specific regions of altered methylation, we used high-throughput sequencing of DNA enriched by pull-down with the methyl binding domain fragment of MBD2 (MBD-Seq). Comparison of summed methylation signal across regions flanking RefSeq transcriptional start sites (TSS) showed the expected decrease in methylation just upstream of the TSS in both groups. However, methylation more distal to the TSS was proportionally lower in PASGΔ75 high than PASGΔ75 low samples (Fig 1a). To evaluate methylation at CpG islands (CGI), UCSC CGI were scaled to 500 bp and MBD-Seq data were summed across 20Kb flanking CGI. While both groups showed the anticipated increase in methylation signal on CGI, methylation in the shore regions immediately flanking CGI was proportionally lower in high PASGΔ75 compared to low PASGΔ75 samples (Fig 1b), further suggesting epigenetic differences between these sample groups. [Display omitted] Because we were unable to identify sequence variants in PASG intron 18 or flanking exons that might explain alternative splicing, we asked whether expression of PASGΔ75 was associated with global changes in transcript splicing. We evaluated gene expression patterns on the Affymetrix HuGene array, with assessment of alternative splicing using Partek software alternative splicing (altsplice) algorithm for quartile-grouped samples. In contrast to comparison of adjacent quartile groups, which showed modest changes in expression and relatively few transcripts with significant altsplice scores, comparison of the highest and lowest quartile samples showed marked changes in gene expression and a large number of alternatively spiced transcripts as assessed by significance of the altsplice score (Fig 2). In addition to splicing changes, these analyses suggested marked differences in gene expression patterns of AML specimens grouped by PASGΔ75 quartile, with clear separation of Q1 and Q4 samples by principal components analysis. Using a conservative Wilcoxon gene sets test and limiting ourselves to small, curated Biocarta pathways, we found expression patterns associated with high deletion variant expression strongly linked to overlapping pathways involving DNA repair, replication, and cell cycle progression. [Display omitted] These data suggest the existence of a previously unrecognized AML subclass characterized by widespread and coordinated changes in RNA expression, alternative transcript splicing, and epigenetic modifications. Disclosures:No relevant conflicts of interest to declare.TablePathways (Biocarta)Pathway ClassBenjamini Hochberg Corrected p-ValueATR/BRCA1/BRCA2DNA Damage Response1.3E-6RB/DNA DamageDNA Damage Response3.5E-3p27 PhosphorylationCell Cycle Progression3.5E-3G2/M CheckpointCell Cycle Progression3.5E-3G1/S CheckpointCell Cycle Progression4.1E-3PLK3Cell Cycle Progression4.1E-3MEF2DApoptosis0.01SRC/PTPaCell Cycle Progression0.02Mitochondrial Acetyl-Co ShuttleMetabolism0.02p53 SignalingDNA Damage Response0.04E2F-1Cell Cycle Progression0.04ATM SignalingDNA Damage Response0.04
Read full abstract