The EMT modulator SNAI1 contributes to AML pathogenesis via its interaction with LSD1

The ETV6-MECOM fusion protein promotes EMT-related properties by repressing the transactivation activity of E-cadherin promoter in K562 leukemia cells

The ETV6-MECOM fusion gene, produced by the rare and recurrent chromosomal translocation t(3; 12) (q26; p13), is associated with high mortality and short survival in myeloid leukemia. However, its function and underlying mechanisms in leukemia progression remain unknown. In this study, leukemia-stable K562 cells expressing the ETV6-MECOM fusion protein were used to investigate the effects of the ETV6-MECOM oncoprotein. K562-ETV6-MECOM cells were undifferentiated and had reduced colony formation, increased cell migration and invasion, and increased sphere number and diameter in a spheroid formation assay, presenting epithelial-to-mesenchymal transition (EMT) traits. The expression of E-cadherin, a hallmark of EMT, was significantly downregulated at the transcriptional and translational level in K562-ETV6-MECOM cells to explore the mechanistic basis of EMT. Stepwise truncation, DNA sequence deletion, mutation analysis for E-cadherin promoter transactivation, and a dual luciferase assay indicated that the regulatory region of ETV6-MECOM is located in the DNA motif -1116 TTAAAA-1111 of E-cadherin promoter. Moreover, a chromatin immunoprecipitation assay showed that this oncoprotein binds to the DNA motif -1116 TTAAAA-1111 with the anti-EVI1 antibody. Although ETV6-MECOM upregulated the expressions of EMT master regulators, including SNAIL, SLUG, ZEB2, and TWIST2, their knockdown had no effect on EMT-related properties. However, overexpression of E-cadherin eliminated EMT traits in the presence of the ETV6-MECOM oncoprotein. These data confirmed that the ETV6-MECOM oncoprotein, not SNAIL, SLUG, ZEB2, or TWIST2, plays a critical role in inducing EMT traits in leukemia K562 cells. ETV6-MECOM induces EMT-related properties by downregulating the transcriptional expression of E-cadherin and repressing its transactivation activity by binding to its core motif -1116TTAAAA-1111 in leukemia K562 cells. These findings could contribute to the development of a therapeutic target for patients with myeloid leukemia characterized by ETV6-MECOM.

A 69 long noncoding RNA signature predicts relapse and acts as independent prognostic factor in pediatric AML

ABSTRACT

Risk stratification using genetics and minimal residual disease has allowed for an increase in the cure rates of pediatric acute myeloid leukemia (pedAML) to up to 70% in contemporary protocols. Nevertheless, ∼30% of patients still experience relapse, indicating a need to optimize stratification strategies. Recently, long noncoding RNA (lncRNA) expression has been shown to hold prognostic power in multiple cancer types. Here, we aimed at refining relapse prediction in pedAML using lncRNA expression. We built a relapse-related lncRNA prognostic signature, named AMLlnc69, using 871 transcriptomes of patients with pedAML obtained from the Therapeutically Applicable Research to Generate Effective Treatments repository. We identified a 69 lncRNA signature AMLlnc69 that is highly predictive of relapse risk (c-index = 0.73), with area under the receiver operating characteristic curve (AUC) values for predicting the 1-, 2-, and 3-year relapse-free survival (RFS) of 0.78, 0.77, and 0.77, respectively. The internal validation using a bootstrap method (resampling times = 1000) resulted in a c-index of 0.72 and AUC values for predicting the 1-, 2-, and 3-year RFS of 0.77, 0.76, and 0.76, respectively. Through a Cox regression analysis, AMLlnc69, nucleophosmin mutation, and white blood cell at diagnosis were identified as independent predictors of RFS. Finally, a nomogram was build using these 2 parameters, showing a c-index of 0.80 and 0.71 after bootstrapping (n = 1000). In conclusion, the identified AMLlnc69 will, after prospective validation, add important information to guide the management of patients with pedAML. The nomogram is a promising tool for easy stratification of patients into a novel scheme of relapse-risk groups.

Epithelial-Mesenchymal Transition in Acute Leukemias

Epithelial-mesenchymal transition (EMT) is a metabolic process that confers phenotypic flexibility to cells and the ability to adapt to new functions. This transition is critical during embryogenesis and is required for the differentiation of many tissues and organs. EMT can also be induced in advanced-stage cancers, leading to further malignant behavior and chemotherapy resistance, resulting in an unfavorable prognosis for patients. Although EMT was long considered and studied only in solid tumors, it has been shown to be involved in the pathogenesis of hematological malignancies, including acute leukemias. Indeed, there is increasing evidence that EMT promotes the progression of acute leukemias, leading to the emergence of a more aggressive phenotype of the disease, and also causes chemotherapy resistance. The current literature suggests that the levels and activities of EMT inducers and markers can be used to predict prognosis, and that targeting EMT in addition to conventional therapies may increase treatment success in acute leukemias.

Combination Therapy and Dual-Target Inhibitors Based on LSD1: New Emerging Tools in Cancer Therapy

Lysine specific demethylase 1 (LSD1), a transcriptional modulator that represses or activates target gene expression, is overexpressed in many cancer and causes imbalance in the expression of normal gene networks. Over two decades, numerous LSD1 inhibitors have been reported, especially some of which have entered clinical trials, including eight irreversible inhibitors (TCP, ORY-1001, GSK-2879552, INCB059872, IMG-7289, ORY-2001, TAK-418, and LH-1802) and two reversible inhibitors (CC-90011 and SP-2577). Most clinical LSD1 inhibitors demonstrated enhanced efficacy in combination with other agents. LSD1 multitarget inhibitors have also been reported, exampled by clinical dual LSD1/histone deacetylases (HDACs) inhibitors 4SC-202 and JBI-802. Herein, we present a comprehensive overview of the combination of LSD1 inhibitors with various antitumor agents, as well as LSD1 multitarget inhibitors. Additionally, the challenges and future research directionsare also discussed, and we hope this review will provide new insight into the development of LSD1-targeted anticancer agents.

Twist family BHLH transcription factor 1 is required for the maintenance of leukemia stem cell in MLL-AF9+ acute myeloid leukemia

