ALBERT

Description: Introduction: Chronic myelomonocytic leukemia (CMML) is a clonal hematopoietic neoplasm characterized by sustained peripheral blood (PB) monocytosis and an inherent risk for leukemic transformation. Clonal origins of the disease can be detected in hematopoietic progenitor cells (CD34+/CD38-), while the complete spectrum of mutational evolution can be seen in circulating monocytes (CD14+). Cell sorting strategies have been employed to select cells in CMML, and while there are adequate monocyte numbers in the PB, there are very few circulating progenitor cells. In addition, attrition related to the selection process significantly depletes primary cells available for biological experiments and multiomics studies such as RNA-seq, ChIP-seq, ATAC-seq, and DIP-seq. While single-cell methods may be able to overcome this challenge, bulk sequencing methods remain a robust and cost-effective approach. We hypothesized that, secondary to the stem cell origin of this disease and significant myeloproliferation, PB mononuclear cells (MNC) would provide comparable results with regards to transcriptomic analysis, in comparison to cell selection procedures. Methods: Peripheral blood obtained from 15 molecularly annotated patients with WHO-defined CMML was ACK-lysed and subjected to a Ficoll procedure for collection of MNC. MNC were left unsorted (n=5) or further selected for CD34+/CD38- (n=5) and CD14+ (n=5) using a fully automated RoboSep-S (StemCell Technologies) protocol. All samples were then subjected to bulk whole transcriptome shotgun sequencing (using Illumina TruSeq and an Illumina HiSeq 4000). After data quality control, counts of detectable transcripts were log2-normalized and Pearson's product-moment correlation coefficients were calculated to evaluate the correlation between the two cell-sorting strategies and unsorted cells in terms of detectable transcripts. To visualize sample differences log2-normalized transcripts counts were centered and scaled per gene for a select number of genes relevant to myeloid biology as well as a number of housekeeping genes. Results: Fifteen patients with WHO-defined CMML, median age 69 years (55-73 years), 66% male, were included. Next generation sequencing for somatic mutations was performed on PB MNC obtained at CMML diagnosis (Figure 1, top heatmap). Considering the small sample size, mutations were evenly distributed among groups with the exception of ASXL1 (higher frequency in CD14+ and CD34+/CD38- cells), ZRSR2 (higher frequency in unsorted cells), and TET2 (lower frequency in CD14+ cells). The three groups were also well matched with regards to other CMML-related variables such as WHO and FAB morphological subtypes, cytogenetic abnormalities, and risk stratification by the Mayo Molecular Model. Transcriptomic analysis revealed a strong positive correlation between the median number of log2-normalized detectable transcripts in unsorted cells and CD34+/CD38- cells (ρ = 0.96, p 〈 0.001, top scatterplot). Likewise, there was a strong positive correlation between the median number of log2-normalized detectable transcripts in unsorted cells and CD14+ cells (ρ = 0.91, p 〈 0.001, bottom scatterplot). The latter correlation was marginally lower, which was explained by increased global gene expression in 3 of the 5 CD14+ samples (bottom heatmap). Increased gene expression in these 3 samples involved key myeloid genes and housekeeping genes known to have stable expression across human tissues alike. In comparison to PB MNC, both cell sorting strategies resulted in significant depletion of primary cells required for other experiments, and for procedures such as ChIP-seq, DIP-seq and ATAC-seq (CD34+/CD38- had greater depletion than CD14+). Additional experiments to assess this strategy for the above mentioned epigenetic studies are currently being planned. Conclusions: Accounting for sample differences, different cell sorting strategies (unsorted, CD34+/CD38- selection, and CD14+ selection) yielded similar results when performing bulk transcriptomic assessments on PB MNC from patients with CMML. For the purpose of gene expression profiling there was no clear advantage with CD34+/CD38- or CD14+ selection. These results support the use of unsorted cells for bulk transcriptomic analysis in CMML. Figure 1 Disclosures Patnaik: Stem Line Pharmaceuticals.: Membership on an entity's Board of Directors or advisory committees.

