Single cell characterization of B-lymphoid differentiation and leukemic cell states during chemotherapy in ETV6-RUNX1-positive pediatric leukemia identifies drug-targetable transcription factor activities

Background Tight regulatory loops orchestrate commitment to B cell fate within bone marrow. Genetic lesions in this gene regulatory network underlie the emergence of the most common childhood cancer, acute lymphoblastic leukemia (ALL). The initial genetic hits, including the common translocation that fuses ETV6 and RUNX1 genes, lead to arrested cell differentiation. Here, we aimed to characterize transcription factor activities along the B-lineage differentiation trajectory as a reference to characterize the aberrant cell states present in leukemic bone marrow, and to identify those transcription factors that maintain cancer-specific cell states for more precise therapeutic intervention. Methods We compared normal B-lineage differentiation and in vivo leukemic cell states using single cell RNA-sequencing (scRNA-seq) and several complementary genomics profiles. Based on statistical tools for scRNA-seq, we benchmarked a workflow to resolve transcription factor activities and gene expression distribution changes in healthy bone marrow lymphoid cell states. We compared these to ALL bone marrow at diagnosis and in vivo during chemotherapy, focusing on leukemias carrying the ETV6-RUNX1 fusion. Results We show that lymphoid cell transcription factor activities uncovered from bone marrow scRNA-seq have high correspondence with independent ATAC- and ChIP-seq data. Using this comprehensive reference for regulatory factors coordinating B-lineage differentiation, our analysis of ETV6-RUNX1-positive ALL cases revealed elevated activity of multiple ETS-transcription factors in leukemic cells states, including the leukemia genome-wide association study hit ELK3. The accompanying gene expression changes associated with natural killer cell inactivation and depletion in the leukemic immune microenvironment. Moreover, our results suggest that the abundance of G1 cell cycle state at diagnosis and lack of differentiation-associated regulatory network changes during induction chemotherapy represent features of chemoresistance. To target the leukemic regulatory program and thereby overcome treatment resistance, we show that inhibition of ETS-transcription factors reduced cell viability and resolved pathways contributing to this using scRNA-seq. Conclusions Our data provide a detailed picture of the transcription factor activities characterizing both normal B-lineage differentiation and those acquired in leukemic bone marrow and provide a rational basis for new treatment strategies targeting the immune microenvironment and the active regulatory network in leukemia

TRIM16 Acts as an E3 Ubiquitin Ligase and Can Heterodimerize with Other TRIM Family Members

The TRIM family of proteins is distinguished by its tripartite motif (TRIM). Typically, TRIM proteins contain a RING finger domain, one or two B-box domains, a coiled-coil domain and the more variable C-terminal domains. TRIM16 does not have a RING domain but does harbour two B-box domains. Here we showed that TRIM16 homodimerized through its coiled-coil domain and heterodimerized with other TRIM family members; TRIM24, Promyelocytic leukaemia (PML) protein and Midline-1 (MID1). Although, TRIM16 has no classic RING domain, three-dimensional modelling of TRIM16 suggested that its B-box domains adopts RING-like folds leading to the hypothesis that TRIM16 acts as an ubiquitin ligase. Consistent with this hypothesis, we demonstrated that TRIM16, devoid of a classical RING domain had auto-polyubiquitination activity and acted as an E3 ubiquitin ligase in vivo and in vitro assays. Thus via its unique structure, TRIM16 possesses both heterodimerization function with other TRIM proteins and also has E3 ubiquitin ligase activity

CDK‐mediated activation of the SCFFBXO28 ubiquitin ligase promotes MYC‐driven transcription and tumourigenesis and predicts poor survival in breast cancer

SCF (Skp1/Cul1/F‐box) ubiquitin ligases act as master regulators of cellular homeostasis by targeting key proteins for ubiquitylation. Here, we identified a hitherto uncharacterized F‐box protein, FBXO28 that controls MYC‐dependent transcription by non‐proteolytic ubiquitylation. SCFFBXO28 activity and stability are regulated during the cell cycle by CDK1/2‐mediated phosphorylation of FBXO28, which is required for its efficient ubiquitylation of MYC and downsteam enhancement of the MYC pathway. Depletion of FBXO28 or overexpression of an F‐box mutant unable to support MYC ubiquitylation results in an impairment of MYC‐driven transcription, transformation and tumourigenesis. Finally, in human breast cancer, high FBXO28 expression and phosphorylation are strong and independent predictors of poor outcome. In conclusion, our data suggest that SCFFBXO28 plays an important role in transmitting CDK activity to MYC function during the cell cycle, emphasizing the CDK‐FBXO28‐MYC axis as a potential molecular drug target in MYC‐driven cancers, including breast cancer

IFN gamma directly counteracts imatinib-induced apoptosis of primary human CD34+CML stem/progenitor cells potentially through the upregulation of multiple key survival factors