Leukemia stem cells (LSC) are a rare population capable of limitless self-renewal and are responsible for the initiation, maintenance, and relapse of leukemia. Elucidation of the mechanisms underlying the regulation of LSC function could provide novel treatment strategies. Here, we show that TWIST1 is extremely highly expressed in the LSC of MLL-AF9+ acute myeloid leukemia (AML), and its upregulation is positively regulated by KDM4C in a H3K9me3 demethylation-dependent manner. We further demonstrate that TWIST1 is essential for the viability, dormancy, and self-renewal capacities of LSC, and that it promotes the initiation and maintenance of MLL-AF9-mediated AML. In addition, TWIST1 directly interacts and collaborates with HOXA9 in inducing AML in mice. Mechanistically, TWIST1 exerts its oncogenic function by activating the WNT5a/RAC1 axis. Collectively, our study uncovers a critical role of TWIST1 in LSC function and provides new mechanistic insights into the pathogenesis of MLL-AF9+ AML.

ALOX5AP is a new prognostic indicator in acute myeloid leukemia

BACKGROUND

The overexpression of ALOX5AP has been observed in many types of cancer and has been identified as an oncogene. However, its role in acute myeloid leukemia (AML) has not been extensively studied. This study aimed to identify the expression and methylation patterns of ALOX5AP in bone marrow (BM) samples of AML patients, and further explore its clinical significance.

METHODS

Eighty-two de novo AML patients and 20 healthy donors were included in the study. Meanwhile, seven public datasets from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) were included to confirm the alteration of ALOX5AP. Receiver operating characteristic (ROC) curve analysis was applied to determine the discriminative capacity of ALOX5AP expression to discriminate AML. The prognostic value of ALOX5AP was identified by the Kaplan-Meier method and log-rank test. It was further validated in four independent cohorts (n = 1186). Significantly different genes associated with ALOX5AP expression were subsequently compared by LinkedOmics, and Metascape database.

RESULTS

The level of ALOX5AP expression was significantly increased in bone marrow cells of AML patients compared with healthy donors (P < 0.05). ROC curve analysis suggested that ALOX5AP expression might be a potential biomarker to discriminate AML from controls. ALOX5AP overexpression was associated with decreased overall survival (OS) in AML according to the TCGA data (P = 0.006), which was validated by other four independent cohorts. DNA methylation levels of ALOX5AP were significantly lower in AML patients compared to normal samples (P < 0.05), as confirmed in the Diseasemeth database and the independent cohort GSE63409. ALOX5AP level was positively associated with genes with proleukemic effects such as PAX2, HOX family, SOX11, H19, and microRNAs that act as oncogenes in leukemia, such as miR125b, miR-93, miR-494, miR-193b, while anti-leukemia-related genes and tumor suppressor microRNAs such as miR-582, miR-9 family and miR-205 were negatively correlated.

CONCLUSION

ALOX5AP overexpression, associated with its hypomethylation, predicts poorer prognosis in AML.

Bone marrow mesenchymal/fibroblastic stromal cells induce a distinctive EMT-like phenotype in AML cells

The development of epithelial-to-mesenchymal transition (EMT) like features is emerging as a critical factor involved in the pathogenesis of acute myeloid leukaemia (AML). However, the extracellular signals and the signalling pathways in AML that may regulate EMT remain largely unstudied. We found that the bone marrow (BM) mesenchymal/fibroblastic cell line HS5 induces an EMT-like migratory phenotype in AML cells. AML cells underwent a strong increase of vimentin (VIM) levels that was not mirrored to the same extent by changes of expression of the other EMT core proteins SNAI1 and SNAI2. We validated these particular pattern of co-expression of core-EMT markers in AML cells by performing an in silico analysis using datasets of human tumours. Our data showed that in AML the expression levels of VIM does not completely correlate with the co-expression of core EMT markers observed in epithelial tumours. We also found that vs epithelial tumours, AML cells display a distinct patterns of co-expression of VIM and the actin binding and adhesion regulatory proteins that regulate F-actin dynamics and integrin-mediated adhesions involved in the invasive migration in cells undergoing EMT. We conclude that the BM stroma induces an EMT related pattern of migration in AML cells in a process involving a distinctive regulation of EMT markers and of regulators of cell adhesion and actin dynamics that should be further investigated. Understanding the tumour specific signalling pathways associated with the EMT process may contribute to the development of new tailored therapies for AML as well as in different types of cancers.

An "unexpected" role for EMT transcription factors in hematological development and malignancy

The epithelial to mesenchymal transition (EMT) is a fundamental developmental process essential for normal embryonic development. It is also important during various pathogenic processes including fibrosis, wound healing and epithelial cancer cell metastasis and invasion. EMT is regulated by a variety of cell signalling pathways, cell-cell interactions and microenvironmental cues, however the key drivers of EMT are transcription factors of the ZEB, TWIST and SNAIL families. Recently, novel and unexpected roles for these EMT transcription factors (EMT-TFs) during normal blood cell development have emerged, which appear to be largely independent of classical EMT processes. Furthermore, EMT-TFs have also begun to be implicated in the development and pathogenesis of malignant hematological diseases such as leukemia and lymphoma, and now present themselves or the pathways they regulate as possible new therapeutic targets within these malignancies. In this review, we discuss the ZEB, TWIST and SNAIL families of EMT-TFs, focusing on what is known about their normal roles during hematopoiesis as well as the emerging and "unexpected" contribution they play during development and progression of blood cancers.

Stromal-induced epithelial-mesenchymal transition induces targetable drug resistance in acute lymphoblastic leukemia

The bone marrow microenvironment (BME) drives drug resistance in acute lymphoblastic leukemia (ALL) through leukemic cell interactions with bone marrow (BM) niches, but the underlying mechanisms remain unclear. Here, we show that the interaction between ALL and mesenchymal stem cells (MSCs) through integrin β1 induces an epithelial-mesenchymal transition (EMT)-like program in MSC-adherent ALL cells, resulting in drug resistance and enhanced survival. Moreover, single-cell RNA sequencing analysis of ALL-MSC co-culture identifies a hybrid cluster of MSC-adherent ALL cells expressing both B-ALL and MSC signature genes, orchestrated by a WNT/β-catenin-mediated EMT-like program. Blockade of interaction between β-catenin and CREB binding protein impairs the survival and drug resistance of MSC-adherent ALL cells in vitro and results in a reduction in leukemic burden in vivo. Targeting of this WNT/β-catenin-mediated EMT-like program is a potential therapeutic approach to overcome cell extrinsically acquired drug resistance in ALL.

