MAPK Signaling Pathway was Dysregulated in Acute Myeloid Leukemia with RUNX1 Mutations

doi:10.21203/rs.3.rs-502740/v1

Download PDF

Research Article

MAPK Signaling Pathway was Dysregulated in Acute Myeloid Leukemia with RUNX1 Mutations

https://doi.org/10.21203/rs.3.rs-502740/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

The purpose of this study was to investigate the relationship between RUNX1 mutations and MAPK signaling pathway in acute myeloid leukemia (AML). In this study, we analyzed miRNA expression with 5 mutant RUNX1 and 9 wild-type RUNX1 cases from TCGA database of AML. Six miRNAs were differently expressed with significance, and three of them were related to overall survival. Predicted target genes of these 3 miRNAs were highly enriched in MAPK signaling pathway by functional enrichment with miRWalk3.0. Besides, genes among RUNX1 associated genes directly regulated by RUNX1 were involved in MAPK signaling pathway, too. Taken together, we demonstrate three DEmiRNAs and three genes correlated to RUNX1 were correlated with prognosis in AML, and RUNX1 modulated MAPK signaling pathway in AML.

Oncology

AML

RUNX1

miRNA

mRNA

MAPK

Acute myeloid leukemia (AML) is a clonal proliferation disease of immature myeloid precursors, caused by series of genetic alterations in hematopoietic stem cells/progenitors(HSPCs)(1). Despite remarkable progress has been made in AML treatment, 5-year survival rates of ∼10%(2). To improve treatment efficiency, it is urgently required to understand the reveal mechanism of AML. Transcription factors (TF) controlling gene expression are vital for hematopoietic stem cells (HSCs) maintenance, differentiation, apoptosis, and aging, such as RUNX1, GATA2, ETV6, and CEBPA(3-5). Bunches of evidence showed changes in these TFs were critical for leukemogenesis(6). Illustrating the dependency of these altered TFs in AML would provide therapeutic targets and improve treatment efficiency of patients.

The Runt-related transcription factor 1 (RUNX1/AML1) gene located on chromosome 21q22 is one of the most frequently altered genes in AML, including chromosomal translocation and loss of function mutations. In FAB classification, M0 is the subtype with the most frequently mutant RUNX1(65.2%), followed by M1(32.4%) and M2(30.2%) (7). According to WHO classification of myeloid neoplasms and acute leukemia, AML with RUNX1 mutations is recognized as a new subtype (or category) of AML. Usually, RUNX1 forming heterodimer with CBF regulates target gene expression by activation or suppression(8, 9). RUNX1 is essential for hematopoietic stem cell generation, and largely dispensable in adult hematopoiesis. But heterozygous germline mutations in RUNX1 cause familial platelet disorder with a predisposition to AML(10). Different from chromosome translocation, loss of function (LOF) mutations in RUNX1 conferred stress-resistance to pre-leukemic stem cells(11). AML patients with RUNX1 mutations usually have poor treatment outcomes, including lower complete remission (CR) rate, inferior disease-free survival (DFS), and overall survival (OS) (12, 13). With the advance in high-through sequencing technology, cytogenetic and molecular genetic abnormalities were identified and co-occurred with RUNX1 mutations in AML, such as TP53, PHF6, DNMT3A, IDH1, TET2, ASXL1, FLT3, NPM1, WT1, RAS, and MLL(14, 15). But the molecular mechanism of chemo-resistance caused by RUNX1 mutations is still elusive. As a transcription factor, RUNX1 can activate or repress its target gene expression, including microRNAs and mRNAs.

MicroRNAs (miRNAs), are small short non-coding RNAs, which can regulate post-transcriptional gene expression through binding to 3’-untranslated regions (3’-UTRs) of target mRNAs(16, 17). As miRNAs are critical regulators in cell proliferation, differentiation, migration, and apoptosis, they have been confirmed that they were involved in tissue development, carcinogenesis, and tumor metastasis. In normal hematopoiesis, the formation of hematopoietic stem cells (HSCs), maintenance of HSCs, and differentiation of HSCs to lineages are regulated by miRNAs. In leukemia, dysregulation of miRNAs also has been demonstrated in pathogenesis. miRNA expression profiles were correlated to prognosis, which can be used for diagnosis, treatment, and prognosis prediction (18-20).

As a master hematopoietic transcription factor, RUNX1 can control miRNA expression, as well as being targeted by miRNAs (21). In megakaryopoiesis, RUNX1 determined the expression of the miRNA-144/451 (22), and its expression was controlled by miR-142 in HSCs(23). Further, many target genes of RUNX1 were screened out, which were essential for hematopoiesis. And lots of evidence of these genes participating in leukemogenesis were accumulated. Hmga2, one of RUNX1 target genes, regulated expansion of myeloid progenitors (24). CSF1R and MPO controlled by RUNX1 were involved in myelogenesis (25, 26). But the relationship of RUNX1 and their targets in AML was not fully illustrated. In this study, we found that RUNX1 specifically modulated MAPK signaling pathway by controlling miRNAs and mRNAs expression.

