Background. Chronic myelomonocytic leukemia (CMML) is a rare myeloid neoplasm that shares both dysplastic and proliferative features and includes patients with highly heterogeneous clinical manifestations and prognosis. Despite known co-mutation patterns that are suggesting of a CMML diagnosis such as TET2 and SRSF2 mutations, the current diagnosis and classification schemes rely solely on morphological criteria that fail to capture the genomic heterogeneity of the disease. Here we perform an unsupervised analysis using only genomic features to generate a novel classification schema of CMML patients with unique clinical features. Methods. A cohort of 2,471 CMML patients defined according to WHO 2016 criteria from European and US centers was analyzed. Bayesian network analysis and hierarchical Dirichlet processes were used to identify genomic associations and subgroups as a basis to define a molecular classification of CMML. Bayesian networks allow to infer the structure of conditional dependencies among mutations. Dirichlet processes were applied to define clusters capturing broad dependencies among all gene mutations and cytogenetic abnormalities. Multivariate logistic regression analysis was applied to compare clinical and hematological characteristics among different groups. ( JCO 2021 39:11, 1223-1233) Results. Baseline clinical and molecular characteristics are summarized in Table 1. We identified 11 clusters defined according to specific genomic features (Table 2). Two clusters were primarily defined by the presence of mutations in splicing genes SRSF2 and ZRSR2, frequently associated with TET2 mutations. Patients assigned to these clusters showed myelodysplastic features according to French American British (FAB) criteria, had a low rate of peripheral blood cytopenias, and low blast counts. Additionally, two clusters were characterized by the same mutations in splicing genes (i.e., SRSF2 or ZRSR2 with associated TET2), together with co-mutations in genes that correlate with a more aggressive phenotype ( RUNX1, ASXL1, EZH2). These groups showed overlapping myelodysplastic and myeloproliferative features by FAB, intermediate-risk per prognostic scores, and decreased survival as compared with the first 2 clusters. A distinct cluster was characterized by SF3B1 gene mutations (frequently associated with TET2 or DNMT3A). The SF3B1 cluster was marked by a dysplastic phenotype and more severe anemia (with a high rate of transfusion-dependence) with respect to the other genomic groups. Four clusters (27.8% of the overall population) were characterized by the presence of mutations in signal transduction and tyrosine kinase pathways genes ( NRAS/KRAS, CBL, JAK2 and SETBP1 often associated with ASXL1/TET2, respectively). Patients within these groups showed a mostly proliferative phenotype and intermediate survival. The presence of JAK2 mutations appeared to confer a better prognosis compared with other proliferative clusters. A minority of patients (4.6%) were assigned to two high-risk clusters, defined by TP53 mutations (often biallelic) with complex karyotype, and AML-like mutations ( NPM1, IDH1/2, FLT3), respectively. Patients in these clusters were predominantly dysplastic by FAB, had a high prevalence of severe cytopenias, higher blast counts, and a dismal outcome. Finally, we identified a cluster of patients defined by the absence of gene mutations and chromosomal abnormalities (4.5%). The majority of patients assigned to this cluster were characterized by a dysplastic phenotype, had low absolute monocyte counts, and favorable outcome. A subset of patients (18.9%) could not be assigned to a specified cluster in this preliminary analysis. An independent validation of this molecular classification is currently being conducted on a US prospective cohort (n=650). Conclusion. Our results provide proof of concept for a molecular classification of CMML based on the identification of homogeneous genomic signatures and associated clinical features. 15% of patients showed clear overlapping features with other myeloid neoplasms, thus providing the rational to refine the boundaries among different clinical entities. A molecular classification of CMML may provide a basis to define a next-generation prognostic tool and improve clinical decision-making.