High Expression of Microtubule-associated Protein TBCB Predicts Adverse Outcome and Immunosuppression in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a devastating blood cancer with high heterogeneity and ill-fated outcome. Despite numerous advances in AML treatment, the prognosis remains poor for a significant proportion of patients. Consequently, it is necessary to accurately and comprehensively identify biomarkers as soon as possible to enhance the efficacy of diagnosis, prognosis and treatment of AML. In this study, we aimed to identify prognostic markers of AML by analyzing the cohorts from TCGA-LAML database and GEO microarray datasets. Interestingly, the transcriptional level of microtubule-associated protein TBCB in AML patients was noticeably increased when compared with normal individuals, and this was verified in two independent cohorts (GSE9476 and GSE13159) and with our AML patients. Furthermore, univariate and multivariate regression analysis revealed that high TBCB expression was an independent poor prognostic factor for AML. GO and GSEA enrichment analysis hinted that immune-related signaling pathways were enriched in up-regulated DEGs between two populations separated by the median expression level of TBCB. By constructing a protein-protein interaction network, we obtained six hub genes, all of which are immune-related molecules, and their expression levels were positively linked to that of TBCB. In addition, the high expression of three hub genes was significantly associated with a poor prognosis in AML. Moreover, we found that the tumor microenvironment in AML with high TBCB expression tended to be infiltrated by NK cells, especially CD56^bright NK cells. The transcriptional levels of NK cell inhibitory receptors and their ligands were positively related to that of TBCB, and their high expression levels also predicted poor prognosis in AML. Notably, we found that the down-regulation of TBCB suppressed cell proliferation in AML cell lines by enhancing the apoptosis and cell cycle arrest. Finally, drug sensitivity prediction illustrated that cells with high TBCB expression were more responsive to ATRA and midostaurin but resistant to cytarabine, dasatinib, and imatinib. In conclusion, our findings shed light on the feasibility of TBCB as a potential predictor of poor outcome and to be an alternative target of treatment in AML.

An optimized approach for multiplexing single-nuclear ATAC-seq using oligonucleotide-conjugated antibodies

BACKGROUND

Single-cell technologies to analyze transcription and chromatin structure have been widely used in many research areas to reveal the functions and molecular properties of cells at single-cell resolution. Sample multiplexing techniques are valuable when performing single-cell analysis, reducing technical variation and permitting cost efficiencies. Several commercially available methods have been used in many scRNA-seq studies. On the other hand, while several methods have been published, multiplexing techniques for single nuclear assay for transposase-accessible chromatin (snATAC)-seq assays remain under development. We developed a simple nucleus hashing method using oligonucleotide-conjugated antibodies recognizing nuclear pore complex proteins, NuHash, to perform snATAC-seq library preparations by multiplexing.

RESULTS

We performed multiplexing snATAC-seq analyses on a mixture of human and mouse cell samples (two samples, 2-plex, and four samples, 4-plex) using NuHash. The analyses on nuclei with at least 10,000 read counts showed that the demultiplexing accuracy of NuHash was high, and only ten out of 9144 nuclei (2-plex) and 150 of 12,208 nuclei (4-plex) had discordant classifications between NuHash demultiplexing and discrimination using reference genome alignments. The differential open chromatin region (OCR) analysis between female and male samples revealed that male-specific OCRs were enriched in chromosome Y (four out of nine). We also found that five female-specific OCRs (20 OCRs) were on chromosome X. A comparative analysis between snATAC-seq and deeply sequenced bulk ATAC-seq on the same samples revealed that the bulk ATAC-seq signal intensity was positively correlated with the number of cell clusters detected in snATAC-seq. Moreover, when we categorized snATAC-seq peaks based on the number of cell clusters in which the peak was present, we observed different distributions over different genomic features between the groups. This result suggests that the peak intensities of bulk ATAC-seq can be used to identify different types of functional loci.

CONCLUSIONS

Our multiplexing method using oligo-conjugated anti-nuclear pore complex proteins, NuHash, permits high-accuracy demultiplexing of samples. The NuHash protocol is straightforward, works on frozen samples, and requires no modifications for snATAC-seq library preparation.

Altered physical phenotypes of leukemia cells that survive chemotherapy treatment

The recurrence of cancer following chemotherapy treatment is a major cause of death across solid and hematologic cancers. In B-cell acute lymphoblastic leukemia (B-ALL), relapse after initial chemotherapy treatment leads to poor patient outcomes. Here we test the hypothesis that chemotherapy-treated versus control B-ALL cells can be characterized based on cellular physical phenotypes. To quantify physical phenotypes of chemotherapy-treated leukemia cells, we use cells derived from B-ALL patients that are treated for 7 days with a standard multidrug chemotherapy regimen of vincristine, dexamethasone, and L-asparaginase (VDL). We conduct physical phenotyping of VDL-treated versus control cells by tracking the sequential deformations of single cells as they flow through a series of micron-scale constrictions in a microfluidic device; we call this method Quantitative Cyclical Deformability Cytometry. Using automated image analysis, we extract time-dependent features of deforming cells including cell size and transit time (TT) with single-cell resolution. Our findings show that VDL-treated B-ALL cells have faster TTs and transit velocity than control cells, indicating that VDL-treated cells are more deformable. We then test how effectively physical phenotypes can predict the presence of VDL-treated cells in mixed populations of VDL-treated and control cells using machine learning approaches. We find that TT measurements across a series of sequential constrictions can enhance the classification accuracy of VDL-treated cells in mixed populations using a variety of classifiers. Our findings suggest the predictive power of cell physical phenotyping as a complementary prognostic tool to detect the presence of cells that survive chemotherapy treatment. Ultimately such complementary physical phenotyping approaches could guide treatment strategies and therapeutic interventions. Insight box Cancer cells that survive chemotherapy treatment are major contributors to patient relapse, but the ability to predict recurrence remains a challenge. Here we investigate the physical properties of leukemia cells that survive treatment with chemotherapy drugs by deforming individual cells through a series of micron-scale constrictions in a microfluidic channel. Our findings reveal that leukemia cells that survive chemotherapy treatment are more deformable than control cells. We further show that machine learning algorithms applied to physical phenotyping data can predict the presence of cells that survive chemotherapy treatment in a mixed population. Such an integrated approach using physical phenotyping and machine learning could be valuable to guide patient treatments.

