Outlier Expression of Isoforms by Targeted or Total RNA Sequencing Identifies Clinically Significant Genomic Variants in Hematolymphoid Tumors

Abstract

Recognition of aberrant gene isoforms due to DNA events can impact risk stratification and molecular classification of hematolymphoid tumors. In myelodysplastic syndromes, KMT2A partial tandem duplication (PTD) was one of the top adverse predictors in the International Prognostic Scoring System-Molecular study. In B-cell acute lymphoblastic leukemia (B-ALL), ERG isoforms have been proposed as markers of favorable-risk DUX4 rearrangements, whereas deletion-mediated IKZF1 isoforms are associated with adverse prognosis and have been extended to the high-risk IKZF1^plus signature defined by codeletions, including PAX5. In this limited study, outlier expression of isoforms as markers of IKZF1 intragenic or 3' deletions, DUX4 rearrangements, or PAX5 intragenic deletions were 92.3% (48/52), 90% (9/10), or 100% (9/9) sensitive, respectively, and 98.7% (368/373), 100% (35/35), or 97.1% (102/105) specific, respectively, by targeted RNA sequencing, and 84.0% (21/25), 85.7% (6/7), or 81.8% (9/11) sensitive, respectively, and 98.2% (109/111), 98.4% (127/129), or 98.7% (78/79) specific, respectively, by total RNA sequencing. Comprehensive split-read analysis identified expressed DNA breakpoints, cryptic splice sites associated with IKZF1 3' deletions, PTD of IKZF1 exon 5 spanning N159Y in B-ALL with mutated IKZF1 N159Y, and truncated KMT2A-PTD isoforms. Outlier isoforms were also effective targeted RNA markers for PAX5 intragenic amplifications (B-ALL), KMT2A-PTD (myeloid malignant cancers), and rare NOTCH1 intragenic deletions (T-cell acute lymphoblastic leukemia). These findings support the use of outlier isoform analysis as a robust strategy for detecting clinically significant DNA events.