BACKGROUND: Natural killer (NK) cell-based therapies represent a promising approach for acute myeloid leukemia (AML) relapse, yet their efficacy is hindered by immunosuppressive factors such as transforming growth factor beta (TGF-β) in the tumor microenvironment. This study investigated the effects of TGF-β on NK cell cytotoxicity and migration using 2D and 3D co-culture models that mimic the leukemic microenvironment. METHODS: TGF-β production was evaluated in AML-derived leukemic cell lines and mesenchymal stromal cells (hTERT-MSCs) using ELISA. Bulk RNA sequencing (RNA-seq) was performed to analyze global gene expression changes in TGF-β-treated primary human NK cells. NK cell cytotoxicity and migration were assessed in 2D monolayer and 3D spheroid co-cultures containing hTERT-MSCs and leukemic cells using flow cytometry and confocal microscopy. RESULTS: Both leukemic cells and MSCs produced TGF-β, with increased levels observed in MSCs after co-culture with primary AML blasts. RNA sequencing revealed that TGF-β altered key gene pathways associated with NK cell cytotoxicity, adhesion, and migration, supporting its immunosuppressive role. In functional assays, TGF-β exposure significantly reduced NK cell-mediated cytotoxicity in a time-dependent manner and impaired NK cell infiltration into 3D spheroids, particularly in models incorporating MSCs. Additionally, MSCs themselves provided a protective environment for leukemic cells, further reducing NK cell effectiveness in 2D co-cultures. CONCLUSION: TGF-β suppresses both NK cell cytotoxicity and migration, limiting their ability to eliminate leukemic cells and infiltrate the bone marrow niche (BMN). These findings provide novel insights into TGF-β-mediated immune evasion mechanisms and provide important insights for the future design of NK-based immunotherapies and clinical trials.
- Publication type
- Journal Article MeSH
BACKGROUND: Myelodysplastic neoplasms (MDS) are heterogeneous hematopoietic disorders characterized by ineffective hematopoiesis and genome instability. Mobilization of transposable elements (TEs) is an important source of genome instability leading to oncogenesis, whereas small PIWI-interacting RNAs (piRNAs) act as cellular suppressors of TEs. However, the roles of TEs and piRNAs in MDS remain unclear. METHODS: In this study, we examined TE and piRNA expression through parallel RNA and small RNA sequencing of CD34+ hematopoietic stem cells from MDS patients. RESULTS: Comparative analysis of TE and piRNA expression between MDS and control samples revealed several significantly dysregulated molecules. However, significant differences were observed between lower-risk MDS (LR-MDS) and higher-risk MDS (HR-MDS) samples. In HR-MDS, we found an inverse correlation between decreased TE levels and increased piRNA expression and these TE and piRNA levels were significantly associated with patient outcomes. Importantly, the upregulation of PIWIL2, which encodes a key factor in the piRNA pathway, independently predicted poor prognosis in MDS patients, underscoring its potential as a valuable disease marker. Furthermore, pathway analysis of RNA sequencing data revealed that dysregulation of the TE‒piRNA axis is linked to the suppression of processes related to energy metabolism, the cell cycle, and the immune response, suggesting that these disruptions significantly affect cellular activity. CONCLUSIONS: Our findings demonstrate the parallel dysregulation of TEs and piRNAs in HR-MDS patients, highlighting their potential role in MDS progression and indicating that the PIWIL2 level is a promising molecular marker for prognosis.
- Publication type
- Journal Article MeSH
Patients with myelodysplastic neoplasms (MDS) are classified according to the risk of acute myeloid leukemia transformation. Some lower-risk MDS patients (LR-MDS) progress rapidly despite expected good prognosis. Using diagnostic samples, we aimed to uncover the mechanisms of this accelerated progression at the transcriptome level. RNAseq was performed on CD34+ ribodepleted RNA samples from 53 LR-MDS patients without accelerated progression (stMDS) and 8 who progressed within 20 months (prMDS); 845 genes were differentially expressed (ІlogFCІ > 1, FDR < 0.01) between these groups. stMDS CD34+ cells exhibited transcriptional signatures of actively cycling, megakaryocyte/erythrocyte lineage-primed progenitors, with upregulation of cell cycle checkpoints and stress pathways, which presumably form a tumor-suppressing barrier. Conversely, cell cycle, DNA damage response (DDR) and energy metabolism-related pathways were downregulated in prMDS samples, whereas cell adhesion processes were upregulated. Also, prMDS samples showed high levels of aberrant splicing and global lncRNA expression that may contribute to the attenuation of DDR pathways. We observed overexpression of multiple oncogenes and diminished differentiation in prMDS; the expression of ZEB1 and NEK3, genes not previously associated with MDS prognosis, might serve as potential biomarkers for LR-MDS progression. Our 19-gene DDR signature showed a significant predictive power for LR-MDS progression. In validation samples (stMDS = 3, prMDS = 4), the key markers and signatures retained their significance. Collectively, accelerated progression of LR-MDS appears to be associated with transcriptome patterns of a quiescent-like cell state, reduced lineage differentiation and suppressed DDR, inherent to CD34+ cells. The attenuation of DDR-related gene-expression signature may refine risk assessment in LR-MDS patients.
- MeSH
- Cell Adhesion MeSH
- Cell Cycle MeSH
- NIMA-Related Kinases genetics metabolism MeSH
- Humans MeSH
- Myelodysplastic Syndromes * genetics MeSH
- Neoplasms * MeSH
- DNA Repair MeSH
- Transcriptome MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
Mutations in the splicing factor 3b subunit 1 (SF3B1) gene are frequent in myelodysplastic neoplasms (MDS). Because the splicing process is involved in the production of circular RNAs (circRNAs), we investigated the impact of SF3B1 mutations on circRNA processing. Using RNA sequencing, we measured circRNA expression in CD34+ bone marrow MDS cells. We defined circRNAs deregulated in a heterogeneous group of MDS patients and described increased circRNA formation in higher-risk MDS. We showed that the presence of SF3B1 mutations did not affect the global production of circRNAs; however, deregulation of specific circRNAs was observed. Particularly, we demonstrated that strong upregulation of circRNAs processed from the zinc finger E-box binding homeobox 1 (ZEB1) transcription factor; this upregulation was exclusive to SF3B1-mutated patients and was not observed in those with mutations in other splicing factors or other recurrently mutated genes, or with other clinical variables. Furthermore, we focused on the most upregulated ZEB1-circRNA, hsa_circ_0000228, and, by its knockdown, we demonstrated that its expression is related to mitochondrial activity. Using microRNA analyses, we proposed miR-1248 as a direct target of hsa_circ_0000228. To conclude, we demonstrated that mutated SF3B1 leads to deregulation of ZEB1-circRNAs, potentially contributing to the defects in mitochondrial metabolism observed in SF3B1-mutated MDS.
- MeSH
- Leukemia, Myeloid, Acute * MeSH
- Phosphoproteins genetics MeSH
- RNA, Circular genetics MeSH
- Humans MeSH
- Mutation genetics MeSH
- Myelodysplastic Syndromes * genetics MeSH
- RNA Splicing Factors genetics MeSH
- Transcription Factors genetics MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
BACKGROUND/AIM: Prediction of response to azacitidine (AZA) treatment is an important challenge in hematooncology. In addition to protein coding genes (PCGs), AZA efficiency is influenced by various noncoding RNAs (ncRNAs), including long ncRNAs (lncRNAs), circular RNAs (circRNAs), and transposable elements (TEs). MATERIALS AND METHODS: RNA sequencing was performed in patients with myelodysplastic syndromes or acute myeloid leukemia before AZA treatment to assess contribution of ncRNAs to AZA mechanisms and propose novel disease prediction biomarkers. RESULTS: Our analyses showed that lncRNAs had the strongest predictive potential. The combined set of the best predictors included 14 lncRNAs, and only four PCGs, one circRNA, and no TEs. Epigenetic regulation and recombinational repair were suggested as crucial for AZA response, and network modeling defined three deregulated lncRNAs (CTC-482H14.5, RP11-419K12.2, and RP11-736I24.4) associated with these processes. CONCLUSION: The expression of various ncRNAs can influence the effect of AZA and new ncRNA-based predictive biomarkers can be defined.
