Multilineage involvement in KMT2A-rearranged B acute lymphoblastic leukaemia: cell-of-origin, biology, and clinical implications.

AIMS: B lymphoblastic leukaemia/lymphoma (B-ALL) is thought to originate from Pro/Pre-B cells and the genetic aberrations largely reside in lymphoid-committed cells. A recent study demonstrated that a proportion of paediatric B-ALL patients have BCR::ABL1 fusion in myeloid cells, suggesting a chronic myeloid leukaemia (CML)-like biology in this peculiar subset of B-ALL, although it is not entirely clear if the CD19-negative precursor compartment is a source of the myeloid cells. Moreover, the observation has not yet been extended to other fusion-driven B-ALLs. METHODS AND RESULTS: In this study we investigated a cohort of KMT2A-rearranged B-ALL patients with a comparison to BCR::ABL1-rearranged B-ALL by performing cell sorting via flow cytometry followed by FISH (fluorescence in situ hybridization) analysis on each of the sorted populations. In addition, RNA sequencing was performed on one of the sorted populations. These analyses showed that (1) multilineage involvement was present in 53% of BCR::ABL1 and 36% of KMT2A-rearranged B-ALL regardless of age, (2) multilineage involvement created pitfalls for residual disease monitoring, and (3) HSPC transcriptome signatures were upregulated in KMT2A-rearranged B-ALL with multilineage involvement. CONCLUSIONS: In summary, multilineage involvement is common in both BCR::ABL1-rearranged and KMT2A-rearranged B-ALL, which should be taken into consideration when interpreting the disease burden during the clinical course.

Jak2V617F Reversible Activation Shows Its Essential Requirement in Myeloproliferative Neoplasms.

Gain-of-function mutations activating JAK/STAT signaling are seen in the majority of patients with myeloproliferative neoplasms (MPN), most commonly JAK2V617F. Although clinically approved JAK inhibitors improve symptoms and outcomes in MPNs, remissions are rare, and mutant allele burden does not substantively change with chronic therapy. We hypothesized this is due to limitations of current JAK inhibitors to potently and specifically abrogate mutant JAK2 signaling. We therefore developed a conditionally inducible mouse model allowing for sequential activation, and then inactivation, of Jak2V617F from its endogenous locus using a combined Dre-rox/Cre-lox dual-recombinase system. Jak2V617F deletion abrogates MPN features, induces depletion of mutant-specific hematopoietic stem/progenitor cells, and extends overall survival to an extent not observed with pharmacologic JAK inhibition, including when cooccurring with somatic Tet2 loss. Our data suggest JAK2V617F represents the best therapeutic target in MPNs and demonstrate the therapeutic relevance of a dual-recombinase system to assess mutant-specific oncogenic dependencies in vivo. SIGNIFICANCE: Current JAK inhibitors to treat myeloproliferative neoplasms are ineffective at eradicating mutant cells. We developed an endogenously expressed Jak2V617F dual-recombinase knock-in/knock-out model to investigate Jak2V617F oncogenic reversion in vivo. Jak2V617F deletion abrogates MPN features and depletes disease-sustaining MPN stem cells, suggesting improved Jak2V617F targeting offers the potential for greater therapeutic efficacy. See related commentary by Celik and Challen, p. 701. This article is featured in Selected Articles from This Issue, p. 695.

Interaction between myelodysplasia-related gene mutations and ontogeny in acute myeloid leukemia.

Accurate classification and risk stratification are critical for clinical decision making in patients with acute myeloid leukemia (AML). In the newly proposed World Health Organization and International Consensus classifications of hematolymphoid neoplasms, the presence of myelodysplasia-related (MR) gene mutations is included as 1 of the diagnostic criteria for AML, AML-MR, based largely on the assumption that these mutations are specific for AML with an antecedent myelodysplastic syndrome. ICC also prioritizes MR gene mutations over ontogeny (as defined in the clinical history). Furthermore, European LeukemiaNet (ELN) 2022 stratifies these MR gene mutations into the adverse-risk group. By thoroughly annotating a cohort of 344 newly diagnosed patients with AML treated at the Memorial Sloan Kettering Cancer Center, we show that ontogeny assignments based on the database registry lack accuracy. MR gene mutations are frequently observed in de novo AML. Among the MR gene mutations, only EZH2 and SF3B1 were associated with an inferior outcome in the univariate analysis. In a multivariate analysis, AML ontogeny had independent prognostic values even after adjusting for age, treatment, allo-transplant and genomic classes or ELN risks. Ontogeny also helped stratify the outcome of AML with MR gene mutations. Finally, de novo AML with MR gene mutations did not show an adverse outcome. In summary, our study emphasizes the importance of accurate ontogeny designation in clinical studies, demonstrates the independent prognostic value of AML ontogeny, and questions the current classification and risk stratification of AML with MR gene mutations.

