Lasho, Terra

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.

Abstract Heterozygous germline variants of DDX41 occur in 2-4% of adult Myelodysplastic Syndrome (MDS) and Acute Myeloid Leukemia (AML) cases, making them the most prevalent cause of inherited predisposition to these diseases. Homozygous deleterious variants are never observed in the germline, and loss of both DDX41 alleles causes apoptosis of proliferative cells, consistent with essentiality. Most disease-associated germline variants cause loss of full-length protein expression (truncating). However, 30-40% are missense (non-truncating) and their role in disease risk compared to truncating variants remains in question. Previous reports have shown that both truncating and non-truncating Pathogenic/Likely Pathogenic (P/LP) variants (defined by ACMG scoring) are associated with DDX41-specific disease features, including co-mutation of the contralateral allele of DDX41. In contrast, variants of unknown significance (VUS) are less associated with these features. Acquired co-mutations of DDX41, of which p.R525H accounts for 50-70%, are restricted to DDX41-mutated patients and are always missense. Our previous studies indicated that p.R525H lacks the essential DDX41 function and cannot support hematopoietic progenitor proliferation. However, the frequency with which this mutation occurs and the absence of acquired truncating mutations suggest that p.R525H has a specific effect on DDX41 function that is distinct from complete loss and promotes pathogenesis specifically in the context of existing P/LP DDX41 mutation. To understand how germline non-truncating variants predispose to disease, we characterized the effect of six recurrent non-truncating P/LP variants (p.G173R, p.R219H, p.P258L, p.Y259C, p.K331del, and p.I396T) and three VUS (p.K9E, p.Y33H, and p.M155H) on DDX41 function. We utilized a genetic complementation approach in murine lineage-negative (Linneg) bone marrow cells with inducible deletion of DDX41. We found that all non-truncating P/LP and VUS forms of DDX41 could fully rescue the colony forming potential of Ddx41-deficient cells to a similar level as wild-type protein, whereas the p.R525H variant completely failed to rescue. To determine the effect of the non-truncating variants on DDX41 function at endogenous expression levels, we utilized CRISPR to modify the DDX41 alleles in MOLM13 human AML cells. We produced viable, proliferating cell lines with homozygous mutations for all 3 VUS and 4/6 P/LP variants. For p.P258L and p.K331del, like p.R525H, we could not produce homozygous mutant cell lines despite screening hundreds of clones and finding several with heterozygosity. To further study the P/LP variants that were functional in the cell line, we generated C57Bl/6 mouse lines harboring p.G173R, p.Y259C, and p.I396T mutations in the endogenous Ddx41 alleles. Heterozygous mutant mice were phenotypically normal in all three strains. However, we could not obtain homozygous mutant mice for p.G173R or p.Y259C, indicating these cannot support embryonic/fetal development. In contrast, p.I396T-homozygous mutant mice were born at Mendelian ratios and were phenotypically normal until at least 15 months of age. Since p.R525H mutation is common in patients with non-truncating P/LP variants, we sought to model the combined effect of these mutations. We utilized the Linneg cells with inducible knockout of Ddx41, co-expressed each non-truncating variant with p.R525H by viral transduction, and monitored the colony formation potential. Remarkably, cultures expressing P/LP variants or M155I combined with p.R525H had dramatically reduced colony formation, while wild-type DDX41 and the other 2 VUS supported normal colony formation when co-expressed with p.R525H. Thus, we conclude p.R525H has a dominant-negative effect on certain non-truncating DDX41 variants. This is consistent with the emerging idea that the presence of an acquired hotspot DDX41 co-mutation is sufficient to classify a non-truncating DDX41 variant as pathogenic.Our results demonstrate that many DDX41 missense variants retain protein function when overexpressed. However, at endogenous expression levels, some P/LP variants are unable to support cell proliferation, indicative of hypomorphic effects on protein function. In patients, acquisition of the p.R525H mutation in the contralateral allele disrupts the functionality of P/LP variants through dominant-negative effects, leading to hematopoietic progenitor defects that are hallmark of the DDX41 pathogenesis.

