• CMML monocytes exhibit a proinflammatory transcriptional signature, contributing to malignant expansion and increased cardiovascular risk.

Chronic myelomonocytic leukemia (CMML) is an aggressive myeloid neoplasm of older individuals characterized by persistent monocytosis. Somatic mutations in CMML are heterogeneous and only partially explain the variability in clinical outcomes. Recent data suggest that cardiovascular morbidity is increased in CMML and contributes to reduced survival. Clonal hematopoiesis of indeterminate potential (CHIP), the presence of mutated blood cells in hematologically normal individuals, is a precursor of age-related myeloid neoplasms and associated with increased cardiovascular risk. To isolate CMML-specific alterations from those related to aging, we performed RNA sequencing and DNA methylation profiling on purified monocytes from CMML patients and from age-matched (old) and young healthy controls. We found that the transcriptional signature of CMML monocytes is highly proinflammatory, with upregulation of multiple inflammatory pathways, including tumor necrosis factor and interleukin (IL)-6 and -17 signaling, whereas age per se does not significantly contribute to this pattern. We observed no consistent correlations between aberrant gene expression and CpG island methylation, suggesting that proinflammatory signaling in CMML monocytes is governed by multiple and complex regulatory mechanisms. We propose that proinflammatory monocytes contribute to cardiovascular morbidity in CMML patients and promote progression by selection of mutated cell clones. Our data raise questions of whether asymptomatic patients with CMML benefit from monocyte-depleting or anti-inflammatory therapies.

Chronic myelomonocytic leukemia (CMML) is a genetically heterogeneous hematopoietic stem cell disorder that combines features of a myelodysplastic syndrome (MDS) and a myeloproliferative neoplasm (MPN) and occurs almost exclusively in older individuals.1,-3  The epigenetic regulators TET2 (40%-60% of patients) and ASXL1 (30%-60% of patients)4  and the splicing factor SRSF2 (30%-50%)5,6  are the most frequently mutated genes in CMML. Mutational activation of the RAS signaling pathway is also common,7  and ∼25% of CMML patients have structural or numerical chromosomal aberrations.8  Multivariate analyses suggest that somatic mutations explain only one quarter of the heterogeneity of CMML outcomes, inferring that leukemia-specific factors other than somatic mutations and/or host factors have a major impact on prognosis.9  For instance, aberrant DNA methylation correlates with outcome and response to hypomethylating agents, irrespective of somatic mutations.10,11 

Aging causes profound changes in the hematopoietic system. Immune cell subsets and cytokine profiles exhibit patterns of enhanced inflammation, and there is a bias toward polyclonal myelomonocytic differentiation.12,,,-16  At the same time, detection of somatic mutations in blood cells of hematologically normal individuals, variably referred to as clonal hematopoiesis of indeterminate potential (CHIP) or age-related clonal hematopoiesis (ARCH), increases steeply with age.17,,,,,,-24  An increasing body of evidence supports an intimate interplay between a systemic inflammatory state and myeloid neoplasms. The risk of myeloid malignancies is significantly elevated in patients with a history of infection or autoimmune disease,25  and a recent report demonstrated that bacterial infection can promote the expansion of Tet2-deficient hematopoietic cells in mice.25,26  The relationship between inflammation and somatic mutations may be particularly important in the context of monocyte-rich myeloid neoplasms, as a recent study found that CMML patients have a more than twofold increase in the rate of cardiovascular events compared with patients with chronic lymphocytic leukemia.27  Given the extreme skewing of the CMML age distribution (median age of onset, 70 years), we hypothesized that separating age-related from CMML-specific alterations will be critical to mechanistically understand CMML oncogenesis and associated CMML-specific risks.28  As a first step in this direction, we performed comparative transcriptional and DNA methylation analyses of CD14+ monocytes from CMML patients, age-matched (old) healthy controls, and young healthy controls.

Study population

CD14+ monocytes were isolated from the blood of untreated CMML patients (n = 12; median age 72 years; range, 57-87), age-matched healthy controls (old controls, n = 12; median age, 68 years; range, 62-74) and young healthy controls (young controls, n = 16; median age, 29 years; range, 24-44) and subjected to RNA sequencing. In addition, we performed DNA methylation profiling for the CMML patients and old controls. The University of Utah Institutional Review Board approved the study protocols, and all participants provided written informed consent before sample collection.

Sample processing, monocyte isolation, and DNA and RNA extraction

Standard methods were used. For details see the supplemental Information.

RNA sequencing

The Illumina TruSeq Stranded Total RNA Sample Prep with Ribo-Zero Gold rRNA removal kit (Illumina, San Diego, CA) was used for RNA library preparation. Sequencing was performed at the High Throughput Genomics and Bioinformatics Analysis Shared Resource of the Huntsman Cancer Institute on an Illumina HiSeq 2500 (Illumina) system. The MultiQC report is included in supplemental Table 1.

Genome-wide DNA methylation profiling

DNA methylation was assessed at whole-genome level in a cohort of 12 CMML patients and 12 age-matched controls with Illumina Infinium Human Methylation 450K BeadChip arrays, which interrogate >485 000 methylation sites per sample at single-nucleotide resolution, as previously described.29 

Analysis pipeline

Human GRCh38 FASTA and GTF files were downloaded from Ensembl,30  release 87, and the reference database was created using STAR,31  version 2.5.2b, with splice junctions optimized for 50-bp reads. Reads were aligned to the reference database by using STAR in 2-pass mode to output a BAM file sorted by coordinates. Mapped reads were assigned to annotated genes in the GTF, using featureCounts,32  version 1.5.1. The output files from FastQC, STAR, and featureCounts were summarized with MultiQC,33  to check for any sample outliers. Samples with unacceptably low mapping rates to coding regions in the MultiQC report were excluded from the analysis. The Huntsman Cancer Institute Bioinformatics Analysis Shared Resource assisted with data analysis.

Identification of differentially expressed genes

Differentially expressed genes were identified by using a 1% false discovery rate (FDR) with DESeq2,34  version 1.16.0. The analysis included 12 CMML patients, 12 age-matched (old) healthy controls, and 16 young controls. Pairwise comparisons were carried out between transcriptomes of the 3 defined groups. The resulting sets of differentially expressed genes were overlapped in a Venn diagram, and changes were identified that were specific among the comparisons. CMML-specific genes were identified by subtracting all genes with significant differences in the young controls vs old controls comparison from the CMML vs old controls comparison.

Identification of DMRs

The DNA methylation arrays from 12 CMML patients and 12 age-matched controls were analyzed with the minfi package.35  The arrays were normalized using the functional normalization algorithm and quality control checked using interactive plots from the shinyMethyl package.36  Probes containing a single-nucleotide polymorphism at the CpG interrogation or at a single-nucleotide extension were removed. The bumphunter package37  was then used to find and annotate differentially methylated regions (DMRs). The DMRs near promoters and significant CMML genes in RNA-Seq data were selected in a subsequent step for analysis. The β-value represents a quantitative measure of the DNA methylation level of specific CpG-sites and ranges from 0 (completely unmethylated) to 1 (completely methylated).

Gene Set Enrichment Analysis