A novel scoring system for the quantitative prediction of prognosis in acute myeloid leukemia

Background

Acute myeloid leukemia (AML) is a heterogeneous hematopoietic malignancy. Patient prognosis cannot be accurately assessed in National Comprehensive Cancer Network (NCCN) risk stratification subgroups based on the current criteria. This study aimed to develop a novel prognostic score model for the quantitative prediction of prognosis in AML.

Results

We developed a prognostic risk scoring model of AML using differentially expressed genes to predict prognosis in patients with AML. Furthermore, we evaluated the effectiveness and clinical significance of this prognostic model in 4 AML cohorts and 905 patients with AML. A prognostic risk scoring model of AML containing eight prognosis-related genes was constructed using a multivariate Cox regression model. The model had a higher predictive value for the prognosis of AML in the training and validation sets. In addition, patients with lower scores had significantly better overall survival (OS) and even-free survival (EFS) than those with higher scores among patients with intermediate-risk AML according to the NCCN guidelines, indicating that the model could be used to further predict the prognosis of the intermediate-risk AML populations. Similarly, patients with high scores had remarkably poor OS and EFS in the normal-karyotype populations, indicating that the scoring model had an excellent predictive performance for patients with AML having normal karyotype.

Conclusions

Our study provided an individualized prognostic risk score model that could predict the prognosis of patients with AML.

A genome-wide CRISPR activation screen identifies SCREEM a novel SNAI1 super-enhancer demarcated by eRNAs

The genome is pervasively transcribed to produce a vast array of non-coding RNAs (ncRNAs). Long non-coding RNAs (lncRNAs) are transcripts of >200 nucleotides and are best known for their ability to regulate gene expression. Enhancer RNAs (eRNAs) are subclass of lncRNAs that are synthesized from enhancer regions and have also been shown to coordinate gene expression. The biological function and significance of most lncRNAs and eRNAs remain to be determined. Epithelial to mesenchymal transition (EMT) is a ubiquitous cellular process that occurs during cellular migration, homeostasis, fibrosis, and cancer-cell metastasis. EMT-transcription factors, such as SNAI1 induce a complex transcriptional program that coordinates the morphological and molecular changes associated with EMT. Such complex transcriptional programs are often subject to coordination by networks of ncRNAs and thus can be leveraged to identify novel functional ncRNA loci. Here, using a genome-wide CRISPR activation (CRISPRa) screen targeting ∼10,000 lncRNA loci we identified ncRNA loci that could either promote or attenuate EMT. We discovered a novel locus that we named SCREEM (SNAI1 cis-regulatory eRNAs expressed in monocytes). The SCREEM locus contained a cluster of eRNAs that when activated using CRISPRa induced expression of the neighboring gene SNAI1, driving concomitant EMT. However, the SCREEM eRNA transcripts themselves appeared dispensable for the induction of SNAI1 expression. Interestingly, the SCREEM eRNAs and SNAI1 were co-expressed in activated monocytes, where the SCREEM locus demarcated a monocyte-specific super-enhancer. These findings suggest a potential role for SNAI1 in monocytes. Exploration of the SCREEM-SNAI axis could reveal novel aspects of monocyte biology.

Integrative analysis reveals histone demethylase LSD1 promotes RNA polymerase II pausing

Lysine-specific demethylase 1 (LSD1) is well-known for its role in decommissioning enhancers during mouse embryonic stem cell (ESC) differentiation. Its role in gene promoters remains poorly understood despite its widespread presence at these sites. Here, we report that LSD1 promotes RNA polymerase II (RNAPII) pausing, a rate-limiting step in transcription regulation, in ESCs. We found the knockdown of LSD1 preferentially affects genes with higher RNAPII pausing. Next, we demonstrate that the co-localization sites of LSD1 and MYC, a factor known to regulate pause-release, are enriched for other RNAPII pausing factors. We show that LSD1 and MYC directly interact and MYC recruitment to genes co-regulated with LSD1 is dependent on LSD1 but not vice versa. The co-regulated gene set is significantly enriched for housekeeping processes and depleted of transcription factors compared to those bound by LSD1 alone. Collectively, our integrative analysis reveals a pleiotropic role of LSD1 in promoting RNAPII pausing.

•

LSD1 promotes RNA polymerase II pausing in mouse embryonic stem cells

•

LSD1 knockdown causes global reduction of RNAPII pausing

•

Co-localized sites of LSD1 and MYC are enriched for RNAPII pausing and releasing factors

•

MYC recruitment to co-regulated genes is dependent on LSD1 but not vice versa

HMG20B stabilizes association of LSD1 with GFI1 on chromatin to confer transcription repression and leukemia cell differentiation block

Pharmacologic inhibition of LSD1 induces molecular and morphologic differentiation of blast cells in acute myeloid leukemia (AML) patients harboring MLL gene translocations. In addition to its demethylase activity, LSD1 has a critical scaffolding function at genomic sites occupied by the SNAG domain transcription repressor GFI1. Importantly, inhibitors block both enzymatic and scaffolding activities, in the latter case by disrupting the protein:protein interaction of GFI1 with LSD1. To explore the wider consequences of LSD1 inhibition on the LSD1 protein complex we applied mass spectrometry technologies. We discovered that the interaction of the HMG-box protein HMG20B with LSD1 was also disrupted by LSD1 inhibition. Downstream investigations revealed that HMG20B is co-located on chromatin with GFI1 and LSD1 genome-wide; the strongest HMG20B binding co-locates with the strongest GFI1 and LSD1 binding. Functional assays demonstrated that HMG20B depletion induces leukemia cell differentiation and further revealed that HMG20B is required for the transcription repressor activity of GFI1 through stabilizing LSD1 on chromatin at GFI1 binding sites. Interaction of HMG20B with LSD1 is through its coiled-coil domain. Thus, HMG20B is a critical component of the GFI1:LSD1 transcription repressor complex which contributes to leukemia cell differentiation block.

