HSPA9/mortalin inhibition disrupts erythroid maturation through a TP53-dependent mechanism in human CD34+ hematopoietic progenitor cells

Myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic stem cell malignancies characterized by abnormal hematopoietic cell maturation, increased apoptosis of bone marrow cells, and anemia. They are the most common myeloid blood cancers in American adults. The full complement of gene mutations that contribute to the phenotypes or clinical symptoms in MDS is not fully understood. Around 10%-25% of MDS patients harbor an interstitial heterozygous deletion on the long arm of chromosome 5 [del(5q)], creating haploinsufficiency for a large set of genes, including HSPA9. The HSPA9 gene encodes for the protein mortalin, a highly conserved heat shock protein predominantly localized in mitochondria. Our prior study showed that knockdown of HSPA9 induces TP53-dependent apoptosis in human CD34+ hematopoietic progenitor cells. In this study, we explored the role of HSPA9 in regulating erythroid maturation using human CD34+ cells. We inhibited the expression of HSPA9 using gene knockdown and pharmacological inhibition and found that inhibition of HSPA9 disrupted erythroid maturation as well as increased expression of p53 in CD34+ cells. To test whether the molecular mechanism of HSPA9 regulating erythroid maturation is TP53-dependent, we knocked down HSPA9 and TP53 individually or in combination in human CD34+ cells. We found that the knockdown of TP53 partially rescued the erythroid maturation defect induced by HSPA9 knockdown, suggesting that the defect in cells with reduced HSPA9 expression is TP53-dependent. Collectively, these findings indicate that reduced levels of HSPA9 may contribute to the anemia observed in del(5q)-associated MDS patients due to the activation of TP53.

Plasma Proteomic Signature Predicts Myeloid Neoplasm Risk

PURPOSE

Clonal hematopoiesis (CH) is thought to be the origin of myeloid neoplasms (MN). Yet our understanding of the mechanisms driving CH progression to MN and clinical risk prediction of MN remains limited. The human proteome reflects complex interactions between genetic and epigenetic regulation of biological systems. We hypothesized that the plasma proteome might predict MN risk and inform our understanding of the mechanisms promoting MN development.

EXPERIMENTAL DESIGN

We jointly characterized CH and plasma proteomic profiles of 46,237 individuals in the UK Biobank at baseline study entry. During 500,036 person-years of follow-up, 115 individuals developed MN. Cox proportional hazard regression was used to test for an association between plasma protein levels and MN risk.

RESULTS

We identified 115 proteins associated with MN risk of which 30% (N=34) were also associated with CH. These were enriched for known regulators of the innate and adaptive immune system. Plasma proteomics improved the prediction of MN risk (AUC=0.85, p=5×10-9) beyond clinical factors and CH (AUC=0.80). In an independent group (N=381,485), we used inherited polygenic risk scores (PRS) for plasma protein levels to validate the relevance of these proteins to MN development. PRS analyses suggest that most MN-associated proteins we identified are not directly causally linked to MN risk, but rather represent downstream markers of pathways regulating the progression of CH to MN.

CONCLUSIONS

These data highlight the role of immune cell regulation in the progression of CH to MN and the promise of leveraging multi-omic characterization of CH to improve MN risk stratification.

R-Loop Accumulation in Spliceosome Mutant Leukemias Confers Sensitivity to PARP1 Inhibition by Triggering Transcription-Replication Conflicts

RNA splicing factor (SF) gene mutations are commonly observed in patients with myeloid malignancies. Here we showed that SRSF2- and U2AF1-mutant leukemias are preferentially sensitive to PARP inhibitors (PARPi), despite being proficient in homologous recombination repair. Instead, SF-mutant leukemias exhibited R-loop accumulation that elicited an R-loop-associated PARP1 response, rendering cells dependent on PARP1 activity for survival. Consequently, PARPi induced DNA damage and cell death in SF-mutant leukemias in an R-loop-dependent manner. PARPi further increased aberrant R-loop levels, causing higher transcription-replication collisions and triggering ATR activation in SF-mutant leukemias. Ultimately, PARPi-induced DNA damage and cell death in SF-mutant leukemias could be enhanced by ATR inhibition. Finally, the level of PARP1 activity at R-loops correlated with PARPi sensitivity, suggesting that R-loop-associated PARP1 activity could be predictive of PARPi sensitivity in patients harboring SF gene mutations. This study highlights the potential of targeting different R-loop response pathways caused by spliceosome gene mutations as a therapeutic strategy for treating cancer.

SIGNIFICANCE

Spliceosome-mutant leukemias accumulate R-loops and require PARP1 to resolve transcription-replication conflicts and genomic instability, providing rationale to repurpose FDA-approved PARP inhibitors for patients carrying spliceosome gene mutations.

Phase I-II Trial of Early Azacitidine after Matched Unrelated Donor Hematopoietic Cell Transplantation

Graft-versus-host disease (GVHD) is a major complication after allogeneic hematopoietic cell transplantation (allo-HCT). The hypomethylating agent azacitidine (AZA) has been shown to be effective in preclinical and clinical studies for the prevention of acute GVHD (aGVHD). We sought to determine the maximum tolerated dose (MTD) of AZA when given on days 1 to 5 of a 28-day cycle for 4 cycles, starting on day +7 after allo-HCT, as well as its impact on aGVHD and chronic GVHD (cGVHD), relapse, and overall survival (OS) in patients undergoing matched unrelated donor allo-HCT. This study was a single-arm, single-center, open-label phase I-II study with a total of 15 and 38 patients enrolled in the phase I and II portions of the trial, respectively. A standard 3+3 study design was used in phase I, and all patients in phase II received AZA at the MTD determined in phase I. The MTD of AZA starting at day +7 post-transplantation was 45 mg/m2. Phase II of the study was halted after enrolling 38 of the planned 46 patients following an interim analysis that suggested futility. Overall, AZA at 45 mg/m2 exhibited a side effect profile consistent with prior reports and had a minimal impact on engraftment. The cumulative incidence of clinically significant aGVHD by day +180 was 39.9% (95% confidence interval [CI], 22% to 53.7%). The incidence of all-grade cGVHD was 61.4% (95% CI, 40.3% to 75%). At 1 year, OS was 73.7% (95% CI, 60.9% to 89.1%), and the disease relapse rate was 11.4% (95% CI, .2% to 21.3%). Our results suggest that early post-allo-HCT AZA has limited efficacy in preventing aGVHD and cGVHD but could have a beneficial effect in preventing disease relapse.

Genomic landscape of TP53-mutated myeloid malignancies

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

Utility of targeted gene sequencing to differentiate myeloid malignancies from other cytopenic conditions

The National Heart, Lung, and Blood Institute-funded National MDS Natural History Study (NCT02775383) is a prospective cohort study enrolling patients with cytopenia with suspected myelodysplastic syndromes (MDS) to evaluate factors associated with disease. Here, we sequenced 53 genes in bone marrow samples harvested from 1298 patients diagnosed with myeloid malignancy, including MDS and non-MDS myeloid malignancy or alternative marrow conditions with cytopenia based on concordance between independent histopathologic reviews (local, centralized, and tertiary to adjudicate disagreements when needed). We developed a novel 2-stage diagnostic classifier based on mutational profiles in 18 of 53 sequenced genes that were sufficient to best predict a diagnosis of myeloid malignancy and among those with a predicted myeloid malignancy, predict whether they had MDS. The classifier achieved a positive predictive value (PPV) of 0.84 and negative predictive value (NPV) of 0.8 with an area under the receiver operating characteristic curve (AUROC) of 0.85 when classifying patients as having myeloid vs no myeloid malignancy based on variant allele frequencies (VAFs) in 17 genes and a PPV of 0.71 and NPV of 0.64 with an AUROC of 0.73 when classifying patients as having MDS vs non-MDS malignancy based on VAFs in 10 genes. We next assessed how this approach could complement histopathology to improve diagnostic accuracy. For 99 of 139 (71%) patients (PPV of 0.83 and NPV of 0.65) with local and centralized histopathologic disagreement in myeloid vs no myeloid malignancy, the classifier-predicted diagnosis agreed with the tertiary pathology review (considered the internal gold standard).

A functional link between lariat debranching enzyme and the intron-binding complex is defective in non-photosensitive trichothiodystrophy

The pre-mRNA life cycle requires intron processing; yet, how intron-processing defects influence splicing and gene expression is unclear. Here, we find that TTDN1/MPLKIP, which is encoded by a gene implicated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to spliceosomal function. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, whereas its N-terminal intrinsically disordered region (IDR) binds the intron-binding complex (IBC). TTDN1 loss, or a mutated IDR, causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. A Ttdn1-deficient mouse model recapitulates intron-processing defects and certain neurodevelopmental phenotypes seen in NP-TTD. Fusing DBR1 to the TTDN1 IDR is sufficient to recruit DBR1 to the IBC and circumvents the functional requirement for TTDN1. Collectively, our findings link RNA lariat processing with splicing outcomes by revealing the molecular function of TTDN1.

