Genome wide mapping of ETV6 binding sites in pre-B leukemic cells

ETV6 represses inflammatory response genes and regulates HSPC function during stress hematopoiesis in mice

ETS variant 6 (ETV6) encodes a transcriptional repressor expressed in hematopoietic stem and progenitor cells (HSPCs), where it is required for adult hematopoiesis. Heterozygous pathogenic germline ETV6 variants are associated with thrombocytopenia 5 (T5), a poorly understood genetic condition resulting in thrombocytopenia and predisposition to hematologic malignancies. To elucidate how germline ETV6 variants affect HSPCs and contribute to disease, we generated a mouse model harboring an Etv6R355X loss-of-function variant, equivalent to the T5-associated variant ETV6R359X. Under homeostatic conditions, all HSPC subpopulations are present in the bone marrow (BM) of Etv6R355X/+ mice; however, these animals display shifts in the proportions and/or numbers of progenitor subtypes. To examine whether the Etv6R355X/+ mutation affects HSPC function, we performed serial competitive transplantation and observed that Etv6R355X/+ lineage-sca1+cKit+ (LSK) cells exhibit impaired reconstitution, with near complete failure to repopulate irradiated recipients by the tertiary transplant. Mechanistic studies incorporating cleavage under target and release under nuclease assay, assay for transposase accessible chromatin sequencing, and high-throughput chromosome conformation capture identify ETV6 binding at inflammatory gene loci, including multiple genes within the tumor necrosis factor (TNF) signaling pathway in ETV6-sufficient mouse and human HSPCs. Furthermore, single-cell RNA sequencing of BM cells isolated after transplantation reveals upregulation of inflammatory genes in Etv6R355X/+ progenitors when compared to Etv6+/+ counterparts. Corroborating these findings, Etv6R355X/+ HSPCs produce significantly more TNF than Etv6+/+ cells post-transplantation. We conclude that ETV6 is required to repress inflammatory gene expression in HSPCs under conditions of hematopoietic stress, and this mechanism may be critical to sustain HSPC function.

ETV6 Deficiency Unlocks ERG-Dependent Microsatellite Enhancers to Drive Aberrant Gene Activation in B-Lymphoblastic Leukemia

Distal enhancers play critical roles in sustaining oncogenic gene-expression programs. We identify aberrant enhancer-like activation of GGAA tandem repeats as a characteristic feature of B-cell acute lymphoblastic leukemia (B-ALL) with genetic defects of the ETV6 transcriptional repressor, including ETV6-RUNX1+ and ETV6-null B-ALL. We show that GGAA repeat enhancers are direct activators of previously identified ETV6-RUNX1+/- like B-ALL "signature" genes, including the likely leukemogenic driver EPOR. When restored to ETV6-deficient B-ALL cells, ETV6 directly binds to GGAA repeat enhancers, represses their acetylation, downregulates adjacent genes, and inhibits B-ALL growth. In ETV6-deficient B-ALL cells, we find that the ETS transcription factor ERG directly binds to GGAA microsatellite enhancers and is required for sustained activation of repeat enhancer-activated genes. Together, our findings reveal an epigenetic gatekeeper function of the ETV6 tumor suppressor gene and establish microsatellite enhancers as a key mechanism underlying the unique gene-expression program of ETV6-RUNX1+/- like B-ALL.

SIGNIFICANCE

We find a unifying mechanism underlying a leukemia subtype-defining gene-expression signature that relies on repetitive elements with poor conservation between humans and rodents. The ability of ETV6 to antagonize promiscuous, nonphysiologic ERG activity may shed light on other roles of these key regulators in hematolymphoid development and human disease. See related commentary by Mercher, p. 2. This article is highlighted in the In This Issue feature, p. 1.

Regulome analysis in B-acute lymphoblastic leukemia exposes Core Binding Factor addiction as a therapeutic vulnerability

The ETV6-RUNX1 onco-fusion arises in utero, initiating a clinically silent pre-leukemic state associated with the development of pediatric B-acute lymphoblastic leukemia (B-ALL). We characterize the ETV6-RUNX1 regulome by integrating chromatin immunoprecipitation- and RNA-sequencing and show that ETV6-RUNX1 functions primarily through competition for RUNX1 binding sites and transcriptional repression. In pre-leukemia, this results in ETV6-RUNX1 antagonization of cell cycle regulation by RUNX1 as evidenced by mass cytometry analysis of B-lineage cells derived from ETV6-RUNX1 knock-in human pluripotent stem cells. In frank leukemia, knockdown of RUNX1 or its co-factor CBFβ results in cell death suggesting sustained requirement for RUNX1 activity which is recapitulated by chemical perturbation using an allosteric CBFβ-inhibitor. Strikingly, we show that RUNX1 addiction extends to other genetic subtypes of pediatric B-ALL and also adult disease. Importantly, inhibition of RUNX1 activity spares normal hematopoiesis. Our results suggest that chemical intervention in the RUNX1 program may provide a therapeutic opportunity in ALL.

Identification of new ETV6 modulators through a high-throughput functional screening

