Cheng, Chia‐Ho

Abstract Despite advances in understanding the genetic abnormalities in myeloproliferative neoplasms (MPN) and the development of JAK2 inhibitors, there is an urgent need to devise new treatment strategies, particularly for patients with triple-negative (TN) myelofibrosis (MF) who lack mutations in the JAK2 kinase pathway and have very poor clinical outcomes. Here we report that MYC copy number gain and increased MYC expression frequently occur in TN-MF and that MYC-directed activation of S100A9, an alarmin protein that plays pivotal roles in inflammation and innate immunity, is necessary and sufficient to drive development and progression of MF. Notably, the MYC-S100A9 circuit provokes a complex network of inflammatory signaling that involves numerous hematopoietic cell types in the bone marrow microenvironment. Accordingly, genetic ablation of S100A9 or treatment with small molecules targeting the MYC-S100A9 pathway effectively ameliorates MF phenotypes, highlighting the MYC–alarmin axis as a novel therapeutic vulnerability for this subgroup of MPNs. Significance: This study establishes that MYC expression is increased in TN-MPNs via trisomy 8, that a MYC-S100A9 circuit manifest in these cases is sufficient to provoke myelofibrosis and inflammation in diverse hematopoietic cell types in the BM niche, and that the MYC-S100A9 circuit is targetable in TN-MPNs.

Primary Myelofibrosis (PMF) is an aggressive myeloproliferative neoplasm (MPN) characterized by constitutional symptoms, cytopenias, splenomegaly, extramedullary hematopoiesis, bone marrow (BM) fibrosis, and a propensity to transform to acute myeloid leukemia (AML). Despite advances in understanding the underlying genetic abnormalities in MPN and the development of JAK2 inhibitors to treat MF, there is an urgent need to devise new treatment strategies, particularly for triple negative myelofibrosis (TN-MF) cases that lack mutations in the JAK2 kinase pathway and who have inferior outcomes. To identify potential oncogenic drivers in TN-MF, we performed cytogenetic analyses and targeted exome sequencing of 98 genes commonly mutated in myeloid malignancies in 584 MF patients identified from the Total Cancer Care database at Moffitt Cancer Center. Notably, there were no significant differences in the somatic mutation profiles (other than JAK2 activating mutations) between TN vs. JAK2/CALR/MPL mutant MF patients. However, trisomy of chromosome 8 occurs more frequently in TN-MF vs. other subtypes (26.7% vs. 6.2%, p=0.0001). To identify oncogenic drivers in trisomy 8+ TN-MF, we performed scRNA-seq analysis from 18,651 cells ( n=15,850 from 3 normal donors and n=2,801 from a trisomy 8+ TN-MF patient) and identified a total of 1,260 genes that are differentially regulated in hematopoietic stem cells (HSCs) of trisomy 8+ TN-MF compared to normal donors. Among these, 131 genes were located on chromosome 8 and MYC was one of the top 5 genes that are significantly upregulated in trisomy 8+ HSCs. Additional immunohistochemistry staining demonstrated that MYC protein levels in BM cells were significantly higher in trisomy 8+ TN-MF patients vs. trisomy 8- patients. To assess whether MYC might be a driver of MPN, we established an Mx1-Cre +/-;Rosa26 LSL-MYC/LSL-MYC transgenic mouse model that inducibly overexpresses MYC in HSCs following treatment of polyinosinic-polycytidylic acid (pIpC). Forced MYC expression induced profound anemia, monocytosis, megakaryocytic atypia, splenomegaly, increased BM collagen deposition/fibrosis, extramedullary hematopoiesis in spleen and liver, and significantly reduced OS vs. Mx1-Cre +/- wild type controls (median OS 258 days vs. not reached (NR), p<0.0001). MYC also promoted expansion of HSCs, myeloid progenitors, and Gr-1 +/CD11b + mature myeloid cells in BM and spleen; thus, MYC was sufficient to provoke MF-like disease independent of JAK2 pathway mutations. Subsequent scRNA-seq analysis of 25,232 cells ( n=13,552 from control and n=11,680 from a MYC transgenic mouse) revealed that S100a9 mRNA levels were increased by MYC in most major hematopoietic cell types. Supporting oncogenic roles of S100A9 in MYC-driven MF, deletion of S100a9 in Mx1-Cre +/-;Rosa26 LSL-MYC/LSL-MYC;S100a9 -/- mice abolished MYC-induced MF phenotypes and significantly improved OS (median OS NR vs. 225 days, p=0.0345). To test whether inhibition of MYC can suppress MF disease progression in vivo, lethally irradiated CD45.1 +/CD45.2 + WT mice were transplanted with BM cells harvested from CD45.2 + Mx1-Cre +/-;Rosa26 LSL-MYC/LSL-MYC mice that were treated with pIpC 20 weeks prior to transplant. MYCi975, a small molecule that inhibits MYC, reduced MYC and S100a9 protein levels in BM and spleen cells. Further, MYCi975 effectively suppressed MYC-driven MF phenotypes and disease progression, and significantly improved OS (median OS NR vs. 252 days, p=0.0109). In summary, our studies are the first to describe an oncogenic role of MYC in MF pathogenesis, where MYC provokes an alarmin-driven inflammatory circuit, and where the MYC-S100A9 axis represents a therapeutic vulnerability for TN-MF patients. Accordingly, our results provide a strong rationale for testing agents targeting MYC or S100A9 in early phase clinical trials in MPNs having increased MYC levels or activity.