Increased clonal hematopoiesis involving DNA damage response genes in patients undergoing lung transplantation

BACKGROUNDCellular stressors influence the development of clonal hematopoiesis (CH). We hypothesized that environmental, inflammatory, and genotoxic stresses drive the emergence of CH in lung transplant recipients. METHODSWe performed a cross-sectional cohort study of 85 lung transplant recipients to characterize CH prevalence. We evaluated somatic variants using duplex error-corrected sequencing and germline variants using whole exome sequencing. We evaluated CH frequency and burden using χ2 and Poisson regression, and we evaluated associations with clinical and demographic variables and clinical outcomes using χ2, logistic regression, and Cox regression. RESULTSCH in DNA damage response (DDR) genes TP53, PPM1D, and ATM was increased in transplant recipients compared with a control group of older adults (28% versus 0%, adjusted OR [aOR], 12.9 [1.7-100.3], P = 0.0002). Age (OR, 1.13 [1.03-1.25], P = 0.014) and smoking history (OR 4.25 [1.02-17.82], P = 0.048) were associated with DDR CH. Germline variants predisposing to idiopathic pulmonary fibrosis were identified but not associated with CH. DDR CH was associated with increased cytomegalovirus viremia versus patients with no (OR, 7.23 [1.95-26.8], P = 0.018) or non-DDR CH (OR, 7.64 [1.77-32.89], P = 0.024) and mycophenolate discontinuation (aOR, 3.8 [1.3-12.9], P = 0.031). CONCLUSIONCH in DDR genes is prevalent in lung transplant recipients and is associated with posttransplant outcomes including cytomegalovirus activation and mycophenolate intolerance. FUNDINGNIH/NHLBI K01HL155231 (LKT), R25HL105400 (LKT), Foundation for Barnes-Jewish Hospital (LKT), Evans MDS Center at Washington University (KAO, MJW), ASH Scholar Award (KAO), NIH K12CA167540 (KAO), NIH P01AI116501 (AEG, DK), NIH R01HL094601 (AEG), and NIH P01CA101937 (DCL).

Impact of U2AF1 mutations on circular RNA expression in myelodysplastic neoplasms

Mutations in U2AF1 are relatively common in myelodysplastic neoplasms (MDS) and are associated with an inferior prognosis, but the molecular mechanisms underlying this are not fully elucidated. Circular RNAs (circRNAs) have been implicated in cancer, but it is unknown how mutations in splicing factors may impact on circRNA biogenesis. Here, we used RNA-sequencing to investigate the effects of U2AF1 mutations on circRNA expression in K562 cells with a doxycycline-inducible U2AF1^S34 mutation, in a mouse model with a doxycycline-inducible U2AF1^S34 mutation, and in FACS-sorted CD34+ bone marrow cells from MDS patients with either U2AF1^S34 or U2AF1^Q157 mutations. In all contexts, we found an increase in global circRNA levels in the U2AF1-mutated setting, which was independent of expression changes in the cognate linear host genes. In patients, the U2AF1^S34 and U2AF1^Q157 mutations were both associated with an overall increased expression of circRNAs. circRNAs generated by a non-Alu-mediated mechanism generally showed the largest increase in expression levels. Several well-described cancer-associated circRNAs, including circZNF609 and circCSNK1G3, were upregulated in MDS patients with U2AF1 mutations compared to U2AF1-wildtype MDS controls. In conclusion, high circRNA expression is observed in association with U2AF1 mutations in three biological systems, presenting an interesting possibility for biomarker and therapeutic investigation.

Clonal Hematopoiesis and Risk of Incident Lung Cancer

PURPOSE

To prospectively examine the association between clonal hematopoiesis (CH) and subsequent risk of lung cancer.

METHODS

Among 200,629 UK Biobank (UKBB) participants with whole-exome sequencing, CH was identified in a nested case-control study of 832 incident lung cancer cases and 3,951 controls (2006-2019) matched on age and year at blood draw, sex, race, and smoking status. A similar nested case-control study (141 cases/652 controls) was conducted among 27,975 participants with whole-exome sequencing in the Mass General Brigham Biobank (MGBB, 2010-2021). In parallel, we compared CH frequency in published data from 5,003 patients with solid tumor (2,279 lung cancer) who had pretreatment blood sequencing performed through Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets.

RESULTS

In UKBB, the presence of CH was associated with increased risk of lung cancer (cases: 12.5% v controls: 8.7%; multivariable-adjusted odds ratio [OR], 1.36; 95% CI, 1.06 to 1.74). The association remained robust after excluding participants with chronic obstructive pulmonary disease. No significant interactions with known risk factors, including polygenic risk score and C-reactive protein, were identified. In MGBB, we observed similar enrichment of CH in lung cancer (cases: 15.6% v controls: 12.7%). The meta-analyzed OR (95% CI) of UKBB and MGBB was 1.35 (1.08 to 1.68) for CH overall and 1.61 (1.19 to 2.18) for variant allele frequencies ≥ 10%. In Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets, CH with a variant allele frequency ≥ 10% was enriched in pretreatment lung cancer compared with other tumors after adjusting for age, sex, and smoking (OR for lung v breast cancer: 1.61; 95% CI, 1.03 to 2.53).

CONCLUSION

Independent of known risk factors, CH is associated with increased risk of lung cancer.

Genomic landscape of TP53 -mutated myeloid malignancies

TP53 -mutated myeloid malignancies are most frequently associated with complex cytogenetics. The presence of complex and extensive structural variants complicates detailed genomic analysis by conventional clinical techniques. We performed whole genome sequencing of 42 AML/MDS cases with paired normal tissue to characterize the genomic landscape of TP53 -mutated myeloid malignancies. The vast majority of cases had multi-hit involvement at the TP53 genetic locus (94%), as well as aneuploidy and chromothripsis. Chromosomal patterns of aneuploidy differed significantly from TP53 -mutated cancers arising in other tissues. Recurrent structural variants affected regions that include ETV6 on chr12p, RUNX1 on chr21, and NF1 on chr17q. Most notably for ETV6 , transcript expression was low in cases of TP53 -mutated myeloid malignancies both with and without structural rearrangements involving chromosome 12p. Telomeric content is increased in TP53 -mutated AML/MDS compared other AML subtypes, and telomeric content was detected adjacent to interstitial regions of chromosomes. The genomic landscape of TP53 -mutated myeloid malignancies reveals recurrent structural variants affecting key hematopoietic transcription factors and telomeric repeats that are generally not detected by panel sequencing or conventional cytogenetic analyses.

Key Points

WGS comprehensively determines TP53 mutation status, resulting in the reclassification of 12% of cases from mono-allelic to multi-hit Chromothripsis is more frequent than previously appreciated, with a preference for specific chromosomes ETV6 is deleted in 45% of cases, with evidence for epigenetic suppression in non-deleted cases NF1 is mutated in 48% of cases, with multi-hit mutations in 17% of these cases TP53 -mutated AML/MDS is associated with altered telomere content compared with other AMLs.

Genomic profiling for clinical decision making in myeloid neoplasms and acute leukemia

Myeloid neoplasms and acute leukemias derive from the clonal expansion of hematopoietic cells driven by somatic gene mutations. Although assessment of morphology plays a crucial role in the diagnostic evaluation of patients with these malignancies, genomic characterization has become increasingly important for accurate diagnosis, risk assessment, and therapeutic decision making. Conventional cytogenetics, a comprehensive and unbiased method for assessing chromosomal abnormalities, has been the mainstay of genomic testing over the past several decades and remains relevant today. However, more recent advances in sequencing technology have increased our ability to detect somatic mutations through the use of targeted gene panels, whole-exome sequencing, whole-genome sequencing, and whole-transcriptome sequencing or RNA sequencing. In patients with myeloid neoplasms, whole-genome sequencing represents a potential replacement for both conventional cytogenetic and sequencing approaches, providing rapid and accurate comprehensive genomic profiling. DNA sequencing methods are used not only for detecting somatically acquired gene mutations but also for identifying germline gene mutations associated with inherited predisposition to hematologic neoplasms. The 2022 International Consensus Classification of myeloid neoplasms and acute leukemias makes extensive use of genomic data. The aim of this report is to help physicians and laboratorians implement genomic testing for diagnosis, risk stratification, and clinical decision making and illustrates the potential of genomic profiling for enabling personalized medicine in patients with hematologic neoplasms.

IL-1β expression in bone marrow dendritic cells is induced by TLR2 agonists and regulates HSC function