Description: Introduction: Chronic myelomonocytic leukemia (CMML), an aggressive myeloid malignancy, can be categorized into two subtypes, proliferative CMML (pCMML) and dysplastic (dCMML), based on a white blood cell (WBC) count of ≥ 13 x 109/L for the former (Arber et al. Blood 2016). While this WBC cut off is somewhat arbitrary, patients with pCMML have unique phenotypic features and a shorter survival. We carried out this study to assess the genomic, transcriptomic and epigenetic landscapes of these two CMML subtypes. Methods: Peripheral blood (PB) and bone marrow (BM) mononuclear cells (MNC) were obtained from WHO-defined CMML patients. Next generation sequencing (NGS) using a 36-gene panel was performed on 350 patients with Illumina HiSeq4000 platform with median read depth of 400X. RNA sequencing (RNA-seq) was performed on 25 patients by bulk whole transcriptome sequencing using Illumina TruSeq. DNA immunoprecipitation sequencing (DIP-seq) was performed on 18 patients using 5-methylcystocine (5mC), 5-hydroxymethylcytosine (5hmC) and bridging monoclonal antibodies with subsequent paired-end sequencing using HiSeq4000. Chromatin immunoprecipitation sequencing (ChIP-seq) was performed on 30 patients with Illumina HiSeq2500 to a depth of 25 million for histone 3 lysine 4 monomethylation (H3K4me1) and histone 3 lysine 4 trimethylation (H3K4me3) and 50 million reads for histone 3 lysine 27 trimethylation (H3K27me3) and Input per sample. Results: Five hundred and seventy-three patients with WHO defined CMML were included; median age 71 years (range 18-95 years), 67% males. 282 patients had pCMML (49%), while 291 (51%) had dCMML. As pre-defined, patients with pCMML were more likely to have higher absolute monocyte counts (p

Description: Introduction: Truncating mutations in the Additional Sex Combs-Like 1 (ASXL1) gene are associated with a proliferative disease phenotype and poor survival outcomes across the spectrum of myeloid malignancies including chronic myelomonocytic leukemia (CMML). ASXL1 is thought to act as a chromatin modifier regulating transcriptional activity, however the exact mechanisms and resulting chromatin states remain controversial. We interrogated the epigenome of 16 patients with ASXL1-mutant and -wildtype CMML using a multiomics approach. Methods: Bone marrow mononuclear cells from patients with CMML (8 ASXL1-mutant, 8 -wildtype) were subjected to targeted NGS of DNA, whole transcriptome shotgun sequencing (RNA-seq), immunoprecipitation (IP) of DNA (hydroxy-)methyl residues (DIP-seq), IP of the histone modifications H3K4me1, H3K4me3, and H3K27me3 (ChIP-seq), and DNA transposase accessibility studies (ATAC-seq). After quality control all samples were sequenced on an Illumina HiSeq 4000 before further processing and data analysis. Global assessments of DNA (hydroxy-)methylation, DNA accessibility, and histone modifications between ASXL1-mutant and -wildtype CMML were performed. Differential gene expression was performed to define the up-regulated genes in ASXL1-mutant disease. The promoter regions of these up-regulated genes (defined as transcription start site ±3kb) were compared using the aforementioned multiomics approach. Epigenomic modification of the promoter region facilitating up-regulation of transcription was defined as the presence of a signal peak in ASXL1-mutant disease (in the absence of a signal peak in -wildtype disease) or 25% increase in a common signal peak (H3K4me1/3, 5hmC), the presence of a unique signal peak in ASXL1-mutant disease (ATAC), or the absence of a signal peak in ASXL1-mutant disease (in the presence of a signal peak in -wildtype disease), or 25% decrease in a common signal peak (H3K27me3, 5mC). Results: Sixteen patients with CMML, median age 69 years (48 - 77), 63% male, were included. Half of the patients had proliferative disease (pCMML) and half of them had truncating frameshift mutations in ASXL1 (heatmap). All ASXL1 variant allele frequencies were compatible with heterozygosity (31 - 48%). The burden of co-mutations was similar between ASXL1-wildtype and ASXL1-mutant disease (21 versus 23 per group; no difference in the median number of co-mutations, p = 0.684). The spectrum of co-mutations was typical for CMML, involving spliceosome components, epigenetic regulators, chromatin regulators, and cell signaling molecules (heatmap). There was a predominant up-regulation of gene expression in ASXL1-mutant patients: 707 genes up- and 124 down-regulated (volcano plot, FDR 〈 0.05 for all genes). Functional annotation of the up-regulated genes showed cell division, mitotic nuclear division, sister chromatid cohesion, DNA replication, and G1/S transition to be the 5 most enriched processes (accounting for 29% of all up-regulated genes, FDR 〈 1x10-10 for all terms). The up-regulated genes included several potential therapeutic targets and HOXA family members (including HOXA9). There were global increases in H3K4me1/3, 5mC, and 5hmC, decreases in H3K27me3, as well as a more relaxed chromatin conformation (bar graphs). Many of these epigenomic changes affected non-coding regions. When focusing on the promoter regions of the 707 up-regulated genes there was evidence of one or more of the interrogated epigenomic mechanisms facilitating transcription for 519 of the genes (73%). The most abundant mechanism was histone modification, followed by changes in DNA (hydroxy-)methylation, and increased chromatin accessibility, with considerable overlap (Venn diagram). For HOXA9, a known driver of leukemogenesis, the data supported a loss of H3K27me3 as the most prominent among the interrogated epigenomic regulatory mechanisms (average signal tracks). Conclusions: The transcriptome and chromatin conformation of ASXL1-mutant CMML are skewed towards proliferation and mirror the aggressive disease phenotype observed in practice. There is evidence of histone modification as well as changes in DNA methylation, and chromatin conformation facilitating transcriptional activity including known leukemogenic drivers. Additional regulatory mechanisms such as gene body methylation and enhancer elements require further exploration. Figure Disclosures Patnaik: Stem Line Pharmaceuticals.: Membership on an entity's Board of Directors or advisory committees.

Description: Introduction: Additional Sex Combs-Like 1 (ASXL1) is a chromatin modifier frequently affected by truncating mutations in myeloid malignancies. These mutations are associated with poor survival outcomes and increased rates of acute leukemic transformation. In chronic myelomonocytic leukemia (CMML), ASXL1 mutations are thought to affect transcriptional activity mainly by modifying histone marks, however additional epigenomic mechanisms have not been fully explored. We interrogated the epigenome of patients with ASXL1-mutant (MT) and -wildtype (WT) CMML using a multiomics approach to define cis-regulatory elements (CREs) such as distal enhancers (DEs). Methods: Bone marrow mononuclear cells from patients with CMML were subjected to targeted NGS of DNA, whole transcriptome shotgun sequencing (RNA-seq), immunoprecipitation (IP) of DNA (hydroxy-)methyl residues (DIP-seq), IP of the histone modifications H3K4me1, H3K4me3, and H3K27me3 (ChIP-seq), and DNA transposase accessibility studies (ATAC-seq). After quality control all samples were sequenced on an Illumina HiSeq 4000 before further processing and data analysis. Global assessments of DNA (hydroxy-)methylation, DNA accessibility, and histone modifications between ASXL1 MT and WT CMML were performed. The samples in the two groups were treated as biological replicates and subjected to a consensus peak calling strategy requiring an overlap of at least 30% between samples and an adjusted p-value 〈 5x10-5 for a signal peak to be considered statistically significant. Differential gene expression was estimated to define the up-regulated genes in ASXL1 MT CMML. Potential CREs were defined as sites with statistically significant signal peaks overlapping in at least two of the three epigenomic marks: H3K4me1, 5hmC, and ATAC. Potential DEs were defined as CREs in non-coding regions outside promoter regions (defined as transcription start site ±3kb) that were annotated in GeneHancer. Annotated DEs only present in ASXL1 MT but not WT CMML (specific DEs) were intersected with the list of up-regulated genes and the ReMap atlas. Results: Sixteen WHO-defined CMML patients were included, median age 69 years (48 - 77), 63% male; of which 8 patients (50%, all truncating frame shift mutations) were ASXL1 MT and 8 (50%) WT. The burden and spectrum of co-mutations was similar between ASXL1 WT and MT CMML (21 versus 23 per group; p = 0.684; heatmap). There was a predominant up-regulation of gene expression in ASXL1 MT CMML: 707 genes up- and 124 down-regulated (volcano plot, FDR 〈 0.050 for all genes). There were 64336 potential CREs, the vast majority (97%) being present in both ASXL1 MT and WT CMML (left Venn diagram). These CREs were most commonly located in introns, promoter regions, and distal non-coding regions (bar graph and pie chart). There were 1303 CREs unique to ASXL1 MT CMML (specific DEs), 1161 (90%) of which were annotated in GeneHancer (left Euler diagram). Of these 1161 annotated specific DEs 859 (74%) were located outside promoter regions and 34 (4%) of them were known to be associated with genes up-regulated in ASXL1 MT CMML (Euler diagrams). These specific DEs were characterized by an increase in H3K4me1 occupancy and DNA accessibility (average signal tracks, purple bars indicating annotated DEs, thin bars below peaks indicating statistical significance). We previously observed epigenomic modification of promoter regions in 519 of the 707 up-regulated genes (73%) facilitating transcriptional activity in ASXL1 MT CMML. For 13 of the up-regulated genes (right Venn diagram, blue genes in volcano plot) the specific DEs were the sole identified mechanism, while for the other 21 genes there were additional mechanisms noted in the promoter region. The top five transcription factor candidates binding the 34 specific DEs included JMJD1C, MYC, KDM5B, RCOR1, and HDAC2 (-log10(E) 〉 40 for all candidates). Conclusions: Using a multiomics approach based on H3K4me1, 5hmC, and ATAC data we identified potential CREs in ASXL1 MT CMML and characterized potential DEs using publicly available annotation data. Specific DEs were associated with up-regulated genes serving as a possible explanation for the observed transcriptional activity, shedding further light on the adverse prognostic impact associated with ASXL1 mutations. Figure 1 Disclosures Patnaik: Stem Line Pharmaceuticals.: Membership on an entity's Board of Directors or advisory committees.