MYC gene amplification and somatic mutations are frequent in both adult and pediatric AML although how MYC drives and contributes to the development and maintenance of AML has not been resolved. Transcription factor EB (TFEB) is a master regulator of genes that control autophagy and lysosome biogenesis, a central catabolic recycling pathway that regulates cell survival. Given the oncogenic effects of MYC in AML and that the induction of autophagy compromises AML cell growth and survival, we tested if the oncogenic effect of MYC depends on its suppression of TFEB transcription programs in AML. In support of this hypothesis, inducible MYC expression in K562 and THP-1 leukemia cells was sufficient to suppress expression of TFEB and its target genes. Further and conversely, MYC knockdown in NB4 AML cells provoked increased expression of TFEB mRNA and protein, as well as increased expression of TFEB target genes. Notably, dose response studies demonstrated that expression of TFEBS211A, constitutively nuclear form of TFEB that is refractory to control by mTORC1 signaling, dramatically impairs proliferation of HL60, OCI-AML2 and OCI-AML3 AML cells. In addition, induction of TFEBS211A provoked rampant apoptosis. Of important, overexpression of TFEBS211A in HL-60 and OCI-AML3 cells was also sufficient to promote monocytic and granulocytic differentiation, as judged by morphological changes and the acquisition of mature monocytic and granulocytic markers including CD11b, Gr1, and CD15. To identify TFEB targets that might contribute to myeloid/granulocytic differentiation, we performed RNA-seq analysis of HL60 leukemia cells engineered to inducibly express the TFEBS211A transgene. Using a cut-off of fold change>4 with q<0.01, a total of 1152 genes were differentially regulated following the induction of TFEBS211A in HL60 cells. As expected, this included the robust induction of nearly all TFEB target genes associated with the autophagy-lysosome pathway, but also of STAT1, KLF4, KLF6, CEBPB, CSF1 and GATA2 genes that are necessary and/or sufficient to provoke terminal monocytic and granulocytic differentiation of AML cells. Surprisingly, among genes induced by TFEB in HL60 cells is IDH1, which catalyzes the production of α-ketoglutarate (α-KG), a required substrate of the TET family of dioxygenases that convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). In particular, TFEBS211A expression provoked increases in levels of α-KG and significant increases in global levels of 5hmC in genomic DNA of HL60 leukemia cells both ex vivo and in vivo. Furthermore, TFEB-mediated induction of IDH1/2 mRNA and protein, and of IDH1 promoter activity, was antagonized by inducible expression of MYC. To assess the global effects of TFEB on 5mC/5hmC landscapes, we performed paired reduced representation bisulfite (BS)- and oxidation bisulfite (oxBS)-sequencing. As predicted, TFEBS211A induced both loss and gains of 5mC, but there were more losses (n = 722) than gains (n = 459) across all 22 chromosomes. Quite remarkably, and consistent with TFEB-provoked increases of 5hmC signals, TFEBS211A exclusively induced 5hmC gains (in a total of 863 genes), and 37% and 36% of these 5hmC gains occurred in promoter regions and CpG islands, and across all 22 chromosomes. Comparison of BS- and oxBS-seq versus RNA-seq analyses of HL60 cells expressing TFEB revealed significant changes in mRNA levels and concomitant differential changes in 5mC and 5hmC marks in KLF4, KLF6, STAT3, TP73, andFOXO1 that have pivotal roles in controlling myeloid cell differentiation and death. Collectively, these findings demonstrate that MYC suppresses TFEB expression and function in AML cells, and that TFEB functions as a tumor suppressor that provokes AML cell differentiation and death. Strikingly, these responses rely on epigenetic control, where TFEB directly induces the transcription of IDH1 and IDH2 to provoke global hydroxylation of 5-methylcytosine and the expression of genes that drive terminal differentiation and apoptosis. Thus, a MYC-TFEB-IDH1/2-TET2 circuit controls AML cell fate. Disclosures Murphy: Merck: Research Funding; Puma Biotech: Research Funding. Ballabio:CASMA Therapeutics: Other: Co-Founder. Kaufmann:Takeda Pharmaceuticals: Research Funding.