Tyrosine kinase inhibitors (TKIs) have dramatically improved the survival in chronic myeloid leukemia (CML), but residual disease typically persists even after prolonged treatment. Several lines of evidence suggest that TKIs administered to CML patients upregulate interferon gamma (IFN gamma) production, which may counteract the anti-tumorigenic effects of the therapy. We now show that activated T cell-conditioned medium (TCM) enhanced proliferation and counteracted imatinib-induced apoptosis of CML cells, and addition of a neutralizing anti-IFN gamma antibody at least partially inhibited the anti-apoptotic effect. Likewise, recombinant IFN gamma also reduced imatinib-induced apoptosis of CML cells. This anti-apoptotic effect of IFN gamma was independent of alternative IFN gamma signaling pathways, but could be notably diminished by STAT1-knockdown. Furthermore, IFN gamma upregulated the expression of several anti-apoptotic proteins, including MCL1, PARP9, and PARP14, both in untreated and imatinib-treated primary human CD34+ CML stem/progenitor cells. Our results suggest that activated T cells in imatinib-treated CML patients can directly rescue CML cells from imatinib-induced apoptosis at least partially through the secretion of IFN gamma, which exerts a rapid, STAT1-dependent anti-apoptotic effect potentially through the simultaneous upregulation of several key hematopoietic survival factors. These mechanisms may have a major clinical impact, when targeting residual leukemic stem/progenitor cells in CML.Funding Agencies|Cancer Research Funds of Radiumhemmet [174283]; Thorsmans stiftelse for preleukemiforskning; Cathrine Everts Research Foundation; Lars Hierta Memorial Foundation; Swedish Research Council [2013-08807]; Karolinska Institute Foundation and Funds; Gunnar Grimfors Gavofond for Hematologisk Forskning; Emil Andersson Fund for Medical Research</p

B-boxes are required for TRIM16’s E3 ligase activity.

<p>(A) TRIM16 <i>in vivo</i> polyubiquitination assay. In HEK293 cells, HA-Ub was co-transfected with various TRIM16-GFP domain deletion expression plasmids. The protein lysate was subjected to immunoprecipitation by GFP antibody, and subjected to Western blot and probed with anti-HA antibody for Ub (right panel) and anti-GFP antibody for TRIM16 (left and middle panel). Two exposure times are shown. GFP antibodies detect both un-ubiquitinated and polyubiquitinated forms of TRIM16. Polyubiquitinated smear is present in the sample transfected with wild-type TRIM16 and shown by anti-GFP and anti-HA antibodies. (B) <i>In vitro</i> ubiquitination assay with myc-His tagged TRIM16 together with a panel of E2 enzymes, showing activity with the UbcH5 family. (C) <i>In vitro</i> ubiquitination assay with full-length TRIM16, TRIM16 domain deletion mutants or empty vector showing extensive polyubiquitination with full-length TRIM16 as detected by Western blot with anti-myc antibodies. Numbers indicate protein size in kDa. (D) Recombinant TRIM16 (Abnova) was evaluated for E3 activity in the presence of recombinant E1, UbcH5b, and HA-Ub as indicated. The capacity to catalyse auto-ubiquitination <i>in vitro</i> was observed only in the presence of ZnCl₂ and ATP. Western Blot (lower panel) with TRIM16 antibody showed amount of TRIM16 protein in each lane.</p

Amino acid sequence comparison of TRIM16 and MID1 and modeling of TRIM16 B-boxes.

<p>The conserved residues in the zinc binding regions are in bold underlined type and are; cysteine; C, histidine; H, alanine; A aspartic acid; D. (<b>A</b>) TRIM16 and MID1 share the zinc binding consensus sequence for B-box1. (B) TRIM16 and MID1 share the Zinc binding consensus sequence for B-box2. (C/D) Modeling of B-boxes reveals zinc binding capability. Superimposition of the alpha-carbon backbone of the B-boxes from the MDM1 NMR structure (purple) and the homology model of TRIM16 (blue-grey). These two structures overlay with an average root-mean-square deviation of 0.4 Å.</p

TRIM16 homodimerizes through its coiled-coil domain

<p>. (A) TRIM16-GFP domain deletion plasmids. (B) Co-transfection of TRIM16-GFP and TRIM16-myc-His in HEK293 cells and subsequent immunoprecipitation by anti-myc antibody (Ab) and Western blot with anti-GFP antibody. Whole cell extract (WCE) used as total input. (C) TRIM16 homodimerizes through its coiled-coil domain. GFP deletion mutants were co-transfected with the TRIM16-myc-His vector. Anti-GFP antibody was used to pull down proteins binding the GFP tagged proteins and the TRIM16-myc-His was used to detect self-association via its different tag (right panel). Transfection efficiency was confirmed (left panel). TRIM16-GFP mutants were efficiently pulled down (middle panel). * non-specific bands, # refer to text.</p

TRIM16 can heterodimerize with MID1, TRIM24 and PML.

