Subunit-specific analysis of cohesin-mutant myeloid malignancies reveals distinct ontogeny and outcomes.

Mutations in the cohesin complex components (STAG2, RAD21, SMC1A, SMC3, and PDS5B) are recurrent genetic drivers in myelodysplastic neoplasm (MDS) and acute myeloid leukemia (AML). Whether the different cohesin subunit mutations share clinical characteristics and prognostic significance is not known. We analyzed 790 cohesin-mutant patients from the Dana-Farber Cancer Institute (DFCI) and the Munich Leukemia Laboratory (MLL), 390 of which had available outcome data, and identified subunit-specific clinical, prognostic, and genetic characteristics suggestive of distinct ontogenies. We found that STAG2 mutations are acquired at MDS stage and are associated with secondary AML, adverse prognosis, and co-occurrence of secondary AML-type mutations. In contrast, mutations in RAD21, SMC1A and SMC3 share features with de novo AML with better prognosis, and co-occurrence with de novo AML-type lesions. The findings show the heterogeneous nature of cohesin complex mutations, and inform clinical and prognostic classification, as well as distinct biology of the cohesin complex.

Cell-type-specific consequences of mosaic structural variants in hematopoietic stem and progenitor cells.

The functional impact and cellular context of mosaic structural variants (mSVs) in normal tissues is understudied. Utilizing Strand-seq, we sequenced 1,133 single-cell genomes from 19 human donors of increasing age, and discovered the heterogeneous mSV landscapes of hematopoietic stem and progenitor cells. While mSVs are continuously acquired throughout life, expanded subclones in our cohort are confined to individuals >60. Cells already harboring mSVs are more likely to acquire additional somatic structural variants, including megabase-scale segmental aneuploidies. Capitalizing on comprehensive single-cell micrococcal nuclease digestion with sequencing reference data, we conducted high-resolution cell-typing for eight hematopoietic stem and progenitor cells. Clonally expanded mSVs disrupt normal cellular function by dysregulating diverse cellular pathways, and enriching for myeloid progenitors. Our findings underscore the contribution of mSVs to the cellular and molecular phenotypes associated with the aging hematopoietic system, and establish a foundation for deciphering the molecular links between mSVs, aging and disease susceptibility in normal tissues.

Treatment with the apoptosis inhibitor Asunercept reduces clone sizes in patients with lower risk Myelodysplastic Neoplasms.

In low-risk Myelodysplastic Neoplasms (MDS), increased activity of apoptosis-promoting factors such as tumor necrosis factor (TNFα) and pro-apoptotic Fas ligand (CD95L) have been described as possible pathomechanisms leading to impaired erythropoiesis. Asunercept (APG101) is a novel therapeutic fusion protein blocking CD95, which has previously shown partial efficacy in reducing transfusion requirement in a clinical phase I trial for low-risk MDS patients (NCT01736436; 2012-11-26). In the current study we aimed to evaluate the effect of Asunercept therapy on the clonal bone marrow composition to identify potential biomarkers to predict response. Bone marrow samples of n = 12 low-risk MDS patients from the above referenced clinical trial were analyzed by serial deep whole exome sequencing in a total of n = 58 time points. We could distinguish a mean of 3.5 molecularly defined subclones per patient (range 2-6). We observed a molecular response defined as reductions of dominant clone sizes by a variant allele frequency (VAF) decrease of at least 10% (mean 20%, range: 10.5-39.2%) in dependency of Asunercept treatment in 9 of 12 (75%) patients. Most of this decline in clonal populations was observed after completion of 12 weeks treatment. Particularly early and pronounced reductions of clone sizes were found in subclones driven by mutations in genes involved in regulation of methylation (n = 1 DNMT3A, n = 1 IDH2, n = 1 TET2). Our results suggest that APG101 could be efficacious in reducing clone sizes of mutated hematopoietic cells in the bone marrow of Myelodysplastic Neoplasms, which warrants further investigation.