Abstract Introduction: TERT promoter (pTERT) variants are among the most frequent genetic alterations in human cancers, where they drive oncogenesis through telomerase reactivation. However, somatic, pathogenic pTERT variants have also been seen in non-cancerous conditions such as telomere biology disorders (TBD), where they attempt to rescue the bone marrow failure phenotype (adaptive compensation). We conducted this study to assess the phenotype and outcomes of pTERT mutations in TBD patients. Methods: After IRB approval, the Mayo Clinic institutional TBD database was queried to identify all patients with somatic pTERT mutations. We have designed a customized error corrected sequencing panel that covers somatic mutations (52 genes) associated with adaptive and maladaptive outcomes in TBD patients. This covers the promoter region for TERT and can detect variants down to 0.5% variant allele fraction (VAF) with 10,000 X sequencing depth. Clinical data was retrospectively abstracted. Results: We identified 7 pTERT variants (3 hotspots) in 6 (two females (33%), mean age 61 years, 4 alive). The three hotspots were: (1) c.-57A>C (2 individuals, mean VAF 17.1%; range 11.6%-22.6%), (2) c.-124C>T (4 individuals, mean VAF 12.9%; range 1.9%-41.3%), and (3) c.-146C>T (1 individual, VAF 1.2%). One individual had both the c.-124C>T (VAF 41.3%) and the c.-146C>T somatic pTERT variants, occurring in the context of a germline TERT variant. The median VAF for all 7 pTERT variants was 6.1% (range 1.2%-41.3%). The hotspot, c.-57A>C, occurred in the context of pathogenic germline TERT (VAF 49.9%) variants with pathogenic somatic co-variants in PPM1D (avg. VAF 1.5%). The second hotspot, c.-124C>T, occurred in the context of germline TERT (avg. VAF 49.6%) and RTEL1 (VAF 47.3%) variants with pathogenic somatic co-variants including ASXL1 (VAF 1%), DNMT3A (VAF 8%), and PPM1D (VAF 4%). Five of six assessable patients with pTERT variants had shortened telomere lengths (TL), with results pending for one individual. Among those with the c.-57A>C variant, granulocyte TLs were at or below the 10th percentile, while lymphocyte TLs were below the 10th percentile in one patient and below the 1st percentile in another. For those with the c.-124C>T variant, all granulocyte TLs were below the 1st percentile. Lymphocyte TLs in this group included one at the 10th percentile and two at or below the 1st percentile. The median follow up for all patients was 2 years. Hotspot 1 (c.-57A>C) variants were associated with normocellular bone marrow biopsies with one patient having normal blood counts and the other having chronic macrocytic anemia (grade 2) and thrombocytopenia (< grade 1). In patients with hotspot 2 and 3 somatic variants, bone marrow cellularity ranged from 15–50% (1 patient had grade 2-3 hypocellularity), with normal bone marrow morphology, except for polytypic plasma cells (10–15%) identified in a patient with IgG kappa multiple myeloma (currently MRD negative on maintenance therapy). This patient also had a history of head and neck cancer and follicular lymphoma. Of these four patients, only one had normal blood counts (co-occurring c.-146C>T variant); the rest had cytopenias (n=3, the worst being grade 2 thrombocytopenia), progressive macrocytosis (n=3), or eosinophilia (n=1). Comorbidities encountered in these patients included myasthenia gravis (1), hypogonadism (1), and sarcoidosis (1). Moreover, four patients (2 each for hotspots 1 and 2) presented with idiopathic pulmonary fibrosis, while liver fibrosis was seen in 4 patients (also 2 each for hotspots 1 and 2). At last follow up, none of these patients had documented progression to myelodysplastic syndrome or acute myeloid leukemia. Conclusions: We describe the spectrum and clinical features of somatic pTERT mutations in TBD patients. pTERT mutations involved the following hotspot regions, c.-124C>T, c.-146C>T, and c.-57A>C, and co-occurred with other somatic variants (i.e., PPM1D and DNMT3A). None of the 6 patients with pTERT mutations had a severe bone marrow failure phenotype or developed myeloid neoplasms at last follow up, exemplifying a unique mechanism of adaptive somatic rescue.