Hematopoietic stem/progenitor cells (HSPCs) reside in localized microenvironments, or niches, in the bone marrow that provide key signals regulating their activity. A fundamental property of hematopoiesis is the ability to respond to environmental cues such as inflammation. How these cues are transmitted to HSPCs within hematopoietic niches is not well established. Here, we show that perivascular bone marrow dendritic cells (DCs) express a high basal level of Toll-like receptor-1 (TLR1) and TLR2. Systemic treatment with a TLR1/2 agonist induces HSPC expansion and mobilization. It also induces marked alterations in the bone marrow microenvironment, including a decrease in osteoblast activity and sinusoidal endothelial cell numbers. TLR1/2 agonist treatment of mice in which Myd88 is deleted specifically in DCs using Zbtb46-Cre show that the TLR1/2-induced expansion of multipotent HPSCs, but not HSPC mobilization or alterations in the bone marrow microenvironment, is dependent on TLR1/2 signaling in DCs. Interleukin-1β (IL-1β) is constitutively expressed in both murine and human DCs and is further induced after TLR1/2 stimulation. Systemic TLR1/2 agonist treatment of Il1r1-/- mice show that TLR1/2-induced HSPC expansion is dependent on IL-1β signaling. Single-cell RNA-sequencing of low-risk myelodysplastic syndrome bone marrow revealed that IL1B and TLR1 expression is increased in DCs. Collectively, these data suggest a model in which TLR1/2 stimulation of DCs induces secretion of IL-1β and other inflammatory cytokines into the perivascular niche, which in turn, regulates multipotent HSPCs. Increased DC TLR1/2 signaling may contribute to altered HSPC function in myelodysplastic syndrome by increasing local IL-1β expression.

Convergent Clonal Evolution of Signaling Gene Mutations Is a Hallmark of Myelodysplastic Syndrome Progression

Progression from myelodysplastic syndromes (MDS) to secondary acute myeloid leukemia (AML) is associated with the acquisition and expansion of subclones. Our understanding of subclone evolution during progression, including the frequency and preferred order of gene mutation acquisition, remains incomplete. Sequencing of 43 paired MDS and secondary AML samples identified at least one signaling gene mutation in 44% of MDS and 60% of secondary AML samples, often below the level of standard sequencing detection. In addition, 19% of MDS and 47% of secondary AML patients harbored more than one signaling gene mutation, almost always in separate, coexisting subclones. Signaling gene mutations demonstrated diverse patterns of clonal evolution during disease progression, including acquisition, expansion, persistence, and loss of mutations, with multiple patterns often coexisting in the same patient. Multivariate analysis revealed that MDS patients who had a signaling gene mutation had a higher risk of AML progression, potentially providing a biomarker for progression.

SIGNIFICANCE

Subclone expansion is a hallmark of progression from MDS to secondary AML. Subclonal signaling gene mutations are common at MDS (often at low levels), show complex and convergent patterns of clonal evolution, and are associated with future progression to secondary AML. See related article by Guess et al., p. 316 (33). See related commentary by Romine and van Galen, p. 270. This article is highlighted in the In This Issue feature, p. 265.

Abstract 841: Mutant Trp53 R172H has gain-of-function or dominant-negative effects in response to hematopoietic stressors in mice

Background: Missense mutations of the TP53 gene frequently occur in patients with myelodysplastic syndrome and acute myeloid leukemia. The functional consequences of mutant p53 expression on hematopoiesis remain incompletely studied. Whether TP53 mutations confer a loss-of-function, gain-of-function, or dominant-negative effect in response to chemotherapies has not been fully explored. One of these hotspot mutations is a G-to-A substitution at nucleotide 525 of TP53, which codes for an Arg175His (R175H) substitution. Thus, we characterized hematopoiesis in Trp53R172H constitutive knock-in mice, corresponding to R175H TP53 mutation in patients, and compared the results to Trp53 knock-out and control mice.

Methods: Five genotypes of mice were analyzed (Trp53+/+ [WT], Trp53R172H/+, Trp53R172H/R172H, Trp53+/-, and Trp53-/-). Frequencies of peripheral blood (PB) and bone marrow (BM) cells were measured by flow cytometry and fluorouracil (5-FU) and N-ethyl-N-nitrosourea (ENU) were used as hematopoietic stressors.

Results: There was no change in the PB counts or lineage distribution of mature PB and BM cells in Trp53R172H/+ and Trp53R172H/R172H mice compared to WT or Trp53 knock-out mice at 8-15 weeks of age. We found that expression of Trp53R172H mutation in mice does not cause hematopoietic stem cell (HSC) expansion under basal conditions, in contrast to heterozygous and homozygous deletion of Trp53, but does confer a functional advantage to HSCs in a competitive repopulation assay. Trp53R172H HSCs expand further following ENU exposure, similar to Trp53-/- cells and consistent with a dominant-negative effect of Trp53R172H. Mechanistically, mutant p53R172H has a dominant-negative effect on p21-induced expression following ENU treatment. Repeated 5-FU injections drive proliferation of progenitors. To test the hematopoietic cell-intrinsic effects of mutant p53R172H on white blood cell count recovery, we exposed mice reconstituted with Trp53+/+, Trp53+/-, Trp53-/- or Trp53R172H/+ BM cells to a single dose of 5-FU and followed mice weekly for a month. Trp53R172H/+ mice were least affected by the 5-FU challenge and showed the fastest and largest recovery of WBC count compared to Trp53+/+, Trp53+/- and Trp53-/- mice. These results demonstrate that mutant p53R172H confers a gain-of-function property following 5-FU exposure that was not seen with deletion of Trp53.

Conclusion: Our data suggest that in vivo expression of mutant Trp53R172H has minimal effects on basal hematopoiesis. However, mutant Trp53R172H induces either a gain-of-function or a dominant-negative effect depending on the type of hematopoietic stress and hematopoietic phenotype being observed. Overall, our results indicate that the functional consequences induced by various Trp53 mutations (e.g., deletion and point mutation) are not identical, but are influenced by the type of hematopoietic stress.