ASXL1 mutations are associated with a response to alvocidib and 5-azacytidine combination in myelodysplastic neoplasms.

Inhibitors of anti-apoptotic BCL-2 family proteins in combination with chemotherapy and hypomethylating agents (HMAs) are promising therapeutic approaches in acute myeloid leukemia (AML) and high-risk myelodysplastic syndromes (MDS). Alvocidib, a cyclin-dependent kinase 9 (CDK9) inhibitor and indirect transcriptional repressor of the anti-apoptotic factor MCL-1, has previously shown clinical activity in AML. Availability of biomarkers for response to the alvocidib + 5- AZA could also extend the rationale of this treatment concept to high-risk MDS. In this study, we performed a comprehensive in vitro assessment of alvocidib and 5-AZA effects in n=45 high-risk MDS patients. Our data revealed additive cytotoxic effects of the combination treatment. Mutational profiling of MDS samples identified ASXL1 mutations as predictors of response. Further, increased response rates were associated with higher gene-expression of the pro-apoptotic factor NOXA in ASXL1 mutated samples. The higher sensitivity of ASXL1 mutant cells to the combination treatment was confirmed in vivo in ASXL1Y588X transgenic mice. Overall, our study demonstrated augmented activity for the alvocidib + 5-AZA combination in higher-risk MDS and identified ASXL1 mutations as a biomarker of response for potential stratification studies.

Significant improvement of bone marrow-derived MSC expansion from MDS patients by defined xeno-free medium.

BACKGROUND: Robust and reliable in vitro and in vivo models of primary cells are necessary to study the pathomechanisms of Myelodysplastic Neoplasms (MDS) and identify novel therapeutic strategies. MDS-derived hematopoietic stem and progenitor cells (HSPCs) are reliant on the support of bone marrow (BM) derived mesenchymal stroma cells (MSCs). Therefore, isolation and expansion of MCSs are essential for successfully modeling this disease. For the clinical use of healthy MSCs isolated from human BM, umbilical cord blood or adipose tissue, several studies showed that xeno-free (XF) culture conditions resulted in superior growth kinetics compared to MSCs cultured in the presence of fetal bovine serum (FBS). In this present study, we investigate, whether the replacement of a commercially available MSC expansion medium containing FBS with a XF medium is beneficial for the expansion of MSCs derived from BM of MDS patients which are often difficult to cultivate. METHODS: MSCs isolated from BM of MDS patients were cultured and expanded in MSC expansion medium with FBS or XF supplement. Subsequently, the impact of culture media on growth kinetics, morphology, immunophenotype, clonogenic potential, differentiation capacity, gene expression profiles and ability to engraft in immunodeficient mouse models was evaluated. RESULTS: Significant higher cell numbers with an increase in clonogenic potential were observed during culture of MDS MSCs with XF medium compared to medium containing FBS. Differential gene expression showed an increase in transcripts associated with MSC stemness after expansion with XF. Furthermore, immunophenotypes of the MSCs and their ability to differentiate into osteoblasts, adipocytes or chondroblasts remained stable. MSCs expanded with XF media were similarly supportive for creating MDS xenografts in vivo as MSCs expanded with FBS. CONCLUSION: Our data indicate that with XF media, higher cell numbers of MDS MSCs can be obtained with overall improved characteristics in in vitro and in vivo experimental models.

Inhibition of lysyl oxidases synergizes with 5-azacytidine to restore erythropoiesis in myelodysplastic and myeloid malignancies.

Limited response rates and frequent relapses during standard of care with hypomethylating agents in myelodysplastic neoplasms (MN) require urgent improvement of this treatment indication. Here, by combining 5-azacytidine (5-AZA) with the pan-lysyl oxidase inhibitor PXS-5505, we demonstrate superior restoration of erythroid differentiation in hematopoietic stem and progenitor cells (HSPCs) of MN patients in 20/31 cases (65%) versus 9/31 cases (29%) treated with 5-AZA alone. This effect requires direct contact of HSPCs with bone marrow stroma components and is dependent on integrin signaling. We further confirm these results in vivo using a bone marrow niche-dependent MN xenograft model in female NSG mice, in which we additionally demonstrate an enforced reduction of dominant clones as well as significant attenuation of disease expansion and normalization of spleen sizes. Overall, these results lay out a strong pre-clinical rationale for efficacy of combination treatment of 5-AZA with PXS-5505 especially for anemic MN.