Data preparation

Use R package ‘TCGAbiolinks-GDCquery’, LAML-somatic mutation, miRNA expression, and gene expression data from GDC TCGA Acute Myeloid Leukemia (LAML) database were downloaded. There were 200 clinically annotated adult cases of de novo AML, including 19 patients with RUNX1 mutations. Patients classified in M0/M1 according to the FAB classification, ‘DECEASED ’statue, the mutations occurred in ‘Runt homology domain’ (amino acid change position<300) (27), no fusion-gene mutations (AML1-ETO or TEL-AML1), and available mRNA and miRNA expression data were used as criteria. There were only 5 cases with RUNX1 mutations according to these filtering criteria. ‘LIVING’ statue, normal karyotype, no fusion-gene mutations, and availability of mRNA and miRNA data were found in 9 cases with wild-type RUNX1. The clinical characteristics of AML samples were summarized in Supplemental 1, sample characteristics diagnosis after filter in AML. The GSE106291 dataset from GEO database was used as the external validation cohort.

Statistical methods

R package ‘DESeq2’ (28) was used to perform the differential analysis. The differently expressed miRNAs with absolute fold changes over 0.5 and p-value less than 0.05 were identified with significance. Spearman rank correlation was used to analyze correlations between RUNX1 with miRNAs or mRNAs, and a p-value less than 0.05 was significant. Cox regression model and log-rank test used for survival analysis by the R packages ‘survival’ and ‘ggplot2’.

Prediction of miRNA target genes and identification of RUNX1 target genes

MiRBase database was used to check miRNA two different matures, (‘-5p’ and ‘-3p’), and confirm previous IDs of differently expressed miRNAs. miRNA target gene was predicted by miRWalk3.0(29-31) which integrates resources of Targetscan, Mirdb, and experimental verification of Mirtarbase. Target region in 3’-UTR position and score more than 0.95 of miRNAs were used as selection criteria in at least one of these three databases.

Functional annotation and pathway enrichment analysis

GSEA genesets reactome pathway, KEGG pathway, and gene ontology can be analyzed by miRWalk3.0. R packages ‘clusterProfiler’ was also used to draw the graph of GO functional annotation and KEGG pathway enrichment. Significantly enriched pathways were sorted by p values less than 0.05. With ENCODE RUNX1 CHIP-Seq on human K562(ENCSR414TYY)(32, 33) optimal and conservative peaks annotation was analyzed by R packages ‘TxDb.Hsapiens. UCSC.hg38. knownGene’ and ‘ChIPseeker’, and ‘annotatePeak’ function were used to define gene promoter, which included 10kb upstream to 0.5kb downstream of the transcription start site (TSS).

RUNX1 regulates miRNA expression in AML

To identify potential roles of RUNX1 in AML, the differentially expressed miRNAs were identified from 14 TCGA AML samples by using DESeq2 algorithms. Compared the expression in AML samples with wildtype RUNX1, in those with mutant RUNX1 28 miRNAs were differently expressed (DEmiRNAs)with significance, in which 16 DEmiRNAs were up-regulated and 12 DEmiRNAs were down-regulated (Fig. 1 and Supplemental 2, significant differentially expression miRNAs and raw miRNAs expression array). Some miRNAs such as hsa-miR-1304, hsa-miR-139, hsa-miR-151a, hsa-miR-1515b, hsa-miR-4473, hsa-miR-873, hsa-miR-874, hsa-miR-877 were differential expressed in patients with NPM1+/FLT3+ AML versus control(34). hsa-miR-126 was involved in AML pathogenesis with cell cycle, proliferation, survival, and differentiation(35, 36).hsa-miR-551b was found to be upregulated in AML(37).

Significant Correlation of RUNX1 miRNAs with the Characteristics and Survival of Patients with AML

To explore the correlation of DEmiRNAs with RUNX1 in AML, the association of miRNAs and RUNX1 in expression was determined by Spearman algorithms. (Supplemental 3, significant correlation miRNAs to RUNX1, p < 0.005). We found 6 miRNAs were significantly correlated to RUNX1 in DEmiRNAs. hsa-miR-551b, has-miR-4473, and has-miR-873 were positive to RUNX1, hsa-miR-6877, has-miR-30c-2, and has-miR-30c-1 were negative to RUNX1. As shown in table I, by Cox regression analyses, the expression of 3 miRNAs (hsa-miR-873, hsa-miR-551b, and hsa-miR-4473) of these six miRNAs was correlated to patients’ overall survival (OS, p < 0.005), and patients with low expression had better treatment outcome (Fig. 2). These findings revealed that these three miRNAs could be involved in chemo-resistance by interfering expression of target mRNAs, and confer poor prognosis in AML.