Citation Format: Tuoen Liu, Tanzir Ahmed, Michael O. Alberti, Brian A. Wadugu, Jin Shao, Sarah Grieb, Grace Przybyl, Matthew Ndonwi, Matthew J. Walter. Mutant Trp53 R172H has gain-of-function or dominant-negative effects in response to hematopoietic stressors in mice [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 841.

Toll-like receptor and cytokine expression throughout the bone marrow differs between patients with low- and high-risk myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders, the pathogenesis of which involves enhanced immune signaling that promotes or selects for mutant hematopoietic stem and progenitor cells (HSPCs). In particular, toll-like receptor (TLR) expression and signaling are enhanced in MDS, and their inhibition is an attractive therapeutic strategy. Although prior studies have reported increased expression of TLR2 and its binding partners TLR1 and TLR6 in the CD34+ cells of patients with MDS (especially those with low-risk disease), TLR expression in other cell types throughout the bone marrow is largely unknown. To address this, we used mass cytometry to assess the expression of TLR1, TLR2, and TLR6 and cytokines in the bone marrow hematopoietic cells of six low/intermediate-risk and six high-risk unmatched MDS bone marrow samples, as well as healthy controls, both at baseline and in response to TLR agonists. We observed several consistent differences between the groups. Most notably, TLR expression was upregulated in multiple cell populations in the low/intermediate-risk, but not high-risk, patients. In addition, many cytokines, including interleukin-6, interleukin-8, tumor necrosis factor α, transforming growth factor β, macrophage inflammatory protein 1β, and granzyme B, were highly expressed from various cell types in low/intermediate-risk patients. However, these same cytokines, with the exception of transforming growth factor β, were expressed at lower levels in high-risk MDS. Together, these findings highlight the differential role of inflammation, and specifically TLR expression, in low/intermediate- versus high-risk MDS, and suggest that elevated TLR expression and cytokine production in multiple cell types likely influences the pathogenesis of MDS in lower-risk patients.

Failure to Detect Mutations in U2AF1 due to Changes in the GRCh38 Reference Sequence

The U2AF1 gene is a core part of mRNA splicing machinery and frequently contains somatic mutations that contribute to oncogenesis in myelodysplastic syndrome, acute myeloid leukemia, and other cancers. A change introduced in the GRCh38 version of the human reference build prevents detection of mutations in this gene, and others, by variant calling pipelines. This study describes the problem in detail and shows that a modified GRCh38 reference build with unchanged coordinates can be used to ameliorate the issue.

Focal disruption of DNA methylation dynamics at enhancers in IDH-mutant AML cells

Recurrent mutations in IDH1 or IDH2 in acute myeloid leukemia (AML) are associated with increased DNA methylation, but the genome-wide patterns of this hypermethylation phenotype have not been comprehensively studied in AML samples. We analyzed whole-genome bisulfite sequencing data from 15 primary AML samples with IDH1 or IDH2 mutations, which identified ~4000 focal regions that were uniquely hypermethylated in IDH^mut samples vs. normal CD34+ cells and other AMLs. These regions had modest hypermethylation in AMLs with biallelic TET2 mutations, and levels of 5-hydroxymethylation that were diminished in IDH and TET-mutant samples, indicating that this hypermethylation results from inhibition of TET-mediated demethylation. Focal hypermethylation in IDH^mut AMLs occurred at regions with low methylation in CD34+ cells, implying that DNA methylation and demethylation are active at these loci. AML samples containing IDH and DNMT3A^R882 mutations were significantly less hypermethylated, suggesting that IDH^mut-associated hypermethylation is mediated by DNMT3A. IDH^mut-specific hypermethylation was highly enriched for enhancers that form direct interactions with genes involved in normal hematopoiesis and AML, including MYC and ETV6. These results suggest that focal hypermethylation in IDH-mutant AML occurs by altering the balance between DNA methylation and demethylation, and that disruption of these pathways at enhancers may contribute to AML pathogenesis.

Genetic and Transcriptional Contributions to Relapse in Normal Karyotype Acute Myeloid Leukemia

To better understand clonal and transcriptional adaptations after relapse in patients with acute myeloid leukemia (AML), we collected presentation and relapse samples from six normal karyotype AML cases. We performed enhanced whole-genome sequencing to characterize clonal evolution, and deep-coverage single-cell RNA sequencing on the same samples, which yielded 142,642 high-quality cells for analysis. Identifying expressed mutations in individual cells enabled us to discriminate between normal and AML cells, to identify coordinated changes in the genome and transcriptome, and to identify subclone-specific cell states. We quantified the coevolution of genetic and transcriptional heterogeneity during AML progression, and found that transcriptional changes were significantly correlated with genetic changes. However, transcriptional adaptation sometimes occurred independently, suggesting that clonal evolution does not represent all relevant biological changes. In three cases, we identified cells at diagnosis that likely seeded the relapse. Finally, these data revealed a conserved relapse-enriched leukemic cell state bearing markers of stemness, quiescence, and adhesion. SIGNIFICANCE: These data enabled us to identify a relapse-enriched leukemic cell state with distinct transcriptional properties. Detailed case-by-case analyses elucidated the complex ways in which the AML genome, transcriptome, and immune microenvironment interact to evade chemotherapy. These analyses provide a blueprint for evaluating these factors in larger cohorts.This article is highlighted in the In This Issue feature, p. 1.

U2af1 is a haplo-essential gene required for hematopoietic cancer cell survival in mice

Somatic mutations in the spliceosome gene U2AF1 are common in patients with myelodysplastic syndromes. U2AF1 mutations that code for the most common amino acid substitutions are always heterozygous, and the retained WT allele is expressed, suggesting that mutant hematopoietic cells may require the residual WT allele to be viable. We show that hematopoiesis and RNA splicing in U2af1 heterozygous knockout mice were similar to those in control mice, but that deletion of the WT allele in U2AF1(S34F) heterozygous mutant-expressing hematopoietic cells (i.e., hemizygous mutant) was lethal. These results confirm that U2AF1 mutant hematopoietic cells are dependent on the expression of WT U2AF1 for survival in vivo and that U2AF1 is a haplo-essential cancer gene. Mutant U2AF1(S34F)-expressing cells were also more sensitive to reduced expression of WT U2AF1 than nonmutant cells. Furthermore, mice transplanted with leukemia cells expressing mutant U2AF1 had significantly reduced tumor burden and improved survival after the WT U2af1 allele was deleted compared with when it was not deleted. These results suggest that selectively targeting the WT U2AF1 allele in heterozygous mutant cells could induce cancer cell death and be a therapeutic strategy for patients harboring U2AF1 mutations.

Mutant U2AF1-induced alternative splicing of H2afy (macroH2A1) regulates B-lymphopoiesis in mice

Somatic mutations in spliceosome genes are found in ∼50% of patients with myelodysplastic syndromes (MDS), a myeloid malignancy associated with low blood counts. Expression of the mutant splicing factor U2AF1(S34F) alters hematopoiesis and mRNA splicing in mice. Our understanding of the functionally relevant alternatively spliced target genes that cause hematopoietic phenotypes in vivo remains incomplete. Here, we demonstrate that reduced expression of H2afy1.1, an alternatively spliced isoform of the histone H2A variant gene H2afy, is responsible for reduced B cells in U2AF1(S34F) mice. Deletion of H2afy or expression of U2AF1(S34F) reduces expression of Ebf1 (early B cell factor 1), a key transcription factor for B cell development, and mechanistically, H2AFY is enriched at the EBF1 promoter. Induced expression of H2AFY1.1 in U2AF1(S34F) cells rescues reduced EBF1 expression and B cells numbers in vivo. Collectively, our data implicate alternative splicing of H2AFY as a contributor to lymphopenia induced by U2AF1(S34F) in mice and MDS.

Bam-readcount - rapid generation of basepair-resolution sequence metrics

Bam-readcount is a utility for generating low-level information about sequencing data at specific nucleotide positions. Originally designed to help filter genomic mutation calls, the metrics it outputs are useful as input for variant detection tools and for resolving ambiguity between variant callers1,2. In addition, it has found broad applicability in diverse fields including tumor evolution, single-cell genomics, climate change ecology, and tracking community spread of SARS-CoV-2.3-6.

Nonsense-Mediated RNA Decay Is a Unique Vulnerability of Cancer Cells Harboring SF3B1 or U2AF1 Mutations

Nonsense-mediated RNA decay (NMD) is recognized as an RNA surveillance pathway that targets aberrant mRNAs with premature translation termination codons (PTC) for degradation, however, its molecular mechanisms and roles in health and disease remain incompletely understood. In this study, we developed a novel reporter system to accurately measure NMD activity in individual cells. A genome-wide CRISPR-Cas9 knockout screen using this reporter system identified novel NMD-promoting factors, including multiple components of the SF3B complex and other U2 spliceosome factors. Interestingly, cells with mutations in the spliceosome genes SF3B1 and U2AF1, which are commonly found in myelodysplastic syndrome (MDS) and cancers, have overall attenuated NMD activity. Compared with wild-type (WT) cells, SF3B1- and U2AF1-mutant cells were more sensitive to NMD inhibition, a phenotype that is accompanied by elevated DNA replication obstruction, DNA damage, and chromosomal instability. Remarkably, the sensitivity of spliceosome mutant cells to NMD inhibition was rescued by overexpression of RNase H1, which removes R-loops in the genome. Together, these findings shed new light on the functional interplay between NMD and RNA splicing and suggest a novel synthetic lethal strategy for the treatment of MDS and cancers with spliceosome mutations. SIGNIFICANCE: This study has developed a novel NMD reporter system and identified a potential therapeutic approach of targeting the NMD pathway to treat cancer with spliceosome gene mutations.

Genome Sequencing as an Alternative to Cytogenetic Analysis in Myeloid Cancers

BACKGROUND

Genomic analysis is essential for risk stratification in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS). Whole-genome sequencing is a potential replacement for conventional cytogenetic and sequencing approaches, but its accuracy, feasibility, and clinical utility have not been demonstrated.

METHODS

We used a streamlined whole-genome sequencing approach to obtain genomic profiles for 263 patients with myeloid cancers, including 235 patients who had undergone successful cytogenetic analysis. We adapted sample preparation, sequencing, and analysis to detect mutations for risk stratification using existing European Leukemia Network (ELN) guidelines and to minimize turnaround time. We analyzed the performance of whole-genome sequencing by comparing our results with findings from cytogenetic analysis and targeted sequencing.

RESULTS

Whole-genome sequencing detected all 40 recurrent translocations and 91 copy-number alterations that had been identified by cytogenetic analysis. In addition, we identified new clinically reportable genomic events in 40 of 235 patients (17.0%). Prospective sequencing of samples obtained from 117 consecutive patients was performed in a median of 5 days and provided new genetic information in 29 patients (24.8%), which changed the risk category for 19 patients (16.2%). Standard AML risk groups, as defined by sequencing results instead of cytogenetic analysis, correlated with clinical outcomes. Whole-genome sequencing was also used to stratify patients who had inconclusive results by cytogenetic analysis into risk groups in which clinical outcomes were measurably different.

CONCLUSIONS

In our study, we found that whole-genome sequencing provided rapid and accurate genomic profiling in patients with AML or MDS. Such sequencing also provided a greater diagnostic yield than conventional cytogenetic analysis and more efficient risk stratification on the basis of standard risk categories. (Funded by the Siteman Cancer Research Fund and others.).

SF3B1-mutant MDS as a distinct disease subtype: a proposal from the International Working Group for the Prognosis of MDS