Lenalidomide-associated B-cell ALL: clinical and pathologic correlates and sensitivity to lenalidomide withdrawal.

Lenalidomide is an effective component of induction and maintenance therapy for multiple myeloma, though with a risk of secondary malignancies, including acute lymphoblastic leukemia (ALL). In contrast to therapy-related myeloid neoplasia, lenalidomide-associated lymphoblastic neoplasia remains poorly characterized. We conducted a dual institution retrospective study of 32 ALL cases that arose after lenalidomide maintenance (all B-lineage, 31/32 BCR::ABL-negative). B-cell ALL (B-ALL) was diagnosed at median 54 months (range, 5-119) after first exposure to lenalidomide and after median 42 months of cumulative lenalidomide exposure (range, 2-114). High incidence of TP53 mutations (9/19 evaluable cases) and low hypodiploidy (8/26 patients) were identified. Despite median age of 65 years and poor-risk B-ALL features observed in the cohort, rates of complete response (CR) or CR with incomplete hematologic recovery were high (25/28 patients receiving treatment). Median event-free survival was 35.4 months among treated patients (not reached among those undergoing allogeneic hematopoietic cell transplantation [HCT]). Sixteen patients remain alive without evidence of B-ALL after HCT or extended maintenance therapy. We also describe regression of B-ALL or immature B-cell populations with B-ALL immunophenotype after lenalidomide discontinuation in 5 patients, suggesting lenalidomide may drive leukemic progression even after initiation of lymphoblastic neoplasia and that lenalidomide withdrawal alone may be an appropriate first-line intervention in selected patients. Monitoring for early B-ALL-like proliferations may offer opportunities for lenalidomide withdrawal to prevent progression. Established combination chemotherapy regimens, newer surface antigen-targeted approaches, and allogeneic HCT are effective in many patients with lenalidomide-associated B-ALL and should be offered to medically fit patients.

Molecular predictors of immunophenotypic measurable residual disease clearance in acute myeloid leukemia.

Measurable residual disease (MRD) is a powerful prognostic factor in acute myeloid leukemia (AML). However, pre-treatment molecular predictors of immunophenotypic MRD clearance remain unclear. We analyzed a dataset of 211 patients with pre-treatment next-generation sequencing who received induction chemotherapy and had MRD assessed by serial immunophenotypic monitoring after induction, subsequent therapy, and allogeneic stem cell transplant (allo-SCT). Induction chemotherapy led to MRD- remission, MRD+ remission, and persistent disease in 35%, 27%, and 38% of patients, respectively. With subsequent therapy, 34% of patients with MRD+ and 26% of patients with persistent disease converted to MRD-. Mutations in CEBPA, NRAS, KRAS, and NPM1 predicted high rates of MRD- remission, while mutations in TP53, SF3B1, ASXL1, and RUNX1 and karyotypic abnormalities including inv (3), monosomy 5 or 7 predicted low rates of MRD- remission. Patients with fewer individual clones were more likely to achieve MRD- remission. Among 132 patients who underwent allo-SCT, outcomes were favorable whether patients achieved early MRD- after induction or later MRD- after subsequent therapy prior to allo-SCT. As MRD conversion with chemotherapy prior to allo-SCT is rarely achieved in patients with specific baseline mutational patterns and high clone numbers, upfront inclusion of these patients into clinical trials should be considered.

A Phase I/II Open-Label Study of Molibresib for the Treatment of Relapsed/Refractory Hematologic Malignancies.