RUNX1 regulates the MAPK signaling pathway through miRNAs

Check the previous IDs of hsa-miR-4473, hsa-miR-873 and hsa- miR-551b in miRBase database, we used hsa-miR-873 as hsa-miR-873-5p, and hsa-miR-551b as hsa-miR-551b-3p. miRNAs regulated gene expression at the post-transcriptional level by binding to 3’ UTR of the target mRNAs. To illustrate the mechanism on chemo-resistance of those 3 miRNAs associated with OS, miRWalk3.0 was applied to predict target genes of hsa-miR-4473, hsa-miR-551b-3p, and hsa-miR-873-5p. By target regions within 3’UTR and more than 0.95 of score, 4633 target genes of these three miRNAs were identified. The KEGG pathway enrichment of these target genes was analyzed with miRWalk3.0 (Supplemental 4, predictions result of DEmiRNAs from miRWalk3.0.). As presented in Fig. 3A-B, with the highest scores, MAPK signaling pathway(hsa04010), Ras signaling pathway(hsa04014), and Rap1 signaling pathway(hsa04015) were found by these enrichment assays. The target genes of three miRNAs were most significantly enriched in the MAPK signaling pathway. The target genes in MAPK pathway signaling were shown in Supplemental 5, graph of miRNA-target genes mapping in the MAPK signaling pathway. Total 95 miRNA-target genes by the red font mapping in the MAPK signaling pathway from KEGG Mapper, and these genes were involved in ERK, JNK, and p38 MAP kinase pathways (red color). Bioinformatics target prediction and pathway analysis results indicated RUNX1 could regulate MAPK signaling pathway through miRNAs, which implied that MAPK signaling could modulate chemo-response of AML patients with RUNX1 mutations.

RUNX1 regulates the MAPK signaling pathway directly

To investigate the control of RUNX1 on mRNA expression, genes associated with RUNX1 expression in AML carrying wild-type RUNX1 were gathered from GEO dataset GSE106291, and the most significant associated genes(p<0.05) were used for KEGG enrichment analysis (Supplemental 6, KEGG enrichment of GSE106291 and row wild-type RUNX1 expression array and Supplemental 7, dot plot of KEGG enrichment of RUNX1 correlate genes in GSE106291). It was noted that the MAPK signaling pathway was the highest enriched pathway, too. One hundred and eighty-two RUNX1-related genes were regulating MAPK signaling pathway (hsa04010, pvalue=0.0071). With RUNX1 ChIP-seq data of human K562 cell line in ENCODE (Fig. 4A), it was found that 35 genes of the optimal and conserved peaks with annotation were enriched to MAPK signaling pathway (Fig. 4B). Meanwhile, in human CD34+ cells there were 9727 genes within whose promoters there were RUNX1 binding motifs (TGTGGT/TGTGTT) (38). Among these genes, 172 genes were participating in MAPK signaling pathway (figure 4B). These results supported RUNX1 could transcriptionally control genes’ expression directly, which were components of MAPK signaling pathway.

As it was shown above, by CHIP-seq data of K562 and human CD34 cells, many genes of MAPK signaling pathway were recognized and potentially controlled by RUNX1. In the RUNX1-related genes, 61.90% genes mediated MAPK signaling pathway. More than 60% genes, which had RUNX1 peaks within 10kb of their promoters (CD34 cell and K562 cell) were involved in MAPK signaling pathway. In comparison, there was less than 30% genes genome-wide with RUNX1 peaks in their promoters. In DemiRNAs target genes, there were 95 genes involved in MAPK signaling pathway, more than 60% of genes (59/95) of which were overlapped with RUNX1-related genes. Further, 63 of 95 target genes of DemiRNAs had RUNX1 binding peaks in their promoters. It meant that RUNX1 could regulate their 63 genes’ expression indirectly and/or directly. Among these four subsets, there were 44 shared genes (Fig. 5 and Supplemental 8, RUNX1 related genes target to MAPK signaling pathway). Log-rank test of these 44 genes found three of them were significantly related to survival (Fig. 5D, p<0.01), including MEF2C, TGFBR2, and CACNB3. It indicated that RUNX1 regulated MAPK signaling pathway indirectly (miRNAs) and directly (mRNA).

RUNX1, as a critical transcription factor, played important roles in hematopoiesis, including maintaining the pool of hematopoietic stem cells, megakaryocyte differentiation, and lymphocyte development. With high-through technology, advances in the research of RUNX1 had been made. However, the mechanism of loss of function mutation in Runx1 required further investigation urgently. In this study, with transcriptome data of miRNAs and mRNAs from GDC-TCGA-LAML, to distinguish targets and uncover the underlying mechanism of RUNX1, differentially expressed targets and enriched pathways of these targets were revealed. hsa-miR-873, hsa-miR-4473, and hsa-miR-551b were suggested to be up-regulated by RUNX1, and their expression levels were positively correlated to OS. MEF2C, TGFBR2, CACNB3 were correlated to RUNX1, and their expression levels were positively correlated to OS, too. And target genes of these miRNAs modulated MAPK signaling pathway.

