Functional analysis of structural variants in single cells using Strand-seq.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Pediatric T-ALL type-1 and type-2 relapses develop along distinct pathways of clonal evolution.

The mechanisms underlying T-ALL relapse remain essentially unknown. Multilevel-omics in 38 matched pairs of initial and relapsed T-ALL revealed 18 (47%) type-1 (defined by being derived from the major ancestral clone) and 20 (53%) type-2 relapses (derived from a minor ancestral clone). In both types of relapse, we observed known and novel drivers of multidrug resistance including MDR1 and MVP, NT5C2 and JAK-STAT activators. Patients with type-1 relapses were specifically characterized by IL7R upregulation. In remarkable contrast, type-2 relapses demonstrated (1) enrichment of constitutional cancer predisposition gene mutations, (2) divergent genetic and epigenetic remodeling, and (3) enrichment of somatic hypermutator phenotypes, related to BLM, BUB1B/PMS2 and TP53 mutations. T-ALLs that later progressed to type-2 relapses exhibited a complex subclonal architecture, unexpectedly, already at the time of initial diagnosis. Deconvolution analysis of ATAC-Seq profiles showed that T-ALLs later developing into type-1 relapses resembled a predominant immature thymic T-cell population, whereas T-ALLs developing into type-2 relapses resembled a mixture of normal T-cell precursors. In sum, our analyses revealed fundamentally different mechanisms driving either type-1 or type-2 T-ALL relapse and indicate that differential capacities of disease evolution are already inherent to the molecular setup of the initial leukemia.

Aging of preleukemic thymocytes drives CpG island hypermethylation in T-cell acute lymphoblastic leukemia.

Cancer cells display DNA hypermethylation at specific CpG islands in comparison to their normal healthy counterparts, but the mechanism that drives this so-called CpG island methylator phenotype (CIMP) remains poorly understood. Here, we show that CpG island methylation in human T-cell acute lymphoblastic leukemia (T-ALL) mainly occurs at promoters of Polycomb Repressor Complex 2 (PRC2) target genes that are not expressed in normal or malignant T-cells and which display a reciprocal association with H3K27me3 binding. In addition, we revealed that this aberrant methylation profile reflects the epigenetic history of T-ALL and is established already in pre-leukemic, self-renewing thymocytes that precede T-ALL development. Finally, we unexpectedly uncover that this age-related CpG island hypermethylation signature in T-ALL is completely resistant to the FDA-approved hypomethylating agent Decitabine. Altogether, we here provide conceptual evidence for the involvement of a pre-leukemic phase characterized by self-renewing thymocytes in the pathogenesis of human T-ALL.

Chromatin accessibility landscape of pediatric T-lymphoblastic leukemia and human T-cell precursors.

We aimed at identifying the developmental stage at which leukemic cells of pediatric T-ALLs are arrested and at defining leukemogenic mechanisms based on ATAC-Seq. Chromatin accessibility maps of seven developmental stages of human healthy T cells revealed progressive chromatin condensation during T-cell maturation. Developmental stages were distinguished by 2,823 signature chromatin regions with 95% accuracy. Open chromatin surrounding SAE1 was identified to best distinguish thymic developmental stages suggesting a potential role of SUMOylation in T-cell development. Deconvolution using signature regions revealed that T-ALLs, including those with mature immunophenotypes, resemble the most immature populations, which was confirmed by TF-binding motif profiles. We integrated ATAC-Seq and RNA-Seq and found DAB1, a gene not related to leukemia previously, to be overexpressed, abnormally spliced and hyper-accessible in T-ALLs. DAB1-negative patients formed a distinct subgroup with particularly immature chromatin profiles and hyper-accessible binding sites for SPI1 (PU.1), a TF crucial for normal T-cell maturation. In conclusion, our analyses of chromatin accessibility and TF-binding motifs showed that pediatric T-ALL cells are most similar to immature thymic precursors, indicating an early developmental arrest.

PDX models recapitulate the genetic and epigenetic landscape of pediatric T-cell leukemia.

We compared 24 primary pediatric T-cell acute lymphoblastic leukemias (T-ALL) collected at the time of initial diagnosis and relapse from 12 patients and 24 matched patient-derived xenografts (PDXs). DNA methylation profile was preserved in PDX mice in 97.5% of the promoters (ρ = 0.99). Similarly, the genome-wide chromatin accessibility (ATAC-Seq) was preserved remarkably well (ρ = 0.96). Interestingly, both the ATAC regions, which showed a significant decrease in accessibility in PDXs and the regions hypermethylated in PDXs, were associated with immune response, which might reflect the immune deficiency of the mice and potentially the incomplete interaction between murine cytokines and human receptors. The longitudinal approach of this study allowed an observation that samples collected from patients who developed a type 1 relapse (clonal mutations maintained at relapse) preserved their genomic composition; whereas in patients who developed a type 2 relapse (subset of clonal mutations lost at relapse), the preservation of the leukemia's composition was more variable. In sum, this study underlines the remarkable genomic stability, and for the first time documents the preservation of the epigenomic landscape in T-ALL-derived PDX models.