Broad Institute Gene Set Enrichment Analysis (GSEA),38  version 3.0 (http://www.broadinstitute.org/gsea/), was used to determine whether defined sets of genes exhibit a statistically significant enrichment in their distribution within our ranked gene list. The rank lists were compared with gene sets in MSigDB 6.1 including pathways from the Kyoto Encyclopedia of Genes and Genomes database.

Ingenuity Pathway Analysis

Pathway-associated gene analyses were generated using commercially available software (Qiagen Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis).

Clinical characteristics of the CMML patients

Consecutive newly diagnosed patients seen at the Hematologic Malignancies Clinic at the Huntsman Cancer Institute were included in our study (8 samples). The remaining 4 samples were obtained from consecutive newly diagnosed CMML patients treated at Oregon Health ’ Science University and were enrolled on a prospective clinical trial of 5-azacitidine (www.clinicaltrials.gov #NCT01350947). Based on the percentage of blast cells in the blood and bone marrow, CMML can be graded from 0-2, with 2 being the most severe. Of the 12 CMML patients studied, 8 were CMML-0, and 3 CMML-2.39  One patient was unclassifiable because no bone marrow biopsy was available. Seven cases were classified as proliferative and 4 as dysplastic CMML. The karyotype was low risk for 9 patients, high risk for 2 patients and unknown for 1 patient. TET2 mutations were detected in 75% of patients by whole exome sequencing. For clinical features see Table 1.

Table 1.

Clinical characteristics at diagnosis

PatientAge, ySexWhite blood cells, 109/LHemoglobin, g/dLPlatelets, 109/LBlood monocytes, % (109/L)Bone marrow blasts, %Circulating immature myeloid cellsWHO subtypeFAB subtypeKaryotypeCytogenetic riskMutated genes
57 Male 29 13.6 97 29 8.4 18 No CMML-2 Proliferative 46,XY[20] Low KRAS, SRSF2, TET2 
73 Male 44.2 11.7 30 3.9 Yes CMML-0 Proliferative 46,XY[20] Low CSF3R, TET2 
68 Female 15.1 13.2 51 41 6.0 No CMML-0 Proliferative 47,XX,+8[6]/46,XX[14] High TET2 
80 Female 15.4 12.3 255 35 5.3 1.4 No CMML-0 Proliferative 46,XX[20] Low PHF6, RUNX1, TET2, CBL 
70 Male 62 12.1 132 22 13.0 18 Yes CMML-2 Proliferative 46,XY[20] Low RUNX1, SRSF2, DNMT3A 
61 Female 20 12.2 353 16 3.2 Yes CMML-0 Proliferative 46,XX[20] Low ASXL1, SETBP1 
67 Male 34 10.7 126 36 12.0 No CMML-0 Dysplastic 46,XX[20] Low ASXL1, NRAS, RUNX1, TET2, ZRSR2 
74 Male 17.3 11.2 21 38 6.5 19 Yes CMML-2 Proliferative 45,XY,-7[3]/46,XY[20] High TET2, RUNX1, CBL 
87 Male 8.4 8.2 652 46 3.8 NA* Yes NA* Dysplastic NA* NA NA 
10 76 Male 3.3 10.5 126 38 1.2 No CMML-0 Dysplastic 46,XY[20] Low TET2, CBL, CSF3R 
11 72 Male 11.6 10.2 18 17 1.9 Yes CMML-0 Dysplastic 46,XY[20] Low TET2, CSF3R, SRSF2 
12 72 Male 5.2 13.6 87 24 1.2 Yes CMML-0 Dysplastic 46,XY[20] Low TET2, RUNX, SRSF2 
PatientAge, ySexWhite blood cells, 109/LHemoglobin, g/dLPlatelets, 109/LBlood monocytes, % (109/L)Bone marrow blasts, %Circulating immature myeloid cellsWHO subtypeFAB subtypeKaryotypeCytogenetic riskMutated genes
57 Male 29 13.6 97 29 8.4 18 No CMML-2 Proliferative 46,XY[20] Low KRAS, SRSF2, TET2 
73 Male 44.2 11.7 30 3.9 Yes CMML-0 Proliferative 46,XY[20] Low CSF3R, TET2 
68 Female 15.1 13.2 51 41 6.0 No CMML-0 Proliferative 47,XX,+8[6]/46,XX[14] High TET2 
80 Female 15.4 12.3 255 35 5.3 1.4 No CMML-0 Proliferative 46,XX[20] Low PHF6, RUNX1, TET2, CBL 
70 Male 62 12.1 132 22 13.0 18 Yes CMML-2 Proliferative 46,XY[20] Low RUNX1, SRSF2, DNMT3A 
61 Female 20 12.2 353 16 3.2 Yes CMML-0 Proliferative 46,XX[20] Low ASXL1, SETBP1 
67 Male 34 10.7 126 36 12.0 No CMML-0 Dysplastic 46,XX[20] Low ASXL1, NRAS, RUNX1, TET2, ZRSR2 
74 Male 17.3 11.2 21 38 6.5 19 Yes CMML-2 Proliferative 45,XY,-7[3]/46,XY[20] High TET2, RUNX1, CBL 
87 Male 8.4 8.2 652 46 3.8 NA* Yes NA* Dysplastic NA* NA NA 
10 76 Male 3.3 10.5 126 38 1.2 No CMML-0 Dysplastic 46,XY[20] Low TET2, CBL, CSF3R 
11 72 Male 11.6 10.2 18 17 1.9 Yes CMML-0 Dysplastic 46,XY[20] Low TET2, CSF3R, SRSF2 
12 72 Male 5.2 13.6 87 24 1.2 Yes CMML-0 Dysplastic 46,XY[20] Low TET2, RUNX, SRSF2 

FAB, French-American-British, NA, not available, WHO, World Health Organization.

*

No bone marrow biopsy and karyotype performed at diagnosis.

Identification of gene expression differences between CMML patients and healthy controls

We performed RNAseq on CD14+ monocytes from treatment-naive CMML patients and healthy controls (old and young) and compared gene expression between CMML vs old controls (Figure 1A,C) and young controls vs old controls (Figure 1B,D), using stringent statistical criteria (FDR < 0.01) as cutoff. Considering only protein-coding transcripts, 1565 genes were upregulated and 915 downregulated in monocytes from CMML patients vs old controls. The 24 genes with the most significant differential expression (sorted by P value: 12 upregulated and 12 downregulated) are summarized in Table 2. The comparison between young and old controls, using the same criteria, identified 1043 upregulated and 308 downregulated transcripts. The 24 genes with the most significant differential expression (sorted by P value: 12 upregulated and 12 downregulated) in the young controls vs old controls comparison are summarized in Table 2. Reynolds et al40  recently reported on age-related transcriptome-wide changes in monocytes collected from 1264 participants (aged 55-94 years) in the Multi-Ethnic Study of Atherosclerosis (MESA) cohort. A total of 2704 genes were differentially expressed with chronological age (FDR ≤ 0.001). Despite the higher age bracket of the MESA cohort, there was statistically significant overlap with the genes differentially expressed in young controls vs old controls in our study (627 of 1351 in common, Fisher’s exact test, P = 1.028 × 10−7), providing independent validation of our control dataset (supplemental Table 2A). In contrast, the overlap with genes differentially expressed between CMML and old controls failed to reach statistical significance (657 of 2840 in common; Fisher’s exact test P = .0614; supplemental Table 2B).

Figure 1.

Identification of differentially expressed genes in CD14+cells from CMML patients and old and young controls. Highly purified CD14+ monocytes were subjected to RNA sequencing. (A) Differentially expressed genes in CMML CD14+ monocytes compared with old healthy controls were graphed according to fold change and FDR-adjusted P value (volcano plot). (B) Differentially expressed genes in CD14+ monocytes of healthy young controls compared with healthy old controls were graphed according to fold change and FDR-adjusted P value (volcano plot). (C) Heat map depicting the top 12 most significantly upregulated genes and the top 12 most significantly downregulated genes in the CMML vs old controls comparison clustered by rlog values. (D) Heat map depicting the top 12 most significantly upregulated genes and the top 12 most significantly downregulated genes in the young controls vs old controls comparison clustered by rlog values.

Figure 1.

Identification of differentially expressed genes in CD14+cells from CMML patients and old and young controls. Highly purified CD14+ monocytes were subjected to RNA sequencing. (A) Differentially expressed genes in CMML CD14+ monocytes compared with old healthy controls were graphed according to fold change and FDR-adjusted P value (volcano plot). (B) Differentially expressed genes in CD14+ monocytes of healthy young controls compared with healthy old controls were graphed according to fold change and FDR-adjusted P value (volcano plot). (C) Heat map depicting the top 12 most significantly upregulated genes and the top 12 most significantly downregulated genes in the CMML vs old controls comparison clustered by rlog values. (D) Heat map depicting the top 12 most significantly upregulated genes and the top 12 most significantly downregulated genes in the young controls vs old controls comparison clustered by rlog values.

Close modal
Table 2.

Summary of the top 12 upregulated and downregulated genes in CMML vs old controls and young controls vs old controls comparisons

CMML vs old controlsYoung controls vs old controls
UpregulatedDownregulatedUpregulatedDownregulated
SNAI1 JSRP1 MME SNX32 
CDKN1A TLDC2 CXCR2 TLDC2 
PHLDA1 SNX22 TNFRSF10C JSRP1 
ATF3 RPA1 FCGR3B SH3D19 
CXCL2 HNRNPUL2 CXCL1 SNX22 
TEX14 UNC119B EPHB1 TYW1B 
OLR1 HIF1AN PHOSPHO1 RPS25 
APBB3 CAMKK2 CHI3L1 RPS24 
C12orf57 SLCO3A1 BAIAP3 FAF1 
HSF4 PTGDR2 AOC3 STX18 
HIST1H3B SH2D1B ASPRV1 SP140L 
PLK2 LYPD2 PLIN4 HK1 
CMML vs old controlsYoung controls vs old controls
UpregulatedDownregulatedUpregulatedDownregulated
SNAI1 JSRP1 MME SNX32 
CDKN1A TLDC2 CXCR2 TLDC2 
PHLDA1 SNX22 TNFRSF10C JSRP1 
ATF3 RPA1 FCGR3B SH3D19 
CXCL2 HNRNPUL2 CXCL1 SNX22 
TEX14 UNC119B EPHB1 TYW1B 
OLR1 HIF1AN PHOSPHO1 RPS25 
APBB3 CAMKK2 CHI3L1 RPS24 
C12orf57 SLCO3A1 BAIAP3 FAF1 
HSF4 PTGDR2 AOC3 STX18 
HIST1H3B SH2D1B ASPRV1 SP140L 
PLK2 LYPD2 PLIN4 HK1 

To identify CMML-specific changes, we cross-referenced the resulting sets of differentially expressed genes in the young controls vs old controls and CMML vs old controls comparisons (Figure 2A). Removing the 360 genes differentially expressed between young and old controls from the list of genes with differential expression between CMML and old controls resulted in a set of 2480 deregulated genes specific to CMML. We next performed unsupervised hierarchical clustering, using the CMML gene signature as input (Figure 2B). CMML patients clearly separated from controls, whereas young and old controls clustered together, as expected after removal of the genes with differential expression between them. The top 25 differentially regulated genes are listed in Figure 2C. For validation, we measured the expression of 4 transcripts by quantitative PCR; all were confirmed (Figure 2D).

Figure 2.

Identification of a CMML transcriptional signature. (A) Venn diagram illustrating the number of differentially expressed genes in CMML vs old controls (blue) compared with young controls vs old controls (red) using an FDR < 0.01. (B) Unsupervised clustering heat map of the 2480 genes that are unique to the CMML vs old controls set, shown in CMML and normal old and young control monocytes. (C) Heat map of the top 25 most significant genes of the CMML signature. (D) RNAseq and qPCR gene expression normalized values for TEX14, PLK2, CAMKK2, RPA1, and CXCL8 transcripts. P values were calculated using the nonparametric Mann-Whitney U test. **P < 0.01; ****P < 0.0001.

Figure 2.

Identification of a CMML transcriptional signature. (A) Venn diagram illustrating the number of differentially expressed genes in CMML vs old controls (blue) compared with young controls vs old controls (red) using an FDR < 0.01. (B) Unsupervised clustering heat map of the 2480 genes that are unique to the CMML vs old controls set, shown in CMML and normal old and young control monocytes. (C) Heat map of the top 25 most significant genes of the CMML signature. (D) RNAseq and qPCR gene expression normalized values for TEX14, PLK2, CAMKK2, RPA1, and CXCL8 transcripts. P values were calculated using the nonparametric Mann-Whitney U test. **P < 0.01; ****P < 0.0001.

Close modal

Functional annotation of the top 25 gene expression signature of CMML monocytes revealed upregulation of several genes involved in inflammation. For instance, activating transcription factor 3 (ATF3) expression is induced by various stress signals and plays an important role in modulating inflammatory response.41 CXCL2 encodes the CXCR2 ligand, macrophage inflammatory protein 2 (MIP2) and is regulated by the granulocyte colony stimulating factor (G-CSF STAT3 signaling cascade.42  Oxidized low-density lipoprotein receptor-1 (OLR-1), is a type II transmembrane receptor of the C-type lectin family that induces oxidative stress, activates tumor necrosis factor α TNFα/nuclear factor-κB (NFκB), upregulates the expression of monocyte chemoattractant protein-1 (MCP1), a key mediator of inflammatory processes, and has a proinflammatory role in cardiovascular disease.43,-45  calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) transcripts were downregulated in CMML CD14+ monocytes. As previous studies have demonstrated, a crucial role of the CAMKK2-5′-AMP–activated protein kinase catalytic subunit-1 (PRKAA1) axis in the differentiation of monocytes into macrophages, our data are consistent with the observation that the differentiation from monocytes to macrophages is compromised in CMML.46,47  Interestingly, we also found upregulation of testis-expressed sequence 14 (TEX14) in CMML monocytes. TEX14 is a tumor-testis antigen whose expression is restricted to malignancy and normal germ cells and may represent a target for immunotherapy.48 

For an unbiased functional pathway annotation, we used the 2480 CMML-specific differentially expressed genes as input into the Ingenuity Pathway Analysis (IPA) application (Figure 3). In addition to enrichment for pathways involved in the regulation of proliferation and cell cycle, inflammatory pathways were highly represented, including lymphotoxin-β receptor signaling, interleukin (IL)-6 signaling, NFκB activation, CD40 signaling, IL-17A signaling and lipopolysaccharide-simulated MAPK signaling. These data suggest that, compared with normal aged-matched controls, CMML monocytes exhibit strongly increased inflammatory signaling. JAK2/STAT signaling was also enriched, in accordance with the granulocyte-macrophage colony-stimulating factor (GM-CSF) hypersensitivity of many CMML cases.49 

Figure 3.

IPA. The 2480-gene, CMML-specific signature was used as input into IPA. Pathways were ranked according to significance.

Figure 3.

IPA. The 2480-gene, CMML-specific signature was used as input into IPA. Pathways were ranked according to significance.

Close modal

A focused comparison of 378 cytokine/chemokine-encoding genes (defined using Gene Ontology Annotations, GO:0005125) identified 46 genes with significantly different expression in the comparison between CMML and old controls or between old and young controls (supplemental Figure 2A-B). Hierarchical clustering, using only this set of cytokine/chemokine genes, clearly distinguished among the 3 groups (supplemental Figure 2C; supplemental Table 3). Interestingly, a group of cytokine/chemokine-encoding genes was upregulated in the young controls compared with both old controls and CMML, including positive (eg, IL-16 and IL-32) and negative (eg, IL-1RN) regulators of inflammation (supplemental Figure 2C). Cytokines/chemokine transcripts with increased expression specific to CMML included IL-1A, IL-1B, IL-10, CSF110, and IL-6, among others (supplemental Figure 2D). CXCL8, IL-6, and VEGFA transcripts overlap with the high levels of cytokines present in CMML patients’ plasma compared with that of healthy controls identified by Niyongere et al.50  Altogether, these data show that CMML monocytes exhibit a unique cytokine/chemokine expression profile that is distinct from both young and old controls.

TET2 is the most commonly mutated CHIP gene and Tet2-deficient monocytes promote atherosclerosis in low-density lipoprotein receptor–deficient mice.51,52  Liu and colleagues53  used the Multiethnic Study of Atherosclerosis (MESA) cohort to establish correlations between the transcriptional signature of monocytes and the risk of human atherosclerosis, as defined by carotid plaque burden and coronary artery calcification. They identified 21 genes with expression associated with carotid plaque score and 104 genes associated with coronary artery calcification (FDR < 0.2). We tested whether the gene expression signature of CMML monocytes may overlap with the atherosclerosis signature of the MESA cohort. Twenty-four differentially expressed genes, detailed in supplemental Table 4, were shared, including ARID5B (expression increased in CMML monocytes and in atherosclerosis). ARID5B is a transcriptional coactivator implicated as a mediator of cardiovascular risk.53  However, as a whole, this overlap did not reach statistical significance (P = .1676).

Methylation differences between CMML and healthy monocytes and integrated analysis

Several studies have shown that CMML mononuclear cells exhibit profound changes in DNA methylation compared with normal controls.11,54,55  However, data specific to monocytes are unavailable. To understand to what extent differences in methylation correlate with differentially expressed genes identified by RNAseq, we compared DNA methylation between CMML and aged-matched normal CD14+ monocytes using high-density methylation arrays. DNA methylation profiling revealed 909 DMRs between CMML and age-matched controls (FDR < 0.05 and methylation difference ≥20%), with most regions being hypermethylated in CMML monocytes (Figure 4A). Of the hypermethylated DMRs, 37% were intronic, 22% were exonic, 14% were in promoter regions, 10% were downstream, 10% were upstream and the remaining 7% were in 5′-untranslated regions.

Figure 4.

DNA methylation analysis of CD14+monocytes from CMML patients and age-matched healthy controls. (A) Comparison of methylation between CMML and old controls. Most regions are hypermethylated. (B) Unsupervised clustering of the 909 DMRs identified in CD14+ monocytes from CMML patients vs healthy age-matched controls, using a P value cutoff of .05. The heat map represents average β-values pseudocolored blue (hypomethylated) to red (hypermethylated). TET2 mutational status is shown. (C) Intersection of methylation differences (909 DMRs in the CMML vs old contrast) and log2 expression fold changes (2480 CMML-specific genes) identified 99 DMRs (94 unique genes, since some genes have 2 or more DMRs). The heat map represents β-values from the 99 DMRs. (D) A total of 124 DMRs localized near gene promoters, and 74 of these genes were also detected in RNA sequencing comparison in monocytes from CMML patients vs old controls. The graphic shows methylation differences in the promoter regions of those 74 DMRs vs fold changes in gene expression. Only 10 of the genes were specific to CMML signature (shown in blue), whereas the remainder represented age-related changes.

Figure 4.

DNA methylation analysis of CD14+monocytes from CMML patients and age-matched healthy controls. (A) Comparison of methylation between CMML and old controls. Most regions are hypermethylated. (B) Unsupervised clustering of the 909 DMRs identified in CD14+ monocytes from CMML patients vs healthy age-matched controls, using a P value cutoff of .05. The heat map represents average β-values pseudocolored blue (hypomethylated) to red (hypermethylated). TET2 mutational status is shown. (C) Intersection of methylation differences (909 DMRs in the CMML vs old contrast) and log2 expression fold changes (2480 CMML-specific genes) identified 99 DMRs (94 unique genes, since some genes have 2 or more DMRs). The heat map represents β-values from the 99 DMRs. (D) A total of 124 DMRs localized near gene promoters, and 74 of these genes were also detected in RNA sequencing comparison in monocytes from CMML patients vs old controls. The graphic shows methylation differences in the promoter regions of those 74 DMRs vs fold changes in gene expression. Only 10 of the genes were specific to CMML signature (shown in blue), whereas the remainder represented age-related changes.

Close modal

As previous studies in unselected CMML blood or bone marrow cells have attributed alterations in DNA methylation to mutations in epigenetic regulators, we initially sought to correlate somatic mutation status with methylation patterns in our cohort.56 TET2 is one of the most frequently mutated genes in CMML,57  and TET2 mutations were detected in 9 of 12 (75%) of our CMML patients. Unsupervised clustering analysis of the 909 DMRs revealed that DNA methylation patterns among all CMML patients differentiate between CMML patients with TET2 mutations and patients with wild-type TET2 (Figure 4B), with most of the regions significantly hypermethylated in CMML monocytes. IPA showed significant enrichment of DMRs in the sonic hedgehog signaling, dopamine-DARPP32 feedback in cAMP signaling, gap junction signaling, and adipogenesis pathways (supplemental Figure 3). Next, we cross-referenced methylation changes to mRNA expression, to identify those genes with coordinated regulation. Using the 909 DMRs and the 2480 CMML-specific genes, we identified 99 shared regions in 94 unique genes, including several genes with 2 or more DMRs. To determine whether the methylation status of this selected gene set would be sufficient to distinguish between CMML and controls, we performed unsupervised clustering analysis using the β-values from the 94 DMRs. The resulting CMML-specific methylation patterns completely separated CMML and normal monocytes (Figure 4C). The 3 TET2 wild-type CMML cases had a methylation profile closer to that of the old controls than to TET2-mutated CMML, resembling a hybrid between old normal and CMML. This is consistent with previously reported methylation patterns of CMML bone marrow mononuclear cells that separated TET2 wild-type from TET2 mutant cases.58 

Methylation of CpG islands in gene promoters is a well-established mechanism of transcriptional silencing in cancer.59  In CMML, hypomethylating agents can induce demethylation of CpG islands within promoters, thereby leading to reexpression of silenced tumor-suppressor genes.55  Therefore, although most of the CMML-specific DMRs are intronic (supplemental Figure 4), we focused specifically on DMRs within genes promoters to understand which CMML-specific genes may be controlled by promoter methylation. We observed significantly coordinated changes in CMML-specific mRNA transcripts and DNA methylation promoter regions for 10 genes (Figure 4D). However, only half of these genes followed the predicted regulation pattern: AOAH, SERINC5, TAF3, and AHCYL1 were hypermethylated and downregulated; TUBA1B was hypomethylated and upregulated. In contrast, MS4A3, TNF, VCAM1, and IFT80 were hypermethylated and upregulated, and PITPNA was hypomethylated and downregulated. Thus, gene expression in CMML monocytes is controlled by complex mechanisms, which may or may not include promoter methylation.

CMML is a genetically heterogeneous MDS/MPN overlap disorder with the unifying phenotypic characteristic of persistent monocytosis.39  Given the steep increase of CMML incidence with advanced age, the myelomonocytic differentiation bias of aged bone marrow, and the overlap of genes mutated in CMML with those mutated in ARCH, CMML could be regarded as the malignant conversion of ARCH.60  However, it is unknown which characteristics of CMML monocytes reflect ageing rather than the disease state of CMML, and the available RNAseq data are limited to 1 small study that provided little detail.61  Gene expression profiles of CMML monocytes and non–age-matched healthy controls using gene expression arrays have been described.62  Identifying features specific for CMML monocytes is critical from the perspective of understanding their biology and for therapeutically targeting CMML over normal monocytes. To approach this question, we performed RNAseq on highly purified blood monocytes from CMML patients, old controls, and young controls. Approximately one fourth of the genes differentially expressed between normal young and old monocytes overlapped with those differentially expressed between CMML and aged-matched monocytes, underscoring the importance of our 3-way comparison to isolate the transcriptional signature of CMML monocytes (Figure 2A). IPA of the CMML monocyte transcriptome identified multiple inflammatory signaling pathways, most of which were dominated by upregulated genes (Figure 3). Specifically, 4 pathways related to the TNF family signaling were enriched in CMML monocytes, including lymphotoxin-β receptor, IL-17A signaling in fibroblasts, NFκB activation by viruses, CD40, and CD27. High-mobility group protein B1 (HMGB1) and IL-6 signaling were additional highly proinflammatory pathways identified in the transcriptome of the CMML monocytes. The latter is consistent with recent reports of increased IL-6 expression in Tet2−/− mouse monocytes and macrophages and with our previous work on cytokine expression in CMML vs normal bone marrow.50,-52  HMGB1 is a nuclear protein that acts as a danger associated molecular pattern (DAMP) and extracellular trigger of inflammation through binding to Toll-like receptors 2 and 4.63  Another proinflammatory pathway characteristic of CMML monocytes is IL-17 signaling (Figure 3). This observation is consistent with our report of increased IL-17 expression in CMML vs normal bone marrow.50  IL-17 has been implicated as a master regulator of pathologic inflammation in psoriasis and other autoimmune diseases, and stimulates the production of proinflammatory cytokines from monocytes (TNFα, IL-1B, IL-6, GM-CSF, and G-CSF).64,65  A second cluster of enriched pathways in CMML monocytes involves cell cycle regulation and DNA damage, with pololike kinases as the most significantly enriched pathway overall, and a third cluster involves proliferative signaling (JAK/STAT and lipopolysaccharide-stimulated MAPK signaling), the latter linking inflammation to proliferation. At the single cytokine level, we also found elevated IL-1A and -1B expression in CMML monocytes, again consistent with data from Tet2−/− mouse monocytes.51,52  In aggregate, these data show that cytokine/chemokine expression patterns in CMML monocytes are highly aberrant. More work is needed to identify the regulatory circuits governing expression of these inflammatory cytokines.

Recent reports have suggested that most myeloid neoplasms, including CMML, originate from age-related precursor lesions that initially manifest as small populations of mutated blood cells in individuals with normal blood counts.17,,,,-22  The term CHIP is typically used in case of mutations with a variant allele frequency of >2%, the sensitivity limit of many next-generation sequencing assays, whereas ARCH refers to the presence of mutated cells at any level.66  When highly sensitive and specific assays are used, DNMT3A and TET2 mutations are detectable at low frequency in almost 100% of healthy middle-aged women.67  It is conceivable that mutated cells are more resistant to an inflammatory environment, providing them with a selective advantage over normal competitors. Indeed, several recent reports support this concept in the case of Tet2−/− mouse hematopoietic stem and progenitor cells.26,68,69  Similarly, we have previously shown that TNFα promotes expansion of JAK2V617F expressing cells in myelofibrosis.70,71  Current thinking holds that multiple mechanisms promote an inflammatory state in the aging bone marrow that manifests as a myelomonocytic differentiation bias (reviewed by Elias et al72 ). Unexpectedly, we found that expression of several inflammatory cytokines such as IL-16 and -32 was higher in young compared with old control monocytes, suggesting that not all aged-related transcriptional changes in monocytes are proinflammatory and that additional cell types are involved in promoting age-associated inflammation. Our study is limited by the lack of clinical information and mutational profiling for the normal controls, and as such we cannot exclude the presence of inflammatory conditions and/or the presence of CHIP/ARCH in the control population. This may explain why 3 of the normal controls grouped with CMML monocytes on unsupervised clustering. Alternatively, mutant cells may promote their own expansion by releasing proinflammatory cytokines into their microenvironment, consistent with previous reports of complex interactions between various CMML cell types. For instance, dysplastic granulocytes generated by the CMML clone secrete large amounts of α-defensin 1 to 3 peptides, which in turn inhibit CSF1-driven differentiation of peripheral blood monocytes into macrophages.47  Mapping the interactions between mutant and residual normal hematopoietic stem and progenitor cells, and mature cells, including monocytes, are critical for understanding their roles in clonal selection. The consequences of CMML-driven inflammation may reach beyond the hematopoietic system. Kyoto Encyclopedia of Genes and Genomes pathway analysis of the CMML-specific deregulated genes identified increased systemic lupus erythematosus, type 1 diabetes mellitus, and Toll-like receptor signaling, whereas the primary immunodeficiency cluster was reduced compared with age-matched controls (supplemental Figure 5). This suggests that monocytes may contribute to the increased incidence of autoimmune disease in patients with CMML.73,-75  In our cohort, 4 CMML patients had history of inflammatory and autoimmune diseases (supplemental Table 5).

CMML is one of the most aggressive of the MDS/MPN overlap disorders, with a median overall survival (OS) of only 20 to 32 months.76,-78  Current therapy is based on hypomethylating agents such as 5-azacytidine, with response rates of ∼50%, but no certain prolongation of OS.79  Although transformation to acute myeloid leukemia and complications of cytopenias are important contributors to poor survival, many patients succumb to comorbid conditions.80  The inclusion of somatic mutations in prognostic scoring systems has failed decisively to improve prognostication, suggesting that additional leukemia and/or host factors influence prognosis.9,76  In contrast to other subtypes of MPN and MDS, accumulation of monocytes is a hallmark of CMML, and our data show that these monocytes are highly proinflammatory compared with normal age-matched controls. As CHIP with a relatively small number of mutated monocytes is sufficient to increase cardiovascular risk, including mortality from ischemic heart disease, it is conceivable that the massive expansion of proinflammatory monocytes in CMML constitutes a considerable vascular risk factor.19,51,52,81  In support of this, a retrospective study found a twofold increased risk for cardiovascular events in CMML compared with patients with chronic lymphocytic leukemia.27  Another study reported a positive correlation between blood monocyte counts and end organ damage in CMML patients, including chronic kidney disease and cardiovascular events, consistent with the reported role of blood monocytes in atherosclerosis and the positive correlation between absolute monocyte counts and cardiovascular risk.82,83  Based on this, we compared the transcriptional signature of atherosclerosis reported from the MESA cohort and our CMML gene signature and identified a set of overlapping genes, including the atherosclerosis regulator ARID5B, suggesting pathogenetic similarities.40,53,84,,-87 

Profound changes in DNA methylation contribute to CMML pathogenesis. These changes are mostly genotype agnostic and correlate with response to 5-azacytidine and prognosis.11  Most differentially methylated regions were intragenic (supplemental Figure 4), consistent with published data.11  However, as cross-referencing of CpG island methylation in promoter regions of genes with CMML-specific gene expression changes revealed relatively little overlap and no consistent regulation pattern (Figure 4D), mechanisms other than promoter methylation must account for the bulk of the CMML-specific transcriptional dysregulation, including the dysregulation of cytokine expression.10,54 

In summary, we report the first systematic analysis of CMML monocyte transcriptional profiles that uses aged-matched and young controls to separate CMML-specific changes from aged-related changes. However, our study comprises a relatively low number of patients (the majority CMML-0 subtype), and a direct comparison with a reactive monocytosis control group is absent. Further validation with larger number of patients and functional studies are needed. Although our cohort is not large enough to fully account for the genetic complexity of CMML, our findings identify CMML monocytes as highly proinflammatory compared with controls. The detrimental effects of proinflammatory monocytes on extrahematopoietic organs, specifically the cardiovascular system, raise the question of whether even asymptomatic and low risk CMML patients would benefit from monocyte-reducing and/or anti-inflammatory therapy such as the IL-1β antibody canakinumab.88 

The authors thank Claire Davis, who provided the original illustrations for the visual abstract.

Research reported in this publication used the High-Throughput Genomics and Bioinformatic Analysis Shared Resource at Huntsman Cancer Institute at the University of Utah and was supported by the National Institutes of Health (NIH), National Cancer Institute Grant P30CA042014. The research was also supported by NIH, National Cancer Institute Grant R01CA178397 (M.W.D. and T.O.) and the V Foundation for Cancer Research (M.W.D. and T.O.). J.S.K. was and D.Y. is supported by a Special Fellow Award from the Leukemia and Lymphoma Society.

The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Contribution: A.F. supervised the study, performed experiments and computational analyses, and cowrote the manuscript; A.D.P supervised collection of the samples (CMML patients, old and young donors), supervised sample purification and analysis, prepared figures and edited the manuscript; D.Y. designed and performed experiments; J.S.K. designed and performed experiments; S.K.T. collected and provided the patients’ clinical information; H.T. contributed to data collection and interpretation; J.M.A. performed the qPCR validation experiments; T.O. supervised the study; and M.W.D initiated and supervised the study and cowrote the manuscript.

Conflict-of-interest disclosure: M.W.D. is a consultant or advisory board member for Novartis, Pfizer, Galena Biopharma, Takeda Pharmaceutical Company, Blueprint Medicines, and Incyte. M.W.D.’s laboratory receives grant funding from Novartis and Pfizer. The remaining authors declare no competing financial interests.

Correspondence: Michael W. Deininger, Huntsman Cancer Institute, The University of Utah, 2000 Circle of Hope, Salt Lake City, UT 84112; e-mail: michael.deininger@hci.utah.edu.

The data reported in this article have been deposited in the Gene Expression Omnibus database (accession number GSE135902).

The full-text version of this article contains a data supplement.

1.
Solary
E
,
Itzykson
R
.
How I treat chronic myelomonocytic leukemia
.
Blood
.
2017
;
130
(
2
):
126
-
136
.
2.
Deininger
MWN
,
Tyner
JW
,
Solary
E
.
Turning the tide in myelodysplastic/myeloproliferative neoplasms
.
Nat Rev Cancer
.
2017
;
17
(
7
):
425
-
440
.
3.
Srour
SA
,
Devesa
SS
,
Morton
LM
, et al
.
Incidence and patient survival of myeloproliferative neoplasms and myelodysplastic/myeloproliferative neoplasms in the United States, 2001-12
.
Br J Haematol
.
2016
;
174
(
3
):
382
-
396
.
4.
Cui
Y
,
Tong
H
,
Du
X
, et al
.
TET2 mutations were predictive of inferior prognosis in the presence of ASXL1 mutations in patients with chronic myelomonocytic leukemia
.
Stem Cell Investig
.
2016
;
3
:
50
.
5.
Haferlach
T
,
Nagata
Y
,
Grossmann
V
, et al
.
Landscape of genetic lesions in 944 patients with myelodysplastic syndromes
.
Leukemia
.
2014
;
28
(
2
):
241
-
247
.
6.
Pellagatti
A
,
Boultwood
J
.
Splicing factor gene mutations in the myelodysplastic syndromes: impact on disease phenotype and therapeutic applications
.
Adv Biol Regul
.
2017
;
63
:
59
-
70
.
7.
Geissler
K
,
Jäger
E
,
Barna
A
, et al
.
Chronic myelomonocytic leukemia patients with RAS pathway mutations show high in vitro myeloid colony formation in the absence of exogenous growth factors
.
Leukemia
.
2016
;
30
(
11
):
2280
-
2281
.
8.
Kohlmann
A
,
Grossmann
V
,
Klein
HU
, et al
.
Next-generation sequencing technology reveals a characteristic pattern of molecular mutations in 72.8% of chronic myelomonocytic leukemia by detecting frequent alterations in TET2, CBL, RAS, and RUNX1
.
J Clin Oncol
.
2010
;
28
(
24
):
3858
-
3865
.
9.
Sallman
D
,
Padron
E
,
Vaupel
C
, et al
Prognostic impact of ASXL1 mutations in MDS and CMML [abstract]. Blood.
2015
;126(23). Abstract 1623.
10.
Palomo
L
,
Malinverni
R
,
Cabezón
M
, et al
.
DNA methylation profile in chronic myelomonocytic leukemia associates with distinct clinical, biological and genetic features
.
Epigenetics
.
2018
;
13
(
1
):
8
-
18
.
11.
Meldi
K
,
Qin
T
,
Buchi
F
, et al
.
Specific molecular signatures predict decitabine response in chronic myelomonocytic leukemia
.
J Clin Invest
.
2015
;
125
(
5
):
1857
-
1872
.
12.
Beerman
I
,
Bhattacharya
D
,
Zandi
S
, et al
.
Functionally distinct hematopoietic stem cells modulate hematopoietic lineage potential during aging by a mechanism of clonal expansion
.
Proc Natl Acad Sci USA
.
2010
;
107
(
12
):
5465
-
5470
.
13.
Challen
GA
,
Boles
NC
,
Chambers
SM
,
Goodell
MA
.
Distinct hematopoietic stem cell subtypes are differentially regulated by TGF-beta1
.
Cell Stem Cell
.
2010
;
6
(
3
):
265
-
278
.
14.
Chambers
SM
,
Shaw
CA
,
Gatza
C
,
Fisk
CJ
,
Donehower
LA
,
Goodell
MA
.
Aging hematopoietic stem cells decline in function and exhibit epigenetic dysregulation
.
PLoS Biol
.
2007
;
5
(
8
):
e201
.
15.
Cho
RH
,
Sieburg
HB
,
Muller-Sieburg
CE
.
A new mechanism for the aging of hematopoietic stem cells: aging changes the clonal composition of the stem cell compartment but not individual stem cells
.
Blood
.
2008
;
111
(
12
):
5553
-
5561
.
16.
King
KY
,
Goodell
MA
.
Inflammatory modulation of HSCs: viewing the HSC as a foundation for the immune response
.
Nat Rev Immunol
.
2011
;
11
(
10
):
685
-
692
.
17.
Busque
L
,
Patel
JP
,
Figueroa
ME
, et al
.
Recurrent somatic TET2 mutations in normal elderly individuals with clonal hematopoiesis
.
Nat Genet
.
2012
;
44
(
11
):
1179
-
1181
.
18.
Genovese
G
,
Kähler
AK
,
Handsaker
RE
, et al
.
Clonal hematopoiesis and blood-cancer risk inferred from blood DNA sequence
.
N Engl J Med
.
2014
;
371
(
26
):
2477
-
2487
.
19.
Jaiswal
S
,
Fontanillas
P
,
Flannick
J
, et al
.
Age-related clonal hematopoiesis associated with adverse outcomes
.
N Engl J Med
.
2014
;
371
(
26
):
2488
-
2498
.
20.
McKerrell
T
,
Park
N
,
Moreno
T
, et al;
Understanding Society Scientific Group
.
Leukemia-associated somatic mutations drive distinct patterns of age-related clonal hemopoiesis
.
Cell Reports
.
2015
;
10
(
8
):
1239
-
1245
.
21.
McKerrell
T
,
Vassiliou
GS
.
Aging as a driver of leukemogenesis
.
Sci Transl Med
.
2015
;
7
(
306
):
306fs38
.
22.
Xie
M
,
Lu
C
,
Wang
J
, et al
.
Age-related mutations associated with clonal hematopoietic expansion and malignancies
.
Nat Med
.
2014
;
20
(
12
):
1472
-
1478
.
23.
Risques
RA
,
Kennedy
SR
.
Aging and the rise of somatic cancer-associated mutations in normal tissues
.
PLoS Genet
.
2018
;
14
(
1
):
e1007108
.
24.
Acuna-Hidalgo
R
,
Sengul
H
,
Steehouwer
M
, et al
.
Ultra-sensitive Sequencing Identifies High Prevalence of Clonal Hematopoiesis-Associated Mutations throughout Adult Life
.
Am J Hum Genet
.
2017
;
101
(
1
):
50
-
64
.
25.
Kristinsson
SY
,
Björkholm
M
,
Hultcrantz
M
,
Derolf
Å
,
Landgren
O
,
Goldin
LR
.
Chronic immune stimulation might act as a trigger for the development of acute myeloid leukemia or myelodysplastic syndromes
.
J Clin Oncol
.
2011
;
29
(
21
):
2897
-
2903
.
26.
Meisel
M
,
Hinterleitner
R
,
Pacis
A
, et al
.
Microbial signals drive pre-leukaemic myeloproliferation in a Tet2-deficient host
.
Nature
.
2018
;
557
(
7706
):
580
-
584
.
27.
Elbæk
MV
,
Sørensen
AL
,
Hasselbalch
HC
.
Cardiovascular disease in chronic myelomonocytic leukemia: do monocytosis and chronic inflammation predispose to accelerated atherosclerosis?
Ann Hematol
.
2019
;
98
(
1
):
101
-
109
.
28.
Guru Murthy
GS
,
Dhakal
I
,
Mehta
P
.
Incidence and survival outcomes of chronic myelomonocytic leukemia in the United States
.
Leuk Lymphoma
.
2017
;
58
(
7
):
1648
-
1654
.
29.
Asmar
F
,
Punj
V
,
Christensen
J
, et al
.
Genome-wide profiling identifies a DNA methylation signature that associates with TET2 mutations in diffuse large B-cell lymphoma
.
Haematologica
.
2013
;
98
(
12
):
1912
-
1920
.
30.
Zerbino
DR
,
Achuthan
P
,
Akanni
W
, et al
.
Ensembl 2018
.
Nucleic Acids Res
.
2018
;
46
(
D1
):
D754
-
D761
.
31.
Dobin
A
,
Davis
CA
,
Schlesinger
F
, et al
.
STAR: ultrafast universal RNA-seq aligner
.
Bioinformatics
.
2013
;
29
(
1
):
15
-
21
.
32.
Liao
Y
,
Smyth
GK
,
Shi
W
.
featureCounts: an efficient general purpose program for assigning sequence reads to genomic features
.
Bioinformatics
.
2014
;
30
(
7
):
923
-
930
.
33.
Ewels
P
,
Magnusson
M
,
Lundin
S
,
Käller
M
.
MultiQC: summarize analysis results for multiple tools and samples in a single report
.
Bioinformatics
.
2016
;
32
(
19
):
3047
-
3048
.
34.
Love
MI
,
Huber
W
,
Anders
S
.
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
.
Genome Biol
.
2014
;
15
(
12
):
550
.
35.
Aryee
MJ
,
Jaffe
AE
,
Corrada-Bravo
H
, et al
.
Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays
.
Bioinformatics
.
2014
;
30
(
10
):
1363
-
1369
.
36.
Fortin
JP
,
Fertig
E
,
Hansen
K
.
shinyMethyl: interactive quality control of Illumina 450k DNA methylation arrays in R
.
F1000 Res
.
2014
;
3
:
175
.
37.
Jaffe
AE
,
Murakami
P
,
Lee
H
, et al
.
Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies
.
Int J Epidemiol
.
2012
;
41
(
1
):
200
-
209
.
38.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
, et al
.
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci USA
.
2005
;
102
(
43
):
15545
-
15550
.
39.
Arber
DA
,
Orazi
A
,
Hasserjian
R
, et al
.
The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia
.
Blood
.
2016
;
127
(
20
):
2391
-
2405
.
40.
Reynolds
LM
,
Ding
J
,
Taylor
JR
, et al
.
Transcriptomic profiles of aging in purified human immune cells
.
BMC Genomics
.
2015
;
16
(
1
):
333
.
41.
Hai
T
,
Wolford
CC
,
Chang
YS
.
ATF3, a hub of the cellular adaptive-response network, in the pathogenesis of diseases: is modulation of inflammation a unifying component?
Gene Expr
.
2010
;
15
(
1
):
1
-
11
.
42.
Nguyen-Jackson
HT
,
Li
HS
,
Zhang
H
,
Ohashi
E
,
Watowich
SS
.
G-CSF-activated STAT3 enhances production of the chemokine MIP-2 in bone marrow neutrophils
.
J Leukoc Biol
.
2012
;
92
(
6
):
1215
-
1225
.
43.
Melgarejo
E
,
Medina
MA
,
Sánchez-Jiménez
F
,
Urdiales
JL
.
Monocyte chemoattractant protein-1: a key mediator in inflammatory processes
.
Int J Biochem Cell Biol
.
2009
;
41
(
5
):
998
-
1001
.
44.
Wang
J
,
Mi
JQ
,
Debernardi
A
, et al
.
A six gene expression signature defines aggressive subtypes and predicts outcome in childhood and adult acute lymphoblastic leukemia
.
Oncotarget
.
2015
;
6
(
18
):
16527
-
16542
.
45.
Arslan
C
,
Bayoglu
B
,
Tel
C
,
Cengiz
M
,
Dirican
A
,
Besirli
K
.
Upregulation of OLR1 and IL17A genes and their association with blood glucose and lipid levels in femoropopliteal artery disease
.
Exp Ther Med
.
2017
;
13
(
3
):
1160
-
1168
.
46.
Obba
S
,
Hizir
Z
,
Boyer
L
, et al
.
The PRKAA1/AMPKα1 pathway triggers autophagy during CSF1-induced human monocyte differentiation and is a potential target in CMML
.
Autophagy
.
2015
;
11
(
7
):
1114
-
1129
.
47.
Droin
N
,
Jacquel
A
,
Hendra
JB
, et al
.
Alpha-defensins secreted by dysplastic granulocytes inhibit the differentiation of monocytes in chronic myelomonocytic leukemia
.
Blood
.
2010
;
115
(
1
):
78
-
88
.
48.
Condomines
M
,
Hose
D
,
Rème
T
, et al
.
Gene expression profiling and real-time PCR analyses identify novel potential cancer-testis antigens in multiple myeloma
.
J Immunol
.
2009
;
183
(
2
):
832
-
840
.
49.
Padron
E
,
Painter
JS
,
Mailloux
AW
, et al
GM-CSF signaling abnormalities in chronic myelomonocytic leukemia [abstract]. Blood.
2011
;118(21):746. Abstract 1713.
50.
Niyongere
S
,
Lucas
N
,
Zhou
JM
, et al
.
Heterogeneous expression of cytokines accounts for clinical diversity and refines prognostication in CMML
.
Leukemia
.
2019
;
33
(
1
):
205
-
216
.
51.
Fuster
JJ
,
MacLauchlan
S
,
Zuriaga
MA
, et al
.
Clonal hematopoiesis associated with TET2 deficiency accelerates atherosclerosis development in mice
.
Science
.
2017
;
355
(
6327
):
842
-
847
.
52.
Jaiswal
S
,
Natarajan
P
,
Silver
AJ
, et al
.
Clonal Hematopoiesis and Risk of Atherosclerotic Cardiovascular Disease
.
N Engl J Med
.
2017
;
377
(
2
):
111
-
121
.
53.
Liu
Y
,
Reynolds
LM
,
Ding
J
, et al
.
Blood monocyte transcriptome and epigenome analyses reveal loci associated with human atherosclerosis
.
Nat Commun
.
2017
;
8
(
1
):
393
.
54.
Yamazaki
J
,
Jelinek
J
,
Lu
Y
, et al
.
TET2 Mutations Affect Non-CpG Island DNA Methylation at Enhancers and Transcription Factor-Binding Sites in Chronic Myelomonocytic Leukemia
.
Cancer Res
.
2015
;
75
(
14
):
2833
-
2843
.
55.
Merlevede
J
,
Droin
N
,
Qin
T
, et al
.
Mutation allele burden remains unchanged in chronic myelomonocytic leukaemia responding to hypomethylating agents
.
Nat Commun
.
2016
;
7
(
1
):
10767
.
56.
Yamazaki
J
,
Taby
R
,
Vasanthakumar
A
, et al
.
Effects of TET2 mutations on DNA methylation in chronic myelomonocytic leukemia
.
Epigenetics
.
2012
;
7
(
2
):
201
-
207
.
57.
Patnaik
MM
,
Zahid
MF
,
Lasho
TL
, et al
.
Number and type of TET2 mutations in chronic myelomonocytic leukemia and their clinical relevance
.
Blood Cancer J
.
2016
;
6
(
9
):
e472
.
58.
Pérez
C
,
Martínez-Calle
N
,
Martín-Subero
JI
, et al
.
TET2 mutations are associated with specific 5-methylcytosine and 5-hydroxymethylcytosine profiles in patients with chronic myelomonocytic leukemia
.
PLoS One
.
2012
;
7
(
2
):
e31605
.
59.
Kulis
M
,
Esteller
M
.
DNA methylation and cancer
.
Adv Genet
.
2010
;
70
:
27
-
56
.
60.
Mason
CC
,
Khorashad
JS
,
Tantravahi
SK
, et al
.
Age-related mutations and chronic myelomonocytic leukemia
.
Leukemia
.
2016
;
30
(
4
):
906
-
913
.
61.
Selimoglu-Buet
D
,
Wagner-Ballon
O
,
Saada
V
, et al;
Francophone Myelodysplasia Group
.
Characteristic repartition of monocyte subsets as a diagnostic signature of chronic myelomonocytic leukemia
.
Blood
.
2015
;
125
(
23
):
3618
-
3626
.
62.
Braun
T
,
Itzykson
R
,
Renneville
A
, et al;
Groupe Francophone des Myélodysplasies
.
Molecular predictors of response to decitabine in advanced chronic myelomonocytic leukemia: a phase 2 trial
.
Blood
.
2011
;
118
(
14
):
3824
-
3831
.
63.
Pilzweger
C
,
Holdenrieder
S
.
Circulating HMGB1 and RAGE as Clinical Biomarkers in Malignant and Autoimmune Diseases
.
Diagnostics (Basel)
.
2015
;
5
(
2
):
219
-
253
.
64.
Beringer
A
,
Noack
M
,
Miossec
P
.
IL-17 in Chronic Inflammation: From Discovery to Targeting
.
Trends Mol Med
.
2016
;
22
(
3
):
230
-
241
.
65.
Miossec
P
,
Korn
T
,
Kuchroo
VK
.
Interleukin-17 and type 17 helper T cells
.
N Engl J Med
.
2009
;
361
(
9
):
888
-
898
.
66.
Steensma
DP
,
Bejar
R
,
Jaiswal
S
, et al
.
Clonal hematopoiesis of indeterminate potential and its distinction from myelodysplastic syndromes
.
Blood
.
2015
;
126
(
1
):
9
-
16
.
67.
Young
AL
,
Challen
GA
,
Birmann
BM
,
Druley
TE
.
Clonal haematopoiesis harbouring AML-associated mutations is ubiquitous in healthy adults
.
Nat Commun
.
2016
;
7
(
1
):
12484
.
68.
Cai
Z
,
Kotzin
JJ
,
Ramdas
B
, et al
.
Inhibition of Inflammatory Signaling in Tet2 Mutant Preleukemic Cells Mitigates Stress-Induced Abnormalities and Clonal Hematopoiesis
.
Cell Stem Cell
.
2018
;
23
(
6
):
P833
-
P849.e5
.
69.
Abegunde
SO
,
Buckstein
R
,
Wells
RA
,
Rauh
MJ
.
An inflammatory environment containing TNFα favors Tet2-mutant clonal hematopoiesis
.
Exp Hematol
.
2018
;
59
:
60
-
65
.
70.
Fleischman
AG
,
Aichberger
KJ
,
Luty
SB
, et al
.
TNFα facilitates clonal expansion of JAK2V617F positive cells in myeloproliferative neoplasms
.
Blood
.
2011
;
118
(
24
):
6392
-
6398
.
71.
Heaton
WL
,
Senina
AV
,
Pomicter
AD
, et al
.
Autocrine Tnf signaling favors malignant cells in myelofibrosis in a Tnfr2-dependent fashion
.
Leukemia
.
2018
;
32
(
11
):
2399
-
2411
.
72.
Elias
HK
,
Bryder
D
,
Park
CY
.
Molecular mechanisms underlying lineage bias in aging hematopoiesis
.
Semin Hematol
.
2017
;
54
(
1
):
4
-
11
.
73.
Isoda
A
,
Sakurai
A
,
Ogawa
Y
, et al
.
Chronic inflammatory demyelinating polyneuropathy accompanied by chronic myelomonocytic leukemia: possible pathogenesis of autoimmunity in myelodysplastic syndrome
.
Int J Hematol
.
2009
;
90
(
2
):
239
-
242
.
74.
Peker
D
,
Padron
E
,
Bennett
JM
, et al
.
A close association of autoimmune-mediated processes and autoimmune disorders with chronic myelomonocytic leukemia: observation from a single institution
.
Acta Haematol
.
2015
;
133
(
2
):
249
-
256
.
75.
Saif
MW
,
Hopkins
JL
,
Gore
SD
.
Autoimmune phenomena in patients with myelodysplastic syndromes and chronic myelomonocytic leukemia
.
Leuk Lymphoma
.
2002
;
43
(
11
):
2083
-
2092
.
76.
Ball
M
,
List
AF
,
Padron
E
.
When clinical heterogeneity exceeds genetic heterogeneity: thinking outside the genomic box in chronic myelomonocytic leukemia
.
Blood
.
2016
;
128
(
20
):
2381
-
2387
.
77.
Duchmann
M
,
Yalniz
FF
,
Sanna
A
, et al
.
Prognostic Role of Gene Mutations in Chronic Myelomonocytic Leukemia Patients Treated With Hypomethylating Agents
.
EBioMedicine
.
2018
;
31
:
174
-
181
.
78.
Itzykson
R
,
Kosmider
O
,
Renneville
A
, et al
.
Prognostic score including gene mutations in chronic myelomonocytic leukemia
.
J Clin Oncol
.
2013
;
31
(
19
):
2428
-
2436
.
79.
Kubasch
AS
,
Platzbecker
U
.
Beyond the Edge of Hypomethylating Agents: Novel Combination Strategies for Older Adults with Advanced MDS and AML
.
Cancers (Basel)
.
2018
;
10
(
6
):
E158
.
80.
Buckstein
R
,
Wells
RA
,
Zhu
N
, et al
.
Patient-related factors independently impact overall survival in patients with myelodysplastic syndromes: an MDS-CAN prospective study
.
Br J Haematol
.
2016
;
174
(
1
):
88
-
101
.
81.
Dorsheimer
L
,
Assmus
B
,
Rasper
T
, et al
.
Association of Mutations Contributing to Clonal Hematopoiesis With Prognosis in Chronic Ischemic Heart Failure
.
JAMA Cardiol
.
2019
;
4
(
1
):
25
-
33
.
82.
Madjid
M
,
Awan
I
,
Willerson
JT
,
Casscells
SW
.
Leukocyte count and coronary heart disease: implications for risk assessment
.
J Am Coll Cardiol
.
2004
;
44
(
10
):
1945
-
1956
.
83.
Naqvi
K
,
Sasaki
K
,
Montalban-Bravo
G
, et al
. Impact of Clonal Hematopoiesis of Indeterminate Potential (CHIP) Associated Mutations and Risk of Comorbidities in Patients with Myelodysplastic Syndrome. Paper presented at American Society of Hematology Annual Meeting; 1-4 December 2018; San Diego, CA.
84.
Ding
J
,
Reynolds
LM
,
Zeller
T
, et al
.
Alterations of a Cellular Cholesterol Metabolism Network Are a Molecular Feature of Obesity-Related Type 2 Diabetes and Cardiovascular Disease
.
Diabetes
.
2015
;
64
(
10
):
3464
-
3474
.
85.
Ley
K
,
Miller
YI
,
Hedrick
CC
.
Monocyte and macrophage dynamics during atherogenesis
.
Arterioscler Thromb Vasc Biol
.
2011
;
31
(
7
):
1506
-
1516
.
86.
Reynolds
LM
,
Taylor
JR
,
Ding
J
, et al
.
Age-related variations in the methylome associated with gene expression in human monocytes and T cells
.
Nat Commun
.
2014
;
5
(
1
):
5366
.
87.
Reynolds
LM
,
Wan
M
,
Ding
J
, et al
.
DNA Methylation of the Aryl Hydrocarbon Receptor Repressor Associations With Cigarette Smoking and Subclinical Atherosclerosis
.
Circ Cardiovasc Genet
.
2015
;
8
(
5
):
707
-
716
.
88.
Ridker
PM
,
Everett
BM
,
Thuren
T
, et al;
CANTOS Trial Group
.
Antiinflammatory Therapy with Canakinumab for Atherosclerotic Disease
.
N Engl J Med
.
2017
;
377
(
12
):
1119
-
1131
.

Supplemental data