Functional analysis of structural variants in single cells using Strand-seq.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Humanized three-dimensional scaffold xenotransplantation models for myelodysplastic syndromes.

Patient-derived xenograft (PDX) models have emerged as versatile preclinical platforms for investigation of functional pathomechanisms in myelodysplastic syndromes (MDS) and other myeloid neoplasms. However, despite increasingly improved methodology, engraftment efficiencies frequently remain low. Humanized three-dimensional scaffold models (ossicle xenotransplantation models) in immunocompromised mice have recently been found to enable improved engraftment rates of healthy and malignant human hematopoiesis. We therefore interrogated the feasibility of using four different three-dimensional ossicle-based PDX models for application with primary MDS samples. In a fully standardized comparison, we evaluated scaffold materials such as Gelfoam, extracellular matrix (ECM), and human or xenogenous bone substance in comparison to intrafemoral (IF) co-injection of bone marrow (BM)-derived mesenchymal stromal cells (MSCs) and CD34+ hematopoietic stem and progenitor cells (HSPCs). Our study included13 primary MDS patient samples transplanted in parallel according to these five different conditions. Engraftment of MDS samples was assessed by flow cytometry, immunohistological staining, and molecular validation. We determined that three-dimensional ossicle-based methods achieved higher relative rates of engraftment and enabled long-term retrievability of patient-derived MSCs from implanted ossicles. In summary, HSPCs and MSCs derived from MDS BM, which did not significantly engraft in NSG mice after intrafemoral injection, were able to colonize humanized scaffold models. Therefore, these models are promising new xenotransplantation techniques for addressing preclinical and functional questions of the interaction between hematopoiesis and the BM niche in MDS.

Preclinical evaluation of eltrombopag in a PDX model of myelodysplastic syndromes.

Preclinical research of myelodysplastic syndromes (MDSs) is hampered by a lack of feasible disease models. Previously, we have established a robust patient-derived xenograft (PDX) model for MDS. Here we demonstrate for the first time that this model is applicable as a preclinical platform to address pending clinical questions by interrogating the efficacy and safety of the thrombopoietin receptor agonist eltrombopag. Our preclinical study included n = 49 xenografts generated from n = 9 MDS patient samples. Substance efficacy was evidenced by FACS-based human platelet quantification and clonal bone marrow evolution was reconstructed by serial whole-exome sequencing of the PDX samples. In contrast to clinical trials in humans, this experimental setup allowed vehicle- and replicate-controlled analyses on a patient-individual level deciphering substance-specific effects from natural disease progression. We found that eltrombopag effectively stimulated thrombopoiesis in MDS PDX without adversely affecting the patients' clonal composition. In conclusion, our MDS PDX model is a useful tool for testing new therapeutic concepts in MDS preceding clinical trials.

Single-cell proteo-genomic reference maps of the hematopoietic system enable the purification and massive profiling of precisely defined cell states.

Single-cell genomics technology has transformed our understanding of complex cellular systems. However, excessive cost and a lack of strategies for the purification of newly identified cell types impede their functional characterization and large-scale profiling. Here, we have generated high-content single-cell proteo-genomic reference maps of human blood and bone marrow that quantitatively link the expression of up to 197 surface markers to cellular identities and biological processes across all main hematopoietic cell types in healthy aging and leukemia. These reference maps enable the automatic design of cost-effective high-throughput cytometry schemes that outperform state-of-the-art approaches, accurately reflect complex topologies of cellular systems and permit the purification of precisely defined cell states. The systematic integration of cytometry and proteo-genomic data enables the functional capacities of precisely mapped cell states to be measured at the single-cell level. Our study serves as an accessible resource and paves the way for a data-driven era in cytometry.