ETV6 transcriptional activity is critical for proper blood cell development in the bone marrow. Despite the accumulating body of evidence linking ETV6 malfunction to hematological malignancies, its regulatory network remains unclear. To uncover genes that modulate ETV6 repressive transcriptional activity, we performed a specifically designed, unbiased genome-wide shRNA screen in pre-B acute lymphoblastic leukemia cells. Following an extensive validation process, we identified 13 shRNAs inducing overexpression of ETV6 transcriptional target genes. We showed that the silencing of AKIRIN1, COMMD9, DYRK4, JUNB, and SRP72 led to an abrogation of ETV6 repressive activity. We identified critical modulators of the ETV6 function which could participate in cellular transformation through the ETV6 transcriptional network.

•

We develop a genome-wide shRNAs screen for ETV6 modulators

•

The screen uncovered 13 novel putative ETV6 modulator genes

•

The modulators demonstrated a broad impact on the ETV6 transcriptional network

•

T-ALL cells results suggest modulators are conserved in other cellular contexts

Molecular basis of ETV6-mediated predisposition to childhood acute lymphoblastic leukemia

There is growing evidence supporting an inherited basis for susceptibility to acute lymphoblastic leukemia (ALL) in children. In particular, we and others reported recurrent germline ETV6 variants linked to ALL risk, which collectively represent a novel leukemia predisposition syndrome. To understand the influence of ETV6 variation on ALL pathogenesis, we comprehensively characterized a cohort of 32 childhood leukemia cases arising from this rare syndrome. Of 34 nonsynonymous germline ETV6 variants in ALL, we identified 22 variants with impaired transcription repressor activity, loss of DNA binding, and altered nuclear localization. Missense variants retained dimerization with wild-type ETV6 with potentially dominant-negative effects. Whole-transcriptome and whole-genome sequencing of this cohort of leukemia cases revealed a profound influence of germline ETV6 variants on leukemia transcriptional landscape, with distinct ALL subsets invoking unique patterns of somatic cooperating mutations. 70% of ALL cases with damaging germline ETV6 variants exhibited hyperdiploid karyotype with characteristic recurrent mutations in NRAS, KRAS, and PTPN11. In contrast, the remaining 30% cases had a diploid leukemia genome and an exceedingly high frequency of somatic copy-number loss of PAX5 and ETV6, with a gene expression pattern that strikingly mirrored that of ALL with somatic ETV6-RUNX1 fusion. Two ETV6 germline variants gave rise to both acute myeloid leukemia and ALL, with lineage-specific genetic lesions in the leukemia genomes. ETV6 variants compromise its tumor suppressor activity in vitro with specific molecular targets identified by assay for transposase-accessible chromatin sequencing profiling. ETV6-mediated ALL predisposition exemplifies the intricate interactions between inherited and acquired genomic variations in leukemia pathogenesis.

ETV6 germline mutations cause HDAC3/NCOR2 mislocalization and upregulation of interferon response genes

ETV6 is an ETS family transcription factor that plays a key role in hematopoiesis and megakaryocyte development. Our group and others have identified germline mutations in ETV6 resulting in autosomal dominant thrombocytopenia and predisposition to malignancy; however, molecular mechanisms defining the role of ETV6 in megakaryocyte development have not been well established. Using a combination of molecular, biochemical, and sequencing approaches in patient-derived PBMCs, we demonstrate abnormal cytoplasmic localization of ETV6 and the HDAC3/NCOR2 repressor complex that led to overexpression of HDAC3-regulated interferon response genes. This transcriptional dysregulation was also reflected in patient-derived platelet transcripts and drove aberrant proplatelet formation in megakaryocytes. Our results suggest that aberrant transcription may predispose patients with ETV6 mutations to bone marrow inflammation, dysplasia, and megakaryocyte dysfunction.

Evidence for a role of RUNX1 as recombinase cofactor for TCRβ rearrangements and pathological deletions in ETV6-RUNX1 ALL

T-cell receptor gene beta (TCRβ) gene rearrangement represents a complex, tightly regulated molecular mechanism involving excision, deletion and recombination of DNA during T-cell development. RUNX1, a well-known transcription factor for T-cell differentiation, has recently been described to act in addition as a recombinase cofactor for TCRδ gene rearrangements. In this work we employed a RUNX1 knock-out mouse model and demonstrate by deep TCRβ sequencing, immunostaining and chromatin immunoprecipitation that RUNX1 binds to the initiation site of TCRβ rearrangement and its homozygous inactivation induces severe structural changes of the rearranged TCRβ gene, whereas heterozygous inactivation has almost no impact. To compare the mouse model results to the situation in Acute Lymphoblastic Leukemia (ALL) we analyzed TCRβ gene rearrangements in T-ALL samples harboring heterozygous Runx1 mutations. Comparable to the Runx1+/- mouse model, heterozygous Runx1 mutations in T-ALL patients displayed no detectable impact on TCRβ rearrangements. Furthermore, we reanalyzed published sequence data from recurrent deletion borders of ALL patients carrying an ETV6-RUNX1 translocation. RUNX1 motifs were significantly overrepresented at the deletion ends arguing for a role of RUNX1 in the deletion mechanism. Collectively, our data imply a role of RUNX1 as recombinase cofactor for both physiological and aberrant deletions.