The 2016 revision of the World Health Organization classification of tumors of hematopoietic and lymphoid tissues is characterized by a closer integration of morphology and molecular genetics. Notwithstanding, the myelodysplastic syndrome (MDS) with isolated del(5q) remains so far the only MDS subtype defined by a genetic abnormality. Approximately half of MDS patients carry somatic mutations in spliceosome genes, with SF3B1 being the most commonly mutated one. SF3B1 mutation identifies a condition characterized by ring sideroblasts (RS), ineffective erythropoiesis, and indolent clinical course. A large body of evidence supports recognition of SF3B1-mutant MDS as a distinct nosologic entity. To further validate this notion, we interrogated the data set of the International Working Group for the Prognosis of MDS (IWG-PM). Based on the findings of our analyses, we propose the following diagnostic criteria for SF3B1-mutant MDS: (1) cytopenia as defined by standard hematologic values, (2) somatic SF3B1 mutation, (3) morphologic dysplasia (with or without RS), and (4) bone marrow blasts <5% and peripheral blood blasts <1%. Selected concomitant genetic lesions represent exclusion criteria for the proposed entity. In patients with clonal cytopenia of undetermined significance, SF3B1 mutation is almost invariably associated with subsequent development of overt MDS with RS, suggesting that this genetic lesion might provide presumptive evidence of MDS in the setting of persistent unexplained cytopenia. Diagnosis of SF3B1-mutant MDS has considerable clinical implications in terms of risk stratification and therapeutic decision making. In fact, this condition has a relatively good prognosis and may respond to luspatercept with abolishment of the transfusion requirement.

Systematic Analysis of Splice-Site-Creating Mutations in Cancer

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

TP53 mutation status divides myelodysplastic syndromes with complex karyotypes into distinct prognostic subgroups

Risk stratification is critical in the care of patients with myelodysplastic syndromes (MDS). Approximately 10% have a complex karyotype (CK), defined as more than two cytogenetic abnormalities, which is a highly adverse prognostic marker. However, CK-MDS can carry a wide range of chromosomal abnormalities and somatic mutations. To refine risk stratification of CK-MDS patients, we examined data from 359 CK-MDS patients shared by the International Working Group for MDS. Mutations were underrepresented with the exception of TP53 mutations, identified in 55% of patients. TP53 mutated patients had even fewer co-mutated genes but were enriched for the del(5q) chromosomal abnormality (p < 0.005), monosomal karyotype (p < 0.001), and high complexity, defined as more than 4 cytogenetic abnormalities (p < 0.001). Monosomal karyotype, high complexity, and TP53 mutation were individually associated with shorter overall survival, but monosomal status was not significant in a multivariable model. Multivariable survival modeling identified severe anemia (hemoglobin < 8.0 g/dL), NRAS mutation, SF3B1 mutation, TP53 mutation, elevated blast percentage (>10%), abnormal 3q, abnormal 9, and monosomy 7 as having the greatest survival risk. The poor risk associated with CK-MDS is driven by its association with prognostically adverse TP53 mutations and can be refined by considering clinical and karyotype features.

The National MDS Natural History Study: design of an integrated data and sample biorepository to promote research studies in myelodysplastic syndromes

Myelodysplastic syndromes (MDS), a spectrum of heterogeneous hematopoietic stem cell diseases, vary in clinical severity, response to therapy, and propensity toward progression to acute myeloid leukemia. These are acquired clonal disorders resulting from somatic mutations within the hematopoietic stem or progenitor cell population. Understanding the natural history and the risk of developing leukemia and other adverse outcomes is dependent on access to well-annotated biospecimens linked to robust clinical and molecular data. To facilitate the acquisition and distribution of MDS biospecimens to the wider scientific community and support scientific discovery in this disease, the National MDS Natural History study was initiated by the National Heart, Lung, and Blood Institute (NHLBI) and is being conducted in collaboration with community hospitals and academic medical centers supported by the National Cancer Institute (NCI). The study will recruit up to 2000 MDS patients or overlapping myeloproliferative neoplasms (MDS/MPN) and up to 500 cases of idiopathic cytopenia of undetermined significance (ICUS). The National MDS Natural History Study (NCT02775383) will offer the world's largest disease-focused tissue biobank linked to longitudinal clinical and molecular data in MDS. Here, we report on the study design features and describe the vanguard phase of 200 cases. The study assembles a comprehensive clinical database, quality of life results, laboratory data, histopathology slides and images, genetic information, hematopoietic and germline tissues representing high-quality biospecimens and data from diverse centers across the United States. These resources will be available to the scientific community for investigator-initiated research.

Immune Escape of Relapsed AML Cells after Allogeneic Transplantation

BACKGROUND

As consolidation therapy for acute myeloid leukemia (AML), allogeneic hematopoietic stem-cell transplantation provides a benefit in part by means of an immune-mediated graft-versus-leukemia effect. We hypothesized that the immune-mediated selective pressure imposed by allogeneic transplantation may cause distinct patterns of tumor evolution in relapsed disease.

METHODS

We performed enhanced exome sequencing on paired samples obtained at initial presentation with AML and at relapse from 15 patients who had a relapse after hematopoietic stem-cell transplantation (with transplants from an HLA-matched sibling, HLA-matched unrelated donor, or HLA-mismatched unrelated donor) and from 20 patients who had a relapse after chemotherapy. We performed RNA sequencing and flow cytometry on a subgroup of these samples and on additional samples for validation.

RESULTS

On exome sequencing, the spectrum of gained and lost mutations observed with relapse after transplantation was similar to the spectrum observed with relapse after chemotherapy. Specifically, relapse after transplantation was not associated with the acquisition of previously unknown AML-specific mutations or structural variations in immune-related genes. In contrast, RNA sequencing of samples obtained at relapse after transplantation revealed dysregulation of pathways involved in adaptive and innate immunity, including down-regulation of major histocompatibility complex (MHC) class II genes ( HLA-DPA1, HLA-DPB1, HLA-DQB1, and HLA-DRB1) to levels that were 3 to 12 times lower than the levels seen in paired samples obtained at presentation. Flow cytometry and immunohistochemical analysis confirmed decreased expression of MHC class II at relapse in 17 of 34 patients who had a relapse after transplantation. Evidence suggested that interferon-γ treatment could rapidly reverse this phenotype in AML blasts in vitro.

CONCLUSIONS

AML relapse after transplantation was not associated with the acquisition of relapse-specific mutations in immune-related genes. However, it was associated with dysregulation of pathways that may influence immune function, including down-regulation of MHC class II genes, which are involved in antigen presentation. These epigenetic changes may be reversible with appropriate therapy. (Funded by the National Cancer Institute and others.).

Mutation Clearance after Transplantation for Myelodysplastic Syndrome

BACKGROUND

Allogeneic hematopoietic stem-cell transplantation is the only curative treatment for patients with myelodysplastic syndrome (MDS). The molecular predictors of disease progression after transplantation are unclear.

METHODS

We sequenced bone marrow and skin samples from 90 adults with MDS who underwent allogeneic hematopoietic stem-cell transplantation after a myeloablative or reduced-intensity conditioning regimen. We detected mutations before transplantation using enhanced exome sequencing, and we evaluated mutation clearance by using error-corrected sequencing to genotype mutations in bone marrow samples obtained 30 days after transplantation. In this exploratory study, we evaluated the association of a mutation detected after transplantation with disease progression and survival.

RESULTS

Sequencing identified at least one validated somatic mutation before transplantation in 86 of 90 patients (96%); 32 of these patients (37%) had at least one mutation with a maximum variant allele frequency of at least 0.5% (equivalent to 1 heterozygous mutant cell in 100 cells) 30 days after transplantation. Patients with disease progression had mutations with a higher maximum variant allele frequency at 30 days than those who did not (median maximum variant allele frequency, 0.9% vs. 0%; P<0.001). The presence of at least one mutation with a variant allele frequency of at least 0.5% at day 30 was associated with a higher risk of progression (53.1% vs. 13.0%; conditioning regimen-adjusted hazard ratio, 3.86; 95% confidence interval [CI], 1.96 to 7.62; P<0.001) and a lower 1-year rate of progression-free survival than the absence of such a mutation (31.3% vs. 59.3%; conditioning regimen-adjusted hazard ratio for progression or death, 2.22; 95% CI, 1.32 to 3.73; P=0.005). The rate of progression-free survival was lower among patients who had received a reduced-intensity conditioning regimen and had at least one persistent mutation with a variant allele frequency of at least 0.5% at day 30 than among patients with other combinations of conditioning regimen and mutation status (P≤0.001). Multivariate analysis confirmed that patients who had a mutation with a variant allele frequency of at least 0.5% detected at day 30 had a higher risk of progression (hazard ratio, 4.48; 95% CI, 2.21 to 9.08; P<0.001) and a lower 1-year rate of progression-free survival than those who did not (hazard ratio for progression or death, 2.39; 95% CI, 1.40 to 4.09; P=0.002).

CONCLUSIONS

