Clonal Evolution Dissection Reveals High MSI2 Level Promotes Chemo-Resistance in T-Cell Acute Lymphoblastic Leukemia

Occurrences of T-cell acute lymphoblastic leukemia (T-ALL) arise in approximately 10-15% of pediatric ALL cases and the outcomes for high-risk T-ALL are inferior. Currently, the precise dissection of clonal evolution with intra-tumoral heterogeneity and to delineate targetable molecular events driving relapse in T-ALL subtypes have become area of interest. We carried out high-throughput droplet-based 5′-single-cell RNA-sequencing (scRNA-seq) and paired T-cell receptor sequencing (scTCR-seq) of isolated CD7 + BMMCs from five patients' Dx_Rel paired samples. To investigate the gene signature variation underlying leukemic progression, we observed two distinct patterns of evolutionary trajectories in the four Dx_Rel nonETP-ALL pairs. From T956 and T723, a significant outgrowth in TCR initially negligible diagnostic subclusters at relapse, to which surrogate TCR repertoires were correspondingly linked and enumerated, suggestive robust dominant “clonal shifts” under continuous chemotherapeutic pressure. Whereas in Dx_Rel pairs T593 and T856, which contained subclusters with distinct gene signatures at diagnosis or relapse respectively, shared an identical TCR genotype (TCR Dx-Rel major clone). The clonal architecture clues gained from simultaneous scRNA-seq and scTCR-seq are strongly indicative of the leukemic “clonal drift” that commonly occurs under chemotherapeutic stress, during which somehow, hypothetically, molecular signatures within the same TCR subclones intrinsically progressed over time ( Fig. A-H). Next, we designed a targeted scDNA-seq library within a custom panel covering 127 amplicons and 102 genes via Tapastri platform to validate these two patterns of clonal evolution during leukemia progression unambiguously. Likewise, we captured that in persisted subclones in T593 and T856 leukemic cells attained reproductive fitness from maintaining a mutational static status with remarkable transcriptomic variation from diagnosis to relapse. To further interpret the transcriptomic signature drift in the persisted TCR-defined major clones from T593 and T856 pairs, differential expression analysis was employed with RNA-binding protein Musashi-2 (MSI2) showing the most prominent differential expression pattern across drifted subclones catching our attention as a potential indicator of intrinsic clonal progression related to chemo-resistance. Comprehensive analysis of RNA-seq data from the CCCG-ALL-2015 cohort (n=72) and TARGET T-ALL cohort (n = 264) revealed