PURPOSE: Molibresib is a selective, small molecule inhibitor of the bromodomain and extra-terminal (BET) protein family. This was an open-label, two-part, Phase I/II study investigating molibresib monotherapy for the treatment of hematological malignancies (NCT01943851). PATIENTS AND METHODS: Part 1 (dose escalation) determined the recommended Phase 2 dose (RP2D) of molibresib in patients with acute myeloid leukemia (AML), Non-Hodgkin lymphoma (NHL), or multiple myeloma. Part 2 (dose expansion) investigated the safety and efficacy of molibresib at the RP2D in patients with relapsed/refractory myelodysplastic syndrome (MDS; as well as AML evolved from antecedent MDS) or cutaneous T-cell lymphoma (CTCL). The primary endpoint in Part 1 was safety and the primary endpoint in Part 2 was objective response rate (ORR). RESULTS: There were 111 patients enrolled (87 in Part 1, 24 in Part 2). Molibresib RP2Ds of 75 mg daily (for MDS) and 60 mg daily (for CTCL) were selected. Most common Grade 3+ adverse events included thrombocytopenia (37%), anemia (15%), and febrile neutropenia (15%). Six patients achieved complete responses [3 in Part 1 (2 AML, 1 NHL), 3 in Part 2 (MDS)], and 7 patients achieved partial responses [6 in Part 1 (4 AML, 2 NHL), 1 in Part 2 (MDS)]. The ORRs for Part 1, Part 2, and the total study population were 10% [95% confidence interval (CI), 4.8-18.7], 25% (95% CI, 7.3-52.4), and 13% (95% CI, 6.9-20.6), respectively. CONCLUSIONS: While antitumor activity was observed with molibresib, use was limited by gastrointestinal and thrombocytopenia toxicities. Investigations of molibresib as part of combination regimens may be warranted.

PRC2-Inactivating Mutations in Cancer Enhance Cytotoxic Response to DNMT1-Targeted Therapy via Enhanced Viral Mimicry.

Polycomb repressive complex 2 (PRC2) has oncogenic and tumor-suppressive roles in cancer. There is clinical success of targeting this complex in PRC2-dependent cancers, but an unmet therapeutic need exists in PRC2-loss cancer. PRC2-inactivating mutations are a hallmark feature of high-grade malignant peripheral nerve sheath tumor (MPNST), an aggressive sarcoma with poor prognosis and no effective targeted therapy. Through RNAi screening in MPNST, we found that PRC2 inactivation increases sensitivity to genetic or small-molecule inhibition of DNA methyltransferase 1 (DNMT1), which results in enhanced cytotoxicity and antitumor response. Mechanistically, PRC2 inactivation amplifies DNMT inhibitor-mediated expression of retrotransposons, subsequent viral mimicry response, and robust cell death in part through a protein kinase R (PKR)-dependent double-stranded RNA sensor. Collectively, our observations posit DNA methylation as a safeguard against antitumorigenic cell-fate decisions in PRC2-loss cancer to promote cancer pathogenesis, which can be therapeutically exploited by DNMT1-targeted therapy. SIGNIFICANCE: PRC2 inactivation drives oncogenesis in various cancers, but therapeutically targeting PRC2 loss has remained challenging. Here we show that PRC2-inactivating mutations set up a tumor context-specific liability for therapeutic intervention via DNMT1 inhibitors, which leads to innate immune signaling mediated by sensing of derepressed retrotransposons and accompanied by enhanced cytotoxicity. See related commentary by Guil and Esteller, p. 2020. This article is highlighted in the In This Issue feature, p. 2007.

Arsenic trioxide therapy predisposes to herpes zoster reactivation despite minimally myelosuppressive therapy.

Acute Promyelocytic Leukemia (APL) is a unique subtype of acute myeloid leukemia that is highly responsive to minimally myelosuppressive therapy with all-trans retinoic acid (ATRA) and arsenic trioxide (ATO). We and others have observed a higher than expected incidence of herpes zoster reactivation in APL patients treated with ATO. Memorial Sloan Kettering Cancer Center (MSKCC) has been using ATO since 1997 in all relapsed APL patients, and more recently has included it in our front-line APL regimens. Here we present a retrospective analysis of the factors contributing to herpes zoster reactivation among APL patients.

Somatic Mutations Drive Specific, but Reversible, Epigenetic Heterogeneity States in AML.