The risk of disease progression was higher among patients with MDS in whom persistent disease-associated mutations were detected in the bone marrow 30 days after transplantation than among those in whom these mutations were not detected. (Funded by the Leukemia and Lymphoma Society and others.).

Spliceosome Mutations Induce R Loop-Associated Sensitivity to ATR Inhibition in Myelodysplastic Syndromes

Heterozygous somatic mutations in spliceosome genes (U2AF1, SF3B1, ZRSR2, or SRSF2) occur in >50% of patients with myelodysplastic syndrome (MDS). These mutations occur early in disease development, suggesting that they contribute to MDS pathogenesis and may represent a unique genetic vulnerability for targeted therapy. Here, we show that RNA splicing perturbation by expression of the U2AF1(S34F) mutant causes accumulation of R loops, a transcription intermediate containing RNA:DNA hybrids and displaced single-stranded DNA, and elicits an ATR response. ATR inhibitors (ATRi) induced DNA damage and cell death in U2AF1(S34F)-expressing cells, and these effects of ATRi were enhanced by splicing modulating compounds. Moreover, ATRi-induced DNA damage was suppressed by overexpression of RNaseH1, an enzyme that specifically removes the RNA in RNA:DNA hybrids, suggesting that the ATRi sensitivity of U2AF1(S34F)-expressing cells arises from R loops. Taken together, our results demonstrate that ATR may represent a novel therapeutic target in patients with MDS carrying the U2AF1(S34F) mutation and potentially other malignancies harboring spliceosome mutations.Significance: This study provides preclinical evidence that patients with MDS or other myeloid malignancies driven by spliceosome mutations may benefit from ATR inhibition to exploit the R loop-associated vulnerability induced by perturbations in splicing. Cancer Res; 78(18); 5363-74. ©2018 AACR.

Abstract 2362: Comprehensive portrait of canonical and non-canonical splicing in cancer

Current annotation methods typically classify mutations as disruptors of splicing if they fall on either the consensus intronic dinucleotide splice donor, GT, or the splice acceptor, AG. As a group, splice site mutations have been presumed to be invariably deleterious because of their disruption of the conserved sequences that are used to identify exon-intron boundaries. While this classification method has been useful, increasing evidence suggests that mutations outside of the canonical splice site can lead to transcriptional changes beyond disruption of the nearest junction. In this study, we have developed the MiSplice pipeline to determine the local and global effects of genomic mutations on splicing across 33 cancer types. To evaluate the local effects of mutations on splicing, we applied MiSplice to identify splice-disrupting mutations (SDMs) and splice-creating mutations (SCMs), genome-wide. We identified 1,964 novel SCMs, of which 26% and 11% were mis-annotated as missense and silent mutations and validated 10 of 11 genes in a mini-gene splicing assay. SDM identification predicted complex splicing patterns associated with canonical splice site mutations and identified a handful of mutations in proximity to the canonical junction that disrupt splicing factor binding sites. Interestingly, further investigation of the novel neoantigens produced by SCMs and SDMs are likely several folds more immunogenic than missense mutations. To explore the global impact of mutations on splicing we focused on recurrent and adjacent mutations disrupting the spliceosomal complex and related splicing factor genes from over 400 curated splice related genes. In addition, we applied HotSpot3D to identify splice-factor mutations (SFMs) that are significantly proximal to one another. Our analysis has identified novel SFMs that disrupt the spliceosomal complex and globally impact downstream splicing targets creating novel peptide sequences and alter key cancer genes. The current study has greatly extended the insight into the transcriptional ramifications of genomic alterations by integrating DNA and RNA sequencing data and painting the portrait of alternative splicing across cancer genomes.

Citation Format: Reyka G. Jayasinghe, Song Cao, Qingsong Gao, Matthew A. Wyczalkowski, Sohini Sengupta, Matthew J. Walter, Christopher Maher, Michael C. Wendl, Feng Chen, Eduardo Eyras, Alexander J. Lazar, Ken Chen, Ilya Shmulevich, Li Ding, The Splicing Analysis Working Group. Comprehensive portrait of canonical and non-canonical splicing in cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2362.

Subclones dominate at MDS progression following allogeneic hematopoietic cell transplant

Allogeneic hematopoietic cell transplantation (alloHCT) is a potentially curative treatment for myelodysplastic syndromes (MDS), but patients who relapse after transplant have poor outcomes. In order to understand the contribution of tumor clonal evolution to disease progression,we applied exome and error-corrected targeted sequencing coupled with copy number analysis to comprehensively define changes in the clonal architecture of MDS in response to therapy using 51 serially acquired tumor samples from 9 patients who progressed after an alloHCT. We show that small subclones before alloHCT can drive progression after alloHCT. Notably, at least one subclone expanded or emerged at progression in all patients. Newly acquired structural variants (SVs) were present in an emergent/expanding subclone in 8 of 9 patients at progression, implicating the acquisition of SVs as important late subclonal progression events. In addition, pretransplant therapy with azacitidine likely influenced the mutation spectrum and evolution of emergent subclones after alloHCT. Although subclone evolution is common, founding clone mutations are always present at progression and could be detected in the bone marrow as early as 30 and/or 100 days after alloHCT in 6 of 8 (75%) patients, often prior to clinical progression. In conclusion, MDS progression after alloHCT is characterized by subclonal expansion and evolution, which can be influenced by pretransplant therapy.

Cellular stressors contribute to the expansion of hematopoietic clones of varying leukemic potential

Hematopoietic clones harboring specific mutations may expand over time. However, it remains unclear how different cellular stressors influence this expansion. Here we characterize clonal hematopoiesis after two different cellular stressors: cytotoxic therapy and hematopoietic transplantation. Cytotoxic therapy results in the expansion of clones carrying mutations in DNA damage response genes, including TP53 and PPM1D. Analyses of sorted populations show that these clones are typically multilineage and myeloid-biased. Following autologous transplantation, most clones persist with stable chimerism. However, DNMT3A mutant clones often expand, while PPM1D mutant clones often decrease in size. To assess the leukemic potential of these expanded clones, we genotyped 134 t-AML/t-MDS samples. Mutations in non-TP53 DNA damage response genes are infrequent in t-AML/t-MDS despite several being commonly identified after cytotoxic therapy. These data suggest that different hematopoietic stressors promote the expansion of distinct long-lived clones, carrying specific mutations, whose leukemic potential depends partially on the mutations they harbor.

Cellular stressors can impact clonal hematopoiesis. Here, the authors explore the impact of cytotoxic therapy and hematopoietic transplantation on clonal expansion, suggesting different stressors can promote expansion of distinct long-lived clones.

Abstract IA11: Spliceosome gene mutations in MDS: Biology and potential therapeutic strategies

Somatic mutations in genes that are involved in pre-mRNA splicing (i.e., spliceosome genes) occur in ~50% of patients with myelodysplastic syndromes (MDS). Commonly mutated spliceosome genes (e.g., SF3B1, SRSF2, and U2AF1) have variable prognostic significance in patients with myeloid malignancies. For example, SRSF2 and U2AF1 mutations are associated with worse and SF3B1 mutations with better clinical outcomes compared to patients without a spliceosome gene mutation. While multiple groups have shown that RNA splicing is altered in patients with spliceosome gene mutations, the contribution of these alterations to disease pathogenesis is largely unknown. Several groups have reported that RNA splicing and hematopoiesis are altered in mutant-expressing mouse models, providing direct evidence that spliceosome gene mutations contribute to some MDS phenotypes. Spliceosome gene mutations have also been detected in the blood cells of healthy elderly individuals with clonal hematopoiesis who are at an increased risk of developing hematopoietic cancer, further implicating these mutations as early genetic drivers in MDS. Consistent with this possibility, spliceosome gene mutations tend to occur in the founding clone (i.e., present in every tumor cell in a patient) of myeloid malignancies. Collectively, these observations indicate that targeting mutant cells could have a large clinical impact. Several observations indicate that spliceosome gene mutations may create a genetic vulnerability in cells that could be exploited therapeutically. Spliceosome gene mutations are typically heterozygous and mutually exclusive of each other in patients, implying either a redundancy in pathogenic function of mutant genes or that a cell cannot tolerate two spliceosome perturbations at once. Accumulating evidence now suggests that heterozygous mutant cells require appropriate expression levels of the nonmutant residual allele or other splicing factors to maintain cell viability. Therefore, we and others hypothesize that cells harboring spliceosome mutations will have increased sensitivity to pharmacologic perturbation of the spliceosome, providing a new treatment approach for patients with these mutations. To examine this, we used sudemycin compounds that bind the SF3B1 protein to modulate pre-mRNA splicing. We show that treatment of U2AF1 mutant mice with sudemycin results in attenuation of mutant-associated hematopoietic progenitor cell expansion compared to control treated mice. Collectively, data from our group and others suggest that cells expressing mutant spliceosome genes have an increased sensitivity to pharmacologic perturbation of splicing. Clinical trials are being designed and conducted to treat spliceosome mutant hematologic cancers with splicing modulator compounds.