<p>(A) Schematic structures of TRIM proteins used in heterodimerization studies. (B) TRIM16 binds MID1. Co-transfection of MID1-GFP and TRIM16-myc-His in HEK293 cells and subsequent immunoprecipitation by anti-myc antibody and Western blot with anti-GFP antibody. (C) MID1 was pulled down via its GFP antibody and a Western blot was performed to detect TRIM16-myc-His in the immunoprecipitated protein complex. (D) Whole cell lysates (WCL) of HEK293 cells transfected with empty vector (EV) or TRIM16-GFP were immunopreciptated with anti-GFP antibody. An anti-PML antibody was used to detect PML as a binding partner of TRIM16. (E) Lysates containing both TRIM16-GFP and TRIM24-His proteins were immunopreciptated with anti-GFP antibody. Anti-His antibody was used to detect the presence of TRIM24 in the TRIM16-associated complex.</p

Single cell characterization of B-lymphoid differentiation and leukemic cell states during chemotherapy in ETV6-RUNX1-positive pediatric leukemia identifies drug-targetable transcription factor activities

Background: Tight regulatory loops orchestrate commitment to B cell fate within bone marrow. Genetic lesions in this gene regulatory network underlie the emergence of the most common childhood cancer, acute lymphoblastic leukemia (ALL). The initial genetic hits, including the common translocation that fuses ETV6 and RUNX1 genes, lead to arrested cell differentiation. Here, we aimed to characterize transcription factor activities along the B-lineage differentiation trajectory as a reference to characterize the aberrant cell states present in leukemic bone marrow, and to identify those transcription factors that maintain cancer-specific cell states for more precise therapeutic intervention. Methods: We compared normal B-lineage differentiation and in vivo leukemic cell states using single cell RNA-sequencing (scRNA-seq) and several complementary genomics profiles. Based on statistical tools for scRNA-seq, we benchmarked a workflow to resolve transcription factor activities and gene expression distribution changes in healthy bone marrow lymphoid cell states. We compared these to ALL bone marrow at diagnosis and in vivo during chemotherapy, focusing on leukemias carrying the ETV6-RUNX1 fusion. Results: We show that lymphoid cell transcription factor activities uncovered from bone marrow scRNA-seq have high correspondence with independent ATAC- and ChIP-seq data. Using this comprehensive reference for regulatory factors coordinating B-lineage differentiation, our analysis of ETV6-RUNX1-positive ALL cases revealed elevated activity of multiple ETS-transcription factors in leukemic cells states, including the leukemia genome-wide association study hit ELK3. The accompanying gene expression changes associated with natural killer cell inactivation and depletion in the leukemic immune microenvironment. Moreover, our results suggest that the abundance of G1 cell cycle state at diagnosis and lack of differentiation-associated regulatory network changes during induction chemotherapy represent features of chemoresistance. To target the leukemic regulatory program and thereby overcome treatment resistance, we show that inhibition of ETS-transcription factors reduced cell viability and resolved pathways contributing to this using scRNA-seq. Conclusions: Our data provide a detailed picture of the transcription factor activities characterizing both normal B-lineage differentiation and those acquired in leukemic bone marrow and provide a rational basis for new treatment strategies targeting the immune microenvironment and the active regulatory network in leukemia.publishedVersionPeer reviewe

Combination of tyrosine kinase inhibitors and the MCL1 inhibitor S63845 exerts synergistic antitumorigenic effects on CML cells

Tyrosine kinase inhibitor (TKI) treatment has dramatically improved the survival of chronic myeloid leukemia (CML) patients, but measurable residual disease typically persists. To more effectively eradicate leukemia cells, simultaneous targeting of BCR-ABL1 and additional CML-related survival proteins has been proposed. Notably, several highly specific myeloid cell leukemia 1 (MCL1) inhibitors have recently entered clinical trials for various hematologic malignancies, although not for CML, reflecting the insensitivity of CML cell lines to single MCL1 inhibition. Here, we show that combining TKI (imatinib, nilotinib, dasatinib, or asciminib) treatment with the small-molecule MCL1 inhibitor S63845 exerted strong synergistic antiviability and proapoptotic effects on CML lines and CD34+ stem/progenitor cells isolated from untreated CML patients in chronic phase. Using wild-type BCR-ABL1-harboring CML lines and their T315I-mutated sublines (generated by CRISPR/Cas9-mediated homologous recombination), we prove that the synergistic proapoptotic effect of the drug combination depended on TKI-mediated BCR-ABL1 inhibition, but not on TKI-related off-target mechanisms. Moreover, we demonstrate that colony formation of CML but not normal hematopoietic stem/progenitor cells became markedly reduced upon combination treatment compared to imatinib monotherapy. Our results suggest that dual targeting of MCL1 and BCR-ABL1 activity may efficiently eradicate residual CML cells without affecting normal hematopoietic stem/progenitors.Funding Agencies|Cancer Research Funds of Radiumhemmet [174283]; Thorsmans Stiftelse for Pre-leukemi Forskning; Gunnar Grimfors Gavofond for Hematologisk Forskning; Cathrine Everts Research Foundation; Lars Hierta Memorial Foundation; Emil Anderson Fund for Medical Research; Swedish Research CouncilSwedish Research CouncilEuropean Commission [2013-08807]; Karolinska Institute Foundation and Funds; Karolinska InstituteKarolinska Institutet</p