MAPKs is highly conserved protein kinases, including nineteen species MAPKKK (also named MAP3K or MEKK) were found, seven kinds of MAPKK (MAP2K) and fourteen kinds of MAPK. They transmit extracellular signals to intracellular to regulate cell growth, differentiation, apoptosis, migration, and so on. There are four distinct MAPK signaling cascades, including extracellular signal-related kinases (ERK1/2 and ERK5)，Jun amino-terminal kinases(JNK1/2/3)，and p38 kinase(39). JNK and p38 MAPK share roles in inflammation, apoptosis, and growth. The kinases of this pathway were dysregulated in various diseases, including cancers so that they could be used as biomarkers or therapeutic targets. MEK inhibitor (PD0325901) was tested in a phase II clinical trial for advanced non-small cell line lung cancer (NSCLC). Moreover, LGX818 is a BRAF inhibitor, which also has shown potential benefit for NSCLC patients. Inhibitor for p35 MAPK improved behavior and synaptic function of mouse with Alzheimer’s disease. Inhibition of JNK pathway provided a potent regulator for brain inflammation in neurological disorders(40, 41).

In AML, MAPK pathway was activated, and ERK phosphorylation was associated with cell growth(42), Flt3 ligand stimulation caused rapid activation of MAPK(43). It has been reported that combinatorial blockade of the MAPK pathways promoted cell death and suppresses proliferation in the majority of primary AML cells(44). Further, activated MAPK signaling contributed to primary therapeutic resistance to an IDH2 inhibitor in AML(45). pan-RAF inhibition, had been proved that it was a promising treatment strategy for AML(46). And in preclinical models, pexmetinib as a p38 MAPK inhibitor showed its potential for MDS/AML treatment (47). In AML, MAPK signaling pathway could be regulated by miR-181a, miR-497, and miR-143. And it has been identified that RUNX1 could affect the activity of MAPK signaling pathway miRNAs and mRNAs(48-50).

In this study, we found Six miRNAs were differently expressed with significance, and three of them were related to overall survival. Predicted target genes of these 3 miRNAs were highly enriched in MAPK signaling pathway by functional enrichment with miRWalk3.0. Besides, genes among RUNX1 associated genes directly regulated by RUNX1 were involved in MAPK signaling pathway, too. Taken together, RUNX1 modulated MAPK signaling pathway in AML. These results could explain the poor treatment outcome with mutant RUNX1, and enhance our understanding of the TF-miRNA regulatory associated with cell proliferation, apoptosis, and survival with MAPK signaling pathway.

Here, in GDC-TCGA-AML, we found new targets, including miRNAs and mRNAs, enriched in MAPK signaling pathway, and presumably proved RUNX1 regulated MAPK signaling pathway. It was concluded that this regulatory network was of great significance to the phenotype of AML with RUNX1 LOF mutations, and this also provided a potential target for clinical treatment in AML. At the same time, there were some shortcomings, the initial samples were relatively small and these research conclusions need to be verified by clinical specimens and animal models.

Acknowledgements

The authors thank Yu Liu for assistance with the bioinformatic analysis, from Shanghai Children's Medical Center, Shanghai Jiaotong University School of Medicine. This research was support by the National Natural Science Foundation of China (Grant No. 8187112228).