Citation Format: Matthew J. Walter. Spliceosome gene mutations in MDS: Biology and potential therapeutic strategies [abstract]. In: Proceedings of the Second AACR Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; May 6-9, 2017; Boston, MA. Philadelphia (PA): AACR; Clin Cancer Res 2017;23(24_Suppl):Abstract nr IA11.

CpG Island Hypermethylation Mediated by DNMT3A Is a Consequence of AML Progression

DNMT3A mutations occur in ∼25% of acute myeloid leukemia (AML) patients. The most common mutation, DNMT3A^R882H, has dominant negative activity that reduces DNA methylation activity by ∼80% in vitro. To understand the contribution of DNMT3A-dependent methylation to leukemogenesis, we performed whole-genome bisulfite sequencing of primary leukemic and non-leukemic cells in patients with or without DNMT3A^R882 mutations. Non-leukemic hematopoietic cells with DNMT3A^R882H displayed focal methylation loss, suggesting that hypomethylation antedates AML. Although virtually all AMLs with wild-type DNMT3A displayed CpG island hypermethylation, this change was not associated with gene silencing and was essentially absent in AMLs with DNMT3A^R882 mutations. Primary hematopoietic stem cells expanded with cytokines were hypermethylated in a DNMT3A-dependent manner, suggesting that hypermethylation may be a response to, rather than a cause of, cellular proliferation. Our findings suggest that hypomethylation is an initiating phenotype in AMLs with DNMT3A^R882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

Knockdown of HSPA9 induces TP53-dependent apoptosis in human hematopoietic progenitor cells

Myelodysplastic syndromes (MDS) are the most common adult myeloid blood cancers in the US. Patients have increased apoptosis in their bone marrow cells leading to low peripheral blood counts. The full complement of gene mutations that contribute to increased apoptosis in MDS remains unknown. Up to 25% of MDS patients harbor and acquired interstitial deletion on the long arm of chromosome 5 [del(5q)], creating haploinsufficiency for a large set of genes including HSPA9. Knockdown of HSPA9 in primary human CD34+ hematopoietic progenitor cells significantly inhibits growth and increases apoptosis. We show here that HSPA9 knockdown is associated with increased TP53 expression and activity, resulting in increased expression of target genes BAX and p21. HSPA9 protein interacts with TP53 in CD34+ cells and knockdown of HSPA9 increases nuclear TP53 levels, providing a possible mechanism for regulation of TP53 by HSPA9 haploinsufficiency in hematopoietic cells. Concurrent knockdown of TP53 and HSPA9 rescued the increased apoptosis observed in CD34+ cells following knockdown of HSPA9. Reduction of HSPA9 below 50% results in severe inhibition of cell growth, suggesting that del(5q) cells may be preferentially sensitive to further reductions of HSPA9 below 50%, thus providing a genetic vulnerability to del(5q) cells. Treatment of bone marrow cells with MKT-077, an HSPA9 inhibitor, induced apoptosis in a higher percentage of cells from MDS patients with del(5q) compared to non-del(5q) MDS patients and normal donor cells. Collectively, these findings indicate that reduced levels of HSPA9 may contribute to TP53 activation and increased apoptosis observed in del(5q)-associated MDS.

TP53 and Decitabine in Acute Myeloid Leukemia and Myelodysplastic Syndromes

BACKGROUND

The molecular determinants of clinical responses to decitabine therapy in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS) are unclear.

METHODS

We enrolled 84 adult patients with AML or MDS in a single-institution trial of decitabine to identify somatic mutations and their relationships to clinical responses. Decitabine was administered at a dose of 20 mg per square meter of body-surface area per day for 10 consecutive days in monthly cycles. We performed enhanced exome or gene-panel sequencing in 67 of these patients and serial sequencing at multiple time points to evaluate patterns of mutation clearance in 54 patients. An extension cohort included 32 additional patients who received decitabine in different protocols.

RESULTS

Of the 116 patients, 53 (46%) had bone marrow blast clearance (<5% blasts). Response rates were higher among patients with an unfavorable-risk cytogenetic profile than among patients with an intermediate-risk or favorable-risk cytogenetic profile (29 of 43 patients [67%] vs. 24 of 71 patients [34%], P<0.001) and among patients with TP53 mutations than among patients with wild-type TP53 (21 of 21 [100%] vs. 32 of 78 [41%], P<0.001). Previous studies have consistently shown that patients with an unfavorable-risk cytogenetic profile and TP53 mutations who receive conventional chemotherapy have poor outcomes. However, in this study of 10-day courses of decitabine, neither of these risk factors was associated with a lower rate of overall survival than the rate of survival among study patients with intermediate-risk cytogenetic profiles.

CONCLUSIONS

Patients with AML and MDS who had cytogenetic abnormalities associated with unfavorable risk, TP53 mutations, or both had favorable clinical responses and robust (but incomplete) mutation clearance after receiving serial 10-day courses of decitabine. Although these responses were not durable, they resulted in rates of overall survival that were similar to those among patients with AML who had an intermediate-risk cytogenetic profile and who also received serial 10-day courses of decitabine. (Funded by the National Cancer Institute and others; ClinicalTrials.gov number, NCT01687400 .).

Comprehensive genomic analysis reveals FLT3 activation and a therapeutic strategy for a patient with relapsed adult B-lymphoblastic leukemia

The genomic events responsible for the pathogenesis of relapsed adult B-lymphoblastic leukemia (B-ALL) are not yet clear. We performed integrative analysis of whole-genome, whole-exome, custom capture, whole-transcriptome (RNA-seq), and locus-specific genomic assays across nine time points from a patient with primary de novo B-ALL. Comprehensive genome and transcriptome characterization revealed a dramatic tumor evolution during progression, yielding a tumor with complex clonal architecture at second relapse. We observed and validated point mutations in EP300 and NF1, a highly expressed EP300-ZNF384 gene fusion, a microdeletion in IKZF1, a focal deletion affecting SETD2, and large deletions affecting RB1, PAX5, NF1, and ETV6. Although the genome analysis revealed events of potential biological relevance, no clinically actionable treatment options were evident at the time of the second relapse. However, transcriptome analysis identified aberrant overexpression of the targetable protein kinase encoded by the FLT3 gene. Although the patient had refractory disease after salvage therapy for the second relapse, treatment with the FLT3 inhibitor sunitinib rapidly induced a near complete molecular response, permitting the patient to proceed to a matched-unrelated donor stem cell transplantation. The patient remains in complete remission more than 4 years later. Analysis of this patient's relapse genome revealed an unexpected, actionable therapeutic target that led to a specific therapy associated with a rapid clinical response. For some patients with relapsed or refractory cancers, this approach may indicate a novel therapeutic intervention that could alter outcome.

Implications of Tumor Clonal Heterogeneity in the Era of Next-Generation Sequencing

Recent whole-genome sequencing (WGS) studies have demonstrated that tumors typically comprise a founding clone and multiple subclones (i.e., clonal heterogeneity is common). The possible combination of mutations in each tumor clone is enormous, making each tumor genetically unique. Clonal heterogeneity likely has a role in cancer progression, relapse, metastasis, and chemoresistance due to functional differences in genetically unique subclones. In current clinical practice, gene mutations are only classified as being present or absent, ignoring the clonal complexity of cancers. In this review, we address how tumor clonality is measured using next-generation sequencing (NGS) data, highlight that clonal heterogeneity is common across multiple tumor types, and discuss the potential clinical implications of tumor clonal heterogeneity.

Abstract B41: Mutant U2AF1 alters hematopoiesis and pre-mRNA splicing in transgenic mice

This abstract is being presented as a short talk. A full abstract is printed in the Proffered Abstracts section (PR06) of the Conference Proceedings.

Citation Format: Cara Lunn Shirai, James N. Ley, Brian White, Justin Tibbitts, Jin Shao, Matthew Ndonwi, Sanghyun Kim, Theresa Okeyo-Owuor, Timothy A. Graubert, Matthew J. Walter. Mutant U2AF1 alters hematopoiesis and pre-mRNA splicing in transgenic mice. [abstract]. In: Proceedings of the AACR Special Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; Sep 20-23, 2014; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(17 Suppl):Abstract nr B41.

Abstract PR06: Mutant U2AF1 alters hematopoiesis and pre-mRNA splicing in transgenic mice

Our group and others discovered recurrent heterozygous mutations in U2AF1, a gene encoding a splicing factor involved in intronic 3'-splice site recognition, in 11% of patients with myelodysplastic syndromes (MDS), suggesting that perturbations in pre-mRNA splicing play a role in MDS pathogenesis. The most common U2AF1 mutation causes a serine to phenylalanine substitution at amino acid 34 (S34F). To study the effects of U2AF1(S34F) expression on hematopoiesis and pre-mRNA splicing in vivo, we created site-specific, single-copy, doxycycline-inducible U2AF1(WT) and U2AF1(S34F) transgenic mice. To examine the cell-autonomous effects of mutant U2AF1(S34F), we transplanted transgenic donor bone marrow into wild type mouse recipients prior to induction of transgene expression.