Abstract MYC oncoproteins regulate transcription of genes directing cell proliferation, metabolism, and tumorigenesis. A variety of alterations drive MYC expression in acute myeloid leukemia (AML), and enforced MYC expression in hematopoietic progenitors is sufficient to induce AML. Here we report that AML and myeloid progenitor cell growth and survival rely on MYC-directed suppression of Transcription Factor EB (TFEB), a master regulator of the autophagy–lysosome pathway. Notably, although originally identified as an oncogene, TFEB functions as a tumor suppressor in AML, where it provokes AML cell differentiation and death. These responses reflect TFEB control of myeloid epigenetic programs by inducing expression of isocitrate dehydrogenase-1 (IDH1) and IDH2, resulting in global hydroxylation of 5-methycytosine. Finally, activating the TFEB–IDH1/IDH2–TET2 axis is revealed as a targetable vulnerability in AML. Thus, epigenetic control by an MYC–TFEB circuit dictates myeloid cell fate and is essential for maintenance of AML. Significance: Alterations in epigenetic control are a hallmark of AML. This study establishes that a MYC–TFEB circuit controls AML differentiation and epigenetic programs by inducing IDH1/IDH2 and hydroxylation of 5-methylcytosine, that TFEB functions as a tumor suppressor in AML, and that this circuit is a targetable vulnerability in AML. See related commentary by Wu and Eisenman, p. 116.

Abstract Motivation Testing SNP–SNP interactions is considered as a key for overcoming bottlenecks of genetic association studies. However, related statistical methods for testing SNP–SNP interactions are underdeveloped. Results We propose the SNP Interaction Pattern Identifier (SIPI), which tests 45 biologically meaningful interaction patterns for a binary outcome. SIPI takes non-hierarchical models, inheritance modes and mode coding direction into consideration. The simulation results show that SIPI has higher power than MDR (Multifactor Dimensionality Reduction), AA_Full, Geno_Full (full interaction model with additive or genotypic mode) and SNPassoc in detecting interactions. Applying SIPI to the prostate cancer PRACTICAL consortium data with approximately 21 000 patients, the four SNP pairs in EGFR-EGFR, EGFR-MMP16 and EGFR-CSF1 were found to be associated with prostate cancer aggressiveness with the exact or similar pattern in the discovery and validation sets. A similar match for external validation of SNP–SNP interaction studies is suggested. We demonstrated that SIPI not only searches for more meaningful interaction patterns but can also overcome the unstable nature of interaction patterns. Availability and Implementation The SIPI software is freely available at http://publichealth.lsuhsc.edu/LinSoftware/. Supplementary information Supplementary data are available at Bioinformatics online.