Epigenetic allele diversity is linked to inferior prognosis in acute myeloid leukemia (AML). However, the source of epiallele heterogeneity in AML is unknown. Herein we analyzed epiallele diversity in a genetically and clinically annotated AML cohort. Notably, AML driver mutations linked to transcription factors and favorable outcome are associated with epigenetic destabilization in a defined set of susceptible loci. In contrast, AML subtypes linked to inferior prognosis manifest greater abundance and highly stochastic epiallele patterning. We report an epiallele outcome classifier supporting the link between epigenetic diversity and treatment failure. Mouse models with TET2 or IDH2 mutations show that epiallele diversity is especially strongly induced by IDH mutations, precedes transformation to AML, and is enhanced by cooperation between somatic mutations. Furthermore, epiallele complexity was partially reversed by epigenetic therapies in AML driven by TET2/IDH2, suggesting that epigenetic therapy might function in part by reducing population complexity and fitness of AMLs. SIGNIFICANCE: We show for the first time that epigenetic clonality is directly linked to specific mutations and that epigenetic allele diversity precedes and potentially contributes to malignant transformation. Furthermore, epigenetic clonality is reversible with epigenetic therapy agents.This article is highlighted in the In This Issue feature, p. 1775.

Leukemia Cell of Origin Influences Apoptotic Priming and Sensitivity to LSD1 Inhibition.

The cell of origin of oncogenic transformation is a determinant of therapeutic sensitivity, but the mechanisms governing cell-of-origin-driven differences in therapeutic response have not been delineated. Leukemias initiating in hematopoietic stem cells (HSC) are less sensitive to chemotherapy and highly express the transcription factor MECOM (EVI1) compared with leukemias derived from myeloid progenitors. Here, we compared leukemias initiated in either HSCs or myeloid progenitors to reveal a novel function for EVI1 in modulating p53 protein abundance and activity. HSC-derived leukemias exhibit decreased apoptotic priming, attenuated p53 transcriptional output, and resistance to lysine-specific demethylase 1 (LSD1) inhibitors in addition to classical genotoxic stresses. p53 loss of function in Evi1 lo progenitor-derived leukemias induces resistance to LSD1 inhibition, and EVI1hi leukemias are sensitized to LSD1 inhibition by venetoclax. Our findings demonstrate a role for EVI1 in p53 wild-type cancers in reducing p53 function and provide a strategy to circumvent drug resistance in chemoresistant EVI1 hi acute myeloid leukemia. SIGNIFICANCE: We demonstrate that the cell of origin of leukemia initiation influences p53 activity and dictates therapeutic sensitivity to pharmacologic LSD1 inhibitors via the transcription factor EVI1. We show that drug resistance could be overcome in HSC-derived leukemias by combining LSD1 inhibition with venetoclax.See related commentary by Gu et al., p. 1445.This article is highlighted in the In This Issue feature, p. 1426.

Rational Targeting of Cooperating Layers of the Epigenome Yields Enhanced Therapeutic Efficacy against AML.

Disruption of epigenetic regulation is a hallmark of acute myeloid leukemia (AML), but epigenetic therapy is complicated by the complexity of the epigenome. Herein, we developed a long-term primary AML ex vivo platform to determine whether targeting different epigenetic layers with 5-azacytidine and LSD1 inhibitors would yield improved efficacy. This combination was most effective in TET2 mut AML, where it extinguished leukemia stem cells and particularly induced genes with both LSD1-bound enhancers and cytosine-methylated promoters. Functional studies indicated that derepression of genes such as GATA2 contributes to drug efficacy. Mechanistically, combination therapy increased enhancer-promoter looping and chromatin-activating marks at the GATA2 locus. CRISPRi of the LSD1-bound enhancer in patient-derived TET2 mut AML was associated with dampening of therapeutic GATA2 induction. TET2 knockdown in human hematopoietic stem/progenitor cells induced loss of enhancer 5-hydroxymethylation and facilitated LSD1-mediated enhancer inactivation. Our data provide a basis for rational targeting of cooperating aberrant promoter and enhancer epigenetic marks driven by mutant epigenetic modifiers. SIGNIFICANCE: Somatic mutations of genes encoding epigenetic modifiers are a hallmark of AML and potentially disrupt many components of the epigenome. Our study targets two different epigenetic layers at promoters and enhancers that cooperate to aberrant gene silencing, downstream of the actions of a mutant epigenetic regulator.This article is highlighted in the In This Issue feature, p. 813.

Loss of plasmacytoid dendritic cell differentiation is highly predictive for post-induction measurable residual disease and inferior outcomes in acute myeloid leukemia.