Following 4 weeks of transgene induction, U2AF1(S34F)-recipient mice have reduced total WBCs in the peripheral blood compared to U2AF1(WT)- and rtTA only-recipient controls (4.3 vs 7.11 and 7.13 K/µl, respectively, p≤0.01), but no significant changes in bone marrow cellularity or spleen size (n=9-11). U2AF1(S34F)-recipient mice have a perturbed mature cell lineage distribution, including reduced monocytes and B cells in both peripheral blood (p≤0.05) and bone marrow (p≤0.01) when compared to control mice (n=9-11). Reduction of bone marrow monocytes occurs as early as 5 days and is associated with increased frequency of Annexin V+ staining (p≤0.05) and phospho-H2AX (p≤0.05) compared to controls, suggesting loss of these cells may be due to apoptosis. In addition, U2AF1(S34F) expression results in alterations of progenitor and stem cell frequency in bone marrow and spleen of mice. U2AF1(S34F)-recipient mice have increased numbers of progenitors in both bone marrow and spleen by CFU-C methylcellulose assay and flow cytometry for c-Kit+ Lineage - (KL) cells, as well as common myeloid progenitors (CMPs), when compared to U2AF1(WT) and rtTA only controls (p≤0.05, n=5-10). U2AF1(S34F)-recipient mice also have an increase in the frequency of bone marrow hematopoietic stem cells (HSCs) measured by flow cytometry for KLS (c-Kit+ Lineage- Sca-1+) cells, short-term HSCs, and dormant long-term HSCs (p≤0.05), as well as an increase in spleen dormant long-term HSCs (p≤0.03) compared to U2AF1(WT)- and rtTA only-recipient controls (n=5-6). The increase in bone marrow KLS cells in U2AF1(S34F)-recipient mice at 5 days is associated with higher levels of intracellular Ki67 (a marker of cell proliferation) in KLS cells compared to controls (p&lt;0.05, n=8-13). Competitive repopulation studies show a disadvantage for bone marrow cells expressing mutant U2AF1(S34F) compared to controls (p≤0.05, n=6-12), suggesting that the increase in KLS cell cycling following U2AF1(S34F) expression may lead to stem cell exhaustion. Collectively, these data indicate U2AF1(S34F) expression alters hematopoiesis in vivo.

Next, we performed unbiased RNA sequencing on sorted bone marrow CMPs following 5 days transgene induction in U2AF1(S34F)- and U2AF1(WT)-transplanted mice. We identified 402 splicing junctions that were differentially expressed in U2AF1(S34F) samples compared to U2AF1(WT) (FDR &lt;10%, &gt;2 |log2| fold change). We also observed a preference of the mutant U2AF1(S34F) to skip exons (n=72) and alternative splice sites (n=45) with a T in the -3 position relative to the AG splice acceptor site. We intersected our junction results with RNA sequencing from AML patient samples with U2AF1 mutations and human CD34+ cells expressing U2AF1(S34F), and we identified homologous dysregulated junctions in all 3 datasets that are currently being tested for their functional significance. Together, these results suggest that mutant U2AF1 expression contributes to the altered pre-mRNA splicing and hematopoiesis observed in patients with U2AF1 mutations.

Citation Format: Cara Lunn Shirai, James N. Ley, Brian White, Justin Tibbitts, Jin Shao, Matthew Ndonwi, Sanghyun Kim, Theresa Okeyo-Owuor, Timothy A. Graubert, Matthew J. Walter. Mutant U2AF1 alters hematopoiesis and pre-mRNA splicing in transgenic mice. [abstract]. In: Proceedings of the AACR Special Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; Sep 20-23, 2014; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(17 Suppl):Abstract nr PR06.

Association Between Mutation Clearance After Induction Therapy and Outcomes in Acute Myeloid Leukemia

IMPORTANCE

Tests that predict outcomes for patients with acute myeloid leukemia (AML) are imprecise, especially for those with intermediate risk AML.

OBJECTIVES

To determine whether genomic approaches can provide novel prognostic information for adult patients with de novo AML.

DESIGN, SETTING, AND PARTICIPANTS

Whole-genome or exome sequencing was performed on samples obtained at disease presentation from 71 patients with AML (mean age, 50.8 years) treated with standard induction chemotherapy at a single site starting in March 2002, with follow-up through January 2015. In addition, deep digital sequencing was performed on paired diagnosis and remission samples from 50 patients (including 32 with intermediate-risk AML), approximately 30 days after successful induction therapy. Twenty-five of the 50 were from the cohort of 71 patients, and 25 were new, additional cases.

EXPOSURES

Whole-genome or exome sequencing and targeted deep sequencing. Risk of identification based on genetic data.

MAIN OUTCOMES AND MEASURES

Mutation patterns (including clearance of leukemia-associated variants after chemotherapy) and their association with event-free survival and overall survival.

RESULTS

Analysis of comprehensive genomic data from the 71 patients did not improve outcome assessment over current standard-of-care metrics. In an analysis of 50 patients with both presentation and documented remission samples, 24 (48%) had persistent leukemia-associated mutations in at least 5% of bone marrow cells at remission. The 24 with persistent mutations had significantly reduced event-free and overall survival vs the 26 who cleared all mutations. Patients with intermediate cytogenetic risk profiles had similar findings. [table: see text].

CONCLUSIONS AND RELEVANCE

The detection of persistent leukemia-associated mutations in at least 5% of bone marrow cells in day 30 remission samples was associated with a significantly increased risk of relapse, and reduced overall survival. These data suggest that this genomic approach may improve risk stratification for patients with AML.

Mutant U2AF1 Expression Alters Hematopoiesis and Pre-mRNA Splicing In Vivo

Heterozygous somatic mutations in the spliceosome gene U2AF1 occur in ∼ 11% of patients with myelodysplastic syndromes (MDS), the most common adult myeloid malignancy. It is unclear how these mutations contribute to disease. We examined in vivo hematopoietic consequences of the most common U2AF1 mutation using a doxycycline-inducible transgenic mouse model. Mice expressing mutant U2AF1(S34F) display altered hematopoiesis and changes in pre-mRNA splicing in hematopoietic progenitor cells by whole transcriptome analysis (RNA-seq). Integration with human RNA-seq datasets determined that common mutant U2AF1-induced splicing alterations are enriched in RNA processing genes, ribosomal genes, and recurrently mutated MDS and acute myeloid leukemia-associated genes. These findings support the hypothesis that mutant U2AF1 alters downstream gene isoform expression, thereby contributing to abnormal hematopoiesis in patients with MDS.

Role of TP53 mutations in the origin and evolution of therapy-related acute myeloid leukaemia

Therapy-related acute myeloid leukaemia (t-AML) and therapy-related myelodysplastic syndrome (t-MDS) are well-recognized complications of cytotoxic chemotherapy and/or radiotherapy. There are several features that distinguish t-AML from de novo AML, including a higher incidence of TP53 mutations, abnormalities of chromosomes 5 or 7, complex cytogenetics and a reduced response to chemotherapy. However, it is not clear how prior exposure to cytotoxic therapy influences leukaemogenesis. In particular, the mechanism by which TP53 mutations are selectively enriched in t-AML/t-MDS is unknown. Here, by sequencing the genomes of 22 patients with t-AML, we show that the total number of somatic single-nucleotide variants and the percentage of chemotherapy-related transversions are similar in t-AML and de novo AML, indicating that previous chemotherapy does not induce genome-wide DNA damage. We identified four cases of t-AML/t-MDS in which the exact TP53 mutation found at diagnosis was also present at low frequencies (0.003-0.7%) in mobilized blood leukocytes or bone marrow 3-6 years before the development of t-AML/t-MDS, including two cases in which the relevant TP53 mutation was detected before any chemotherapy. Moreover, functional TP53 mutations were identified in small populations of peripheral blood cells of healthy chemotherapy-naive elderly individuals. Finally, in mouse bone marrow chimaeras containing both wild-type and Tp53(+/-) haematopoietic stem/progenitor cells (HSPCs), the Tp53(+/-) HSPCs preferentially expanded after exposure to chemotherapy. These data suggest that cytotoxic therapy does not directly induce TP53 mutations. Rather, they support a model in which rare HSPCs carrying age-related TP53 mutations are resistant to chemotherapy and expand preferentially after treatment. The early acquisition of TP53 mutations in the founding HSPC clone probably contributes to the frequent cytogenetic abnormalities and poor responses to chemotherapy that are typical of patients with t-AML/t-MDS.