Low expression of the metabolism-related gene SLC25A21 predicts unfavourable prognosis in patients with acute myeloid leukaemia

Acute myeloid leukaemia (AML) is a heterogeneous disease associated with poor outcomes. To identify AML-specific genes with prognostic value, we analysed transcriptome and clinical information from The Cancer Genome Atlas (TCGA) database, Gene Expression Omnibus (GEO) datasets, and Genotype-Tissue Expression (GTEx) project. The metabolism-related gene, SLC25A21 was found to be significantly downregulated in AML, and was associated with high white blood cell (WBC) counts, high pretrial blood (PB) and bone marrow (BM) blast abundance, FLT3 mutation, NPM1 mutation, and death events (all p value <0.05). We validated the expression of SLC25A21 in our clinical cohort, and found that SLC25A21 was downregulated in AML. Moreover, we identified low expression of SLC25A21 as an independent prognostic factor by univariate Cox regression (hazard ratio [HR]: 0.550; 95% Confidence interval [CI]: 0.358–0.845; p value = 0.006) and multivariate Cox regression analysis (HR: 0.341; 95% CI: 0.209–0.557; p value <0.05). A survival prediction nomogram was established with a C-index of 0.735, which indicated reliable prognostic prediction. Subsequently, based on the median SLC25A21 expression level, patients in the TCGA-LAML cohort were divided into low- and high-expression groups. Gene ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs highlighted growth factor binding, extracellular structure organization, cytokine‒cytokine receptor interaction, etc. The results of gene set enrichment analysis (GSEA) indicated that the epithelial-mesenchymal transition, KRAS signalling, oxidative phosphorylation, and reactive oxygen species pathways were enriched. Through gene coexpression and protein‒protein interaction (PPI) network analysis, we identified two hub genes, EGFR and COL1A2, which were linked to worse clinical outcomes. Furthermore, we found that lower SLC25A21 expression was closely associated with a significant reduction in the levels of infiltrating immune cells, which might be associated with immune escape of AML cells. A similar trend was observed for the expression of checkpoint genes (CTLA4, LAG3, TIGIT, and HAVCR2). Finally, drug sensitivity testing suggested that the low-expression SLC25A21 group is sensitive to doxorubicin, mitomycin C, linifanib but resistant to JQ1, belinostat, and dasatinib. Hence, our study demonstrated that a low expression level of SLC25A21 predicts an unfavourable prognosis in patients with AML.

SMGR: a joint statistical method for integrative analysis of single-cell multi-omics data

Unravelling the regulatory programs from single-cell multi-omics data has long been one of the major challenges in genomics, especially in the current emerging single-cell field. Currently there is a huge gap between fast-growing single-cell multi-omics data and effective methods for the integrative analysis of these inherent sparse and heterogeneous data. In this study, we have developed a novel method, Single-cell Multi-omics Gene co-Regulatory algorithm (SMGR), to detect coherent functional regulatory signals and target genes from the joint single-cell RNA-sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data obtained from different samples. Given that scRNA-seq and scATAC-seq data can be captured by zero-inflated Negative Binomial distribution, we utilize a generalized linear regression model to identify the latent representation of consistently expressed genes and peaks, thus enables the identification of co-regulatory programs and the elucidation of regulating mechanisms. Results from both simulation and experimental data demonstrate that SMGR outperforms the existing methods with considerably improved accuracy. To illustrate the biological insights of SMGR, we apply SMGR to mixed-phenotype acute leukemia (MPAL) and identify the MPAL-specific regulatory program with significant peak-gene links, which greatly enhance our understanding of the regulatory mechanisms and potential targets of this complex tumor.

Comprehensive analysis of ferroptosis-related gene signatures as a potential therapeutic target for acute myeloid leukemia: A bioinformatics analysis and experimental verification

Background

Ferroptosis plays an important role in the development of acute myeloid leukemia (AML); however, the exact role of ferroptosis-related genes in the prognosis of AML patients is unclear.

Methods

RNA sequencing data and the clinicopathological characteristics of AML patients were obtained from The Cancer Genome Atlas database, and ferroptosis-related genes were obtained from the FerrDb database. Cox regression analysis and least absolute shrinkage and selection operator analysis were performed to identify ferroptosis-related gene signatures. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and single-sample gene set enrichment analysis (ssGSEA) were performed to explore the biological functions of the ferroptosis-related genes. Finally, ferroptosis of AML cells was induced by erastin and sulfasalazine to detect the changes in the expression of relevant prognostic genes and explore the underlying mechanisms using quantitative real-time polymerase chain reaction (qRT-PCR).

Results

Seven ferroptosis-related gene signatures (SOCS1, ACSF2, MYB, EIF2AK4, AIFM2, SLC7A11, and GPX4) were identified in the training group. Kaplan-Meier and Cox regression analyses confirmed that risk score was an independent prognostic predictor of AML in the training and validation groups (P<0.05). Further, functional enrichment analysis revealed that seven ferroptosis-related genes were associated with many immune-related biological processes. Most importantly, erastin and sulfasalazine can induce the ferroptosis of AML cells. Overall, SLC7A11 and the SLC7A11/xCT-GSH-GPX4 pathway may be the respective key gene and potential regulatory pathway in erastin- and sulfasalazine-induced ferroptosis of AML cells.

Conclusions

A novel signature involving seven ferroptosis-related genes that could accurately predict AML prognosis was identified. Further, the Food and Drug Administration-approved drug, sulfasalazine, was demonstrated for the first time to induce the ferroptosis of AML cells. SLC7A11 and the SLC7A11/xCT-GSH-GPX4 pathway may be the respective key gene and underlying mechanism in this process, ultimately providing new insights into the strategies for the development of new AML therapies.

