Abstract
Long-range oncogenic enhancers play an important role in cancer. Yet, whether similar regulation of tumor suppressor genes is relevant remains unclear. Loss of expression of PTEN is associated with the pathogenesis of various cancers, including T-cell acute lymphoblastic leukemia (T-ALL). Here, we identify a highly conserved distal enhancer (PE) that interacts with the PTEN promoter in multiple hematopoietic populations, including T cells, and acts as a hub of relevant transcription factors in T-ALL. Consistently, loss of PE leads to reduced PTEN levels in T-ALL cells. Moreover, PE-null mice show reduced Pten levels in thymocytes and accelerated development of NOTCH1-induced T-ALL. Furthermore, secondary loss of PE in established leukemias leads to accelerated progression and a gene expression signature driven by Pten loss. Finally, we uncovered recurrent deletions encompassing PE in T-ALL, which are associated with decreased PTEN levels. Altogether, our results identify PE as the first long-range tumor suppressor enhancer directly implicated in cancer.
Here, we identify a PTEN enhancer that is recurrently deleted in patients with T-ALL. Loss of this enhancer leads to reduced PTEN levels in T cells together with accelerated generation and progression of NOTCH1-induced leukemia in vivo. These results uncover long-range regulation of tumor suppressor genes as a relevant mechanism in cancer.
This article is highlighted in the In This Issue feature, p. 1
Introduction
Enhancers are distal DNA regulatory regions that act in cis, and in an orientation-independent manner, to regulate gene expression through three-dimensional physical interaction with gene promoters (1–4). In recent years, enhancer malfunction has been shown to play a key role in cancer (5). Specifically, cancer cells can use enhancers via epigenetic activation or via noncoding mutations that create enhancers de novo, increasing the expression of a variety of oncogenes. In addition, we and others have described numerous tissue-specific oncogenic enhancers as recurrently duplicated in cancer (5, 6). In this setting, MYC is the best example to date of an oncogene regulated by a plethora of enhancer regions, all of which can be dysregulated in different cancer types (7). However, whether long-range enhancer regions might control the expression of tumor suppressor genes (TSG), or whether these might also play a role in cancer development, remains largely unknown. T-cell acute lymphoblastic leukemia (T-ALL) is a NOTCH1-driven aggressive hematologic malignancy that presents in both pediatric and adult patients and requires treatment with intensified chemotherapy (8). Importantly, 15% to 20% of patients with T-ALL show loss of expression of the PTEN TSG (9, 10), which leads to resistance to anti-NOTCH1 therapies with γ-secretase inhibitors (GSI) in T-ALL in vivo (11). However, not all of the T-ALL cases with loss of PTEN expression can be explained by prototypical mutations or deletions in its coding region, by defective splicing, or by epigenetic silencing (10), highlighting our incomplete understanding of the mechanisms mediating loss of expression of PTEN. In this context, we hypothesized that PTEN expression in T-ALL might be controlled by yet unidentified enhancer regions and that loss of these putative enhancers might help explain the lack of PTEN expression in certain T-ALL cases. To test this hypothesis, we integrated chromosome conformation and epigenetic profiling analyses and identified a bona fide enhancer of PTEN, whose functionality we further dissected using CRISPR/Cas9-mutant T-ALL cell lines and novel mouse models in vivo.
Results
Identification of PE, a PTEN Enhancer in T-ALL
To reveal the global pattern of chromosomal interactions of the PTEN promoter in the leukemia genome and identify potential regions that might act as enhancers of PTEN, we performed chromosome conformation capture 4C-seq experiments (12, 13) in three independent T-ALL cell lines (DND41, HPB-ALL, and JURKAT) using the PTEN promoter as the viewpoint. These analyses revealed that the PTEN promoter interacts at high frequency with several regions along most of its topologically associating domain (TAD; ref. 14; Fig. 1A, upper tracks). Among these, we identified a region approximately 550 kb downstream of PTEN (hereby named PE, for PTEN enhancer), located in an intronic sequence of the neighboring RNLS gene, that shows high levels of interaction with the PTEN promoter. Moreover, epigenetic profiling analyses revealed concomitant bona fide enhancer marks in this region, including high levels of H3K27ac and H3K4me1, together with binding of CTCF, BRD4, and ZNF143 (Fig. 1A, middle tracks). Importantly, the CTCF-binding sites in the PTEN promoter and the PE enhancer are located in convergent orientation (Fig. 1A, CTCF track), suggesting that cohesin-mediated loops might be responsible for this interaction (15–18). In addition, GRO-seq analyses of T-ALL cells (19) uncovered bidirectional transcription from this region (Fig. 1A, bottom tracks), which has also been suggested to mark active enhancer regulatory elements (20, 21). Finally, analyses of publicly available H3K27ac Hi-C chromatin immunoprecipitation (Hi-ChIP) data in CUTLL1 T-ALL cells (19) revealed that, among the different interacting regions with the PTEN promoter, the interaction with the PE enhancer is the only one detected at FDR < 1E-15 (Fig. 1A, top). Indeed, the PTEN–PE interaction is detected at an extremely low FDR value < 1E-49 (Fig. 1A, top), overall suggesting a prominent role for PE in the regulation of PTEN.
Following up on these findings, we next performed reciprocal 4C-seq analyses using this putative PE enhancer region as the viewpoint. Notably, these analyses confirmed its interaction with the PTEN promoter in all three cell lines (Fig. 1A, upper tracks) regardless of PTEN mutational status, as JURKAT cells (PTEN null; ref. 9) show similar levels of interaction as compared with HPB-ALL and DND41 cells (PTEN positive; ref. 9). This is also in line with chromatin immunoprecipitation sequencing (ChIP-seq) data showing H3K27ac in the PE enhancer both in PTEN-positive and PTEN-null cell lines (Supplementary Fig. S1A). Further inspection of ENCODE data (22) in this area in DND41 T-ALL cells revealed that the PE region is an approximately 10 kb superenhancer with higher H3K27ac levels along an approximately 4 kb core sequence. However, reduced levels of H3K27ac are seen in the GM12878 B-lymphoblastoid cells or primary CD20+ B cells (Supplementary Fig. S1B), whereas greatly decreased or completely absent levels of H3K27ac were observed in a variety of solid tumor cell lines, primary CD14+ monocytes, or other primary cells derived from a variety of nonhematopoietic tissues (Supplementary Fig. S1B). Moreover, H3K27ac was also absent in PTEN-positive prostate (23) and thyroid (24) cancer cells (Supplementary Fig. S1C), which is of particular relevance given the well-described role of PTEN loss in the development of both prostate and thyroid cancers (25–27). Overall, these data suggest that PE might be a hematopoietic-specific enhancer prominently active in the lymphoid compartment.
Next, to investigate what transcription factors might regulate PE, we mined publicly available ChIP-seq data in T-ALL cells and identified a wide range of key T-cell transcription factors that bind to the PE enhancer, including MED1, UTX, NOTCH1, RBPJ, ETS1, RUNX1, MYC, HOXA9, MYB, LMO1, TLX1, TCF12, and CDK7 (Fig. 1B). Of note, the cobinding of NOTCH1 and CDK7 to enhancer regions has been previously reported in T-ALL (19). Indeed, some of these NOTCH1/CDK7-bound enhancers are GSI sensitive, whereas the rest of them are GSI insensitive but sensitive to CDK7 inhibition (19). In the case of PE, analyses of CUTLL1 T-ALL cells treated with either a GSI or the THZ1 CDK7 inhibitor revealed lower H3K27ac levels only after THZ1 treatment, highlighting that PE likely belongs to the GSI-insensitive but CDK7-sensitive group of enhancers (Supplementary Fig. S1D). On the other hand, we detected minimal or no binding of TAL1 or TCF3, and we also failed to detect significant binding of GATA3 binding (Fig. 1B), which we previously showed to be a critical pioneering factor for the N-Me enhancer (1), underscoring that other pioneering factors might be more relevant in the case of PE. Finally, we also performed epigenetic profiling in human primary T-ALLs. These analyses revealed a robust H3K27ac mark (Fig. 1C) together with open chromatin via Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq; Fig. 1D). Interestingly, even if both marks show variable intensity in the PE region, they are present in all of the samples tested and correlate with the H3K27ac and ATAC-seq signal observed in the PTEN promoter (Fig. 1C and D). Consistently, inspection of Hi-C data from human normal T cells and primary T-ALLs (19) revealed a clear interaction between the PTEN promoter and PE in both normal and transformed T cells (Supplementary Fig. S1E and S1F). Overall, these results suggest that PE is a very relevant enhancer hub in human T-ALL controlled by inputs from a plethora of important T-cell transcription factors and highlight a relevant role for CDK7 in maintaining its activity.
Notably, multispecies DNA sequence alignment of the PE sequence revealed remarkable conservation of this region in mice and other mammals, with its central sequence also being conserved in reptiles but not in birds, fish, or amphibians (Supplementary Fig. S2). Thus, we then performed 4C-seq analyses in NOTCH1-induced mouse primary T-ALLs (either driven by a potent ΔE-NOTCH1 construct or a weaker HDΔP-NOTCH1 construct) to uncover the global pattern of chromosomal interactions of the mouse Pten promoter in the leukemia genome. These experiments revealed that the mouse Pten promoter shows more promiscuous interactions inside its own TAD as compared with the human PTEN promoter (Fig. 1E). However, and most notably, among these contacts, we observed again high levels of interaction with the PE homologous region in mice, which also shows high levels of H3K27ac (Fig. 1E). In addition, reciprocal 4C-seq analyses using the mouse PE enhancer as the viewpoint confirmed its interaction with the Pten promoter. Finally, and similar to the results observed in human cells, ChIP-seq in mouse Th1 cells showed CTCF binding in both the Pten promoter and PE enhancer with CTCF-binding sites in convergent orientation, whereas ATAC-seq analyses from mouse T-ALL cells revealed that the PE enhancer is enriched for open chromatin in mouse leukemias (Fig. 1E). Overall, our results identify the PE enhancer as a long-range regulatory region that might be involved in the control of PTEN in T-ALL.
Functional Characterization of the PE Enhancer
To functionally characterize the potential role of this regulatory region in the control of PTEN expression, we first performed reporter luciferase assays to test the transactivation capability of PE. These experiments showed strong, orientation-independent activation of reporter constructs containing the core ∼4 kb sequence of both human and mouse PE (Fig. 2A). Next, we hypothesized that PE loss in T-ALL cells might lead to reduced levels of PTEN. To test this, we generated DND41 cells with CRISPR/Cas9-induced heterozygous loss of PE (Fig. 2B). Notably, qRT-PCR analyses showed reduced PTEN levels in two independent clones with heterozygous loss of PE (Fig. 2C). Consistently, intracellular FACS staining revealed reduced levels of PTEN protein expression in both clones (Fig. 2D). This effect was further confirmed by Western blot analysis, which showed reduced levels of PTEN together with a concomitant activation of AKT (Fig. 2E), consistent with the well-known role of PTEN as a negative regulator of the PI3K pathway (28). Finally, and in line with the well-described role of PTEN as a TSG, both clones of DND41 cells with heterozygous loss of PE displayed a proliferative advantage as compared with nontargeted DND41 control cells (Fig. 2F). Together, these results demonstrate that PE behaves as a bona fide enhancer of PTEN.
Germline Loss of PE Leads to Reduced Pten Levels during Thymocyte Development
Next, and to test the relevance of the PE enhancer in T-cell development and transformation, we genetically engineered mice with PE deletion using CRISPR/Cas9. To this end, C57BL/6J mouse zygotes were injected with Cas9 protein complexed with synthetic single-guide RNAs (sgRNA) that targeted two sites flanking the ∼4 kb core PE enhancer region with high levels of H3K27ac in the mouse (Supplementary Fig. S3A). Importantly, we obtained three different founders that showed specific deletion of the PE enhancer, with only 1 to 3 bp differences in the size of the deleted fragment (Supplementary Fig. S3A and S3B). PE-null mice show normal postnatal development without any apparent abnormality in behavior or body size, are fertile, and are born with the expected Mendelian ratios (Supplementary Table S1). However, analyses of thymi from approximately 6-week-old mice revealed a marked reduction in Pten mRNA and protein expression levels in PE-knockout mice from the three independent strains as compared with littermate controls (Fig. 3A–C), whereas Pten expression remained unchanged across a comprehensive list of nonthymic tissues (Supplementary Fig. S4A). Of note, PTEN levels seem to be partially reduced in PE heterozygous mice, suggesting a dose-dependent effect for the PE enhancer (Fig. 3A–C). In addition, we did not detect any decrease in Rnls mRNA levels in these same samples, suggesting that the deletion of the PE enhancer does not affect Rnls expression, which is also consistent with the lack of interaction between the PE enhancer and the Rnls promoter in our previous 4C-seq analyses (Fig. 1E; Supplementary Fig. S4B). Finally, immunohistochemistry (IHC) analyses of these thymi revealed a generally preserved architecture and confirmed reduced PTEN levels in PE heterozygous mice, which were further and drastically reduced in PE-knockout mice (Fig. 3D and E).
Following up on these results, we next performed detailed immunophenotypic analyses of thymic populations in approximately 6-week-old mice with different PE allelic dosages. These experiments revealed that mice lacking PE show normal thymus weight and display a mostly normal thymocyte development across the different strains, with no consistent changes in CD4−CD8− double-negative (DN), CD4+CD8+ double-positive (DP), or single-positive (SP) CD4+ and CD8+ cells as compared with wild-type (WT) control mice (Fig. 4A–C; Supplementary Fig. S5A and S5B), and further analyses also failed to show noticeable differences in the spleens from these mice (Supplementary Fig. S5C and S5D). Of note, this is consistent with the phenotypes previously reported in T-cell–specific Pten-knockout mice (29). Given the drastic reduction in PTEN levels observed in total thymus lysates from PE-knockout mice (Fig. 3), we next performed single-cell RNA sequencing (scRNA-seq) analyses in thymi from PE WT or knockout mice to dissect the specific stage/s of T-cell development in which the loss of PE might result in reduced levels of Pten. These analyses revealed that Pten levels in PE-knockout mice are significantly reduced at the CD4−CD8− DN, CD8+ immature SP (ISP), and CD4+CD8+ DP stages of thymocyte development (Fig. 4D–F). Moreover, a statistically nonsignificant trend of decreased Pten levels was also observed in CD4+ and CD8+ SP mature cells (Fig. 4D–F). Interestingly, Rnls levels were not detectable by scRNA-seq in the thymus, and other neighboring genes contained in the TADs either upstream or downstream of the Pten TAD showed no differences or were also undetectable by scRNA-seq (Supplementary Fig. S6), underscoring that PE might be primarily controlling the expression of Pten. We next investigated PE chromatin status throughout mouse hematopoietic development by mining ATAC-seq data available from the ImmGen Consortium (30). These data showed open chromatin already present in long-term hematopoietic stem cells, which is maintained throughout the development of different bone marrow and thymocyte populations (Supplementary Fig. S7A). Interestingly, even if the ATAC-seq signal is variable between DN–ISP–DP–SP stages, we detected significant differences in Pten levels in DN–ISP–DP cells only (Fig. 4F). In addition, analyses of tagHi-C data from different early mouse hematopoietic populations (31) confirmed the interaction of the Pten promoter and PE (Supplementary Fig. S7B). Similarly, epigenetic profiling of different human hematologic and thymic populations along the T-cell differentiation trajectory showed open chromatin and H3K27ac mark in PE across the different populations (Supplementary Fig. S8A and S8B; ref. 32). Consistently, analyses of publicly available capture-HiC data in human hematologic populations (33) showed strong interaction between the PTEN promoter and the PE enhancer in fetal thymus, total CD4 cells, total CD8 cells, total B cells, and megakaryocytes, with lower interactions detected in neutrophils, erythroblasts, monocytes, or endothelial precursors (Supplementary Fig. S9).
Overall, our results demonstrate that deletion of the PE enhancer leads to reduced expression of PTEN throughout thymic T-cell development in mice and highlight that the PE enhancer might also contribute to regulate PTEN expression in other hematopoietic populations.
Loss of PE Results in Acceleration of NOTCH1-Induced T-ALL
Given the fundamental importance of NOTCH1 in T-ALL pathogenesis (34), where approximately 60% of patients with T-ALL show activating mutations of NOTCH1 (8), and given that approximately 20% of patients with T-ALL show loss of PTEN, which might co-occur with NOTCH1 mutations (9), we hypothesized that reduced levels of PTEN due to loss of the PE enhancer might accelerate NOTCH1-induced leukemogenesis. To test this hypothesis, we retrovirally expressed a constitutively active oncogenic mutant form of the NOTCH1 receptor (lacking its extracellular portion but still tethered to the membrane and requiring γ-secretase cleavage) with concomitant expression of the GFP marker (ΔE-NOTCH1-GFP; ref. 6) in bone marrow progenitors obtained from PE WT, heterozygous, or homozygous knockout mice, followed by transplantation into lethally irradiated recipient mice (Fig. 5A). Enforced expression of constitutively active NOTCH1 leads to the exclusive development of T-cell lineage leukemia with CD4+CD8+ DP features accompanied by acute leukemic infiltration of the spleen, in line with the critical role of NOTCH1 in T-cell specification (34, 35). Notably, mice typically develop a CD4+CD8+ DP wave of preleukemic cells in peripheral blood only 3 weeks after transplantation of NOTCH1-infected cells, and this preleukemic DP wave is known to be a good predictor of mouse survival (6). Indeed, our results showed a significantly higher number of DP cells in mice transplanted with either PE heterozygous or PE homozygous knockout NOTCH1–infected bone marrow progenitor cells as compared with mice transplanted with PE WT progenitor cells (Supplementary Fig. S10). Consistently, although mice transplanted with ΔE-NOTCH1–expressing bone marrow progenitor cells from PE WT mice developed T-ALL in 60 to 80 days (Fig. 5B), mice transplanted with ΔE-NOTCH1–expressing bone marrow progenitor cells from mice lacking one copy of PE developed T-ALL with significantly faster kinetics, which was further accelerated by biallelic deletion of PE (strain #6, 20% lower median survival; strain #15, 13% lower median survival; strain #18, 19% lower median survival), underscoring a dose-dependent tumorsuppressor role for PE in T-cell leukemia generation (Fig. 5B).
Next, to test the effects of PE loss in leukemia progression and maintenance, we generated PE conditional knockout mice (Supplementary Fig. S11A and S11B). These mice were subsequently bred with a mouse strain harboring a tamoxifen-inducible Cre from the Rosa26 locus (6). In this context, we generated PE conditional homozygous knockout leukemias following a similar strategy (Fig. 6A). Transplant of these already generated leukemias in a secondary cohort of mice led to rapid disease progression; however, tamoxifen-induced isogenic deletion of PE still resulted in a significant acceleration (32% median lower survival) of T-ALL progression (Fig. 6B). Of note, this acceleration is in line with what we previously observed and reported upon secondary complete loss of Pten in already established leukemias, which led to a 35% lower median survival (PTEN WT = 20 days; PTEN- knockout = 13 days; ref. 11). Analyses of these leukemic cells after tamoxifen-induced isogenic deletion of PE revealed markedly reduced Pten mRNA and protein levels (Fig. 6C and D). Furthermore, global gene expression profiling analyses upon inactivation of PE in these leukemias unveiled Pten as one of the most significantly downregulated genes (Fig. 6E; Supplementary Fig. S12A), whereas Rnls levels were undetectable, and other genes located in the TADs neighboring the Pten TAD did not show significant differences either (Supplementary Fig. S12B). In addition, gene set enrichment analyses (GSEA; ref. 36) against our previously published gene expression signature obtained upon isogenic deletion of Pten in T-ALL using Pten conditional knockout leukemias (11) revealed that the transcriptional program elicited by deletion of PE is driven by Pten loss (Fig. 6F), and GSEA against the C2 database revealed 342 gene sets significantly downregulated (nominal P < 0.01) upon loss of PE (Supplementary Table S2). Altogether, our results demonstrate that loss of the PE enhancer in vivo results in an acceleration of T-ALL generation and progression characterized by reduced expression levels of Pten.
Recurrent Deletions Encompassing PE in Human T-ALL
Finally, to test whether genetic lesions in PE might contribute to T-ALL pathogenesis in human patients, we analyzed a panel of 398 T-ALL cases via single-nucleotide polymorphism (SNP) arrays (n = 398), whole-genome sequencing (WGS; n = 73), and whole-transcriptome sequencing (WTS; n = 360). These analyses uncovered five cases (1.25%) with focal deletions encompassing PE without affecting the PTEN coding sequence (CDS; Fig. 6G). Interestingly, two of these samples show homozygous deletions of the PE enhancer (Supplementary Fig. S13A). However, four out of five of these cases show additional simultaneous deletions targeting the CDS of PTEN independently (Supplementary Fig. S13A). In addition, three out of five of these cases show mutations in NOTCH1 or FBXW7 (Supplementary Fig. S13A). As expected, cases with heterozygous deletions in the PTEN CDS show lower expression levels of PTEN as compared with cases without deletions (Fig. 6H). In addition, further reduced levels are seen when these deletions affect specifically the transcriptional start site (TSS) of PTEN (Fig. 6H). Among cases with loss of PE, there is only one sample that does not concomitantly present with deletions in the PTEN CDS, therefore precluding any relevant comparison in this group. The other four cases with PE loss present with heterozygous deletions affecting the PTEN TSS. However, these four cases show markedly and significantly reduced PTEN levels when compared with other samples similarly having heterozygous deletions in the PTEN TSS but without concomitant loss of PE (Fig. 6H). Furthermore, we performed comparative analyses in a cohort of 1,415 B-cell acute lymphoblastic leukemia (B-ALL) cases. Interestingly, unlike in T-ALL, there was no case with focal deletion of PE among these 1,415 B-ALLs (Fig. 6I), underscoring the relevance of PE in T-ALL. Lastly, we also analyzed the methylation status of the PE DNA sequence in T-ALL. We found four CpG sites in the PE DNA sequence, two of which are hypomethylated and two of which are hypermethylated in normal T cells (Supplementary Fig. S13B). However, both T-ALL cell lines and primary T-ALL cells maintain the hypermethylation in one of these CpG sites, whereas the other three are hypomethylated (Supplementary Fig. S13C and S13D). Overall, these results support that the loss of the PE enhancer does contribute to reduce the expression levels of PTEN specifically in T-ALL and suggest that DNA hypermethylation is likely not involved in the dysregulation of the PE enhancer in T-ALL.
Discussion
PTEN, a key negative regulator of the PI3K/AKT pathway, is the second most frequently mutated TSG (37) and is deregulated in a wide range of cancer types (38). Specifically, loss of PTEN expression is seen in 15% to 20% of patients with T-ALL (9, 10), where it leads to secondary resistance to NOTCH1 inhibitory therapies in vivo (11). In addition, several studies suggest that genetic lesions driving loss of PTEN expression are associated with poor prognosis and outcome in both adult and pediatric patients with T-ALL (39), highlighting its relevance in this hematologic malignancy. However, and most notably, approximately 5% of patients with T-ALL show loss of PTEN expression that cannot be currently explained by the typical mutations or deletions of its CDS (10), underscoring that the mechanisms leading to loss of PTEN expression in cancer remain incompletely understood. Here, we identify a previously unknown long-range enhancer region of PTEN located in a highly conserved intronic region of the neighboring RNLS gene, approximately 550 kb downstream of its TSS. We further prove that this enhancer region, named PE for PTEN enhancer, acts as a hub for a plethora of major T-cell transcription factors, shows bona fide enhancer marks and open chromatin in human T-ALL, physically interacts with the PTEN promoter in normal and transformed T cells, and drives increased PTEN transactivation in reporter assays. In addition, we demonstrate that CRISPR-Cas9–induced deletion of PE in human T-ALL cells leads to reduced PTEN levels and a proliferative advantage. Still, even if PE seems to play a prominent role in the control of PTEN expression, the PTEN promoter interacts with additional regions along the PTEN TAD that also show several enhancer features, suggesting that some of these interacting regions might also play relevant roles in the regulation of PTEN.
To analyze PE effects in T-cell development and transformation in vivo, we generated mice with germinal or conditional deletion of PE. These novel mouse models unveiled that loss of PE is associated with reduced expression of Pten in developing thymocytes, whereas a wide variety of other solid tissues obtained from PE-knockout mice failed to show any difference in Pten levels, suggesting that PE might behave as a hematopoietic-specific enhancer prominently active in the lymphoid compartment. In addition, thymocyte scRNA-seq analyses demonstrated that Pten levels in PE-knockout mice are reduced in CD4−CD8− DN, CD8+ ISP, and CD4+CD8+ DP cells; however, this reduction does not lead to increased cellularity or altered T-cell development in the thymus. Moreover, in line with the majority of cells in the thymus being CD4+CD8+ DP cells, and consistent with the striking difference in PTEN protein levels observed by IHC, scRNA-seq analyses revealed that CD4+CD8+ DP thymocytes express the highest levels of Pten across all thymocyte populations, and it is in this specific population where the most significant decrease in Pten is observed in PE-knockout thymocytes. Still, even if PE clearly plays a prominent role in T cells and T-ALL, epigenetic and chromosome conformation analyses via ATAC-seq, ChIP-seq, and Hi-C in earlier and/or non-T hematopoietic populations in both mice and humans highlight that PE might play a relevant role in the control of PTEN expression beyond the T-cell lineage, as PE shows H3K27ac mark and open chromatin in multiple hematologic populations, and this is also accompanied by physical interaction with the PTEN promoter in non-T populations. Overall, these data suggest that the activity of the PE enhancer might be very prominently controlled by transcription factors that could be specifically present/absent in different stages of development, in addition to (or rather than) by how open its chromatin is. Subsequent follow-up studies are therefore warranted to thoroughly dissect its effects in each specific population and lineage. In addition, it is also worth mentioning that multiple reports have shown abnormalities in different hematopoietic lineages upon loss of PTEN using either an Mx-Cre allele (40–42), an Scl-CreERT allele (39), or a CD4-Cre allele (43, 44), whereas we failed to detect noticeable abnormalities in T-cell development in PE-knockout mice. However, all of these studies investigated the effects of complete loss of PTEN (i.e., using a Ptenflox/flox background, not a Ptenflox/+ background). In our case, even after complete loss of PE, we still detect PTEN expression even if at lower levels, which is what would be expected from the loss of an enhancer region that modulates the expression of its target gene. Thus, the different levels of reduction in PTEN expression might help explain the discordant phenotypes observed. In addition, specific loss of PTEN in T cells using an Lck-Cre allele showed no differences in T-cell development at 6 weeks of age (29), similar to what we observed in PE-knockout mice. However, even if PE-null mice do not show complete loss of PTEN, we cannot rule out at this point that PE loss might result in T-cell abnormalities at later time points beyond 6 weeks of age or might also lead to some developmental defects in earlier stages of hematopoietic development and/or other hematologic lineages beyond T cells.
Regarding its effects in cancer development, it is well- known that PTEN acts as a haploinsufficient TSG, and even minor changes in its expression can lead to measurable effects in cancer development and progression (38). Consistently, germinal loss of one copy of the PE enhancer led to a significant acceleration in the development of NOTCH1-induced T-ALL in vivo, which was further accelerated by deletion of both allelic copies of PE. Furthermore, secondary deletion of PE in mice harboring already established PE conditional knockout T-ALL led to accelerated disease progression in vivo, reduced levels of Pten expression, and a gene expression signature driven by Pten loss. Overall, our results underscore a strong dose-dependent effect of the PE enhancer in T-cell leukemia generation and progression.
We have previously demonstrated the importance of long-range regulatory sequences in T-ALL pathogenesis in vivo, as MYC expression is exquisitely dependent on the NOTCH1-driven N-Me MYC enhancer in T cells (1, 6). Indeed, enhancer deregulation is now considered to be a common mechanism contributing to cancer development (5). However, most of the cases described to date are related to oncogenic enhancers (5–7, 45–47). Recent studies have suggested that a regulatory region of the tumor suppressor PAX5 might be affected in B-cell leukemia (48). Similarly, deleterious mutations in an enhancer region of ACTRT1 have been described in inherited and sporadic basal cell carcinoma (49), and a noncoding regulatory variant in an enhancer region of PER2 has recently been shown to affect the growth of acute myeloid leukemia cells (50). Finally, and more specifically related to PTEN, a different study suggested that p53 might regulate PTEN expression through a 5′ proximal regulatory region in epithelial cells. However, this region falls in the promoter of the neighboring CFL1P1 gene and also seems to affect its expression, as well as the expression of other neighboring genes (51). Moreover, whether this region physically interacts with the PTEN promoter or if it has any relevance in vivo remains unknown (51). In this context, to the best of our knowledge, our study is the first demonstration of the existence of an enhancer region of a major TSG and its relevance in cancer development in vivo. A fascinating question that still needs to be addressed is how the PE enhancer might be developmentally regulated, including what pioneering factor/s controls its chromatin opening and what transcription factors control its effects on PTEN transactivation, both positively or negatively. PTEN regulation is notoriously complex, involving both transcriptional and posttranscriptional mechanisms (38). Key recent findings in the regulation of PTEN include the discovery of alternative translation start sites leading to the expression of PTEN-β, which negatively regulates rDNA transcription and ribosomal biogenesis in the nucleus (52), and PTEN-Long, which is secreted and can exert tumor-suppressive effects in a paracrine manner (53). However, the transcriptional regulation of PTEN is still poorly understood, with only a handful of transcription factors that have been suggested to transcriptionally activate or repress PTEN in certain contexts (54). Follow-up studies to dissect the mechanisms controlling the activity of PE will shed light on the dynamic processes and factors involved in the activation of this enhancer region.
Finally, analyses of a panel of 398 T-ALL cases uncovered 5 cases with deletions encompassing PE. Of note, T-ALL cases with loss of PE show reduced PTEN levels, underscoring its relevance in the pathophysiology of human T-ALL. Although this is a relatively small number of cases (1.25%), it is in a similar range to what has been previously described for copy-number variations affecting other enhancer regions such as the MYC N-Me enhancer (5% of cases; ref. 7), the MYC BENC enhancer (3% of cases; ref. 7), the MYC-LASE enhancer (2% of cases; ref. 7), or the MYC-ECSE enhancer (4% of cases; ref. 7). In addition, no case with loss of PE was detected among a cohort of 1,415 B-ALL cases, highlighting the specific relevance of PE in T-ALL. Loss of PTEN expression is associated with resistance to NOTCH1 inhibition (9, 11) and to glucocorticoids (55), together with increased risk of relapse and long-term adverse outcomes in patients with T-ALL (56, 57). Thus, PE status could be clinically used in the short-term as an additional biomarker to refine risk stratification and to select risk-adapted therapeutic strategies in patients. In line with this, it is also worth noting that Mcm2-deficient mice have been previously shown to develop T-ALL harboring genetic lesions affecting Pten, and WGS of eight Mcm2-deficient T-ALLs uncovered one case (12.5%) harboring a deletion encompassing PE (58), which might have facilitated leukemia development. Still, additional mechanisms beyond genetic deletions may also lead to impaired activity of the PE enhancer. Although we did not find evidence that DNA hypermethylation might play a role in the dysregulation of this enhancer, mutations in the PE sequence or SNPs affecting the binding of important transcription factors controlling the activity of PE might conceivably result in decreased levels of PTEN. Similarly, epigenetic silencing of the PE enhancer resulting in closed chromatin along this region might also occur. In addition, chromosome conformation changes leading to alterations in interactions within TADs, as well as their boundaries, have also recently been shown to affect enhancer activity (19); thus, another potential mechanism driving impaired PE activity might involve the establishment of a TAD boundary gain separating the PTEN promoter from the PE enhancer. It is also tempting to speculate that additional regions yet to be identified might control PTEN expression in other types of cancer. Further studies to elucidate how widespread PE activity might be impaired in T-ALL, and whether PTEN is regulated by additional enhancer regions in other tumor types, are therefore warranted. Overall, our study clearly demonstrates that impaired activity of the PE is a novel mechanism that cooperates in leukemia development and, more broadly, paves the way for the discovery and dissection of enhancers of other TSGs in additional malignancies.
Methods
Cell Culture
We performed cell culture of cell lines in standard conditions in a humidified atmosphere at 37°C under 5% CO2. DND41 (#ACC-525; RRID:CVCL_2022), JURKAT (#ACC-282; RRID:CVCL_0065), and HPB-ALL (#ACC-483; RRID:CVCL_1820) cells were obtained from Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ). We purchased HEK293T cells from the ATCC (#CRL-3216; RRID:CVCL_0063). No cell authentication was performed. T-ALL cells were cultured in HyClone RPMI 1640 Media (SH3002701, Fisher Scientific) supplemented with 10% FBS (900-108, Gemini Bio-Products) and 1% penicillin–streptomycin (45000-652, VWR).
Generation of DND41 Cells with CRISPR/Cas9-Induced Loss of PE
DND41 PE+/− cells were generated following a previously described protocol (59). Briefly, custom primers (Sigma-Aldrich; Supplementary Table S3) were used to generate sgRNA DNA templates. Evaluation for potential off targets was performed through CRISPRscan (https://www.crisprscan.org). In vitro transcription of sgRNAs from templates was performed following the protocol guidelines. Note that 500 ng of each purified sgRNA was independently incubated with 500 ng of Cas9 (1074181, IDT) for 20 minutes at room temperature. Subsequently, 1 million DND41 cells were electroporated with 1 μg of each complex per transfection using the Neon Transfection System (MPK5000, Thermo Fisher Scientific) and 100 μL electroporating tips (MPK10025, Thermo Fisher Scientific). Electroporation conditions were as follows: pulse voltage 1350 V, pulse width 10 ms, pulse number 3, and cell density 107/mL. Single-cell clones were sorted using the Influx High Speed Sorter (BD Biosciences), subsequently grown in tissue culture, and screened for PE loss by PCR using REDTaq ReadyMix (R2523-100RXN, Sigma-Aldrich) following standard conditions and using custom primers (Supplementary Table S3).
Genomic Analysis of Primary T-ALL and B-ALL Samples
Genomic profiling of patient samples was performed by SNP array (Affymetrix SNP 6.0: T-ALL, n = 398; B-ALL, n = 1415), WGS (n = 73), and WTS (n = 360). Sequencing libraries for WGS were prepared from genomic DNA sheared on LE220 ultrasonicator (Covaris) with HyperPrep Library Preparation Kits (Kapa Biosystems). WTS libraries were generated by using TruSeq Stranded Total RNA Library Prep Kit (Illumina). Libraries were analyzed for insert size distribution with a 2100 BioAnalyzer High Sensitivity kit (Agilent Technologies) or Caliper LabChip GX DNA High Sensitivity Reagent kit (PerkinElmer), and were quantified by using the Quant-iT PicoGreen dsDNA assay (Life Technologies). Paired-end sequencing was performed with HiSeq2000/2500 or NovaSeq 6000 (Illumina).
The detection of copy-number variations on PTEN CDS and PE regions was performed by using SNP array data as previously described (60). For WGS, reads were mapped to the GRCh37 reference human genome assembly using BWA, sorted, and deduplicated using Samtools (https://www.htslib.org/). Multiple softwares on Genomon Project (https://github.com/Genomon-Project) were applied to call single-nucleotide variants (SNV) and indels (GenomonFisher, GenomonMutationFilter, and EBFilter) or structural variants (GenomonSV). Expression of PTEN was calculated from gene-level read counts generated by HTseq-count (61). After normalization with DESeq2 Bioconductor R package (62), a regularized log-transformed (rlog) value was calculated by variance-stabilizing transformation of DESeq2. All statistical analyses were performed using R (v3.6.0) software (www.r-project.org).
Statistical significance on PTEN expression was assessed by t test (ggpubr R package; https://rpkgs.datanovia.com/ggpubr/), and values were considered statistically significant at P < 0.05. Associations between categorical values were examined using the Fisher exact test. Data are available at the St. Jude PeCan Data Portal (ref. 63; https://pecan.stjude.cloud/), the Therapeutically Applicable Research to Generate Effective Treatment (TARGET) data portal (https://ocg.cancer.gov/programs/target), the database of Genotypes and Phenotypes (dbGaP), and the European Genome-phenome Archive, as previously described (64, 65).
DNA Methylation Analysis
DNA methylation status of PE was determined using the Illumina MethylationEPIC BeadChip 850K microarray (66). Briefly, 600 ng of DNA of the studied samples (T-cell normal donor samples, T-ALL cell lines, and primary malignancies) was used to hybridize the BeadChip and scanned using HiScan SQ system (Illumina). Raw signal intensity data were initially QC'd and preprocessed from resulting idat files in R statistical environment (v3.6.2) using minfi Bioconductor package (v1.32.0). A number of quality control steps were applied to minimize errors and remove erratic probe signals, such as failed probes (detection P value > 0.01), cross-reacting probes, and probes that overlapped SNVs within ±1 bp of CpG sites, followed by background correction and dye-based normalization using ssNoob algorithm (single-sample normal-exponential out-of-band). DNA methylation scores for each CpG were represented as β values ranging between standard 0 and 1, where 1 represents fully methylated CpGs and 0 represents fully unmethylated. All downstream analyses were performed under R statistical environment (v3.6.1). Unsupervised hierarchical clustering analysis using Euclidean distance with complete-linkage clustering and heatmaps of methylation CpG β values in the studied samples were used to represent the DNA methylation data. Raw DNA methylation data are available upon request.
ChIP-seq and GRO-seq Analysis
We analyzed PE occupancy of chromatin marks and epigenetic and transcription factors using the following T-ALL publicly available ChIP-seq datasets from Gene Expression Omnibus (GEO): GSE29611, GSE54379, GSE54380, GSE29600, GSE72300, GSE59657, GSE58406, GSE83777, GSE17954, GSE85524, GSE83777, GSE109653, GSE94391, GSE62264, GSE29180, GSE50622, GSE76783, GSE115893, GSE60482, and GSE104992. We analyzed nascent RNA around the PE locus using the following T-ALL publicly available GRO-seq data set from GEO: GSE115894.
ChIPmentation Analysis
Fifty to 100,000 CD45+ sorted T-ALL cells were fixed using 1% formaldehyde (Thermo Fisher Scientific 28906) and quenched by glycin (125 mmol/L). Next, cell pellets were lysed in 100 μL of short-term complete lysis buffer (50 mmol/L Tris-HCl, pH 8.0, 10 mmol/L EDTA, 0.25% SDS, 20 mmol/L NaBu histone deacetylase inhibitor, 1× complete protease inhibitors cocktail EDTA free (Roche, 5056489001)]. Chromatin was sheared on the Bioruptor Pico (Diagenode) using a 15-second on/30-second off, 7-cycle regimen. Sheared chromatin was magnetically immunoprecipitated using a polyclonal anti-H3K27ac antibody (Supplementary Table S4) and subsequently tagmentated using the Auto ChIPmentation Kit for Histones (Diagenode, C01011010) on the IP-Star Compact Automated System (Diagenode, B03000002), according to the manufacturer's instructions. Input DNA was de–cross-linked, purified (MinElute, QIAGEN), and tagmentated (Nextera DNA Library Prep, Illumina). Stripping, end repair, and library amplification were performed according to the ChIPmentation Kit guidelines. Libraries were sequenced using the NextSeq500 (SR75, High Output). Reads were trimmed by Trimmomatic and aligned to hg38 with Bowtie2 using the parameters -N 1 -k 1. Peaks were called with MACS2 with the respective input control for each patient sample.
We analyzed H3K27ac mark around the PE sequence in normal human thymocytes using previously published ChIPmentation results (GSE151078).
ATAC-seq Analysis
Analyses of open chromatin regions around the PE sequence in a mouse NOTCH1-induced T-ALL cell line (previously described in refs. 1 and 6) were done by Active Motif, following well-established protocols (67). We analyzed open chromatin regions around the PE sequence in normal human thymocytes and human T-ALL primary samples using publicly available ATAC-seq datasets from GEO (GSE151075 and GSE124223).
Conservation Analysis
The conservation track in Supplementary Fig. S2 shows a subset of the vertebrates' genomic sequence alignment as produced by the UCSC Genome Browser based on the MULTIZ algorithm (68).
4C-seq Analyses
We performed 4C-seq analysis in JURKAT, HPB-ALL, and DND41 human T-ALL cells and HDΔP-NOTCH1–induced or ΔE-NOTCH1–induced mouse primary T-ALL lymphoblasts (previously described in ref. 11) following published protocols (69) and using the restriction enzymes HindIII (R3104M, NEB) and DpnII (R0543M, NEB). We constructed sequencing libraries from 4C DNA including barcoded Illumina adapters to the 5′ end of each PCR primer (Supplementary Table S3). We sequenced pooled libraries using a HiSeq 2500 sequencer (Illumina). Processing of 4C-seq data was done as previously described (12, 69). For visualization purposes, signal was averaged across three replicates using 25 kb bins with 90% overlap for smoothening.
Hi-C, Hi-ChIP, and Capture-HiC Analyses
We analyzed Hi-C of human normal CD4 cells or T-ALLs, as well as Hi-ChIP data in CUTTL1 T-ALL cells using available datasets from GEO (GSE115896). We analyzed tagHi-C data of mouse hematopoietic populations using an available dataset from GEO (GSE152918). Hi-C data were represented using Juicebox (70). Capture Hi-C data were represented using Capture Hi-C plotter (71).
Luciferase Reporter Assays
We performed reporter assays using a pGL3-Promoter Vector (E1761, Promega) luciferase construct alone and coupled with human and mouse PE enhancer sequences (hPE: hg19; chr10:90,144,435–90,149, 109; mPE: mm10; chr19:33,220,867–33,224,929; both sequences were synthesized by Genewiz, RRID:SCR_003177), cloned in the forward and reverse orientations. In these assays, we electroporated JURKAT cells with a Neon Transfection System Device (MPK5000, Thermo Fisher Scientific) using 100 μL tips (MPK10025, Thermo Fisher Scientific). Electroporation conditions were as follows: pulse voltage 1,350 V, pulse width 10 ms, pulse number 3, and cell density 107/mL. PE constructs were transfected together with a plasmid driving the expression of the Renilla luciferase gene (pCMV-Renilla) used as an internal control. We measured luciferase activity 42 hours after electroporation with the Dual-Luciferase Reporter Assay kit (E1980, Promega).
Mice and Animal Procedures
All animals were maintained in specific pathogen–free facilities at New Brunswick Rutgers Biomedical and Health Sciences Campus. All animal procedures were approved by the Rutgers Institutional Animal Care and Use Committee. PE-knockout and PE conditional knockout mice were generated by the Rutgers Cancer Institute of New Jersey Genome Editing Shared Resource.
To generate PE germinal knockout mice, CRISPR sgRNAs flanking the PE region (Supplementary Table S3) and HiFiCas9 (1081060, IDT) were used. sgRNAs were complexed with Cas9 protein and electroporated into C57BL/6J embryos (JAX:000664; RRID:IMSR_JAX:000664). Three correctly targeted founders were identified and bred to C57BL/6NTac (Taconic) to confirm germline transmission.
To generate PE conditional knockout mice, we introduced LoxP sites at mouse genomic positions GRCm38/mm10 chr19: 33,220,867–33,224,930, flanking the mouse PE enhancer. CRISPR sgRNAs flanking the PE region (Supplementary Table S3) and HiFiCas9 (1081060, IDT) were used to insert LoxP sites at the specified locations. sgRNAs were complexed with Cas9 protein, mixed with a donor oligo–containing LoxP sites, and electroporated into C57BL/6J embryos. LoxP sites were inserted sequentially with the 5′ LoxP site being done first. A male mouse with a 5′ LoxP site correctly inserted was used to perform in vitro fertilization with WT oocytes, and the 3′ LoxP was inserted in these embryos. Both LoxP sites were confirmed in founders and bred to C57BL/6NTac mice to confirm the presence of both LoxP sites on the same allele.
To generate conditional inducible PE-knockout mice, we bred PE conditional knockout mice with animals harboring a tamoxifen-inducible form of the Cre recombinase from the ubiquitous Rosa26 locus (72). Phenotypic characterization of PE-knockout mice included male and female mice and isogenic age- and gender-matched controls without blinding. For secondary deletion of PE using tamoxifen-induced activation of Cre recombinase in vivo, animals were randomly assigned to the vehicle (corn oil, C8267, Sigma-Aldrich) or tamoxifen (T5648, Sigma-Aldrich; 3 mg per mouse in corn oil) groups, and disease progression was followed without blinding.
IHC
Freshly isolated thymi from 6-week-old PE WT, heterozygous, or homozygous knockout littermates were fixed overnight in 3.7% zinc formalin solution (175, Anatech). We performed IHC stainings following standard protocols using hematoxylin and eosin stain or a monoclonal antibody against PTEN (Supplementary Table S4).
Generation of NOTCH1-Induced Leukemias
We generated NOTCH1-induced T-ALL tumors in mice as previously described (6, 11). Briefly, we performed retroviral transduction of bone marrow cells enriched in lineage-negative cells isolated using magnetic beads (130-090-858, Miltenyi Biotec) with an activated form of the NOTCH1 oncogene (δE-NOTCH1). Primary mouse lineage–negative bone marrow cells were cultured in Opti-MEM media (51985091, Life Technologies) supplemented with 10% FBS, penicillin G (100 U/mL), streptomycin (100 μg/mL), β-mercaptoethanol (55 μmol/L), IL3 (10 ng/mL; 213-13, PeproTech), IL6 (10 ng/mL; 216-16, PeproTech), IL7 (25 ng/mL; 217-17, PeproTech), SCF (50 ng/mL; 250-03, PeproTech), and Flt3L (50 ng/mL; 250-31L, PeproTech). NOTCH1-infected bone marrow progenitor cells were then transplanted by intravenous injection into lethally irradiated female isogenic recipients (C57BL/6NTac). No blinding was performed.
Flow Cytometry Analysis of Hematopoietic Populations
All flow cytometry data were collected on an Attune NxT Flow Cytometer (Thermo Fisher Scientific) and were later analyzed using FlowJo Software (BD; RRID:SCR_008520).
Single-cell suspensions of total thymocytes or splenocytes from 6-week-old PE WT, heterozygous, and homozygous knockout littermates were prepared by pressing tissues through a 70-μm filter. Red cells in spleen samples were removed by incubation with red blood cell lysis buffer (155 mmol/L NH4Cl, 12 mmol/L KHCO3, and 0.1 mmol/L EDTA) for 5 minutes at room temperature.
For flow cytometry–based analysis of thymocyte and splenocyte populations and discrete stages of T-cell development, cells were stained with anti-mouse fluorochrome-conjugated antibodies (Supplementary Table S4). T-cell subsets in the thymus were as follows: DN1 (CD4−CD8−CD44+CD25−), DN2 (CD4−CD8−CD44+CD25+), DN3 (CD4−CD8−CD44−CD25+) and DN4 (CD4−CD8−CD44−CD25−), DP (CD4+CD8+), SP CD4+ (CD4+CD8−), and SP CD8+ (CD4−CD8+).
For PTEN intracellular staining, DND41 control cells with heterozygous loss of PE were fixed and permeabilized with BD Cytofix/Cytoperm solution (554722, BD), washed with BD Perm/Wash buffer (554723, BD), and stained with a mouse anti-PTEN or a mouse isotype control antibody (Supplementary Table S4) following datasheet guidelines.
RNA-seq Gene Expression Profiling
PE conditional knockout NOTCH1–induced T-ALL–bearing mice were treated with vehicle only (corn oil, C8267, Sigma-Aldrich) or tamoxifen (T5648, Sigma-Aldrich; 3 mg per mouse in corn oil) to induce isogenic loss of PE via i.p. injection. Single-cell suspensions of total leukemic splenocytes were prepared by pressing leukemic spleens through a 70-μm filter. We removed red cells in spleen samples by incubation with red blood cell lysis buffer (155 mmol/L NH4Cl, 12 mmol/L KHCO3, and 0.1 mmol/L EDTA) for 5 minutes at room temperature. RNA was extracted using QIAshredder (#79656) and Rneasy Mini (#74106) QIAGEN kits. RNA library preparations and next-generation sequencing were performed at Genewiz (RRID:SCR_003177) using an Illumina Next-Seq platform (Illumina). We calculated gene counts for each sample using kallisto (73), normalized the counts, and evaluated differential expression between the samples in R using DESeq (74) and DESeq2 (62), respectively. We used GSEA to assess enrichment of PTEN signature genes as previously described in ref. 11, based on 10,000 permutations of the gene list (36).
scRNA-seq
Single-cell suspensions of total thymocytes from 6-week-old WT and PE homozygous knockout littermates were prepared by pressing thymic lobes through a 70-μm filter. This single-cell suspension was enriched in CD4−CD3− immature thymocyte progenitor cells by staining with biotinylated antibodies against CD3 and CD4 (Supplementary Table S4), followed by magnetic depletion of antibody-labeled cells using streptavidin microbeads (Miltenyi Biotec 130-048-101). CD4−CD3−-enriched thymocytes were then mixed 1:1 with single-cell suspensions from total unenriched thymocytes. We then loaded thymocyte suspensions (10,000 cells per sample) isolated from 6-week-old mice on a 10x Chromium instrument (10x Genomics) and prepared scRNA-seq libraries using the Chromium Single Cell 3′ v3 Reagent Kit (10x Genomics PN-1000092) and sequenced them at the Princeton University Genomics Core Facility using NovaSeq instruments (Illumina). We demultiplexed the reads and aligned them to the mm10 reference using 10x Cell Ranger version 3.0.2. We loaded single-cell data into count matrices and merged the PE WT and PE-knockout thymocytes using Seurat v3.0 (75), and calculated each cell's hybrid doublet score using scds (76). We retained cells that expressed 200 to 4,000 genes, contained less than 12% of transcripts mapping to mitochondrial genes, and had a hybrid doublet score ≤ 1. We then scaled and log-normalized gene expression measurements for each cell and performed principal component analysis using the top 2,000 genes selected by variance stabilizing transformation approach as implemented in Seurat. We used Harmony to integrate PE WT and knockout thymocyte populations and to remove dataset-specific characteristics (77). We used UMAP to visualize the distribution of cells in the projection of the significant principal components. We then extracted markers for the different thymocyte populations from data produced by the ImmGen project (GSE127267; ref. 78), and annotated cells accordingly, using singleR (79). The percentiles for each population were compared with population percentages determined by flow cytometry analysis in each genotype as ground truth. Cells with multiple or non–T-cell population calls were excluded from the final analysis. We statistically evaluated differential expression between the PE WT and PE-knockout populations within each cell type using the Wilcoxon rank-sum test and adjusted for the multiple hypotheses using the Bonferroni correction.
Quantitative Real-time PCR
We performed reverse transcription reactions with the High-Capacity cDNA Reverse Transcription Kit (4368814, Applied Biosystems) and analyzed the resulting cDNA products by quantitative real-time PCR using PowerUp SYBR Green Master Mix (A25778, Applied Biosystems) using a QuantStudio 3 Real-Time PCR System (Applied Biosystems). Relative expression levels were normalized using Gapdh or Actin as a reference control (Supplementary Table S3).
Western Blotting
Western blot analysis was performed using standard procedures. Antibodies used are listed in Supplementary Table S4.
Statistical Analyses
Unless otherwise noted, we evaluated statistical significance using two-tailed Student t test, assuming normality and equal distribution of variance between the different groups analyzed. Survival in mouse experiments was represented with Kaplan–Meier curves, and significance was estimated with the log-rank test. We performed all statistical analyses using Prism GraphPad (RRID:SCR_002798).
Data Availability
All generated data were deposited in the NCBI Gene Expression Omnibus (GEO; RRID:SCR_005012) under the following accession numbers: 4C-seq data, GSE160427; mouse PE conditional knockout T-ALL RNA-seq data, GSE150894; mouse thymus scRNA-seq data, GSE156070; human H3K27ac T-ALL ChIPmentation data, GSE155332; and mouse T-ALL ATAC-seq data, GSE160620.
Authors' Disclosures
M.A. Lawlor reports personal fees from C4 Therapeutics (consulting) outside the submitted work. J. Soulier reports grants from ERC and ANR during the conduct of the study. T. Taghon reports grants from Fund for Scientific Research Flanders and Cancer Research Institute Ghent during the conduct of the study. C.G. Mullighan reports personal fees from Illumina (scientific advisory board) during the conduct of the study, as well as grants from Loxo Oncology, Pfizer, and AbbVie and personal fees and other from Amgen (stock) outside the submitted work. H. Khiabanian reports grants from NIH and V Foundation during the conduct of the study. D. Herranz reports grants from the NIH, the American Cancer Society, Alex's Lemonade Stand Foundation, the Leukemia Research Foundation, the Children's Leukemia Research Association, and Gabrielle's Angel Foundation for Cancer Research during the conduct of the study. No disclosures were reported by the other authors.
Authors' Contributions
L. Tottone: Data curation, formal analysis, investigation, writing–review and editing. O. Lancho: Data curation, formal analysis, investigation, writing–review and editing. J.-W. Loh: Data curation, formal analysis, investigation. A. Singh: Data curation, formal analysis, investigation. S. Kimura: Data curation, formal analysis, visualization, writing–review and editing. J. Roels: Data curation, formal analysis, investigation, writing–review and editing. A. Kuchmiy: Data curation, formal analysis, investigation. S. Strubbe: Data curation, formal analysis, supervision, funding acquisition, investigation, writing–review and editing. M.A. Lawlor: Data curation, formal analysis, investigation, writing–review and editing. V. da Silva-Diz: Data curation, formal analysis, visualization, writing–review and editing. S. Luo: Data curation, formal analysis, investigation, writing–review and editing. S. Gachet: Data curation, investigation. C.A. García-Prieto: Data curation, formal analysis, investigation, writing–review and editing. R. Hagelaar: Data curation, formal analysis. M. Esteller: Data curation, formal analysis, supervision, writing–review and editing. J.P.P. Meijerink: Formal analysis, supervision. J. Soulier: Formal analysis, supervision. T. Taghon: Data curation, formal analysis, supervision, funding acquisition. P. Van Vlierberghe: Data curation, formal analysis, supervision, funding acquisition, writing–review and editing. C.G. Mullighan: Data curation, formal analysis, supervision, funding acquisition. H. Khiabanian: Data curation, formal analysis, supervision, funding acquisition, writing–review and editing. P.P. Rocha: Data curation, formal analysis, supervision, writing–review and editing. D. Herranz: Conceptualization, data curation, formal analysis, supervision, funding acquisition, writing–original draft, project administration, writing–review and editing.
Acknowledgments
The authors are grateful to Adolfo A. Ferrando (Columbia University Medical Center) and Antonio Maraver (IRCM, Montpellier) for their constructive criticism and support. They thank Peter Romanienko and Ghassan Yehia (Rutgers University) for outstanding technical assistance in the generation of PE germinal and conditional knockout mice. The authors thank Jessica Wiggins, Wei Wang and the Princeton University Genomics Core Facility for the single-cell library preparation and sequencing. They thank Iannis Aifantis and Aristotelis Tsirigos (New York University School of Medicine) for sharing already processed HiC files with them. The authors also thank Lorenzo Brunetti (University of Perugia) for his advice with the CRISPR/Cas9 electroporation protocol. They also thank everyone involved with JuanLord for their support. Figures 5A and 6A were created using BioRender.com.
Work in the laboratory of D. Herranz is supported by the NIH (R00CA197869 and R01CA236936), the American Cancer Society (RSG-19-161-01-TBE), Alex's Lemonade Stand Foundation, the Leukemia Research Foundation, the Children's Leukemia Research Association, Gabrielle's Angel Foundation for Cancer Research, and the Rutgers Cancer Institute of New Jersey. Work in the laboratory of H. Khiabanian is supported by the NIH (R01CA233662) and the V Foundation (T2019-012), as well as by Rutgers Office of Advanced Research Computing (NIH 1S10OD012346-01A1). In addition, Rutgers Cancer Institute of New Jersey Shared Resources, which are supported in part by the National Cancer Institute Cancer Center Support Grant P30CA072720, were instrumental for this project, including Biomedical Informatics Shared Resource (P30CA072720-5917), Flow Cytometry and Cell Sorting Shared Resource (P30CA072720-5921), and Genome Editing Shared Resource (P30CA072720-5922). Work in the lab of P.P. Rocha is supported by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Development (NICHD). The authors also thank the NICHD Molecular Genomics Core and computational resources of the NIH HPC Biowulf cluster (https://hpc.nih.gov). Work in the labs of P. Van Vlierberghe and T. Taghon is supported by the European Research Council (grant no. StG-639784), the Fund for Scientific Research Flanders (FWO), and the Cancer Research Institute Ghent (CRIG). Fellowships from the New Jersey Commission on Cancer Research support the work of L. Tottone (DCHS20PPC010), J.-W. Loh (DFS18PPC017), and V. da Silva-Diz (DCHS19PPC008).