Prevalence, clinical correlates and clonal inferences of concurrent splicing factor mutations in patients with myeloid neoplasms

Abstract Introduction: Somatic mutations in splicing factors (SF) - SF3B1, SRSF2, U2AF1, and ZRSR2 – represent one of the most common drivers of premalignant clonal hematopoesis. While SF mutations are potent drivers of clonal expansion (Fabre et al., Nature 2022), multiple SF mutations are rarely observed in the same individual (negative epistasis). Prior studies have identified rare cases of concurrent SF mutations (<4%) and described specific hotspot combinations that escape negative epistasis (Taylor et al., Blood 2020). We conducted this study to further characterize the frequency and clonal relationship between concurrent SF mutations using a larger patient cohort. Methods: After IRB approval, we assessed 3,369 myeloid next generation sequencing (NGS) results that identified somatic variants in SF3B1, SRSF2, U2AF1, or ZRSR2, in 2,500 adult patients from 08/2015 to 07/2025. Patient demographics, ICD-9/10 codes, and complete blood counts with differentials were obtained, retaining only those datapoints assigned 30 days before or after NGS testing. To allow comparison of clone size between autosomal and the X-linked SF ZRSR2, VAFs were multiplied by two to estimate the cancer cell fraction (CCF), except for ZRSR2 mutations in males. For patients harboring concurrent SF mutations, we categorized them as co-dominant if the difference in respective CCFs was < 10% of the greater CCF. Specific SF variants were categorized into common hotspots (n=13), except for ZRSR2 variants, which do not occur in defined hotspots. Fisher's exact testing was used to assess the enrichment of SF hotspot variants in patients with concurrent SF mutations, as well as combinations of SF variants in patients with co-dominant SF mutations. Results: We identified 2,331 patients with at least 1 pathogenic SF mutation, of which 6.7% (n=156) had > 1 SF mutation (multi-SF MT). The total cohort comprised cases of AML (n=549, 24%), MDS (n=977, 42%), MDS/MPN (n=233, 10%), MPN (n=228,10%), and uncategorized ICD codes (“Other”, n=326), of which the largest category was clonal cytopenia(s). Notably, the prevalence of multi-SF cases was higher than previously described and our cohort included 5 patients with 3 SF mutations and 1 with 4 SF mutations. Of the multi-SF MT patients, we identified 42 (27%) with co-dominant CCFs among their respective SF mutations. Thirty-three multi-SF MT patients had multiple NGS sampling timepoints (median 2, range 2-4) spanning a median of 0.97 years (range 0.06-5.69 years); we found that 17 of these 33 demonstrated correlative SF CCFs over multiple timepoints, suggesting their SF mutations occurred in the same clone (“co-clonal”). Of these 17 patients, 82% had co-dominant CCFs at the single-timepoint level, whereas none of the 16 non-co-clonal demonstrated co-dominance. Thus, we used SF co-dominance as a surrogate for cases likely to be co-clonal. ZRSR2 (p=8.8x10-29), U2AF1Q157 (p=9.4x10-9), SF3B1K666 (p=1.4x10-4) and infrequent SF3B1 (“SF3B1-Other”, p=.0029) mutations were overrepresented in the multi-SF MT cohort compared to patients with a single SF mutation. Within the multi-SF MT cohort, 3 SF mutation combinations showed overrepresentation in codominant cases – SRSF2P95R/H/L+ ZRSR2 (p=0.014), SRSF2P95R/H/L+ SF3B1K666 (p=0.044), and U2AF1S34 + U2AF1Q157 (p=0.024), although not significant after multiple hypothesis correction. Dual ZRSR2 mutations were the most common combination seen in the multi-SF MT cohort but were underrepresented in codominant cases (p=0.025). Although there was no disease subtype uniquely associated with total or co-dominant multi-SF MT cases, occurrence of multiple SF mutations was more common in males (8% vs. 3% in females), an effect driven by ZRSR2 mutations. Hemoglobin and WBC were lower in multi-SF MT cases, with the lower hemoglobin strongly driven by the codominant multi-SF MT cases (median 8.3 g/dL vs 9.4 g/dL, p=0.00046). Conclusions: These findings suggest that SF mutation co-occurrence in myeloid neoplasms is more prevalent than previously estimated (~7%). Furthermore, we identified 6 individuals harboring > 2 SF mutations, which has not been previously described. We validate specific SF variant pairs that appear to exist in the same clone and are working to confirm the clonal architecture and characterize global splicing aberration in selected cases with single-cell studies to fully uncover the mechanism by which negative epistasis is overcome.