Mechanisms of carcinogenic activity triggered by lysine-specific demethylase 1A

Epigenetics has emerged as a prime focus area in the field of cancer research. Lysine-specific demethylase 1A (LSD1), the first discovered histone demethylase, is mainly responsible for catalysing demethylation of histone 3 lysine 4 (H3K4) and H3K9 to activate or inhibit gene transcription. LSD1 is abnormally expressed in various cancers and participates in cancer proliferation, apoptosis, metastasis, invasion, drug resistance and other processes by interacting with regulatory factors. Therefore, it may serve as a potential therapeutic target for cancer. This review summarises the major oncogenic mechanisms mediated by LSD1 and provides a reference for developing novel and efficient anticancer strategies targeting LSD1.

A comprehensive comparative study on LSD1 in different cancers and tumor specific LSD1 inhibitors

LSD1 was significantly over-expressed in several cancer types, and its aberrant overexpression was revealed to play a crucial role in the initiation and progression of cancer. Several LSD1 inhibitors that were discovered and developed so far were found to be effective in attenuating tumor growth in both in vivo and in vitro studies. However, the major challenge associated with the development of cancer therapies is personalized treatment. Therefore, it is essential to look in detail at how LSD1 plays its part in carcinogenesis and whether there are any different expression levels of LSD1 in different tumors. Here in this review, fresh insight into a list of function correlated LSD1 binding proteins are provided, and we tried to figure out the role of LSD1 in different cancer types, including hematological malignancies and solid tumors. A critical description of mutation preference for LSD1 in different tumors was also discussed. Recent research findings clearly showed that the abrogation of LSD1 demethylase activity via LSD1 inhibitors markedly reduced the growth of cancer cells. But there are still many ambiguities regarding the role of LSD1 in different cancers. Therefore, targeting LSD1 for treating different cancers is still reductionist, and many challenges need to be met to improve the therapeutic outcomes of LSD1 inhibitors.

CDK7/12/13 inhibition targets an oscillating leukemia stem cell network and synergizes with venetoclax in acute myeloid leukemia

The heterogeneous response of acute myeloid leukemia (AML) to current anti‐leukemic therapies is only partially explained by mutational heterogeneity. We previously identified GPR56 as a surface marker associated with poor outcome across genetic groups, which characterizes two leukemia stem cell (LSC)‐enriched compartments with different self‐renewal capacities. How these compartments self‐renew remained unclear. Here, we show that GPR56⁺ LSC compartments are promoted in a complex network involving epithelial‐to‐mesenchymal transition (EMT) regulators besides Rho, Wnt, and Hedgehog (Hh) signaling. Unexpectedly, Wnt pathway inhibition increased the more immature, slowly cycling GPR56⁺CD34⁺ fraction and Hh/EMT gene expression, while Wnt activation caused opposite effects. Our data suggest that the crucial role of GPR56 lies in its ability to co‐activate these opposing signals, thus ensuring the constant supply of both LSC subsets. We show that CDK7 inhibitors suppress both LSC‐enriched subsets in vivo and synergize with the Bcl‐2 inhibitor venetoclax. Our data establish reciprocal transition between LSC compartments as a novel concept underlying the poor outcome in GPR56^high AML and propose combined CDK7 and Bcl‐2 inhibition as LSC‐directed therapy in this disease.

Targeting Thyrointegrin αvβ3 Using Fluorobenzyl Polyethylene Glycol Conjugated Tetraiodothyroacetic Acid (NP751) in Acute Myeloid Leukemia

Background

Acute myeloid leukemia (AML) is associated with poor long-term survival, even with newer therapeutic agents. Here, we show the results of our preclinical study, in which we evaluated the efficacy of a new thyrointegrin αvβ3 antagonist, named fluorobenzyl polyethylene glycol conjugated tetraiodothyroacetic acid (fb-PMT).

Methods and Results

fb-PMT (NP751) is a potent αvβ3 antagonist of molecular weight of 2,478.9 Da. it represents a conjugate of tetraiodothyroacetic acid (tetrac) and monodisperse polyethylene glycol (PEG36), with a 4-fluorobenzyl group capping the other end of the PEG. fb-PMT effectively suppresses the malignant growth of human acute myeloid leukemia (AML) after successful engraftment in transgenic NSG-S xenograft mouse models of either established human AML cell line or primary AML cells. Daily treatment with fb-PMT (1–10 mg/kg body weight) subcutaneously (s.c.) for 3–4 weeks was associated with marked regression of leukemogenesis and extended survival in both models. The efficiency of the fb-PMT therapy was verified using in vivo imaging system (IVIS) imaging, flow cytometry, and histopathological examination to monitor the engraftment of leukemic cells in the bone marrow and other organs. fb-PMT therapy for 3–4 weeks at 3 and 10 mg/kg daily doses exhibited significant reduction (p < 0.0001) of leukemic cell burden of 74% and >95%, respectively. All fb-PMT-treated mice in the 10 mg/kg treatment arm successfully maintained remission after discontinuing the daily treatment. Comprehensive fb-PMT safety assessments demonstrated excellent safety and tolerability at multiple folds above the anticipated human therapeutic doses. Lastly, our genome-wide microarray screens demonstrated that fb-PMT works through the molecular interference mechanism with multiple signaling pathways contributing to growth and survival of leukemic cells.

Conclusion

Our preclinical findings of the potent anticancer activities of fb-PMT and its favorable safety profiles warrant its clinical investigation for the effective and safe management of AML.

Targeting lysine-specific demethylase 1A inhibits renal epithelial-mesenchymal transition and attenuates renal fibrosis

