Single cell long-read genotyping of transcriptomes reveals discrete mechanisms of clonal evolution in post-myeloproliferative neoplasm acute myeloid leukemia.

Myeloproliferative neoplasms (MPNs) are caused by acquired mutations in hematopoietic stem and progenitor cells (HSPCs). The acquisition of additional mutations like TP53 and the overall mutational burden influence a patients risk of disease progression toward lethal post-MPN acute myeloid leukemia (AML). Recent technological advancements in linking single-cell gene expression with genotype have improved our understanding of tumor heterogeneity. However, current methodologies have limitations in simultaneously genotyping low-expression genes (such as JAK2) alongside other pathogenic loci. To address this, we developed a novel long read genotyping pipeline of cDNA transcripts called LOTR-Seq, which can genotype the full length of expressed transcripts of 30 genes at once. Using LOTR-Seq, we genotyped HSPCs at the JAK2V617 locus in 9,075 single cells from eight patients with chronic phase MPN (CP-MPN) and in 5,016 cells from four patients with post-MPN AML. We then linked the mutations to the single cell transcriptome of 29,712 JAK2V617F-driven CP-MPN cells and 16,895 post-MPN AML cells. In our analysis of post-MPN AMLs, we identified nine mutated loci across six genes (JAK2, IDH1/2, TP53, SRSF2, U2AF1) and linked these mutations to specific transcriptional phenotypes. Overall, LOTR-Seq provides novel insights into the evolution of post-MPN AML.