Measurable residual disease is associated with inferior outcomes in patients with acute myeloid leukemia (AML). Measurable residual disease monitoring enhances risk stratification and may guide therapeutic intervention. The European LeukemiaNet working party recently came to a consensus recommendation incorporating leukemia associated immunophenotype-based different from normal approach by multi-color flow cytometry for measurable residual disease evaluation. However, the analytical approach is highly expertise-dependent and difficult to standardize. Here we demonstrate that loss of plasmacytoid dendritic cell differentiation after 7+3 induction in AML is highly specific for measurable residual disease positivity (specificity 97.4%) in a uniformly treated patient cohort. Moreover, loss of plasmacytoid dendritic cell differentiation as determined by a blast-to-plasmacytoid dendritic cell ratio >10 was strongly associated with inferior overall and relapse-free survival (RFS) [Hazard ratio 2.79, 95% confidence interval (95%CI): 0.98-7.97; P=0.077) and 3.83 (95%CI: 1.51-9.74; P=0.007), respectively), which is similar in magnitude to measurable residual disease positivity. Importantly, measurable residual disease positive patients who reconstituted plasmacytoid dendritic cell differentiation (blast/ plasmacytoid dendritic cell ratio <10) showed a higher rate of measurable residual disease clearance at later pre-transplant time points compared to patients with loss of plasmacytoid dendritic cell differentiation (blast/ plasmacytoid dendritic cell ratio <10) (6 of 12, 50% vs 2 of 18, 11%; P=0.03). Furthermore pre-transplant plasmacytoid dendritic cell recovery was associated with superior outcome in measurable residual disease positive patients. Our study provides a novel, simple, broadly applicable, and quantitative multi-color flow cytometry approach to risk stratification in AML.

Epigenetic Identity in AML Depends on Disruption of Nonpromoter Regulatory Elements and Is Affected by Antagonistic Effects of Mutations in Epigenetic Modifiers.

We performed cytosine methylation sequencing on genetically diverse patients with acute myeloid leukemia (AML) and found leukemic DNA methylation patterning is primarily driven by nonpromoter regulatory elements and CpG shores. Enhancers displayed stronger differential methylation than promoters, consisting predominantly of hypomethylation. AMLs with dominant hypermethylation featured greater epigenetic disruption of promoters, whereas those with dominant hypomethylation displayed greater disruption of distal and intronic regions. Mutations in IDH and DNMT3A had opposing and mutually exclusive effects on the epigenome. Notably, co-occurrence of both mutations resulted in epigenetic antagonism, with most CpGs affected by either mutation alone no longer affected in double-mutant AMLs. Importantly, this epigenetic antagonism precedes malignant transformation and can be observed in preleukemic LSK cells from Idh2R140Q or Dnmt3aR882H single-mutant and Idh2R140Q/Dnmt3aR882H double-mutant mice. Notably, IDH/DNMT3A double-mutant AMLs manifested upregulation of a RAS signaling signature and displayed unique sensitivity to MEK inhibition ex vivo as compared with AMLs with either single mutation.Significance: AML is biologically heterogeneous with subtypes characterized by specific genetic and epigenetic abnormalities. Comprehensive DNA methylation profiling revealed that differential methylation of nonpromoter regulatory elements is a driver of epigenetic identity, that gene mutations can be context-dependent, and that co-occurrence of mutations in epigenetic modifiers can result in epigenetic antagonism. Cancer Discov; 7(8); 868-83. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 783.

Late-replicating heterochromatin is characterized by decreased cytosine methylation in the human genome.

Heterochromatin is believed to be associated with increased levels of cytosine methylation. With the recent availability of genome-wide, high-resolution molecular data reflecting chromatin organization and methylation, such relationships can be explored systematically. As well-defined surrogates for heterochromatin, we tested the relationship between DNA replication timing and DNase hypersensitivity with cytosine methylation in two human cell types, unexpectedly finding the later-replicating, more heterochromatic regions to be less methylated than early replicating regions. When we integrated gene-expression data into the study, we found that regions of increased gene expression were earlier replicating, as previously identified, and that transcription-targeted cytosine methylation in gene bodies contributes to the positive correlation with early replication. A self-organizing map (SOM) approach was able to identify genomic regions with early replication and increased methylation, but lacking annotated transcripts, loci missed in simple two variable analyses, possibly encoding unrecognized intergenic transcripts. We conclude that the relationship of cytosine methylation with heterochromatin is not simple and depends on whether the genomic context is tandemly repetitive sequences often found near centromeres, which are known to be heterochromatic and methylated, or the remaining majority of the genome, where cytosine methylation is targeted preferentially to the transcriptionally active, euchromatic compartment of the genome.