Lysine-specific histone demethylase 1 (LSD1) as the first identified histone/lysine demethylase regulates gene expression and protein functions in diverse diseases. In this study, we show that the expression of LSD1 is increased in mouse kidneys with unilateral ureteral obstruction (UUO) and in cultured NRK-52E cells undergoing TGF-β1-induced epithelial-mesenchymal transition (EMT). Inhibition of LSD1 with its specific inhibitor ORY1001 attenuated renal EMT and fibrosis, which was associated with decreased the deposition of extracellular matrix proteins and the expression of fibrotic markers, including α-smooth muscle actin (α-SMA) and fibronectin, and the recovery of E-cadherin expression and decrease of N-cadherin expression in UUO kidneys and in NRK-52E cells induced with TGF-β1. Targeting LSD1 also decreased the expression of Snail family transcriptional repressor 1 (Snail-1) and its interaction with LSD1 in UUO kidneys and in NRK-52E cells treated with TGF-β1. In addition, we identified a novel LSD1-14-3-3ζ-PKCα axis in the regulation of the activation of AKT and Stat3 and then the activation of fibroblasts. This study suggests that LSD1 plays a critical role in regulation of renal EMT and fibrosis through activation of diverse signaling pathways and places an emphasis that LSD1 has potential as a therapeutic target for the treatment of renal fibrosis.

Aspects of the Epigenetic Regulation of EMT Related to Cancer Metastasis

Epithelial to mesenchymal transition (EMT) occurs during the pathological process associated with tumor progression and is considered to influence and promote the metastatic cascade. Characterized by loss of cell adhesion and apex base polarity, EMT enhances cell motility and metastasis. The key markers of the epithelial to mesenchymal transition are proteins characteristic of the epithelial phenotype, e.g., E-cadherin, cytokeratins, occludin, or desmoplakin, the concentration and activity of which are reduced during this process. On the other hand, as a result of acquiring the characteristics of mesenchymal cells, an increased amount of N-cadherin, vimentin, fibronectin, or vitronectin is observed. Importantly, epithelial cells undergo partial EMT where some of the cells show both epithelial and mesenchymal characteristics. The significant influence of epigenetic regulatory mechanisms is observed in the gene expression involved in EMT. Among the epigenetic modifications accompanying incorrect genetic reprogramming in cancer are changes in the level of DNA methylation within the CpG islands and posttranslational covalent changes of histone proteins. All observed modifications, which are stable but reversible changes, affect the level of gene expression leading to the development and progression of the disease, and consequently affect the uncontrolled growth of the population of cancer cells.

Establishment of a prognostic ferroptosis-related gene profile in acute myeloid leukaemia

Acute myeloid leukaemia (AML) is a heterogeneous disease with a difficult to predict prognosis. Ferroptosis, an iron‐induced programmed cell death, is a promising target for cancer therapy. Nevertheless, not much is known about the relationship between ferroptosis‐related genes and AML prognosis. Herein, we retrieved RNA profile and corresponding clinical data of AML patients from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Univariate Cox analysis was employed to identify ferroptosis‐related genes significantly associated with AML prognosis. Next, the least absolute shrinkage and selection operator (LASSO) regression was employed to establish a prognostic ferroptosis‐related gene profile. 12 ferroptosis‐related genes were screened to generate a prognostic model, which stratified patients into a low‐ (LR) or high‐risk (HR) group. Using Kaplan‐Meier analysis, we demonstrated that the LR patients exhibited better prognosis than HR patients. Moreover, receiver operating characteristic (ROC) curve analysis confirmed that the prognostic model showed good predictability. Functional enrichment analysis indicated that the infiltration of regulatory T cells (Treg) differed vastly between the LR and HR groups. Our prognostic model can offer guidance into the accurate prediction of AML prognosis and selection of personalized therapy in clinical practice.

RNA m6A Demethylase ALKBH5 Protects Against Pancreatic Ductal Adenocarcinoma via Targeting Regulators of Iron Metabolism

Although RNA m⁶A regulators have been implicated in the tumorigenesis of several different types of tumors, including pancreatic cancer, their clinical relevance and intrinsic regulatory mechanism remain elusive. This study analyzed eight m⁶A regulators (METTL3, METTL14, WTAP, FTO, ALKBH5, and YTHDF1-3) in pancreatic ductal adenocarcinoma (PDAC) and found that only RNA m⁶A demethylase ALKBH5 serves as an independent favorable prognostic marker for this tumor. To better understand the molecular mechanism underlying the protective effect conferred by ALKBH5 against pancreatic tumorigenesis, we performed a transcriptome-wide analysis of m⁶A methylation, gene expression, and alternative splicing (AS) using the MIA PaCa-2 stable cell line with ALKBH5 overexpression. We demonstrated that ALKBH5 overexpression induced a reduction in RNA m⁶A levels globally. Furthermore, mRNAs encoding ubiquitin ligase FBXL5, and mitochondrial iron importers SLC25A28 and SLC25A37, were identified as substrates of ALKBH5. Mechanistically, the RNA stabilities of FBXL5 and SLC25A28, and the AS of SLC25A37 were affected, which led to their upregulation in pancreatic cancer cell line. Particularly, we observed that downregulation of FBXL5 in tumor samples correlated with shorter survival time of patients. Owing to FBXL5-mediated degradation, ALKBH5 overexpression incurred a significant reduction in iron-regulatory protein IRP2 and the modulator of epithelial-mesenchymal transition (EMT) SNAI1. Notably, ALKBH5 overexpression led to a significant reduction in intracellular iron levels as well as cell migratory and invasive abilities, which could be rescued by knocking down FBXL5. Overall, our results reveal a previously uncharacterized mechanism of ALKBH5 in protecting against PDAC through modulating regulators of iron metabolism and underscore the multifaceted role of m⁶A in pancreatic cancer.

Interplay between the EMT transcription factors ZEB1 and ZEB2 regulates hematopoietic stem and progenitor cell differentiation and hematopoietic lineage fidelity