Abstract Background: DNMT3A and TET2 are epigenetic regulator genes commonly mutated in age-related clonal hematopoiesis (CH). Despite having opposed epigenetic functions, these mutations are associated with increased all-cause mortality and a low risk for hematological neoplasms. In this study, we make use of a natural inflammatory response occurring during COVID-19 to define the impact of these mutations on inflammation and outcomes. Methods: A cohort of 243 community-based patients (Olmsted County, MN) with COVID-19 were identified. CH status was determined on PBMC using error corrected sequencing (lower detection limit 0.5% VAF). In a subset of patients, we use multi-omics interrogation of the genome, transcriptome, and epigenome at single cell (sc) resolution along with DNA methylation, in PBMC, from patients with COVID-19 and CH. Results: We detected 97 CH mutations in 72 (29.6%) of 243 patients with COVID-19. The most frequent being DNMT3A (n=30, 30%) and TET2 (n=26, 28%). The presence of CH negatively impacted overall survival when restricted to TET2 and DNMT3A mutations. However, only DNMT3Amt CH retained age and comorbidity independent prognostication (HR 2.84, p = 0.022). We then assessed for changes in DNA methylation using Illumina MethylationEPIC array on patients with COVID-19 and DNMT3A (n=3) or TET2mt (n=4) CH. Site specific DNA methylation changes revealed significant loss of methylation in DNMT3Amt CH in comparison to TET2mt patients, with 10,944 hypomethylated and 1,160 hypermethylated sites. We found that actively transcribed states (Tx, TxWk) were more commonly hypomethylated in DNMT3Amt CH while enhancers (Enh) and promoters (TssA, TssAFlnk) were more commonly hypomethylated in TET2mt CH. We then conducted scDNA-seq in DNMT3Amt (n=3) and TET2mt (n=1) CH in COVID-19 and observed that while in TET2mt CH, mutations were largely restricted to classical and intermediate monocytes, in DNMT3Amt CH, mutations were seen in myeloid and lymphoid lineage cells. We performed scRNA-seq on patients with COVID-19 and DNMT3A (n=3) or TET2mt (n=6) CH. Differential gene expression analysis identified 1,569 upregulated genes in DNMT3Amt patients and 205 downregulated. Pathway analysis demonstrated an upregulation of genes involved in T cell survival, function, and inflammation in DNMT3Amt patients compared to TET2mt patients with a significant overexpression of GIMAP1, GIMAP4, IRF2, IL32 in multiple cell types. We correlated transcriptomic changes with chromatin accessibility using the 10X Genomics Multiome platform to profile both gene expression and open chromatin from the same samples. Analysis of global distribution of cut sites and differentially accessible peaks showed increased chromatin accessibility in DNMT3Amt CH, especially CD4+ T lymphocytes and NK cells. Conclusion: We demonstrate an age, sex, and comorbidity-independent adverse impact on mortality associated with DNMT3Amt CH in the context of COVID-19 and highlight transcriptomic and epigenetic differences between DNMT3A and TET2mt CH, with single cell resolution. Citation Format: Jenna A Fernandez, Wazim Mohammed Ismail, Moritz Binder, Terra Lasho, Susan M. Geyer, Amelia Mazzone, Christy M. Finke, Abhishek Mangaonkar, Jeong-Heon Lee, Liguo Wang, Vernadette A. Simon, Fariborz Rakhshan Rohakthar, Amik Munankarmy, Jonathan J. Harrington, Melissa R. Snyder, Keith D. Robertson, Alexandre Gaspar-Maia, Mrinal Patnaik. Divergent impacts of DNMT3A and TET2 mutant clonal hematopoiesis on COVID-19 related outcomes [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Methylation, Clonal Hematopoiesis, and Cancer; 2025 Feb 1-4; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2025;85(3 Suppl):Abstract nr A012.

Abstract Background: Age-related clonal hematopoiesis (ARCH) is a precursor state to blood cancers and is associated with epigenetic age acceleration (EpiAgeAccel). While ARCH is associated with increased all-cause mortality and cardiovascular disease, not much is known about these risks with therapy-related clonal hematopoiesis (t-CH) – CH occurring in the context of prior cytotoxic chemotherapy/radiation exposure. Methods: Patients were prospectively recruited for CH and EpiAge assessment from Mayo Clinic between 2022 and 2024. CH status was determined using error corrected sequencing on peripheral blood (PB) DNA (sensitivity 0.5% VAF). Patients with CH were classified as t-CH if they had received chemotherapy or radiation prior to CH detection. All other CH cases were classifed as ARCH. For EpiAge analysis, PB DNA was analyzed using the Illumina MethylationEPIC array. The resulting beta values were used to obtain epigentic age estimates using Horvath’s clock, Hannum’s Clock, PhenoAge, DunedinPACE, and estimated telomere lengths (DNAmTL). Results: 106 (62%) of 170 patients assessed were found to have CH. CH patients were classified into ARCH-63 (59%) and t-CH-43 (41%). While the median chronological age (69 vs 67 years) was similar between ARCH and t-CH, t-CH patients had lower Hgb [10.9g/dL vs. 13g/dL (p=0.0001)], lower platelet [146x109/L vs. 192x109/L (p=0.0132)] and lower WBC [4.4x109/L vs. 5.2x109/L (p=0.0345)]. The 3 most common ARCH mutations were TET2 (20%), DNMT3A (19%), and ASXL1 (14%), while PPM1D (27%), TET2 (18%), and DNMT3A (14%) were most frequent in t-CH. There was no significant difference in median EpiAge via Horvath (68 vs 67.5years p=0.2826), Hannum (65 vs 66 years p=0.9333), or PhenoAge (55 vs 56 years p=0.618) between ARCH and t-CH. However, DunedinPACE of aging was higher in t-CH (1.27 vs 1.18 bioyear/chronoyear p=0.0148) and median DNAmTL was shorter in t-CH (6.35 vs 6.49kb p=0.0247), compared to ARCH (Note patients with telomere biology disorders were excluded). Multivariable regression analysis adjusting for chronologic age, smoking status, and male sex showed that both ARCH and t-CH had EpiAgeAccel via Horvath(5.8years p=0.0001, 3.4 years p=0.051), Hannum(5years p=0.001, 4.5 years p=0.007), and PhenoAge (6.7years p=0.002, 7.6 years p=0.002) clocks. t-CH more than doubled DunedinPACE (0.2 bioyear/chronoyear, p<0.001) of aging compared to ARCH (0.08 bioyear/chronoyear, p=0.026). t-CH also demonstrated a significant decrease in DNAmTL (0.34kb, p<0.001) in comparison to ARCH (0.17kb, p=0.015). With a median follow-up of 20 months, t-CH was associated with worse overall survival in comparison to ARCH, and non-CH patients (p=0.013, 93% vs. 90% alive at 10 months, 87 vs. 72% alive at 20 months). Conclusion: Our prospective CH cohort study validates the association between CH and EpiAgeAccel using several different methylation clocks and for the first time highlights an increased pace of aging and shorter DNAmTL in the context of inferior survival outcomes in t-CH. Citation Format: Marissa Li, Jenna Fernandez, Kristina Kirschner, Tamir Chandra, Kristen McCullough, Sarah Aug, Terra Lasho, Christy Finke, James Foran, Cecilia Arana Yi, Yael Kusne, Talha Badar, Abhishek Mangaonkar, Naseema Gangat, Rachel Simon, Alejandro Ferrer, Mrinal Patnaik. A prospective cohort study of biological aging and epigenetic age acceleration in 106 patients with age and therapy-related clonal hematopoiesis [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Methylation, Clonal Hematopoiesis, and Cancer; 2025 Feb 1-4; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2025;85(3 Suppl):Abstract nr A015.

Abstract Background DNMT3A and TET2 are epigenetic regulator genes commonly mutated in age-related clonal hematopoiesis (CH). Despite having opposing epigenetic functions, these mutations are associated with increased all-cause mortality and a low risk for progression to hematological neoplasms. While individual impacts on the epigenome have been described using different model systems, the phenotypic complexity in humans remains to be elucidated. Results Here we make use of a natural inflammatory response occurring during coronavirus disease 2019 (COVID-19), to understand the association of these mutations with inflammatory morbidity and mortality. We demonstrate the age-independent, negative impact of DNMT3A mutant CH on COVID-19-related cytokine release severity and mortality. Using single cell proteogenomics we show that DNMT3A mutations involve myeloid and lymphoid cells. Using single cell multiomics sequencing, we identify cell-specific gene expression changes associated with DNMT3A mutations, along with significant epigenomic deregulation affecting enhancer accessibility, resulting in overexpression of IL32, a proinflammatory cytokine that can result in inflammasome activation in monocytes and macrophages. Finally, we show with single cell resolution that the loss of function of DNMT3A is directly associated with increased chromatin accessibility in mutant cells. Conclusions We demonstrate the negative prognostic impact of DNMT3Amt CH on COVID-19 related inflammatory morbidity and mortality. DNMT3Amt CH involves myeloid and lymphoid cells and in the context of COVID-19, was associated with inflammatory transcriptional priming, resulting in overexpression of IL32. This overexpression was secondary to increased chromatic accessibility, specific to DNMT3Amt CH cells. DNMT3Amt CH can serve as a potential biomarker for adverse inflammatory outcomes.

Abstract Clonal cytopenia of undetermined significance (CCUS) is defined by a myeloid driver mutation in the context of otherwise unexplained cytopenia. CCUS has an inherent risk of progressing to myeloid neoplasm. However, it is unknown how exposure to previous cytotoxic therapy may impact the risk of progression and survival. We stratified patients with CCUS by prior exposure to DNA-damaging therapy. Of 151 patients, 46 (30%) had received cytotoxic therapy and were classified as therapy-related CCUS (t-CCUS), whereas 105 (70%) had de novo CCUS. A lower proportion of t-CCUS had hypercellular marrows (17.8% vs 44.8%, P = .002) but had higher median bone marrow blast percentages. After a median follow-up of 2.2 years, t-CCUS had significantly shorter progression-free survival (PFS, 1.8 vs 6.3 years; hazard ratio [HR], 2.1; P = .007) and median overall survival (OS; 3.6 years vs not reached; HR, 2.3; P = .007) compared with CCUS. Univariable and multivariable time-to-event analyses showed that exposure to cytotoxic therapy independently accounted for inferior PFS and OS. Despite the similarities in clinical presentation between CCUS and t-CCUS, we show that exposure to prior cytotoxic therapies was an independent risk factor for inferior outcomes. This suggests that t-CCUS represents a unique clinical entity that needs more stringent monitoring or earlier intervention strategies.

Abstract VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome, caused by somatic mutations in UBA1, is an autoinflammatory disorder with diverse systemic manifestations. Thrombosis is a prominent clinical feature of VEXAS syndrome. The risk factors and frequency of thrombosis in VEXAS syndrome are not well described, due to the disease’s recent discovery and the paucity of large databases. We evaluated 119 patients with VEXAS syndrome for venous and arterial thrombosis and correlated their presence with clinical outcomes and survival. Thrombosis occurred in 49% of patients, mostly venous thromboembolism (VTE; 41%). Almost two-thirds of VTEs were unprovoked, 41% were recurrent, and 20% occurred despite anticoagulation. The cumulative incidence of VTE was 17% at 1 year from symptom onset and 40% by 5 years. Cardiac and pulmonary inflammatory manifestations were associated with time to VTE. M41L was positively associated specifically with pulmonary embolism by univariate (odds ratio [OR]: 4.58, confidence interval [CI] 1.28-16.21, P = .02) and multivariate (OR: 16.94, CI 1.99-144.3, P = .01) logistic regression. The cumulative incidence of arterial thrombosis was 6% at 1 year and 11% at 5 years. The overall survival of the entire patient cohort at median follow-up time of 4.8 years was 88%, and there was no difference in survival between patients with or without thrombosis (P = .8). Patients with VEXAS syndrome are at high risk of VTE; thromboprophylaxis should administered be in high-risk settings unless strongly contraindicated.