Cai SF and Levine RL: Genetic and epigenetic determinants of aml pathogenesis. Semin Hematol 56(2): 84-89, 2019. PMID, DOI: 10.1053/j.seminhematol.2018.08.001
DeWolf S and Tallman MS: How i treat relapsed or refractory aml. Blood 136(9): 1023-1032, 2020. PMID: PMC7453152, DOI: 10.1182/blood.2019001982
Mangaonkar AA and Patnaik MM: Hereditary predisposition to hematopoietic neoplasms: When bloodline matters for blood cancers. Mayo Clin Proc 95(7): 1482-1498, 2020. PMID, DOI: 10.1016/j.mayocp.2019.12.013
Babushok DV, Bessler M and Olson TS: Genetic predisposition to myelodysplastic syndrome and acute myeloid leukemia in children and young adults. Leuk Lymphoma 57(3): 520-536, 2016. PMID: PMC4798888, DOI: 10.3109/10428194.2015.1115041
Duployez N, Lejeune S, Renneville A and Preudhomme C: Myelodysplastic syndromes and acute leukemia with genetic predispositions: A new challenge for hematologists. Expert Rev Hematol 9(12): 1189-1202, 2016. PMID, DOI: 10.1080/17474086.2016.1257936
Hu J, Zhou GB, Wang ZY, Chen SJ and Chen Z: Mutant transcription factors and tyrosine kinases as therapeutic targets for leukemias: From acute promyelocytic leukemia to chronic myeloid leukemia and beyond. Adv Cancer Res 98(191-220, 2007. PMID, DOI: 10.1016/S0065-230X(06)98006-3
Schnittger S, Dicker F, Kern W, Wendland N, Sundermann J, Alpermann T, Haferlach C and Haferlach T: Runx1 mutations are frequent in de novo aml with noncomplex karyotype and confer an unfavorable prognosis. Blood 117(8): 2348-2357, 2011. PMID, DOI: 10.1182/blood-2009-11-255976
Zhen T, Kwon EM, Zhao L, Hsu J, Hyde RK, Lu Y, Alemu L, Speck NA and Liu PP: Chd7 deficiency delays leukemogenesis in mice induced by cbfb-myh11. Blood 130(22): 2431-2442, 2017. PMID: PMC5709785, DOI: 10.1182/blood-2017-04-780106
Zhen T, Cao Y, Ren G, Zhao L, Hyde RK, Lopez G, Feng D, Alemu L, Zhao K and Liu PP: Runx1 and cbfbeta-smmhc transactivate target genes together in abnormal myeloid progenitors for leukemia development. Blood 136(21): 2373-2385, 2020. PMID: PMC7685208, DOI: 10.1182/blood.2020007747
Harada Y and Harada H: Molecular mechanisms that produce secondary mds/aml by runx1/aml1 point mutations. J Cell Biochem 112(2): 425-432, 2011. PMID, DOI: 10.1002/jcb.22974
Cai X, Gaudet JJ, Mangan JK, Chen MJ, De Obaldia ME, Oo Z, Ernst P and Speck NA: Runx1 loss minimally impacts long-term hematopoietic stem cells. PLoS One 6(12): e28430, 2011. PMID: PMC3228772, DOI: 10.1371/journal.pone.0028430
Fu Y, Zhu H, Wu W, Xu J, Chen T, Xu B, Qian S, Li J and Liu P: Clinical significance of lymphoid enhancer-binding factor 1 expression in acute myeloid leukemia. Leuk Lymphoma 55(2): 371-377, 2014. PMID, DOI: 10.3109/10428194.2013.805759
Greif PA, Konstandin NP, Metzeler KH, Herold T, Pasalic Z, Ksienzyk B, Dufour A, Schneider F, Schneider S, Kakadia PM, Braess J, Sauerland MC, Berdel WE, Buchner T, Woermann BJ, Hiddemann W, Spiekermann K and Bohlander SK: Runx1 mutations in cytogenetically normal acute myeloid leukemia are associated with a poor prognosis and up-regulation of lymphoid genes. Haematologica 97(12): 1909-1915, 2012. PMID: PMC3590097, DOI: 10.3324/haematol.2012.064667
Cancer Genome Atlas Research N, Ley TJ, Miller C, Ding L, Raphael BJ, Mungall AJ, Robertson A, Hoadley K, Triche TJ, Jr., Laird PW, Baty JD, Fulton LL, Fulton R, Heath SE, Kalicki-Veizer J, Kandoth C, Klco JM, Koboldt DC, Kanchi KL, Kulkarni S, Lamprecht TL, Larson DE, Lin L, Lu C, McLellan MD, McMichael JF, Payton J, Schmidt H, Spencer DH, Tomasson MH, Wallis JW, Wartman LD, Watson MA, Welch J, Wendl MC, Ally A, Balasundaram M, Birol I, Butterfield Y, Chiu R, Chu A, Chuah E, Chun HJ, Corbett R, Dhalla N, Guin R, He A, Hirst C, Hirst M, Holt RA, Jones S, Karsan A, Lee D, Li HI, Marra MA, Mayo M, Moore RA, Mungall K, Parker J, Pleasance E, Plettner P, Schein J, Stoll D, Swanson L, Tam A, Thiessen N, Varhol R, Wye N, Zhao Y, Gabriel S, Getz G, Sougnez C, Zou L, Leiserson MD, Vandin F, Wu HT, Applebaum F, Baylin SB, Akbani R, Broom BM, Chen K, Motter TC, Nguyen K, Weinstein JN, Zhang N, Ferguson ML, Adams C, Black A, Bowen J, Gastier-Foster J, Grossman T, Lichtenberg T, Wise L, Davidsen T, Demchok JA, Shaw KR, Sheth M, Sofia HJ, Yang L, Downing JR and Eley G: Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med 368(22): 2059-2074, 2013. PMID: PMC3767041, DOI: 10.1056/NEJMoa1301689
Patel JP, Gonen M, Figueroa ME, Fernandez H, Sun Z, Racevskis J, Van Vlierberghe P, Dolgalev I, Thomas S, Aminova O, Huberman K, Cheng J, Viale A, Socci ND, Heguy A, Cherry A, Vance G, Higgins RR, Ketterling RP, Gallagher RE, Litzow M, van den Brink MR, Lazarus HM, Rowe JM, Luger S, Ferrando A, Paietta E, Tallman MS, Melnick A, Abdel-Wahab O and Levine RL: Prognostic relevance of integrated genetic profiling in acute myeloid leukemia. N Engl J Med 366(12): 1079-1089, 2012. PMID: PMC3545649, DOI: 10.1056/NEJMoa1112304
Sticht C, De La Torre C, Parveen A and Gretz N: Mirwalk: An online resource for prediction of microrna binding sites. Plos One 13(10), 2018. PMID, DOI: ARTN e0206239 10.1371/journal.pone.0206239
Dykes IM and Emanueli C: Transcriptional and post-transcriptional gene regulation by long non-coding rna. Genomics Proteomics Bioinformatics 15(3): 177-186, 2017. PMID: PMC5487525, DOI: 10.1016/j.gpb.2016.12.005
Lee YS and Dutta A: Micrornas in cancer. Annu Rev Pathol 4(199-227, 2009. PMID: PMC2769253, DOI: 10.1146/annurev.pathol.4.110807.092222
Mishra S, Yadav T and Rani V: Exploring mirna based approaches in cancer diagnostics and therapeutics. Crit Rev Oncol Hematol 98(12-23, 2016. PMID, DOI: 10.1016/j.critrevonc.2015.10.003
Ha M and Kim VN: Regulation of microrna biogenesis. Nat Rev Mol Cell Biol 15(8): 509-524, 2014. PMID, DOI: 10.1038/nrm3838
Rossetti S and Sacchi N: Runx1: A microrna hub in normal and malignant hematopoiesis. Int J Mol Sci 14(1): 1566-1588, 2013. PMID: PMC3565335, DOI: 10.3390/ijms14011566
Kohrs N, Kolodziej S, Kuvardina ON, Herglotz J, Yillah J, Herkt S, Piechatzek A, Salinas Riester G, Lingner T, Wichmann C, Bonig H, Seifried E, Platzbecker U, Medyouf H, Grez M and Lausen J: Mir144/451 expression is repressed by runx1 during megakaryopoiesis and disturbed by runx1/eto. PLoS Genet 12(3): e1005946, 2016. PMID: PMC4798443, DOI: 10.1371/journal.pgen.1005946
Lu X, Li X, He Q, Gao J, Gao Y, Liu B and Liu F: Mir-142-3p regulates the formation and differentiation of hematopoietic stem cells in vertebrates. Cell Res 23(12): 1356-1368, 2013. PMID: PMC3847575, DOI: 10.1038/cr.2013.145
Lam K, Muselman A, Du R, Harada Y, Scholl AG, Yan M, Matsuura S, Weng S, Harada H and Zhang DE: Hmga2 is a direct target gene of runx1 and regulates expansion of myeloid progenitors in mice. Blood 124(14): 2203-2212, 2014. PMID: PMC4183983, DOI: 10.1182/blood-2014-02-554543
Hoogenkamp M, Lichtinger M, Krysinska H, Lancrin C, Clarke D, Williamson A, Mazzarella L, Ingram R, Jorgensen H, Fisher A, Tenen DG, Kouskoff V, Lacaud G and Bonifer C: Early chromatin unfolding by runx1: A molecular explanation for differential requirements during specification versus maintenance of the hematopoietic gene expression program. Blood 114(2): 299-309, 2009. PMID: PMC2714206, DOI: 10.1182/blood-2008-11-191890
Austin GE, Zhao WG, Regmi A, Lu JP and Braun J: Identification of an upstream enhancer containing an aml1 site in the human myeloperoxidase (mpo) gene. Leuk Res 22(11): 1037-1048, 1998. PMID, DOI: 10.1016/s0145-2126(98)00105-2
Tang JL, Hou HA, Chen CY, Liu CY, Chou WC, Tseng MH, Huang CF, Lee FY, Liu MC, Yao M, Huang SY, Ko BS, Hsu SC, Wu SJ, Tsay W, Chen YC, Lin LI and Tien HF: Aml1/runx1 mutations in 470 adult patients with de novo acute myeloid leukemia: Prognostic implication and interaction with other gene alterations. Blood 114(26): 5352-5361, 2009. PMID, DOI: 10.1182/blood-2009-05-223784
Love MI, Huber W and Anders S: Moderated estimation of fold change and dispersion for rna-seq data with deseq2. Genome Biol 15(12): 550, 2014. PMID: PMC4302049, DOI: 10.1186/s13059-014-0550-8
Dweep H and Gretz N: Mirwalk2.0: A comprehensive atlas of microrna-target interactions. Nat Methods 12(8): 697, 2015. PMID, DOI: 10.1038/nmeth.3485
Parveen A, Gretz N and Dweep H: Obtaining mirna-target interaction information from mirwalk2.0. Curr Protoc Bioinformatics 55(12 15 11-12 15 27, 2016. PMID, DOI: 10.1002/cpbi.14
Sticht C, De La Torre C, Parveen A and Gretz N: Mirwalk: An online resource for prediction of microrna binding sites. PLoS One 13(10): e0206239, 2018. PMID: PMC6193719, DOI: 10.1371/journal.pone.0206239
Consortium EP: An integrated encyclopedia of DNA elements in the human genome. Nature 489(7414): 57-74, 2012. PMID: PMC3439153, DOI: 10.1038/nature11247
Consortium EP: A user's guide to the encyclopedia of DNA elements (encode). PLoS Biol 9(4): e1001046, 2011. PMID: PMC3079585, DOI: 10.1371/journal.pbio.1001046
Pandita A, Ramadas P, Poudel A, Saad N, Anand A, Basnet A, Wang D, Middleton F and Gilligan DM: Differential expression of mirnas in acute myeloid leukemia quantified by nextgen sequencing of whole blood samples. PLoS One 14(3): e0213078, 2019. PMID: PMC6426230, DOI: 10.1371/journal.pone.0213078
Dorrance AM, Neviani P, Ferenchak GJ, Huang X, Nicolet D, Maharry KS, Ozer HG, Hoellarbauer P, Khalife J, Hill EB, Yadav M, Bolon BN, Lee RJ, Lee LJ, Croce CM, Garzon R, Caligiuri MA, Bloomfield CD and Marcucci G: Targeting leukemia stem cells in vivo with antagomir-126 nanoparticles in acute myeloid leukemia. Leukemia 29(11): 2143-2153, 2015. PMID: PMC4633325, DOI: 10.1038/leu.2015.139
Wallace JA and O'Connell RM: Micrornas and acute myeloid leukemia: Therapeutic implications and emerging concepts. Blood 130(11): 1290-1301, 2017. PMID: PMC5600138, DOI: 10.1182/blood-2016-10-697698
Zhang H, Zhang C, Feng R, Zhang H, Gao M and Ye L: Investigating the microrna-mrna regulatory network in acute myeloid leukemia. Oncol Lett 14(4): 3981-3988, 2017. PMID: PMC5620483, DOI: 10.3892/ol.2017.6686
Beck D, Thoms JA, Perera D, Schutte J, Unnikrishnan A, Knezevic K, Kinston SJ, Wilson NK, O'Brien TA, Gottgens B, Wong JW and Pimanda JE: Genome-wide analysis of transcriptional regulators in human hspcs reveals a densely interconnected network of coding and noncoding genes. Blood 122(14): e12-22, 2013. PMID, DOI: 10.1182/blood-2013-03-490425
Sun Y, Liu WZ, Liu T, Feng X, Yang N and Zhou HF: Signaling pathway of mapk/erk in cell proliferation, differentiation, migration, senescence and apoptosis. J Recept Signal Transduct Res 35(6): 600-604, 2015. PMID, DOI: 10.3109/10799893.2015.1030412
Pradhan R, Singhvi G, Dubey SK, Gupta G and Dua K: Mapk pathway: A potential target for the treatment of non-small-cell lung carcinoma. Future Med Chem 11(8): 793-795, 2019. PMID, DOI: 10.4155/fmc-2018-0468
Kim EK and Choi EJ: Pathological roles of mapk signaling pathways in human diseases. Biochim Biophys Acta 1802(4): 396-405, 2010. PMID, DOI: 10.1016/j.bbadis.2009.12.009
Ricciardi MR, McQueen T, Chism D, Milella M, Estey E, Kaldjian E, Sebolt-Leopold J, Konopleva M and Andreeff M: Quantitative single cell determination of erk phosphorylation and regulation in relapsed and refractory primary acute myeloid leukemia. Leukemia 19(9): 1543-1549, 2005. PMID, DOI: 10.1038/sj.leu.2403859
Hayakawa F, Towatari M, Kiyoi H, Tanimoto M, Kitamura T, Saito H and Naoe T: Tandem-duplicated flt3 constitutively activates stat5 and map kinase and introduces autonomous cell growth in il-3-dependent cell lines. Oncogene 19(5): 624-631, 2000. PMID, DOI: 10.1038/sj.onc.1203354
Han L, Zhang Q, Dail M, Shi C, Cavazos A, Ruvolo VR, Zhao Y, Kim E, Rahmani M, Mak DH, Jin SS, Chen J, Phillips DC, Koller PB, Jacamo R, Burks JK, DiNardo C, Daver N, Jabbour E, Wang J, Kantarjian HM, Andreeff M, Grant S, Leverson JD, Sampath D and Konopleva M: Concomitant targeting of bcl2 with venetoclax and mapk signaling with cobimetinib in acute myeloid leukemia models. Haematologica 105(3): 697-707, 2020. PMID: PMC7049339, DOI: 10.3324/haematol.2018.205534
Amatangelo MD, Quek L, Shih A, Stein EM, Roshal M, David MD, Marteyn B, Farnoud NR, de Botton S, Bernard OA, Wu B, Yen KE, Tallman MS, Papaemmanuil E, Penard-Lacronique V, Thakurta A, Vyas P and Levine RL: Enasidenib induces acute myeloid leukemia cell differentiation to promote clinical response. Blood 130(6): 732-741, 2017. PMID: PMC5553578, DOI: 10.1182/blood-2017-04-779447
Tambe M, Karjalainen E, Vaha-Koskela M, Bulanova D, Gjertsen BT, Kontro M, Porkka K, Heckman CA and Wennerberg K: Pan-raf inhibition induces apoptosis in acute myeloid leukemia cells and synergizes with bcl2 inhibition. Leukemia 34(12): 3186-3196, 2020. PMID, DOI: 10.1038/s41375-020-0972-0
Bachegowda L, Morrone K, Winski SL, Mantzaris I, Bartenstein M, Ramachandra N, Giricz O, Sukrithan V, Nwankwo G, Shahnaz S, Bhagat T, Bhattacharyya S, Assal A, Shastri A, Gordon-Mitchell S, Pellagatti A, Boultwood J, Schinke C, Yu Y, Guha C, Rizzi J, Garrus J, Brown S, Wollenberg L, Hogeland G, Wright D, Munson M, Rodriguez M, Gross S, Chantry D, Zou Y, Platanias L, Burgess LE, Pradhan K, Steidl U and Verma A: Pexmetinib: A novel dual inhibitor of tie2 and p38 mapk with efficacy in preclinical models of myelodysplastic syndromes and acute myeloid leukemia. Cancer Res 76(16): 4841-4849, 2016. PMID: PMC5398415, DOI: 10.1158/0008-5472.CAN-15-3062
Hartmann JU, Brauer-Hartmann D, Kardosova M, Wurm AA, Wilke F, Schodel C, Gerloff D, Katzerke C, Krakowsky R, Namasu CY, Bill M, Schwind S, Muller-Tidow C, Niederwieser D, Alberich-Jorda M and Behre G: Microrna-143 targets erk5 in granulopoiesis and predicts outcome of patients with acute myeloid leukemia. Cell Death Dis 9(8): 814, 2018. PMID: PMC6062564, DOI: 10.1038/s41419-018-0837-x
Kang H, Tong C, Li C and Luo J: Mir-497 plays a key role in tanshinone iia-attenuated proliferation in oci-aml3 cells via the mapk/erk1/2 pathway. Cytotechnology 72(3): 427-432, 2020. PMID: PMC7225231, DOI: 10.1007/s10616-020-00389-5
Huang X, Schwind S, Santhanam R, Eisfeld AK, Chiang CL, Lankenau M, Yu B, Hoellerbauer P, Jin Y, Tarighat SS, Khalife J, Walker A, Perrotti D, Bloomfield CD, Wang H, Lee RJ, Lee LJ and Marcucci G: Targeting the ras/mapk pathway with mir-181a in acute myeloid leukemia. Oncotarget 7(37): 59273-59286, 2016. PMID: PMC5312311, DOI: 10.18632/oncotarget.11150

Table I: Univariate Cox regression analyses in the TCGA filtered cohorts.

type	beta	HR (95% CI for HR)	wald.test	p.value
hsa-mir-551b	0.027	1 (1-1)	6.3	0.012
hsa-mir-4473	0.56	1.7 (1.2-2.6)	7	0.008
hsa-mir-873	2.5	12 (1.6-85)	5.9	0.015
hsa-mir-6877	-1.7	0.19 (0.027-1.3)	2.8	0.096
hsa-mir-30c-2	-0.00077	1 (1-1)	0.87	0.35
hsa-mir-30c-1	-0.00064	1 (1-1)	0.68	0.41

Download PDF

Version 1

posted

You are reading this latest preprint version

MAPK Signaling Pathway was Dysregulated in Acute Myeloid Leukemia with RUNX1 Mutations

Status:

Version 1

Abstract

Figures

Introduction

Materials and methods

Data preparation

Statistical methods

Prediction of miRNA target genes and identification of RUNX1 target genes

Functional annotation and pathway enrichment analysis

Results

RUNX1 regulates miRNA expression in AML

Significant Correlation of RUNX1 miRNAs with the Characteristics and Survival of Patients with AML

RUNX1 regulates the MAPK signaling pathway through miRNAs

RUNX1 regulates the MAPK signaling pathway directly

Discussion

Conclusions

Declarations

Acknowledgements

References

Tables

Supplementary Files

Status:

Version 1