The ZEB2 transcription factor has been demonstrated to play important roles in hematopoiesis and leukemic transformation. ZEB1 is a close family member of ZEB2 but has remained more enigmatic concerning its roles in hematopoiesis. Here, we show using conditional loss-of-function approaches and bone marrow (BM) reconstitution experiments that ZEB1 plays a cell-autonomous role in hematopoietic lineage differentiation, particularly as a positive regulator of monocyte development in addition to its previously reported important role in T-cell differentiation. Analysis of existing single-cell (sc) RNA sequencing (RNA-seq) data of early hematopoiesis has revealed distinctive expression differences between Zeb1 and Zeb2 in hematopoietic stem and progenitor cell (HSPC) differentiation, with Zeb2 being more highly and broadly expressed than Zeb1 except at a key transition point (short-term HSC [ST-HSC]➔MPP1), whereby Zeb1 appears to be the dominantly expressed family member. Inducible genetic inactivation of both Zeb1 and Zeb2 using a tamoxifen-inducible Cre-mediated approach leads to acute BM failure at this transition point with increased long-term and short-term hematopoietic stem cell numbers and an accompanying decrease in all hematopoietic lineage differentiation. Bioinformatics analysis of RNA-seq data has revealed that ZEB2 acts predominantly as a transcriptional repressor involved in restraining mature hematopoietic lineage gene expression programs from being expressed too early in HSPCs. ZEB1 appears to fine-tune this repressive role during hematopoiesis to ensure hematopoietic lineage fidelity. Analysis of Rosa26 locus–based transgenic models has revealed that Zeb1 as well as Zeb2 cDNA-based overexpression within the hematopoietic system can drive extramedullary hematopoiesis/splenomegaly and enhance monocyte development. Finally, inactivation of Zeb2 alone or Zeb1/2 together was found to enhance survival in secondary MLL-AF9 acute myeloid leukemia (AML) models attesting to the oncogenic role of ZEB1/2 in AML.

This study shows that the closely related transcription factors ZEB1 and ZEB2 cooperate to restrain myeloid and lymphoid differentiation programs in hematopoietic stem and progenitor cells, ensuring fidelity of differentiation in multiple lineages.

The novel LSD1 inhibitor ZY0511 suppresses diffuse large B-cell lymphoma proliferation by inducing apoptosis and autophagy

Lysine-specific demethylase 1 (LSD1, also known as KDM1A) is an attractive agent for treatment of cancer. However, the anti-tumor effect of LSD1 inhibitors against diffuse large B-cell lymphoma (DLBCL) and the underlying mechanism are still unclear. Here, we report that KDM1A is overexpressed in human DLBCL tissues and negatively related to overall survival rate of DLBCL patients. ZY0511, a novel and potent LSD1 inhibitor developed by our group, inhibited the proliferation of human DLBCL cells. ZY0511 interacted with LSD1, induced methylation level of histone 3 lysine 4 and histone 3 lysine 9 in DLBCL cells. Mechanistically, transcriptome sequencing results indicated that ZY0511 induced the genes enrichment significantly related to cell cycle, autophagy, and apoptosis signaling pathways. Further study confirmed that ZY0511 blocked cell cycle at G0/G1 phase and expression of CDK4 and cyclin D1. ZY0511 decreased mitochondrial membrane potential and induced apoptosis, which can be reverted by a pan-caspase inhibitor, Z-VAD-FMK. Moreover, ZY0511 treatment significantly increased autophagy-associated marker proteins and autophagosomes formation in DLBCL cells. In vivo xenograft experiments confirmed that intraperitoneal administration of ZY0511 significantly suppressed SU-DHL-6 xenograft tumor growth in vivo. In conclusion, our findings identify that ZY0511 inhibits DLBCL growth both in vitro and in vivo via the induction of apoptosis and autophagy, and LSD1 inhibitor might be a promising strategy for treating DLBCL.

miR-1287-5p upregulation inhibits the EMT and pro-inflammatory cytokines in LPS-induced human nasal epithelial cells (HNECs)

BACKGROUND

Chronic rhinosinusitis is an intractable symptom that influences daily lives of patients. miR-1287-5p was discovered to play a suppressive role in cervical cancer and HBV-related infection.

PURPOSE

This study investigated the potential role of miR-1287-5p in the in-vitro model of chronic rhinosinusitis.

METHODS

GSE169376 dataset was analyzed and differential miRNAs in nasal mucosa tissues in the chronic rhinosinusitis group were screened out. LPS was used to treat HNECs for 12h, 24h and 48h. Cells underwent LPS treatment after SNAI1 downregulation, miR-1287-5p upregulation or pretreatment of the HMGB1 inhibitor, Glycyrrhizin. RT-PCR was used to measure the RNA expression of miR-1287-5p, SNAI1 and HMGB1. ELISA was used for the detection of IL-6, IL-8, TNF-α changes. Targetscan and starBase were used to predict the targets (SNAI1 and HMGB1) of miR-1287-5p. Dual-luciferase reporter assays were applied to validate this. Western blot was used to analyze the protein changes of Snai1, Vimentin, E-cadherin and HMGB1.

RESULTS

miR-1287-5p was downregulated in the chronic rhinosinusitis group and decreased after LPS treatment in HNECs. The upregulation of miR-1287-5p inhibited IL-6, IL-8, TNF-α and EMT. miR-1287-5p targeted and inhibited SNAI1 and HMGB1. SNAI1 downregulation led to inhibition in EMT while loss of HMGB1 contributed to the decrease in pro-inflammatory cytokines. Knockdown of SNAI1 decreased HMGB1, resulting in the reduction of pro-inflammatory cytokines while HMGB1 inhibitor reduced SNAI1 and thus suppressed the EMT process.

CONCLUSION

miR-1287-5p downregulation was associated with chronic rhinosinusitis and its upregulation inhibited the EMT and inflammation in LPS-induced HNECs through Snai1/HMGB1 pathway.

Roles of lysine-specific demethylase 1 (LSD1) in homeostasis and diseases

Lysine-specific demethylase 1 (LSD1) targets mono- or di-methylated histone H3K4 and H3K9 as well as non-histone substrates and functions in the regulation of gene expression as a transcriptional repressor or activator. This enzyme plays a pivotal role in various physiological processes, including development, differentiation, inflammation, thermogenesis, neuronal and cerebral physiology, and the maintenance of stemness in stem cells. LSD1 also participates in pathological processes, including cancer as the most representative disease. It promotes oncogenesis by facilitating the survival of cancer cells and by generating a pro-cancer microenvironment. In this review, we discuss the role of LSD1 in several aspects of cancer, such as hypoxia, epithelial-to-mesenchymal transition, stemness versus differentiation of cancer stem cells, as well as anti-tumor immunity. Additionally, the current understanding of the involvement of LSD1 in various other pathological processes is discussed.