LINE-1 activity in facultative heterochromatin formation during X chromosome inactivation.

During X chromosome inactivation (XCI), Xist RNA coats and silences one of the two X chromosomes in female cells. Little is known about how XCI spreads across the chromosome, although LINE-1 elements have been proposed to play a role. Here we show that LINEs participate in creating a silent nuclear compartment into which genes become recruited. A subset of young LINE-1 elements, however, is expressed during XCI, rather than being silenced. We demonstrate that such LINE expression requires the specific heterochromatic state induced by Xist. These LINEs often lie within escape-prone regions of the X chromosome, but close to genes that are subject to XCI, and are associated with putative endo-siRNAs. LINEs may thus facilitate XCI at different levels, with silent LINEs participating in assembly of a heterochromatic nuclear compartment induced by Xist, and active LINEs participating in local propagation of XCI into regions that would otherwise be prone to escape.

CG dinucleotide periodicities recognized by the Dnmt3a-Dnmt3L complex are distinctive at retroelements and imprinted domains.

The Dnmt3a and Dnmt3L genes are critical mediators of cytosine methylation during gametogenesis, with major actions noted at transposable elements and imprinted loci. The Dnmt3a-Dnmt3L complex was recently described to have preferential activity at CG dinucleotides located 8-10 bp apart. Because cytosine methylation is heterogeneously distributed in the genome, we tested whether this relative sequence preference explains the effects of mutation of the Dnmt3a and Dnmt3L genes using bioinformatic analysis. We found that the human and mouse genomes are significantly enriched in a CG dinucleotide periodicity of 2 bp, leading to an increased frequency of CGs spaced 8 bp apart that represent widespread targets for this protein complex. When we broke down the human and mouse genomes by annotation, we found that this significant 2-bp periodicity and increased 8-bp periodicity are maintained in Alu SINEs in both species. The 8-bp periodicity was mapped genome-wide, identifying enrichment at the promoters of both paternally and maternally methylated imprinted genes and at CG dinucleotide-enriched sequences. We conclude that CG dinucleotide periodicity helps to explain some but not all of the relative sequence specificity of mutations of Dnmt3a or Dnmt3L in the establishment of germline cytosine methylation patterns.

High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers.

Many genome-wide assays involve the generation of a subset (or representation) of the genome following restriction enzyme digestion. The use of enzymes sensitive to cytosine methylation allows high-throughput analysis of this epigenetic regulatory process. We show that the use of a dual-adapter approach allows us to generate genomic representations that includes fragments of <200 bp in size, previously not possible when using the standard approach of using a single adapter. By expanding the representation to smaller fragments using HpaII or MspI, we increase the representation by these isoschizomers to more than 1.32 million loci in the human genome, representing 98.5% of CpG islands and 91.1% of refSeq promoters. This advance allows the development of a new, high-resolution version of our HpaII-tiny fragment Enrichment by Ligation-mediated PCR (HELP) assay to study cytosine methylation. We also show that the MspI representation generates information about copy-number variation, that the assay can be used on as little as 10 ng of DNA and that massively parallel sequencing can be used as an alternative to microarrays to read the output of the assay, making this a powerful discovery platform for studies of genomic and epigenomic abnormalities.

CG dinucleotide clustering is a species-specific property of the genome.

Cytosines at cytosine-guanine (CG) dinucleotides are the near-exclusive target of DNA methyltransferases in mammalian genomes. Spontaneous deamination of methylcytosine to thymine makes methylated cytosines unusually susceptible to mutation and consequent depletion. The loci where CG dinucleotides remain relatively enriched, presumably due to their unmethylated status during the germ cell cycle, have been referred to as CpG islands. Currently, CpG islands are solely defined by base compositional criteria, allowing annotation of any sequenced genome. Using a novel bioinformatic approach, we show that CG clusters can be identified as an inherent property of genomic sequence without imposing a base compositional a priori assumption. We also show that the CG clusters co-localize in the human genome with hypomethylated loci and annotated transcription start sites to a greater extent than annotations produced by prior CpG island definitions. Moreover, this new approach allows CG clusters to be identified in a species-specific manner, revealing a degree of orthologous conservation that is not revealed by current base compositional approaches. Finally, our approach is able to identify methylating genomes (such as Takifugu rubripes) that lack CG clustering entirely, in which it is inappropriate to annotate CpG islands or CG clusters.