Bone marrow derived stromal cells from myelodysplastic syndromes are altered but not clonally mutated in vivo.

The bone marrow (BM) stroma in myeloid neoplasms is altered and it is hypothesized that this cell compartment may also harbor clonal somatically acquired mutations. By exome sequencing of in vitro expanded mesenchymal stromal cells (MSCs) from n = 98 patients with myelodysplastic syndrome (MDS) and n = 28 healthy controls we show that these cells accumulate recurrent mutations in genes such as ZFX (n = 8/98), RANK (n = 5/98), and others. MDS derived MSCs display higher mutational burdens, increased replicative stress, senescence, inflammatory gene expression, and distinct mutational signatures as compared to healthy MSCs. However, validation experiments in serial culture passages, chronological BM aspirations and backtracking of high confidence mutations by re-sequencing primary sorted MDS MSCs indicate that the discovered mutations are secondary to in vitro expansion but not present in primary BM. Thus, we here report that there is no evidence for clonal mutations in the BM stroma of MDS patients.

Cohesin mutations in myeloid malignancies.

Cohesin is a multisubunit protein complex that forms a ring-like structure around DNA. It is essential for sister chromatid cohesion, chromatin organization, transcriptional regulation, and DNA damage repair and plays a major role in dynamically shaping the genome architecture and maintaining DNA integrity. The core complex subunits STAG2, RAD21, SMC1, and SMC3, as well as its modulators PDS5A/B, WAPL, and NIPBL, have been found to be recurrently mutated in hematologic and solid malignancies. These mutations are found across the full spectrum of myeloid neoplasia, including pediatric Down syndrome-associated acute megakaryoblastic leukemia, myelodysplastic syndromes, chronic myelomonocytic leukemia, and de novo and secondary acute myeloid leukemias. The mechanisms by which cohesin mutations act as drivers of clonal expansion and disease progression are still poorly understood. Recent studies have described the impact of cohesin alterations on self-renewal and differentiation of hematopoietic stem and progenitor cells, which are associated with changes in chromatin and epigenetic state directing lineage commitment, as well as genomic integrity. Herein, we review the role of the cohesin complex in healthy and malignant hematopoiesis. We discuss clinical implications of cohesin mutations in myeloid malignancies and discuss opportunities for therapeutic targeting.

Genome-wide DNA methylation analysis pre- and post-lenalidomide treatment in patients with myelodysplastic syndrome with isolated deletion (5q).

Myelodysplastic syndrome (MDS) with isolated deletion of chromosome 5q (MDS del5q) is a distinct subtype of MDS with quite favorable prognosis and excellent response to treatment with lenalidomide. Still, a relevant percentage of patients do not respond to lenalidomide and even experience progression to acute myeloid leukemia (AML). In this study, we aimed to investigate whether global DNA methylation patterns could predict response to lenalidomide. Genome-wide DNA methylation analysis using Illumina 450k methylation arrays was performed on n=51 patients with MDS del5q who were uniformly treated with lenalidomide in a prospective multicenter trial of the German MDS study group. To study potential direct effects of lenalidomide on DNA methylation, 17 paired samples pre- and post-treatment were analyzed. Our results revealed no relevant effect of lenalidomide on methylation status. Furthermore, methylation patterns prior to therapy could not predict lenalidomide response. However, methylation clustering identified a group of patients with a trend towards inferior overall survival. These patients showed hypermethylation of several interesting target genes, including genes of relevant signaling pathways, potentially indicating the evaluation of novel therapeutic targets.

Identification of leukemic and pre-leukemic stem cells by clonal tracking from single-cell transcriptomics.

Cancer stem cells drive disease progression and relapse in many types of cancer. Despite this, a thorough characterization of these cells remains elusive and with it the ability to eradicate cancer at its source. In acute myeloid leukemia (AML), leukemic stem cells (LSCs) underlie mortality but are difficult to isolate due to their low abundance and high similarity to healthy hematopoietic stem cells (HSCs). Here, we demonstrate that LSCs, HSCs, and pre-leukemic stem cells can be identified and molecularly profiled by combining single-cell transcriptomics with lineage tracing using both nuclear and mitochondrial somatic variants. While mutational status discriminates between healthy and cancerous cells, gene expression distinguishes stem cells and progenitor cell populations. Our approach enables the identification of LSC-specific gene expression programs and the characterization of differentiation blocks induced by leukemic mutations. Taken together, we demonstrate the power of single-cell multi-omic approaches in characterizing cancer stem cells.

High erythroferrone expression in CD71+ erythroid progenitors predicts superior survival in myelodysplastic syndromes.

Ineffective erythropoiesis and iron overload are common in myelodysplastic syndromes (MDS). Erythroferrone (ERFE) and growth/differentiation factor 15 (GDF15) are two regulators of iron homeostasis produced by erythroid progenitors. Elevated systemic levels of ERFE and GDF15 in MDS are associated with dysregulated iron metabolism and iron overload, which is especially pronounced in MDS with SF3B1 gene mutations. However, the role of ERFE and GDF15 in MDS pathogenesis and their influence on disease progression are largely unknown. Here, we analyzed the expression of ERFE and GDF15 in CD71+ erythroid progenitors of n = 111 MDS patients and assessed their effects on patient survival. The expression of ERFE and GDF15 in MDS was highly aberrant. Unexpectedly, ERFE expression in erythroprogenitors was highly relevant for MDS prognosis and independent of International Prognostic Scoring System (IPSS) stratification. Although ERFE expression was increased in patients with SF3B1 mutations, it predicted overall survival (OS) in both the SF3B1wt and SF3B1mut subgroups. Of note, ERFE overexpression predicted superior OS in the IPSS low/Int-1 subgroup and in patients with normal karyotype. Similar observations were made for GDF15, albeit not reaching statistical significance. In summary, our results revealed a strong association between ERFE expression and MDS outcome, suggesting a possible involvement of ERFE in molecular MDS pathogenesis.

Replication stress signaling is a therapeutic target in myelodysplastic syndromes with splicing factor mutations.

Somatic mutations in genes coding for splicing factors, e.g. SF3B1, U2AF1, SRSF2, and others are found in approximately 50% of patients with Myelodysplastic Syndromes (MDS). These mutations have been predicted to frequently occur early in the mutational hierarchy of the disease therefore making them particularly attractive potential therapeutic targets. Recent studies in cell lines engineered to carry splicing factor mutations have revealed a strong association with elevated levels of DNA:RNA intermediates (R-loops) and a dependency on proper ATR function. However, data confirming this hypothesis in a representative cohort of primary MDS patient samples have so far been missing. Using CD34+ cells isolated from MDS patients with and without splicing factor mutations as well as healthy controls we show that splicing factor mutation-associated R-loops lead to elevated levels of replication stress and ATR pathway activation. Moreover, splicing factor mutated CD34+ cells are more susceptible to pharmacological inhibition of ATR resulting in elevated levels of DNA damage, cell cycle blockade, and cell death. This can be enhanced by combination treatment with low-dose splicing modulatory compound Pladienolide B. We further confirm the direct association of R-loops and ATR sensitivity with the presence of a splicing factor mutation using lentiviral overexpression of wild-type and mutant SRSF2 P95H in cord blood CD34+ cells. Collectively, our results from n=53 MDS patients identify replication stress and associated ATR signaling to be critical pathophysiological mechanisms in primary MDS CD34+ cells carrying splicing factor mutations, and provide a preclinical rationale for targeting ATR signaling in these patients.

RNA-sequencing of acute promyelocytic leukemia primary blasts reveals novel molecular biomarkers of early death events.

Although acute promyelocytic leukemia (APL) has evolved to the AML entity with the best prognosis, typical 'early death' (ED) events still account for mortality rates of ∼20% in population-based studies. To investigate this poorly understood issue we performed whole transcriptome analysis of n = 7 APL ED cases compared to n = 7 APL cases with long term remission. We discovered the proteins S100A8/S100A9 and EFEMP1 as the most differentially expressed factors. In an independent cohort of n = 58 APL patients EFEMP1 over-expression was associated with a worse overall survival. Furthermore, a subgroup analysis of ED caused by hemorrhagic complications revealed an association of metallothioneins (MT1G/MT1E) with higher bleeding rates, ED events and negative prognostic effects on overall survival. Finally, we identified a novel TPM4-KLF2 fusion transcripts in 44/64 APL samples. In summary, we report a comprehensive transcriptomic analysis and novel potential biomarkers of ED biology, which highlight novel pathways in ED events in APL.

Characterization of Stem-like Circulating Tumor Cells in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) is the fourth most frequent cause of death from cancer. Circulating tumor cells (CTCs) with stem-like characteristics lead to distant metastases and thus contribute to the dismal prognosis of PDAC. Our purpose is to investigate the role of stemness in CTCs derived from a genetically engineered mouse model of PDAC and to further explore the potential molecular mechanisms. The publically available RNA sequencing dataset GSE51372 was analyzed, and CTCs with (CTC-S) or without (CTC-N) stem-like features were discriminated based on a principal component analysis (PCA). Differentially expressed genes, weighted gene co-expression network analysis (WGCNA), and further functional enrichment analyses were performed. The prognostic role of the candidate gene (CTNNB1) was assessed in a clinical PDAC patient cohort. Overexpression of the pluripotency marker Klf4 (Krüppel-like factor 4) in CTC-S cells positively correlates with Ctnnb1 (ß-Catenin) expression, and their interaction presumably happens via protein-protein binding in the nucleus. As a result, the adherens junction pathway is significantly enriched in CTC-S. Furthermore, the overexpression of Ctnnb1 is a negative prognostic factor for progression-free survival (PFS) and relapse-free survival (RFS) in human PDAC cohort. Overexpression of Ctnnb1 may thus promote the metastatic capabilities of CTCs with stem-like properties via adherens junctions in murine PDAC.

Deep sequencing of bone marrow microenvironments of patients with del(5q) myelodysplastic syndrome reveals imprints of antigenic selection as well as generation of novel T-cell clusters as a response pattern to lenalidomide.

In myelodysplastic syndromes with a partial deletion of the long arm of chromosome 5, del(5q), lenalidomide is believed to reverse anergic T-cell immunity in the bone marrow resulting in suppression of the del(5q) clone. In this study we used next-generation sequencing of immunoglobulin heavy chain (IGH) and T-cell receptor beta (TRB) rearrangements in bone marrow-residing and peripheral blood-circulating lymphocytes of patients with del(5q) myelodysplastic syndromes to assess the immune architecture and track adaptive immune responses during treatment with lenalidomide. The baseline bone marrow B-cell space in patients was comparable to that of age-matched healthy controls in terms of gene usage and IGH clonality, but showed a higher percentage of hypermutated IGH sequences, indicating an expanded number of antigen-experienced B lineage cells. Bone marrow B lineage clonality decreased significantly and hypermutated IGH clones normalized upon lenalidomide treatment, well in line with the proliferative effect on healthy antigen-inexperienced B-cell precursors previously described for this drug. The T-cell space in bone marrow of patients with del(5q) myelodysplastic syndromes showed higher TRB clonality compared to that of healthy controls. Upon lenalidomide treatment, myelodysplastic syndrome-specific T-cell clusters with low to medium spontaneous generation probabilities emerged; these clusters were shared across patients, indicating a common antigen-driven T-cell response pattern. Hence, we observed B lineage diversification and generation of new, antigen-dependent T-cell clusters, compatible with a model of adaptive immunity induced against the del(5q) clone by lenalidomide. Overall, this supports the concept that lenalidomide not only alters the functional T-cell state, but also the composition of the T- and B-cell repertoires in del(5q) myelodysplastic syndromes.